Investigating the chemolithoautotrophic and formate metabolism of Nitrospira moscoviensis by constraint-based metabolic modeling and 13C-tracer analysis

Nitrite-oxidizing bacteria belonging to the genus Nitrospira mediate a key step in nitrification and play important roles in the biogeochemical nitrogen cycle and wastewater treatment. While these organisms have recently been shown to exhibit metabolic flexibility beyond their chemolithoautotrophic lifestyle, including the use of simple organic compounds to fuel their energy metabolism, the metabolic networks controlling their autotrophic and mixotrophic growth remain poorly understood. Here, we reconstructed a genome-scale metabolic model for Nitrospira moscoviensis (iNmo686) and used constraint-based analysis to evaluate the metabolic networks controlling autotrophic and formatotrophic growth on nitrite and formate, respectively. Subsequently, proteomic analysis and 13C-tracer experiments with bicarbonate and formate coupled to metabolomic analysis were performed to experimentally validate model predictions. Our findings support that N. moscoviensis uses the reductive tricarboxylic acid cycle for CO2 fixation. We also show that N. moscoviensis can indirectly use formate as a carbon source by oxidizing it first to CO2 followed by reassimilation, rather than direct incorporation via the reductive glycine pathway. Our study offers the first measurements of Nitrospira’s in vivo central carbon metabolism and provides a quantitative tool that can be used for understanding and predicting their metabolic processes.ImportanceNitrospira are globally abundant nitrifying bacteria in soil and aquatic ecosystems and wastewater treatment plants, where they control the oxidation of nitrite to nitrate. Despite their critical contribution to nitrogen cycling across diverse environments, detailed understanding of their metabolic network and prediction of their function under different environmental conditions remains a major challenge. Here, we provide the first constraint-based metabolic model of N. moscoviensis representing the ubiquitous Nitrospira lineage II and subsequently validate this model using proteomics and 13C-tracers combined with intracellular metabolomic analysis. The resulting genome-scale model will serve as a knowledge base of Nitrospira metabolism and lays the foundation for quantitative systems biology studies of these globally important nitrite- oxidizing bacteria.

Functional diversity of microbial communities in pristine aquifers inferred by PLFA- and sequencing-based approaches

Microorganisms in groundwater play an important role in aquifer biogeochemical cycles and water quality. However, the mechanisms linking the functional diversity of microbial populations and the groundwater physico-chemistry are still not well understood due to the complexity of interactions between surface and subsurface. Within the framework of Hainich (north-western Thuringia, central Germany) Critical Zone Exploratory of the Collaborative Research Centre AquaDiva, we used the relative abundances of phospholipid-derived fatty acids (PLFAs) to link specific biochemical markers within the microbial communities to the spatio-temporal changes of the groundwater physico-chemistry. The functional diversities of the microbial communities were mainly correlated with groundwater chemistry, including dissolved O2, Fet and NH4+ concentrations. Abundances of PLFAs derived from eukaryotes and potential nitrite-oxidizing bacteria (11Me16:0 as biomarker for Nitrospira moscoviensis) were high at sites with elevated O2 concentration where groundwater recharge supplies bioavailable substrates. In anoxic groundwaters more rich in Fet, PLFAs abundant in sulfate-reducing bacteria (SRB), iron-reducing bacteria and fungi increased with Fet and HCO3− concentrations, suggesting the occurrence of active iron reduction and the possible role of fungi in meditating iron solubilization and transport in those aquifer domains. In more NH4+-rich anoxic groundwaters, anammox bacteria and SRB-derived PLFAs increased with NH4+ concentration, further evidencing the dependence of the anammox process on ammonium concentration and potential links between SRB and anammox bacteria. Additional support of the PLFA-based bacterial communities was found in DNA- and RNA-based Illumina MiSeq amplicon sequencing of bacterial 16S rRNA genes, which showed high predominance of nitrite-oxidizing bacteria Nitrospira, e.g. Nitrospira moscoviensis, in oxic aquifer zones and of anammox bacteria in more NH4+-rich anoxic groundwater. Higher relative abundances of sequence reads in the RNA-based datasets affiliated with iron-reducing bacteria in more Fet-rich groundwater supported the occurrence of active dissimilatory iron reduction. The functional diversity of the microbial communities in the biogeochemically distinct groundwater assemblages can be largely attributed to the redox conditions linked to changes in bioavailable substrates and input of substrates with the seepage. Our results demonstrate the power of complementary information derived from PLFA-based and sequencing-based approaches.

Functional diversity of microbial communities in pristine aquifers inferred by PLFA – and sequencing – based approaches

Microorganisms in groundwater play an important role in aquifer biogeochemical cycles and water quality. However, the mechanisms linking the functional diversity of microbial populations and the groundwater physicochemistry are still not well understood due to the complexity of interactions between surface and subsurface. Here, we used phospholipid fatty acids (PLFAs) relative abundances to link specific biochemical markers within the microbial communities to the spatio-temporal changes of the groundwater physicochemistry. PLFAs were isolated from groundwater of two physicochemically distinct aquifer assemblages in central Germany (Thuringia). The functional diversities of the microbial communities were mainly correlated with groundwater chemistry, including dissolved O2, Fet and NH4+ concentrations. Abundances of PLFAs derived from eukaryotes and potential nitrite oxidizing bacteria (11MeC16:0 as biomarker for Nitrospira moscoviensis) were high at sites with elevated O2 concentration where groundwater recharge supplies both bioavailable organic substrates and NH4+ needed to sustain heterotrophic growth and nitrification processes. In anoxic groundwaters more rich in Fet, PLFAs abundant in sulphate reducing bacteria (SRB), iron-reducing bacteria and fungi increased with Fet and HCO3− concentrations suggesting the occurrence of active iron-reduction and the possible role of fungi in meditating iron solubilisation and transport in those aquifer domains. In NH4+ richer anoxic groundwaters, anammox bacteria and SRB- derived PLFAs increased with NH4+ concentration further evidencing the dependence of the anammox process on ammonium concentration and potential links between SRB and anammox bacteria. Additional support of the PLFA-based bacterial communities was found in DNA and RNA-based Illumina MiSeq amplicon sequencing of bacterial 16S rRNA genes, which evidenced high predominance of nitrite-oxidizing bacteria Nitrospira e.g. Nitrospira moscoviensis in oxic zones of the aquifers and of anammox bacteria in NH4+ richer anoxic groundwater. Higher relative abundances of sequence reads in the RNA-based data sets affiliated with iron-reducing bacteria in Fet richer groundwater supported the occurrence of active dissimilatory iron-reduction. The functional diversity of the microbial communities in these biogeochemically distinct groundwater assemblages can be largely attributed to the redox conditions linked to changes in bioavailable substrates and input of substrates with the seepage. Our results demonstrate the power of complementary information derived from PLFA-based and sequencing-based approaches.