Search

Public health emergencies like SARS, MERS, and COVID-19 have prioritized surveillance of zoonotic coronaviruses, resulting in extensive genomic characterization of coronavirus diversity in bats. Sequencing viral genomes directly from animal specimens remains a laboratory challenge, however, and most bat coronaviruses have been characterized solely by PCR amplification of small regions from the best-conserved gene. This has resulted in limited phylogenetic resolution and left viral genetic factors relevant to threat assessment undescribed. In this study, we evaluated whether a technique called hybridization probe capture can achieve more extensive genome recovery from surveillance specimens. Using a custom panel of 20,000 probes, we captured and sequenced coronavirus genomic material in 21 swab specimens collected from bats in the Democratic Republic of the Congo. For 15 of these specimens, probe capture recovered more genome sequence than had been previously generated with standard amplicon sequencing protocols, providing a median 6.1-fold improvement (ranging up to 69.1-fold). Probe capture data also identified five novel alpha- and betacoronaviruses in these specimens, and their full genomes were recovered with additional deep sequencing. Based on these experiences, we discuss how probe capture could be effectively operationalized alongside other sequencing technologies for high-throughput, genomics-based discovery and surveillance of bat coronaviruses., Competing Interests: KK, KL, DS, JR, AS, ML, NL, JM, NP, AC No competing interests declared, CK, FM, PM, IN, FN, JA, JA, CM, ER, DM were employed by Metabiota Inc, MM, KS, CL are employees of Labyrinth Global Health Inc and were employed by Metabiota Inc, AG is an employee of Development Alternatives Inc and was employed by Metabiota Inc, DJ is an employee of Nyati Health Consulting and was employed by Metabiota Inc, NW is an employee of Metabiota Inc, (© 2022, Kuchinski et al.)

Zoonotic spillover of animal viruses into human populations is a continuous and increasing public health risk. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) highlights the global impact of emergence. Considering the history and diversity of coronaviruses (CoVs), especially in bats, SARS-CoV-2 will likely not be the last to spillover from animals into human populations. We sampled and tested wildlife in the Central African country Cameroon to determine which CoVs are circulating and how they relate to previously detected human and animal CoVs. We collected animal and ecological data at sampling locations and used family-level consensus PCR combined with amplicon sequencing for virus detection. Between 2003 and 2018, samples were collected from 6,580 animals of several different orders. CoV RNA was detected in 175 bats, a civet, and a shrew. The CoV RNAs detected in the bats represented 17 different genetic clusters, coinciding with alpha ( n  = 8) and beta ( n  = 9) CoVs. Sequences resembling human CoV-229E (HCoV-229E) were found in 40 Hipposideridae bats. Phylogenetic analyses place the human-derived HCoV-229E isolates closest to those from camels in terms of the S and N genes but closest to isolates from bats for the envelope, membrane, and RNA-dependent RNA polymerase genes. The CoV RNA positivity rate in bats varied significantly (P < 0.001) between the wet (8.2 per cent) and dry seasons (4.5 per cent). Most sampled species accordingly had a wet season high and dry season low, while for some the opposite was found. Eight of the suspected CoV species of which we detected RNA appear to be entirely novel CoV species, which suggests that CoV diversity in African wildlife is still rather poorly understood. The detection of multiple different variants of HCoV-229E-like viruses supports the bat reservoir hypothesis for this virus, with the phylogenetic results casting some doubt on camels as an intermediate host. The findings also support the previously proposed influence of ecological factors on CoV circulation, indicating a high level of underlying complexity to the viral ecology. These results indicate the importance of investing in surveillance activities among wild animals to detect all potential threats as well as sentinel surveillance among exposed humans to determine emerging threats., (© The Author(s) 2021. Published by Oxford University Press.)

Enteroviruses infect humans and animals and can cause disease, and some may be transmitted across species barriers. We tested Central African wildlife and found Enterovirus RNA in primates (17) and rodents (2). Some sequences were very similar, while others were dissimilar to known species, highlighting the underexplored enterovirus diversity in wildlife.

Coronaviruses play an important role as pathogens of humans and animals, and the emergence of epidemics like SARS, MERS and COVID-19 is closely linked to zoonotic transmission events primarily from wild animals. Bats have been found to be an important source of coronaviruses with some of them having the potential to infect humans, with other animals serving as intermediate or alternate hosts or reservoirs. Host diversity may be an important contributor to viral diversity and thus the potential for zoonotic events. To date, limited research has been done in Africa on this topic, in particular in the Congo Basin despite frequent contact between humans and wildlife in this region. We sampled and, using consensus coronavirus PCR-primers, tested 3,561 wild animals for coronavirus RNA. The focus was on bats (38%), rodents (38%), and primates (23%) that posed an elevated risk for contact with people, and we found coronavirus RNA in 121 animals, of which all but two were bats. Depending on the taxonomic family, bats were significantly more likely to be coronavirus RNA-positive when sampled either in the wet (Pteropodidae and Rhinolophidae) or dry season (Hipposideridae, Miniopteridae, Molossidae, and Vespertilionidae). The detected RNA sequences correspond to 15 alpha- and 6 betacoronaviruses, with some of them being very similar (>95% nucleotide identities) to known coronaviruses and others being more unique and potentially representing novel viruses. In seven of the bats, we detected RNA most closely related to sequences of the human common cold coronaviruses 229E or NL63 (>80% nucleotide identities). The findings highlight the potential for coronavirus spillover, especially in regions with a high diversity of bats and close human contact, and reinforces the need for ongoing surveillance., Competing Interests: Metabiota (https://metabiota.com), Labyrinth Global Health (https://www.labyrinthgh.com), and Mosaic (https://mosaic.cm) are private contractors who receive funds from international donor organizations to conduct technical assistance and research. Metabiota employs or employed CK, FM, PMK, INL, FNK, JAL, PM, MM, UT, AG, BSS, DJM, JAA, CM, DOJ, KS, NDW, EMR, and CEL, Mosaic employs ML, and Labyrinth Global Health employs KS, MM and CEL. There are no patents, products in development or marketed products associated with this research to declare. The association with commercial entities that some of the authors have does not alter our adherence to PLOS ONE policies on sharing data and materials.

Adenoviruses (AdVs) are diverse pathogens of humans and animals, with several dozen bat AdVs already identified. Considering that over 100 human AdVs are known, and the huge diversity of bat species, many bat AdVs likely remain undiscovered. To learn more about AdV prevalence, diversity and evolution, we sampled and tested bats in Cameroon using several PCR assays for viral and host DNA. AdV DNA was detected in 14 % of the 671 sampled animals belonging to 37 different bat species. There was a correlation between species roosting in larger groups and AdV DNA detection. The detected AdV DNA belonged to between 28 and 44 different, mostly previously unknown, mastadenovirus species. The novel isolates are phylogenetically diverse and while some cluster with known viruses, others appear to form divergent new clusters. The phylogenetic tree of novel and previously known bat AdVs does not mirror that of the various host species, but does contain structures consistent with a degree of virus-host co-evolution. Given that closely related isolates were found in different host species, it seems likely that at least some bat AdVs have jumped species barriers, probably in the more recent past; however, the tree is also consistent with such events having taken place throughout bat AdV evolution. AdV diversity was highest in bat species roosting in large groups. The study significantly increased the diversity of AdVs known to be harboured by bats, and suggests that host behaviours, such as roosting size, may be what limits some AdVs to one species rather than an inability of AdVs to infect other related hosts.

Coronaviruses can become zoonotic, as in the case of COVID-19, and hunting, sale, and consumption of wild animals in Southeast Asia increases the risk for such incidents. We sampled and tested rodents (851) and other mammals and found betacoronavirus RNA in 12 rodents. The sequences belong to two separate genetic clusters and are closely related to those of known rodent coronaviruses detected in the region and distantly related to those of human coronaviruses OC43 and HKU1. Considering the close human-wildlife contact with many species in and beyond the region, a better understanding of virus diversity is urgently needed for the mitigation of future risks.

Dengue fever is an understudied disease in many parts of Africa and little is known about its prevalence in Cameroon. We tested blood from 629 individuals from the South Region of Cameroon, collected over the course of one year, for flavivirus RNA using conventional broad range PCR. Flavivirus RNA corresponding to dengue virus (DENV) serotype 1 was identified in two individuals who were also diagnosed with malaria. This finding confirms previous reports that indicate the presence of low-level circulation of DENV in Cameroon and supports the concern that dengue fever may be underdiagnosed due to more prevalent diseases that have similar symptomology and insufficient diagnostic capacity.

It is generally believed that exchange of secondary metabolite biosynthetic gene clusters (BGCs) among closely related bacteria is an important driver of BGC evolution and diversification. Applying this idea may help researchers efficiently connect many BGCs to their products and characterize the products' roles in various environments. However, existing genetic tools support only a small fraction of these efforts. Here, we present the development of chassis-independent recombinase-assisted genome engineering (CRAGE), which enables single-step integration of large, complex BGC constructs directly into the chromosomes of diverse bacteria with high accuracy and efficiency. To demonstrate the efficacy of CRAGE, we expressed three known and six previously identified but experimentally elusive non-ribosomal peptide synthetase (NRPS) and NRPS-polyketide synthase (PKS) hybrid BGCs from Photorhabdus luminescens in 25 diverse γ-Proteobacteria species. Successful activation of six BGCs identified 22 products for which diversity and yield were greater when the BGCs were expressed in strains closely related to the native strain than when they were expressed in either native or more distantly related strains. Activation of these BGCs demonstrates the feasibility of exploiting their underlying catalytic activity and plasticity, and provides evidence that systematic approaches based on CRAGE will be useful for discovering and identifying previously uncharacterized metabolites.

Rodent adenoviruses are important models for human disease. In contrast to the over 70 adenovirus types isolated from humans, few rodent adenoviruses are known, despite the vast diversity of rodent species. PCR and Sanger sequencing were used to investigate adenovirus diversity in wild rodents and shrews in Cameroon. Adenovirus DNA was detected in 13.8% of animals (n = 218). All detected sequences differ from known adenovirus types by more than 10% at the amino acid level, thus indicating up to 14 novel adenovirus species. These results highlight the diversity of rodent adenoviruses, their phylogeny, and opportunities for studying alternative adenovirus rodent models.

In the version of this Article originally published, the bat species for 12 individuals were incorrectly identified in Supplementary Table 1 and 2. After resequencing the MT-CytB and MT-CO1 segments and reviewing the data, the authors have corrected the errors for these 12 animals. In the amended version of the Supplementary Information, Supplementary Tables 1 and 2 have been replaced to include the corrected host species information. None of the 12 bats affected were positive for the Bombali virus, and the conclusions of the study are therefore unchanged.

Tomato variety Hawaii 7996 is resistant to the soil-borne pathogen Ralstonia solanacearum, whereas the Moneymaker variety is susceptible to the pathogen. To evaluate whether plant-associated microorganisms have a role in disease resistance, we analyzed the rhizosphere microbiomes of both varieties in a mesocosm experiment. Microbiome structures differed between the two cultivars. Transplantation of rhizosphere microbiota from resistant plants suppressed disease symptoms in susceptible plants. Comparative analyses of rhizosphere metagenomes from resistant and susceptible plants enabled the identification and assembly of a flavobacterial genome that was far more abundant in the resistant plant rhizosphere microbiome than in that of the susceptible plant. We cultivated this flavobacterium, named TRM1, and found that it could suppress R. solanacearum-disease development in a susceptible plant in pot experiments. Our findings reveal a role for native microbiota in protecting plants from microbial pathogens, and our approach charts a path toward the development of probiotics to ameliorate plant diseases.

Here we describe the complete genome of a new ebolavirus, Bombali virus (BOMV) detected in free-tailed bats in Sierra Leone (little free-tailed (Chaerephon pumilus) and Angolan free-tailed (Mops condylurus)). The bats were found roosting inside houses, indicating the potential for human transmission. We show that the viral glycoprotein can mediate entry into human cells. However, further studies are required to investigate whether exposure has actually occurred or if BOMV is pathogenic in humans.

Bocaparvoviruses are members of the family Parvovirinae and human bocaviruses have been found to be associated with respiratory and gastrointestinal disease. There are four known human bocaviruses, as well as several distinct ones in great apes. The goal of the presented study was to detect other non-human primate (NHP) bocaviruses in NHP species in the Democratic Republic of the Congo using conventional broad-range PCR. We found bocavirus DNA in blood and tissues samples in 6 out of 620 NHPs, and all isolates showed very high identity (>97 %) with human bocaviruses 2 or 3. These findings suggest cross-species transmission of bocaviruses between humans and NHPs.

Objective: Herpesviruses belong to a diverse order of large DNA viruses that can cause diseases in humans and animals. With the goal of gathering information about the distribution and diversity of herpesviruses in wild rodent and shrew species in central Africa, animals in Cameroon and the Democratic Republic of the Congo were sampled and tested by PCR for the presence of herpesvirus DNA., Methods: A broad range PCRs targeting either the Polymerase or the terminase gene were used for virus detection. Amplified products from PCR were sequenced and isolates analysed for phylogenetic placement., Results: Overall, samples of 1,004 animals of various rodent and shrew species were tested and 24 were found to be positive for herpesvirus DNA. Six of these samples contained strains of known viruses, while the other positive samples revealed DNA sequences putatively belonging to 11 previously undescribed herpesviruses. The new isolates are beta- and gammaherpesviruses and the shrew isolates appear to form a separate cluster within the Betaherpesvirinae subfamily., Conclusion: The diversity of viruses detected is higher than in similar studies in Europe and Asia. The high diversity of rodent and shrew species occurring in central Africa may be the reason for a higher diversity in herpesviruses in this area., (© 2018 S. Karger AG, Basel.)

Methods for assembly, taxonomic profiling and binning are key to interpreting metagenome data, but a lack of consensus about benchmarking complicates performance assessment. The Critical Assessment of Metagenome Interpretation (CAMI) challenge has engaged the global developer community to benchmark their programs on highly complex and realistic data sets, generated from ∼700 newly sequenced microorganisms and ∼600 novel viruses and plasmids and representing common experimental setups. Assembly and genome binning programs performed well for species represented by individual genomes but were substantially affected by the presence of related strains. Taxonomic profiling and binning programs were proficient at high taxonomic ranks, with a notable performance decrease below family level. Parameter settings markedly affected performance, underscoring their importance for program reproducibility. The CAMI results highlight current challenges but also provide a roadmap for software selection to answer specific research questions.

Background: Ruminants are important contributors to global methane emissions via microbial fermentation in their reticulo-rumens. This study is part of a larger program, characterising the rumen microbiomes of sheep which vary naturally in methane yield (g CH 4 /kg DM/day) and aims to define differences in microbial communities, and in gene and transcript abundances that can explain the animal methane phenotype., Methods: Rumen microbiome metagenomic and metatranscriptomic data were analysed by Gene Set Enrichment, sparse partial least squares regression and the Wilcoxon Rank Sum test to estimate correlations between specific KEGG bacterial pathways/genes and high methane yield in sheep. KEGG genes enriched in high methane yield sheep were reassembled from raw reads and existing contigs and analysed by MEGAN to predict their phylogenetic origin. Protein coding sequences from Succinivibrio dextrinosolvens strains were analysed using Effective DB to predict bacterial type III secreted proteins. The effect of S. dextrinosolvens strain H5 growth on methane formation by rumen methanogens was explored using co-cultures., Results: Detailed analysis of the rumen microbiomes of high methane yield sheep shows that gene and transcript abundances of bacterial type III secretion system genes are positively correlated with methane yield in sheep. Most of the bacterial type III secretion system genes could not be assigned to a particular bacterial group, but several genes were affiliated with the genus Succinivibrio, and searches of bacterial genome sequences found that strains of S. dextrinosolvens were part of a small group of rumen bacteria that encode this type of secretion system. In co-culture experiments, S. dextrinosolvens strain H5 showed a growth-enhancing effect on a methanogen belonging to the order Methanomassiliicoccales, and inhibition of a representative of the Methanobrevibacter gottschalkii clade., Conclusions: This is the first report of bacterial type III secretion system genes being associated with high methane emissions in ruminants, and identifies these secretions systems as potential new targets for methane mitigation research. The effects of S. dextrinosolvens on the growth of rumen methanogens in co-cultures indicate that bacteria-methanogen interactions are important modulators of methane production in ruminant animals.

We report the identification of novel tRNA species with 12-base pair amino-acid acceptor branches composed of longer acceptor stem and shorter T-stem. While canonical tRNAs have a 7/5 configuration of the branch, the novel tRNAs have either 8/4 or 9/3 structure. They were found during the search for selenocysteine tRNAs in terabytes of genome, metagenome and metatranscriptome sequences. Certain bacteria and their phages employ the 8/4 structure for serine and histidine tRNAs, while minor cysteine and selenocysteine tRNA species may have a modified 8/4 structure with one bulge nucleotide. In Acidobacteria, tRNAs with 8/4 and 9/3 structures may function as missense and nonsense suppressor tRNAs and/or regulatory noncoding RNAs. In δ-proteobacteria, an additional cysteine tRNA with an 8/4 structure mimics selenocysteine tRNA and may function as opal suppressor. We examined the potential translation function of suppressor tRNA species in Escherichia coli; tRNAs with 8/4 or 9/3 structures efficiently inserted serine, alanine and cysteine in response to stop and sense codons, depending on the identity element and anticodon sequence of the tRNA. These findings expand our view of how tRNA, and possibly the genetic code, is diversified in nature., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The evolution of body shape is thought to be tightly coupled to changes in regulatory sequences, but specific molecular events associated with major morphological transitions in vertebrates have remained elusive. We identified snake-specific sequence changes within an otherwise highly conserved long-range limb enhancer of Sonic hedgehog (Shh). Transgenic mouse reporter assays revealed that the in vivo activity pattern of the enhancer is conserved across a wide range of vertebrates, including fish, but not in snakes. Genomic substitution of the mouse enhancer with its human or fish ortholog results in normal limb development. In contrast, replacement with snake orthologs caused severe limb reduction. Synthetic restoration of a single transcription factor binding site lost in the snake lineage reinstated full in vivo function to the snake enhancer. Our results demonstrate changes in a regulatory sequence associated with a major body plan transition and highlight the role of enhancers in morphological evolution. PAPERCLIP., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Background: Enteric fermentation by farmed ruminant animals is a major source of methane and constitutes the second largest anthropogenic contributor to global warming. Reducing methane emissions from ruminants is needed to ensure sustainable animal production in the future. Methane yield varies naturally in sheep and is a heritable trait that can be used to select animals that yield less methane per unit of feed eaten. We previously demonstrated elevated expression of hydrogenotrophic methanogenesis pathway genes of methanogenic archaea in the rumens of high methane yield (HMY) sheep compared to their low methane yield (LMY) counterparts. Methane production in the rumen is strongly connected to microbial hydrogen production through fermentation processes. In this study, we investigate the contribution that rumen bacteria make to methane yield phenotypes in sheep., Results: Using deep sequence metagenome and metatranscriptome datasets in combination with 16S rRNA gene amplicon sequencing from HMY and LMY sheep, we show enrichment of lactate-producing Sharpea spp. in LMY sheep bacterial communities. Increased gene and transcript abundances for sugar import and utilisation and production of lactate, propionate and butyrate were also observed in LMY animals. Sharpea azabuensis and Megasphaera spp. act as important drivers of lactate production and utilisation according to phylogenetic analysis and read mappings., Conclusions: Our findings show that the rumen microbiome in LMY animals supports a rapid heterofermentative growth, leading to lactate production. We postulate that lactate is subsequently metabolised mainly to butyrate in LMY animals, producing 2 mol of hydrogen and 0.5 mol of methane per mol hexose, which represents 24 % less than the 0.66 mol of methane formed from the 2.66 mol of hydrogen produced if hexose fermentation was directly to acetate and butyrate. These findings are consistent with the theory that a smaller rumen size with a higher turnover rate, where rapid heterofermentative growth would be an advantage, results in lower hydrogen production and lower methane formation. Together with previous methanogen gene expression data, this builds a strong concept of how animal traits and microbial communities shape the methane phenotype in sheep.

Whole-genome sequencing is identifying growing numbers of non-coding variants in human disease studies, but the lack of accurate functional annotations prevents their interpretation. We describe the genome-wide landscape of distant-acting enhancers active in the developing and adult human heart, an organ whose impairment is a predominant cause of mortality and morbidity. Using integrative analysis of >35 epigenomic data sets from mouse and human pre- and postnatal hearts we created a comprehensive reference of >80,000 putative human heart enhancers. To illustrate the importance of enhancers in the regulation of genes involved in heart disease, we deleted the mouse orthologs of two human enhancers near cardiac myosin genes. In both cases, we observe in vivo expression changes and cardiac phenotypes consistent with human heart disease. Our study provides a comprehensive catalogue of human heart enhancers for use in clinical whole-genome sequencing studies and highlights the importance of enhancers for cardiac function.

Selenocysteine (Sec or U) is encoded by UGA, a stop codon reassigned by a Sec-specific elongation factor and a distinctive RNA structure. To discover possible code variations in extant organisms we analyzed 6.4 trillion base pairs of metagenomic sequences and 24 903 microbial genomes for tRNA(Sec) species. As expected, UGA is the predominant Sec codon in use. We also found tRNA(Sec) species that recognize the stop codons UAG and UAA, and ten sense codons. Selenoprotein synthesis programmed by UAG in Geodermatophilus and Blastococcus, and by the Cys codon UGU in Aeromonas salmonicida was confirmed by metabolic labeling with (75) Se or mass spectrometry. Other tRNA(Sec) species with different anticodons enabled E. coli to synthesize active formate dehydrogenase H, a selenoenzyme. This illustrates the ease by which the genetic code may evolve new coding schemes, possibly aiding organisms to adapt to changing environments, and show the genetic code is much more flexible than previously thought., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

DNA methylation acts in concert with restriction enzymes to protect the integrity of prokaryotic genomes. Studies in a limited number of organisms suggest that methylation also contributes to prokaryotic genome regulation, but the prevalence and properties of such non-restriction-associated methylation systems remain poorly understood. Here, we used single molecule, real-time sequencing to map DNA modifications including m6A, m4C, and m5C across the genomes of 230 diverse bacterial and archaeal species. We observed DNA methylation in nearly all (93%) organisms examined, and identified a total of 834 distinct reproducibly methylated motifs. This data enabled annotation of the DNA binding specificities of 620 DNA Methyltransferases (MTases), doubling known specificities for previously hard to study Type I, IIG and III MTases, and revealing their extraordinary diversity. Strikingly, 48% of organisms harbor active Type II MTases with no apparent cognate restriction enzyme. These active 'orphan' MTases are present in diverse bacterial and archaeal phyla and show motif specificities and methylation patterns consistent with functions in gene regulation and DNA replication. Our results reveal the pervasive presence of DNA methylation throughout the prokaryotic kingdoms, as well as the diversity of sequence specificities and potential functions of DNA methylation systems.

Analysis of the increasing wealth of metagenomic data collected from diverse environments can lead to the discovery of novel branches on the tree of life. Here we analyse 5.2 Tb of metagenomic data collected globally to discover a novel bacterial phylum ('Candidatus Kryptonia') found exclusively in high-temperature pH-neutral geothermal springs. This lineage had remained hidden as a taxonomic 'blind spot' because of mismatches in the primers commonly used for ribosomal gene surveys. Genome reconstruction from metagenomic data combined with single-cell genomics results in several high-quality genomes representing four genera from the new phylum. Metabolic reconstruction indicates a heterotrophic lifestyle with conspicuous nutritional deficiencies, suggesting the need for metabolic complementarity with other microbes. Co-occurrence patterns identifies a number of putative partners, including an uncultured Armatimonadetes lineage. The discovery of Kryptonia within previously studied geothermal springs underscores the importance of globally sampled metagenomic data in detection of microbial novelty, and highlights the extraordinary diversity of microbial life still awaiting discovery.

Yakutia, Sakha Republic, in the Siberian Far East, represents one of the coldest places on Earth, with winter record temperatures dropping below -70 °C. Nevertheless, Yakutian horses survive all year round in the open air due to striking phenotypic adaptations, including compact body conformations, extremely hairy winter coats, and acute seasonal differences in metabolic activities. The evolutionary origins of Yakutian horses and the genetic basis of their adaptations remain, however, contentious. Here, we present the complete genomes of nine present-day Yakutian horses and two ancient specimens dating from the early 19th century and ∼5,200 y ago. By comparing these genomes with the genomes of two Late Pleistocene, 27 domesticated, and three wild Przewalski's horses, we find that contemporary Yakutian horses do not descend from the native horses that populated the region until the mid-Holocene, but were most likely introduced following the migration of the Yakut people a few centuries ago. Thus, they represent one of the fastest cases of adaptation to the extreme temperatures of the Arctic. We find cis-regulatory mutations to have contributed more than nonsynonymous changes to their adaptation, likely due to the comparatively limited standing variation within gene bodies at the time the population was founded. Genes involved in hair development, body size, and metabolic and hormone signaling pathways represent an essential part of the Yakutian horse adaptive genetic toolkit. Finally, we find evidence for convergent evolution with native human populations and woolly mammoths, suggesting that only a few evolutionary strategies are compatible with survival in extremely cold environments.

To understand patterns of geochemical cycling in pristine versus contaminated groundwater ecosystems, pristine shallow groundwater (FW301) and contaminated groundwater (FW106) samples from the Oak Ridge Integrated Field Research Center (OR-IFRC) were sequenced and compared to each other to determine phylogenetic and metabolic difference between the communities. Proteobacteria (e.g., Burkholderia, Pseudomonas) are the most abundant lineages in the pristine community, though a significant proportion ( >55%) of the community is composed of poorly characterized low abundance (individually <1%) lineages. The phylogenetic diversity of the pristine community contributed to a broader diversity of metabolic networks than the contaminated community. In addition, the pristine community encodes redundant and mostly complete geochemical cycles distributed over multiple lineages and appears capable of a wide range of metabolic activities. In contrast, many geochemical cycles in the contaminated community appear truncated or minimized due to decreased biodiversity and dominance by Rhodanobacter populations capable of surviving the combination of stresses at the site. These results indicate that the pristine site contains more robust and encodes more functional redundancy than the stressed community, which contributes to more efficient nutrient cycling and adaptability than the stressed community.

The domestication of the horse ∼ 5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski's horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the "cost of domestication" hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.

Harnessing the biotechnological potential of the large number of proteins available in sequence databases requires scalable methods for functional characterization. Here we propose a workflow to address this challenge by combining phylogenomic guided DNA synthesis with high-throughput mass spectrometry and apply it to the systematic characterization of GH1 β-glucosidases, a family of enzymes necessary for biomass hydrolysis, an important step in the conversion of lignocellulosic feedstocks to fuels and chemicals. We synthesized and expressed 175 GH1s, selected from over 2000 candidate sequences to cover maximum sequence diversity. These enzymes were functionally characterized over a range of temperatures and pHs using nanostructure-initiator mass spectrometry (NIMS), generating over 10,000 data points. When combined with HPLC-based sugar profiling, we observed GH1 enzymes active over a broad temperature range and toward many different β-linked disaccharides. For some GH1s we also observed activity toward laminarin, a more complex oligosaccharide present as a major component of macroalgae. An area of particular interest was the identification of GH1 enzymes compatible with the ionic liquid 1-ethyl-3-methylimidazolium acetate ([C2mim][OAc]), a next-generation biomass pretreatment technology. We thus searched for GH1 enzymes active at 70 °C and 20% (v/v) [C2mim][OAc] over the course of a 24-h saccharification reaction. Using our unbiased approach, we identified multiple enzymes of different phylogentic origin with such activities. Our approach of characterizing sequence diversity through targeted gene synthesis coupled to high-throughput screening technologies is a broadly applicable paradigm for a wide range of biological problems.

Short non-coding transcripts can be transcribed from distant-acting transcriptional enhancer loci, but the prevalence of such enhancer RNAs (eRNAs) within the transcriptome, and the association of eRNA expression with tissue-specific enhancer activity in vivo remain poorly understood. Here, we investigated the expression dynamics of tissue-specific non-coding RNAs in embryonic mouse tissues via deep RNA sequencing. Overall, approximately 80% of validated in vivo enhancers show tissue-specific RNA expression that correlates with tissue-specific enhancer activity. Globally, we identified thousands of tissue-specifically transcribed non-coding regions (TSTRs) displaying various genomic hallmarks of bona fide enhancers. In transgenic mouse reporter assays, over half of tested TSTRs functioned as enhancers with reproducible activity in the predicted tissue. Together, our results demonstrate that tissue-specific eRNA expression is a common feature of in vivo enhancers, as well as a major source of extragenic transcription, and that eRNA expression signatures can be used to predict tissue-specific enhancers independent of known epigenomic enhancer marks.

Ruminant livestock represent the single largest anthropogenic source of the potent greenhouse gas methane, which is generated by methanogenic archaea residing in ruminant digestive tracts. While differences between individual animals of the same breed in the amount of methane produced have been observed, the basis for this variation remains to be elucidated. To explore the mechanistic basis of this methane production, we measured methane yields from 22 sheep, which revealed that methane yields are a reproducible, quantitative trait. Deep metagenomic and metatranscriptomic sequencing demonstrated a similar abundance of methanogens and methanogenesis pathway genes in high and low methane emitters. However, transcription of methanogenesis pathway genes was substantially increased in sheep with high methane yields. These results identify a discrete set of rumen methanogens whose methanogenesis pathway transcription profiles correlate with methane yields and provide new targets for CH4 mitigation at the levels of microbiota composition and transcriptional regulation., (© 2014 Shi et al.; Published by Cold Spring Harbor Laboratory Press.)

The SMARCA4 (also known as BRG1 in humans) chromatin remodeling factor is critical for establishing lineage-specific chromatin states during early mammalian development. However, the role of SMARCA4 in tissue-specific gene regulation during embryogenesis remains poorly defined. To investigate the genome-wide binding landscape of SMARCA4 in differentiating tissues, we engineered a Smarca4(FLAG) knock-in mouse line. Using ChIP-seq, we identified ∼51,000 SMARCA4-associated regions across six embryonic mouse tissues (forebrain, hindbrain, neural tube, heart, limb, and face) at mid-gestation (E11.5). The majority of these regions was distal from promoters and showed dynamic occupancy, with most distal SMARCA4 sites (73%) confined to a single or limited subset of tissues. To further characterize these regions, we profiled active and repressive histone marks in the same tissues and examined the intersection of informative chromatin states and SMARCA4 binding. This revealed distinct classes of distal SMARCA4-associated elements characterized by activating and repressive chromatin signatures that were associated with tissue-specific up- or down-regulation of gene expression and relevant active/repressed biological pathways. We further demonstrate the predicted active regulatory properties of SMARCA4-associated elements by retrospective analysis of tissue-specific enhancers and direct testing of SMARCA4-bound regions in transgenic mouse assays. Our results indicate a dual active/repressive function of SMARCA4 at distal regulatory sequences in vivo and support its role in tissue-specific gene regulation during embryonic development., (© 2014 Attanasio et al.; Published by Cold Spring Harbor Laboratory Press.)

The canonical genetic code is assumed to be deeply conserved across all domains of life with very few exceptions. By scanning 5.6 trillion base pairs of metagenomic data for stop codon reassignment events, we detected recoding in a substantial fraction of the >1700 environmental samples examined. We observed extensive opal and amber stop codon reassignments in bacteriophages and of opal in bacteria. Our data indicate that bacteriophages can infect hosts with a different genetic code and demonstrate phage-host antagonism based on code differences. The abundance and diversity of genetic codes present in environmental organisms should be considered in the design of engineered organisms with altered genetic codes in order to preclude the exchange of genetic information with naturally occurring species., (Copyright © 2014, American Association for the Advancement of Science.)

Epigenetic information is available from contemporary organisms, but is difficult to track back in evolutionary time. Here, we show that genome-wide epigenetic information can be gathered directly from next-generation sequence reads of DNA isolated from ancient remains. Using the genome sequence data generated from hair shafts of a 4000-yr-old Paleo-Eskimo belonging to the Saqqaq culture, we generate the first ancient nucleosome map coupled with a genome-wide survey of cytosine methylation levels. The validity of both nucleosome map and methylation levels were confirmed by the recovery of the expected signals at promoter regions, exon/intron boundaries, and CTCF sites. The top-scoring nucleosome calls revealed distinct DNA positioning biases, attesting to nucleotide-level accuracy. The ancient methylation levels exhibited high conservation over time, clustering closely with modern hair tissues. Using ancient methylation information, we estimated the age at death of the Saqqaq individual and illustrate how epigenetic information can be used to infer ancient gene expression. Similar epigenetic signatures were found in other fossil material, such as 110,000- to 130,000-yr-old bones, supporting the contention that ancient epigenomic information can be reconstructed from a deep past. Our findings lay the foundation for extracting epigenomic information from ancient samples, allowing shifts in epialleles to be tracked through evolutionary time, as well as providing an original window into modern epigenomics.

Enhancers are distal regulatory elements that can activate tissue-specific gene expression and are abundant throughout mammalian genomes. Although substantial progress has been made toward genome-wide annotation of mammalian enhancers, their temporal activity patterns and global contributions in the context of developmental in vivo processes remain poorly explored. Here we used epigenomic profiling for H3K27ac, a mark of active enhancers, coupled to transgenic mouse assays to examine the genome-wide utilization of enhancers in three different mouse tissues across seven developmental stages. The majority of the ∼90,000 enhancers identified exhibited tightly temporally restricted predicted activity windows and were associated with stage-specific biological functions and regulatory pathways in individual tissues. Comparative genomic analysis revealed that evolutionary conservation of enhancers decreases following midgestation across all tissues examined. The dynamic enhancer activities uncovered in this study illuminate rapid and pervasive temporal in vivo changes in enhancer usage that underlie processes central to development and disease., (Copyright © 2013 Elsevier Inc. All rights reserved.)

The shape of the human face and skull is largely genetically determined. However, the genomic basis of craniofacial morphology is incompletely understood and hypothesized to involve protein-coding genes, as well as gene regulatory sequences. We used a combination of epigenomic profiling, in vivo characterization of candidate enhancer sequences in transgenic mice, and targeted deletion experiments to examine the role of distant-acting enhancers in craniofacial development. We identified complex regulatory landscapes consisting of enhancers that drive spatially complex developmental expression patterns. Analysis of mouse lines in which individual craniofacial enhancers had been deleted revealed significant alterations of craniofacial shape, demonstrating the functional importance of enhancers in defining face and skull morphology. These results demonstrate that enhancers are involved in craniofacial development and suggest that enhancer sequence variation contributes to the diversity of human facial morphology.

Genome sequencing enhances our understanding of the biological world by providing blueprints for the evolutionary and functional diversity that shapes the biosphere. However, microbial genomes that are currently available are of limited phylogenetic breadth, owing to our historical inability to cultivate most microorganisms in the laboratory. We apply single-cell genomics to target and sequence 201 uncultivated archaeal and bacterial cells from nine diverse habitats belonging to 29 major mostly uncharted branches of the tree of life, so-called 'microbial dark matter'. With this additional genomic information, we are able to resolve many intra- and inter-phylum-level relationships and to propose two new superphyla. We uncover unexpected metabolic features that extend our understanding of biology and challenge established boundaries between the three domains of life. These include a novel amino acid use for the opal stop codon, an archaeal-type purine synthesis in Bacteria and complete sigma factors in Archaea similar to those in Bacteria. The single-cell genomes also served to phylogenetically anchor up to 20% of metagenomic reads in some habitats, facilitating organism-level interpretation of ecosystem function. This study greatly expands the genomic representation of the tree of life and provides a systematic step towards a better understanding of biological evolution on our planet.

The mammalian telencephalon plays critical roles in cognition, motor function, and emotion. Though many of the genes required for its development have been identified, the distant-acting regulatory sequences orchestrating their in vivo expression are mostly unknown. Here, we describe a digital atlas of in vivo enhancers active in subregions of the developing telencephalon. We identified more than 4,600 candidate embryonic forebrain enhancers and studied the in vivo activity of 329 of these sequences in transgenic mouse embryos. We generated serial sets of histological brain sections for 145 reproducible forebrain enhancers, resulting in a publicly accessible web-based data collection comprising more than 32,000 sections. We also used epigenomic analysis of human and mouse cortex tissue to directly compare the genome-wide enhancer architecture in these species. These data provide a primary resource for investigating gene regulatory mechanisms of telencephalon development and enable studies of the role of distant-acting enhancers in neurodevelopmental disorders., (Copyright © 2013 Elsevier Inc. All rights reserved.)

The cyanobacterial phylum encompasses oxygenic photosynthetic prokaryotes of a great breadth of morphologies and ecologies; they play key roles in global carbon and nitrogen cycles. The chloroplasts of all photosynthetic eukaryotes can trace their ancestry to cyanobacteria. Cyanobacteria also attract considerable interest as platforms for "green" biotechnology and biofuels. To explore the molecular basis of their different phenotypes and biochemical capabilities, we sequenced the genomes of 54 phylogenetically and phenotypically diverse cyanobacterial strains. Comparison of cyanobacterial genomes reveals the molecular basis for many aspects of cyanobacterial ecophysiological diversity, as well as the convergence of complex morphologies without the acquisition of novel proteins. This phylum-wide study highlights the benefits of diversity-driven genome sequencing, identifying more than 21,000 cyanobacterial proteins with no detectable similarity to known proteins, and foregrounds the diversity of light-harvesting proteins and gene clusters for secondary metabolite biosynthesis. Additionally, our results provide insight into the distribution of genes of cyanobacterial origin in eukaryotic nuclear genomes. Moreover, this study doubles both the amount and the phylogenetic diversity of cyanobacterial genome sequence data. Given the exponentially growing number of sequenced genomes, this diversity-driven study demonstrates the perspective gained by comparing disparate yet related genomes in a phylum-wide context and the insights that are gained from it.

The Deepwater Horizon oil spill in the Gulf of Mexico resulted in a deep-sea hydrocarbon plume that caused a shift in the indigenous microbial community composition with unknown ecological consequences. Early in the spill history, a bloom of uncultured, thus uncharacterized, members of the Oceanospirillales was previously detected, but their role in oil disposition was unknown. Here our aim was to determine the functional role of the Oceanospirillales and other active members of the indigenous microbial community using deep sequencing of community DNA and RNA, as well as single-cell genomics. Shotgun metagenomic and metatranscriptomic sequencing revealed that genes for motility, chemotaxis and aliphatic hydrocarbon degradation were significantly enriched and expressed in the hydrocarbon plume samples compared with uncontaminated seawater collected from plume depth. In contrast, although genes coding for degradation of more recalcitrant compounds, such as benzene, toluene, ethylbenzene, total xylenes and polycyclic aromatic hydrocarbons, were identified in the metagenomes, they were expressed at low levels, or not at all based on analysis of the metatranscriptomes. Isolation and sequencing of two Oceanospirillales single cells revealed that both cells possessed genes coding for n-alkane and cycloalkane degradation. Specifically, the near-complete pathway for cyclohexane oxidation in the Oceanospirillales single cells was elucidated and supported by both metagenome and metatranscriptome data. The draft genome also included genes for chemotaxis, motility and nutrient acquisition strategies that were also identified in the metagenomes and metatranscriptomes. These data point towards a rapid response of members of the Oceanospirillales to aliphatic hydrocarbons in the deep sea.

In the process of clone-based genome sequencing, initial assemblies frequently contain cloning gaps that can be resolved using cloning-independent methods, but the reason for their occurrence is largely unknown. By analyzing 9,328,693 sequencing clones from 393 microbial genomes, we systematically mapped more than 15,000 genes residing in cloning gaps and experimentally showed that their expression products are toxic to the Escherichia coli host. A subset of these toxic sequences was further evaluated through a series of functional assays exploring the mechanisms of their toxicity. Among these genes, our assays revealed novel toxins and restriction enzymes, and new classes of small, non-coding toxic RNAs that reproducibly inhibit E. coli growth. Further analyses also revealed abundant, short, toxic DNA fragments that were predicted to suppress E. coli growth by interacting with the replication initiator DnaA. Our results show that cloning gaps, once considered the result of technical problems, actually serve as a rich source for the discovery of biotechnologically valuable functions, and suggest new modes of antimicrobial interventions.

Development and function of the human heart depend on the dynamic control of tissue-specific gene expression by distant-acting transcriptional enhancers. To generate an accurate genome-wide map of human heart enhancers, we used an epigenomic enhancer discovery approach and identified ∼6,200 candidate enhancer sequences directly from fetal and adult human heart tissue. Consistent with their predicted function, these elements were markedly enriched near genes implicated in heart development, function and disease. To further validate their in vivo enhancer activity, we tested 65 of these human sequences in a transgenic mouse enhancer assay and observed that 43 (66%) drove reproducible reporter gene expression in the heart. These results support the discovery of a genome-wide set of noncoding sequences highly enriched in human heart enhancers that is likely to facilitate downstream studies of the role of enhancers in development and pathological conditions of the heart.

Permafrost contains an estimated 1672 Pg carbon (C), an amount roughly equivalent to the total currently contained within land plants and the atmosphere. This reservoir of C is vulnerable to decomposition as rising global temperatures cause the permafrost to thaw. During thaw, trapped organic matter may become more accessible for microbial degradation and result in greenhouse gas emissions. Despite recent advances in the use of molecular tools to study permafrost microbial communities, their response to thaw remains unclear. Here we use deep metagenomic sequencing to determine the impact of thaw on microbial phylogenetic and functional genes, and relate these data to measurements of methane emissions. Metagenomics, the direct sequencing of DNA from the environment, allows the examination of whole biochemical pathways and associated processes, as opposed to individual pieces of the metabolic puzzle. Our metagenome analyses reveal that during transition from a frozen to a thawed state there are rapid shifts in many microbial, phylogenetic and functional gene abundances and pathways. After one week of incubation at 5 °C, permafrost metagenomes converge to be more similar to each other than while they are frozen. We find that multiple genes involved in cycling of C and nitrogen shift rapidly during thaw. We also construct the first draft genome from a complex soil metagenome, which corresponds to a novel methanogen. Methane previously accumulated in permafrost is released during thaw and subsequently consumed by methanotrophic bacteria. Together these data point towards the importance of rapid cycling of methane and nitrogen in thawing permafrost.

The paucity of enzymes that efficiently deconstruct plant polysaccharides represents a major bottleneck for industrial-scale conversion of cellulosic biomass into biofuels. Cow rumen microbes specialize in degradation of cellulosic plant material, but most members of this complex community resist cultivation. To characterize biomass-degrading genes and genomes, we sequenced and analyzed 268 gigabases of metagenomic DNA from microbes adherent to plant fiber incubated in cow rumen. From these data, we identified 27,755 putative carbohydrate-active genes and expressed 90 candidate proteins, of which 57% were enzymatically active against cellulosic substrates. We also assembled 15 uncultured microbial genomes, which were validated by complementary methods including single-cell genome sequencing. These data sets provide a substantially expanded catalog of genes and genomes participating in the deconstruction of cellulosic biomass.

Modern methods to develop microbe-based biomass conversion processes require a system-level understanding of the microbes involved. Clostridium species have long been recognized as ideal candidates for processes involving biomass conversion and production of various biofuels and other industrial products. To expand the knowledge base for clostridial species relevant to current biofuel production efforts, we have sequenced the genomes of 20 species spanning multiple genera. The majority of species sequenced fall within the class III cellulosome-encoding Clostridium and the class V saccharolytic Thermoanaerobacteraceae. Species were chosen based on representation in the experimental literature as model organisms, ability to degrade cellulosic biomass either by free enzymes or by cellulosomes, ability to rapidly ferment hexose and pentose sugars to ethanol, and ability to ferment synthesis gas to ethanol. The sequenced strains significantly increase the number of noncommensal/nonpathogenic clostridial species and provide a key foundation for future studies of biomass conversion, cellulosome composition, and clostridial systems biology.

Accurate control of tissue-specific gene expression plays a pivotal role in heart development, but few cardiac transcriptional enhancers have thus far been identified. Extreme noncoding-sequence conservation has successfully predicted enhancers that are active in many tissues but has failed to identify substantial numbers of heart-specific enhancers. Here, we used ChIP-Seq with the enhancer-associated protein p300 from mouse embryonic day 11.5 heart tissue to identify over 3,000 candidate heart enhancers genome wide. Compared to enhancers active in other tissues we studied at this time point, most candidate heart enhancers were less deeply conserved in vertebrate evolution. Nevertheless, transgenic mouse assays of 130 candidate regions revealed that most function reproducibly as enhancers active in the heart, irrespective of their degree of evolutionary constraint. These results provide evidence for a large population of poorly conserved heart enhancers and suggest that the evolutionary conservation of embryonic enhancers can vary depending on tissue type.