A bacterial sensor taxonomy across earth ecosystems for machine learning applications

Microbial communities have evolved to colonize all ecosystems of the planet, from the deep sea to the human gut. Microbes survive by sensing, responding, and adapting to immediate environmental cues. This process is driven by signal transduction proteins such as histidine kinases, which use their sensing domains to bind or otherwise detect environmental cues and "transduce" signals to adjust internal processes. We hypothesized that an ecosystem's unique stimuli leave a sensor "fingerprint," able to identify and shed insight on ecosystem conditions. To test this, we collected 20,712 publicly available metagenomes from Host-associated, Environmental, and Engineered ecosystems across the globe. We extracted and clustered the collection's nearly 18M unique sensory domains into 113,712 similar groupings with MMseqs2. We built gradient-boosted decision tree machine learning models and found we could classify the ecosystem type (accuracy: 87%) and predict the levels of different physical parameters (R2 score: 83%) using the sensor cluster abundance as features. Feature importance enables identification of the most predictive sensors to differentiate between ecosystems which can lead to mechanistic interpretations if the sensor domains are well annotated. To demonstrate this, a machine learning model was trained to predict patient's disease state and used to identify domains related to oxygen sensing present in a healthy gut but missing in patients with abnormal conditions. Moreover, since 98.7% of identified sensor domains are uncharacterized, importance ranking can be used to prioritize sensors to determine what ecosystem function they may be sensing. Furthermore, these new predictive sensors can function as targets for novel sensor engineering with applications in biotechnology, ecosystem maintenance, and medicine.IMPORTANCEMicrobes infect, colonize, and proliferate due to their ability to sense and respond quickly to their surroundings. In this research, we extract the sensory proteins from a diverse range of environmental, engineered, and host-associated metagenomes. We trained machine learning classifiers using sensors as features such that it is possible to predict the ecosystem for a metagenome from its sensor profile. We use the optimized model's feature importance to identify the most impactful and predictive sensors in different environments. We next use the sensor profile from human gut metagenomes to classify their disease states and explore which sensors can explain differences between diseases. The sensors most predictive of environmental labels here, most of which correspond to uncharacterized proteins, are a useful starting point for the discovery of important environment signals and the development of possible diagnostic interventions.

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.

Pangenomics reveals alternative environmental lifestyles among chlamydiae

Chlamydiae are highly successful strictly intracellular bacteria associated with diverse eukaryotic hosts. Here we analyzed metagenome-assembled genomes of the "Genomes from Earth's Microbiomes" initiative from diverse environmental samples, which almost double the known phylogenetic diversity of the phylum and facilitate a highly resolved view at the chlamydial pangenome. Chlamydiae are defined by a relatively large core genome indicative of an intracellular lifestyle, and a highly dynamic accessory genome of environmental lineages. We observe chlamydial lineages that encode enzymes of the reductive tricarboxylic acid cycle and for light-driven ATP synthesis. We show a widespread potential for anaerobic energy generation through pyruvate fermentation or the arginine deiminase pathway, and we add lineages capable of molecular hydrogen production. Genome-informed analysis of environmental distribution revealed lineage-specific niches and a high abundance of chlamydiae in some habitats. Together, our data provide an extended perspective of the variability of chlamydial biology and the ecology of this phylum of intracellular microbes.

A genomic catalog of Earth's microbiomes

The reconstruction of bacterial and archaeal genomes from shotgun metagenomes has enabled insights into the ecology and evolution of environmental and host-associated microbiomes. Here we applied this approach to >10,000 metagenomes collected from diverse habitats covering all of Earth's continents and oceans, including metagenomes from human and animal hosts, engineered environments, and natural and agricultural soils, to capture extant microbial, metabolic and functional potential. This comprehensive catalog includes 52,515 metagenome-assembled genomes representing 12,556 novel candidate species-level operational taxonomic units spanning 135 phyla. The catalog expands the known phylogenetic diversity of bacteria and archaea by 44% and is broadly available for streamlined comparative analyses, interactive exploration, metabolic modeling and bulk download. We demonstrate the utility of this collection for understanding secondary-metabolite biosynthetic potential and for resolving thousands of new host linkages to uncultivated viruses. This resource underscores the value of genome-centric approaches for revealing genomic properties of uncultivated microorganisms that affect ecosystem processes.

Insights into the dynamics between viruses and their hosts in a hot spring microbial mat

Our current knowledge of host-virus interactions in biofilms is limited to computational predictions based on laboratory experiments with a small number of cultured bacteria. However, natural biofilms are diverse and chiefly composed of uncultured bacteria and archaea with no viral infection patterns and lifestyle predictions described to date. Herein, we predict the first DNA sequence-based host-virus interactions in a natural biofilm. Using single-cell genomics and metagenomics applied to a hot spring mat of the Cone Pool in Mono County, California, we provide insights into virus-host range, lifestyle and distribution across different mat layers. Thirty-four out of 130 single cells contained at least one viral contig (26%), which, together with the metagenome-assembled genomes, resulted in detection of 59 viruses linked to 34 host species. Analysis of single-cell amplification kinetics revealed a lack of active viral replication on the single-cell level. These findings were further supported by mapping metagenomic reads from different mat layers to the obtained host-virus pairs, which indicated a low copy number of viral genomes compared to their hosts. Lastly, the metagenomic data revealed high layer specificity of viruses, suggesting limited diffusion to other mat layers. Taken together, these observations indicate that in low mobility environments with high microbial abundance, lysogeny is the predominant viral lifestyle, in line with the previously proposed "Piggyback-the-Winner" theory.

IMG/M v.5.0: an integrated data management and comparative analysis system for microbial genomes and microbiomes

The Integrated Microbial Genomes & Microbiomes system v.5.0 (IMG/M: https://img.jgi.doe.gov/m/) contains annotated datasets categorized into: archaea, bacteria, eukarya, plasmids, viruses, genome fragments, metagenomes, cell enrichments, single particle sorts, and metatranscriptomes. Source datasets include those generated by the DOE's Joint Genome Institute (JGI), submitted by external scientists, or collected from public sequence data archives such as NCBI. All submissions are typically processed through the IMG annotation pipeline and then loaded into the IMG data warehouse. IMG's web user interface provides a variety of analytical and visualization tools for comparative analysis of isolate genomes and metagenomes in IMG. IMG/M allows open access to all public genomes in the IMG data warehouse, while its expert review (ER) system (IMG/MER: https://img.jgi.doe.gov/mer/) allows registered users to access their private genomes and to store their private datasets in workspace for sharing and for further analysis. IMG/M data content has grown by 60% since the last report published in the 2017 NAR Database Issue. IMG/M v.5.0 has a new and more powerful genome search feature, new statistical tools, and supports metagenome binning.

Mapping metabolic activity at single cell resolution in intact volcanic fumarole sediment

Interactions among microorganisms and their mineralogical substrates govern the structure, function and emergent properties of microbial communities. These interactions are predicated on spatial relationships, which dictate metabolite exchange and access to key substrates. To quantitatively assess links between spatial relationships and metabolic activity, this study presents a novel approach to map all organisms, the metabolically active subset and associated mineral grains, all while maintaining spatial integrity of an environmental microbiome. We applied this method at an outgassing fumarole of Vanuatu's Marum Crater, one of the largest point sources of several environmentally relevant gaseous compounds, including H2O, CO2 and SO2. With increasing distance from the sediment-air surface and from mineral grain outer boundaries, organism abundance decreased but the proportion of metabolically active organisms often increased. These protected niches may provide more stable conditions that promote consistent metabolic activity of a streamlined community. Conversely, exterior surfaces accumulate more organisms that may cover a wider range of preferred conditions, implying that only a subset of the community will be active under any particular environmental regime. More broadly, the approach presented here allows investigators to see microbial communities ‘as they really are’ and explore determinants of metabolic activity across a range of microbiomes.

Divergent methyl-coenzyme M reductase genes in a deep-subseafloor Archaeoglobi

The methyl-coenzyme M reductase (MCR) complex is a key enzyme in archaeal methane generation and has recently been proposed to also be involved in the oxidation of short-chain hydrocarbons including methane, butane, and potentially propane. The number of archaeal clades encoding the MCR continues to grow, suggesting that this complex was inherited from an ancient ancestor, or has undergone extensive horizontal gene transfer. Expanding the representation of MCR-encoding lineages through metagenomic approaches will help resolve the evolutionary history of this complex. Here, a near-complete Archaeoglobi metagenome-assembled genome (MAG; Ca. Polytropus marinifundus gen. nov. sp. nov.) was recovered from the deep subseafloor along the Juan de Fuca Ridge flank that encodes two divergent McrABG operons similar to those found in Ca. Bathyarchaeota and Ca. Syntrophoarchaeum MAGs. Ca. P. marinifundus is basal to members of the class Archaeoglobi, and encodes the genes for β-oxidation, potentially allowing an alkanotrophic metabolism similar to that proposed for Ca. Syntrophoarchaeum. Ca. P. marinifundus also encodes a respiratory electron transport chain that can potentially utilize nitrate, iron, and sulfur compounds as electron acceptors. Phylogenetic analysis suggests that the Ca. P. marinifundus MCR operons were horizontally transferred, changing our understanding of the evolution and distribution of this complex in the Archaea.

Expanded diversity of microbial groups that shape the dissimilatory sulfur cycle

A critical step in the biogeochemical cycle of sulfur on Earth is microbial sulfate reduction, yet organisms from relatively few lineages have been implicated in this process. Previous studies using functional marker genes have detected abundant, novel dissimilatory sulfite reductases (DsrAB) that could confer the capacity for microbial sulfite/sulfate reduction but were not affiliated with known organisms. Thus, the identity of a significant fraction of sulfate/sulfite-reducing microbes has remained elusive. Here we report the discovery of the capacity for sulfate/sulfite reduction in the genomes of organisms from 13 bacterial and archaeal phyla, thereby more than doubling the number of microbial phyla associated with this process. Eight of the 13 newly identified groups are candidate phyla that lack isolated representatives, a finding only possible given genomes from metagenomes. Organisms from Verrucomicrobia and two candidate phyla, Candidatus Rokubacteria and Candidatus Hydrothermarchaeota, contain some of the earliest evolved dsrAB genes. The capacity for sulfite reduction has been laterally transferred in multiple events within some phyla, and a key gene potentially capable of modulating sulfur metabolism in associated cells has been acquired by putatively symbiotic bacteria. We conclude that current functional predictions based on phylogeny significantly underestimate the extent of sulfate/sulfite reduction across Earth's ecosystems. Understanding the prevalence of this capacity is integral to interpreting the carbon cycle because sulfate reduction is often coupled to turnover of buried organic carbon. Our findings expand the diversity of microbial groups associated with sulfur transformations in the environment and motivate revision of biogeochemical process models based on microbial community composition.

Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea

Standards for sequencing the microbial 'uncultivated majority', namely bacterial and archaeal single-cell genome sequences, and genome sequences from metagenomic datasets, are proposed.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.

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.

Genomic comparisons of a bacterial lineage that inhabits both marine and terrestrial deep subsurface systems

It is generally accepted that diverse, poorly characterized microorganisms reside deep within Earth’s crust. One such lineage of deep subsurface-dwelling bacteria is an uncultivated member of the Firmicutes phylum that can dominate molecular surveys from both marine and continental rock fracture fluids, sometimes forming the sole member of a single-species microbiome. Here, we reconstructed a genome from basalt-hosted fluids of the deep subseafloor along the eastern Juan de Fuca Ridge flank and used a phylogenomic analysis to show that, despite vast differences in geographic origin and habitat, it forms a monophyletic clade with the terrestrial deep subsurface genome of “Candidatus Desulforudis audaxviator” MP104C. While a limited number of differences were observed between the marine genome of “Candidatus Desulfopertinax cowenii” modA32 and its terrestrial relative that may be of potential adaptive importance, here it is revealed that the two are remarkably similar thermophiles possessing the genetic capacity for motility, sporulation, hydrogenotrophy, chemoorganotrophy, dissimilatory sulfate reduction, and the ability to fix inorganic carbon via the Wood-Ljungdahl pathway for chemoautotrophic growth. Our results provide insights into the genetic repertoire within marine and terrestrial members of a bacterial lineage that is widespread in the global deep subsurface biosphere, and provides a natural means to investigate adaptations specific to these two environments.

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.

Nanocalorimetric Characterization of Microbial Activity in Deep Subsurface Oceanic Crustal Fluids

Although fluids within the upper oceanic basaltic crust harbor a substantial fraction of the total prokaryotic cells on Earth, the energy needs of this microbial population are unknown. In this study, a nanocalorimeter (sensitivity down to 1.2 nW ml(-1)) was used to measure the enthalpy of microbially catalyzed reactions as a function of temperature in samples from two distinct crustal fluid aquifers. Microorganisms in unamended, warm (63°C) and geochemically altered anoxic fluids taken from 292 meters sub-basement (msb) near the Juan de Fuca Ridge produced 267.3 mJ of heat over the course of 97 h during a step-wise isothermal scan from 35.5 to 85.0°C. Most of this heat signal likely stems from the germination of thermophilic endospores (6.66 × 10(4) cells ml(-1) FLUID) and their subsequent metabolic activity at temperatures greater than 50°C. The average cellular energy consumption (5.68 pW cell(-1)) reveals the high metabolic potential of a dormant community transported by fluids circulating through the ocean crust. By contrast, samples taken from 293 msb from cooler (3.8°C), relatively unaltered oxic fluids, produced 12.8 mJ of heat over the course of 14 h as temperature ramped from 34.8 to 43.0°C. Corresponding cell-specific energy turnover rates (0.18 pW cell(-1)) were converted to oxygen uptake rates of 24.5 nmol O2 ml(-1) FLUID d(-1), validating previous model predictions of microbial activity in this environment. Given that the investigated fluids are characteristic of expansive areas of the upper oceanic crust, the measured metabolic heat rates can be used to constrain boundaries of habitability and microbial activity in the oceanic crust.

Novel microbial assemblages inhabiting crustal fluids within mid-ocean ridge flank subsurface basalt

Although little is known regarding microbial life within our planet's rock-hosted deep subseafloor biosphere, boreholes drilled through deep ocean sediment and into the underlying basaltic crust provide invaluable windows of access that have been used previously to document the presence of microorganisms within fluids percolating through the deep ocean crust. In this study, the analysis of 1.7 million small subunit ribosomal RNA genes amplified and sequenced from marine sediment, bottom seawater and basalt-hosted deep subseafloor fluids that span multiple years and locations on the Juan de Fuca Ridge flank was used to quantitatively delineate a subseafloor microbiome comprised of distinct bacteria and archaea. Hot, anoxic crustal fluids tapped by newly installed seafloor sampling observatories at boreholes U1362A and U1362B contained abundant bacterial lineages of phylogenetically unique Nitrospirae, Aminicenantes, Calescamantes and Chloroflexi. Although less abundant, the domain Archaea was dominated by unique, uncultivated lineages of marine benthic group E, the Terrestrial Hot Spring Crenarchaeotic Group, the Bathyarchaeota and relatives of cultivated, sulfate-reducing Archaeoglobi. Consistent with recent geochemical measurements and bioenergetic predictions, the potential importance of methane cycling and sulfate reduction were imprinted within the basalt-hosted deep subseafloor crustal fluid microbial community. This unique window of access to the deep ocean subsurface basement reveals a microbial landscape that exhibits previously undetected spatial heterogeneity.

Activity and phylogenetic diversity of sulfate-reducing microorganisms in low-temperature subsurface fluids within the upper oceanic crust

The basaltic ocean crust is the largest aquifer system on Earth, yet the rates of biological activity in this environment are unknown. Low-temperature (<100°C) fluid samples were investigated from two borehole observatories in the Juan de Fuca Ridge (JFR) flank, representing a range of upper oceanic basement thermal and geochemical properties. Microbial sulfate reduction rates (SRR) were measured in laboratory incubations with ³⁵S-sulfate over a range of temperatures and the identity of the corresponding sulfate-reducing microorganisms (SRM) was studied by analyzing the sequence diversity of the functional marker dissimilatory (bi)sulfite reductase (dsrAB) gene. We found that microbial sulfate reduction was limited by the decreasing availability of organic electron donors in higher temperature, more altered fluids. Thermodynamic calculations indicate energetic constraints for metabolism, which together with relatively higher cell-specific SRR reveal increased maintenance requirements, consistent with novel species-level dsrAB phylotypes of thermophilic SRM. Our estimates suggest that microbially-mediated sulfate reduction may account for the removal of organic matter in fluids within the upper oceanic crust and underscore the potential quantitative impact of microbial processes in deep subsurface marine crustal fluids on marine and global biogeochemical carbon cycling.

Phylogenetic diversity of microorganisms in subseafloor crustal fluids from Holes 1025C and 1026B along the Juan de Fuca Ridge flank

To expand investigations into the phylogenetic diversity of microorganisms inhabiting the subseafloor biosphere, basalt-hosted crustal fluids were sampled from Circulation Obviation Retrofit Kits (CORKs) affixed to Holes 1025C and 1026B along the Juan de Fuca Ridge (JdFR) flank using a clean fluid pumping system. These boreholes penetrate the crustal aquifer of young ocean crust (1.24 and 3.51 million years old, respectively), but differ with respect to borehole depth and temperature at the sediment-basement interface (147 m and 39°C vs. 295 m and 64°C, respectively). Cloning and sequencing of PCR-amplified small subunit ribosomal RNA genes revealed that fluids retrieved from Hole 1025C were dominated by relatives of the genus Desulfobulbus of the Deltaproteobacteria (56% of clones) and Candidatus Desulforudis of the Firmicutes (17%). Fluids sampled from Hole 1026B also contained plausible deep subseafloor inhabitants amongst the most abundant clone lineages; however, both geochemical analysis and microbial community structure reveal the borehole to be compromised by bottom seawater intrusion. Regardless, this study provides independent support for previous observations seeking to identify phylogenetic groups of microorganisms common to the deep ocean crustal biosphere, and extends previous observations by identifying additional lineages that may be prevalent in this unique environment.