Advancing environmental biotechnology with microbial community modeling rooted in functional 'omics

Emerging biotechnologies that solve pressing environmental and climate emergencies will require harnessing the vast functional diversity of the underlying microbiomes driving such engineered processes. Modeling is a critical aspect of process engineering that informs system design as well as aids diagnostic optimization of performance. 'Conventional' bioprocess models assume homogenous biomass within functional guilds and thus fail to predict emergent properties of diverse microbial physiologies, such as product specificity and community interactions. Yet, recent advances in functional 'omics-based approaches can provide a 'lens' through which we can probe and measure in situ ecophysiologies of environmental microbiomes. Here, we overview microbial community modeling approaches that incorporate functional 'omics data, which we posit can advance our ability to design and control new environmental biotechnologies going forward.

Metatranscriptomics-guided genome-scale metabolic reconstruction reveals the carbon flux and trophic interaction in methanogenic communities

BACKGROUND

Despite rapid advances in genomic-resolved metagenomics and remarkable explosion of metagenome-assembled genomes (MAGs), the function of uncultivated anaerobic lineages and their interactions in carbon mineralization remain largely uncertain, which has profound implications in biotechnology and biogeochemistry.

RESULTS

In this study, we combined long-read sequencing and metatranscriptomics-guided metabolic reconstruction to provide a genome-wide perspective of carbon mineralization flow from polymers to methane in an anaerobic bioreactor. Our results showed that incorporating long reads resulted in a substantial improvement in the quality of metagenomic assemblies, enabling the effective recovery of 132 high-quality genomes meeting stringent criteria of minimum information about a metagenome-assembled genome (MIMAG). In addition, hybrid assembly obtained 51% more prokaryotic genes in comparison to the short-read-only assembly. Metatranscriptomics-guided metabolic reconstruction unveiled the remarkable metabolic flexibility of several novel Bacteroidales-affiliated bacteria and populations from Mesotoga sp. in scavenging amino acids and sugars. In addition to recovering two circular genomes of previously known but fragmented syntrophic bacteria, two newly identified bacteria within Syntrophales were found to be highly engaged in fatty acid oxidation through syntrophic relationships with dominant methanogens Methanoregulaceae bin.74 and Methanothrix sp. bin.206. The activity of bin.206 preferring acetate as substrate exceeded that of bin.74 with increasing loading, reinforcing the substrate determinantal role.

CONCLUSION

Overall, our study uncovered some key active anaerobic lineages and their metabolic functions in this complex anaerobic ecosystem, offering a framework for understanding carbon transformations in anaerobic digestion. These findings advance the understanding of metabolic activities and trophic interactions between anaerobic guilds, providing foundational insights into carbon flux within both engineered and natural ecosystems. Video Abstract.

Quantifying functional redundancy in polysaccharide-degrading prokaryotic communities

BACKGROUND

Functional redundancy (FR) is widely present, but there is no consensus on its formation process and influencing factors. Taxonomically distinct microorganisms possessing genes for the same function in a community lead to within-community FR, and distinct assemblies of microorganisms in different communities playing the same functional roles are termed between-community FR. We proposed two formulas to respectively quantify the degree of functional redundancy within and between communities and analyzed the FR degrees of carbohydrate degradation functions in global environment samples using the genetic information of glycoside hydrolases (GHs) encoded by prokaryotes.

RESULTS

Our results revealed that GHs are each encoded by multiple taxonomically distinct prokaryotes within a community, and the enzyme-encoding prokaryotes are further distinct between almost any community pairs. The within- and between-FR degrees are primarily affected by the alpha and beta community diversities, respectively, and are also affected by environmental factors (e.g., pH, temperature, and salinity). The FR degree of the prokaryotic community is determined by deterministic factors.

CONCLUSIONS

We conclude that the functional redundancy of GHs is a stabilized community characteristic. This study helps to determine the FR formation process and influencing factors and provides new insights into the relationships between prokaryotic community biodiversity and ecosystem functions. Video Abstract.

Enzyme modified biodegradable plastic preparation and performance in anaerobic co-digestion with food waste

A modified biodegradable plastic (PLA/PBAT) was developed by through covalent bonding with proteinase K, porcine pancreatic lipase, or amylase, and was then investigated in anaerobic co-digestion mixed with food waste. Fluorescence microscope validated that enzymes could remain stable in modified the plastic, even after co-digestion. The results of thermophilic anaerobic co-digestion showed that, degradation of the plastic modified with Proteinase K increased from 5.21 ± 0.63 % to 29.70 ± 1.86 % within 30 days compare to blank. Additionally, it was observed that the cumulative methane production increased from 240.9 ± 0.5 to 265.4 ± 1.8 mL/gVS, and the methane production cycle was shortened from 24 to 20 days. Interestingly, the kinetic model suggested that the modified the plastic promoted the overall hydrolysis progression of anaerobic co-digestion, possibly as a result of the enhanced activities of Bacteroidota and Thermotogota. In conclusion, under anaerobic co-digestion, the modified the plastic not only achieved effective degradation but also facilitated the co-digestion process.

New insights on ecological roles of waste activated sludge in nutrient-stressed co-digestion

There have been extensive applications of waste activated sludge (WAS) in anaerobic co-digestion (AcoD). Nonetheless, mechanisms through which AcoD systems maintain stability, particularly under nutrient-stressed conditions, are under-appreciated. In this study, the role of WAS in a nutrient-stressed WAS-food waste AcoD system was re-evaluated. Our findings demonstrated that WAS-based co-digestion increased methane production (by 20-60%) as WAS bolsters such systems' resilience via establishing a core niche-based microbial balance. The carbon utilization investigation suggested a microbial niche balance is attainable if two conditions are satisfied: 1) hydrolysis efficiency is greater than 50%; and 2) both the acidogenesis-to-hydrolysis and acetogenesis-to-hydrolysis efficiencies surpass 0.5. Metagenomic assembly genome (MAG) analysis indicated that the versatile metabolic characteristics strengthened the microbial niche balance, rendering the system resilient and efficient through a syntrophic mode, contributing to both acidogenesis and acetogenesis. The findings of this study provide new insights into the ecological effects of WAS on AcoD.

Microbial community structure in a constructed wetland based on a recirculating aquaculture system: Exploring spatio-temporal variations and assembly mechanisms

The diversity, composition and performance of microbial communities within constructed wetlands (CW) were markedly influenced by spatio-temporal variations. A pilot-scale integrated vertical-flow constructed wetland (IVCW) as the biological purification unit within a recirculating aquaculture system (RAS) was established and monitored in this study. The investigation aimed to elucidate the responses of community structure, co-occurrence networks, and assembly mechanisms of the microbial community to spatial and temporal changes. Spatially, all a-diversity indices and microbial networks complexity were significantly higher in the upstream pool of the IVCW than in the downstream pool. Temporally, the richness increased over time, while the evenness showed a decreasing trend. The number of nodes and edges of microbial networks increased over time. Notably, the stable pollutant removal efficiencies were observed during IVCW operations, despite a-diversity and bacterial community networks exhibited significant variations across time. Functional redundancy emerged as a likely mechanism contributing to the stability of microbial ecosystem functions. Null model and neutral model analyses revealed the dominance of deterministic processes shaping microbial communities over time, with deterministic influences being more pronounced at lower a-diversity levels. DO and inorganic nitrogen emerged as the principal environmental factor influencing microbial community dynamics. This study provides a theoretical foundation for the regulation of microbial communities and environmental factors within the context of IVCW.

Effects of operational parameters on bacterial communities in Hong Kong and global wastewater treatment plants

Wastewater treatment plants (WWTPs) are indispensable biotechnology facilities for modern cities and play an essential role in modern urban infrastructure by employing microorganisms to remove pollutants in wastewater, thus protecting public health and the environment. This study conducted a 13-month bacterial community survey of six full-scale WWTPs in Hong Kong with samples of influent, activated sludge (AS), and effluent to explore their synchronism and asynchronism of bacterial community. Besides, we compared AS results of six Hong Kong WWTPs with data from 1,186 AS amplicon data in 269 global WWTPs and a 9-year metagenomic sequencing survey of a Hong Kong WWTP. Our results showed the compositions of bacterial communities varied and the bacterial community structure of AS had obvious differences across Hong Kong WWTPs. The co-occurrence analysis identified 40 pairs of relationships that existed among Hong Kong WWTPs to show solid associations between two species and stochastic processes took large proportions for the bacterial community assembly of six WWTPs. The abundance and distribution of the functional bacteria in worldwide and Hong Kong WWTPs were examined and compared, and we found that ammonia-oxidizing bacteria had more diversity than nitrite-oxidizing bacteria. Besides, Hong Kong WWTPs could make great contributions to the genome mining of microbial dark matter in the global "wanted list." Operational parameters had important effects on OTUs' abundance, such as the temperature to the genera of Tetrasphaera, Gordonia and Nitrospira. All these results obtained from this study can deepen our understanding of the microbial ecology in WWTPs and provide foundations for further studies.

IMPORTANCE

Wastewater treatment plants (WWTPs) are an indispensable component of modern cities, as they can remove pollutants in wastewater to prevent anthropogenic activities. Activated sludge (AS) is a fundamental wastewater treatment process and it harbors a highly complex microbial community that forms the main components and contains functional groups. Unveiling "who is there" is a long-term goal of the research on AS microbiology. High-throughput sequencing provides insights into the inventory diversity of microbial communities to an unprecedented level of detail. At present, the analysis of communities in WWTPs usually comes from a specific WWTP and lacks comparisons and verification among different WWTPs. The wide-scale and long-term sampling project and research in this study could help us evaluate the AS community more accurately to find the similarities and different results for different WWTPs in Hong Kong and other regions of the world.

Metabolism of novel potential syntrophic acetate-oxidizing bacteria in thermophilic methanogenic chemostats

Acetate is a major intermediate in the anaerobic digestion of organic waste to produce CH4. In methanogenic systems, acetate degradation is carried out by either acetoclastic methanogenesis or syntrophic degradation by acetate oxidizers and hydrogenotrophic methanogens. Due to challenges in the isolation of syntrophic acetate-oxidizing bacteria (SAOB), the diversity and metabolism of SAOB and the mechanisms of their interactions with methanogenic partners are not fully characterized. In this study, the in situ activity and metabolic characteristics of potential SAOB and their interactions with methanogens were elucidated through metagenomics and metatranscriptomics. In addition to the reported SAOB classified in the genera Tepidanaerobacter, Desulfotomaculum, and Thermodesulfovibrio, we identified a number of potential SAOB that are affiliated with Clostridia, Thermoanaerobacteraceae, Anaerolineae, and Gemmatimonadetes. The potential SAOB possessing the glycine-mediated acetate oxidation pathway dominates SAOB communities. Moreover, formate appeared to be the main product of the acetate degradation by the most active potential SAOB. We identified the methanogen partner of these potential SAOB in the acetate-fed chemostat as Methanosarcina thermophila. The dominated potential SAOB in each chemostat had similar metabolic characteristics, even though they were in different fatty-acid-fed chemostats. These novel syntrophic lineages are prevalent and may play critical roles in thermophilic methanogenic reactors. This study expands our understanding of the phylogenetic diversity and in situ biological functions of uncultured syntrophic acetate degraders and presents novel insights into how they interact with methanogens.IMPORTANCECombining reactor operation with omics provides insights into novel uncultured syntrophic acetate degraders and how they perform in thermophilic anaerobic digesters. This improves our understanding of syntrophic acetate degradation and contributes to the background knowledge necessary to better control and optimize anaerobic digestion processes.

Core taxa, co-occurrence pattern, diversity, and metabolic pathways contributing to robust anaerobic biodegradation of chlorophenol

It is hard to achieve robustness in anaerobic biodegradation of trichlorophenol (TCP). We hypothesized that specific combinations of environmental factors determine phylogenetic diversity and play important roles in the decomposition and stability of TCP-biodegrading bacteria. The anaerobic bioreactor was operated at 35 °C (H condition) or 30 °C (L condition) and mainly fed with TCP (from 28 μM to 180 μM) and organic material. Metagenome sequencing was combined with 16S rRNA gene amplicon sequencing for the microbial community analysis. The results exhibited that the property of robustness occurred in specific conditions. The corresponding co-occurrence and diversity patterns suggest high collectivization, degree and evenness for robust communities. Two types of core functional taxa were recognized: dechlorinators (unclassified Anaerolineae, Thermanaerothrix and Desulfovibrio) and ring-opening members (unclassified Proteobacteria, Methanosarcina, Methanoperedens, and Rubrobacter). The deterministic process of the expansion of niche of syntrophic bacteria at higher temperatures was confirmed. The reductive and hydrolytic dechlorination mechanisms jointly lead to C-Cl bond cleavage. H ultimately adapted to the stress of high TCP loading, with more abundant ring-opening enzyme (EC 3.1.1.45, ∼55%) and hydrolytic dechlorinase (EC 3.8.1.5, 26.5%) genes than L (∼47%, 10.5%). The functional structure (based on KEGG) in H was highly stable despite the high loading of TCP (up to 60 μM), but not in L. Furthermore, an unknown taxon with multiple functions (dechlorinating and ring-opening) was found based on genetic sequencing; its functional contribution of EC 3.8.1.5 in H (26.5%) was higher than that in L (10.5%), and it possessed a new metabolic pathway for biodegradation of halogenated aromatic compounds. This new finding is supplementary to the robust mechanisms underlying organic chlorine biodegradation, which can be used to support the engineering, regulation, and design of synthetic microbiomes.

Riboflavin-loaded carbon cloth aids the anaerobic digestion of cow dung by promoting direct interspecies electron transfer

Cow dung generates globally due to increased beef and milk consumption, but its treatment efficiency remains low. Previous studies have shown that riboflavin-loaded conductive materials can improve anaerobic digestion through enhance direct interspecies electron transfer (DIET). However, its effect on the practical anaerobic digestion of cow dung remained unclear. In this study, carbon cloth loaded with riboflavin (carbon cloth-riboflavin) was added into an anaerobic digester treating cow dung. The carbon cloth-riboflavin reactor showed a better performance than other two reactors. The metagenomic analysis revealed that Methanothrix on the surface of the carbon cloth predominantly utilized the CO2 reduction for methane production, further enhanced after riboflavin addition, while Methanothrix in bulk sludge were using the acetate decarboxylation pathway. Furthermore, the carbon cloth-riboflavin enriched various major methanogenic pathways and activated a large number of enzymes associated with DIET. Riboflavin's presence altered the microbial communities and the abundance of functional genes relate to DIET, ultimately leading to a better performance of anaerobic digestion for cow dung.

Analysis on methane production from various coal slime fermentations based on metagenomics

Metagenomic sequencing technology was applied to evaluate differences in the anaerobic fermentation process of coal slimes by analyzing microbial diversity, functional activity structure, and cooperative relationship during the anaerobic fermentation of coal slimes with different coal ranks. The obtained results showed that the production of biomethane from coal slime was decreased by increasing metamorphism degree. Internal reason was higher abundance of microbial community in low rank coal slimes compared to that in high rank coal which had higher activity in the gene expression of key steps such as hydrolysis and acidification, methanation and the production of hydrogen and acetic acid. Acetic acid decarboxylation and CO₂ reduction are two key pathways of methanation process. At the same time, K11261 (formylmethanofuran dehydrogenase subunit) and K01499 (methenyltetrahydromethanopterin cyclohydrolase) genes were further enriched in low rank slime systems, which enhanced the proportion of CO₂ reduction in methanation pathway and was beneficial to biomethane production. Research revealed the roles of different coal slime ranks in biomethane production process and is considered as an important reference significance for further exploration of coal slime resource utilization.

Potential effect of key soil bacterial taxa on the increase of rice yield under milk vetch rotation

Legume crop rotation is often adopted in rice cultivation to improve soil productivity. However, little is known about the role of microbes under legume rotation in affecting soil productivity. To elucidate this, a long-term paddy cropping experiment was set up to study the relationship between crop yield, soil chemical properties, and key microbial taxa under a double-rice and milk vetch rotation. Milk vetch rotation significantly improved soil chemical properties compared to no fertilization treatment, and soil phosphorus was a major factor correlated with crop yield. Long-term legume rotation increased soil bacterial alpha diversity and changed soil bacterial community. After milk vetch rotation, the relative abundances of Bacteroidota, Desulfobacterota, Firmicutes, and Proteobacteria increased while those of Acidobacteriota, Chloroflexi, and Planctomycetota decreased. Moreover, milk vetch rotation increased the relative abundance of phosphorus-related gene K01083 (bpp), which was significantly correlated with soil phosphorus content and crop yield. Network analysis showed that taxa of Vicinamibacterales were positively correlated with total phosphorus and available phosphorus, which was a potential taxon contributing to the availability of soil phosphorus stock. Our results indicated that milk vetch rotation could enrich key taxa with latent phosphate-solubilizing ability, increase the content of soil available phosphorus, and finally enhance crop yield. This could provide scientific guidance for better crop production.

Genome-centric metagenomics revealed functional traits in high-solids anaerobic co-digestion of restaurant food waste, household food waste and rice straw

High-solids anaerobic co-digestion (HS-AcoD) of food waste (FW) and other organic wastes is an effective option to improve the biogas production and system stability compared to mono-digestion. However, the clean and sustainable HS-AcoD strategy for FW and associated microbial functional traits have not been well explored. Here, HS-AcoD of restaurant food waste (RFW), household food waste (HFW) and rice straw (RS) were performed. Results showed that the maximum synergy index (SI) of 1.28 were achieved when the volatile solids ratio of RFW, HFW and RS was 0.45:0.45:0.1. HS-AcoD alleviated the acidification process by regulating metabolism associated with hydrolysis and volatile fatty acids formation. The synergistic relationship between syntrophic bacteria and Methanothrix sp., and the enhanced metabolic capacity associated with the acetotrophic and hydrogenotrophic pathways dominated by Methanothrix sp., provided a further explanation of the synergistic mechanism. These findings advance the knowledge about microbial mechanisms underlying the synergistic effect of HS-AcoD.

Sustainability for wastewater treatment: bioelectricity generation and emission reduction

This review covers the technological measures of a self-sustainable anaerobic up-flow sludge blanket (UASB) system compared with an aerobic activated sludge process (ASP) for wastewater treatment plants (WWTPs). The ASP requires a huge amount of electricity and chemicals and also results in the emission of carbon. The UASB system, instead, is based on greenhouse gas (GHG) emission reduction and is associated with biogas production for cleaner electricity. WWTPs including the ASP system are not sustainable due to the massive financial power required for clean wastewater. When the ASP system was used, the amount of production was estimated to be 10658.98 tonnes CO₂eq-d^- of carbon dioxide. Whereas it was 239.19 tonnes CO₂eq-d^-1 with the UASB. The UASB system is advantageous over the ASP system as it has a high production of biogas, needs low maintenance, yields a low amount of sludge, and is also a source of electricity that can be used as a power source for the WWTPs. Also, the UASB system produces less biomass, and this helps in reducing costs and maintaining work. Moreover, the aeration tank of the ASP needs 60% of energy distribution; on the other hand, the UASB consumes less energy, approximately 3-11%.

Trait biases in microbial reference genomes

Common culturing techniques and priorities bias our discovery towards specific traits that may not be representative of microbial diversity in nature. So far, these biases have not been systematically examined. To address this gap, here we use 116,884 publicly available metagenome-assembled genomes (MAGs, completeness ≥80%) from 203 surveys worldwide as a culture-independent sample of bacterial and archaeal diversity, and compare these MAGs to the popular RefSeq genome database, which heavily relies on cultures. We compare the distribution of 12,454 KEGG gene orthologs (used as trait proxies) in the MAGs and RefSeq genomes, while controlling for environment type (ocean, soil, lake, bioreactor, human, and other animals). Using statistical modeling, we then determine the conditional probabilities that a species is represented in RefSeq depending on its genetic repertoire. We find that the majority of examined genes are significantly biased for or against in RefSeq. Our systematic estimates of gene prevalences across bacteria and archaea in nature and gene-specific biases in reference genomes constitutes a resource for addressing these issues in the future.

Spatiotemporal Dynamics of Bacterial Taxonomic and Functional Profiles in Estuarine Intertidal Soils of China Coastal Zone

Bacteria play an important role in regulating carbon (C), nitrogen (N), and sulfur (S) in estuarine intertidal wetlands. To gain insights into the ecological and metabolic modes possessed by bacteria in estuarine intertidal wetlands, a total of 78 surface soil samples were collected from China's coastal intertidal wetlands to examine the spatial and seasonal variations of bacterial taxonomic composition, assembly processes, and ecological system functions through shotgun metagenomic and 16S rRNA gene sequencing. Obvious spatiotemporal dynamic patterns in the bacterial community structure were identified, with more pronounced seasonal rather than spatial variations. Dispersion limitation was observed to act as a critical factor affecting community assembly, explaining approximately half of the total variation in the bacterial community. Functional bacterial community structure exhibited a more significant latitudinal change than seasonal variability, highlighting that functional stability of the bacterial communities differed with their taxonomic variability. Identification of biogeochemically related links between C, N, and S cycles in the soils showed the adaptive routed metabolism of the bacterial communities and the strong interactions between coupled metabolic pathways. Our study broadens the insights into the taxonomic and functional profiles of bacteria in China's estuarine intertidal soils and helps us understand the effects exerted by environmental factors on the ecological health and microbial diversity of estuarine intertidal flats.

Comparison of prokaryotes between Mount Everest and the Mariana Trench

BACKGROUND

Mount Everest and the Mariana Trench represent the highest and deepest places on Earth, respectively. They are geographically separated, with distinct extreme environmental parameters that provide unique habitats for prokaryotes. Comparison of prokaryotes between Mount Everest and the Mariana Trench will provide a unique perspective to understanding the composition and distribution of environmental microbiomes on Earth.

RESULTS

Here, we compared prokaryotic communities between Mount Everest and the Mariana Trench based on shotgun metagenomic analysis. Analyzing 25 metagenomes and 1176 metagenome-assembled genomes showed distinct taxonomic compositions between Mount Everest and the Mariana Trench, with little taxa overlap, and significant differences in genome size, GC content, and predicted optimal growth temperature. However, community metabolic capabilities exhibited striking commonality, with > 90% of metabolic modules overlapping among samples of Mount Everest and the Mariana Trench, with the only exception for CO2 fixations (photoautotrophy in Mount Everest but chemoautotrophy in the Mariana Trench). Most metabolic pathways were common but performed by distinct taxa in the two extreme habitats, even including some specialized metabolic pathways, such as the versatile degradation of various refractory organic matters, heavy metal metabolism (e.g., As and Se), stress resistance, and antioxidation. The metabolic commonality indicated the overall consistent roles of prokaryotes in elemental cycling and common adaptation strategies to overcome the distinct stress conditions despite the intuitively huge differences in Mount Everest and the Mariana Trench.

CONCLUSION

Our results, the first comparison between prokaryotes in the highest and the deepest habitats on Earth, may highlight the principles of prokaryotic diversity: although taxa are habitat-specific, primary metabolic functions could be always conserved. Video abstract.

Interactions between Bacterial Inoculants and Native Soil Bacterial Community: the Case of Spore-forming Bacillus spp

Microbial diversity can restrict the invasion and impact of alien microbes into soils via resource competition. However, this theory has not been tested on various microbial invaders with different ecological traits, particularly spore-forming bacteria. Here we investigated the survival capacity of two introduced spore-forming bacteria, Bacillus mycoides (BM) and B. pumillus (BP) and their impact on the soil microbiome niches with low and high diversity. We hypothesized that higher soil bacterial diversity would better restrict Bacillus survival via resource competition, and the invasion would alter the resident bacterial communities' niches only if inoculants do not escape competition with the soil community (e.g. through sporulation). Our findings showed that BP could not survive as viable propagules and transiently impacted the bacterial communities' niche structure. This may be linked to its poor resource usage and low growth rate. Having better resource use capacities, BM better survived in soil, though its survival was weakly related to the remaining resources left for them by the soil community. BM strongly affected the community niche structure, ultimately in less diverse communities. These findings show that the inverse diversity-invasibility relationship can be valid for some spore-forming bacteria, but only when they have sufficient resource use capacity.

Absolute quantification and genome-centric analyses elucidate the dynamics of microbial populations in anaerobic digesters

Anaerobic digestion (AD) relies on myriads of functions performed by complex microbial communities in customized settings, thus, a comprehensive investigation on the AD microbiome is central to the fine-tuned control. Most current AD microbiome studies are based on relative abundance, which hinders the interpretation of microbes' dynamics and inter-sample comparisons. Here, we developed an absolute quantification (AQ) approach that integrated cellular spike-ins with metagenomic sequencing to elucidate microbial community variations and population dynamics in four anaerobic digesters. Using this method, 253 microbes were defined as decaying populations with decay rates ranging from -0.05 to -5.85 d^-1, wherein, a population from Flavobacteriaceae family decayed at the highest rates of -3.87 to -5.85 d^-1 in four digesters. Meanwhile, 25 microbes demonstrated the growing trend in the AD processes with growth rates ranging from 0.11 to 1.77 d^-1, and genome-centric analysis assigned some of the populations to the functional niches of hydrolysis, short-chain fatty acids metabolism, and methane generation. Additionally, we observed that the specific activity of methanogens was lower in the prolonged digestion stage, and redundancy analysis revealed that the feedstock composition and the digestion duration were the two key parameters in governing the AD microbial compositions.

Integrating Genome-Resolved Metagenomics with Trait-Based Process Modeling to Determine Biokinetics of Distinct Nitrifying Communities within Activated Sludge

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.

Temperature mediated the balance between stochastic and deterministic processes and reoccurrence of microbial community during treating aniline wastewater

Seasonal temperature changes significantly affect microbial community diversity, composition, and performance in wastewater treatment plants. However, the community assembly mechanisms under seasonal temperature variations remain unclear. Here, we carried out temperature cycling experiments (30 °C, 35 °C, 37 °C, 40 °C, 42 °C, 45 °C, 40 °C, and 30 °C) to investigate how temperature impacts microbial performance and co-occurrence network and how assembly processes determine the structure and function of microbial communities during treating aniline wastewater. During the 195-day operation, the system achieved an efficient and stable aniline removal of 99%. Interestingly, α-diversity and network complexity were negatively correlated with temperature but could be recovered when the temperature was returned to 30 °C. The results showed that functional redundancy was probably responsible for the excellent microbial performance during the whole process. Null model analyses presented that deterministic process dominated the community when the temperature was 30 °C, and stochasticity dominated the assembly process when the temperature was over 30 °C. Overall, the balance between stochastic and deterministic processes in the treatment of aniline wastewater mediated the reoccurrence of microbial community and co-occurrence network at different temperatures. This study provides new insights into microbial community reoccurrence under seasonal temperature changes and a theoretical basis for regulating microbial communities in wastewater treatment plants.

Impacts of food waste to sludge ratios on microbial dynamics and functional traits in thermophilic digesters

A self-stabilizing microbial community lays the foundation of the efficient biochemical reactions of the anaerobic digestion (AD) process. Despite extensive profiling of microbial community dynamics under varying operating parameters, the effects of food waste (FW) to feeding sewage sludge (FSS) ratios on the microbial assembly, functional traits, and syntrophic interspecies interactions in thermophilic microbial consortia remain poorly understood. Here, we investigated the long-term impacts of the FW: FSS ratio on the thermophilic AD microbiome using genome-centric metagenomics. Both the short reads (SRs) assembly, and the iterative hybrid assembly (IHA) of SRs and nanopore long reads (LRs) were used to reconstruct metagenome-assembled genomes (MAGs) and four microbial clusters were identified, demonstrating different microbial dynamics patterns in response to varying FW:FSS ratios. Cluster C1-C3 were comprised of full functional members with genetic potentials in fulfilling empirical AD biochemical reactions, wherein, syntrophic decarboxylating acetogens could interact with methanogens, and some microbes could be energized by the electron bifurcation mechanism to drive thermodynamics unfavorable reactions. We found the co-existence of both acetogenic and hydrogenotrophic methanogens in the AD microbiome, and they altered their trophic groups to scavenge the methanogenic substrates in ensuring the methane generation in digesters with different FW:FSS ratios. Another interesting observation was that two phylogenetically close Thermotogota species showed a possible strong competition on carbon source inferred by the nearly complete genetic overlap of their relevant pathways.

Impact of substrate digestibility on microbial community stability in methanogenic digestors: The mechanism and solution

This study investigated the temporal dynamics of digestion efficiency and community stability in digesters fed with waste activated sludge (WAS), straw (STR-AD), food waste (FW-AD) and mixture of straw-and-food waste (STR-FW-AD). Results showed that carbon removals of recalcitrant substrates (i.e., 48.2 ± 3.9% in WAS-AD and 57.8 ± 4.9% in STR-AD) were lower than that of labile substrates (i.e., 70.7 ± 4.0% in FW-AD). Nonetheless, carbon removal of recalcitrant substrates was largely improved through co-digestion (70.3 ± 3.2% in STR-FW-AD). In contrast to monopoly communities (e.g., the highly enriched Paludibacter) fed with the labile substrates, recalcitrant substrates supported highly diverse communities. Accordingly, the medians of negative/positive cohesions of communities in WAS-AD, STR-AD, STR-FW-AD and FW-AD decreased from 0.86 to 0.63, suggesting their decreasing community stability. Microbial source tracking analyses showed the major contribution of the STR-AD community to the co-digestion community. This study provided unprecedented mechanistic insight into stability improvement of substrate co-digestion on the methanogenic digestion microbiome.

The Coupling Response between Different Bacterial Metabolic Functions in Water and Sediment Improve the Ability to Mitigate Climate Change

Extreme climatic events, such as heat wave and large temperature fluctuations, are predicted to increase in frequency and intensity during the next hundred years, which may rapidly alter the composition and function of lake bacterial communities. Here, we conducted a year-long experiment to explore the effect of warming on bacterial metabolic function of lake water and sediment. Predictions of the metabolic capabilities of these communities were performed with FAPROTAX using 16S rRNA sequencing data. The results indicated that the increase in temperature changed the structure of bacterial metabolic functional groups in water and sediment. During periods of low temperature, the carbon degradation pathway decreased, and the synthesis pathway increased, under the stimulation of warming, especially under the conditions temperature fluctuation. We also observed that nitrogen fixation ability was especially important in the warming treatments during the summer season. However, an elevated temperature significantly led to reduced nitrogen fixation abilities in winter. Compared with the water column, the most predominant functional groups of nitrogen cycle in sediment were nitrite oxidation and nitrification. Variable warming significantly promoted nitrite oxidation and nitrification function in winter, and constant warming was significantly inhibited in spring, with control in sediments. Co-occurrence network results showed that warming, especially variable warming, made microbial co-occurrence networks larger, more connected and less modular, and eventually functional groups in the water column and sediment cooperated to resist warming. We concluded that warming changed bacterial functional potentials important to the biogeochemical cycling in the experimental mesocosms in winter and spring with low temperature. The effect of different bacteria metabolism functions in water column and sediment may change the carbon and nitrogen fluxes in aquatic ecosystems. In conclusion, the coupling response between different bacterial metabolic functions in water and sediment may improve the ability to mitigate climate change.

Current understanding and perspectives in anaerobic digestion based on genome-resolved metagenomic approaches

Anaerobic digestion (AD) is a technique that can be used to treat high concentrations of various organic wastes using a consortium of functionally diverse microorganisms under anaerobic conditions. Methane gas, a beneficial by-product of the AD process, is a renewable energy source that can replace fossil fuels following purification. However, detailed functional roles and metabolic interactions between microbial populations involved in organic waste removal and methanogenesis are yet to be known. Recent metagenomic approaches based on advanced high-throughput sequencing techniques have enabled the exploration of holistic microbial taxonomy and functionality of complex microbial populations involved in the AD process. Gene-centric and genome-centric analyses based on metagenome-assembled genomes are a platform that can be used to study the composition of microbial communities and their roles during AD. This review looks at how these up-to-date metagenomic analyses can be applied to promote our understanding and improved the development of the AD process.

Extracellular hydrolytic potential drives microbiome shifts during anaerobic co-digestion of sewage sludge and food waste

Bacterial community structure and dynamics in anaerobic digesters are primarily influenced by feedstock composition. It is therefore important to unveil microbial traits that explain microbiome variations in response to substrate changes. Here, gene and genome-centric metagenomics were used to examine microbiome dynamics in four laboratory-scale reactors, in which sewage sludge was co-digested with increasing amounts of food waste. A co-occurrence network revealed microbiome shifts in response to changes in substrate composition and concentration. Food waste concentration correlated with extracellular enzymes and metagenome-assembled genomes (MAGs) involved in the degradation of complex carbohydrates commonly found in fruits and plant cell walls as well as with the abundance of hydrolytic MAGs. A key role was attributed to Proteiniphillum for being the only bacteria that encoded the complete pectin degradation pathway. These results suggest that changes of feedstock composition establish new microbial niches for bacteria with the capacity to degrade newly added substrates.

Convergent Microbial Community Formation in Replicate Anaerobic Reactors Inoculated from Different Sources and Treating Ersatz Crew Waste

Future manned space travel will require efficient recycling of nutrients from organic waste back into food production. Microbial systems are a low-energy, efficient means of nutrient recycling, but their use in a life support system requires predictability and reproducibility in community formation and reactor performance. To assess the reproducibility of microbial community formation in fixed-film reactors, we inoculated replicate anaerobic reactors from two methanogenic inocula: a lab-scale fixed-film, plug-flow anaerobic reactor and an acidic transitional fen. Reactors were operated under identical conditions, and we assessed reactor performance and used 16s rDNA amplicon sequencing to determine microbial community formation. Reactor microbial communities were dominated by similar groups, but differences in community membership persisted in reactors inoculated from different sources. Reactor performance overlapped, suggesting a convergence of both reactor communities and organic matter mineralization. The results of this study suggest an optimized microbial community could be preserved and used to start new, or restart failed, anaerobic reactors in a life support system with predictable reactor performance.

Metagenomic insights into the effect of thermal hydrolysis pre-treatment on microbial community of an anaerobic digestion system

Thermal hydrolysis process (THP) is an effective pre-treatment method to reduce solids volume and improve biogas production during anaerobic digestion (AD) via increasing the biodegradability of waste activated sludge (WAS). However, the effects of THP pre-treated sludge on microbial diversity, interspecies interactions, and metabolism in AD systems remain largely unknown. We therefore setup and operated an anaerobic digester during a long-term period to shed light on the effect of THP pre-treatment on AD microbial ecology in comparison to conventional AD via Illumina based 16S rRNA gene amplicon sequencing and genome-centric metagenomics analysis. Results showed THP sludge significantly reduced the microbial diversity, shaped the microbial community structure, and resulted in more intense microbial interactions. Compared to WAS as the feed sludge, THP sludge shaped the core functional groups, but functional redundancy ensured the system's stability. The metabolic interactions between methanogens and syntrophic bacteria as well as the specific metabolic pathways were further elucidated. Hydrogenotrophic methanogens, Methanospirillum sp. and Methanolinea sp., were the primary contributors for methane production when treating THP and WAS, respectively, which also have potential for acetate oxidation to methane. Collectively, this study provides in-depth information on the interspecies interactions to better understand how THP pre-treatment influences AD microbial community.

Elevated pCO2 alters the interaction patterns and functional potentials of rearing seawater microbiota

Mean oceanic CO₂ values have already risen and are expected to rise further on a global scale. Elevated pCO₂ (eCO₂) changes the bacterial community in seawater. However, the ecological association of seawater microbiota and related geochemical functions are largely unknown. We provide the first evidence that eCO₂ alters the interaction patterns and functional potentials of microbiota in rearing seawater of the swimming crab, Portunus trituberculatus. Network analysis showed that eCO₂ induced a simpler and more modular bacterial network in rearing seawater, with increased negative associations and distinct keystone taxa. Using the quantitative microbial element cycling method, nitrogen (N) and phosphorus (P) cycling genes exhibited the highest increase after one week of eCO₂ stress and were significantly associated with keystone taxa. However, the functional potential of seawater bacteria was decoupled from their taxonomic composition and strongly coupled with eCO₂ levels. The changed functional potential of seawater bacteria contributed to seawater N and P chemistry, which was highlighted by markedly decreased NH₃, NH₄⁺-N, and PO₄^3--P levels and increased NO₂^--N and NO₃^--N levels. This study suggests that eCO₂ alters the interaction patterns and functional potentials of seawater microbiota, which lead to the changes of seawater chemical parameters. Our findings provide new insights into the mechanisms underlying the effects of eCO₂ on marine animals from the microbial ecological perspective.

Polysaccharide-Bacteria Interactions From the Lens of Evolutionary Ecology

Microbes have the unique ability to break down the complex polysaccharides that make up the bulk of organic matter, initiating a cascade of events that leads to their recycling. Traditionally, the rate of organic matter degradation is perceived to be limited by the chemical and physical structure of polymers. Recent advances in microbial ecology, however, suggest that polysaccharide persistence can result from non-linear growth dynamics created by the coexistence of alternate degradation strategies, metabolic roles as well as by ecological interactions between microbes. This complex "landscape" of degradation strategies and interspecific interactions present in natural microbial communities appears to be far from evolutionarily stable, as frequent gene gain and loss reshape enzymatic repertoires and metabolic roles. In this perspective, we discuss six challenges at the heart of this problem, ranging from the evolution of genetic repertoires, phenotypic heterogeneity in clonal populations, the development of a trait-based ecology, and the impact of metabolic interactions and microbial cooperation on degradation rates. We aim to reframe some of the key questions in the study of polysaccharide-bacteria interactions in the context of eco-evolutionary dynamics, highlighting possible research directions that, if pursued, would advance our understanding of polysaccharide degraders at the interface between biochemistry, ecology and evolution.

Linkage of community composition and function over short response time in anaerobic digestion systems with food fermentation wastewater

We investigated the short-term dynamics of microbial composition and function in bioreactors with inocula collected from full-scale and laboratory-based anaerobic digestion (AD) systems. The Bray-Curtis dissimilarity of both inocula was approximately 10% of the predicted Kyoto Encyclopedia of Genes and Genomes pathway and 40% of the taxonomic composition and yet resulted in a similar performance in methane production, implying that the variation of community composition may be decoupled from performance. However, the significant correlation of volatile fatty acids with taxonomic variation suggested that the pathways of AD could be different because of the varying genus. The predicted function of the significantly varying genus was mostly related to fermentation, which strengthened the conclusion that most microbial variation occurred within the fermentative species and led to alternative routes to result in similar methane production in methanogenic bioreactors. This finding sheds some light on the understanding of AD community regulation, which depends on the aims to recover intermediates or methane.

Trophic Interactions and the Drivers of Microbial Community Assembly

Despite numerous surveys of gene and species content in heterotrophic microbial communities, such as those found in animal guts, oceans, or soils, it is still unclear whether there are generalizable biological or ecological processes that control their dynamics and function. Here, we review experimental and theoretical advances to argue that networks of trophic interactions, in which the metabolic excretions of one species are the primary resource for another, constitute the central drivers of microbial community assembly. Trophic interactions emerge from the deconstruction of complex forms of organic matter into a wealth of smaller metabolic intermediates, some of which are released to the environment and serve as a nutritional buffet for the community. The structure of the emergent trophic network and the rate at which primary resources are supplied control many features of microbial community assembly, including the relative contributions of competition and cooperation and the emergence of alternative community states. Viewing microbial community assembly through the lens of trophic interactions also has important implications for the spatial dynamics of communities as well as the functional redundancy of taxonomic groups. Given the ubiquity of trophic interactions across environments, they impart a common logic that can enable the development of a more quantitative and predictive microbial community ecology.

Profiling temporal dynamics of acetogenic communities in anaerobic digesters using next-generation sequencing and T-RFLP

Acetogens play a key role in anaerobic degradation of organic material and in maintaining biogas process efficiency. Profiling this community and its temporal changes can help evaluate process stability and function, especially under disturbance/stress conditions, and avoid complete process failure. The formyltetrahydrofolate synthetase (FTHFS) gene can be used as a marker for acetogenic community profiling in diverse environments. In this study, we developed a new high-throughput FTHFS gene sequencing method for acetogenic community profiling and compared it with conventional terminal restriction fragment length polymorphism of the FTHFS gene, 16S rRNA gene-based profiling of the whole bacterial community, and indirect analysis via 16S rRNA profiling of the FTHFS gene-harbouring community. Analyses and method comparisons were made using samples from two laboratory-scale biogas processes, one operated under stable control and one exposed to controlled overloading disturbance. Comparative analysis revealed satisfactory detection of the bacterial community and its changes for all methods, but with some differences in resolution and taxonomic identification. FTHFS gene sequencing was found to be the most suitable and reliable method to study acetogenic communities. These results pave the way for community profiling in various biogas processes and in other environments where the dynamics of acetogenic bacteria have not been well studied.

Organic Matter Type Defines the Composition of Active Microbial Communities Originating From Anoxic Baltic Sea Sediments

Carbon cycling in anoxic marine sediments is dependent on uncultured microbial communities. Niches of heterotrophic microorganisms are defined by organic matter (OM) type and the different phases in OM degradation. We investigated how OM type defines microbial communities originating from organic-rich, anoxic sediments from the Baltic Sea. We compared changes in the sediment microbial community, after incubation with different stable isotope labeled OM types [i.e., particulate algal organic matter (PAOM), protein, and acetate], by using DNA stable isotope probing (DNA-SIP). Incorporation of 13C and/or 15N label was predominantly detected in members of the phyla Planctomycetes and Chloroflexi, which also formed the majority (>50%) of the original sediment community. While these phylum-level lineages incorporated label from all OM types, phylogenetic analyses revealed a niche separation at the order level. Members of the MSBL9 (Planctomycetes), the Anaerolineales (Chloroflexi), and the class Bathyarchaeota, were identified as initial degraders of carbohydrate-rich OM, while other uncultured orders, like the CCM11a and Phycisphaerales (Planctomycetes), Dehalococcoidia, and JG30-KF-CM66 (Chloroflexi), incorporated label also from protein and acetate. Our study highlights the importance of initial fermentation of complex carbon pools in shaping anoxic sediment microbial communities and reveals niche specialization at the order level for the most important initial degraders in anoxic sediments.

Overlooked Ecological Roles of Influent Wastewater Microflora in Improving Biological Phosphorus Removal in an Anoxic/Aerobic MBR Process

The ecological roles of influent microflora in activated sludge communities have not been well investigated. Herein, parallel lab-scale anoxic/aerobic (A/O) membrane bioreactors (MBRs), which were fed with raw (MBR-C) and sterilized (MBR-T) municipal wastewater, were operated. The MBRs showed comparable nitrogen removal but superior phosphorus removal in MBR-C than MBR-T over the long-term operation. The MBR-C sludge community had higher diversity and deterministic assembly than the MBR-T sludge community as revealed by 16S rRNA gene sequencing and null model analysis. Moreover, the MBR-C sludge community had higher abundance of polyphosphate accumulating organisms (PAOs) and hydrolytic/fermentative bacteria (HFB) but lower abundance of glycogen-accumulating organisms (GAOs), in comparison with MBR-T sludge. Intriguingly, the results of both the net growth rate and Sloan's neutral model demonstrated that HFB in the sludge community were generally slow-growing or nongrowing and their consistent presence in activated sludge was primarily attributed to the HFB immigration from influent microflora. Positive correlations between PAOs and HFB and potential competitions between HFB and GAOs were observed, as revealed by the putative species-species associations in the ecological networks. Taken together, this work deciphers the positive ecological roles of influent microflora, particularly HFB, in system functioning and highlights the necessity of incorporating influent microbiota for the design and modeling of A/O MBR plants.

Microbial Communities in Full-Scale Wastewater Treatment Systems Exhibit Deterministic Assembly Processes and Functional Dependency over Time

Microbial communities constitute the core component of biological wastewater treatment processes. We conducted a meta-analysis based on the 16S rRNA gene of temporal samples obtained from diverse full-scale activated sludge and anaerobic digestion systems treating municipal and industrial wastewater (collected in this study and published previously) to investigate their community assembly mechanism and functional traits over time, which are not currently well understood. The influent composition was found to be the main driver of the microbial community's composition, and relatively large proportions of specialist (26.1% and 18.6%) and transient taxa (67.2% and 68.1%) were estimated in both systems. Deterministic processes, especially homogeneous selection events (accounting for >53.8% of assembly events), were consistently identified as the dominant microbial community assembly mechanisms in both systems over time. Significant and strong correlations (Pearson's r = 0.51-0.92) were detected between the dynamics of the temporal community and the functional compositions in both systems, which suggests functional dependency. In contrast, the occurrence of sludge bulking and foaming in the activated sludge system led to an increase in stochastic assembly processes (i.e., limited dispersal and undominated events), a shift toward functional redundancy and less community diversity, a decreased community niche breadth index, and a more compact co-association network. This study illustrates that the mechanism of microbial community assembly and functional traits over time can be used to diagnose system performance and provide information on potential system malfunction.

Profiling microbial community structures and functions in bioremediation strategies for treating 1,4-dioxane-contaminated groundwater

Microbial community compositions and functional profiles were analyzed in microcosms established using aquifer materials from a former automobile factory site, where 1,4-dioxane was identified as the primary contaminant of concern. Propane or oxygen biostimulation resulted in limited 1,4-dioxane degradation, which was markedly enhanced with the addition of nutrients, resulting in abundant Mycobacterium and Methyloversatilis taxa and high expressions of propane monooxygenase gene, prmA. In bioaugmented treatments, Pseudonocardia dioxanivorans CB1190 or Rhodococcus ruber ENV425 strains dominated immediately after augmentation and degraded 1,4-dioxane rapidly which was consistent with increased representation of xenobiotic and lipid metabolism-related functions. Although the bioaugmented microbes decreased due to insufficient growth substrates and microbial competition, they did continue to degrade 1,4-dioxane, presumably by indigenous propanotrophic and heterotrophic bacteria, inducing similar community structures across bioaugmentation conditions. In various treatments, functional redundancy acted as buffer capacity to ensure a stable microbiome, drove the restoration of the structure and microbial functions to original levels, and induced the decoupling between basic metabolic functions and taxonomy. The results of this study provided valuable information for design and decision-making for ex-situ bioreactors and in-situ bioremediation applications. A metagenomics-based understanding of the treatment process will enable efficient and accurate adjustments when encountering unexpected issues in bioremediation.

Sludge digestibility and functionally active microorganisms in methanogenic sludge digesters revealed by E. coli-fed digestion and microbial source tracking

Methanogenic sludge digestion plays a pivotal role in attenuating and hygienizing the massively-produced waste activated sludge (WAS), which is predominantly composed of microbial cells and extracellular polymeric substances (EPS). The efficient sludge digestion requires a variety of functionally active microorganisms working together closely to convert sludge organic matter into biogas. Nonetheless, the digestion efficiency (or digestibility quantified as carbon removal efficiency) of major sludge constituents (i.e., microbial cells and EPS) and associated functionally active microorganisms in sludge digesters remain elusive. In this study, we identified the digestibility of sludge microbial cells and the associated functionally active microorganisms by using Escherichia coli (E. coli)-fed digestion and microbial source tracking. The average carbon removals in four digesters fed with fresh WAS (WAS-AD), thermal pretreated WAS (Thermal-WAS-AD), E. coli cells (E.coli-AD) and thermal pretreated E. coli cells (Thermal-E.coli-AD) were 30.6 ± 3.4%, 45.8 ± 2.9%, 69.0 ± 3.4% and 68.9 ± 4.6%, respectively. Compared to WAS-AD and Thermal-WAS-AD, the significantly higher carbon removals in E. coli-AD and Thermal-E. coli-AD suggested the remarkably higher digestibility of microbial cells than EPS, and releasing organic matter from EPS might be a rate-limiting step in sludge digestion. Functionally active microorganisms for microbial cell digestion predominantly included fermenters (e.g., Petrimonas and Lentimicrobium), syntrophic acetogens (e.g., Synergistaceae) and methanogens (e.g., Methanosaeta and Methanosarcina). Microbial source tracking estimation showed that the microbial cell-digesting populations accounted for 35.6 ± 9.1% and 70.3 ± 10.1% of total microbial communities in the WAS-AD and Thermal-WAS-AD, respectively. Accordingly, the functionally active microorganisms for digestion of both microbial cells and EPS accounted for 64.5 ± 12.1% and 97.3 ± 2.0% of total digestion sludge microbiome in WAS-AD and Thermal-WAS-AD, respectively. By contrast, feeding WAS-derived microorganisms accounted for 23.2 ± 4.4% and 2.3 ± 1.2% of total microbial communities in the WAS-AD and Thermal-WAS-AD, respectively.

Genomic and functional analyses of fungal and bacterial consortia that enable lignocellulose breakdown in goat gut microbiomes

The herbivore digestive tract is home to a complex community of anaerobic microbes that work together to break down lignocellulose. These microbiota are an untapped resource of strains, pathways and enzymes that could be applied to convert plant waste into sugar substrates for green biotechnology. We carried out more than 400 parallel enrichment experiments from goat faeces to determine how substrate and antibiotic selection influence membership, activity, stability and chemical productivity of herbivore gut communities. We assembled 719 high-quality metagenome-assembled genomes (MAGs) that are unique at the species level. More than 90% of these MAGs are from previously unidentified herbivore gut microorganisms. Microbial consortia dominated by anaerobic fungi outperformed bacterially dominated consortia in terms of both methane production and extent of cellulose degradation, which indicates that fungi have an important role in methane release. Metabolic pathway reconstructions from MAGs of 737 bacteria, archaea and fungi suggest that cross-domain partnerships between fungi and methanogens enabled production of acetate, formate and methane, whereas bacterially dominated consortia mainly produced short-chain fatty acids, including propionate and butyrate. Analyses of carbohydrate-active enzyme domains present in each anaerobic consortium suggest that anaerobic bacteria and fungi employ mostly complementary hydrolytic strategies. The division of labour among herbivore anaerobes to degrade plant biomass could be harnessed for industrial bioprocessing.

Genome-Centric Metagenomic Insights into the Impact of Alkaline/Acid and Thermal Sludge Pretreatment on the Microbiome in Digestion Sludge

Pretreatment of waste-activated sludge (WAS) is an effective way to destabilize sludge floc structure and release organic matter for improving sludge digestion efficiency. Nonetheless, information on the impact of WAS pretreatment on digestion sludge microbiomes, as well as mechanistic insights into how sludge pretreatment improves digestion performance, remains elusive. In this study, a genome-centric metagenomic approach was employed to investigate the digestion sludge microbiome in four sludge digesters with different types of feeding sludge: WAS pretreated with 0.25 mol/liter alkaline/acid (APAD), WAS pretreated with 0.8 mol/liter alkaline/acid (HS-APAD), thermally pretreated WAS (thermal-AD), and fresh WAS (control-AD). We retrieved 254 metagenome-assembled genomes (MAGs) to identify the key functional populations involved in the methanogenic digestion process. These MAGs span 28 phyla, including 69 yet-to-be-cultivated lineages, and 30 novel lineages were characterized with metabolic potential associated with hydrolysis and fermentation. Interestingly, functional populations involving carbohydrate digestion were enriched in APAD and HS-APAD, while lineages related to protein and lipid fermentation were enriched in thermal-AD, corroborating the idea that different substrates are released from alkaline/acid and thermal pretreatments. Among the major functional populations (i.e., fermenters, syntrophic acetogens, and methanogens), significant correlations between genome sizes and abundance of the fermenters were observed, particularly in APAD and HS-APAD, which had improved digestion performance.IMPORTANCE Wastewater treatment generates large amounts of waste-activated sludge (WAS), which consists mainly of recalcitrant microbial cells and particulate organic matter. Though WAS pretreatment is an effective way to release sludge organic matter for subsequent digestion, detailed information on the impact of the sludge pretreatment on the digestion sludge microbiome remains scarce. Our study provides unprecedented genome-centric metagenomic insights into how WAS pretreatments change the digestion sludge microbiomes, as well as their metabolic networks. Moreover, digestion sludge microbiomes could be a unique source for exploring microbial dark matter. These results may inform future optimization of methanogenic sludge digestion and resource recovery.

Proposal to reclassify the proteobacterial classes Deltaproteobacteria and Oligoflexia, and the phylum Thermodesulfobacteria into four phyla reflecting major functional capabilities

The class Deltaproteobacteria comprises an ecologically and metabolically diverse group of bacteria best known for dissimilatory sulphate reduction and predatory behaviour. Although this lineage is the fourth described class of the phylum Proteobacteria , it rarely affiliates with other proteobacterial classes and is frequently not recovered as a monophyletic unit in phylogenetic analyses. Indeed, one branch of the class Deltaproteobacteria encompassing Bdellovibrio-like predators was recently reclassified into a separate proteobacterial class, the Oligoflexia . Here we systematically explore the phylogeny of taxa currently assigned to these classes using 120 conserved single-copy marker genes as well as rRNA genes. The overwhelming majority of markers reject the inclusion of the classes Deltaproteobacteria and Oligoflexia in the phylum Proteobacteria . Instead, the great majority of currently recognized members of the class Deltaproteobacteria are better classified into four novel phylum-level lineages. We propose the names Desulfobacterota phyl. nov. and Myxococcota phyl. nov. for two of these phyla, based on the oldest validly published names in each lineage, and retain the placeholder name SAR324 for the third phylum pending formal description of type material. Members of the class Oligoflexia represent a separate phylum for which we propose the name Bdellovibrionota phyl. nov. based on priority in the literature and general recognition of the genus Bdellovibrio. Desulfobacterota phyl. nov. includes the taxa previously classified in the phylum Thermodesulfobacteria , and these reclassifications imply that the ability of sulphate reduction was vertically inherited in the Thermodesulfobacteria rather than laterally acquired as previously inferred. Our analysis also indicates the independent acquisition of predatory behaviour in the phyla Myxococcota and Bdellovibrionota, which is consistent with their distinct modes of action. This work represents a stable reclassification of one of the most taxonomically challenging areas of the bacterial tree and provides a robust framework for future ecological and systematic studies.

Metagenomic and Metatranscriptomic Analyses Revealed Uncultured Bacteroidales Populations as the Dominant Proteolytic Amino Acid Degraders in Anaerobic Digesters

Current understanding of amino acid (AA) degraders in anaerobic digesters is mainly based on cultured species, whereas microorganisms that play important roles in a complex microbial community remain poorly characterized. This study investigated short-term enrichments degrading single AAs using metagenomics and metatranscriptomics. Metagenomic analysis revealed that populations related to cultured AA degraders had an abundance <2.5% of the sequences. In contrast, metagenomic-assembled bins related to uncultured Bacteroidales collectively accounted for >35% of the sequences. Phylogenetic analyses suggested that these Bacteroidales populations represented a yet-to-be characterized family lineage, i.e., Bacteroidetes vadinHA17. The bins possessed the genetic capacity related to protein degradation, including surface adhesion (3–7 genes), secreted peptidase (52–77 genes), and polypeptide-specific transporters (2–5 genes). Furthermore, metatranscriptomics revealed that these Bacteroidales populations expressed the complete metabolic pathways for degrading 16 to 17 types of AAs in enrichments fed with respective substrates. These characteristics were distinct from cultured AA degraders including Acidaminobacter and Peptoclostridium, suggesting the uncultured Bacteroidales were the major protein-hydrolyzing and AA-degrading populations. These uncultured Bacteroidales were further found to be dominant and active in full-scale anaerobic digesters, indicating their important ecological roles in the native habitats. “Candidatus Aminobacteroidaceae” was proposed to represent the previously uncharted family Bacteroidetes vadinHA17.

Functional redundancy imparts process stability to acidic Fe(II)-oxidizing microbial reactors

In previous work, lab-scale reactors designed to study microbial Fe(II) oxidation rates at low pH were found to have stable rates under a wide range of pH and Fe(II) concentrations. Since the stirred reactor environment eliminates many of the temporal and spatial variations that promote high diversity among microbial populations in nature, we were surprised that the reactors supported multiple taxa presumed to be autotrophic Fe(II) oxidizers based on their phylogeny. Metagenomic analyses of the reactor communities revealed differences in the metabolic potential of these taxa with respect to Fe(II) oxidation and carbon fixation pathways, acquisition of potentially growth-limiting substrates and the ability to form biofilms. Our findings support the hypothesis that the long-term co-existence of multiple autotrophic Fe(II)-oxidizing populations in the reactors are due to distinct metabolic potential that supports differential growth in response to limiting resources such as nitrogen, phosphorus and oxygen. Our data also highlight the role of biofilms in creating spatially distinct geochemical niches that enable the co-existence of multiple taxa that occupy the same apparent metabolic niche when the system is viewed in bulk. The distribution of key metabolic functions across different co-existing taxa supported functional redundancy and imparted process stability to these reactors.

Effects of forced taxonomic transitions on metabolic composition and function in microbial microcosms

Surveys of microbial systems indicate that in many situations taxonomy and function may constitute largely independent ('decoupled') axes of variation. However, this decoupling is rarely explicitly tested experimentally, partly because it is hard to directly induce taxonomic variation without affecting functional composition. Here we experimentally evaluate this paradigm using microcosms resembling lake sediments and subjected to two different levels of salinity (0 and 19) and otherwise similar environmental conditions. We used DNA sequencing for taxonomic and functional profiling of bacteria and archaea and physicochemical measurements to monitor metabolic function, over 13 months. We found that the taxonomic composition of the saline systems gradually but strongly diverged from the fresh systems. In contrast, the metabolic composition (in terms of proportions of various genes) remained nearly identical across treatments and over time. Oxygen consumption rates and methane concentrations were substantially lower in the saline treatment, however, their similarity either increased (for oxygen) or did not change significantly (for methane) between the first and last sampling time, indicating that the lower metabolic activity in the saline treatments was directly and immediately caused by salinity rather than the gradual taxonomic divergence. Our experiment demonstrates that strong taxonomic shifts need not directly affect metabolic rates.

Transition of microbial communities and degradation pathways in anaerobic digestion at decreasing retention time

Tuning of operational variables is a common practice to control the anaerobic digestion process and, in advanced applications, to promote the accumulation of fermentation products. However, process variables are interrelated. In this study, the hydraulic retention time (HRT) was decoupled from the organic loading rate (OLR) in order to isolate the effect of HRT as a selective pressure on: process performance, metabolic rates (hydrolytic, acetogenic, and methanogenic) and the microbial community. Four mesophilic anaerobic digesters were subjected to a sequential decrease in HRT from 15 to 8, 4 and 2 days while keeping the OLR constant at chemical oxygen demand of 1 gCOD L _r^-1 d^-1. The results showed that HRT alone was insufficient to washout methanogens from the digesters, which in turn prevented the accumulation of volatile fatty acids (VFA). Methanosaeta was the dominant genus in the four digesters at all HRTs. Metabolic rates showed that process performance was controlled by hydrolysis, with a clear shift in acetogenic rates, from butyrate and propionate degradation to ethanol degradation at 4 and 2d HRT. The change in acetogenic pathways was attributed to a shift in the fermentation pathways co-current with changes in fermentative bacteria. At 2d HRT, biofilm was formed on the walls and paddles of the digesters, probably as a survival strategy. Most of the taxa in the biofilm were also present in the digester media. Overall, it is the combination of HRT with other operational parameters which promotes the washout of methanogens and the accumulation of VFA.

Genome-centric microbiome analysis reveals solid retention time (SRT)-shaped species interactions and niche differentiation in food waste and sludge co-digesters

Co-digestion of food waste with sewage sludge is widely applied for waste stabilization and energy recovery around the world. However, the effect of solid retention time (SRT) on the microbial population dynamics, metabolism and interspecies interaction have not been fully elucidated. Here, the influence of SRTs (5-25 days) on the performance of the co-digestion system was investigated and state-of-the-art genome-centric metagenomic analysis was employed to uncover the dynamics and metabolic network of the key players underlying the well-functioned and poorly-functioned co-digestion microbial communities. The results of the microbial analyses indicated that SRT largely shaped microbial community structure by enriching the syntrophic specialist Syntrophomonas and CO₂/H₂ ( formate)-using methanogen Methanocorpusculum in the well-functioned co-digester operated at SRT of 25 days, while selecting acid-tolerant populations Lactobacillus at SRT of 5 days. The metagenome assembled genomes (MAGs) of key players, such as Syntrophomonadaceae, Methanocorpusculum, and Mesotoga, were retrieved, additionally, the syntrophic acetate oxidation plus hydrogenotrophic methanogenesis (SAO-HM) were proposed as the dominant pathway for methane production. The metabolic interaction in the co-digestion microbial consortia was profiled by assigning MAGs into functional guilds. Functional redundancy was found in the bacterial groups in hydrolysis step, and the members in these groups reduced the direct competition by niche differentiation.

Disentangling Community Structure of Ecological System in Activated Sludge: Core Communities, Functionality, and Functional Redundancy

The microbial ecosystems of the sludge were characterized in terms of the core community structure, functional pathways, and functional redundancy through Illumina MiSeq sequencing and PICRUSt analysis on the activated sludge (AS) samples from an extended activated aeration process. Based on the identified OTU distribution, we identified 125 core community genera, including 3 abundant core genera and 21 intermittent abundant core genera. Putative genera Nitrosomonas, Nitrotoga, Zoogloea, Novosphingobium, Thermomonas, Amaricoccus, Tetrasphaera, Candidatus Microthrix, and Haliscomenobacter, which are associated with functions of nitrifying, denitrifying, phosphorus accumulating, and bulking and foaming, were found to present as the core community organisms in the AS sampled from the conventional extended aeration AS processes. The high-abundant nitrogen metabolic pathways were associated with nitrate reduction to ammonium (DNRA and ANRA), denitrification, and nitrogen fixation, while the ammonia oxidation-related genes (amo) were rarely annotated in the AS samples. Strict functional redundancy was not found with the AS ecosystem as it showed a high correlation between the community composition similarity and function similarity. In addition, the classified dominant core genera community was found to be sufficient to characterize the functionality of AS, which could invigorate applications of 16S rDNA MiSeq sequencing and PICRUSt for the prediction of functions of AS ecosystems.

Syntrophic acetate oxidation replaces acetoclastic methanogenesis during thermophilic digestion of biowaste

BACKGROUND

Anaerobic digestion (AD) is a globally important technology for effective waste and wastewater management. In AD, microorganisms interact in a complex food web for the production of biogas. Here, acetoclastic methanogens and syntrophic acetate-oxidizing bacteria (SAOB) compete for acetate, a major intermediate in the mineralization of organic matter. Although evidence is emerging that syntrophic acetate oxidation is an important pathway for methane production, knowledge about the SAOB is still very limited.

RESULTS

A metabolic reconstruction of metagenome-assembled genomes (MAGs) from a thermophilic solid state biowaste digester covered the basic functions of the biogas microbial community. Firmicutes was the most abundant phylum in the metagenome (53%) harboring species that take place in various functions ranging from the hydrolysis of polymers to syntrophic acetate oxidation. The Wood-Ljungdahl pathway for syntrophic acetate oxidation and corresponding genes for energy conservation were identified in a Dethiobacteraceae MAG that is phylogenetically related to known SAOB. 16S rRNA gene amplicon sequencing and enrichment cultivation consistently identified the uncultured Dethiobacteraceae together with Syntrophaceticus, Tepidanaerobacter, and unclassified Clostridia as members of a potential acetate-oxidizing core community in nine full-scare digesters, whereas acetoclastic methanogens were barely detected.

CONCLUSIONS

Results presented here provide new insights into a remarkable anaerobic digestion ecosystem where acetate catabolism is mainly realized by Bacteria. Metagenomics and enrichment cultivation revealed a core community of diverse and novel uncultured acetate-oxidizing bacteria and point to a particular niche for them in dry fermentation of biowaste. Their genomic repertoire suggests metabolic plasticity besides the potential for syntrophic acetate oxidation. Video Abstract.

New insights from the biogas microbiome by comprehensive genome-resolved metagenomics of nearly 1600 species originating from multiple anaerobic digesters

BACKGROUND

Microorganisms in biogas reactors are essential for degradation of organic matter and methane production. However, a comprehensive genome-centric comparison, including relevant metadata for each sample, is still needed to identify the globally distributed biogas community members and serve as a reliable repository.

RESULTS

Here, 134 publicly available metagenomes derived from different biogas reactors were used to recover 1635 metagenome-assembled genomes (MAGs) representing different biogas bacterial and archaeal species. All genomes were estimated to be > 50% complete and nearly half ≥ 90% complete with ≤ 5% contamination. In most samples, specialized microbial communities were established, while only a few taxa were widespread among the different reactor systems. Metabolic reconstruction of the MAGs enabled the prediction of functional traits related to biomass degradation and methane production from waste biomass. An extensive evaluation of the replication index provided an estimation of the growth dynamics for microbes involved in different steps of the food chain.

CONCLUSIONS

The outcome of this study highlights a high flexibility of the biogas microbiome, allowing it to modify its composition and to adapt to the environmental conditions, including temperatures and a wide range of substrates. Our findings enhance our mechanistic understanding of the AD microbiome and substantially extend the existing repository of genomes. The established database represents a relevant resource for future studies related to this engineered ecosystem.