Universally distributed single-copy genes indicate a constant rate of horizontal transfer

Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity

The potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored the potential for rapid evolution in Festuca rubra, a long-lived grass with extensive clonal growth dominating in alpine grasslands. We used a field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from five independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored the genetic profiles of established seedlings after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important than increasing temperature or the interaction of the two climatic factors in selecting pressure. This can indicate negative genetic interaction in response to the combined effects or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions, likely affecting the function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity), suggesting these changes may have a consequence on plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.

Discovery and characterization of the α-amylases cDNAs from Enchytraeus albidus shed light on the evolution of "Enchytraeus-Eisenia type" Amy homologs in Annelida

Although enchytraeids have gained popularity in scientific research, fundamental questions regarding their feeding ecology and biology remain largely unexplored. This study investigates α-amylases, major digestive enzymes responsible for hydrolyzing starch and similar polysaccharides into sugars, in Enchytraeus albidus. Genetic data related to α-amylases is currently lacking for the family Enchytraeidae but also for the entire Annelida. To detect and identify coding sequences of the expressed α-amylase genes in COI-monohaplotype culture (PL-A strain) of E. albidus, we used classical "gene fishing" and transcriptomic approaches. We also compared coding sequence variants of α-amylase retrieved from transcriptomic data related to freeze-tolerant strains. Our results reveal that E. albidus possesses two distinct α-amylase genes (Amy I and Amy II) that are homologs to earthworm Eisenia fetida Ef-Amy genes. Different strains of E. albidus possess distinctive alleles of α-amylases with unique SNP patterns specific to a particular strain. Unlike Amy II, Amy I seems to be a highly polymorphic and multicopy gene. The domain architecture of the putative Amy proteins was found the same as for classical animal α-amylases with ABC-domains. A characteristic feature of Amy II is the lack of GHGA motif in the flexible loop region, similarly to many insect amylases. We identified "Enchytraeus-Eisenia type" α-amylase homologs in other clitellates and polychaetes, indicating the ancestral origin of Amy I/II proteins in Annelida. This study provides the first insight into the endogenous non-proteolytic digestive enzyme genes in potworms, discusses the evolution of Amy α-amylases in Annelida, and explores phylogenetic implications.

Generation of accurate, expandable phylogenomic trees with uDance

Phylogenetic trees provide a framework for organizing evolutionary histories across the tree of life and aid downstream comparative analyses such as metagenomic identification. Methods that rely on single-marker genes such as 16S rRNA have produced trees of limited accuracy with hundreds of thousands of organisms, whereas methods that use genome-wide data are not scalable to large numbers of genomes. We introduce updating trees using divide-and-conquer (uDance), a method that enables updatable genome-wide inference using a divide-and-conquer strategy that refines different parts of the tree independently and can build off of existing trees, with high accuracy and scalability. With uDance, we infer a species tree of roughly 200,000 genomes using 387 marker genes, totaling 42.5 billion amino acid residues.

A metagenomic view of novel microbial and metabolic diversity found within the deep terrestrial biosphere at DeMMO: A microbial observatory in South Dakota, USA

The deep terrestrial subsurface is a large and diverse microbial habitat and vast repository of biomass. However, in relation to its size and physical heterogeneity we have limited understanding of taxonomic and metabolic diversity in this realm. Here we present a detailed metagenomic analysis of samples from the Deep Mine Microbial Observatory (DeMMO) spanning depths from the surface to 1.5 km into the crust. From eight geochemically and spatially distinct fluid samples we reconstructed ~600 partial to near-complete metagenome-assembled genomes (MAGs), representing 50 distinct phyla and including 18 candidate phyla. These novel clades include members of the candidate phyla radiation, two new MAGs from OLB16, a phylum originally identified in DeMMO fluids and for which only one other MAG is currently available, and new MAGs from the Eisenbacteria, Omnitrophota, and Edwardsbacteria. We find that microbes spanning this expansive phylogenetic diversity and physical subsurface space gain a competitive edge by maintaining a wide variety of functional pathways, are often capable of numerous dissimilatory energy metabolisms and poised to take advantage of nutrients as they become available in isolated fracture fluids. Our results support and expand on emerging themes of tight nutrient cycling and genomic plasticity in deep subsurface biosphere taxa.

VBCG: 20 validated bacterial core genes for phylogenomic analysis with high fidelity and resolution

BACKGROUND

Phylogenomic analysis has become an inseparable part of studies of bacterial diversity and evolution, and many different bacterial core genes have been collated and used for phylogenomic tree reconstruction. However, these genes have been selected based on their presence and single-copy ratio in all bacterial genomes, leaving out the gene's 'phylogenetic fidelity' unexamined.

RESULTS

From 30,522 complete genomes covering 11,262 species, we examined 148 bacterial core genes that have been previously used for phylogenomic analysis. In addition to the gene presence and single-copy rations, we evaluated the gene's phylogenetic fidelity by comparing each gene's phylogeny with its corresponding 16S rRNA gene tree. Out of the 148 bacterial genes, 20 validated bacterial core genes (VBCG) were selected as the core gene set with the highest bacterial phylogenetic fidelity. Compared to the larger gene set, the 20-gene core set resulted in more species having all genes present and fewer species with missing data, thereby enhancing the accuracy of phylogenomic analysis. Using Escherichia coli strains as examples of prominent bacterial foodborne pathogens, we demonstrated that the 20 VBCG produced phylogenies with higher fidelity and resolution at species and strain levels while 16S rRNA gene tree alone could not.

CONCLUSION

The 20 validated core gene set improves the fidelity and speed of phylogenomic analysis. Among other uses, this tool improves our ability to explore the evolution, typing and tracking of bacterial strains, such as human pathogens. We have developed a Python pipeline and a desktop graphic app (available on GitHub) for users to perform phylogenomic analysis with high fidelity and resolution. Video Abstract.

Horizontal Gene Transfer and Drug Resistance Involving Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) acquires drug resistance at a rate comparable to that of bacterial pathogens that replicate much faster and have a higher mutation rate. One explanation for this rapid acquisition of drug resistance in Mtb is that drug resistance may evolve in other fast-replicating mycobacteria and then be transferred to Mtb through horizontal gene transfer (HGT). This paper aims to address three questions. First, does HGT occur between Mtb and other mycobacterial species? Second, what genes after HGT tend to survive in the recipient genome? Third, does HGT contribute to antibiotic resistance in Mtb? I present a conceptual framework for detecting HGT and analyze 39 ribosomal protein genes, 23S and 16S ribosomal RNA genes, as well as several genes targeted by antibiotics against Mtb, from 43 genomes representing all major groups within Mycobacterium. I also included mgtC and the insertion sequence IS6110 that were previously reported to be involved in HGT. The insertion sequence IS6110 shows clearly that the Mtb complex participates in HGT. However, the horizontal transferability of genes depends on gene function, as was previously hypothesized. HGT is not observed in functionally important genes such as ribosomal protein genes, rRNA genes, and other genes chosen as drug targets. This pattern can be explained by differential selection against functionally important and unimportant genes after HGT. Functionally unimportant genes such as IS6110 are not strongly selected against, so HGT events involving such genes are visible. For functionally important genes, a horizontally transferred diverged homologue from a different species may not work as well as the native counterpart, so the HGT event involving such genes is strongly selected against and eliminated, rendering them invisible to us. In short, while HGT involving the Mtb complex occurs, antibiotic resistance in the Mtb complex arose from mutations in those drug-targeted genes within the Mtb complex and was not gained through HGT.

OrthoVenn3: an integrated platform for exploring and visualizing orthologous data across genomes

Advancements in comparative genomics research have led to a growing interest in studying species evolution and genetic diversity. To facilitate this research, OrthoVenn3 has been developed as a powerful, web-based tool that enables users to efficiently identify and annotate orthologous clusters and infer phylogenetic relationships across a range of species. The latest upgrade of OrthoVenn includes several important new features, including enhanced orthologous cluster identification accuracy, improved visualization capabilities for numerous sets of data, and wrapped phylogenetic analysis. Furthermore, OrthoVenn3 now provides gene family contraction and expansion analysis to support researchers better understanding the evolutionary history of gene families, as well as collinearity analysis to detect conserved and variable genomic structures. With its intuitive user interface and robust functionality, OrthoVenn3 is a valuable resource for comparative genomics research. The tool is freely accessible at https://orthovenn3.bioinfotoolkits.net.

Metagenomics revealing molecular profiles of microbial community structure and metabolic capacity in Bamucuo lake, Tibet

Microorganisms play critical ecological roles in the global biogeochemical cycles. However, extensive information on the microbial communities in Qinghai-Tibet Plateau (QTP), which is the highest plateau in the world, is still lacking, particularly in high elevation locations above 4500 m. Here, we performed a survey of th e soil and water microbial communities in Bamucuo Lake, Tibet, by using shotgun metagenomic methods. In the soil and water samples, we reconstructed 75 almost complete metagenomic assembly genomes, and 74 of the metagenomic assembly genomes from the water sample represented novel species. Proteobacteria and Actinobacteria were found to be the dominant bacterial phyla, while Euryarchaeota was the dominant archaeal phylum. The largest virus, Pandoravirus salinus, was found in the soil microbial community. We concluded that the microorganisms in Bamucuo Lake are most likely to fix carbon mainly through the 3-hydroxypropionic bi-cycle pathway. This study, for the first time, characterized the microbial community composition and metabolic capacity in QTP high-elevation locations with 4555 m, confirming that QTP is a vast and valuable resource pool, in which many microorganisms can be used to develop new bioactive substances and new antibiotics to which pathogenic microorganisms have not yet developed resistance.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

A revisit to universal single-copy genes in bacterial genomes

Universal single-copy genes (USCGs) are widely used for species classification and taxonomic profiling. Despite many studies on USCGs, our understanding of USCGs in bacterial genomes might be out of date, especially how different the USCGs are in different studies, how well a set of USCGs can distinguish two bacterial species, whether USCGs can separate different strains of a bacterial species, to name a few. To fill the void, we studied USCGs in the most updated complete bacterial genomes. We showed that different USCG sets are quite different while coming from highly similar functional categories. We also found that although USCGs occur once in almost all bacterial genomes, each USCG does occur multiple times in certain genomes. We demonstrated that USCGs are reliable markers to distinguish different species while they cannot distinguish different strains of most bacterial species. Our study sheds new light on the usage and limitations of USCGs, which will facilitate their applications in evolutionary, phylogenomic, and metagenomic studies.

Comprehensive phylogenomics of Methylobacterium reveals four evolutionary distinct groups and underappreciated phyllosphere diversity

Methylobacterium is a group of methylotrophic microbes associated with soil, fresh water, and particularly the phyllosphere, the aerial part of plants that has been well-studied in terms of physiology but whose evolutionary history and taxonomy are unclear. Recent work has suggested that Methylobacterium is much more diverse than thought previously, questioning its status as an ecologically and phylogenetically coherent taxonomic genus. However, taxonomic and evolutionary studies of Methylobacterium have mostly been restricted to model species, often isolated from habitats other than the phyllosphere, and have yet to utilize comprehensive phylogenomic methods to examine gene trees, gene content, or synteny. By analyzing 189 Methylobacterium genomes from a wide range of habitats, including the phyllosphere, we inferred a robust phylogenetic tree while explicitly accounting for the impact of horizontal gene transfer. We showed that Methylobacterium contains four evolutionarily distinct groups of bacteria (namely A, B, C, D), characterized by different genome size, GC content, gene content and genome architecture, revealing the dynamic nature of Methylobacterium genomes. In addition to recovering 59 described species, we identified 45 candidate species, mostly phyllosphere-associated, stressing the significance of plants as a reservoir of Methylobacterium diversity. We inferred an ancient transition from a free-living lifestyle to association with plant roots in Methylobacteriaceae ancestor, followed by phyllosphere association of three of the major groups (A, B, D), whose early branching in Methylobacterium history has been heavily obscured by HGT. Together, our work lays the foundations for a thorough redefinition of Methylobacterium taxonomy, beginning with the abandonment of Methylorubrum.

An expanded reference map of the human gut microbiome reveals hundreds of previously unknown species

The gut is the richest ecosystem of microbes in the human body and has great influence on our health. Despite many efforts, the set of microbes inhabiting this environment is not fully known, limiting our ability to identify microbial content and to research it. In this work, we combine new microbial metagenomic assembled genomes from 51,052 samples, with previously published genomes to produce a curated set of 241,118 genomes. Based on this set, we procure a new and improved human gut microbiome reference set of 3594 high quality species genomes, which successfully matches 83.65% validation samples’ reads. This improved reference set contains 310 novel species, including one that exists in 19% of validation samples. Overall, this study provides a gut microbial genome reference set that can serve as a valuable resource for further research.

Here, Leviatan et al. produce 241,118 genome assemblies to produce a new human gut microbiome reference set of 3,594 species genomes, of which 310 represent previously undescribed species, making the catalog a valuable resource for further research.

Paraburkholderia bengalensis sp. nov. isolated from roots of Oryza sativa, IR64

Paraburkholderia bengalensis sp. nov. strain IR64_4_BI was isolated from rice roots cultivated in Madhyamgram field station of Bose Institute, West Bengal, India. IR64_4_BI is a Gram-negative, motile, nitrate-reducing, nitrogen-fixing bacterium. Whole-cell fatty acid analyses of IR64_4_BI show C16:0, summed feature 8 (comprising C18:1ω7c and/or C18:1 ω 6c) and summed feature 3(C16:1 w7c/C16:1 w6c or C16:1 ω 7c/C16:1 ω 6c) were the predominant fatty acids. 16S rRNA phylogeny showed that it was most similar to P. phymatum STM815^T (98.5% identity), P. terrae KMY02^T (98.44% identity) and P. hospita LMG 20598T (98.32% identity). The Average Nucleotide Identity-BLAST (ANIb) of P. bengalensis IR64_4_BI with P. hospita DSM 17164^T, P. terrae DSM 17804^T, P. phymatum STM815^T and P. hospita LMG 20598T was 83.11, 83.52, 84.5 and 83.12% respectively. Comparison of genome sequence of IR64_4_BI with other species of Paraburkholderia using the Multi-locus species tree software show that P. bengalensis IR64_4_BI is a novel species. The ability of P. bengalensis IR64_4_BI to survive on nitrogen-free medium under microaerophilic conditions and the abundance of nitrogen metabolism-related genes makes this strain a potential candidate for developing a nitrogen-fixing system in rice. Based on genotypic, phenotypic and chemotaxonomic studies, we propose that IR64_4_BI (= MTCC 13051 = JCM 34777) is a new species of Paraburkholderia which has been assigned as Paraburkholderia bengalensis sp.nov.

Landscape of mobile genetic elements and their antibiotic resistance cargo in prokaryotic genomes

Prokaryotic Mobile Genetic Elements (MGEs) such as transposons, integrons, phages and plasmids, play important roles in prokaryotic evolution and in the dispersal of cargo functions like antibiotic resistance. However, each of these MGE types is usually annotated and analysed individually, hampering a global understanding of phylogenetic and environmental patterns of MGE dispersal. We thus developed a computational framework that captures diverse MGE types, their cargos and MGE-mediated horizontal transfer events, using recombinases as ubiquitous MGE marker genes and pangenome information for MGE boundary estimation. Applied to ∼84k genomes with habitat annotation, we mapped 2.8 million MGE-specific recombinases to six operational MGE types, which together contain on average 13% of all the genes in a genome. Transposable elements (TEs) dominated across all taxa (∼1.7 million occurrences), outnumbering phages and phage-like elements (<0.4 million). We recorded numerous MGE-mediated horizontal transfer events across diverse phyla and habitats involving all MGE types, disentangled and quantified the extent of hitchhiking of TEs (17%) and integrons (63%) with other MGE categories, and established TEs as dominant carriers of antibiotic resistance genes. We integrated all these findings into a resource (proMGE.embl.de), which should facilitate future studies on the large mobile part of genomes and its horizontal dispersal.

Growth promotion and antibiotic induced metabolic shifts in the chicken gut microbiome

Antimicrobial growth promoters (AGP) have played a decisive role in animal agriculture for over half a century. Despite mounting concerns about antimicrobial resistance and demand for antibiotic alternatives, a thorough understanding of how these compounds drive performance is missing. Here we investigate the functional footprint of microbial communities in the cecum of chickens fed four distinct AGP. We find relatively few taxa, metabolic or antimicrobial resistance genes similarly altered across treatments, with those changes often driven by the abundances of core microbiome members. Constraints-based modeling of 25 core bacterial genera associated increased performance with fewer metabolite demands for microbial growth, pointing to altered nitrogen utilization as a potential mechanism of narasin, the AGP with the largest performance increase in our study. Untargeted metabolomics of narasin treated birds aligned with model predictions, suggesting that the core cecum microbiome might be targeted for enhanced performance via its contribution to host-microbiota metabolic crosstalk.

This study compares the functional profiles of the cecal microbiome among chickens fed four different antimicrobial growth promoters. Chickens receiving narasin exhibited the largest performance increase via apparent nitrogen recycling by the core cecal microbiome.

An estimate of the deepest branches of the tree of life from ancient vertically-evolving genes

Core gene phylogenies provide a window into early evolution, but different gene sets and analytical methods have yielded substantially different views of the tree of life. Trees inferred from a small set of universal core genes have typically supported a long branch separating the archaeal and bacterial domains. By contrast, recent analyses of a broader set of non-ribosomal genes have suggested that Archaea may be less divergent from Bacteria, and that estimates of inter-domain distance are inflated due to accelerated evolution of ribosomal proteins along the inter-domain branch. Resolving this debate is key to determining the diversity of the archaeal and bacterial domains, the shape of the tree of life, and our understanding of the early course of cellular evolution. Here, we investigate the evolutionary history of the marker genes key to the debate. We show that estimates of a reduced Archaea-Bacteria (AB) branch length result from inter-domain gene transfers and hidden paralogy in the expanded marker gene set. By contrast, analysis of a broad range of manually curated marker gene datasets from an evenly sampled set of 700 Archaea and Bacteria reveals that current methods likely underestimate the AB branch length due to substitutional saturation and poor model fit; that the best-performing phylogenetic markers tend to support longer inter-domain branch lengths; and that the AB branch lengths of ribosomal and non-ribosomal marker genes are statistically indistinguishable. Furthermore, our phylogeny inferred from the 27 highest-ranked marker genes recovers a clade of DPANN at the base of the Archaea and places the bacterial Candidate Phyla Radiation (CPR) within Bacteria as the sister group to the Chloroflexota.

Novel nitrite reductase domain structure suggests a chimeric denitrification repertoire in the phylum Chloroflexi

Denitrification plays a central role in the global nitrogen cycle, reducing and removing nitrogen from marine and terrestrial ecosystems. The flux of nitrogen species through this pathway has a widespread impact, affecting ecological carrying capacity, agriculture, and climate. Nitrite reductase (Nir) and nitric oxide reductase (NOR) are the two central enzymes in this pathway. Here we present a previously unreported Nir domain architecture in members of phylum Chloroflexi. Phylogenetic analyses of protein domains within Nir indicate that an ancestral horizontal transfer and fusion event produced this chimeric domain architecture. We also identify an expanded genomic diversity of a rarely reported NOR subtype, eNOR. Together, these results suggest a greater diversity of denitrification enzyme arrangements exist than have been previously reported.

Assembling a Reference Phylogenomic Tree of Bacteria and Archaea by Summarizing Many Gene Phylogenies

Phylogenomics is the inference of phylogenetic trees based on multiple marker genes sampled in the genomes of interest. An important challenge in phylogenomics is the potential incongruence among the evolutionary histories of individual genes, which can be widespread in microorganisms due to the prevalence of horizontal gene transfer. This protocol introduces the procedures for building a phylogenetic tree of a large number of microbial genomes using a broad sampling of marker genes that are representative of whole-genome evolution. The protocol highlights the use of a gene tree summary method, which can effectively reconstruct the species tree while accounting for the topological conflicts among individual gene trees. The pipeline described in this protocol is scalable to tens of thousands of genomes while retaining high accuracy. We discussed multiple software tools, libraries, and scripts to enable convenient adoption of the protocol. The protocol is suitable for microbiology and microbiome studies based on public genomes and metagenomic data.

Phylogenetic Signal, Congruence, and Uncertainty across Bacteria and Archaea

Reconstruction of the Tree of Life is a central goal in biology. Although numerous novel phyla of bacteria and archaea have recently been discovered, inconsistent phylogenetic relationships are routinely reported, and many inter-phylum and inter-domain evolutionary relationships remain unclear. Here, we benchmark different marker genes often used in constructing multidomain phylogenetic trees of bacteria and archaea and present a set of marker genes that perform best for multidomain trees constructed from concatenated alignments. We use recently-developed Tree Certainty metrics to assess the confidence of our results and to obviate the complications of traditional bootstrap-based metrics. Given the vastly disparate number of genomes available for different phyla of bacteria and archaea, we also assessed the impact of taxon sampling on multidomain tree construction. Our results demonstrate that biases between the representation of different taxonomic groups can dramatically impact the topology of resulting trees. Inspection of our highest-quality tree supports the division of most bacteria into Terrabacteria and Gracilicutes, with Thermatogota and Synergistota branching earlier from these superphyla. This tree also supports the inclusion of the Patescibacteria within the Terrabacteria as a sister group to the Chloroflexota instead of as a basal-branching lineage. For the Archaea, our tree supports three monophyletic lineages (DPANN, Euryarchaeota, and TACK/Asgard), although we note the basal placement of the DPANN may still represent an artifact caused by biased sequence composition. Our findings provide a robust and standardized framework for multidomain phylogenetic reconstruction that can be used to evaluate inter-phylum relationships and assess uncertainty in conflicting topologies of the Tree of Life.

Genome-resolved metagenomics provides insights into the functional complexity of microbial mats in Blue Holes, Shark Bay

The present study describes for the first time the community composition and functional potential of the microbial mats found in the supratidal, gypsum-rich, and hypersaline region of Blue Holes, Shark Bay. This was achieved via high throughput metagenomic sequencing of total mat community DNA and complementary analyses using hyperspectral confocal microscopy. Mat communities were dominated by Proteobacteria (29%), followed by Bacteroidetes/Chlorobi Group (11%), and Planctomycetes (10%). These mats were found to also harbor a diverse community of potentially novel microorganisms including members from the DPANN, Asgard archaea, and Candidate Phyla Radiation, with highest diversity found in the lower regions (∼14-20 mm depth) of the mat. In addition to pathways for major metabolic cycles, a range of putative rhodopsins with previously uncharacterized motifs and functions were identified along with heliorhodopsins and putative schizorhodopsins. Critical microbial interactions were also inferred, and from 117 medium-to-high quality metagenome-assembled genomes (MAGs), viral defense mechanisms (CRISPR, BREX, and DISARM), elemental transport, osmoprotection, heavy metal and UV resistance were also detected. These analyses have provided a greater understanding of these distinct mat systems in Shark Bay, including key insights into adaptive responses and proposing that photoheterotrophy may be an important lifestyle in Blue Holes.

Phylogenetic systematics of Butyrivibrio and Pseudobutyrivibrio genomes illustrate vast taxonomic diversity, open genomes and an abundance of carbohydrate-active enzyme family isoforms

Butyrivibrio and Pseudobutyrivibrio dominate in anaerobic gastrointestinal microbiomes, particularly the rumen, where they play a key role in harvesting dietary energy. Within these genera, five rumen species have been classified (Butyrivibrio fibrisolvens, Butyrivibrio hungatei, Butyrivibrio proteoclasticus, Pseudobutyrivibrio ruminis and Pseudobutyrivibrio xylanivorans) and more recently an additional Butyrivibrio sp. group was added. Given the recent increase in available genomes, we re-investigated the phylogenetic systematics and evolution of Butyrivibrio and Pseudobutyrivibrio. Across 71 genomes, we show using 16S rDNA and 40 gene marker phylogenetic trees that the current six species designations (P. ruminis, P. xylanivorans, B. fibrisolvens, Butyrivibrio sp., B. hungatei and B. proteclasticus) are found. However, pangenome analysis showed vast genomic variation and a high abundance of accessory genes (91.50–99.34 %), compared with core genes (0.66–8.50 %), within these six taxonomic groups, suggesting incorrectly assigned taxonomy. Subsequent pangenome accessory genomes under varying core gene cut-offs (%) and average nucleotide identity (ANI) analysis suggest the existence of 42 species within 32 genera. Pangenome analysis of those that still group within B. fibrisolvens, B. hungatei and P. ruminis, based on revised ANI phylogeny, also showed possession of very open genomes, illustrating the diversity that exists even within these groups. All strains of both Butyrivibrio and Pseudobutyrivibrio also shared a broad range of clusters of orthologous genes (COGs) (870), indicating recent evolution from a common ancestor. We also demonstrate that the carbohydrate-active enzymes (CAZymes) predominantly belong to glycosyl hydrolase (GH)2, 3, 5, 13 and 43, with numerous within family isoforms apparent, likely facilitating metabolic plasticity and resilience under dietary perturbations. This study provides a major advancement in our functional and evolutionary understanding of these important anaerobic bacteria.

Evidence for Selection in the Abundant Accessory Gene Content of a Prokaryote Pangenome

A pangenome is the complete set of genes (core and accessory) present in a phylogenetic clade. We hypothesize that a pangenome's accessory gene content is structured and maintained by selection. To test this hypothesis, we interrogated the genomes of 40 Pseudomonas species for statistically significant coincident (i.e., co-occurring/avoiding) gene patterns. We found that 86.7% of common accessory genes are involved in ≥1 coincident relationship. Further, genes that co-occur and/or avoid each other-but are not vertically inherited-are more likely to share functional categories, are more likely to be simultaneously transcribed, and are more likely to produce interacting proteins, than would be expected by chance. These results are not due to coincident genes being adjacent to one another on the chromosome. Together, these findings suggest that the accessory genome is structured into sets of genes that function together within a given strain. Given the similarity of the Pseudomonas pangenome with open pangenomes of other prokaryotic species, we speculate that these results are generalizable.

Reductive evolution and unique predatory mode in the CPR bacterium Vampirococcus lugosii

The Candidate Phyla Radiation (CPR) constitutes a large group of mostly uncultured bacterial lineages with small cell sizes and limited biosynthetic capabilities. They are thought to be symbionts of other organisms, but the nature of this symbiosis has been ascertained only for cultured Saccharibacteria, which are epibiotic parasites of other bacteria. Here, we study the biology and the genome of Vampirococcus lugosii, which becomes the first described species of Vampirococcus, a genus of epibiotic bacteria morphologically identified decades ago. Vampirococcus belongs to the CPR phylum Absconditabacteria. It feeds on anoxygenic photosynthetic gammaproteobacteria, fully absorbing their cytoplasmic content. The cells divide epibiotically, forming multicellular stalks whose apical cells can reach new hosts. The genome is small (1.3 Mbp) and highly reduced in biosynthetic metabolism genes, but is enriched in genes possibly related to a fibrous cell surface likely involved in interactions with the host. Gene loss has been continuous during the evolution of Absconditabacteria, and generally most CPR bacteria, but this has been compensated by gene acquisition by horizontal gene transfer and de novo evolution. Our findings support parasitism as a widespread lifestyle of CPR bacteria, which probably contribute to the control of bacterial populations in diverse ecosystems.

Microbial ecology of the newly discovered serpentinite-hosted Old City hydrothermal field (southwest Indian ridge)

Lost City (mid-Atlantic ridge) is a unique oceanic hydrothermal field where carbonate-brucite chimneys are colonized by a single phylotype of archaeal Methanosarcinales, as well as sulfur- and methane-metabolizing bacteria. So far, only one submarine analog of Lost City has been characterized, the Prony Bay hydrothermal field (New Caledonia), which nonetheless shows more microbiological similarities with ecosystems associated with continental ophiolites. This study presents the microbial ecology of the ‘Lost City’-type Old City hydrothermal field, recently discovered along the southwest Indian ridge. Five carbonate-brucite chimneys were sampled and subjected to mineralogical and geochemical analyses, microimaging, as well as 16S rRNA-encoding gene and metagenomic sequencing. Dominant taxa and metabolisms vary between chimneys, in conjunction with the predicted redox state, while potential formate- and CO-metabolizing microorganisms as well as sulfur-metabolizing bacteria are always abundant. We hypothesize that the variable environmental conditions resulting from the slow and diffuse hydrothermal fluid discharge that currently characterizes Old City could lead to different microbial populations between chimneys that utilize CO and formate differently as carbon or electron sources. Old City discovery and this first description of its microbial ecology opens up attractive perspectives for understanding environmental factors shaping communities and metabolisms in oceanic serpentinite-hosted ecosystems.

Culture-enriched metagenomic sequencing enables in-depth profiling of the cystic fibrosis lung microbiota

Amplicon sequencing (for example, of the 16S rRNA gene) identifies the presence and relative abundance of microbial community members. However, metagenomic sequencing is needed to identify the genetic content and functional potential of a community. Metagenomics is challenging in samples dominated by host DNA, such as those from the skin, tissue and respiratory tract. Here, we combine advances in amplicon and metagenomic sequencing with culture-enriched molecular profiling to study the human microbiota. Using the cystic fibrosis lung as an example, we cultured an average of 82.13% of the operational taxonomic units representing 99.3% of the relative abundance identified in direct sequencing of sputum samples; importantly, culture enrichment identified 63.3% more operational taxonomic units than direct sequencing. We developed the PLate Coverage Algorithm (PLCA) to determine a representative subset of culture plates on which to conduct culture-enriched metagenomics, resulting in the recovery of greater taxonomic diversity-including of low-abundance taxa-with better metagenome-assembled genomes, longer contigs and better functional annotations when compared to culture-independent methods. The PLCA is also applied as a proof of principle to a previously published gut microbiota dataset. Culture-enriched molecular profiling can be used to better understand the role of the human microbiota in health and disease.

Phylogenomics of 10,575 genomes reveals evolutionary proximity between domains Bacteria and Archaea

Rapid growth of genome data provides opportunities for updating microbial evolutionary relationships, but this is challenged by the discordant evolution of individual genes. Here we build a reference phylogeny of 10,575 evenly-sampled bacterial and archaeal genomes, based on a comprehensive set of 381 markers, using multiple strategies. Our trees indicate remarkably closer evolutionary proximity between Archaea and Bacteria than previous estimates that were limited to fewer "core" genes, such as the ribosomal proteins. The robustness of the results was tested with respect to several variables, including taxon and site sampling, amino acid substitution heterogeneity and saturation, non-vertical evolution, and the impact of exclusion of candidate phyla radiation (CPR) taxa. Our results provide an updated view of domain-level relationships.

Optimizing sequencing protocols for leaderboard metagenomics by combining long and short reads

As metagenomic studies move to increasing numbers of samples, communities like the human gut may benefit more from the assembly of abundant microbes in many samples, rather than the exhaustive assembly of fewer samples. We term this approach leaderboard metagenome sequencing. To explore protocol optimization for leaderboard metagenomics in real samples, we introduce a benchmark of library prep and sequencing using internal references generated by synthetic long-read technology, allowing us to evaluate high-throughput library preparation methods against gold-standard reference genomes derived from the samples themselves. We introduce a low-cost protocol for high-throughput library preparation and sequencing.

AzeR, a transcriptional regulator that responds to azelaic acid in Pseudomonas nitroreducens

Azelaic acid is a dicarboxylic acid that has recently been shown to play a role in plant-bacteria signalling and also occurs naturally in several cereals. Several bacteria have been reported to be able to utilize azelaic acid as a unique source of carbon and energy, including Pseudomonas nitroreducens. In this study, we utilize P. nitroreducens as a model organism to study bacterial degradation of and response to azelaic acid. We report genetic evidence of azelaic acid degradation and the identification of a transcriptional regulator that responds to azelaic acid in P. nitroreducens DSM 9128. Three mutants possessing transposons in genes of an acyl-CoA ligase, an acyl-CoA dehydrogenase and an isocitrate lyase display a deficient ability in growing in azelaic acid. Studies on transcriptional regulation of these genes resulted in the identification of an IclR family repressor that we designated as AzeR, which specifically responds to azelaic acid. A bioinformatics survey reveals that AzeR is confined to a few proteobacterial genera that are likely to be able to degrade and utilize azelaic acid as the sole source of carbon and energy.

Persistent Core Populations Shape the Microbiome Throughout the Water Column in the North Pacific Subtropical Gyre

Marine microbial communities are responsible for many important ecosystem processes in the oceans. Their variability across time and depths is well recognized, but mostly at a coarse-grained taxonomic resolution. To gain a deeper perspective on ecological patterns of bacterioplankton diversity in the North Pacific Subtropical Gyre, we characterized bacterioplankton communities throughout the water column at a fine-grained taxonomic level with a focus on temporally persistent (core) populations. Considerable intra-clade microdiversity was evident in virtually every microbial clade examined. While some of the most abundant populations comprised only a small fraction of the intra-clade microdiversity, they formed a temporally persistent core within a more diverse array of less abundant ephemeral populations. The depth-stratified population structure within many phylogenetically disparate clades suggested that ecotypic variation was the rule among most planktonic bacterial and archaeal lineages. Our results suggested that the abundant, persistent core populations comprised the bulk of the biomass within any given clade. As such, we postulate that these core populations are largely responsible for microbially driven ecosystem processes, and so represent ideal targets for elucidating key microbial processes in the open-ocean water column.

Metabolic versatility in a modern lineage of cyanobacteria from terrestrial hot springs

The extent of oxygenated environments on the early Earth was much lower than today, and cyanobacteria were critical players in Earth's shift from widespread anoxia to oxygenated surface environments. Extant cyanobacteria that aggregate into cones, tufts and ridges are used to understand the long record of photosynthesis and microbe-mineral interactions during times when oxygen was much lower, i.e., the Archean and the Proterozoic. To better understand the metabolic versatility and physiological properties of these organisms, we examined publicly available genomes of cyanobacteria from modern terrestrial hydrothermal systems and a newly sequenced genome of a cyanobacterium isolated from conical and ridged microbialites that grow in occasionally sulfidic hydrothermal springs in Yellowstone National Park, USA. Phylogenomic analyses reveal that cyanobacteria from globally distributed terrestrial and shallow marine hydrothermal systems form a monophyletic clade within the Cyanobacteria phylum. Comparative genomics of this clade reveals the genetic capacity for oxygenic photosynthesis that uses photosystems I and II, and anoxygenic photosynthesis that uses a putative sulfide quinone reductase to oxidize sulfide and bypass photosystem II. Surprisingly large proportions of the newly sequenced genome from Yellowstone National Park are also dedicated to secondary metabolite production (15.1-15.6%), of which ∼6% can be attributed to antibiotic production and resistance genes. All this may be advantageous to benthic, mat-forming photosynthesizers that have to compete for light and nutrients in sporadically or permanently sulfidic environments, and may have also improved the tolerance of ancient counterparts of these cyanobacteria to sulfidic conditions in benthic communities that colonized the coastal margins in the Archean and the Proterozoic.

Unexpected host dependency of Antarctic Nanohaloarchaeota

We demonstrate that Candidatus Nanohaloarchaeum antarcticus requires Halorubrum lacusprofundi for growth, illustrating that Nanohaloarchaeota require a host rather than being free living as previously proposed. Developing the means of cultivating Nanohaloarchaeota in the laboratory provides the capacity to advance understanding of how archaea interact and the factors that control their symbiotic relationship (e.g. mutualism, commensalism, antagonism). Our findings amplify the view that Antarctic lakes are a treasure trove for the discovery of microbes with previously unknown properties.

In hypersaline environments, Nanohaloarchaeota (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, Nanohaloarchaeota [DPANN] superphylum) are thought to be free-living microorganisms. We report cultivation of 2 strains of Antarctic Nanohaloarchaeota and show that they require the haloarchaeon Halorubrum lacusprofundi for growth. By performing growth using enrichments and fluorescence-activated cell sorting, we demonstrated successful cultivation of Candidatus Nanohaloarchaeum antarcticus, purification of Ca. Nha. antarcticus away from other species, and growth and verification of Ca. Nha. antarcticus with Hrr. lacusprofundi; these findings are analogous to those required for fulfilling Koch’s postulates. We use fluorescent in situ hybridization and transmission electron microscopy to assess cell structures and interactions; metagenomics to characterize enrichment taxa, generate metagenome assembled genomes, and interrogate Antarctic communities; and proteomics to assess metabolic pathways and speculate about the roles of certain proteins. Metagenome analysis indicates the presence of a single species, which is endemic to Antarctic hypersaline systems that support the growth of haloarchaea. The presence of unusually large proteins predicted to function in attachment and invasion of hosts plus the absence of key biosynthetic pathways (e.g., lipids) in metagenome assembled genomes of globally distributed Nanohaloarchaeota indicate that all members of the lineage have evolved as symbionts. Our work expands the range of archaeal symbiotic lifestyles and provides a genetically tractable model system for advancing understanding of the factors controlling microbial symbiotic relationships.

Horizontal gene transfer in human-associated microorganisms inferred by phylogenetic reconstruction and reconciliation

Horizontal gene transfer (HGT) is widespread in the evolution of prokaryotes, especially those associated with the human body. Here, we implemented large-scale gene-species phylogenetic tree reconstructions and reconciliations to identify putative HGT-derived genes in the reference genomes of microbiota isolated from six major human body sites by the NIH Human Microbiome Project. Comparisons with a control group representing microbial genomes from diverse natural environments indicated that HGT activity increased significantly in the genomes of human microbiota, which is confirmatory of previous findings. Roughly, more than half of total genes in the genomes of human-associated microbiota were transferred (donated or received) by HGT. Up to 60% of the detected HGTs occurred either prior to the colonization of the human body or involved bacteria residing in different body sites. The latter could suggest 'genetic crosstalk' and movement of bacterial genes within the human body via hitherto poorly understood mechanisms. We also observed that HGT activity increased significantly among closely-related microorganisms and especially when they were united by physical proximity, suggesting that the 'phylogenetic effect' can significantly boost HGT activity. Finally, we identified several core and widespread genes least influenced by HGT that could become useful markers for building robust 'trees of life' and address several outstanding technical challenges to improve the phylogeny-based genome-wide HGT detection method for future applications.

Functional shifts in microbial mats recapitulate early Earth metabolic transitions

Phototrophic microbial mats dominated terrestrial ecosystems for billions of years, largely causing, through cyanobacterial oxygenic photosynthesis, but also undergoing, the Great Oxidation Event approximately 2.5 billion years ago. Taking a space-for-time approach based on the universality of core metabolic pathways expressed at ecosystem level, we studied gene content and co-occurrence networks in high-diversity metagenomes from spatially close microbial mats along a steep redox gradient. The observed functional shifts suggest that anoxygenic photosynthesis was present but not predominant under early Precambrian conditions, being accompanied by other autotrophic processes. Our data also suggest that, in contrast to general assumptions, anoxygenic photosynthesis largely expanded in parallel with the subsequent evolution of oxygenic photosynthesis and aerobic respiration. Finally, our observations might represent space-for-time evidence that the Wood-Ljungdahl carbon fixation pathway dominated phototrophic mats in early ecosystems, whereas the Calvin cycle probably evolved from pre-existing variants before becoming the dominant contemporary form of carbon fixation.

Widespread Antibiotic, Biocide, and Metal Resistance in Microbial Communities Inhabiting a Municipal Waste Environment and Anthropogenically Impacted River

Landfill leachate is a persistent contamination threat for terrestrial waters. Microbial metabolism in landfills transforms contaminants and contributes to greenhouse gas emissions. A better understanding of landfill-associated microbial communities will inform bioremediation of solid waste environments and improve pathogen monitoring. We leveraged shotgun metagenomics to investigate the microbial communities of the Riverton City dump and the adjoining Duhaney River near Kingston City, Jamaica. We identified no overlap between the microbial communities inhabiting the Riverton City dump leachate and the Duhaney River. Both communities are predicted to degrade aromatic compounds, which are ubiquitous environmental pollutants. Adversely, microbes in both environments are predicted to withstand widely used antibiotics, antiseptics, and metal contamination. The absence of evidence for microbial transfer from the leachate to the river is encouraging; however, the Duhaney River contained several organisms with predicted pathogenic lifestyles, indicating that the river represents a human health risk regardless of impact from the dump.

The Riverton City dump is Jamaica’s largest solid waste disposal site, but it lacks engineered protection for leachate containment and treatment. Shotgun metagenomics was used to survey the microbial communities in the Riverton City dump leachate and in surface waters of the Duhaney River, an urban waterway abutting the dump. The community within the leachate pond was taxonomically distinct from that found in the surface waters of the Duhaney River. Higher microbial diversity was observed within the dump leachate, with members of the Bacteroidetes, Firmicutes, Gammaproteobacteria, Deltaproteobacteria, and Tenericutes being the most abundant, while the river community was dominated by Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria. The microbial communities exhibit similar functional potential profiles, including chemoorganoheterotrophy as the dominant metabolism, and the potential to degrade aromatic compounds. From reconstruction of metagenome-assembled genomes (MAGs), organisms within both environments are predicted to survive in the presence of multiple antibiotics, antiseptics, biocides, and metals. Strong virulence potential coincided with the most diverse multiple resistance profiles in 1 of 5 leachate MAGs and 5 of 33 river MAGs. Unexpectedly, the microbial resistance profiles were more varied and widespread in the river populations, where we had expected the chemical composition of the leachate to select and enrich for resistance characteristics. This study provides valuable insights into the total functional potential of a landfill leachate microbial community and identifies possible human health hazards within the Duhaney River and Riverton City dump, urban environments with the potential to impact human populations.

IMPORTANCE Landfill leachate is a persistent contamination threat for terrestrial waters. Microbial metabolism in landfills transforms contaminants and contributes to greenhouse gas emissions. A better understanding of landfill-associated microbial communities will inform bioremediation of solid waste environments and improve pathogen monitoring. We leveraged shotgun metagenomics to investigate the microbial communities of the Riverton City dump and the adjoining Duhaney River near Kingston City, Jamaica. We identified no overlap between the microbial communities inhabiting the Riverton City dump leachate and the Duhaney River. Both communities are predicted to degrade aromatic compounds, which are ubiquitous environmental pollutants. Adversely, microbes in both environments are predicted to withstand widely used antibiotics, antiseptics, and metal contamination. The absence of evidence for microbial transfer from the leachate to the river is encouraging; however, the Duhaney River contained several organisms with predicted pathogenic lifestyles, indicating that the river represents a human health risk regardless of impact from the dump.

Enhanced antibacterial effect of the novel T4-like bacteriophage KARL-1 in combination with antibiotics against multi-drug resistant Acinetobacter baumannii

The continuing rise of infections caused by multi-drug resistant bacteria has led to a renewed interest in bacteriophage therapy. Here we characterize phage vB_AbaM-KARL-1 with lytic activity against multi-drug resistant clinical isolates of Acinetobacter baumannii (AB). Besides genomic and phenotypic phage analysis, the objective of our study was to investigate the antibacterial outcome when the phage acts in concert with distinct antibiotics. KARL-1 belongs to the family of Myoviridae and is able to lyse 8 of 20 (40%) tested clinical isolates. Its double-stranded DNA genome consists of 166,560 bp encoding for 253 open reading frames. Genome wide comparison suggests that KARL-1 is a novel species within the subfamily Tevenvirinae, sharing 77% nucleotide identity (coverage 58%) with phage ZZ1. The antibacterial efficacy at various multiplicities of infection (MOI) was monitored either alone or in combination with meropenem, ciprofloxacin, and colistin. A complete clearance of liquid cultures was achieved with KARL-1 at an MOI of 10^-1 and meropenem (>128 mg/l). KARL-1 was still effective at an MOI of 10^-7, but antibacterial activity was significantly augmented with meropenem. While ciprofloxacin did generally not support phage activity, the application of KARL-1 at an MOI of 10^-7 and therapeutic doses of colistin significantly elevated bacterial suppression. Hence, KARL-1 represents a novel candidate for use against multi-drug resistant AB and the therapeutic outcome may be positively influenced by the addition of traditional antibiotics.

Geological and Geochemical Constraints on the Origin and Evolution of Life

The traditional tree of life from molecular biology with last universal common ancestor (LUCA) branching into bacteria and archaea (though fuzzy) is likely formally valid enough to be a basis for discussion of geological processes on the early Earth. Biologists infer likely properties of nodal organisms within the tree and, hence, the environment they inhabited. Geologists both vet tenuous trees and putative origin of life scenarios for geological and ecological reasonability and conversely infer geological information from trees. The latter approach is valuable as geologists have only weakly constrained the time when the Earth became habitable and the later time when life actually existed to the long interval between ∼4.5 and ∼3.85 Ga where no intact surface rocks are known. With regard to vetting, origin and early evolution hypotheses from molecular biology have recently centered on serpentinite settings in marine and alternatively land settings that are exposed to ultraviolet sunlight. The existence of these niches on the Hadean Earth is virtually certain. With regard to inferring geological environment from genomics, nodes on the tree of life can arise from true bottlenecks implied by the marine serpentinite origin scenario and by asteroid impact. Innovation of a very useful trait through a threshold allows the successful organism to quickly become very abundant and later root a large clade. The origin of life itself, that is, the initial Darwinian ancestor, the bacterial and archaeal roots as free-living cellular organisms that independently escaped hydrothermal chimneys above marine serpentinite or alternatively from shallow pore-water environments on land, the Selabacteria root with anoxygenic photosynthesis, and the Terrabacteria root colonizing land are attractive examples that predate the geological record. Conversely, geological reasoning presents likely events for appraisal by biologists. Asteroid impacts may have produced bottlenecks by decimating life. Thermophile roots of bacteria and archaea as well as a thermophile LUCA are attractive.

Genomic Description of 'Candidatus Abyssubacteria,' a Novel Subsurface Lineage Within the Candidate Phylum Hydrogenedentes

The subsurface biosphere is a massive repository of fixed carbon, harboring approximately 90% of Earth’s microbial biomass. These microbial communities drive transformations central to Earth’s biogeochemical cycles. However, there is still much we do not understand about how complex subterranean microbial communities survive and how they interact with these cycles. Recent metagenomic investigation of deeply circulating terrestrial subsurface fluids revealed the presence of several novel lineages of bacteria. In one particular example, phylogenomic analyses do not converge on any one previously identified taxon; here we describe the first full genomic sequences of a new bacterial lineage within the candidate phylum Hydrogenedentes, ‘Candidatus Abyssubacteria.’ A global survey revealed that members of this proposed lineage are widely distributed in both marine and terrestrial subsurface environments, but their physiological and ecological roles have remained unexplored. Two high quality metagenome assembled genomes (SURF_5: 97%, 4%; SURF_17: 91% and 4% completeness and contamination, respectively) were reconstructed from fluids collected 1.5 kilometers below surface in the former Homestake gold mine—now the Sanford Underground Research Facility (SURF)—in Lead, South Dakota, United States. Metabolic reconstruction suggests versatile metabolic capability, including possible nitrogen reduction, sulfite oxidation, sulfate reduction and homoacetogenesis. This first glimpse into the metabolic capabilities of these cosmopolitan bacteria suggests that they are involved in key geochemical processes, including sulfur, nitrogen, and carbon cycling, and that they are adapted to survival in the dark, often anoxic, subsurface biosphere.

Nitrogen-fixing populations of Planctomycetes and Proteobacteria are abundant in surface ocean metagenomes

Nitrogen fixation in the surface ocean impacts global marine nitrogen bioavailability and thus microbial primary productivity. Until now, cyanobacterial populations have been viewed as the main suppliers of bioavailable nitrogen in this habitat. Although PCR amplicon surveys targeting the nitrogenase reductase gene have revealed the existence of diverse non-cyanobacterial diazotrophic populations, subsequent quantitative PCR surveys suggest that they generally occur in low abundance. Here, we use state-of-the-art metagenomic assembly and binning strategies to recover nearly one thousand non-redundant microbial population genomes from the TARA Oceans metagenomes. Among these, we provide the first genomic evidence for non-cyanobacterial diazotrophs inhabiting surface waters of the open ocean, which correspond to lineages within the Proteobacteria and, most strikingly, the Planctomycetes. Members of the latter phylum are prevalent in aquatic systems, but have never been linked to nitrogen fixation previously. Moreover, using genome-wide quantitative read recruitment, we demonstrate that the discovered diazotrophs were not only widespread but also remarkably abundant (up to 0.3% of metagenomic reads for a single population) in both the Pacific Ocean and the Atlantic Ocean northwest. Our results extend decades of PCR-based gene surveys, and substantiate the importance of heterotrophic bacteria in the fixation of nitrogen in the surface ocean.

UBCG: Up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction

Genome-based phylogeny plays a central role in the future taxonomy and phylogenetics of Bacteria and Archaea by replacing 16S rRNA gene phylogeny. The concatenated core gene alignments are frequently used for such a purpose. The bacterial core genes are defined as single-copy, homologous genes that are present in most of the known bacterial species. There have been several studies describing such a gene set, but the number of species considered was rather small. Here we present the up-to-date bacterial core gene set, named UBCG, and software suites to accommodate necessary steps to generate and evaluate phylogenetic trees. The method was successfully used to infer phylogenomic relationship of Escherichia and related taxa and can be used for the set of genomes at any taxonomic ranks of Bacteria. The UBCG pipeline and file viewer are freely available at https://www.ezbiocloud.net/tools/ubcg and https://www.ezbiocloud.net/tools/ubcg_viewer , respectively.

Whole-Genome Characterization of Bacillus cereus Associated with Specific Disease Manifestations

Bacillus cereus remains an important cause of infections, particularly in immunocompromised hosts. While typically associated with enteric infections, disease manifestations can be quite diverse and include skin infections, bacteremia, pneumonia, and meningitis. Whether there are any genetic correlates of bacterial strains with particular clinical manifestations remains unknown. To address this gap in understanding, we undertook whole-genome analysis of B. cereus strains isolated from patients with a range of disease manifestations, including noninvasive colonizing disease, superficial skin infections, and invasive bacteremia. Interestingly, strains involved in skin infection tended to form a distinct genetic cluster compared to isolates associated with invasive disease. Other disease manifestations, despite not being exclusively clustered, nonetheless had unique genetic features. The unique features associated with the specific types of infections ranged from traditional virulence determinants to metabolic pathways and gene regulators. These data represent the largest genetic analysis to date of pathogenic B. cereus isolates with associated clinical parameters.

Subspecies in the global human gut microbiome

Population genomics of prokaryotes has been studied in depth in only a small number of primarily pathogenic bacteria, as genome sequences of isolates of diverse origin are lacking for most species. Here, we conducted a large-scale survey of population structure in prevalent human gut microbial species, sampled from their natural environment, with a culture-independent metagenomic approach. We examined the variation landscape of 71 species in 2,144 human fecal metagenomes and found that in 44 of these, accounting for 72% of the total assigned microbial abundance, single-nucleotide variation clearly indicates the existence of sub-populations (here termed subspecies). A single subspecies (per species) usually dominates within each host, as expected from ecological theory. At the global scale, geographic distributions of subspecies differ between phyla, with Firmicutes subspecies being significantly more geographically restricted. To investigate the functional significance of the delineated subspecies, we identified genes that consistently distinguish them in a manner that is independent of reference genomes. We further associated these subspecies-specific genes with properties of the microbial community and the host. For example, two of the three Eubacterium rectale subspecies consistently harbor an accessory pro-inflammatory flagellum operon that is associated with lower gut community diversity, higher host BMI, and higher blood fasting insulin levels. Using an additional 676 human oral samples, we further demonstrate the existence of niche specialized subspecies in the different parts of the oral cavity. Taken together, we provide evidence for subspecies in the majority of abundant gut prokaryotes, leading to a better functional and ecological understanding of the human gut microbiome in conjunction with its host.

High reactivity of deep biota under anthropogenic CO2 injection into basalt

Basalts are recognized as one of the major habitats on Earth, harboring diverse and active microbial populations. Inconsistently, this living component is rarely considered in engineering operations carried out in these environments. This includes carbon capture and storage (CCS) technologies that seek to offset anthropogenic CO₂ emissions into the atmosphere by burying this greenhouse gas in the subsurface. Here, we show that deep ecosystems respond quickly to field operations associated with CO₂ injections based on a microbiological survey of a basaltic CCS site. Acidic CO₂-charged groundwater results in a marked decrease (by ~ 2.5–4) in microbial richness despite observable blooms of lithoautotrophic iron-oxidizing Betaproteobacteria and degraders of aromatic compounds, which hence impact the aquifer redox state and the carbon fate. Host-basalt dissolution releases nutrients and energy sources, which sustain the growth of autotrophic and heterotrophic species whose activities may have consequences on mineral storage.

The impacts of carbon capture and storage (CCS) on subsurface microorganisms are poorly understood. Here, the authors show that deep ecosystems respond quickly to CO2 injections and that the environmental consequences of their metabolic activities need to be properly assessed for sustainable CCS in basalt.

A Massively Parallel Sequence Similarity Search for Metagenomic Sequencing Data

Sequence similarity searches have been widely used in the analyses of metagenomic sequencing data. Finding homologous sequences in a reference database enables the estimation of taxonomic and functional characteristics of each query sequence. Because current metagenomic sequencing data consist of a large number of nucleotide sequences, the time required for sequence similarity searches account for a large proportion of the total time. This time-consuming step makes it difficult to perform large-scale analyses. To analyze large-scale metagenomic data, such as those found in the human oral microbiome, we developed GHOST-MP (Genome-wide HOmology Search Tool on Massively Parallel system), a parallel sequence similarity search tool for massively parallel computing systems. This tool uses a fast search algorithm based on suffix arrays of query and database sequences and a hierarchical parallel search to accelerate the large-scale sequence similarity search of metagenomic sequencing data. The parallel computing efficiency and the search speed of this tool were evaluated. GHOST-MP was shown to be scalable over 10,000 CPU (Central Processing Unit) cores, and achieved over 80-fold acceleration compared with mpiBLAST using the same computational resources. We applied this tool to human oral metagenomic data, and the results indicate that the oral cavity, the oral vestibule, and plaque have different characteristics based on the functional gene category.

DESMAN: a new tool for de novo extraction of strains from metagenomes

We introduce DESMAN for De novo Extraction of Strains from Metagenomes. Large multi-sample metagenomes are being generated but strain variation results in fragmentary co-assemblies. Current algorithms can bin contigs into metagenome-assembled genomes but are unable to resolve strain-level variation. DESMAN identifies variants in core genes and uses co-occurrence across samples to link variants into haplotypes and abundance profiles. These are then searched for against non-core genes to determine the accessory genome of each strain. We validated DESMAN on a complex 50-species 210-genome 96-sample synthetic mock data set and then applied it to the Tara Oceans microbiome.

Evidence for H2 consumption by uncultured Desulfobacterales in coastal sediments

Molecular hydrogen (H₂ ) is the key intermediate in the anaerobic degradation of organic matter. Its removal by H₂ -oxidizing microorganisms is essential to keep anaerobic degradation energetically favourable. Sulfate-reducing microorganisms (SRM) are known as the main H₂ scavengers in anoxic marine sediments. Although the community of marine SRM has been extensively studied, those consuming H₂ in situ are completely unknown. We combined metagenomics, PCR-based clone libraries, single-amplified genomes (SAGs) and metatranscriptomics to identify potentially H₂ -consuming SRM in anoxic coastal sediments. The vast majority of SRM-related H₂ ase sequences were assigned to group 1b and 1c [NiFe]-H₂ ases of the deltaproteobacterial order Desulfobacterales. Surprisingly, the same sequence types were similarly highly expressed in spring and summer, suggesting that these are stable and integral members of the H₂ -consuming community. Notably, one sequence cluster from the SRM group 1 consistently accounted for around half of all [NiFe]-H₂ ase transcripts. Using SAGs, we could link this cluster with the 16S rRNA genes of the uncultured Sva0081-group of the family Desulfobacteraceae. Sequencing of 16S rRNA gene amplicons and H₂ ase gene libraries suggested consistently high in situ abundance of the Sva0081 group also in other marine sediments. Together with other Desulfobacterales these likely are important H₂ -scavengers in marine sediments.

Energy and carbon metabolisms in a deep terrestrial subsurface fluid microbial community

The terrestrial deep subsurface is a huge repository of microbial biomass, but in relation to its size and physical heterogeneity, few sites have been investigated in detail. Here, we applied a culture-independent metagenomic approach to characterize the microbial community composition in deep (1500 meters below surface) terrestrial fluids. Samples were collected from a former gold mine in Lead, South Dakota, USA, now Sanford Underground Research Facility (SURF). We reconstructed 74 genomes from metagenomes (MAGs), enabling the identification of common metabolic pathways. Sulfate and nitrate/nitrite reduction were the most common putative energy metabolisms. Complete pathways for autotrophic carbon fixation were found in more than half of the MAGs, with the reductive acetyl-CoA pathway by far the most common. Nearly 40% (29 of 74) of the recovered MAGs belong to bacterial phyla without any cultivated members-microbial dark matter. Three of our MAGs constitute two novel phyla previously only identified in 16 S rRNA gene surveys. The uniqueness of this data set-its physical depth in the terrestrial subsurface, the relative abundance and completeness of microbial dark matter genomes and the overall diversity of this physically deep, dark, community-make it an invaluable addition to our knowledge of deep subsurface microbial ecology.

Draft Genome Sequence and Assembly of a Lysobacter enzymogenes Strain with Biological Control Activity against Root Knot Nematodes

Lysobacter enzymogenes strain B25, an isolate from an agricultural field, acts as a biological control agent against root knot nematodes in tomato plants. B25 also controls several fungal diseases and promotes plant growth under abiotic stress. We hereby report on the draft genome sequence and assembly of B25.

Metagenome sequencing and 98 microbial genomes from Juan de Fuca Ridge flank subsurface fluids

The global deep subsurface biosphere is one of the largest reservoirs for microbial life on our planet. This study takes advantage of new sampling technologies and couples them with improvements to DNA sequencing and associated informatics tools to reconstruct the genomes of uncultivated Bacteria and Archaea from fluids collected deep within the Juan de Fuca Ridge subseafloor. Here, we generated two metagenomes from borehole observatories located 311 meters apart and, using binning tools, retrieved 98 genomes from metagenomes (GFMs). Of the GFMs, 31 were estimated to be >90% complete, while an additional 17 were >70% complete. Phylogenomic analysis revealed 53 bacterial and 45 archaeal GFMs, of which nearly all were distantly related to known cultivated isolates. In the GFMs, abundant Bacteria included Chloroflexi, Nitrospirae, Acetothermia (OP1), EM3, Aminicenantes (OP8), Gammaproteobacteria, and Deltaproteobacteria, while abundant Archaea included Archaeoglobi, Bathyarchaeota (MCG), and Marine Benthic Group E (MBG-E). These data are the first GFMs reconstructed from the deep basaltic subseafloor biosphere, and provide a dataset available for further interrogation.

Divergent functional isoforms drive niche specialisation for nutrient acquisition and use in rumen microbiome

Many microbes in complex competitive environments share genes for acquiring and utilising nutrients, questioning whether niche specialisation exists and if so, how it is maintained. We investigated the genomic signatures of niche specialisation in the rumen microbiome, a highly competitive, anaerobic environment, with limited nutrient availability determined by the biomass consumed by the host. We generated individual metagenomic libraries from 14 cows fed an ad libitum diet of grass silage and calculated functional isoform diversity for each microbial gene identified. The animal replicates were used to calculate confidence intervals to test for differences in diversity of functional isoforms between microbes that may drive niche specialisation. We identified 153 genes with significant differences in functional isoform diversity between the two most abundant bacterial genera in the rumen (Prevotella and Clostridium). We found Prevotella possesses a more diverse range of isoforms capable of degrading hemicellulose, whereas Clostridium for cellulose. Furthermore, significant differences were observed in key metabolic processes indicating that isoform diversity plays an important role in maintaining their niche specialisation. The methods presented represent a novel approach for untangling complex interactions between microorganisms in natural environments and have resulted in an expanded catalogue of gene targets central to rumen cellulosic biomass degradation.