Toward an Intensive Longitudinal Understanding of Activated Sludge Bacterial Assembly and Dynamics

Longitudinal dynamics and cross-domain interactions of eukaryotic populations in wastewater treatment plants

Activated sludge is a large reservoir of novel microorganisms and microbial genetic diversity. While much attention has been given to the profile and functions of prokaryotes, the eukaryotic diversity remains largely unexplored. In this study, we analysed longitudinal activated sludge samples spanning 13 years from the largest secondary wastewater treatment plants in Hong Kong, unveiling a wealth of eukaryotic taxa and 681 856 non-redundant protein-coding genes, the majority (416 044) of which appeared novel. Ciliophora was the most dominant phylum with a significant increase after a transient intervention (bleaching event). Our metagenomic analysis reveals close linkage and covariation of eukaryotes, prokaryotes, and prokaryotic viruses (phages), indicating common responses to environmental changes such as transient intervention and intermittent fluctuations. Furthermore, high-resolution cross-domain relationships were interpreted by S-map, demonstrating a predatory role of Arthropoda, Ascomycota, Mucoromycota, and Rotifera. This high-resolution profile of microbial dynamics expands our knowledge on yet-to-be-cultured populations and their cross-domain interactions and highlights the ecological importance of eukaryotes in the activated sludge ecosystem.

High-throughput single-cell sequencing of activated sludge microbiome

Wastewater treatment plants (WWTPs) represent one of biotechnology's largest and most critical applications, playing a pivotal role in environmental protection and public health. In WWTPs, activated sludge (AS) plays a major role in removing contaminants and pathogens from wastewater. While metagenomics has advanced our understanding of microbial communities, it still faces challenges in revealing the genomic heterogeneity of cells, uncovering the microbial dark matter, and establishing precise links between genetic elements and their host cells as a bulk method. These issues could be largely resolved by single-cell sequencing, which can offer unprecedented resolution to show the unique genetic information. Here we show the high-throughput single-cell sequencing to the AS microbiome. The single-amplified genomes (SAGs) of 15,110 individual cells were clustered into 2,454 SAG bins. We find that 27.5% of the genomes in the AS microbial community represent potential novel species, highlighting the presence of microbial dark matter. Furthermore, we identified 1,137 antibiotic resistance genes (ARGs), 10,450 plasmid fragments, and 1,343 phage contigs, with shared plasmid and phage groups broadly distributed among hosts, indicating a high frequency of horizontal gene transfer (HGT) within the AS microbiome. Complementary analysis using 1,529 metagenome-assembled genomes from the AS samples allowed for the taxonomic classification of 98 SAG bins, which were previously unclassified. Our study establishes the feasibility of single-cell sequencing in characterizing the AS microbiome, providing novel insights into its ecological dynamics, and deepening our understanding of HGT processes, particularly those involving ARGs. Additionally, this valuable tool could monitor the distribution, spread, and pathogenic hosts of ARGs both within AS environments and between AS and other environments, which will ultimately contribute to developing a health risk evaluation system for diverse environments within a One Health framework.

Deep insights into the assembly mechanisms, co-occurrence patterns, and functional roles of microbial community in wastewater treatment plants

The understanding of activated sludge microbial status and roles is imperative for improving and enhancing the performance of wastewater treatment plants (WWTPs). In this study, we conducted a deep analysis of activated sludge microbial communities across five compartments (inflow, effluent, and aerobic, anoxic, anaerobic tanks) over temporal scales, employing high-throughput sequencing of 16S rRNA amplicons and metagenome data. Clearly discernible seasonal patterns, exhibiting cyclic variations, were observed in microbial diversity, assembly, co-occurrence network, and metabolic functions. Notably, summer samples exhibited higher α-diversity and were distinctly separated from winter samples. Our analysis revealed that microbial community assembly is influenced by both stochastic processes (66%) and deterministic processes (34%), with winter samples demonstrating more random assembly compared to summer. Co-occurrence patterns were predominantly mutualistic, with over 96% positive correlations, and summer networks were more organized than those in winter. These variations were significantly correlated with temperature, total phosphorus and sludge volume index. However, no significant differences were found among microbial community across five compartments in terms of β diversity. A core community of keystone taxa was identified, playing key roles in eight nitrogen and eleven phosphorus cycling pathways. Understanding the assembly mechanisms, co-occurrence patterns, and functional roles of microbial communities is essential for the design and optimization of biotechnological treatment processes in WWTPs.

Assembly mechanism of microbial community under different seasons in Shantou sea area

Coastal areas are often affected by a variety of climates, and microbial composition patterns are conducive to adaptation to these environments. In this study, the composition and pattern of microbial communities in the Shantou sea from four seasons were analyzed. The diversity of microbial community was significant differences under different seasons (p < 0.01). Meanwhile, dissolved oxygen levels, temperature were key factors to shift microbial communities. The assembly mechanism of microbial communities was constructed by the iCAMP (Infer community assembly mechanism by the phylogenetic bin-based null). Interestingly, the analyses revealed that drift was the predominant driver of this process (44.5 %), suggesting that microbial community assembly in this setting was dominated by stochastic processes. For example, Vibrio was found to be particularly susceptible to stochastic processes, indicating that the pattern of bacterial community was governed by stochastic processes. Thus, these results offering novel insight into the regulation of microbial ecology in marine environments.

Filamentous bacteria-induced sludge bulking can alter antibiotic resistance gene profiles and increase potential risks in wastewater treatment systems

Sludge bulking caused by filamentous bacteria is a prevalent issue in wastewater treatment systems. While previous studies have primarily concentrated on controlling sludge bulking, the biological risks associated with it have been overlooked. This study demonstrates that excessive growth of filamentous bacteria during sludge bulking can significantly increase the abundance of antibiotic resistance genes (ARGs) in activated sludge. Through metagenomic analysis, we identified specific ARGs carried by filamentous bacteria, such as Sphaerotilus and Thiothrix, which are responsible for bulking. Additionally, by examining over 1,000 filamentous bacterial genomes, we discovered a diverse array of ARGs across different filamentous bacteria derived from wastewater treatment systems. Our findings indicate that 74.84% of the filamentous bacteria harbor at least one ARG, with the occurrence frequency of ARGs in these bacteria being approximately 1.5 times higher than that in the overall bacterial population in activated sludge. Furthermore, genomic and metagenomic analyses have shown that the ARGs in filamentous bacteria are closely linked to mobile genetic elements and are frequently found in potentially pathogenic bacteria, highlighting potential risks posed by these filamentous bacteria. These insights enhance our understanding of ARGs in activated sludge and underscore the importance of risk management in wastewater treatment systems.

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.

Sulfur autotrophic denitrification as an efficient nitrogen removals method for wastewater treatment towards lower organic requirement: A review

The sulfur autotrophic denitrification (SADN) process is an organic-free denitrification process that utilizes reduced inorganic sulfur compounds (RISCs) as the electron donor for nitrate reduction. It has been proven to be a cost-effective and environment-friendly approach to achieving carbon neutrality in wastewater treatment plants. However, there is no consensus on whether SADN can become a dominant denitrification process to treat domestic wastewater or industrial wastewater if organic carbon is desired to be saved. Through a comprehensive summary of the SADN process and extensive discussion of state-of-the-art SADN-based technologies, this review provides a systematic overview of the potential of the SADN process as a sustainable alternative for the heterotrophic denitrification (HD) process (organic carbons as electron donor). First, we introduce the mechanism of the SADN process that is different from the HD process, including its transformation pathways based on different RISCs as well as functional bacteria and key enzymes. The SADN process has unique theoretical advantages (e.g., economy and carbon-free, less greenhouse gas emissions, and a great potential for coupling with novel autotrophic processes), even if there are still some potential issues (e.g., S intermediates undesired production, and relatively slow growth rate of sulfur-oxidizing bacteria [SOB]) for wastewater treatment. Then we present the current representative SADN-based technologies, and propose the outlooks for future research in regards to SADN process, including implement of coupling of SADN with other nitrogen removal processes (e.g., HD, and sulfate-dependent anaerobic ammonium oxidation), and formation of SOB-enriched biofilm. This review will provide guidance for the future applications of the SADN process to ensure a robust-performance and chemical-saving denitrification for wastewater treatment.

The microbial dark matter and "wanted list" in worldwide wastewater treatment plants

BACKGROUND

Wastewater treatment plants (WWTPs) are one of the largest biotechnology applications in the world and are of critical importance to modern urban societies. An accurate evaluation of the microbial dark matter (MDM, microorganisms whose genomes remain uncharacterized) proportions in WWTPs is of great value, while there is no such research yet. This study conducted a global meta-analysis of MDM in WWTPs with 317,542 prokaryotic genomes from the Genome Taxonomy Database and proposed a "wanted list" for priority targets in further investigations of activated sludge.

RESULTS

Compared with the Earth Microbiome Project data, WWTPs had relatively lower genome-sequenced proportions of prokaryotes than other ecosystems, such as the animal related environments. Analysis showed that the median proportions of the genome-sequenced cells and taxa (100% identity and 100% coverage in 16S rRNA gene region) in WWTPs reached 56.3% and 34.5% for activated sludge, 48.6% and 28.5% for aerobic biofilm, and 48.3% and 28.5% for anaerobic digestion sludge, respectively. This result meant MDM had high proportions in WWTPs. Besides, all of the samples were occupied by a few predominant taxa, and the majority of the sequenced genomes were from pure cultures. The global-scale "wanted list" for activated sludge contained four phyla that have few representatives and 71 operational taxonomic units with the majority of them having no genome or isolate yet. Finally, several genome mining methods were verified to successfully recover genomes from activated sludge such as hybrid assembly of the second- and third-generation sequencing.

CONCLUSIONS

This work elucidated the proportion of MDM in WWTPs, defined the "wanted list" of activated sludge for future investigations, and certified potential genome recovery methods. The proposed methodology of this study can be applied to other ecosystems and improve understanding of ecosystem structure across diverse habitats. Video Abstract.

Deterministic mechanisms drive bacterial communities assembly in industrial wastewater treatment system

Microbial communities are responsible for biological treatment of many industrial wastewater, but our knowledge of their diversity, assembly patterns, and function is still poor. Here, we analyzed the bacterial communities of wastewater and activated sludge samples taken from 11 full-scale industrial wastewater treatment plants (IWWTPs) characterized by the same process design but different wastewater types and WWTP compartments. We found significantly different diversity and compositions of bacterial assemblages among distinct wastewater types and IWWTPs compartments. IWWTPs bacterial communities exhibited a clear species abundance distribution. The dispersal-driven process was weak in shaping IWWTP communities. Meanwhile, environmental and operating conditions were important factors in regulating the structure of the activated sludge community and pollutants removal, indicating that bacterial community was largely driven by deterministic mechanisms. The core microbial community in IWWTPs was different from that in municipal wastewater treatment plants (MWWTPs), and many taxa (e.g. the genus Citreitalea) rarely were detected before, indicating IWWTPs harbored unique core bacterial communities. Furthermore, we found that bacterial community compositions were strongly linked to activated sludge function. These findings are important to both microbial ecologists and environmental engineers, who may optimize the operation strategies jointly for maintaining biodiversity, which in turn may promote a more stable performance of the IWWTP. Overall, our study enhances the mechanistic understanding of the IWWTP microbial community diversity, assembly patterns, and function, and provides important implications for microbial ecology and wastewater treatment processes.

State-Space-Based Framework for Predicting Microbial Interaction Variability in Wastewater Treatment Plants

Substantial attempts have been made to control microbial communities for environmental integrity, biosystem performance, and human health. However, it is difficult to manipulate microbial communities in practice due to the varying and nonlinear nature of interspecific interaction networks. Here, we develop a manifold-based framework to investigate the patterns of microbial interaction variability in wastewater treatment plants using manifold geometric properties and design a simple control strategy to manipulate the microbes in nonlinear communities. We validate our framework using the readily available and nonsequential microbiome profiles of wastewater treatment plants. Our results show that some microbes in the activated sludge and anammox communities display deterministic rival or cooperative relationships and constitute a stable subnetwork within the whole nonlinear community network. We further use a simulation to demonstrate that these microbes can be used to drive a microbe in a target direction regardless of the community dynamics. Overall, our framework can provide a time-efficient solution to select effective control inputs for reliable manipulation in varying microbial networks, opening up new possibilities across a range of biological fields, including wastewater treatment plants.

Causality and correlation analysis for deciphering the microbial interactions in activated sludge

Time series data has been considered to be a massive information provider for comprehending more about microbial dynamics and interaction, leading to a causality inference in a complex microbial community. Granger causality and correlation analysis have been investigated and applied for the construction of a microbial causal correlation network (MCCN) and efficient prediction of the ecological interaction within activated sludge, which thereby exhibited ecological interactions at the OTU-level. Application of MCCN to a time series of activated sludge data revealed that the hub species OTU56, classified as the one belonging to the genus Nitrospira, was responsible for nitrification in activated sludge and interaction with Proteobacteria and Bacteroidetes in the form of amensal and commensal relationships, respectively. The phylogenetic tree suggested a mutualistic relationship between Nitrospira and denitrifiers. Zoogloea displayed the highest ncf value within the classified OTUs of the MCCN, indicating that it could be a foundation for activated sludge through the formation of characteristic cell aggregate matrices where other organisms embed during floc formation. Inclusively, the research outcomes of this study have provided a deep insight into the ecological interactions within the communities of activated sludge.

Nationwide biogeography and health implications of bacterial communities in household drinking water

Drinking water at the point of use harbors microorganisms that may pose potential risks to human health. However, the microbial diversity and health impacts of household drinking water are poorly understood, since culture-based methods only target on specific microorganisms and low biomass of drinking water hinders a high-throughput profiling. Here, we used an optimized workflow to efficiently collect microorganisms from low-biomass drinking water and performed deep sequencing of 16S rRNA genes to profile the bacterial diversity and biogeography of 110 household drinking water samples covering 38 cities of 29 provinces/regions in China, and further explored environmental drivers and potential health implications. Our analyses revealed a diverse drinking water community comprising a total of 22,771 operational taxonomic units (OTUs). The spatial turnover of drinking water communities is scale-dependent and appears to be driven largely by rainfall and water source river. The identified potential pathogenic species may have the possibility of causing health risks. Our novel insights enhance the current understanding of the diversity and biogeography of drinking water bacterial communities within a theoretical ecological framework and have further important implications for safe drinking water management and public health protection.

Effects of influent immigration and environmental factors on bacterial assembly of activated sludge microbial communities

The functional mechanism of microbial assembly of activated sludge (AS) in urban wastewater treatment plants (UWTPs) remains unclear. A comprehensive quantitative evaluation of the contribution of influent immigration and environmental factors to AS community composition requires investigation. In this study, the microbial characteristics of six full-scale UWTPs with different influent compositions and environmental factors (altitude, temperature, dissolved oxygen (DO), pH, chemical oxygen demand (COD), total nitrogen (TN), ammonia nitrogen (NH₄⁺-N), and total phosphorus (TP)) were analyzed to determine the main forces affecting the bacterial assembly of AS microbial communities. Abundant and core taxa were screened out based on the abundance and frequency of operational taxonomic units (OTUs) occurrence in all samples. Abundant OTUs (18.7% occurrence) accounted for 87.7% of the total 16S rRNA sequences, while rare OTUs (71.7% occurrence) accounted for only 7.8% of the total 16S rRNA sequences. A total of 135 OTUs were identified as core taxa, accounting for 14.6-26.2% of the total reads, of which 83 OTUs belonged to abundant taxa. The richness and uniformity of the influent community were significantly lower than those of the AS system. The community composition in influent varied from that in AS. Moreover, about 89.7% (86.5% of 16S rRNA sequences) OTUs in AS samples showed positive growth rates, indicating that immigration of influent communities had a limited effect on the microbial composition of AS. Redundancy analysis (RDA) combined with co-occurrence network showed that the bacterial assembly of microbial communities was significantly correlated with altitude, pH, and TN (P < 0.05), and these three parameters could explain 23.3%, 21.1%, and 17.7% of the bacterial assembly of AS microbial communities in UWTPs, respectively.

Predation increases multiple components of microbial diversity in activated sludge communities

Protozoan predators form an essential component of activated sludge communities that is tightly linked to wastewater treatment efficiency. Nonetheless, very little is known how protozoan predation is channelled via bacterial communities to affect ecosystem functioning. Therefore, we experimentally manipulated protozoan predation pressure in activated-sludge communities to determine its impacts on microbial diversity, composition and putative functionality. Different components of bacterial diversity such as taxa richness, evenness, genetic diversity and beta diversity all responded strongly and positively to high protozoan predation pressure. These responses were non-linear and levelled off at higher levels of predation pressure, supporting predictions of hump-shaped relationships between predation pressure and prey diversity. In contrast to predation intensity, the impact of predator diversity had both positive (taxa richness) and negative (evenness and phylogenetic distinctiveness) effects on bacterial diversity. Furthermore, predation shaped the structure of bacterial communities. Reduction in top-down control negatively affected the majority of taxa that are generally associated with increased treatment efficiency, compromising particularly the potential for nitrogen removal. Consequently, our findings highlight responses of bacterial diversity and community composition as two distinct mechanisms linking protozoan predation with ecosystem functioning in activated sludge communities.

Microbial communities across activated sludge plants show recurring species-level seasonal patterns

Microbial communities in activated sludge (AS) are the core of sanitation in wastewater treatment plants (WWTPs). Microbial communities in AS have shown seasonal changes, however, long-term experiments (>2 years) are rarely conducted, limiting our understanding of the true seasonal dynamics in WWTPs. In this study, we resolved the microbial seasonal dynamics at the species level in four municipal full-scale WWTPs, sampled every 7-10 days, during 3-5 consecutive years. By applying a new time-series analysis approach, we revealed that the seasonal pattern was species-specific, where species belonging to the same functional guild or genus may show different seasonal dynamics. Species could be grouped into cohorts according to their seasonal patterns, where seasonal cohorts showed repeatable annual dynamics across years and plants. Species were also grouped according to their net growth rate in the AS (i.e., growing species and disappearing species). Growing species were more prevailing in spring and autumn cohorts, while disappearing species, which were only present due to the continuous immigration from influent wastewater, were mostly associated with winter and spring cohorts. Most known process-critical species, such as nitrifiers, polyphosphate accumulating organisms and filamentous organisms, showed distinct species-specific patterns. Overall, our study showed that overarching seasonal patterns affected microbial species in full-scale AS plants, with similar seasonal patterns across plants for many dominant species. These recurrent seasonal variations should be taken into account in the operation, understanding and management of the WWTPs.

Candidatus Kaistella beijingensis sp. nov., Isolated from a Municipal Wastewater Treatment Plant, Is Involved in Sludge Foaming

Biological foaming (or biofoaming) is a frequently occurring problem in wastewater treatment plants (WWTPs) and is attributed to the overwhelming growth of filamentous bulking and foaming bacteria (BFB). Biological foaming has been intensively investigated, with BFB like Microthrix and Skermania having been identified from WWTPs and implicated in foaming. Nevertheless, studies are still needed to improve our understanding of the microbial diversity of WWTP biofoams and how microbial activities contribute to foaming. In this study, sludge foaming at the Qinghe WWTP of China was monitored, and sludge foams were investigated using culture-dependent and culture-independent microbiological methods. The foam microbiomes exhibited high abundances of Skermania, Mycobacterium, Flavobacteriales, and Kaistella. A previously unknown bacterium, Candidatus Kaistella beijingensis, was cultivated from foams, its genome was sequenced, and it was phenotypically characterized. Ca. K. beijingensis exhibits hydrophobic cell surfaces, produces extracellular polymeric substances (EPS), and metabolizes lipids. Ca. K. beijingensis abundances were proportional to EPS levels in foams. Several proteins encoded by the Ca. K. beijingensis genome were identified from EPS that was extracted from sludge foams. Ca. K. beijingensis populations accounted for 4 to 6% of the total bacterial populations in sludge foam samples within the Qinghe WWTP, although their abundances were higher in spring than in other seasons. Cooccurrence analysis indicated that Ca. K. beijingensis was not a core node among the WWTP community network, but its abundances were negatively correlated with those of the well-studied BFB Skermania piniformis among cross-season Qinghe WWTP communities. IMPORTANCE Biological foaming, also known as scumming, is a sludge separation problem that has become the subject of major concern for long-term stable activated sludge operation in decades. Biological foaming was considered induced by foaming bacteria. However, the occurrence and deterioration of foaming in many WWTPs are still not completely understood. Cultivation and characterization of the enriched bacteria in foaming are critical to understand their genetic, physiological, phylogenetic, and ecological traits, as well as to improve the understanding of their relationships with foaming and performance of WWTPs.

Sequencing Batch Integrated Fixed-Film Activated Sludge Membrane Process for Treatment of Tapioca Processing Wastewater

Tapioca processing industries are very popular in the rural community to produce a variety of foods as the end products. Due to their small scales and scattered locations, they require robust modular systems to operate at low capacity with minimum supervision. This study explores the application of a novel sequencing batch-integrated fixed-film activated sludge membrane (SB-IFASM) process to treat tapioca processing wastewater for reuse purposes. The SB-IFASM employed a gravity-driven system and utilizes biofilm to enhance biodegradation without requiring membrane cleaning. The SB-IFASM utilizes the biofilm as a secondary biodegradation stage to enhance the permeate quality applicable for reuse. A lab-scale SB-IFASM was developed, preliminarily assessed, and used to treat synthetic tapioca processing industry wastewater. The results of short-term filtration tests showed the significant impact of hydrostatic pressure on membrane compaction and instant cake layer formation. Increasing the pressure from 2.2 to 10 kPa lowered the permeability of clean water and activated sludge from 720 to 425 and from 110 to 50 L/m²·h bar, respectively. The unsteady-state operation of the SB-IFASM showed the prominent role of the bio-cake in removing the organics reaching the permeate quality suitable for reuse. High COD removals of 63-98% demonstrated the prominence contribution of the biofilm in enhancing biological performance and ultimate COD removals of >93% make it very attractive for application in small-scale tapioca processing industries. However, the biological ecosystem was unstable, as shown by foaming that deteriorated permeability and was detrimental to the organic removal. Further developments are still required, particularly to address the biological stability and low permeability.

Microbial Community Composition in Municipal Wastewater Treatment Bioreactors Follows a Distance Decay Pattern Primarily Controlled by Environmental Heterogeneity

Understanding spatiotemporal patterns in microbial community composition is a central goal of microbial ecology. The objective of this study was to better understand the biogeography of activated sludge microbial communities, which are important for the protection of surface water quality. Monthly samples were collected from 20 facilities (25 bioreactors) within 442 km of each other for 1 year. Microbial community composition was characterized by sequencing of PCR-amplified 16S rRNA gene fragments. Statistically significant distance decay of community similarity was observed in these bioreactors independent of clustering method (operational taxonomic units [OTUs] at 97% similarity, genus-level phylotypes) and community dissimilarity metric (Sørensen, Bray-Curtis, and weighted Unifrac). Universal colonizers (i.e., detected in all samples) and ubiquitous genus-level phylotypes (i.e., detected in every facility at least once) also exhibited a significant distance decay relationship. Variation partitioning analysis of community composition showed that environmental characteristics (temperature, influent characteristics, etc.) explained more of the variance in community composition than geographic distance did, suggesting that environmental heterogeneity is more important than dispersal limitation as a mechanism for determining microbial community composition. Distance decay relationships also became stronger with increasing distance between facilities. Seasonal variation in community composition was also observed from selected bioreactors, but there was no clear seasonal pattern in the distance decay relationships.

IMPORTANCE Understanding the spatiotemporal patterns of biodiversity is a central goal of ecology. The distance decay of community similarity is one of the spatial scaling patterns observed in many forms of life, including plants, animals, and microbial communities. Municipal wastewater treatment relies on microorganisms to prevent the release of excessive quantities of nutrients and other pollutants, but relatively few studies have explored distance decay relationships in wastewater treatment bioreactors. Our results demonstrate a strong distance decay pattern in wastewater treatment bioreactors, regardless of the sequence clustering method or the community dissimilarity metric. Our results suggest that microbial communities in wastewater treatment bioreactors are not randomly assembled but rather exhibit a statistically significant spatial pattern.

Charting the complexity of the activated sludge microbiome through a hybrid sequencing strategy

BACKGROUND

Long-read sequencing has shown its tremendous potential to address genome assembly challenges, e.g., achieving the first telomere-to-telomere assembly of a gapless human chromosome. However, many issues remain unresolved when leveraging error-prone long reads to characterize high-complexity metagenomes, for instance, complete/high-quality genome reconstruction from highly complex systems.

RESULTS

Here, we developed an iterative haplotype-resolved hierarchical clustering-based hybrid assembly (HCBHA) approach that capitalizes on a hybrid (error-prone long reads and high-accuracy short reads) sequencing strategy to reconstruct (near-) complete genomes from highly complex metagenomes. Using the HCBHA approach, we first phase short and long reads from the highly complex metagenomic dataset into different candidate bacterial haplotypes, then perform hybrid assembly of each bacterial genome individually. We reconstructed 557 metagenome-assembled genomes (MAGs) with an average N50 of 574 Kb from a deeply sequenced, highly complex activated sludge (AS) metagenome. These high-contiguity MAGs contained 14 closed genomes and 111 high-quality (HQ) MAGs including full-length rRNA operons, which accounted for 61.1% of the microbial community. Leveraging the near-complete genomes, we also profiled the metabolic potential of the AS microbiome and identified 2153 biosynthetic gene clusters (BGCs) encoded within the recovered AS MAGs.

CONCLUSION

Our results established the feasibility of an iterative haplotype-resolved HCBHA approach to reconstruct (near-) complete genomes from highly complex ecosystems, providing new insights into "complete metagenomics". The retrieved high-contiguity MAGs illustrated that various biosynthetic gene clusters (BGCs) were harbored in the AS microbiome. The high diversity of BGCs highlights the potential to discover new natural products biosynthesized by the AS microbial community, aside from the traditional function (e.g., organic carbon and nitrogen removal) in wastewater treatment. Video Abstract.

Direct evidence for deterministic assembly of bacterial communities in full-scale municipal wastewater treatment facilities

In this study we investigated whether bacterial community composition in full-scale wastewater treatment bioreactors can be better explained by niche- or neutral- based theory (deterministic or stochastic) and whether bioreactor design (continuous-flow vs. fill-and-draw) affected community assembly. Four wastewater treatment facilities (one with quadruplicated continuous-flow bioreactors, two with one continuous-flow bioreactor each, one with triplicated fill-and-draw bioreactors) were investigated. Bioreactor community composition was characterized by sequencing of PCR-amplified 16S rRNA gene fragments. Replicate bioreactors at the same wastewater treatment facility had largely reproducible (i.e., deterministic) bacterial community composition, although bacterial community composition in continuous-flow bioreactors was significantly more reproducible (P < 0.001) than in fill-and-draw bioreactors (Bray-Curtis dissimilarity: μ = 0.48 ± 0.06 vs. 0.58 ± 0.08). Next, we compared our results to previously-used indirect methods for distinguishing between deterministic and stochastic community assembly mechanisms. Synchronicity was observed in the bacterial community composition among bioreactors within the same metropolitan region, consistent with deterministic community assembly. Similarly, a null model-based analysis also indicated that all wastewater bioreactor communities were controlled by deterministic factors and that continuous-flow bioreactors were significantly more deterministic (P < 0.001) than fill-and-draw bioreactors (nearest-taxon index: μ = 3.3 ± 0.6 vs. 2.7 ± 0.8). Our results indicate that bacterial community composition in wastewater treatment bioreactors is better explained by deterministic community assembly theory; simultaneously, our results validate previously-used but indirect methods to quantify whether microbial communities were assembled via deterministic or stochastic mechanisms. IMPORTANCE Understanding the mechanisms of bacterial community assembly is one of the grand challenges of microbial ecology. In environmental systems, this challenge is exacerbated because "replicate" experiments are typically impossible; that is, microbial ecologists cannot fabricate multiple field-scale experiments of identical, natural ecosystems. Our results directly demonstrate that deterministic mechanisms are more prominent than stochastic mechanisms in the assembly of wastewater treatment bioreactor communities. Our results also suggest that wastewater treatment bioreactor design is pertinent, such that the imposition of feast-famine conditions (i.e., fill-and-draw bioreactors) nudge bacterial community assembly more towards stochastic mechanisms compared to the imposition of stringent nutrient limitation (i.e., continuous-flow bioreactors). Our research also validates the previously-used indirect methods (synchronous community dynamics and an application of a null-model) for characterizing the relative importance of deterministic versus stochastic mechanisms of community assembly.

Seasonal dynamics of the microbial community in two full-scale wastewater treatment plants: Diversity, composition, phylogenetic group based assembly and co-occurrence pattern

The optimal operation and functional stability of a wastewater treatment plant (WWTP) strongly depend on the properties of its microbial community. However, a knowledge gap remains regarding the seasonal dynamics of microbial community properties, especially phylogenetic group based assembly and co-occurrence patterns. Accordingly, in this study, influent and activated sludge (AS) samples were weekly collected from 2 full-scale WWTPs for one year (89 influent and 103 AS samples in total) and examined by high-throughput Illumina-MiSeq sequencing. The results suggested that the microbial community diversity and composition in the influent fluctuated substantially with season, while those in the AS had a relatively more stable pattern throughout the year. The phylogenetic group based assembly mechanisms of AS community were identified by using "Infer Community Assembly Mechanisms by Phylogenetic-bin-based null model (iCAMP)". The results showed that drift accounted for the largest proportion (52.8%), while homogeneous selection (18.2%) was the most important deterministic process. Deterministic processes dominated in 47 microbial groups (bins), which were also found (~40%) in the AS core taxa dataset. Moreover, the results suggested that Nitrospira were more susceptible to stochastic processes in winter, which may provide a possible explanation for nitrification failure in winter. Network analysis results suggested that the network structure of the AS community could be more stable in summer and autumn. In addition, there were no identical keystone taxa found in different networks (constructed from different plants, sources, and seasons), which supported the context dependency theory. The results of this study deepened our understanding of the microbial ecology in AS systems and provided a foundation for further studies on the community regulation strategy of WWTPs.

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.

Seasonal Prevalence of Ammonia-Oxidizing Archaea in a Full-Scale Municipal Wastewater Treatment Plant Treating Saline Wastewater Revealed by a 6-Year Time-Series Analysis

Although several molecular-based studies have demonstrated the involvement of ammonia-oxidizing archaea (AOA) in ammonia oxidation in wastewater treatment plants (WWTPs), factors affecting the persistence and growth of AOA in these engineered systems have not been resolved. Here, we show a seasonal prevalence of AOA in a full-scale WWTP (Shatin, Hong Kong SAR) over a 6-year period of observation, even outnumbering ammonia-oxidizing bacteria in the seasonal peaks in 3 years, which may be due to the high bioavailable copper concentrations. Comparative analysis of three metagenome-assembled genomes of group I.1a AOA obtained from the activated sludge and 16S rRNA gene sequences recovered from marine sediments suggested that the seawater used for toilet flushing was the primary source of the WWTP AOA. A rare AOA population in the estuarine source water became transiently abundant in the WWTP with a metagenome-based relative abundance of up to 1.3% over three seasons of observation. Correlation-based network analysis revealed a robust co-occurrence relationship between these AOA and organisms potentially active in nitrite oxidation. Moreover, a strong correlation between the dominant AOA and an abundant proteobacterial organism suggested that capacity for extracellular polymeric substance production by the proteobacterium could provide a niche for AOA within bioaggregates. Together, the study highlights the importance of long-term observation in identifying biotic and abiotic factors governing population dynamics in open systems such as full-scale WWTPs.

Tropical and temperate wastewater treatment plants assemble different and diverse microbiomes

The diversity and assembly of activated sludge microbiomes play a key role in the performances of municipal wastewater treatment plants (WWTPs), which are the most widely applied biotechnological process systems. In this study, we investigated the microbiomes of municipal WWTPs in Bangkok, Wuhan, and Beijing that respectively represent tropical, subtropical, and temperate climate regions, and also explored how microbiomes assembled in these municipal WWTPs. Our results showed that the microbiomes from these municipal WWTPs were significantly different. The assembly of microbiomes in municipal WWTPs followed deterministic and stochastic processes governed by geographical location, temperature, and nutrients. We found that both taxonomic and phylogenetic α-diversities of tropical Bangkok municipal WWTPs were the highest and were rich in yet-to-be-identified microbial taxa. Nitrospirae and β-Proteobacteria were more abundant in tropical municipal WWTPs, but did not result in better removal efficiencies of ammonium and total nitrogen. Overall, these results suggest that tropical and temperate municipal WWTPs harbored diverse and unique microbial resources, and the municipal WWTP microbiomes were assembled with different processes. Implications of these findings for designing and running tropical municipal WWTPs were discussed. KEY POINTS: • Six WWTPs of tropical Thailand and subtropical and temperate China were investigated. • Tropical Bangkok WWTPs had more diverse and yet-to-be-identified microbial taxa. • Microbiome assembly processes were associated with geographical location.

Bacterial community evolution along full-scale municipal wastewater treatment processes

Sewage pollution is a major threat to public health because sewage is always accompanied by pathogens. Generally, wastewater treatment plants (WWTP) receive and treat sewage to control pathogenic risks and improve environmental health. This study investigated the changes in the bacterial community over the course of treatment by a WWTP. Illumina MiSeq high-throughput sequencing was performed to characterize the bacterial communities in the WWTP. This study found that potential pathogens in the WWTP, especially the genera Arcobacter and Acinetobacter, were greatly reduced. In addition, high chemical oxygen demand levels provided excessive growth substrates for the genera Hyphomicrobium and Rhodoplanes, the abundance of which could exceed autotrophic bacteria, increasing the ammonium removal. According to the network analysis, the bacterial assemblage was not randomly arranged in the WWTP, and various defined processes led to higher intra-phylum (such as Proteobacteria) coexistence than expected. Moreover, the metabolic functions of bacterial communities significantly improved in the WWTP compared with the influent. Together, the data in this study emphasize the need to understand the bacterial community of WWTPs better. When analyzing the risks of WWTP drainage systems to the environment and human health, these data should be considered.

Casimicrobium huifangae gen. nov., sp. nov., a Ubiquitous "Most-Wanted" Core Bacterial Taxon from Municipal Wastewater Treatment Plants

Microorganisms in wastewater treatment plants (WWTPs) play a key role in the removal of pollutants from municipal and industrial wastewaters. A recent study estimated that activated sludge from global municipal WWTPs harbors 1 × 10⁹ to 2 × 10⁹ microbial species, the majority of which have not yet been cultivated, and 28 core taxa were identified as "most-wanted" ones (L. Wu, D. Ning, B. Zhang, Y. Li, et al., Nat Microbiol 4:1183-1195, 2019, https://doi.org/10.1038/s41564-019-0426-5). Cultivation and characterization of the "most-wanted" core bacteria are critical to understand their genetic, physiological, phylogenetic, and ecological traits, as well as to improve the performance of WWTPs. In this study, we isolated a bacterial strain, designated SJ-1, that represents a novel cluster within Betaproteobacteria and corresponds to OTU_16 within the 28 core taxa in the "most-wanted" list. Strain SJ-1 was identified and nominated as Casimicrobium huifangae gen. nov., sp. nov., of a novel family, Casimicrobiaceae. C. huifangae is ubiquitously distributed and is metabolically versatile. In addition to mineralizing various carbon sources (including carbohydrates, aromatic compounds, and short-chain fatty acids), C. huifangae is capable of nitrate reduction and phosphorus accumulation. The population of C. huifangae accounted for more than 1% of the bacterial population of the activated sludge microbiome from the Qinghe WWTP, which showed seasonal dynamic changes. Cooccurrence analysis suggested that C. huifangae was an important module hub in the bacterial network of Qinghe WWTP.IMPORTANCE The activated sludge process is the most widely applied biotechnology and is one of the best ecosystems to address microbial ecological principles. Yet, the cultivation of core bacteria and the exploration of their physiology and ecology are limited. In this study, the core and novel bacterial taxon C. huifangae was cultivated and characterized. This study revealed that C. huifangae functioned as an important module hub in the activated sludge microbiome, and it potentially plays an important role in municipal wastewater treatment plants.

Characterizing community dynamics and exploring bacterial assemblages in two activated sludge systems

Bacterial communities in the activated sludge (AS) determine the wastewater treatment performance in the municipal wastewater treatment plants (WWTPs). Aiming at identifying the affecting factors and the variation patterns of the bacterial assemblages in AS, a 2-year time-series AS samples were collected from two separated WWTPs and metagenomic sequencing was conducted. Obvious seasonal shift and succession of the bacterial community were observed in both WWTPs on the genus and species levels, especially for the persistent taxa, implying that temperature was a decisive factor for maintaining bacterial assemblage patterns in long-term period. Taxa abundance distribution (TAD) concerning occurrence frequency and average abundance were found fitting for exponential formulations, and the approximately equal total abundance of persistent taxa suggested that stable and high abundance (~ 90%) of core functional bacterial groups would help to maintain wastewater treatment performance. Drastic changes of environmental factors were found causing temporally significant bacterial structure variation, while the innate correlations between bacterial species could recover the community gradually and maintain relative stability of the AS system. Delayed correlations between environmental factors and abundant (persistent or intermittent) bacterial species were observed widely, while synchronous biotic interactions were identified more frequently. Besides, bacterial species with similar functions were prone to cluster together and shared the same seasonal variation pattern, implicating that the cooperation of functional correlated taxa played the most dominant role in shaping the bacterial assemblages. Furthermore, rare bacterial groups were to be explored for removing emerging pollutants with lower concentrations. The results of this study would assist dealing with operational defect and optimize the treatment system in WWTPs.

Insights into the effects of carbon source on sequencing batch reactors: Performance, quorum sensing and microbial community

Effects of carbon source on the performance, quorum sensing (QS) and microbial communities in the sequencing batch reactors were investigated in this work. Among the chosen carbon source, sodium acetate (R1), glucose (R2), starch (R3) and Tween 80 (R4), sodium acetate was the best carbon source for nutrient removal, while starch was favorable for inducing the sludge bulking, and Tween 80 was beneficial to the production of extracellular polymeric substances (EPS) and proliferation of Microthrix parvicella. Additionally, the R² value of linear correlation between sludge settleability and particle size in four reactors followed an order of R1 > R2 > R3 > R4. Moreover, Person correlation analysis showed that various significant correlations were observed in reactors fed with different carbon sources and the QS mainly mediated the production and component of EPS. High-throughput sequencing analysis revealed that the carbon source affected microbial communities and the Canonical correspondence analysis results indicated that QS related to microbial communities. It was inferred that the interactions between microbial communities and QS affected system performance.

Effect of low-intensity direct current electric field on microbial nitrate removal in coal pyrolysis wastewater with low COD to nitrogen ratio

The coupling of bioelectrochemical system with the biological denitrification process (R1) was constructed for nitrate removal in coal pyrolysis wastewater (CPW) and the effect of low-intensity direct current electric field was investigated. Compared with control reactor (R2), the total nitrogen (TN) removal efficiency in R1 at the optimized voltage of 0.9 V was 94.20 ± 2.14%, which was 14.42% higher than that in R2. The biofilm on the cathode of R1 enhanced the nitrate reducing, however, nitrite was only reduced by bacteria in suspended activated sludge, which was confirmed by cyclic voltammetry measurement (CV). Microbial community network analysis revealed that exoelectrogenic bacteria of Pseudomonas was enriched on the anode of R1, and the "small world", including Zoogloea, Pseudomonas and Arenimonas, was established under the stimulation of voltage. Therefore, Pseudomonas transferred electron to anode, and Arenimonas could utilize electron from anode to reduce nitrate, which enhanced TN removal in R1.

Bacterial assembly in the bio-cake of membrane bioreactors: Stochastic vs. deterministic processes

Much about assembly processes dictating bio-cake microbiota remains uncertain, leading to poor understanding of membrane biofouling in membrane bioreactors (MBRs). This work aimed to reveal the underlying mechanisms driving bio-cake community during the biofouling process under different flux conditions. On the basis of 16S rRNA sequences, the results showed that bacterial diversity decreased with increasing fouling. Additionally, low-flux bio-cake (8 LMH) communities harbored much lower diversity than high-flux (16 LMH) bio-cake microbiomes. Ecological null model analyses and phylogenetic molecular ecological networks (pMENs) revealed that environmental filtering deterministically governed low-flux bio-cake communities. In contrast, high-flux bio-cake communities were mainly shaped in a stochastic manner. This is likely due to the higher stochastic deposition of bacterial taxa from bulk sludge because of the presence of a stronger drag force. Moreover, by lowering the flux, the interactions between bacterial lineages were enhanced; this is evidenced by the greater number of links, the higher average degree, and the higher average clustering coefficients within the pMENs in low-flux bio-cakes than those in high-flux bio-cakes. Most keystone fouling-related taxa in low-flux bio-cakes were motile and involved in nitrate reduction and polysaccharide/protein metabolism. This corroborated the important role of environmental filtering in the assembly process dictating low-flux bio-cake formation. Some low-abundance taxa were observed to be key fouling-related bacteria under both flux conditions, indicating that a few populations play paramount ecological roles in triggering biofouling. In summary, our findings clearly indicate distinct bio-cake community assembly patterns under different operational conditions and highlight the importance of developing specialized strategies for fouling control in individual MBR systems.

High variations of methanogenic microorganisms drive full-scale anaerobic digestion process

Anaerobic digestion is one of the most successful waste management strategies worldwide, wherein microorganisms play an essential role in reducing organic pollutants and producing renewable energy. However, variations of microbial community in full-scale anaerobic digesters, particularly functional groups relevant to biogas production, remain elusive. Here, we examined microbial community in a year-long monthly time series of 3 full-scale anaerobic digesters. We observed substantial diversification in community composition, with only a few abundant OTUs (e.g. Clostridiales, Anaerolineaceae and Methanosaeta) persistently present across different samples. Similarly, there were high variations in relative abundance of methanogenic archaea and methanogenic genes, which were positively correlated (r² = 0.530, P < 0.001). Variations of methanogens explained 55.7% of biogas producing rates, much higher than the explanatory percentage of environmental parameters (16.4%). Hydrogenotrophic methanogens, especially abundant Methanomicrobiales taxa, were correlated with biogas production performance (r = 0.665, P < 0.001) and nearly all methanogenic genes (0.430 < r < 0.735, P < 0.012). Given that methanogenic archaea or genes are well established for methanogenesis, we conclude that high variations in methanogenic traits (e.g. taxa or genes) are responsible for biogas production variations in full-scale anaerobic digesters.