Metabolic model for the filamentous 'Candidatus Microthrix parvicella' based on genomic and metagenomic analyses

Seasonal Changes in Bacterial Communities Cause Foaming in a Wastewater Treatment Plant

Bio-foaming is a major problem in solid separation in activated sludge (AS) wastewater treatment systems. Understanding the changes in bacterial communities during sludge foaming is vital for explaining foam formation. Changes in bacterial communities in the foam, corresponding foaming AS, and non-foaming AS in a seasonal foaming wastewater treatment plant (WWTP) in Northern China were investigated by high-throughput pyrosequencing and molecular quantification-based approaches. We found that bacterial communities of the foam and the corresponding foaming AS were similar but markedly different from those of the non-foaming AS. Actinobacteria was the predominant phylum in the foam and the corresponding foaming AS, whereas Proteobacteria was predominant in the non-foaming AS. Similar to the results of most previous studies, our results showed that Candidatus "Microthrix parvicella" was the predominant filamentous bacteria in the foam and the corresponding foaming AS and was significantly enriched in the foam compared to the corresponding foaming AS. Its abundance decreased gradually with a slow disappearance of sludge foaming, indicating that its overgrowth had a direct relationship with sludge foaming. In addition to Candidatus M. parvicella, Tetrasphaera and Trichococcus might play a role in sludge foaming, because they supported the changes in AS microbial ecology for foam formation. The effluent water quality of the surveyed plant remained stable during the period of sludge foaming, but the microbial consortia responsible for nitrogen and phosphorus transformation and removal markedly changed compared to that in the non-foaming AS. This study adds to the previous understanding of bacterial communities causing foaming in WWTPs.

Genomic and in situ investigations of the novel uncultured Chloroflexi associated with 0092 morphotype filamentous bulking in activated sludge

Overgrowth of filamentous bacteria in activated sludge wastewater treatment plants (WWTPs) leads to impaired sludge settleability, a condition known as bulking, which is a common operational problem worldwide. Filaments with the Eikelboom 0092 morphotype are commonly associated with such bulking episodes. Members of the uncultured B45 phylotype, which is embraced within the phylum Chloroflexi, were recently shown to exhibit this morphology. Although these organisms are among the most abundant populations recorded in activated sludge processes, nothing is known about their metabolic characteristics. In this study, a genome sequence, representing the B45 phylotype, was retrieved from a metagenome generated from an activated sludge WWTP. The genome consisted of two chromosomes and one plasmid, which were 4.0, 1.0 and 0.04 Mbps in size, respectively. A metabolic model was constructed for this organism, based on annotation of its genome, showing its ability to generate energy by respiration, utilizing oxygen, nitrite or nitrous oxide as electron acceptors, or by fermentation of sugars. The ability of B45 members to ferment sugars under anaerobic conditions was validated in situ with microautoradiography-fluorescence in situ hybridization. The provisional name of 'Candidatus Promineofilum breve' is proposed for this species. This study represents the first detailed information on an uncultured genus of filamentous organisms from activated sludge.

Genomic, Proteomic, and Metabolite Characterization of Gemfibrozil-Degrading Organism Bacillus sp. GeD10

Gemfibrozil is a widely used hypolipidemic and triglyceride lowering drug. Excess of the drug is excreted and discharged into the environment primarily via wastewater treatment plant effluents. Bacillus sp. GeD10, a gemfibrozil-degrader, was previously isolated from activated sludge. It is the first identified bacterium capable of degrading gemfibrozil. Gemfibrozil degradation by Bacillus sp. GeD10 was here studied through genome sequencing, quantitative proteomics and metabolite analysis. From the bacterial proteome of Bacillus sp. GeD10 1974 proteins were quantified, of which 284 proteins were found to be overabundant by more than 2-fold (FDR corrected p-value ≤0.032, fold change (log2) ≥ 1) in response to gemfibrozil exposure. Metabolomic analysis identified two hydroxylated intermediates as well as a glucuronidated hydroxyl-metabolite of gemfibrozil. Overall, gemfibrozil exposure in Bacillus sp. GeD10 increased the abundance of several enzymes potentially involved in gemfibrozil degradation as well as resulted in the production of several gemfibrozil metabolites. The potential catabolic pathway/modification included ring-hydroxylation preparing the substrate for subsequent ring cleavage by a meta-cleaving enzyme. The identified genes may allow for monitoring of potential gemfibrozil-degrading organisms in situ and increase the understanding of microbial processing of trace level contaminants. This study represents the first omics study on a gemfibrozil-degrading bacterium.

Factors affecting population of filamentous bacteria in wastewater treatment plants with nutrients removal

Filamentous population in activated sludge and key operational parameters of full-scale municipal wastewater treatment plants (WWTPs) with bulking problems representative for Poland were investigated with quantitative fluorescence in situ hybridization. Statistical analyses revealed few relationships between operational parameters and biovolume of filamentous bacteria. Sludge age was not only positively correlated with abundance of Chloroflexi (parametric correlation and principal component analysis (PCA)), but also differentiated Microthrix population (analysis of variance (ANOVA)). Phylum Chloroflexi and pH presented a negative relation during the study (PCA). ANOVA showed that pH of influent and sludge volume index (SVI) differentiated abundance of types 0803 and 1851 of Chloroflexi and candidate division TM7. SVI increased along with higher abundance of Microthrix (positive parametric and non-parametric correlations and positive relation in PCA). Biovolumes of morphotypes 0803 and 1851 of Chloroflexi were differentiated by organic matter in influent, also by nutrients in the case of Chloroflexi type 1851. Chemical and biological oxygen demands (COD and BOD5, respectively) were negatively correlated with Microthrix. COD also differentiated the abundance of Haliscomenobacter hydrossis. Results of the study can be used to prevent WWTPs from excessive proliferation of filamentous bacteria and operational problems caused by them--bulking and foaming of activated sludge.

Evaluating the influence of wastewater composition on the growth of Microthrix parvicella by GCxGC/qMS and real-time PCR

This study underlines the significance of long chain fatty acid (LCFA) content in wastewater influents as an influencing factor promoting the growth of Candidatus 'Microthrix parvicella' (M. parvicella), the most common filamentous bacteria causing foam in activated sludge systems worldwide. Quantification of M. parvicella by real-time polymerase chain reaction (real-time PCR) and analysis of LCFAs by means of two-dimensional gas chromatography coupled with mass spectrometry (GCxGC/qMS), involving solid phase micro-extraction (SPME) to enhance sensitivity, were combined for the first time as a monitoring tool. The results indicate a highly significant correlation between the abundance of M. parvicella and the total LCFA loading (r = 0.96) and linolenic acid C18:3 (r = 0.98) in particular. Additionally, comparison of slope values for the direct correlations of all significant LCFAs found in the analyses showed that the influence of LCFAs on M. parvicella growth increases with an increasing degree of unsaturation of carbon chains. These findings suggest that by removing lipid compounds from the incoming waters, substrate availability would be limited for M. parvicella.

Identification of genes coding for putative wax ester synthase/diacylglycerol acyltransferase enzymes in terrestrial and marine environments

Synthesis of neutral lipids such as triacylglycerols (TAG) and wax esters (WE) is catalyzed in bacteria by wax ester synthase/diacylglycerol acyltransferase enzymes (WS/DGAT). We investigated the diversity of genes encoding this enzyme in contrasting natural environments from Patagonia (Argentina). The content of petroleum hydrocarbons in samples collected from oil-producing areas was measured. PCR-based analysis covered WS/DGAT occurrence in marine sediments and soil. No product was obtained in seawater samples. All clones retrieved from marine sediments affiliated with gammaproteobacterial sequences and within them, most phylotypes formed a unique cluster related to putative WS/DGAT belonging to marine OM60 clade. In contrast, soils samples contained phylotypes only related to actinomycetes. Among them, phylotypes affiliated with representatives largely or recently reported as oleaginous bacteria, as well as with others considered as possible lipid-accumulating bacteria based on the analysis of their annotated genomes. Our study shows for the first time that the environment could contain a higher variety of ws/dgat than that reported from bacterial isolates. The results of this study highlight the relevance of the environment in a natural process such as the synthesis and accumulation of neutral lipids. Particularly, both marine sediments and soil may serve as a useful source for novel WS/DGAT with biotechnological interest.

In situ phenotypic heterogeneity among single cells of the filamentous bacterium Candidatus Microthrix parvicella

Microorganisms in biological wastewater treatment plants require adaptive strategies to deal with rapidly fluctuating environmental conditions. At the population level, the filamentous bacterium Candidatus Microthrix parvicella (Ca. M. parvicella) has been found to fine-tune its gene expression for optimized substrate assimilation. Here we investigated in situ substrate assimilation by single cells of Ca. M. parvicella using nano-scale secondary-ion mass spectrometry (nanoSIMS). NanoSIMS imaging highlighted phenotypic heterogeneity among Ca. M. parvicella cells of the same filament, whereby ¹³C-oleic acid and ¹³C-glycerol-3-phosphate assimilation occurred in ≈21–55% of cells, despite non-assimilating cells being intact and alive. In response to alternating aerobic–anoxic regimes, ¹³C-oleic acid assimilation occurred among subpopulations of Ca. M. parvicella cells (≈3–28% of cells). Furthermore, Ca. M. parvicella cells exhibited two temperature optima for ¹³C-oleic acid assimilation and associated growth rates. These results suggest that phenotypic heterogeneity among Ca. M. parvicella cells allows the population to adapt rapidly to fluctuating environmental conditions facilitating its widespread occurrence in biological wastewater treatment plants.

Comparative integrated omics: identification of key functionalities in microbial community-wide metabolic networks

BACKGROUND

Mixed microbial communities underpin important biotechnological processes such as biological wastewater treatment (BWWT). A detailed knowledge of community structure and function relationships is essential for ultimately driving these systems towards desired outcomes, e.g., the enrichment in organisms capable of accumulating valuable resources during BWWT.

METHODS

A comparative integrated omic analysis including metagenomics, metatranscriptomics and metaproteomics was carried out to elucidate functional differences between seasonally distinct oleaginous mixed microbial communities (OMMCs) sampled from an anoxic BWWT tank. A computational framework for the reconstruction of community-wide metabolic networks from multi-omic data was developed. These provide an overview of the functional capabilities by incorporating gene copy, transcript and protein abundances. To identify functional genes, which have a disproportionately important role in community function, we define a high relative gene expression and a high betweenness centrality relative to node degree as gene-centric and network topological features, respectively.

RESULTS

Genes exhibiting high expression relative to gene copy abundance include genes involved in glycerolipid metabolism, particularly triacylglycerol lipase, encoded by known lipid accumulating populations, e.g., CandidatusMicrothrix parvicella. Genes with a high relative gene expression and topologically important positions in the network include genes involved in nitrogen metabolism and fatty acid biosynthesis, encoded by Nitrosomonas spp. and Rhodococcus spp. Such genes may be regarded as 'keystone genes' as they are likely to be encoded by keystone species.

CONCLUSION

The linking of key functionalities to community members through integrated omics opens up exciting possibilities for devising prediction and control strategies for microbial communities in the future.

Rare taxa have potential to make metabolic contributions in enhanced biological phosphorus removal ecosystems

Enhanced biological phosphorus removal (EBPR) relies on diverse but specialized microbial communities to mediate the cycling and ultimate removal of phosphorus from municipal wastewaters. However, little is known about microbial activity and dynamics in relation to process fluctuations in EBPR ecosystems. Here, we monitored temporal changes in microbial community structure and potential activity across each bioreactor zone in a pilot-scale EBPR treatment plant by examining the ratio of small subunit ribosomal RNA (SSU rRNA) to SSU rRNA gene (rDNA) over a 120 day study period. Although the majority of operational taxonomic units (OTUs) in the EBPR ecosystem were rare, many maintained high potential activities based on SSU rRNA : rDNA ratios, suggesting that rare OTUs contribute substantially to protein synthesis potential in EBPR ecosystems. Few significant differences in OTU abundance and activity were observed between bioreactor redox zones, although differences in temporal activity were observed among phylogenetically cohesive OTUs. Moreover, observed temporal activity patterns could not be explained by measured process parameters, suggesting that other ecological drivers, such as grazing or viral lysis, modulated community interactions. Taken together, these results point towards complex interactions selected for within the EBPR ecosystem and highlight a previously unrecognized functional potential among low abundance microorganisms in engineered ecosystems.

Integrated production of polyhydroxyalkanoates (PHAs) with municipal wastewater and sludge treatment at pilot scale

A pilot-scale process was operated over 22 months at the Brussels North Wastewater Treatment Plant (WWTP) in order to evaluate polyhydroxyalkanoate (PHA) production integration with services of municipal wastewater and sludge management. Activated sludge was produced with PHA accumulation potential (PAP) by applying feast-famine selection while treating the readily biodegradable COD from influent wastewater (average removals of 70% COD, 60% CODsol, 24% nitrogen, and 46% phosphorus). The biomass PAP was evaluated to be in excess of 0.4gPHA/gVSS. Batch fermentation of full-scale WWTP sludge at selected temperatures (35, 42 and 55 °C) produced centrate (6-9.4 gCODVFA/L) of consistent VFA composition, with optimal fermentation performance at 42 °C. Centrate was used to accumulate PHA up to 0.39 gPHA/gVSS. The centrate nutrients are a challenge to the accumulation process but producing a biomass with 0.5 gPHA/gVSS is considered to be realistically achievable within the typically available carbon flows at municipal waste management facilities.

Community-integrated omics links dominance of a microbial generalist to fine-tuned resource usage

Microbial communities are complex and dynamic systems that are primarily structured according to their members’ ecological niches. To investigate how niche breadth (generalist versus specialist lifestyle strategies) relates to ecological success, we develop and apply an integrative workflow for the multi-omic analysis of oleaginous mixed microbial communities from a biological wastewater treatment plant. Time- and space-resolved coupled metabolomic and taxonomic analyses demonstrate that the community-wide lipid accumulation phenotype is associated with the dominance of the generalist bacterium Candidatus Microthrix spp. By integrating population-level genomic reconstructions (reflecting fundamental niches) with transcriptomic and proteomic data (realised niches), we identify finely tuned gene expression governing resource usage by Candidatus Microthrix parvicella over time. Moreover, our results indicate that the fluctuating environmental conditions constrain the accumulation of genetic variation in Candidatus Microthrix parvicella likely due to fitness trade-offs. Based on our observations, niche breadth has to be considered as an important factor for understanding the evolutionary processes governing (microbial) population sizes and structures in situ.

Within microbial communities, microorganisms adopt different lifestyle strategies to use the available resources. Here, the authors use an integrated ‘multi-omic’ approach to study niche breadth (generalist versus specialist lifestyles) in oleaginous microbial assemblages from an anoxic wastewater treatment tank.

A hundred years of activated sludge: time for a rethink

Biological wastewater treatment plants (BWWTPs) based on the activated sludge (AS) process have dramatically improved worldwide water sanitation despite increased urbanization and industrialization. However, current AS-based operations are considered economically and environmentally unsustainable. In this Perspective, we discuss our current understanding of microbial populations and their metabolic transformations in AS-based BWWTPs in view of developing more sustainable processes in the future. In particular, much has been learned over the course of the past 25 years about specialized microorganisms, which could be more comprehensively leveraged to recover energy and/or nutrients from wastewater streams. To achieve this, we propose a bottom-up design approach, focused around the concept of a "wastewater biorefinery column", which would rely on the engineering of distinct ecological niches into a BWWTP in order to guarantee the targeted enrichment of specific organismal groups which in turn will allow the harvest of high-value resources from wastewater. This concept could be seen as a possible grand challenge to microbial ecologists and engineers alike at the centenary of the discovery of the AS process.

'Candidatus Competibacter'-lineage genomes retrieved from metagenomes reveal functional metabolic diversity

The glycogen-accumulating organism (GAO) 'Candidatus Competibacter' (Competibacter) uses aerobically stored glycogen to enable anaerobic carbon uptake, which is subsequently stored as polyhydroxyalkanoates (PHAs). This biphasic metabolism is key for the Competibacter to survive under the cyclic anaerobic-'feast': aerobic-'famine' regime of enhanced biological phosphorus removal (EBPR) wastewater treatment systems. As they do not contribute to phosphorus (P) removal, but compete for resources with the polyphosphate-accumulating organisms (PAO), thought responsible for P removal, their proliferation theoretically reduces the EBPR capacity. In this study, two complete genomes from Competibacter were obtained from laboratory-scale enrichment reactors through metagenomics. Phylogenetic analysis identified the two genomes, 'Candidatus Competibacter denitrificans' and 'Candidatus Contendobacter odensis', as being affiliated with Competibacter-lineage subgroups 1 and 5, respectively. Both have genes for glycogen and PHA cycling and for the metabolism of volatile fatty acids. Marked differences were found in their potential for the Embden-Meyerhof-Parnas and Entner-Doudoroff glycolytic pathways, as well as for denitrification, nitrogen fixation, fermentation, trehalose synthesis and utilisation of glucose and lactate. Genetic comparison of P metabolism pathways with sequenced PAOs revealed the absence of the Pit phosphate transporter in the Competibacter-lineage genomes--identifying a key metabolic difference with the PAO physiology. These genomes are the first from any GAO organism and provide new insights into the complex interaction and niche competition between PAOs and GAOs in EBPR systems.

Modular community structure suggests metabolic plasticity during the transition to polar night in ice-covered Antarctic lakes

High-latitude environments, such as the Antarctic McMurdo Dry Valley lakes, are subject to seasonally segregated light-dark cycles, which have important consequences for microbial diversity and function on an annual basis. Owing largely to the logistical difficulties of sampling polar environments during the darkness of winter, little is known about planktonic microbial community responses to the cessation of photosynthetic primary production during the austral sunset, which lingers from approximately February to April. Here, we hypothesized that changes in bacterial, archaeal and eukaryotic community structure, particularly shifts in favor of chemolithotrophs and mixotrophs, would manifest during the transition to polar night. Our work represents the first concurrent molecular characterization, using 454 pyrosequencing of hypervariable regions of the small-subunit ribosomal RNA gene, of bacterial, archaeal and eukaryotic communities in permanently ice-covered lakes Fryxell and Bonney, before and during the polar night transition. We found vertically stratified populations that varied at the community and/or operational taxonomic unit-level between lakes and seasons. Network analysis based on operational taxonomic unit level interactions revealed nonrandomly structured microbial communities organized into modules (groups of taxa) containing key metabolic potential capacities, including photoheterotrophy, mixotrophy and chemolithotrophy, which are likely to be differentially favored during the transition to polar night.