Microbial biogeography across a full-scale wastewater treatment plant transect: evidence for immigration between coupled processes

Validation of 16S rRNA gene sequencing and metagenomics for evaluating microbial immigration in a methanogenic bioreactor

To quantitatively evaluate the impact of microbial immigration from an upstream community on the microbial assembly of a downstream community, an ecological genomics (ecogenomics)-based mass balance (EGMB) model coupled with 16S rRNA gene sequencing was previously developed. In this study, a mock community was used to further validate the EGMB models and demonstrate the feasibility of using metagenome-based EGMB model to reveal both microbial activity and function. The mock community consisting of Aeromonas, Escherichia, and Pseudomonas was fed into a lab-scale methanogenic bioreactor together with dissolved organic substrate. Using qPCR, 16S rRNA gene, 16S rRNA gene copy number normalization (GCN), and metagenome, results showed highly comparable community profiles in the feed. In the bioreactor, Aeromonas and Pseudomonas exhibited negative growth rates throughout the experiment by all approaches. Escherichia's growth rate was negative by most biomarkers but was slightly positive by 16S rRNA gene. Still, all approaches showed a decreasing trend toward negative in the growth rate of Escherichia as reactor operation time increased. Uncultivated populations of phyla Desulfobacterota, Chloroflexi, Actinobacteriota, and Spirochaetota were observed to increase in abundance, suggesting their contribution in degrading the feed biomass. Based on metabolic reconstruction of metagenomes, these populations possessed functions of hydrolysis, fermentation, fatty acid degradation, or acetate oxidation. Overall results supported the application of both 16S rRNA gene- and metagenome-based EGMB models to measure the growth rate of microbes in the bioreactor, and the latter had advantage in providing insights into the microbial functions of uncultivated populations.

Disturbance-based management of ecosystem services and disservices in partial nitritation-anammox biofilms

The resistance and resilience provided by functional redundancy, a common feature of microbial communities, is not always advantageous. An example is nitrite oxidation in partial nitritation-anammox (PNA) reactors designed for nitrogen removal in wastewater treatment, where suppression of nitrite oxidizers like Nitrospira is sought. In these ecosystems, biofilms provide microhabitats with oxygen gradients, allowing the coexistence of aerobic and anaerobic bacteria. We designed a disturbance experiment where PNA biofilms, treating water from a high-rate activated sludge process, were constantly or intermittently exposed to anaerobic sidestream wastewater, which has been proposed to inhibit nitrite oxidizers. With increasing sidestream exposure we observed decreased abundance, alpha-diversity, functional versatility, and hence functional redundancy, among Nitrospira in the PNA biofilms, while the opposite patterns were observed for anammox bacteria within Brocadia. At the same time, species turnover was observed for aerobic ammonia-oxidizing Nitrosomonas populations. The different exposure regimens were associated with metagenomic assembled genomes of Nitrosomonas, Nitrospira, and Brocadia, encoding genes related to N-cycling, substrate usage, and osmotic stress response, possibly explaining the three different patterns by niche differentiation. These findings imply that disturbances can be used to manage the functional redundancy of biofilm microbiomes in a desirable direction, which should be considered when designing operational strategies for wastewater treatment.

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.

Layered viscoelastic properties of granular biofilms

Granular biofilms are dense spherical complex biological systems composed mainly of multi-microbial cells, water, and extracellular polymeric substances (EPS). They facilitate efficient purification and settling of activated sludge in wastewater treatment processes. The viscoelastic properties of these complex biofilm systems are important characteristics that control their growth and dictate how they respond to hydrodynamic forces and chemical stimuli. However, the viscoelastic properties of granular biofilms are poorly understood. In this paper, we study granular biofilms' viscoelastic properties using optical coherence elastography (OCE), a nondestructive method that integrates optical coherence tomography (OCT) with elastic wave propagation. While quantitative viscoelastic characterization of granular biofilms is challenging due to their heterogeneous properties, we show that elastic waves are suitable for this purpose. First, we employ guided elastic waves in a thin section of a granular biofilm to reveal a two-layered profile for the viscoelastic properties. Next, we utilize circumferential elastic waves that propagate near the surface of a non-sectioned spherical biofilm to quantify the layered system's viscoelastic properties. To the best of our knowledge, this work is the first quantitative study that characterizes the layered viscoelastic properties of granular biofilms. The measurement approach may provide a platform to study the interplay between the viscoelastic properties and other characteristics of granular biofilms such as the complex microbial system, morphology, and oxygen distribution.

Mass-immigration determines the assembly of activated sludge microbial communities

Wastewater treatment plants are engineering technologies used worldwide to protect the environment and human health. Microbial communities sustain these plants, so it is crucial to know the key factors responsible for the community assembly. We show, in contrast to existing understanding, that microbial immigration largely controls the community structure in these plants and that the fate (growth or death) of immigrating species in the plants is controlled by local factors. The community structure was quantitatively predicted by the immigrating microbial community, highlighting the need to revise the way we today understand, design, and manage microbial communities in wastewater treatment plants.

The assembly of bacterial communities in wastewater treatment plants (WWTPs) is affected by immigration via wastewater streams, but the impact and extent of bacterial immigrants are still unknown. Here, we quantify the effect of immigration at the species level in 11 Danish full-scale activated sludge (AS) plants. All plants have different source communities but have very similar process design, defining the same overall environmental growth conditions. The AS community composition in each plant was strongly reflected by the corresponding influent wastewater (IWW) microbial composition. Most species in AS across the plants were detected and quantified in the corresponding IWW, allowing us to identify their fate in the AS: growing, disappearing, or surviving. Most of the abundant species in IWW disappeared in AS, so their presence in the AS biomass was only due to continuous mass-immigration. In AS, most of the abundant growing species were present in the IWW at very low abundances. We predicted the AS species abundances from their abundance in IWW by using a partial least square regression model. Some species in AS were predicted by their own abundance in IWW, while others by multiple species abundances. Detailed analyses of functional guilds revealed different prediction patterns for different species. We show, in contrast to the present understanding, that the AS microbial communities were strongly controlled by the IWW source community and could be quantitatively predicted by taking into account immigration. This highlights a need to revise the way we understand, design, and manage the microbial communities in WWTPs.

Urban wastewater bacterial communities assemble into seasonal steady states

BACKGROUND

Microorganisms in urban sanitary sewers exhibit community properties that suggest sewers are a novel ecosystem. Sewer microorganisms present both an opportunity as a control point for wastewater treatment and a risk to human health. If treatment processes are to be improved and health risks quantified, then it is necessary to understand microbial distributions and dynamics within this community. Here, we use 16S rRNA gene sequencing to characterize raw influent wastewater bacterial communities in a 5-year time series from two wastewater treatment plants in Milwaukee, WI; influent wastewater from 77 treatment plants across the USA; and wastewater in 12 Milwaukee residential sewers.

RESULTS

In Milwaukee, we find that in transit from residences to treatment plants, the human bacterial component of wastewater decreases in proportion and exhibits stochastic temporal variation. In contrast, the resident sewer community increases in abundance during transit and cycles seasonally according to changes in wastewater temperature. The result is a bacterial community that assembles into two distinct community states each year according to the extremes in wastewater temperature. Wastewater bacterial communities from other northern US cities follow temporal trends that mirror those in Milwaukee, but southern US cities have distinct community compositions and differ in their seasonal patterns.

CONCLUSIONS

Our findings provide evidence that environmental conditions associated with seasonal change and climatic differences related to geography predictably structure the bacterial communities residing in below-ground sewer pipes. Video abstract.

Core activated sludge communities are influenced little by immigration: Case study of a membrane bioreactor plant

Microbial immigrants arriving with influent wastewater may influence activated sludge (AS) ecosystems. However, the extent to which immigration impacts AS communities is still debated. To explore the intensity of immigration impact, we used sequencing technology to track the raw wastewater and AS communities from a membrane bioreactor plant over a 12-month period. We first distinguished core populations from peripheral ones in both raw wastewater and AS based on their occurrence frequency and abundance. The results showed that core OTUs (≥ 80% occurrence frequency) made up a large fraction (> 90%) of total sequences, while peripheral OTUs composed the majority of all detected OTUs but merely occupied a few sequences. A significant difference in core communities between the influent and AS was found, as well as between the compositions of core and peripheral populations. Additionally, the persistent functional bacteria of AS, although not numerically dominant, accounted for 96.24% of the total sequences related to nutrient turnover, suggesting the presence of a small number of longstanding and core functional bacteria in the AS ecosystem. Importantly, 64% of the 5188 OTUs in AS, which accounted for 91.51% of the sequences, exhibited positive growth rates, which suggested that their apparent abundances were due to growth within the plant, not from immigration. Taken together, these results demonstrated that the impact of influent populations on core AS communities was limited. Overall, this work provides quantitative insights into the impact of immigration, which is expected to advance our understanding of the AS community assembly.

A framework for P-cycle assessment in wastewater treatment plants

Phosphorus (P) in wastewater has a variety of negative effects and is usually permanently lost as a non-renewable resource. To mitigate future P shortage, P must be recovered from wastewater, preferably by bio-based technologies to avoid toxic side streams. A standardized procedure for the determination of P types and P concentrations in all liquid and solid process stages was established, which is applicable to all full-scale wastewater treatment plants (WWTPs). Based on this, an equally universal calculation framework for P-cycle assessment based on volume flow and mass load rates was designed to identify the most promising process streams for biological P recovery. As an example, in 16 process streams of a typical WWTP, concentrations of free, bound and total P were calculated and microbial communities were analyzed by flow cytometry over 748 days. The most promising process streams for the recovery of free P were anaerobic digester sludge, centrate and the water-extracts of the biosolids with 0.510 kg P m^-3, 0.075 kg P m^-3 and 1.023 kg P m^-3, while the best process streams for the recovery of bound P were return sludge, excess sludge, anaerobic digester sludge, and the solids of the biosolids with 0.300 kg P m^-3, 0.268 kg P m^-3, 0.213 kg P m^-3 and 1.336 kg P m^-3, respectively. Microorganisms capable of P accumulation were active in all process stages and it was observed that chemical P precipitation antagonizes biological P removal. The framework for P-cycle assessment was able to identify process streams that are economically viable to make future in-stream technologies for biological P removal feasible.

A conceptual framework modeling of functional microbial communities in wastewater treatment electro-bioreactors

Understanding the microbial ecology of a system allows linking members of the community and their metabolic functions to the performance of the wastewater bioreactor. This study provided a comprehensive conceptual framework for microbial communities in wastewater treatment electro-bioreactors (EBRs). The model was based on data acquired from monitoring the effect of altering different bioreactor operational parameters, such as current density and hydraulic retention time, on the microbial communities of an EBR and its nutrient removal efficiency. The model was also based on the 16S rRNA gene high-throughput sequencing data analysis and bioreactor efficiency data. The collective data clearly demonstrated that applying various electric currents affected the microbial community composition and stability and the reactor efficiency in terms of chemical oxygen demand, N and P removals. Moreover, a schematic that recommends operating conditions that are tailored to the type of wastewater that needs to be treated based on the functional microbial communities enriched at specific operating conditions was suggested. In this study, a conceptual model as a simplified representation of the behavior of microbial communities in EBRs was developed. The proposed conceptual model can be used to predict how biological treatment of wastewater in EBRs can be improved by varying several operating conditions.

Metagenomics for taxonomy profiling: tools and approaches

The study of metagenomics is an emerging field that identifies the total genetic materials in an organism along with the set of all genetic materials like deoxyribonucleic acid and ribose nucleic acid, which play a key role with the maintenance of cellular functions. The best part of this technology is that it gives more flexibility to environmental microbiologists to instantly pioneer the immense genetic variability of microbial communities. However, it is intensively complex to identify the suitable sequencing measures of any specific gene that can exclusively indicate the involvement of microbial metagenomes and be able to advance valuable results about these communities. This review provides an overview of the metagenomic advancement that has been advantageous for aggregation of more knowledge about specific genes, microbial communities and its metabolic pathways. More specific drawbacks of metagenomes technology mainly depend on sequence-based analysis. Therefore, this 'targeted based metagenomics' approach will give comprehensive knowledge about the ecological, evolutionary and functional sequence of significantly important genes that naturally exist in living beings either human, animal and microorganisms from distinctive ecosystems.

Insitu kinetics of human pharmaceutical conjugates and the impact of transformation, deconjugation, and sorption on persistence in wastewater batch bioreactors

The fate of selected common pharmaceuticals and four of their major conjugates in wastewater batch bioreactors was evaluated to determine how treatment plant parameters such as addition of air, and the presence of waste activated sludge (WAS) could influence the removal of parent compounds and conjugates. Under a realistic hydraulic residence time (HRT) for each treatment sub-process of approximately 2 h, acetaminophen and its sulfate metabolite were both rapidly degraded (>99%). Propranolol was sulfated and concurrently removed. Deconjugation of N-acetylsulfamethoxazole and sulfamethoxazole-glucuronide contributed to increases of the parent sulfamethoxazole. Thyroxine was resistant to degradation, while thyroxine-glucuronide was rapidly deconjugated (>90% in <2 h). In the absence of WAS, sorption to suspended solids was another major removal mechanism for acetaminophen, propranolol, sulfamethoxazole, and thyroxine. However, with WAS, concentrations associated with suspended solids decreased for all analytes within 24 h. These results indicate that both conjugation and back-transformation are compound-specific and dependent on parameters such as HRT, addition of microbial content, and suspended solids levels. Therefore, conjugation-deconjugation processes may strongly influence the speciation of pharmaceuticals and their fate in wastewater treatment plant effluents.

The seasonal changes of core bacterial community decide sewage purification in sub-plateau municipal sewage treatment plants

The decline of sewage purification efficiency in winter is a frequent problem in sub-plateau municipal sewage treatment plants (MSTPs). Understanding the links between activated sludge (AS) bacterial community and sewage purification is crucial for exploring the cause of this problem. In this study, Illumina high-throughput sequencing technology was applied to investigate the seasonal changes of AS bacterial community in sub-plateau MSTPs. The sequencing result indicates that the bacterial community OTU number, diversity, and relative abundance in winter are significantly lower than that in summer samples. The discriminant linear effect size analysis (LEfSe) reveals that Proteobacteria and Chloroflexi members were enriched in summer AS, while Actinobacteria and Firmicutes were enriched in winter AS. The results indicate that different core bacterial community assembly was developed in summer and winter, respectively. The changes in bacterial community may be the reasons for the lower sewage purification efficiency in winter. Furthermore, redundancy analysis (RDA) shows that temperature and dissolved oxygen (DO) are the principal factors that drive the seasonal changes in the core bacterial community diversity, richness and structure in sub-plateau MSTPs. Thus, the sub-plateau AS selects for a unique community assembly pattern and shapes the particular AS ecosystem. These results expand previous understanding and provide insight into the relationship between bacterial community and performance of sub-plateau MSTPs.

Ecological Processes Affecting Long-Term Eukaryote and Prokaryote Biofilm Persistence in Nitrogen Removal from Sewage

The factors affecting long-term biofilm stability in sewage treatment remain largely unexplored. We therefore analyzed moving bed bioreactors (MBBRs) biofilm composition and function two years apart from four reactors in a nitrogen-removal sewage treatment plant. Multivariate ANOVA revealed a similar prokaryote microbiota composition on biofilm carriers from the same reactors, where reactor explained 84.6% of the variance, and year only explained 1.5%. Eukaryotes showed a less similar composition with reactor explaining 56.8% of the variance and year 9.4%. Downstream effects were also more pronounced for eukaryotes than prokaryotes. For prokaryotes, carbon source emerged as a potential factor for deterministic assembly. In the two reactors with methanol as a carbon source, the bacterial genus Methylotenera dominated, with M. versatilis as the most abundant species. M. versatilis showed large lineage diversity. The lineages mainly differed with respect to potential terminal electron acceptor usage (nitrogen oxides and oxygen). Searches in the Sequence Read Archive (SRA) database indicate a global distribution of the M. versatilis strains, with methane-containing sediments as the main habitat. Taken together, our results support long-term prokaryote biofilm persistence, while eukaryotes were less persistent.

Microbial and Viral Communities and Their Antibiotic Resistance Genes Throughout a Hospital Wastewater Treatment System

Antibiotic resistance poses a serious threat to global public health, and antibiotic resistance determinants can enter natural aquatic systems through discharge of wastewater effluents. Hospital wastewater in particular is expected to contain high abundances of antibiotic resistance genes (ARGs) compared to municipal wastewater because it contains human enteric bacteria that may include antibiotic-resistant organisms originating from hospital patients, and can also have high concentrations of antibiotics and antimicrobials relative to municipal wastewater. Viruses also play an important role in wastewater treatment systems since they can influence the bacterial community composition through killing bacteria, facilitating transduction of genetic material between organisms, and modifying the chromosomal content of bacteria as prophages. However, little is known about the fate and connections between ARGs, viruses, and their associated bacteria in hospital wastewater systems. To address this knowledge gap, we characterized the composition and persistence of ARGs, dsDNA viruses, and bacteria from influent to effluent in a pilot-scale hospital wastewater treatment system in Israel using shotgun metagenomics. Results showed that ARGs, including genes conferring resistance to antibiotics of high clinical relevance, were detected in all sampling locations throughout the pilot-scale system, with only 16% overall depletion of ARGs per genome equivalent between influent and effluent. The most common classes of ARGs detected throughout the system conferred resistance to aminoglycoside, cephalosporin, macrolide, penam, and tetracycline antibiotics. A greater proportion of total ARGs were associated with plasmid-associated genes in effluent compared to in influent. No strong associations between viral sequences and ARGs were identified in viral metagenomes from the system, suggesting that phage may not be a significant vector for ARG transfer in this system. The majority of viruses in the pilot-scale system belonged to the families Myoviridae, Podoviridae, and Siphoviridae. Gammaproteobacteria was the dominant class of bacteria harboring ARGs and the most common putative viral host in all samples, followed by Bacilli and Betaproteobacteria. In the total bacterial community, the dominant class was Betaproteobacteria for each sample. Overall, we found that a variety of different types of ARGs and viruses were persistent throughout this hospital wastewater treatment system, which can be released to the environment through effluent discharge.

Quantifying the contribution of microbial immigration in engineered water systems

Immigration is a process that can influence the assembly of microbial communities in natural and engineered environments. However, it remains challenging to quantitatively evaluate the contribution of this process to the microbial diversity and function in the receiving ecosystems. Currently used methods, i.e., counting shared microbial species, microbial source tracking, and neutral community model, rely on abundance profile to reveal the extent of overlapping between the upstream and downstream communities. Thus, they cannot suggest the quantitative contribution of immigrants to the downstream community function because activities of individual immigrants are not considered after entering the receiving environment. This limitation can be overcome by using an approach that couples a mass balance model with high-throughput DNA sequencing, i.e., ecogenomics-based mass balance. It calculates the net growth rate of individual microbial immigrants and partitions the entire community into active populations that contribute to the community function and inactive ones that carry minimal function. Linking activities of immigrants to their abundance further provides quantification of the contribution from an upstream environment to the downstream community. Considering only active populations can improve the accuracy of identifying key environmental parameters dictating process performance using methods such as machine learning.

Viral composition and context in metagenomes from biofilm and suspended growth municipal wastewater treatment plants

Wastewater treatment plants (WWTPs) contain high density and diversity of viruses which can significantly impact microbial communities in aquatic systems. While previous studies have investigated viruses in WWTP samples that have been specifically concentrated for viruses and filtered to exclude bacteria, little is known about viral communities associated with bacterial communities throughout wastewater treatment systems. Additionally, differences in viral composition between attached and suspended growth wastewater treatment bioprocesses are not well characterized. Here, shotgun metagenomics was used to analyse wastewater and biomass from transects through two full-scale WWTPs for viral composition and associations with bacterial hosts. One WWTP used a suspended growth activated sludge bioreactor and the other used a biofilm reactor (trickling filter). Myoviridae, Podoviridae and Siphoviridae were the dominant viral families throughout both WWTPs, which are all from the order Caudovirales. Beta diversity analysis of viral sequences showed that samples clustered significantly both by plant and by specific sampling location. For each WWTP, the overall bacterial community structure was significantly different than community structure of bacterial taxa associated with viral sequences. These findings highlight viral community composition in transects through different WWTPs and provide context for dsDNA viral sequences in bacterial communities from these systems.

Microbial immigration in wastewater treatment systems: analytical considerations and process implications

Microbial immigration from sewers to wastewater treatment systems is attracting increasing attention for understanding community assembly mechanisms, and improving process modeling and operation. While there is no consensus on approaches to analyze immigration, we suggest to classify them as relevant to either rare (non-observable) diffusive immigration or to time-continuous high-rate mass flow immigration (i.e. mass effects). When analyzed by a mass flow approach, heterotrophs appear to be strongly influenced by deterministic selection, suggesting that the heterotrophs should be subdivided into several functional guilds when assessing their assembly mechanisms. Conversely, nitrifiers appear to transfer neutrally from sewer to activated sludge, and this immigration can restore full nitrification in otherwise non-nitrifying reactors. With further refinement, these findings could be included in predictive process models with various objectives.

Wastewater microbial community structure and functional traits change over short timescales

Wastewater contains microorganisms coming from various sources, e.g. feces discharges, soil infiltrations and sewer biofilms and sediments. The primary objective of this work was to determine if end-of-pipe wastewater microbial community structures exhibits short-timescale variation, and assess possible microbial origins. To this end, we measured hourly physicochemical characteristics of wastewater influent for 2 days and analyzed the microbial community at 4-h intervals using 16S rRNA gene amplicon sequencing. Results showed large variations in the microbial community composition at phylum and genus levels, i.e. Proteobacteria ranged from 44 to 63% of the total relative abundance and Arcobacter ranged from 11 to 22%. Diurnal patterns were observed in the alpha-diversity, beta-diversity and the prevalence of several taxa. Wastewater physicochemical characteristics explained 61% of the total microbial community variance by Canonical Correspondence Analysis (CCA), with flow rate being the main explanatory variable exhibiting a clear diurnal profile. Comparison with public databases using closed reference OTUs revealed that only 7.3% of the sequences were shared with human gut microbiota and 21.7% with soil microbiota, the majority being from the sewer biofilms and sediments. The functional trait, weighted average ribosomal RNA operon (rrn) copy number per genome, was found to be relatively high in the wastewater microbiota (average 3.6, soil 2.1, and human gut 2.6) and significantly correlated with flow, inferring active microbial enrichments in the sewer. The prevalence of Methylophilaceae, methanol oxidation genes and denitrification genes were related to high influent methanol and NO₃^- concentration in the influent wastewater. These functional organisms and genes indicate important carbon and nutrient removal related functions in the sewer. Together, the observed temporal patterns of the microbial community and functional traits suggest that high wastewater flow causes greater transport of active sewer microorganisms which are functionally important.

Coupling growth kinetics modeling with machine learning reveals microbial immigration impacts and identifies key environmental parameters in a biological wastewater treatment process

BACKGROUND

Ubiquitous in natural and engineered ecosystems, microbial immigration is one of the mechanisms shaping community assemblage. However, quantifying immigration impact remains challenging especially at individual population level. The activities of immigrants in the receiving community are often inadequately considered, leading to potential bias in identifying the relationship between community composition and environmental parameters.

RESULTS

This study quantified microbial immigration from an upstream full-scale anaerobic reactor to downstream activated sludge reactors. A mass balance was applied to 16S rRNA gene amplicon sequencing data to calculate the net growth rates of individual populations in the activated sludge reactors. Among the 1178 observed operational taxonomic units (OTUs), 582 had a positive growth rate, including all the populations with abundance > 0.1%. These active populations collectively accounted for 99% of the total sequences in activated sludge. The remaining 596 OTUs with a growth rate ≤ 0 were classified as inactive populations. All the abundant populations in the upstream anaerobic reactor were inactive in the activated sludge process, indicating a negligible immigration impact. We used a supervised learning regressor to predict environmental parameters based on community composition and compared the prediction accuracy based on either the entire community or the active populations. Temperature was the most predictable parameter, and the prediction accuracy was improved when only active populations were used to train the regressor.

CONCLUSIONS

Calculating growth rate of individual microbial populations in the downstream system provides an effective approach to determine microbial activity and quantify immigration impact. For the studied biological process, a marginal immigration impact was observed, likely due to the significant differences in the growth environments between the upstream and downstream processes. Excluding inactive populations as a result of immigration further enhanced the prediction of key environmental parameters affecting process performance.

Reactor staging influences microbial community composition and diversity of denitrifying MBBRs- Implications on pharmaceutical removal

The subdivision of biofilm reactor in two or more stages (i.e., reactor staging) represents an option for process optimisation of biological treatment. In our previous work, we showed that the gradient of influent organic substrate availability (induced by the staging) can influence the microbial activity (i.e., denitrification and pharmaceutical biotransformation kinetics) of a denitrifying three-stage Moving Bed Biofilm Reactor (MBBR) system. However, it is unclear whether staging and thus the long-term exposure to varying organic carbon type and loading influences the microbial community structure and diversity. In this study, we investigated biofilm structure and diversity in the three-stage MBBR system (S) compared to a single-stage configuration (U) and their relationship with microbial functions. Results from 16S rRNA amplicon libraries revealed a significantly higher microbial richness in the staged MBBR (at 99% sequence similarity) compared to single-stage MBBR. A more even and diverse microbial community was selected in the last stage of S (S3), likely due to exposure to carbon limitation during continuous-flow operation. A core of OTUs was shared in both systems, consisting of Burkholderiales, Xanthomonadales, Flavobacteriales and Sphingobacteriales, while MBBR staging selected for specific taxa (i.e., Candidate division WS6 and Deinococcales). Results from quantitative PCR (qPCR) showed that S3 exhibited the lowest abundance of 16S rRNA but the highest abundance of atypical nosZ, suggesting a selection of microbes with more diverse N-metabolism (i.e., incomplete denitrifiers) in the stage exposed to the lowest carbon availability. A positive correlation (p < 0.05) was observed between removal rate constants of several pharmaceuticals with abundance of relevant denitrifying genes, but not with biodiversity. Despite the previously suggested positive relationship between microbial diversity and functionality in macrobial and microbial ecosystems, this was not observed in the current study, indicating a need to further investigate structure-function relationships for denitrifying systems.

Assessment of Microbial Community Structure and Function in Serially Passaged Wastewater Electro-Bioreactor Sludge: An Approach to Enhance Sludge Settleability

Several studies have been carried out to understand bulking phenomena and the importance of environmental factors on sludge settling characteristics. The main objective of this study was to carry out functional characterization of microbial community structure of wastewater electro-bioreactor sludge as it undergoes serial passaging in the presence or absence of a current density over 15 days. Illumina MiSeq sequencing and QIIME were used to assess sludge microbial community shifts over time. (α) and (β) diversity analysis were conducted to assess the microbial diversity in electro-bioreactors. A phylogeny-based weighted UniFrac distance analysis was used to compare between bacterial communities while BIO-ENV trend and Spearman’s rank correlation analysis were performed to investigate how reactor operational parameters correlated with bacterial community diversity. Results showed that the removal efficiency of soluble chemical oxygen demand (sCOD) ranged from 91–97%, while phosphorous (PO₄³⁻-P) removal was approximately 99%. Phylogenetic analysis revealed stark differences in the development of sludge microbial communities in the control and treatment reactor. There was no mention of any studies aimed at characterizing functional microbial communities under electric field and the results communicated here are the first, to our knowledge, that address this gap in the literature.

Microbial rRNA gene expression and co-occurrence profiles associate with biokinetics and elemental composition in full-scale anaerobic digesters

This study examined whether the abundance and expression of microbial 16S rRNA genes were associated with elemental concentrations and substrate conversion biokinetics in 20 full-scale anaerobic digesters, including seven municipal sewage sludge (SS) digesters and 13 industrial codigesters. SS digester contents had higher methane production rates from acetate, propionate and phenyl acetate compared to industrial codigesters. SS digesters and industrial codigesters were distinctly clustered based on their elemental concentrations, with higher concentrations of NH₃ -N, Cl, K and Na observed in codigesters. Amplicon sequencing of 16S rRNA genes and reverse-transcribed 16S rRNA revealed divergent grouping of microbial communities between mesophilic SS digesters, mesophilic codigesters and thermophilic digesters. Higher intradigester distances between Archaea 16S rRNA and rRNA gene profiles were observed in mesophilic codigesters, which also had the lowest acetate utilization biokinetics. Constrained ordination showed that microbial rRNA and rRNA gene profiles were significantly associated with maximum methane production rates from acetate, propionate, oleate and phenyl acetate, as well as concentrations of NH₃ -N, Fe, S, Mo and Ni. A co-occurrence network of rRNA gene expression confirmed the three main clusters of anaerobic digester communities based on active populations. Syntrophic and methanogenic taxa were highly represented within the subnetworks, indicating that obligate energy-sharing partnerships play critical roles in stabilizing the digester microbiome. Overall, these results provide new evidence showing that different feed substrates associate with different micronutrient compositions in anaerobic digesters, which in turn may influence microbial abundance, activity and function.

Bacterial community structure transformed after thermophilically composting human waste in Haiti

Recycling human waste for beneficial use has been practiced for millennia. Aerobic (thermophilic) composting of sewage sludge has been shown to reduce populations of opportunistically pathogenic bacteria and to inactivate both Ascaris eggs and culturable Escherichia coli in raw waste, but there is still a question about the fate of most fecal bacteria when raw material is composted directly. This study undertook a comprehensive microbial community analysis of composting material at various stages collected over 6 months at two composting facilities in Haiti. The fecal microbiota signal was monitored using a high-density DNA microarray (PhyloChip). Thermophilic composting altered the bacterial community structure of the starting material. Typical fecal bacteria classified in the following groups were present in at least half the starting material samples, yet were reduced below detection in finished compost: Prevotella and Erysipelotrichaceae (100% reduction of initial presence), Ruminococcaceae (98–99%), Lachnospiraceae (83–94%, primarily unclassified taxa remained), Escherichia and Shigella (100%). Opportunistic pathogens were reduced below the level of detection in the final product with the exception of Clostridium tetani, which could have survived in a spore state or been reintroduced late in the outdoor maturation process. Conversely, thermotolerant or thermophilic Actinomycetes and Firmicutes (e.g., Thermobifida, Bacillus, Geobacillus) typically found in compost increased substantially during the thermophilic stage. This community DNA-based assessment of the fate of human fecal microbiota during thermophilic composting will help optimize this process as a sanitation solution in areas where infrastructure and resources are limited.

Regional synchrony in full-scale activated sludge bioreactors due to deterministic microbial community assembly

Seasonal community structure and regionally synchronous population dynamics have been observed in natural microbial ecosystems, but have not been well documented in wastewater treatment bioreactors. Few studies of community dynamics in full-scale activated sludge systems facing similar meteorological conditions have been done to compare the importance of deterministic and neutral community assembly mechanisms. We subjected weekly activated sludge samples from six regional full-scale bioreactors at four wastewater treatment plants obtained over 1 year to Illumina sequencing of 16S ribosomal RNA genes, resulting in a library of over 17 million sequences. All samples derived from reactors treating primarily municipal wastewater. Despite variation in operational characteristics and location, communities displayed temporal synchrony at the individual operational taxonomic unit (OTU), broad phylogenetic affiliation and community-wide scale. Bioreactor communities were dominated by 134 abundant and highly regionally synchronized OTU populations that accounted for over 50% of the total reads. Non-core OTUs displayed abundance-dependent population synchrony. Alpha diversity varied by reactor, but showed a highly reproducible and synchronous seasonal fluctuation. Community similarity was dominated by seasonal changes, but individual reactors maintained minor stable differences after 1 year. Finally, the impacts of mass migration driven by direct biomass transfers between reactors was investigated, but had no significant effect on community similarity or diversity in the sink community. Our results show that population dynamics in activated sludge bioreactors are consistent with niche-driven assembly guided by seasonal temperature fluctuations.

Microbial Community Structure of Activated Sludge in Treatment Plants with Different Wastewater Compositions

Activated sludge (AS) plays a crucial role in the treatment of domestic and industrial wastewater. AS is a biocenosis of microorganisms capable of degrading various pollutants, including organic compounds, toxicants, and xenobiotics. We performed 16S rRNA gene sequencing of AS and incoming sewage in three wastewater treatment plants (WWTPs) responsible for processing sewage with different origins: municipal wastewater, slaughterhouse wastewater, and refinery sewage. In contrast to incoming wastewater, the taxonomic structure of AS biocenosis was found to become stable in time, and each WWTP demonstrated a unique taxonomic pattern. Most pathogenic microorganisms (Streptococcus, Trichococcus, etc.), which are abundantly represented in incoming sewage, were significantly decreased in AS of all WWTPs, except for the slaughterhouse wastewater. Additional load of bioreactors with influent rich in petroleum products and organic matter was associated with the increase of bacteria responsible for AS bulking and foaming. Here, we present a novel approach enabling the prediction of the metabolic potential of bacterial communities based on their taxonomic structures and MetaCyc database data. We developed a software application, XeDetect, to implement this approach. Using XeDetect, we found that the metabolic potential of the three bacterial communities clearly reflected the substrate composition. We revealed that the microorganisms responsible for AS bulking and foaming (most abundant in AS of slaughterhouse wastewater) played a leading role in the degradation of substrates such as fatty acids, amino acids, and other bioorganic compounds. Moreover, we discovered that the chemical, rather than the bacterial composition of the incoming wastewater was the main factor in AS structure formation. XeDetect (freely available: https://sourceforge.net/projects/xedetect) represents a novel powerful tool for the analysis of the metabolic capacity of bacterial communities. The tool will help to optimize bioreactor performance and avoid some most common technical problems.

Ionic liquid biodegradability depends on specific wastewater microbial consortia

Complete biodegradation of a newly-synthesized chemical in a wastewater treatment plant (WWTP) eliminates the potential for novel environmental pollutants. However, differences within- and between-WWTP microbial communities may alter expectations for biodegradation. WWTP communities can also serve as a source of unique consortia that, when enriched, can metabolize chemicals that tend to resist degradation, but are otherwise promising green alternatives. We tested the biodegradability of three ionic liquids (ILs): 1-octyl-3-methylpyridinium bromide (OMP), 1-butyl-3-methylpyridinium bromide (BMP) and 1-butyl-3-methylimidazolium chloride (BMIM). We performed tests using communities from two WWTPs at three time points. Site-specific and temporal variation both influenced community composition, which impacted the success of OMP biodegradability. Neither BMP nor BMIM degraded in any test, suggesting that these ILs are unlikely to be removed by traditional treatment. Following standard biodegradation assays, we enriched for three consortia that were capable of quickly degrading OMP, BMP and BMIM. Our results indicate WWTPs are not functionally redundant with regard to biodegradation of specific ionic liquids. However, consortia can be enriched to degrade chemicals that fail biodegradability assays. This information can be used to prepare pre-treatment procedures and prevent environmental release of novel pollutants.

Comparison of bacterial diversity in full scale anammox bioreactors operated under different conditions

Bacterial community structure of full-scale anammox bioreactor is still mainly unknown. It has never been analyzed whether different anammox bioreactor configurations might result in the development of different bacterial community structures among these systems. In this work, the bacterial community structure of six full-scale autotrophic nitrogen removal bioreactors located in The Netherlands and China operating under three different technologies and with different influent wastewater characteristics was studied by the means of pyrotag sequencing evaluation of the bacterial assemblage yielded a great diversity in all systems. The most represented phyla were the Bacteroidetes and the Proteobacteria, followed by the Planctomycetes. 14 OTUs were shared by all bioreactors, but none of them belonged to the Brocadiales order. Statistical analysis at OTU level showed that differences in the microbial communities were high, and that the main driver of the bacterial assemblage composition was different for the distinct phyla identified in the six bioreactors, depending on bioreactor technology or influent wastewater characteristics.

Microbiological Analysis of Necrosols Collected from Urban Cemeteries in Poland

Decomposition of organic matter is the primary function in the soil ecosystem, which involves bacteria and fungi. Soil microbial content depends on many factors, and secondary biological and chemical contaminations change and affect environmental feedback. Little work has been done to estimate the microbiological risk for cemetery employees and visitors. The potential risk of infection for people in the cemetery is primarily associated with injury and wound contamination during performing the work. The aim of this study was to analyze the microbiota of cemetery soil obtained from cemeteries and bacterial composition in selected soil layers encountered by gravediggers and cemetery caretakers. The most common bacterial pathogens were Enterococcus spp. (80.6%), Bacillus spp. (77.4%), and E. coli (45.1%). The fungi Penicillium spp. and Aspergillus spp. were isolated from 51% and 6.4% of samples, respectively. Other bacterial species were in the ground cemetery relatively sparse. Sampling depth was not correlated with bacterial growth (p > 0.05), but it was correlated with several differences in microbiota composition (superficial versus deep layer).

Influence of influent wastewater communities on temporal variation of activated sludge communities

Continuously feeding influent wastewater containing diverse bacterial species to a wastewater treatment activated sludge bioreactor may influence the activated sludge bacterial community temporal dynamics. To explore this possibility, this study tracked influent wastewater and activated sludge bacterial communities by pyrosequencing 16S rRNA genes from four full-scale wastewater treatment plants over a 9-month period. The activated sludge communities showed significantly higher richness and evenness than the influent wastewater communities. Furthermore, the two communities were different in composition and temporal dynamics. These results demonstrate that the impact of the influent wastewater communities on the activated sludge communities was weak. Nevertheless, 4.3-9.3% of the operational taxonomic units (OTUs) detected in the activated sludge were shared with the influent wastewater, implying contribution from influent wastewater communities to some extent. However, the relative OTU abundance of the influent wastewater was not maintained in the activated sludge communities (i.e., weak neutral assembly). In addition, the variability of the communities of the shared OTUs was moderately correlated with abiotic factors imposed to the bioreactors. Taken together, temporal dynamics of activated sludge communities appear to be predominantly explained by species sorting processes in response to influent wastewater communities.