Nitrogen removal from activated sludge wastewater treatment systems is an energy-intensive process due to the large aeration requirement for nitrification. This energy footprint could be minimized with engineering control strategies that wash out nitrite-oxidizing bacteria (NOB) to limit oxygen demands.

Conventional bioprocess models for wastewater treatment are based on aggregated bulk biomass concentrations and do not incorporate microbial physiological diversity. Such a broad aggregation of microbial functional groups can fail to predict ecosystem dynamics when high levels of physiological diversity exist within trophic guilds. For instance, functional diversity among nitrite-oxidizing bacteria (NOB) can obfuscate engineering strategies for their out-selection in activated sludge (AS), which is desirable to promote energy-efficient nitrogen removal. Here, we hypothesized that different NOB populations within AS can have different physiological traits that drive process performance, which we tested by estimating biokinetic growth parameters using a combination of highly replicated respirometry, genome-resolved metagenomics, and process modeling. A lab-scale AS reactor subjected to a selective pressure for over 90 days experienced resilience of NOB activity. We recovered three coexisting Nitrospira population genomes belonging to two sublineages, which exhibited distinct growth strategies and underwent a compositional shift following the selective pressure. A trait-based process model calibrated at the NOB genus level better predicted nitrite accumulation than a conventional process model calibrated at the NOB guild level. This work demonstrates that trait-based modeling can be leveraged to improve our prediction, control, and design of functionally diverse microbiomes driving key environmental biotechnologies.

Using genome-centric metagenomics, we discovered a new Methanothrix harundinacea subspecies that appears to be a halotolerant acetoclastic methanogen with the flexibility for adaptation in the anaerobic digestion process both at low (5 g/L Na + ) and high salinity conditions (20 g/L Na + ). Annotation of the recovered M. harundinacea genome revealed salinity stress-related functions, including the modification of EPS glycoconjugates and the production of compatible solutes.

Stable isotope probing (SIP) facilitates culture-independent identification of active microbial populations within complex ecosystems through isotopic enrichment of nucleic acids. Many SIP studies rely on 16S rRNA sequences to identify active taxa but connecting these sequences to specific bacterial genomes is often challenging. Here, we describe a standardized laboratory and analysis framework to quantify isotopic enrichment on a per-genome basis using shotgun metagenomics instead of 16S rRNA sequencing. To develop this framework, we explored various sample processing and analysis approaches using a designed microbiome where the identity of labeled genomes, and their level of isotopic enrichment, were experimentally controlled. With this ground truth dataset, we empirically assessed the accuracy of different analytic models for identifying active taxa, and examined how sequencing depth impacts the detection of isotopically labeled genomes. We also demonstrate that using synthetic DNA internal standards to measure absolute genome abundances in SIP density fractions improves estimates of isotopic enrichment. In addition, our study illustrates the utility of internal standards to reveal anomalies in sample handling that could negatively impact SIP metagenomic analyses if left undetected. Finally, we present SIPmg, an R package to facilitate the estimation of absolute abundances and perform statistical analyses for identifying labeled genomes within SIP metagenomic data. This experimentally validated analysis framework strengthens the foundation of DNA-SIP metagenomics as a tool for accurately measuring the in situ activity of environmental microbial populations and assessing their genomic potential.

No abstract