Diversity and dynamics of dominant and rare bacterial taxa in replicate sequencing batch reactors operated under different solids retention time

Influence of high-load shocks on achieving mainstream partial nitrification: Microbial community succession

Driving microbial community succession through the regulation of operational strategies is crucial for achieving partial nitrification (PN) in municipal wastewater. However, at present, there is a decoupling between the strategic regulation of PN systems and the succession characteristics of the microbial community. This study examined the correlation between microbial community succession and PN performance under two high-load shocks (HLS1 and HLS2) treating actual sewage. During HLS1, the influent organic loading rate (OLR) and nitrogen loading rate (NLR) increased from 116.7 ± 37.7 to 219.7 ± 24.7 mg COD/(g VSS·d) and 0.21±0.02 to 0.33±0.02 kg N/m3/d respectively, with the nitrite concentration and nitrite accumulation ratio only reaching 11.7 ± 2.7 mg/L and 49.3 ± 13.9 %, respectively. During HLS2, the influent OLR and NLR increased from 123.5 ± 17.2 to 300.3 ± 49.2 mg COD/(g VSS·d) and 0.19±0.03 to 0.32±0.03 kg N/m3/d respectively, resulting in a nitrite accumulation ratio of 89.4 ± 10.7 %. The system achieved efficient PN performance and sustained for 124 days. High-throughput sequencing results showed that community diversity remained consistently high, and the community composition returned to its initial state following a minor succession during HLS1. During HLS2, the high-load shock reduced the richness and evenness of the microbial community. The community underwent succession in a new direction, leading to community composition and function changes. The results indicate that the realization, stabilization, and disruption of PN are influenced not only by operational parameters but also by microbial community structure.

Stochastic processes shape microbial community assembly in biofilters: Hidden role of rare taxa

Stochastic and deterministic processes are the major themes governing microbial community assembly; however, their roles in bioreactors are poorly understood. Herein, the mechanisms underlying microbial assembly and the effect of rare taxa were studied in biofilters. Phylogenetic tree analysis revealed differences in microbial communities at various stages. Null model analysis showed that stochastic processes shaped the community assembly, and deterministic processes emerged only in the inoculated activated sludge after domestication. This finding indicates the dominant role of stochastic factors (biofilm formation, accumulation, and aging). The Sloan neutral model corroborated the advantages of stochastic processes and mainly attributed these advantages to rare taxa. Cooccurrence networks revealed the importance of rare taxa, which accounted for more than 85% of the keystones. Overall, these results provide good foundations for understanding community assembly, especially the role of rare taxa, and offer theoretical support for future community design and reactor regulation.

Resistance of aerobic granular sludge microbiomes to periodic loss of biomass

Granular sludge is a biofilm process used for wastewater treatment which is currently being implemented worldwide. It is important to understand how disturbances affect the microbial community and performance of reactors. Here, two acetate-fed replicate reactors were inoculated with acclimatized sludge and the reactor performance, and the granular sludge microbial community succession were studied for 149 days. During this time, the microbial community was challenged by periodically removing half of the reactor biomass, subsequently increasing the food-to-microorganism (F/M) ratio. Diversity analysis together with null models show that overall, the microbial communities were resistant to the disturbances, observing some minor effects on polyphosphate-accumulating and denitrifying microbial communities and their associated reactor functions. Community turnover was driven by drift and random granule loss, and stochasticity was the governing ecological process for community assembly. These results evidence the aerobic granular sludge process as a robust system for wastewater treatment.

Influence of solids and hydraulic retention times on microbial diversity and removal of estrogens and nonylphenols in a pilot-scale activated sludge plant

The removal of EDCs in activated sludge processes can be enhanced by increasing solid and hydraulic retention times (SRT and HRT); it has been suggested that the improvement in removal is due to changes in microbial community structure (MCS). Though the influence of SRT and HRT on chemical removal and MCS has been studied in isolation, their synergistic impact on MCS and the removal of estrogens and nonylphenols in activated sludge remains unknown. Hence, we investigated how both parameters influence MCS in activated sludge processes and their ulterior effect on EDC removal. In our study, an activated sludge pilot-plant was fed with domestic sewage fortified with 100 and 1000 ng/L nonylphenols or 2 and 15 ng/L estrogens and operated at 3, 10 and 27 d SRT (constant HRT) and at 8, 16 and 24 h HRT (constant SRT). The MCS was assessed by phospholipid fatty acids (PLFA) analysis, and the archaeal and bacterial diversities were determined by 16S rRNA analysis. From the PLFA, the microbial abundance ranked as follows: Gram-negative > fungi > Gram-positive > actinomycetes whilst 16S rRNA analysis revealed Proteobacteria > Bacteroidetes > Others. Both PLFA and 16S rRNA analysis detected changes in MCS as SRT and HRT were increased. An SRT increment from 3 to 10 d resulted in higher estrone (E1) removal from 19 to 93% and nonylphenol-4-exthoxylate (NP4EO) from 44 to 73%. These findings demonstrate that EDC-removal in activated sludge plants can be optimised where longer SRT (>10 d) and HRT (>8 h) are suitable. We have also demonstrated that PLFA can be used for routine monitoring of changes in MCS in activated sludge plants.

Deciphering community assembly and succession in sequencing batch moving bed biofilm reactor: Differentiation between attached and suspended communities

Elucidating community assembly and succession is crucial to understanding the ecosystem functioning. Herein, the ecological processes underpinning community assembly and succession were studied to uncover the respective ecological functions of attached biofilms and suspended biomass in a sequencing batch moving bed biofilm reactor. Compared with suspended biomass, attached biofilms presented higher relative abundances of Nitrospira (2.94 %) and Nitrosomonas (1.25 %), and contributed to 66.89 ± 11.37 % and 68.11 ± 12.72 % of nitrification and denitrification activities, respectively. The microbial source tracking result demonstrated that early formation of suspended biomass was dominated by the seeding effect of detached biofilms in the start-up period (days 0-30), while self-growth of previous suspended biomass was eventually outcompeted the seeding effect when the reactor stabilized (days 31-120). Null model and ecological network analysis further suggested distinctive ecological processes underpinning the differentiation between attached and suspended communities in the same reactor. Specifically, in the start-up period, positive interactions facilitated early formation of attached (73.84 %) and suspended communities (59.41 %), while homogenous selection (88.89 %) and homogenizing dispersal (65.71 %) governed assembly of attached and suspended communities, respectively. When the reactor stabilized, attached and suspended communities showed low composition turnover as reflected by dominant homogenizing dispersal, while they presented distinctive trends of interspecies interactions. This study sheds light on discrepant ecological processes governing community differentiation of attached biofilms and suspended biomass, which would provide ecological insights into the regulation of hybrid ecosystems.

The counteraction of anammox community to long-term nitrite stress: Crucial roles of rare subcommunity

Understanding the temporal dynamics and recovery of anammox community under nitrite stress is critical for successful application of anammox-related processes. Here, the response behaviors of anammox community were investigated to characterize the reactor performance and ecological function under varied levels of nitrite stress (changing from 0, 50, 100, 200 to 0 mg-N/L) across a large temporal scale (588 days). The nitrogen removal rates decreased from 0.51 ± 0.02 to 0.16 ± 0.04 kg-N/(m³·d) under nitrite stress from 0 to 200 mg-N/L, while it was recovered to 0.29 ± 0.06 kg-N/(m³·d) as nitrite stress terminated. A strong community succession was driven by the initial nitrite stress of 50 mg-N/L, while the community dissimilarity mainly resulted from the increased beta diversity of rare subcommunity. Meanwhile, the rare subcommunity with high functional redundancy likely warranted the functional resilience of anammox community across the nitrite stress gradients. Moreover, the increased positive interactions between anammox bacteria and side populations supported the resilience of anammox after discontinuing nitrite stress, which facilitated the recovery of nitrogen removal efficiency. This study deciphers the interspecies interactions and functional redundancy of rare subcommunity in shaping the robustness and resilience of anammox-related processes when treating nitrite fluctuated wastewater.

Deciphering the succession dynamics of dominant and rare genera in biofilm development process

Deciphering the succession dynamics of dominant and rare taxa is crucial to understand the stability and ecosystem functions of biofilm communities. However, the essential laws of the succession dynamics based on dominant and rare taxa were still unenlightened. Herein, we investigated the succession dynamics of dominant and rare genera in multi-species biofilms developed in flow cells fed with 10 and 40 mg-TOC/L LB broth. The relative abundance of dominant genera (Enterobacteria and Acinetobacter) decreased remarkably (from 94.63% to 73.22%) in 10 mg-TOC/L LB broth, whereas they kept relatively steady (93.75 ± 4.23%) along with the cultivation time in 40 mg-TOC/L LB broth. Fluorescence in situ hybridization showed that rare genera tended to form clusters at both concentrations, while weaker dispersal of dominant genera caused patchier biofilm structures in 10 mg-TOC/L LB broth compared to that in 40 mg-TOC/L LB broth. Null model analyses further demonstrated that the stochastic ecological drift was more pronounced in the community assembly of biofilms in 10 mg-TOC/L LB broth (73.33%) than those in 40 mg-TOC/L LB broth (60.95%), weakening the competitive superiority of dominant taxa in the patchier biofilms. In addition, the co-occurrence network reflected that the positive interactions among rare genera contributed to exclude dominant genera in 10 mg-TOC/L LB broth, whereas negative interactions only occurred between the dominant Enterobacter and Acinetobacter or rare Comamonas in 40 mg-TOC/L LB broth. This study highlighted the distinctive succession dynamics of dominant and rare genera in biofilms at different substrate concentrations, which would advance our understanding of the biofilm communities in biofilm-related process.

Press Disturbance Alters Community Structure and Assembly Mechanisms of Bacterial Taxa and Functional Genes in Mesocosm-Scale Bioreactors

Press disturbances are of interest in microbial ecology, as they can drive microbial communities to alternative stable states. However, the effect of press disturbances in community assembly mechanisms, particularly with regard to taxa and functional genes at different levels of abundance (i.e., common and rare), remains largely unknown. Here, we tested the effect of a continuous alteration in substrate feeding scheme on the structure, function, and assembly of bacterial communities. Two sets of replicate 5-liter sequencing batch reactors were operated at two different organic carbon loads for a period of 74 days, following 53 days of acclimation after inoculation with sludge from a full-scale treatment plant. Temporal dynamics of community taxonomic and functional gene structure were derived from metagenomics and 16S rRNA gene metabarcoding data. Disturbed reactors exhibited different community function, structure, and assembly compared to undisturbed reactors. Bacterial taxa and functional genes showed dissimilar α-diversity and community assembly patterns. Deterministic assembly mechanisms were generally stronger in disturbed reactors and in common fractions compared to rare ones. Function quickly recovered after the disturbance was removed, but community structure did not. Our results highlight that functional gene data from metagenomics can indicate patterns of community assembly that differ from those obtained from taxon data. This study reveals how a joint evaluation of assembly mechanisms and community structure of bacterial taxa and functional genes as well as ecosystem function can unravel the response of complex microbial systems to a press disturbance.IMPORTANCE Ecosystem management must be viewed in the context of increasing frequencies and magnitudes of various disturbances that occur at different scales. This work provides a glimpse of the changes in assembly mechanisms found in microbial communities exposed to sustained changes in their environment. These mechanisms, deterministic or stochastic, can cause communities to reach a similar or variable composition and function. For a comprehensive view, we use a joint evaluation of temporal dynamics in assembly mechanisms and community structure for both bacterial taxa and their functional genes at different abundance levels, in both disturbed and undisturbed states. We further reverted the disturbance state to contrast recovery of function with community structure. Our findings are relevant, as very few studies have employed such an approach, while there is a need to assess the relative importance of assembly mechanisms for microbial communities across different spatial and temporal scales, environmental gradients, and types of disturbance.

Comparison of bacterial community structure and function under different petroleum hydrocarbon degradation conditions

Bioremediation methods have been successfully applied to the removal of organic pollutants for decades, but the responses of the microbial community to environmental factors remain less well known. In this work, the degradation rates of petroleum hydrocarbons (PHs) reached up to 50.11% ± 2.74% after optimizing the degradation conditions. Under the influence of the optimized degradation conditions, the diversity of the bacterial community gradually increased. Meanwhile, the dominant bacterial genera, encompassing Burkholderia-Paraburkholderia, Luteibacter, and Acinetobacter, remained stable. Moreover, statistical analysis indicated that the genera Bacterium, Burkholderia-Paraburkholderia, Luteibacter, and Acinetobacter contributed the most to PHs degradation. Additionally, the functional modules of amino acid metabolism, carbohydrate metabolism, as well as global and overview maps played a vital role in the metabolization of PHs. Therefore, understanding the changes of the microbial community structure and function can provide valuable guidance to further improve the degradation rate of organic waste via bioremediation methods.

Importance of Species Sorting and Immigration on the Bacterial Assembly of Different-Sized Aggregates in a Full-Scale Aerobic Granular Sludge Plant

In aerobic granular sludge (AGS) systems, different-sized microbial aggregates having different solids retention time (SRT) coexist in the same reactor compartment and are subjected to the same influent wastewater. Thus, the AGS system provides a unique ecosystem to study the importance of local (species sorting) and regional (immigration) processes in bacterial community assembly. The microbial communities of different-sized aggregates (flocs <0.2 mm, small granules (0.2-1.0 mm) and large granules >1.0 mm), influent wastewater, excess sludge and effluent of a full-scale AGS plant were characterized over a steady-state operation period of 6 months. Amplicon sequencing was integrated with mass balance to determine the SRT and net growth rate of operational taxonomic units (OTUs). We found strong evidence of species sorting as opposed to immigration, which was significantly higher at short SRT (i.e., flocs and small granules) than that at long SRT (large granules). Rare OTUs in wastewater belonging to putative functional groups responsible for nitrogen and phosphorus removal were progressively enriched with an increase in microbial aggregates size. In contrast, fecal- and sewage infrastructure-derived microbes progressively decreased in relative abundance with increase in microbial aggregate size. These findings highlight the importance of AGS as a unique model ecosystem to study fundamental microbial ecology concepts.

Combined Deterministic and Stochastic Processes Control Microbial Succession in Replicate Granular Biofilm Reactors

Granular sludge is an efficient and compact biofilm process for wastewater treatment. However, the ecological factors involved in microbial community assembly during the granular biofilm formation are poorly understood, and little is known about the reproducibility of the process. Here, three replicate bioreactors were used to investigate microbial succession during the formation of granular biofilms. We identified three successional phases. During the initial phase, the successional turnover was high and α-diversity decreased as a result of the selection of taxa adapted to grow on acetate and form aggregates. Despite these dynamic changes, the microbial communities in the replicate reactors were similar. The second successional phase occurred when the settling time was rapidly decreased to selectively retain granules in the reactors. The influence of stochasticity on succession increased and new niches were created as granules emerged, resulting in temporarily increased α-diversity. The third successional phase occurred when the settling time was kept stable and granules dominated the biomass. Turnover was low, and selection resulted in the same abundant taxa in the reactors, but drift, which mostly affected low-abundant community members, caused the community in one reactor to diverge from the other two. Even so, performance was stable and similar between reactors.

Impact of solids residence time on community structure and nutrient dynamics of mixed phototrophic wastewater treatment systems

Suspended growth, mixed community phototrophic wastewater treatment systems (including high-rate algal ponds and photobioreactors) have the potential to achieve biological nitrogen and phosphorus recovery with effluent nutrient concentrations below the current limit-of-technology. In order to achieve reliable and predictive performance, it is necessary to establish a thorough understanding of how design and operational decisions influence the complex community structure governing nutrient recovery in these systems. Solids residence time (SRT), a critical operational parameter governing growth rate, was leveraged as a selective pressure to shape microbial community structure in laboratory-scale photobioreactors fed secondary effluent from a local wastewater treatment plant. In order to decouple the effects of SRT and hydraulic retention time (HRT), nutrient loading was fixed across all experimental conditions and the effect of changing SRT on microbial community structure, diversity, and stability, as well as its impact on nutrient recovery, was characterized. Reactors were operated at distinct SRTs (5, 10, and 15 days) with diurnal lighting over long-term operation (>6 SRTs), and in-depth examination of the eukaryotic and bacterial community structure was performed using amplicon-based sequencing of the 18S and 16S rRNA genes, respectively. In order to better represent the microalgal community structure, this study leveraged improved 18S rRNA gene primers that have been shown to provide a more accurate representation of the wastewater process-relevant algal community members. Long-term operation resulted in distinct eukaryotic communities across SRTs, independent of the relative abundance of Operational Taxonomic Units (OTUs) in the inoculum. The longest SRT (15 days, SRT 15) resulted in a more stable algal community along with stable bacterial nitrification, while the shortest SRT (5 days, SRT 5) resulted in a less stable, more dynamic community. Although SRT was not strongly associated with overall bacterial diversity, the eukaryotic community of SRT 15 was significantly less diverse and less even than SRT 5, with a few dominant OTUs making up a majority of the eukaryotic community structure in the former. Overall, although longer SRTs promote stable bacterial nitrification, short SRTs promote higher eukaryotic diversity, increased functional stability, and better total N removal via biomass assimilation. These results indicate that SRT may be a key factor in not only controlling microalgal community membership, but community diversity and functional stability as well. Ultimately, the efficacy and reliability of NH₄⁺ removal may be in tension with TN removal in mixed phototrophic systems given that lower SRTs may achieve better total N removal (via biomass assimilation) through increased eukaryotic diversity, biomass productivity, and functional stability.

Microbial residence time is a controlling parameter of the taxonomic composition and functional profile of microbial communities

A remaining challenge within microbial ecology is to understand the determinants of richness and diversity observed in environmental microbial communities. In a range of systems, including activated sludge bioreactors, the microbial residence time (MRT) has been previously shown to shape the microbial community composition. However, the physiological and ecological mechanisms driving this influence have remained unclear. Here, this relationship is explored by analyzing an activated sludge system fed with municipal wastewater. Using a model designed in this study based on Monod-growth kinetics, longer MRTs were shown to increase the range of growth parameters that enable persistence, resulting in increased richness and diversity in the modeled community. In laboratory experiments, six sequencing batch reactors treating domestic wastewater were operated in parallel at MRTs between 1 and 15 days. The communities were characterized using both 16S ribosomal RNA and non-target messenger RNA sequencing (metatranscriptomic analysis), and model-predicted monotonic increases in richness were confirmed in both profiles. Accordingly, taxonomic Shannon diversity also increased with MRT. In contrast, the diversity in enzyme class annotations resulting from the metatranscriptomic analysis displayed a non-monotonic trend over the MRT gradient. Disproportionately high abundances of transcripts encoding for rarer enzymes occur at longer MRTs and lead to the disconnect between taxonomic and functional diversity profiles.

Community structure of rare methanogenic archaea: insight from a single functional group

The rare biosphere, the low abundant microbial populations, is suggested to be a conserved way of microbial life. Here we conducted a molecular survey of rare methanogenic archaea in the environment targeting the mcrA gene in order to test if general concepts associated with the structure of the rare bacterial biosphere also apply to single functional groups. Similar to what is known about rare bacterial communities, the contribution of rare methanogens to the alpha diversity is much larger than to Bray-Curtis measures. Moreover, a similar core group of methanogens harbored by the abundant and rare communities suggests similar sources and environmental controls of both groups. Among the communities of different levels of rarity, the conditionally rare methanogenic taxa largely account for the overall community dynamics of the rare biosphere and likely enter the dominant community under favorable environmental conditions. In addition, we observed a positive correlation between the alpha diversity and the production of methane when the rare taxa were taken into account. This supports the concept that increasing microbial biodiversity enhances ecological function. The composition and environmental associations of the rare methanogenic biosphere allow us to conclude that rarity is a conserved way also for single functional groups.

Membrane biofilm communities in full-scale membrane bioreactors are not randomly assembled and consist of a core microbiome

Finding efficient biofouling control strategies requires a better understanding of the microbial ecology of membrane biofilm communities in membrane bioreactors (MBRs). Studies that characterized the membrane biofilm communities in lab-and pilot-scale MBRs are numerous, yet similar studies in full-scale MBRs are limited. Also, most of these studies have characterized the mature biofilm communities with very few studies addressing early biofilm communities. In this study, five full-scale MBRs located in Seattle (Washington, U.S.A.) were selected to address two questions concerning membrane biofilm communities (early and mature): (i) Is the assembly of biofilm communities (early and mature) the result of random immigration of species from the source community (i.e. activated sludge)? and (ii) Is there a core membrane biofilm community in full-scale MBRs? Membrane biofilm (early and mature) and activated sludge (AS) samples were collected from the five MBRs, and 16S rRNA gene sequencing was applied to investigate the bacterial communities of AS and membrane biofilms (early and mature). Alpha and beta diversity measures revealed clear differences in the bacterial community structure between the AS and biofilm (early and mature) samples in the five full-scale MBRs. These differences were mainly due to the presence of large number of unique but rare operational taxonomic units (∼13% of total reads in each MBR) in each sample. In contrast, a high percentage (∼87% of total reads in each MBR) of sequence reads was shared between AS and biofilm samples in each MBR, and these shared sequence reads mainly belong to the dominant taxa in these samples. Despite the large fraction of shared sequence reads between AS and biofilm samples, simulated biofilm communities from random sampling of the respective AS community revealed that biofilm communities differed significantly from the random assemblages (P < 0.001 for each MBR), indicating that the biofilm communities (early and mature) are unlikely to represent a random sample of the AS community. In addition to the presence of unique operational taxonomic units in each biofilm sample (early or mature), comparative analysis of operational taxonomic units and genera revealed the presence of a core biofilm community in the five full-scale MBRs. These findings provided insight into the membrane biofilm communities in full-scale MBRs. More comparative studies are needed in the future to elucidate the factors shaping the core and unique biofilm communities in full-scale MBRs.

High-rate activated sludge communities have a distinctly different structure compared to low-rate sludge communities, and are less sensitive towards environmental and operational variables

High-rate activated sludge processes allow for the recovery of organics and energy from wastewaters. These systems are operated at a short sludge retention time and high sludge-specific loading rates, which results in a higher sludge yield and better digestibility than conventional, low-rate activated sludge. Little is known about the microbial ecology of high-rate systems. In this work, we address the need for a fundamental understanding of how high-rate microbial communities differ from low-rate communities. We investigated the high-rate and low-rate communities in a sewage treatment plant in relation to environmental and operational variables over a period of ten months. We demonstrated that (1) high-rate and low-rate communities are distinctly different in terms of richness, evenness and composition, (2) high-rate community dynamics are more variable and less shaped by deterministic factors compared to low-rate communities, (3) sub-communities of continuously core and transitional members are more shaped by deterministic factors than the continuously rare members, both in high-rate and low-rate communities, and (4) high-rate community members showed a co-occurrence pattern similar to that of low-rate community members, but were less likely to be correlated to environmental and operational variables. These findings provide a basis for further optimization of high-rate systems, in order to facilitate resource recovery from wastewater.