Genomic resolution of a cold subsurface aquifer community provides metabolic insights for novel microbes adapted to high CO2 concentrations

Microbiological and Geochemical Survey of CO2-Dominated Mofette and Mineral Waters of the Cheb Basin, Czech Republic

The Cheb Basin (NW Bohemia, Czech Republic) is a shallow, neogene intracontinental basin. It is a non-volcanic region which features frequent earthquake swarms and large-scale diffuse degassing of mantle-derived CO₂ at the surface that occurs in the form of CO₂-rich mineral springs and wet and dry mofettes. So far, the influence of CO₂ degassing onto the microbial communities has been studied for soil environments, but not for aquatic systems. We hypothesized, that deep-trenching CO₂ conduits interconnect the subsurface with the surface. This admixture of deep thermal fluids should be reflected in geochemical parameters and in the microbial community compositions. In the present study four mineral water springs and two wet mofettes were investigated through an interdisciplinary survey. The waters were acidic and differed in terms of organic carbon and anion/cation concentrations. Element geochemical and isotope analyses of fluid components were used to verify the origin of the fluids. Prokaryotic communities were characterized through quantitative PCR and Illumina 16S rRNA gene sequencing. Putative chemolithotrophic, anaerobic and microaerophilic organisms connected to sulfur (e.g., Sulfuricurvum, Sulfurimonas) and iron (e.g., Gallionella, Sideroxydans) cycling shaped the core community. Additionally, CO₂-influenced waters form an ecosystem containing many taxa that are usually found in marine or terrestrial subsurface ecosystems. Multivariate statistics highlighted the influence of environmental parameters such as pH, Fe²⁺ concentration and conductivity on species distribution. The hydrochemical and microbiological survey introduces a new perspective on mofettes. Our results support that mofettes are either analogs or rather windows into the deep biosphere and furthermore enable access to deeply buried paleo-sediments.

Geochemical Influence on Microbial Communities at CO2-Leakage Analog Sites

Microorganisms influence the chemical and physical properties of subsurface environments and thus represent an important control on the fate and environmental impact of CO₂ that leaks into aquifers from deep storage reservoirs. How leakage will influence microbial populations over long time scales is largely unknown. This study uses natural analog sites to investigate the long-term impact of CO₂ leakage from underground storage sites on subsurface biogeochemistry. We considered two sites with elevated CO₂ levels (sample groups I and II) and one control site with low CO₂ content (group III). Samples from sites with elevated CO₂ had pH ranging from 6.2 to 4.5 and samples from the low-CO₂ control group had pH ranging from 7.3 to 6.2. Solute concentrations were relatively low for samples from the control group and group I but high for samples from group II, reflecting varying degrees of water-rock interaction. Microbial communities were analyzed through clone library and MiSeq sequencing. Each 16S rRNA analysis identified various bacteria, methane-producing archaea, and ammonia-oxidizing archaea. Both bacterial and archaeal diversities were low in groundwater with high CO₂ content and community compositions between the groups were also clearly different. In group II samples, sequences classified in groups capable of methanogenesis, metal reduction, and nitrate reduction had higher relative abundance in samples with relative high methane, iron, and manganese concentrations and low nitrate levels. Sequences close to Comamonadaceae were abundant in group I, while the taxa related to methanogens, Nitrospirae, and Anaerolineaceae were predominant in group II. Our findings provide insight into subsurface biogeochemical reactions that influence the carbon budget of the system including carbon fixation, carbon trapping, and CO₂ conversion to methane. The results also suggest that monitoring groundwater microbial community can be a potential tool for tracking CO₂ leakage from geologic storage sites.

Members of the Candidate Phyla Radiation are functionally differentiated by carbon- and nitrogen-cycling capabilities

BACKGROUND

The Candidate Phyla Radiation (CPR) is a recently described expansion of the tree of life that represents more than 15% of all bacterial diversity and potentially contains over 70 different phyla. Despite this broad phylogenetic variation, these microorganisms appear to feature little functional diversity, with members generally characterized as obligate fermenters. Additionally, much of the data describing CPR phyla has been generated from a limited number of environments, constraining our knowledge of their functional roles and biogeographical distribution. To expand our understanding of subsurface CPR microorganisms, we sampled four separate groundwater wells over 2 years across three Ohio counties.

RESULTS

Samples were analyzed using 16S rRNA gene amplicon and shotgun metagenomic sequencing. Amplicon results indicated that CPR members comprised between 2 and 20% of the microbial communities with relative abundances stable through time in Athens and Greene samples but dynamic in Licking groundwater. Shotgun metagenomic analyses generated 71 putative CPR genomes, representing roughly 32 known phyla and 2 putative new phyla, Candidatus Brownbacteria and Candidatus Hugbacteria. While these genomes largely mirrored metabolic characteristics of known CPR members, some features were previously uncharacterized. For instance, nitrite reductase, encoded by nirK, was found in four of our Parcubacteria genomes and multiple CPR genomes from other studies, indicating a potentially undescribed role for these microorganisms in denitrification. Additionally, glycoside hydrolase (GH) family profiles for our 71 genomes and over 2000 other CPR genomes were analyzed to characterize their carbon-processing potential. Although common trends were present throughout the radiation, differences highlighted potential mechanisms that could allow microorganisms across the CPR to occupy various subsurface niches. For example, members of the Microgenomates superphylum appear to potentially degrade a wider range of carbon substrates than other CPR phyla.

CONCLUSIONS

CPR members are present across a range of environments and often constitute a significant fraction of the microbial population in groundwater systems, particularly. Further sampling of such environments will resolve this portion of the tree of life at finer taxonomic levels, which is essential to solidify functional differences between members that populate this phylogenetically broad region of the tree of life.

Microbial Community Structure and the Persistence of Cyanobacterial Populations in Salt Crusts of the Hyperarid Atacama Desert from Genome-Resolved Metagenomics

Although once thought to be devoid of biology, recent studies have identified salt deposits as oases for life in the hyperarid Atacama Desert. To examine spatial patterns of microbial species and key nutrient sources, we genomically characterized 26 salt crusts from three sites along a fog gradient. The communities are dominated by a large variety of Halobacteriales and Bacteroidetes, plus a few algal and Cyanobacterial species. CRISPR locus analysis suggests the distribution of a single Cyanobacterial population among all sites. This is in stark contrast to the extremely high sample specificity of most other community members. Only present at the highest moisture site is a genomically characterized Thermoplasmatales archaeon (Marine Group II) and six Nanohaloarchaea, one of which is represented by a complete genome. Parcubacteria (OD1) and Saccharibacteria (TM7), not previously reported from hypersaline environments, were found at low abundances. We found no indication of a N2 fixation pathway in the communities, suggesting acquisition of bioavailable nitrogen from atmospherically derived nitrate. Samples cluster by site based on bacterial and archaeal abundance patterns and photosynthetic capacity decreases with increasing distance from the ocean. We conclude that moisture level, controlled by coastal fog intensity, is the strongest driver of community membership.

dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication

The number of microbial genomes sequenced each year is expanding rapidly, in part due to genome-resolved metagenomic studies that routinely recover hundreds of draft-quality genomes. Rapid algorithms have been developed to comprehensively compare large genome sets, but they are not accurate with draft-quality genomes. Here we present dRep, a program that reduces the computational time for pairwise genome comparisons by sequentially applying a fast, inaccurate estimation of genome distance, and a slow, accurate measure of average nucleotide identity. dRep achieves a 28 × increase in speed with perfect recall and precision when benchmarked against previously developed algorithms. We demonstrate the use of dRep for genome recovery from time-series datasets. Each metagenome was assembled separately, and dRep was used to identify groups of essentially identical genomes and select the best genome from each replicate set. This resulted in recovery of significantly more and higher-quality genomes compared to the set recovered using co-assembly.

Simulation of Deepwater Horizon oil plume reveals substrate specialization within a complex community of hydrocarbon degraders

The Deepwater Horizon (DWH) accident released an estimated 4.1 million barrels of oil and 10¹⁰ mol of natural gas into the Gulf of Mexico, forming deep-sea plumes of dispersed oil droplets and dissolved gases that were largely degraded by bacteria. During the course of this 3-mo disaster a series of different bacterial taxa were enriched in succession within deep plumes, but the metabolic capabilities of the different populations that controlled degradation rates of crude oil components are poorly understood. We experimentally reproduced dispersed plumes of fine oil droplets in Gulf of Mexico seawater and successfully replicated the enrichment and succession of the principal oil-degrading bacteria observed during the DWH event. We recovered near-complete genomes, whose phylogeny matched those of the principal biodegrading taxa observed in the field, including the DWH Oceanospirillales (now identified as a Bermanella species), multiple species of Colwellia, Cycloclasticus, and other members of Gammaproteobacteria, Flavobacteria, and Rhodobacteria. Metabolic pathway analysis, combined with hydrocarbon compositional analysis and species abundance data, revealed substrate specialization that explained the successional pattern of oil-degrading bacteria. The fastest-growing bacteria used short-chain alkanes. The analyses also uncovered potential cooperative and competitive relationships, even among close relatives. We conclude that patterns of microbial succession following deep ocean hydrocarbon blowouts are predictable and primarily driven by the availability of liquid petroleum hydrocarbons rather than natural gases.

Microbial potential for carbon and nutrient cycling in a geogenic supercritical carbon dioxide reservoir

Microorganisms catalyze carbon cycling and biogeochemical reactions in the deep subsurface and thus may be expected to influence the fate of injected supercritical (sc) CO2 following geological carbon sequestration (GCS). We hypothesized that natural subsurface scCO2 reservoirs, which serve as analogs for the long-term fate of sequestered scCO2 , harbor a 'deep carbonated biosphere' with carbon cycling potential. We sampled subsurface fluids from scCO2 -water separators at a natural scCO2 reservoir at McElmo Dome, Colorado for analysis of 16S rRNA gene diversity and metagenome content. Sequence annotations indicated dominance of Sulfurospirillum, Rhizobium, Desulfovibrio and four members of the Clostridiales family. Genomes extracted from metagenomes using homology and compositional approaches revealed diverse mechanisms for growth and nutrient cycling, including pathways for CO2 and N2 fixation, anaerobic respiration, sulfur oxidation, fermentation and potential for metabolic syntrophy. Differences in biogeochemical potential between two production well communities were consistent with differences in fluid chemical profiles, suggesting a potential link between microbial activity and geochemistry. The existence of a microbial ecosystem associated with the McElmo Dome scCO2 reservoir indicates that potential impacts of the deep biosphere on CO2 fate and transport should be taken into consideration as a component of GCS planning and modelling.

RNA-Dependent Cysteine Biosynthesis in Bacteria and Archaea

The diversity of the genetic code systems used by microbes on earth is yet to be elucidated. It is known that certain methanogenic archaea employ an alternative system for cysteine (Cys) biosynthesis and encoding; tRNA^Cys is first acylated with phosphoserine (Sep) by O-phosphoseryl-tRNA synthetase (SepRS) and then converted to Cys-tRNA^Cys by Sep-tRNA:Cys-tRNA synthase (SepCysS). In this study, we searched all genomic and metagenomic protein sequence data in the Integrated Microbial Genomes (IMG) system and at the NCBI to reveal new clades of SepRS and SepCysS proteins belonging to diverse archaea in the four major groups (DPANN, Euryarchaeota, TACK, and Asgard) and two groups of bacteria (“Candidatus Parcubacteria” and Chloroflexi). Bacterial SepRS and SepCysS charged bacterial tRNA^Cys species with cysteine in vitro. Homologs of SepCysE, a scaffold protein facilitating SepRS⋅SepCysS complex assembly in Euryarchaeota class I methanogens, are found in a few groups of TACK and Asgard archaea, whereas the C-terminally truncated homologs exist fused or genetically coupled with diverse SepCysS species. Investigation of the selenocysteine (Sec)- and pyrrolysine (Pyl)-utilizing traits in SepRS-utilizing archaea and bacteria revealed that the archaea carrying full-length SepCysE employ Sec and that SepRS is often found in Pyl-utilizing archaea and Chloroflexi bacteria. We discuss possible contributions of the SepRS-SepCysS system for sulfur assimilation, methanogenesis, and other metabolic processes requiring large amounts of iron-sulfur enzymes or Pyl-containing enzymes.

Comprehensive analyses of all genomic and metagenomic protein sequence data in public databases revealed the distribution and evolution of an alternative cysteine-encoding system in diverse archaea and bacteria. The finding that the SepRS-SepCysS-SepCysE- and the selenocysteine-encoding systems are shared by the Euryarchaeota class I methanogens, the Crenarchaeota AK8/W8A-19 group, and an Asgard archaeon suggests that ancient archaea may have used both systems. In contrast, bacteria may have obtained the SepRS-SepCysS system from archaea. The SepRS-SepCysS system sometimes coexists with a pyrrolysine-encoding system in both archaea and bacteria. Our results provide additional bioinformatic evidence for the contribution of the SepRS-SepCysS system for sulfur assimilation and diverse metabolisms which require vast amounts of iron-sulfur enzymes and proteins. Among these biological activities, methanogenesis, methylamine metabolism, and organohalide respiration may have local and global effects on earth. Taken together, uncultured bacteria and archaea provide an expanded record of the evolution of the genetic code.

Poorly known microbial taxa dominate the microbiome of permafrost thaw ponds

In the transition zone of the shifting permafrost border, thaw ponds emerge as hotspots of microbial activity, processing the ancient carbon freed from the permafrost. We analyzed the microbial succession across a gradient of recently emerged to older ponds using three molecular markers: one universal, one bacterial and one fungal. Age was a major modulator of the microbial community of the thaw ponds. Surprisingly, typical freshwater taxa comprised only a small fraction of the community. Instead, thaw ponds of all age classes were dominated by enigmatic bacterial and fungal phyla. Our results on permafrost thaw ponds lead to a revised perception of the thaw pond ecosystem and their microbes, with potential implications for carbon and nutrient cycling in this increasingly important class of freshwaters.

Potential for microbial H2 and metal transformations associated with novel bacteria and archaea in deep terrestrial subsurface sediments

Geological sequestration in deep underground repositories is the prevailing proposed route for radioactive waste disposal. After the disposal of radioactive waste in the subsurface, H₂ may be produced by corrosion of steel and, ultimately, radionuclides will be exposed to the surrounding environment. To evaluate the potential for microbial activities to impact disposal systems, we explored the microbial community structure and metabolic functions of a sediment-hosted ecosystem at the Horonobe Underground Research Laboratory, Hokkaido, Japan. Overall, we found that the ecosystem hosted organisms from diverse lineages, including many from the phyla that lack isolated representatives. The majority of organisms can metabolize H₂, often via oxidative [NiFe] hydrogenases or electron-bifurcating [FeFe] hydrogenases that enable ferredoxin-based pathways, including the ion motive Rnf complex. Many organisms implicated in H₂ metabolism are also predicted to catalyze carbon, nitrogen, iron and sulfur transformations. Notably, iron-based metabolism is predicted in a novel lineage of Actinobacteria and in a putative methane-oxidizing ANME-2d archaeon. We infer an ecological model that links microorganisms to sediment-derived resources and predict potential impacts of microbial activity on H₂ consumption and retardation of radionuclide migration.

New CRISPR-Cas systems from uncultivated microbes

CRISPR-Cas systems provide microbes with adaptive immunity by employing short sequences, termed spacers, that guide Cas proteins to cleave foreign DNA^1,2. Class 2 CRISPR-Cas systems are streamlined versions in which a single Cas protein bound to RNA recognizes and cleaves targeted sequences^3,4. The programmable nature of these minimal systems has enabled their repurposing as a versatile technology that is broadly revolutionizing biological and clinical research⁵. However, current CRISPR-Cas technologies are based solely on systems from isolated bacteria, leaving untapped the vast majority of enzymes from organisms that have not been cultured. Metagenomics, the sequencing of DNA extracted from natural microbial communities, provides access to the genetic material of a huge array of uncultivated organisms^6,7. Here, using genome-resolved metagenomics, we identified novel CRISPR-Cas systems, including the first reported Cas9 in the archaeal domain of life. This divergent Cas9 protein was found in little-studied nanoarchaea as part of an active CRISPR-Cas system. In bacteria, we discovered two previously unknown systems, CRISPR-CasX and CRISPR-CasY, which are among the most compact systems yet identified. Notably, all required functional components were identified by metagenomics, enabling validation of robust in vivo RNA-guided DNA interference activity in E. coli. Interrogation of environmental microbial communities combined with in vivo experiments allows access to an unprecedented diversity of genomes whose content will expand the repertoire of microbe-based biotechnologies.

Proteogenomic analyses indicate bacterial methylotrophy and archaeal heterotrophy are prevalent below the grass root zone

Annually, half of all plant-derived carbon is added to soil where it is microbially respired to CO₂. However, understanding of the microbiology of this process is limited because most culture-independent methods cannot link metabolic processes to the organisms present, and this link to causative agents is necessary to predict the results of perturbations on the system. We collected soil samples at two sub-root depths (10–20 cm and 30–40 cm) before and after a rainfall-driven nutrient perturbation event in a Northern California grassland that experiences a Mediterranean climate. From ten samples, we reconstructed 198 metagenome-assembled genomes that represent all major phylotypes. We also quantified 6,835 proteins and 175 metabolites and showed that after the rain event the concentrations of many sugars and amino acids approach zero at the base of the soil profile. Unexpectedly, the genomes of novel members of the Gemmatimonadetes and Candidate Phylum Rokubacteria phyla encode pathways for methylotrophy. We infer that these abundant organisms contribute substantially to carbon turnover in the soil, given that methylotrophy proteins were among the most abundant proteins in the proteome. Previously undescribed Bathyarchaeota and Thermoplasmatales archaea are abundant in deeper soil horizons and are inferred to contribute appreciably to aromatic amino acid degradation. Many of the other bacteria appear to breakdown other components of plant biomass, as evidenced by the prevalence of various sugar and amino acid transporters and corresponding hydrolyzing machinery in the proteome. Overall, our work provides organism-resolved insight into the spatial distribution of bacteria and archaea whose activities combine to degrade plant-derived organics, limiting the transport of methanol, amino acids and sugars into underlying weathered rock. The new insights into the soil carbon cycle during an intense period of carbon turnover, including biogeochemical roles to previously little known soil microbes, were made possible via the combination of metagenomics, proteomics, and metabolomics.

Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system

The subterranean world hosts up to one-fifth of all biomass, including microbial communities that drive transformations central to Earth's biogeochemical cycles. However, little is known about how complex microbial communities in such environments are structured, and how inter-organism interactions shape ecosystem function. Here we apply terabase-scale cultivation-independent metagenomics to aquifer sediments and groundwater, and reconstruct 2,540 draft-quality, near-complete and complete strain-resolved genomes that represent the majority of known bacterial phyla as well as 47 newly discovered phylum-level lineages. Metabolic analyses spanning this vast phylogenetic diversity and representing up to 36% of organisms detected in the system are used to document the distribution of pathways in coexisting organisms. Consistent with prior findings indicating metabolic handoffs in simple consortia, we find that few organisms within the community can conduct multiple sequential redox transformations. As environmental conditions change, different assemblages of organisms are selected for, altering linkages among the major biogeochemical cycles.

Culture Independent Genomic Comparisons Reveal Environmental Adaptations for Altiarchaeales

The recently proposed candidatus order Altiarchaeales remains an uncultured archaeal lineage composed of genetically diverse, globally widespread organisms frequently observed in anoxic subsurface environments. In spite of 15 years of studies on the psychrophilic biofilm-producing Candidatus Altiarchaeum hamiconexum and its close relatives, very little is known about the phylogenetic and functional diversity of the widespread free-living marine members of this taxon. From methanogenic sediments in the White Oak River Estuary, NC, USA, we sequenced a single cell amplified genome (SAG), WOR_SM1_SCG, and used it to identify and refine two high-quality genomes from metagenomes, WOR_SM1_79 and WOR_SM1_86-2, from the same site. These three genomic reconstructions form a monophyletic group, which also includes three previously published genomes from metagenomes from terrestrial springs and a SAG from Sakinaw Lake in a group previously designated as pMC2A384. A synapomorphic mutation in the Altiarchaeales tRNA synthetase β subunit, pheT, caused the protein to be encoded as two subunits at non-adjacent loci. Consistent with the terrestrial spring clades, our estuarine genomes contained a near-complete autotrophic metabolism, H2 or CO as potential electron donors, a reductive acetyl-CoA pathway for carbon fixation, and methylotroph-like NADP(H)-dependent dehydrogenase. Phylogenies based on 16S rRNA genes and concatenated conserved proteins identified two distinct sub-clades of Altiarchaeales, Alti-1 populated by organisms from actively flowing springs, and Alti-2 which was more widespread, diverse, and not associated with visible mats. The core Alti-1 genome suggested Alti-1 is adapted for the stream environment with lipopolysaccharide production capacity and extracellular hami structures. The core Alti-2 genome suggested members of this clade are free-living with distinct mechanisms for energy maintenance, motility, osmoregulation, and sulfur redox reactions. These data suggested that the hamus structures found in Candidatus Altiarchaeum hamiconexum are not present outside of stream-adapted Altiarchaeales. Homologs to a Na(+) transporter and membrane bound coenzyme A disulfide reductase that were unique to the brackish sediment Alti-2 genomes, could indicate adaptations to the estuarine, sulfur-rich environment.