Metabolic Diversity and Aero-Tolerance in Anammox Bacteria from Geochemically Distinct Aquifers

Metagenomic evidence of a novel anammox community in a cold aquifer with high nitrogen pollution

The process of anaerobic ammonium oxidation by nitrite (anammox) is a globally essential part of N cycle. To date, 8 Candidatus genera and more than 22 species of anammox bacteria have been discovered in various anthropogenic and natural habitats, including nitrogen-polluted aquifers. In this work, anammox bacteria were detected for the first time in the groundwater ecosystem with high anthropogenic nitrogen pollution (up to 1760 mg NO3--N/L and 280 mg NH4+-N/L) and low year-round temperature (7-8 °C) in the zone of a uranium sludge repository. Further metagenomic analysis resulted in retrieval of metagenome-assembled genomes of 4 distinct anammox bacteria: a new genus named Ca. Frigussubterria, new species in Ca. Kuenenia, and two strains of a new species in Ca. Scalindua. Analysis of the genomes revealed essential genes involved in anammox metabolism. Both strains of Ca. Scalindua chemeplantae had a high copy number of genes encoding the cold shock proteins CspA/B, which can also function as an antifreeze protein (CspB). Ca. Kuenenia glazoviensis and Ca. Frigussubterria udmurtiae were abundant in less N-polluted site, while Ca. Scalindua chemeplantae inhabited both sites. Genes for urea utilization, reduction of insoluble Fe2O3 or MnO2, assimilatory sulfate reduction, reactive oxygen detoxification, nitrate reduction to ammonium, and putatively arsenate respiration were found. These findings enrich knowledge of the functional and phylogenetic diversity of anammox bacteria and improve understanding of the nitrogen cycle in polluted aquifers.

Metatranscriptomic Analysis Reveals Synergistic Activities of Comammox and Anammox Bacteria in Full-Scale Attached Growth Nitrogen Removal System

Leveraging comammox Nitrospira and anammox bacteria for shortcut nitrogen removal can drastically lower the carbon footprint of wastewater treatment facilities by decreasing aeration energy, carbon, alkalinity, and tank volume requirements while also potentially reducing nitrous oxide emissions. However, their co-occurrence as dominant nitrifying bacteria is rarely reported in full-scale wastewater treatment. As a result, there is a poor understanding of how operational parameters, in particular, dissolved oxygen, impact their activity and synergistic behavior. Here, we report the impact of dissolved oxygen concentration (DO = 2, 4, 6 mg/L) on the microbial community's transcriptomic expression in a full-scale integrated fixed film activated sludge (IFAS) municipal wastewater treatment facility where nitrogen removal is predominantly performed by comammox Nitrospira and anammox bacterial populations. 16S rRNA transcript compositions revealed anammox bacteria and Nitrospira were significantly more active in IFAS biofilms compared to suspended sludge biomass. In IFAS biofilms, anammox bacteria significantly increased hzo expression at lower dissolved oxygen concentrations and this increase was highly correlated with the amoA expression levels of comammox bacteria. Interestingly, the genes involved in nitrite oxidation by comammox bacteria were significantly more upregulated, relative to the genes involved in ammonia oxidation with decreasing dissolved oxygen concentrations. Ultimately, our findings suggest that comammox Nitrospira supplies anammox bacteria with nitrite via ammonia oxidation and that this synergistic behavior is dependent on dissolved oxygen concentrations.

Predicted metabolic roles and stress responses provide insights into candidate phyla Hydrogenedentota and Sumerlaeota as members of the rare biosphere in biofilms from various environments

Pustular mats from Shark Bay, Western Australia, host complex microbial communities bound within an organic matrix. These mats harbour many poorly characterized organisms with low relative abundances (<1%), such as candidate phyla Hydrogenedentota and Sumerlaeota. Here, we aim to constrain the metabolism and physiology of these candidate phyla by analyzing two representative metagenome-assembled genomes (MAGs) from a pustular mat. Metabolic reconstructions of these MAGs suggest facultatively anaerobic, chemoorganotrophic lifestyles of both organisms and predict that both MAGs can metabolize a diversity of carbohydrate substrates. Ca. Sumerlaeota possesses genes involved in degrading chitin, cellulose and other polysaccharides, while Ca. Hydrogenedentota can metabolize cellulose derivatives in addition to glycerol, fatty acids and phosphonates. Both Ca. phyla can respond to nitrosative stress and participate in nitrogen metabolism. Metabolic comparisons of MAGs from Shark Bay and those from various polyextreme environments (i.e., hot springs, hydrothermal vents, subsurface waters, anaerobic digesters, etc.) reveal similar metabolic capabilities and adaptations to hypersalinity, oxidative stress, antibiotics, UV radiation, nitrosative stress, heavy metal toxicity and life in surface-attached communities. These adaptations and capabilities may account for the widespread nature of these organisms and their contributions to biofilm communities in a range of extreme surface and subsurface environments.

High niche specificity and host genetic diversity of groundwater viruses

Viruses are key members of microbial communities that exert control over host abundance and metabolism, thereby influencing ecosystem processes and biogeochemical cycles. Aquifers are known to host taxonomically diverse microbial life, yet little is known about viruses infecting groundwater microbial communities. Here, we analysed 16 metagenomes from a broad range of groundwater physicochemistries. We recovered 1571 viral genomes that clustered into 468 high-quality viral operational taxonomic units. At least 15% were observed to be transcriptionally active, although lysis was likely constrained by the resource-limited groundwater environment. Most were unclassified (95%), and the remaining 5% were Caudoviricetes. Comparisons with viruses inhabiting other aquifers revealed no shared species, indicating substantial unexplored viral diversity. In silico predictions linked 22.4% of the viruses to microbial host populations, including to ultra-small prokaryotes, such as Patescibacteria and Nanoarchaeota. Many predicted hosts were associated with the biogeochemical cycling of carbon, nitrogen, and sulfur. Metabolic predictions revealed the presence of 205 putative auxiliary metabolic genes, involved in diverse processes associated with the utilization of the host's intracellular resources for biosynthesis and transformation reactions, including those involved in nucleotide sugar, glycan, cofactor, and vitamin metabolism. Viruses, prokaryotes overall, and predicted prokaryotic hosts exhibited narrow spatial distributions, and relative abundance correlations with the same groundwater parameters (e.g. dissolved oxygen, nitrate, and iron), consistent with host control over viral distributions. Results provide insights into underexplored groundwater viruses, and indicate the large extent to which viruses may manipulate microbial communities and biogeochemistry in the terrestrial subsurface.

Implications for nitrogen and sulphur cycles: phylogeny and niche-range of Nitrospirota in terrestrial aquifers

Increasing evidence suggests Nitrospirota are important contributors to aquatic and subsurface nitrogen and sulphur cycles. We determined the phylogenetic and ecological niche associations of Nitrospirota colonizing terrestrial aquifers. Nitrospirota compositions were determined across 59 groundwater wells. Distributions were strongly influenced by oxygen availability in groundwater, marked by a trade-off between aerobic (Nitrospira, Leptospirillum) and anaerobic (Thermodesulfovibrionia, unclassified) lineages. Seven Nitrospirota metagenome-assembled genomes (MAGs), or populations, were recovered from a subset of wells, including three from the recently designated class 9FT-COMBO-42-15. Most were relatively more abundant and transcriptionally active in dysoxic groundwater. These MAGs were analysed with 743 other Nitrospirota genomes. Results illustrate the predominance of certain lineages in aquifers (e.g. non-nitrifying Nitrospiria, classes 9FT-COMBO-42-15 and UBA9217, and Thermodesulfovibrionales family UBA1546). These lineages are characterized by mechanisms for nitrate reduction and sulphur cycling, and, excluding Nitrospiria, the Wood-Ljungdahl pathway, consistent with carbon-limited, low-oxygen, and sulphur-rich aquifer conditions. Class 9FT-COMBO-42-15 is a sister clade of Nitrospiria and comprises two families spanning a transition in carbon fixation approaches: f_HDB-SIOIB13 encodes rTCA (like Nitrospiria) and f_9FT-COMBO-42-15 encodes Wood-Ljungdahl CO dehydrogenase (like Thermodesulfovibrionia and UBA9217). The 9FT-COMBO-42-15 family is further differentiated by its capacity for sulphur oxidation (via DsrABEFH and SoxXAYZB) and dissimilatory nitrate reduction to ammonium, and gene transcription indicated active coupling of nitrogen and sulphur cycles by f_9FT-COMBO-42-15 in dysoxic groundwater. Overall, results indicate that Nitrospirota are widely distributed in groundwater and that oxygen availability drives the spatial differentiation of lineages with ecologically distinct roles related to nitrogen and sulphur metabolism.

Metabolic and phylogenetic diversity in the phylum Nitrospinota revealed by comparative genome analyses

The most abundant known nitrite-oxidizing bacteria in the marine water column belong to the phylum Nitrospinota. Despite their importance in marine nitrogen cycling and primary production, there are only few cultured representatives that all belong to the class Nitrospinia. Moreover, although Nitrospinota were traditionally thought to be restricted to marine environments, metagenome-assembled genomes have also been recovered from groundwater. Over the recent years, metagenomic sequencing has led to the discovery of several novel classes of Nitrospinota (UBA9942, UBA7883, 2-12-FULL-45-22, JACRGO01, JADGAW01), which remain uncultivated and have not been analyzed in detail. Here, we analyzed a nonredundant set of 98 Nitrospinota genomes with focus on these understudied Nitrospinota classes and compared their metabolic profiles to get insights into their potential role in biogeochemical element cycling. Based on phylogenomic analysis and average amino acid identities, the highly diverse phylum Nitrospinota could be divided into at least 33 different genera, partly with quite distinct metabolic capacities. Our analysis shows that not all Nitrospinota are nitrite oxidizers and that members of this phylum have the genomic potential to use sulfide and hydrogen for energy conservation. This study expands our knowledge of the phylogeny and potential ecophysiology of the phylum Nitrospinota and offers new avenues for the isolation and cultivation of these elusive bacteria.

Nitrogen fixation and diazotroph diversity in groundwater systems

Biological nitrogen fixation (BNF), the conversion of N2 into bioavailable nitrogen (N), is the main process for replenishing N loss in the biosphere. However, BNF in groundwater systems remains poorly understood. In this study, we examined the activity, abundance, and community composition of diazotrophs in groundwater in the Hetao Plain of Inner Mongolia using 15N tracing methods, reverse transcription qPCR (RT-qPCR), and metagenomic/metatranscriptomic analyses. 15N2 tracing incubation of near in situ groundwater (9.5-585.4 nmol N L-1 h-1) and N2-fixer enrichment and isolates (13.2-1728.4 nmol N g-1 h-1, as directly verified by single-cell resonance Raman spectroscopy), suggested that BNF is a non-negligible source of N in groundwater in this region. The expression of nifH genes ranged from 3.4 × 103 to 1.2 × 106 copies L-1 and was tightly correlated with dissolved oxygen (DO), Fe(II), and NH4+. Diazotrophs in groundwater were chiefly aerobes or facultative anaerobes, dominated by Stutzerimonas, Pseudomonas, Paraburkholderia, Klebsiella, Rhodopseudomonas, Azoarcus, and additional uncultured populations. Active diazotrophs, which prefer reducing conditions, were more metabolically diverse and potentially associated with nitrification, sulfur/arsenic mobilization, Fe(II) transport, and CH4 oxidation. Our results highlight the importance of diazotrophs in subsurface geochemical cycles.

Ultra-small bacteria and archaea exhibit genetic flexibility towards groundwater oxygen content, and adaptations for attached or planktonic lifestyles

Aquifers are populated by highly diverse microbial communities, including unusually small bacteria and archaea. The recently described Patescibacteria (or Candidate Phyla Radiation) and DPANN radiation are characterized by ultra-small cell and genomes sizes, resulting in limited metabolic capacities and probable dependency on other organisms to survive. We applied a multi-omics approach to characterize the ultra-small microbial communities over a wide range of aquifer groundwater chemistries. Results expand the known global range of these unusual organisms, demonstrate the wide geographical range of over 11,000 subsurface-adapted Patescibacteria, Dependentiae and DPANN archaea, and indicate that prokaryotes with ultra-small genomes and minimalistic metabolism are a characteristic feature of the terrestrial subsurface. Community composition and metabolic activities were largely shaped by water oxygen content, while highly site-specific relative abundance profiles were driven by a combination of groundwater physicochemistries (pH, nitrate-N, dissolved organic carbon). We provide insights into the activity of ultra-small prokaryotes with evidence that they are major contributors to groundwater community transcriptional activity. Ultra-small prokaryotes exhibited genetic flexibility with respect to groundwater oxygen content, and transcriptionally distinct responses, including proportionally greater transcription invested into amino acid and lipid metabolism and signal transduction in oxic groundwater, along with differences in taxa transcriptionally active. Those associated with sediments differed from planktonic counterparts in species composition and transcriptional activity, and exhibited metabolic adaptations reflecting a surface-associated lifestyle. Finally, results showed that groups of phylogenetically diverse ultra-small organisms co-occurred strongly across sites, indicating shared preferences for groundwater conditions.

Characterization of Enrichment Cultures of Anammox, Nitrifying and Denitrifying Bacteria Obtained from a Cold, Heavily Nitrogen-Polluted Aquifer

Anammox bacteria related to Candidatus Scalindua were recently discovered in a cold (7.5 °C) aquifer near sludge repositories containing solid wastes of uranium and processed polymetallic concentrate. Groundwater has a very high level of nitrate and ammonia pollution (up to 10 and 0.5 g/L, respectively) and a very low content of organic carbon (2.5 mg/L). To assess the potential for bioremediation of polluted groundwater in situ, enrichment cultures of anammox, nitrifying, and denitrifying bacteria were obtained and analyzed. Fed-batch enrichment of anammox bacteria was not successful. Stable removal of ammonium and nitrite (up to 100%) was achieved in a continuous-flow reactor packed with a nonwoven fabric at 15 °C, and enrichment in anammox bacteria was confirmed by FISH and qPCR assays. The relatively low total N removal efficiency (up to 55%) was due to nonstoichiometric nitrate buildup. This phenomenon can be explained by a shift in the metabolism of anammox bacteria towards the production of more nitrates and less N₂ at low temperatures compared to the canonical stoichiometry. In addition, the too high an estimate of specific anammox activity suggests that N cycle microbial groups other than anammox bacteria may have contributed significantly to N removal. Stable nitrite production was observed in the denitrifying enrichment culture, while no "conventional" nitrifiers were found in the corresponding enrichment cultures. Xanthomonadaceae was a common taxon for all microbial communities, indicating its exclusive role in this ecosystem. This study opens up new knowledge about the metabolic capabilities of N cycle bacteria and potential approaches for sustainable bioremediation of heavily N-polluted cold ecosystems.

Novel and unusual genes for nitrogen and metal cycling in Planctomycetota- and KSB1-affiliated metagenome-assembled genomes reconstructed from a marine subsea tunnel

The Oslofjord subsea road tunnel is a unique environment in which the typically anoxic marine deep subsurface is exposed to oxygen. Concrete biodeterioration and steel corrosion in the tunnel have been linked to the growth of iron- and manganese-oxidizing biofilms in areas of saline water seepage. Surprisingly, previous 16S rRNA gene surveys of biofilm samples revealed microbial communities dominated by sequences affiliated with nitrogen-cycling microorganisms. This study aimed to identify microbial genomes with metabolic potential for novel nitrogen- and metal-cycling reactions, representing biofilm microorganisms that could link these cycles and play a role in concrete biodeterioration. We reconstructed 33 abundant, novel metagenome-assembled genomes (MAGs) affiliated with the phylum Planctomycetota and the candidate phylum KSB1. We identified novel and unusual genes and gene clusters in these MAGs related to anaerobic ammonium oxidation, nitrite oxidation, and other nitrogen-cycling reactions. Additionally, 26 of 33 MAGs also had the potential for iron, manganese, and arsenite cycling, suggesting that bacteria represented by these genomes might couple these reactions. Our results expand the diversity of microorganisms putatively involved in nitrogen and metal cycling, and contribute to our understanding of potential biofilm impacts on built infrastructure.

Nitrogen cycling and microbial cooperation in the terrestrial subsurface

The nitrogen cycle plays a major role in aquatic nitrogen transformations, including in the terrestrial subsurface. However, the variety of transformations remains understudied. To determine how nitrogen cycling microorganisms respond to different aquifer chemistries, we sampled groundwater with varying nutrient and oxygen contents. Genes and transcripts involved in major nitrogen-cycling pathways were quantified from 55 and 26 sites, respectively, and metagenomes and metatranscriptomes were analyzed from a subset of oxic and dysoxic sites (0.3-1.1 mg/L bulk dissolved oxygen). Nitrogen-cycling mechanisms (e.g. ammonia oxidation, denitrification, dissimilatory nitrate reduction to ammonium) were prevalent and highly redundant, regardless of site-specific physicochemistry or nitrate availability, and present in 40% of reconstructed genomes, suggesting that nitrogen cycling is a core function of aquifer communities. Transcriptional activity for nitrification, denitrification, nitrite-dependent anaerobic methane oxidation and anaerobic ammonia oxidation (anammox) occurred simultaneously in oxic and dysoxic groundwater, indicating the availability of oxic-anoxic interfaces. Concurrent activity by these microorganisms indicates potential synergisms through metabolite exchange across these interfaces (e.g. nitrite and oxygen). Fragmented denitrification pathway encoding and transcription was widespread among groundwater bacteria, although a considerable proportion of associated transcriptional activity was driven by complete denitrifiers, especially under dysoxic conditions. Despite large differences in transcription, the capacity for the final steps of denitrification was largely invariant to aquifer conditions, and most genes and transcripts encoding N₂O reductases were the atypical Sec-dependant type, suggesting energy-efficiency prioritization. Results provide insights into the capacity for cooperative relationships in groundwater communities, and the richness and complexity of metabolic mechanisms leading to the loss of fixed nitrogen.