Metatranscriptomic insights on gene expression and regulatory controls in Candidatus Accumulibacter phosphatis

Candidatus Thiothrix phosphatis SCUT-1: A novel polyphosphate-accumulating organism abundant in the enhanced biological phosphorus removal system

A novel coccus Thiothrix-related polyphosphate-accumulating organism (PAO) was enriched in an acetate-fed enhanced biological phosphorus removal (EBPR) system. High EBPR performance was achieved for an extended period (>100 days). A high-quality draft genome (completeness 97.2 %, contamination 3.26 %) was retrieved, representing a novel Thiothrix species (with similarity<93.2 % to known Thiothrix species), and was denoted as 'Candidatus Thiothrix phosphatis SCUT-1'. Its acetate uptake rate (6.20 mmol C/g VSS/h) surpassed most Ca. Accumulibacter and known glycogen-accumulating organisms (GAOs), conferring their predominance in the acetate-fed system. Metatranscriptomic analysis suggested that Ca. Thiothrix phosphatis SCUT-1 employed both low- and high-affinity pathways for acetate activation, and both the conventional (PhaABC) pathway and the fatty acid β-oxidation pathway for PHA synthesis; additionally, a much more efficient FAD-dependent malate: quinone oxidoreductase (MQO) were encoded and employed than the traditional malate dehydrogenase (MDH) to oxidize malate to oxaloacetate in the TCA and glyoxylate cycle, collectively contributing to a higher acetate utilization and processing rate of this microorganism. Batch tests further demonstrated the versatile ability of this PAO in using VFA (acetate, propionate, and butyrate), lactate, amino acids (aspartate and glutamate), and glucose as carbon sources for EBPR, showing a partially overlapped but unique ecological niche of this microorganism comparing to Ca. Accumulibacter and known GAOs. A metabolic model was built for Ca. Thiothrix phosphatis SCUT-1 using the above-mentioned carbon sources for EBPR. Overall, this study represents the first comprehensive characterization of the physiology and metabolic characteristics of representative coccus Thiothrix-related PAOs, which are expected to provide new insights into PAO microbiology in EBPR systems.

Integrated genomics provides insights into the evolution of the polyphosphate accumulation trait of Ca. Accumulibacter

Candidatus Accumulibacter, a prominent polyphosphate-accumulating organism (PAO) in wastewater treatment, plays a crucial role in enhanced biological phosphorus removal (EBPR). The genetic underpinnings of its polyphosphate accumulation capabilities, however, remain largely unknown. Here, we conducted a comprehensive genomic analysis of Ca. Accumulibacter-PAOs and their relatives within the Rhodocyclaceae family, identifying 124 core genes acquired via horizontal gene transfer (HGT) at its least common ancestor. Metatranscriptomic analysis of an enrichment culture of Ca. Accumulibacter revealed active transcription of 44 of these genes during an EBPR cycle, notably including the polyphosphate kinase 2 (PPK2) gene instead of the commonly recognized polyphosphate kinase 1 (PPK1) gene. Intriguingly, the phosphate regulon (Pho) genes showed minimal transcriptions, pointing to a distinctive fact of Pho dysregulation, where PhoU, the phosphate signaling complex protein, was not regulating the high-affinity phosphate transport (Pst) system, resulting in continuous phosphate uptake. To prevent phosphate toxicity, Ca. Accumulibacter utilized the laterally acquired PPK2 to condense phosphate into polyphosphate, resulting in the polyphosphate-accumulating feature. This study provides novel insights into the evolutionary emergence of the polyphosphate-accumulating trait in Ca. Accumulibacter, offering potential advancements in understanding the PAO phenotype in the EBPR process.

Strategies to attenuate ciprofloxacin inhibition on enhanced biological phosphorus removal from wastewater and its recoverability

The inhibiting effects of ciprofloxacin (CIP) on enhanced biological phosphorus removal (EBPR) were investigated with no change in reactor operation and with increased aeration rate and sludge retention time (SRT) to explore inhibition-alleviating solutions. Additionally, performance recoverability was evaluated. The results showed that the phosphorus removal efficiency in the presence of 0.002-0.092 mg/L CIP for 7 days was only 12.5%. Increasing the aeration rate relieved inhibition (33.5% phosphorus removal efficiency on Day 7), and increasing SRT slowed EBPR performance deterioration. The EBPR performance recovered from CIP inhibition and increases in the aeration rate and SRT resulted in different recovery phenomena. The maximum PO43--P release rate continued to decrease in the first 2 days of the recovery stage and then gradually increased. However, the maximum PO43--P uptake rate immediately increased at different rates among reactors, which might be attributed to variations in the microbial community structure, decreased poly-P content, and enhanced abundances of ABC transporters and quorum sensing. It was found that some microorganisms associated with phosphorus removal were more tolerant to CIP than glycogen accumulating organisms. Moreover, the increased relative abundance of the qepA gene indicated that the microorganisms in the EBPR system had strong antibiotic resistance capacity. The bacterial community structure was significantly affected by CIP and could not recover to the initial structure. The results help to provide technical support for the operation of the EBPR process in the presence of CIP and to increase the understanding of system recoverability.

From metagenomes to metabolism: Systematically assessing the metabolic flux feasibilities for "Candidatus Accumulibacter" species during anaerobic substrate uptake

With the rapid growing availability of metagenome assembled genomes (MAGs) and associated metabolic models, the identification of metabolic potential in individual community members has become possible. However, the field still lacks an unbiassed systematic evaluation of the generated metagenomic information to uncover not only metabolic potential, but also feasibilities of these models under specific environmental conditions. In this study, we present a systematic analysis of the metabolic potential in species of "Candidatus Accumulibacter", a group of polyphosphate-accumulating organisms (PAOs). We constructed a metabolic model of the central carbon metabolism and compared the metabolic potential among available MAGs for "Ca. Accumulibacter" species. By combining Elementary Flux Modes Analysis (EFMA) with max-min driving force (MDF) optimization, we obtained all possible flux distributions of the metabolic network and calculated their individual thermodynamic feasibility. Our findings reveal significant variations in the metabolic potential among "Ca. Accumulibacter" MAGs, particularly in the presence of anaplerotic reactions. EFMA revealed 700 unique flux distributions in the complete metabolic model that enable the anaerobic uptake of acetate and its conversion into polyhydroxyalkanoates (PHAs), a well-known phenotype of "Ca. Accumulibacter". However, thermodynamic constraints narrowed down this solution space to 146 models that were stoichiometrically and thermodynamically feasible (MDF > 0 kJ/mol), of which only 8 were strongly feasible (MDF > 7 kJ/mol). Notably, several novel flux distributions for the metabolic model were identified, suggesting putative, yet unreported, functions within the PAO communities. Overall, this work provides valuable insights into the metabolic variability among "Ca. Accumulibacter" species and redefines the anaerobic metabolic potential in the context of phosphate removal. More generally, the integrated workflow presented in this paper can be applied to any metabolic model obtained from a MAG generated from microbial communities to objectively narrow the expected phenotypes from community members.

Two new clades recovered at high temperatures provide novel phylogenetic and genomic insights into Candidatus Accumulibacter

Candidatus Accumulibacter, a key genus of polyphosphate-accumulating organisms, plays key roles in lab- and full-scale enhanced biological phosphorus removal (EBPR) systems. A total of 10 high-quality Ca. Accumulibacter genomes were recovered from EBPR systems operated at high temperatures, providing significantly updated phylogenetic and genomic insights into the Ca. Accumulibacter lineage. Among these genomes, clade IIF members SCELSE-3, SCELSE-4, and SCELSE-6 represent the to-date known genomes encoding a complete denitrification pathway, suggesting that Ca. Accumulibacter alone could achieve complete denitrification. Clade IIC members SSA1, SCUT-1, SCELCE-2, and SCELSE-8 lack the entire set of denitrifying genes, representing to-date known non-denitrifying Ca. Accumulibacter. A pan-genomic analysis with other Ca. Accumulibacter members suggested that all Ca. Accumulibacter likely has the potential to use dicarboxylic amino acids. Ca. Accumulibacter aalborgensis AALB and Ca. Accumulibacter affinis BAT3C720 seemed to be the only two members capable of using glucose for EBPR. A heat shock protein Hsp20 encoding gene was found exclusively in genomes recovered at high temperatures, which was absent in clades IA, IC, IG, IIA, IIB, IID, IIG, and II-I members. High transcription of this gene in clade IIC members SCUT-2 and SCUT-3 suggested its role in surviving high temperatures for Ca. Accumulibacter. Ambiguous clade identity was observed for newly recovered genomes (SCELSE-9 and SCELSE-10). Five machine learning models were developed using orthogroups as input features. Prediction results suggested that they belong to a new clade (IIK). The phylogeny of Ca. Accumulibacter was re-evaluated based on the laterally derived polyphosphokinase 2 gene, showing improved resolution in differentiating different clades.

Microbial community dynamics and metagenomics reveal the potential role of unconventional functional microorganisms in nitrogen and phosphorus removal biofilm system

The conventional functional microorganisms for nitrogen and phosphorus removal, such as Nitrosomonas, Nitrobacter, Nitrospira and Candidatus Accumulibacter, were hotspots in past research. However, the role of diverse unconventional functional microorganisms was neglected. In this study, a biofilm system was developed to explore the potential role of unconventional functional microorganisms in nutrients removal. According to the results of microbial community dynamics and metagenomics, complete ammonia oxidizing (comammox) bacteria was 20 times more abundant than ammonia-oxidizing bacteria (AOB) at day 121 and its abundance of amoA gene was almost the same as AOB. Although Nitrospira dominated the nitrite-oxidizing bacteria (NOB), diverse unconventional nxrB-containing microorganisms, particularly Chloroflexi, also significantly contributed to the nitrite oxidation. Binning analysis showed that Myxococcota-affiliated Haliangium had the necessary genes owns by phosphorus-accumulating organisms (PAO) and was likely to be the primary PAO since its abundance (6.38 %) was much higher than other conventional PAO (0.70 %). Comparing metagenome-assembled genomes of comammox bacteria with AOB and ammonia-oxidizing archaea (AOA), it possessed potential metabolic versatility in hydrogen and phosphorus, which may be the primary reason for the positive effect of the alternating anaerobic and aerobic conditions on the enrichment of comammox bacteria. Collectively, our findings broaden the understanding on the microbial mechanism of nitrogen and phosphorus removal in biofilm system.

New insight into the granule formation in the reactor for enhanced biological phosphorus removal

While granulated activated sludge exhibits high productivity, the processes of granule formation are incompletely studied. The processes of granule formation and succession of communities were investigated in a laboratory sequencing batch reactor (SBR) under conditions for enhanced biological phosphorus removal (EBPR) using microbiological and molecular techniques. Active consumption of acetate, primarily by the phosphate-accumulating organisms (PAO), commenced at day 150 of cultivation. This was indicated by the high ratio of molar P-released/acetate uptake (0.73-0.77 P-mol/C-mol), characteristic of PAO. During this period, two types of granule-like aggregates formed spontaneously out of the activated sludge flocs. The aggregates differed in morphology and microbial taxonomic composition. While both aggregate types contained phosphorus-enriched bacterial cells, PAO prevailed in those of morphotype I, and glycogen-accumulating organisms (GAOs) were predominant in the aggregates of morphotype II. After 250 days, the elimination of the morphotype II aggregates from the reactor was observed. The subsequent selection of the community was associated with the development of the morphotype I aggregates, in which the relative abundance of PAO increased significantly, resulting in higher efficiency of phosphorus removal. Metagenomic analysis revealed a predominance of the organisms closely related to Candidatus Accumulibacter IС and IIС and of Ca. Accumulibacter IIB among the PAO. Based on the content of the genes of the key metabolic pathways, the genomes of potential PAO belonging to the genera Amaricoccus, Azonexus, Thauera, Zoogloea, Pinisolibacter, and Siculibacillus were selected. The patterns of physicochemical processes and the microbiome structure associated with granule formation and succession of the microbial communities were revealed.

Candidatus Accumulibacter use fermentation products for enhanced biological phosphorus removal

Previous research suggested that two major groups of polyphosphate-accumulating organisms (PAOs), i.e., Ca. Accumulibacter and Tetrasphaera, play cooperative roles in enhanced biological phosphorus removal (EBPR). The fermentation of complex organic compounds by Tetrasphaera provides carbon sources for Ca. Accumulibacter. However, the viability of the fermentation products (e.g., lactate, succinate, alanine) as carbon sources for Ca. Accumulibacter and their potential effects on the metabolism of Ca. Accumulibacter were largely unknown. This work for the first time investigated the capability and metabolic details of Ca. Accumulibacter cognatus clade IIC strain SCUT-2 (enriched in a lab-scale reactor with a relative abundance of 42.8%) in using these fermentation products for EBPR. The enrichment culture was able to assimilate lactate and succinate with the anaerobic P release to carbon uptake ratios of 0.28 and 0.36 P mol/C mol, respectively. In the co-presence of acetate, the uptake of lactate was strongly inhibited, since two substrates shared the same transporter as suggested by the carbon uptake bioenergetic analysis. When acetate and succinate were fed at the same time, Ca. Accumulibacter assimilated two carbon sources simultaneously. Proton motive force (PMF) was the key driving force (up to 90%) for the uptake of lactate and succinate by Ca. Accumulibacter. Apart from the efflux of proton in symport with phosphate via the inorganic phosphate transport system, translocation of proton via the activity of fumarate reductase contributed to the generation of PMF, which agreed with the fact that PHV was a major component of PHA when lactate and succinate were used as carbon sources, involving the succinate-propionate pathway. Metabolic models for the usage of lactate and succinate by Ca. Accumulibacter for EBPR were built based on the combined physiological, biochemical, metagenomic, and metatranscriptomic analyses. Alanine was shown as an invalid carbon source for Ca. Accumulibacter. Instead, it significantly and adversely affected Ca. Accumulibacter-mediated EBPR. Phosphate release was observed without alanine uptake. Significant inhibitions on the aerobic phosphate uptake was also evident. Overall, this study suggested that there might not be a simply synergic relationship between Ca. Accumulibacter and Tetrasphaera. Their interactions would largely be determined by the kind of fermentation products released by the latter.

Frequency of change determines effectiveness of microbial response strategies

Nature challenges microbes with change at different frequencies and demands an effective response for survival. Here, we used controlled laboratory experiments to investigate the effectiveness of different response strategies, such as post-translational modification, transcriptional regulation, and specialized versus adaptable metabolisms. For this, we inoculated replicated chemostats with an enrichment culture obtained from sulfidic stream microbiomes 16 weeks prior. The chemostats were submitted to alternatingly oxic and anoxic conditions at three frequencies, with periods of 1, 4 and 16 days. The microbial response was recorded with 16S rRNA gene amplicon sequencing, shotgun metagenomics, transcriptomics and proteomics. Metagenomics resolved provisional genomes of all abundant bacterial populations, mainly affiliated with Proteobacteria and Bacteroidetes. Almost all these populations maintained a steady growth rate under both redox conditions at all three frequencies of change. Our results supported three conclusions: (1) Oscillating oxic/anoxic conditions selected for generalistic species, rather than species specializing in only a single condition. (2) A high frequency of change selected for strong codon usage bias. (3) Alignment of transcriptomes and proteomes required multiple generations and was dependent on a low frequency of change.

Effect of different types of volatile fatty acids on the performance and bacterial population in a batch reactor performing biological nutrient removal

Different ratios of four volatile fatty acids (VFAs) were used as the primary feed to a laboratory scale biological nutrient reactor during four operational stages. The reactor performed efficiently over 500 days of operation with over 90% dissolved phosphorus and over 98% ammonium-nitrogen (NH4+-N) removal. Through in the first experimental phase, acetate and propionate were present in a significant proportion as carbon sources, the relative abundance of Candidatus Accumulibacter, a potential polyphosphate accumulating organism, increased from 10% to 57% and the Defluviicoccus genus, a known glycogen accumulating organism (GAO), decreased from 41% to 5%. Further tests indicated the presence of denitrifying phosphorus accumulating organisms (DPAO) belonging to Clade IIC, that could use nitrite as the electron acceptor during P-uptake. In general, VFAs favored the increase of the genus Defluviicoccus and Candidatus Accumulibacter. High relative abundance of Defluviicoccus did not affect the stability and the performance of the BNR process.

Microbiome assembly mechanism and functional potential in enhanced biological phosphorus removal system enriched with Tetrasphaera-related polyphosphate accumulating organisms

Tetrasphaera-related polyphosphate accumulating organisms (PAOs) are the key functional guilds for enhanced biological phosphorus removal (EBPR) systems. Their biomass enrichment can be enhanced by the nitrification inhibitor allylthiourea (ATU). However, the underlying assembly mechanism and the functional potential of the EBPR microbiome regulated by ATU are unclear. This study investigates the effect of ATU on microbiome assembly and functional potential by closely following the microbiota dynamics in an EBPR system enriched with Tetrasphaera-related PAOs for 288-days before, during and after ATU addition. The results showed that ATU addition increased microbiota structural similarity and compositional convergence, and enhanced determinism in the assembly of EBPR microbiome. During exposure to ATU, Tetrasphaera-related PAOs were governed by homogeneous selection and the dominant species revealed by 16S rRNA gene-based phylogenetic analysis shifted from clade III to clade I. Meanwhile, ATU supply promoted significant enrichment of functional genes involved in phosphate transport (pit) and polyphosphate synthesis and degradation (ppk1 and ppk2), whereas both Nitrosomonas and ammonia monooxygenase-encoding genes (amoA/B/C) assignable to this group of nitrifying bacteria decreased. Moreover, ATU addition relieved the significant abundance correlation between filamentous bacteria Ca. Promineofilum and denitrifying Brevundimonas (FDR-adjusted P < 0.01), damaging their potential synergic or cooperative interactions, thus weakening their competitiveness against Tetrasphaera-related PAOs. Notably, ATU withdrawn created opportunistic conditions for the unexpected explosive growth and predominance of Thiothrix filaments, leading to a serious bulking event. Our study provides new insights into the microbial ecology of Tetrasphaera-related PAOs in EBPR system, which could guide the establishment of an efficient microbiota for EBPR.

Expanding the Diversity of Accumulibacter with a Novel Type and Deciphering the Transcriptional and Morphological Features among Co-Occurring Strains

"Candidatus Accumulibacter" is the major polyphosphate-accumulating organism (PAO) in global wastewater treatment systems, and its phylogenetic and functional diversity have expanded in recent years. In addition to the widely recognized type I and II sublineages, we discovered a novel type enriched in laboratory bioreactors. Core gene and machine learning-based gene feature profiling supported the assertion that type III "Ca. Accumulibacter" is a potential PAO with the unique function of using dimethyl sulfoxide as an electron acceptor. Based on the correlation between ppk1 and genome similarity, the species-level richness of Accumulibacter was estimated to be over 100, suggesting that the currently recognized species are only the tip of the iceberg. Meanwhile, the interstrain transcriptional and morphological features of multiple "Ca. Accumulibacter" strains co-occurring in a bioreactor were investigated. Metatranscriptomics of seven co-occurring strains indicated that the expression level and interphasic dynamics of PAO phenotype-related genes had minimal correlation with their phylogeny. In particular, the expression of denitrifying and polyphosphate (poly-P) metabolism genes exhibited higher interstrain and interphasic divergence than expression of glycogen and polyhydroxyalkanoate metabolic genes. A strategy of cloning rRNA genes from different strains based on similar genomic synteny was successfully applied to differentiate their morphology via fluorescence in situ hybridization. Our study further expands the phylogenetic and functional diversity of "Ca. Accumulibacter" and proposes that deciphering the function and capability of certain "Ca. Accumulibacter" should be tailored to the environment and population in question. IMPORTANCE In the last 2 decades, "Ca. Accumulibacter" has garnered significant attention as the core functional but uncultured taxon for enhanced biological phosphorus removal due to its phylogenetic and functional diversity and intragenus niche differentiation. Since 2002, it has been widely known that this genus has two sublineages (type I and II). However, in this study, a metagenomic approach led to the discovery of a novel type (type III) with proposed novel functional features. By comparing the average nucleotide identity of "Ca. Accumulibacter" genomes and the similarity of ppk1, a phylogenetic biomarker largely deposited in databases, the global species-level richness of "Ca. Accumulibacter" was estimated for the first time to be over 100. Furthermore, we observed the co-occurrence of multiple "Ca. Accumulibacter" strains in a single bioreactor and found the simultaneous transcriptional divergence of these strains intriguing with regard to their niche differentiation within a single community. Our results indicated a decoupling feature between transcriptional pattern and phylogeny for co-occurring strains.

Impact of the anaerobic feeding mode on substrate distribution in aerobic granular sludge

There is a growing interest to implement aerobic granular sludge (AGS) in existing conventional activated sludge (CAS) systems with a continuous flow-through configuration. The mode of anaerobic contact of raw sewage with the sludge is an important aspect in the adaptation of CAS systems to accommodate AGS. It remains unclear how the distribution of substrate over the sludge by a conventional anaerobic selector compares to the distribution via bottom-feeding applied in sequencing batch reactors (SBRs). This study investigated the effect of the anaerobic contact mode on the substrate (and storage) distribution by operating two lab-scale SBRs; one with the traditional bottom-feeding through a settled sludge bed similar to full-scale AGS systems, and one where the synthetic wastewater was fed as a pulse at the start of the anaerobic phase while the reactor was mixed through sparging of nitrogen gas (mimicking a plug-flow anaerobic selector in continuous flow-through systems). The distribution of the substrate over the sludge particle population was quantified via PHA analysis, combined with the obtained granule size distribution. Bottom-feeding was found to primarily direct substrate towards the large granular size classes (i.e. large volume and close to the bottom), while completely mixed pulse-feeding gives a more equal distribution of substrate over all granule sizes (i.e. surface area dependant). The anaerobic contact mode directly controls the substrate distribution over the different granule sizes, irrespective of the solids retention time of a granule as an entity. Preferential feeding of the larger granules will enhance and stabilise the granulation compared to pulse-feeding, certainly under less advantageous conditions imposed by real sewage.

Putative metabolism of Ca. Accumulibacter via the utilization of glucose

Ca. Accumulibacter was the predominant microorganism (relative FISH bio-abundance of 67 ± 5%) in a lab-scale sequential batch reactor that accomplished enhanced biological phosphorus removal (EBPR) while using glucose and acetate as the carbon sources (1:1 COD-based ratio). Both organic compounds were completely anaerobically consumed. The reactor's performance in terms of P/C ratio, phosphorous release and uptake, and overall kinetic and stoichiometric parameters were on the high end of the reported spectrum for EBPR systems (100:9.3 net mg phosphate removal per mg COD consumed when using glucose and acetate in a 1:1 ratio). The batch tests showed that, to the best of our knowledge, this is the first time a reactor enriched with Ca. Accumulibacter can putatively utilize glucose as the sole carbon source to biologically remove phosphate (COD:P (mg/mg) removal ratio of 100:6.3 when using only glucose). Thus, this research proposes that Ca. Accumulibacter directly anaerobically stored the fed glucose primarily as glycogen by utilizing the ATP provided via the hydrolysis of poly-P and secondarily as PHA by balancing its ATP utilization (glycogen generation) and formation (PHA storage). Alternative hypotheses are also discussed. The reported findings could challenge the conventional theories of glucose assimilation by Ca. Accumulibacter, and can be of significance for the biological removal of phosphorus from wastewaters with high contents of fermentable compounds or low VFAs.

Predicting the impact of temperature on metabolic fluxes using resource allocation modelling: Application to polyphosphate accumulating organisms

The understanding of microbial communities and the biological regulation of its members is crucial for implementation of novel technologies using microbial ecology. One poorly understood metabolic principle of microbial communities is resource allocation and biosynthesis. Resource allocation theory in polyphosphate accumulating organisms (PAOs) is limited as a result of their slow imposed growth rate (typical sludge retention times of at least 4 days) and limitations to quantify changes in biomass components over a 6 hours cycle (less than 10% of their growth). As a result, there is no direct evidence supporting that biosynthesis is an exclusive aerobic process in PAOs that alternate continuously between anaerobic and aerobic phases. Here, we apply resource allocation metabolic flux analysis to study the optimal phenotype of PAOs over a temperature range of 4 °C to 20 °C. The model applied in this research allowed to identify optimal metabolic strategies in a core metabolic model with limited constraints based on biological principles. The addition of a constraint limiting biomass synthesis to be an exclusive aerobic process changed the metabolic behaviour and improved the predictability of the model over the studied temperature range by closing the gap between prediction and experimental findings. The results validate the assumption of limited anaerobic biosynthesis in PAOs, specifically "Candidatus Accumulibacter" related species. Interestingly, the predicted growth yield was lower, suggesting that there are mechanistic barriers for anaerobic growth not yet understood nor reflected in the current models of PAOs. Moreover, we identified strategies of resource allocation applied by PAOs at different temperatures as a result of the decreased catalytic efficiencies of their biochemical reactions. Understanding resource allocation is paramount in the study of PAOs and their currently unknown complex metabolic regulation, and metabolic modelling based on biological first principles provides a useful tool to develop a mechanistic understanding.

Advanced simultaneous nitrogen and phosphorus removal for non-sterile wastewater through a novel coupled yeast-sludge system: Performance, microbial interaction, and mechanism

A novel coupled yeast-sludge system (CYSS) was constructed by the yeast Candida sp. PNY integrated with activated sludge to treat non-sterile mainstream wastewater. After 240-day cultivation, compared with single activated sludge, simultaneous removal efficiency of total organic carbon (TOC), nitrogen and phosphorus increased by 19.5% (176.34 mg TOC g^-1 d^-1), 21.3% (11.25 mg TN g^-1 d^-1) and 15.0% (6.95 mg TP g^-1 d^-1), respectively, while the amount of sludge reduced by 50%. Amplicon sequencing analysis showed that the abundance of Nitrosomonas, Nitrospira, Zoogloea, Dechloromonas, and Candidatus Accumulibacter significantly decreased to 0% on Day 200. Abundance of nirS and nirK for denitrification significantly decreased in CYSS by quantitative PCR (qPCR), and the copies of nirS and nirK were 3.37-fold and 1.71-fold decrease from Day 0 to Day 240, respectively. The results of Fluorescence in situ hybridization and co-occurrence network showed that Candida sp. PNY predominated its distribution in CYSS, and strongly connected with environmental variables based on network analysis. Furthermore, this study reconstructed the carbon, nitrogen and phosphorus metabolic pathways of the CYSS based on metagenomics.

TbasCO: trait-based comparative 'omics identifies ecosystem-level and niche-differentiating adaptations of an engineered microbiome

A grand challenge in microbial ecology is disentangling the traits of individual populations within complex communities. Various cultivation-independent approaches have been used to infer traits based on the presence of marker genes. However, marker genes are not linked to traits with complete fidelity, nor do they capture important attributes, such as the timing of gene expression or coordination among traits. To address this, we present an approach for assessing the trait landscape of microbial communities by statistically defining a trait attribute as a shared transcriptional pattern across multiple organisms. Leveraging the KEGG pathway database as a trait library and the Enhanced Biological Phosphorus Removal (EBPR) model microbial ecosystem, we demonstrate that a majority (65%) of traits present in 10 or more genomes have niche-differentiating expression attributes. For example, while many genomes containing high-affinity phosphorus transporter pstABCS display a canonical attribute (e.g. up-regulation under phosphorus starvation), we identified another attribute shared by many genomes where transcription was highest under high phosphorus conditions. Taken together, we provide a novel framework for unravelling the functional dynamics of uncultivated microorganisms by assigning trait-attributes through genome-resolved time-series metatranscriptomics.

Reevaluation of the Phylogenetic Diversity and Global Distribution of the Genus "Candidatus Accumulibacter"

“Candidatus Accumulibacter” was the first microorganism identified as a polyphosphate-accumulating organism (PAO) important for phosphorus removal from wastewater. Members of this genus are diverse, and the current phylogeny and taxonomic framework appear complicated, with most publicly available genomes classified as “Candidatus Accumulibacter phosphatis,” despite notable phylogenetic divergence. The ppk1 marker gene allows for a finer-scale differentiation into different “types” and “clades”; nevertheless, taxonomic assignments remain inconsistent across studies. Therefore, a comprehensive reevaluation is needed to establish a common understanding of this genus, in terms of both naming and basic conserved physiological traits. Here, we provide this reassessment using a comparison of genome, ppk1, and 16S rRNA gene-based approaches from comprehensive data sets. We identified 15 novel species, along with “Candidatus Accumulibacter phosphatis,” “Candidatus Accumulibacter delftensis,” and “Candidatus Accumulibacter aalborgensis.” To compare the species in situ, we designed new species-specific fluorescence in situ hybridization (FISH) probes and revealed their morphology and arrangement in activated sludge. Based on the MiDAS global survey, “Ca. Accumulibacter” species were widespread in wastewater treatment plants (WWTPs) with phosphorus removal, indicating process design as a major driver for their abundance. Genome mining for PAO-related pathways and FISH-Raman microspectroscopy confirmed the potential for PAO metabolism in all “Ca. Accumulibacter” species, with detection in situ of the typical PAO storage polymers. Genome annotation further revealed differences in the nitrate/nitrite reduction pathways. This provides insights into the niche differentiation of these lineages, potentially explaining their coexistence in the same ecosystem while contributing to overall phosphorus and nitrogen removal.

IMPORTANCE “Candidatus Accumulibacter” is the most studied PAO, with a primary role in biological nutrient removal. However, the species-level taxonomy of this lineage is convoluted due to the use of different phylogenetic markers or genome sequencing approaches. Here, we redefined the phylogeny of these organisms, proposing a comprehensive approach which could be used to address the classification of other diverse and uncultivated lineages. Using genome-resolved phylogeny, compared to phylogeny based on the 16S rRNA gene and other phylogenetic markers, we obtained a higher-resolution taxonomy and established a common understanding of this genus. Furthermore, genome mining of genes and pathways of interest, validated in situ by application of a new set of FISH probes and Raman microspectroscopy, provided additional high-resolution metabolic insights into these organisms.

Enzymatic Depletion of Mitochondrial Inorganic Polyphosphate (polyP) Increases the Generation of Reactive Oxygen Species (ROS) and the Activity of the Pentose Phosphate Pathway (PPP) in Mammalian Cells

Inorganic polyphosphate (polyP) is an ancient biopolymer that is well preserved throughout evolution and present in all studied organisms. In mammals, it shows a high co-localization with mitochondria, and it has been demonstrated to be involved in the homeostasis of key processes within the organelle, including mitochondrial bioenergetics. However, the exact extent of the effects of polyP on the regulation of cellular bioenergetics, as well as the mechanisms explaining these effects, still remain poorly understood. Here, using HEK293 mammalian cells under Wild-type (Wt) and MitoPPX (cells enzymatically depleted of mitochondrial polyP) conditions, we show that depletion of polyP within mitochondria increased oxidative stress conditions. This is characterized by enhanced mitochondrial O2- and intracellular H2O2 levels, which may be a consequence of the dysregulation of oxidative phosphorylation (OXPHOS) that we have demonstrated in MitoPPX cells in our previous work. These findings were associated with an increase in basal peroxiredoxin-1 (Prx1), superoxide dismutase-2 (SOD2), and thioredoxin (Trx) antioxidant protein levels. Using 13C-NMR and immunoblotting, we assayed the status of glycolysis and the pentose phosphate pathway (PPP) in Wt and MitoPPX cells. Our results show that MitoPPX cells display a significant increase in the activity of the PPP and an increase in the protein levels of transaldolase (TAL), which is a crucial component of the non-oxidative phase of the PPP and is involved in the regulation of oxidative stress. In addition, we observed a trend towards increased glycolysis in MitoPPX cells, which corroborates our prior work. Here, for the first time, we show the crucial role played by mitochondrial polyP in the regulation of mammalian redox homeostasis. Moreover, we demonstrate a significant effect of mitochondrial polyP on the regulation of global cellular bioenergetics in these cells.

Disentangle genus microdiversity within a complex microbial community by using a multi-distance long-read binning method: example of Candidatus Accumulibacter

Complete genomes can be recovered from metagenomes by assembling and binning DNA sequences into metagenome assembled genomes (MAGs). Yet, the presence of microdiversity can hamper the assembly and binning processes, possibly yielding chimeric, highly fragmented and incomplete genomes. Here, the metagenomes of four samples of aerobic granular sludge bioreactors containing Candidatus (Ca.) Accumulibacter, a phosphate-accumulating organism of interest for wastewater treatment, were sequenced with both PacBio and Illumina. Different strategies of genome assembly and binning were investigated, including published protocols and a binning procedure adapted to the binning of long contigs (MuLoBiSC). Multiple criteria were considered to select the best strategy for Ca. Accumulibacter, whose multiple strains in every sample represent a challenging microdiversity. In this case, the best strategy relies on long-read only assembly and a custom binning procedure including MuLoBiSC in metaWRAP. Several high-quality Ca. Accumulibacter MAGs, including a novel species, were obtained independently from different samples. Comparative genomic analysis showed that MAGs retrieved in different samples harbour genomic rearrangements in addition to accumulation of point mutations. The microdiversity of Ca. Accumulibacter, likely driven by mobile genetic elements, causes major difficulties in recovering MAGs, but it is also a hallmark of the panmictic lifestyle of these bacteria.

Denitrification performance and mechanism of denitrification biofilm reactor based on carbon-nitrate counter-diffusional

This study researched denitrification performance and mechanism of denitrification biofilm reactor with different HRTs and carbon sources dosages. Experimental group (EG) had better nitrate and COD removal performance than control group (CG) with different HRTs or carbon doses, and the maximum nitrate-to-nitrite transformation ratio (NTR) of them reached 7.91 ± 1.60% and 17.50 ± 1.92%, respectively. Because organic carbon sources were added to the carrier's interior in EG, forming high local concentrations in biofilms and counter-diffusional with nitrate. By contrast, carbon sources and nitrate were provided from the aqueous phase in CG. Thus, the EG system has more active regions of the biofilm than CG. In addition, EG had higher proportions of microorganisms and enzymes related to denitrification and carbon metabolism. The most dominant phylum, genus, and species were Proteobacteria, Thaurea, and Thauera_sp._27, respectively. The transcript of acetyl-CoA synthetase (K01895) and denitrification (M00529) was mainly originated from unclassified_g__Pseudomonas and unclassified_g__Thauera, respectively.

Glycine adversely affects enhanced biological phosphorus removal

Enhanced biological phosphorus removal (EBPR) is used extensively in full-scale wastewater treatment plants for the removal of phosphorus. Despite previous evidence showing that glycine is a carbon source for a certain lineage of polyphosphate accumulating organisms (PAOs) such as Tetrasphaera, it is still unknown whether glycine can support EBPR. We observed an overall adverse effect of glycine on EBPR using activated sludge from both full-scale wastewater treatment plants and lab-scale reactors harboring distant and diverse PAOs and glycogen accumulating organisms (GAOs), including Candidatus Accumulibacter, Thiothrix, Tetrasphaera, Dechloromonas, Ca. Competibacter, and Defluviicoccus, among others. Glycine induced phosphorus (P) release under anaerobic conditions without being effectively taken up by cells. The induced P release rate correlated with glycine concentration in the range of 10 to 50 mg C/L. PAOs continued to release P in the presence of glycine under aerobic conditions without any evident P uptake. Under mixed carbon conditions, the occurrence of glycine did not seem to affect acetate uptake; however, it significantly reduced the rate of P uptake in the aerobic phase. Overall, glycine did not appear to be an effective carbon source for a majority of PAOs and GAOs in full-scale and lab-scale systems, and neither did other community members utilize glycine under anaerobic or aerobic conditions. Metatranscriptomic analysis showed the transcription of glycine cleavage T, P and H protein genes, but not of the L protein or the downstream genes in the glycine cleavage pathway, suggesting barriers to metabolizing glycine. The high transcription of a gene encoding a drug/metabolite transporter suggests a potential efflux mechanism, where glycine transported into the cells is in turn exported at the expense of ATP, resulting in P release without affecting the glycine concentration in solution. The ability of glycine to induce P release without cellular uptake suggests a way to effectively recover P from P-enriched waste sludge.

"Candidatus Dechloromonas phosphoritropha" and "Ca. D. phosphorivorans", novel polyphosphate accumulating organisms abundant in wastewater treatment systems

Members of the genus Dechloromonas are often abundant in enhanced biological phosphorus removal (EBPR) systems and are recognized putative polyphosphate accumulating organisms (PAOs), but their role in phosphate removal is still unclear. Here, we used 16S rRNA gene sequencing and fluorescence in situ hybridization (FISH) to investigate the abundance and distribution of Dechloromonas spp. in Danish and global wastewater treatment plants. The two most abundant species worldwide revealed in situ dynamics of important intracellular storage polymers, measured by FISH-Raman in activated sludge from four full-scale EBPR plants and from a lab-scale reactor fed with different substrates. Moreover, seven distinct Dechloromonas species were determined from a set of ten high-quality metagenome-assembled genomes (MAGs) from Danish EBPR plants, each encoding the potential for polyphosphate (poly-P), glycogen, and polyhydroxyalkanoates (PHA) accumulation. The two species exhibited an in situ phenotype in complete accordance with the metabolic information retrieved by the MAGs, with dynamic levels of poly-P, glycogen, and PHA during feast-famine anaerobic-aerobic cycling, legitimately placing these microorganisms among the important PAOs. They are potentially involved in denitrification showing niche partitioning within the genus and with other important PAOs. As no isolates are available for the two species, we propose the names Candidatus Dechloromonas phosphoritropha and Candidatus Dechloromonas phosphorivorans.

Oligotyping and metagenomics reveal distinct Candidatus Accumulibacter communities in side-stream versus conventional full-scale enhanced biological phosphorus removal (EBPR) systems

Candidatus Accumulibacter phosphatis (CAP) and its clade-level micro-diversity has been associated with and implicated in functional differences in phosphorus removal performance in enhanced biological phosphorus removal (EBPR) systems. Side-stream EBPR (S2EBPR) is an emerging process that has been shown to present a suite of advantages over the conventional EBPR design, however, large knowledge gaps remain in terms of its underlying ecological mechanisms. Here, we compared and revealed the higher-resolution differences in microbial ecology of CAP between a full-scale side-stream EBPR configuration and a conventional A2O EBPR process that were operated in parallel and with the same influent feed. Even though the relative abundance of CAP, revealed by 16S rRNA gene amplicon sequencing, was similar in both treatment trains, a clade-level analysis, using combined 16S rRNA-gene based amplicon sequencing and oligotyping analysis and metagenomics analysis, revealed the distinct CAP microdiversity between the S2EBPR and A2O configurations that likely attributed to the improved performance in S2EBPR in comparison to conventional EBPR. Furthermore, genome-resolved metagenomics enabled extraction of three metagenome-assembled genomes (MAGs) belonging to CAP clades IIB (RCAB4-2), IIC (RC14) and II (RC18), from full-scale EBPR sludge for the first time, including a distinct Ca. Accumulibacter clade that is dominant and associated only with the S2EBPR configuration. The results also revealed the temporally increasing predominance of RC14, which belonged to Clade IIC, during the implementation of the S2EBPR configuration. Finally, we also show the existence of previously uncharacterized diversity of clades of CAP, namely the clades IIB and as yet unidentified clade of type II, in full-scale EBPR communities, highlighting the unknown diversity of CAP communities in full-scale EBPR systems.

Prospects for multi-omics in the microbial ecology of water engineering

Advances in high-throughput sequencing technologies and bioinformatics approaches over almost the last three decades have substantially increased our ability to explore microorganisms and their functions - including those that have yet to be cultivated in pure isolation. Genome-resolved metagenomic approaches have enabled linking powerful functional predictions to specific taxonomical groups with increasing fidelity. Additionally, related developments in both whole community gene expression surveys and metabolite profiling have permitted for direct surveys of community-scale functions in specific environmental settings. These advances have allowed for a shift in microbiome science away from descriptive studies and towards mechanistic and predictive frameworks for designing and harnessing microbial communities for desired beneficial outcomes. Water engineers, microbiologists, and microbial ecologists studying activated sludge, anaerobic digestion, and drinking water distribution systems have applied various (meta)omics techniques for connecting microbial community dynamics and physiologies to overall process parameters and system performance. However, the rapid pace at which new omics-based approaches are developed can appear daunting to those looking to apply these state-of-the-art practices for the first time. Here, we review how modern genome-resolved metagenomic approaches have been applied to a variety of water engineering applications from lab-scale bioreactors to full-scale systems. We describe integrated omics analysis across engineered water systems and the foundations for pairing these insights with modeling approaches. Lastly, we summarize emerging omics-based technologies that we believe will be powerful tools for water engineering applications. Overall, we provide a framework for microbial ecologists specializing in water engineering to apply cutting-edge omics approaches to their research questions to achieve novel functional insights. Successful adoption of predictive frameworks in engineered water systems could enable more economically and environmentally sustainable bioprocesses as demand for water and energy resources increases.

Metabolic Differentiation of Co-occurring Accumulibacter Clades Revealed through Genome-Resolved Metatranscriptomics

Natural microbial communities consist of closely related taxa that may exhibit phenotypic differences and inhabit distinct niches. However, connecting genetic diversity to ecological properties remains a challenge in microbial ecology due to the lack of pure cultures across the microbial tree of life. "Candidatus Accumulibacter phosphatis" (Accumulibacter) is a polyphosphate-accumulating organism that contributes to the enhanced biological phosphorus removal (EBPR) biotechnological process for removing excess phosphorus from wastewater and preventing eutrophication from downstream receiving waters. Distinct Accumulibacter clades often coexist in full-scale wastewater treatment plants and laboratory-scale enrichment bioreactors and have been hypothesized to inhabit distinct ecological niches. However, since individual strains of the Accumulibacter lineage have not been isolated in pure culture to date, these predictions have been made solely on genome-based comparisons and enrichments with varying strain compositions. Here, we used genome-resolved metagenomics and metatranscriptomics to explore the activity of coexisting Accumulibacter strains in an engineered bioreactor environment. We obtained four high-quality genomes of Accumulibacter strains that were present in the bioreactor ecosystem, one of which is a completely contiguous draft genome scaffolded with long Nanopore reads. We identified core and accessory genes to investigate how gene expression patterns differed among the dominating strains. Using this approach, we were able to identify putative pathways and functions that may confer distinct functions to Accumulibacter strains and provide key functional insights into this biotechnologically significant microbial lineage. IMPORTANCE "Candidatus Accumulibacter phosphatis" is a model polyphosphate-accumulating organism that has been studied using genome-resolved metagenomics, metatranscriptomics, and metaproteomics to understand the EBPR process. Within the Accumulibacter lineage, several similar but diverging clades are defined by the shared sequence identity of the polyphosphate kinase (ppk1) locus. These clades are predicted to have key functional differences in acetate uptake rates, phage defense mechanisms, and nitrogen-cycling capabilities. However, such hypotheses have largely been made based on gene content comparisons of sequenced Accumulibacter genomes, some of which were obtained from different systems. Here, we performed time series genome-resolved metatranscriptomics to explore gene expression patterns of coexisting Accumulibacter clades in the same bioreactor ecosystem. Our work provides an approach for elucidating ecologically relevant functions based on gene expression patterns between closely related microbial populations.

Low Global Diversity of Candidatus Microthrix, a Troublesome Filamentous Organism in Full-Scale WWTPs

Candidatus Microthrix is one of the most common bulking filamentous microorganisms found in activated sludge wastewater treatment plants (WWTPs) across the globe. One species, Ca. M. parvicella, is frequently observed, but global genus diversity, as well as important aspects of its ecology and physiology, are still unknown. Here, we use the MiDAS ecosystem-specific 16S rRNA gene database in combination with amplicon sequencing of Danish and global WWTPs to investigate Ca. Microthrix spp. diversity, distribution, and factors affecting their global presence. Only two species were abundant across the world confirming low diversity of the genus: the dominant Ca. M. parvicella and an unknown species typically present along with Ca. M. parvicella, although usually in lower abundances. Both species were mostly found in Europe at low-to-moderate temperatures and their growth was favored in municipal WWTPs with advanced process designs. As no isolate is available for the novel species, we propose the name "Candidatus Microthrix subdominans." Ten high-quality metagenome-assembled genomes recovered from Danish WWTPs, including 6 representing the novel Ca. M. subdominans, demonstrated high genetic similarity between the two species with a likely preference for lipids, a putative capability to reduce nitrate and nitrite, and the potential to store lipids and poly-P. Ca. M. subdominans had a potentially more versatile metabolism including additional sugar transporters, higher oxygen tolerance, and the potential to use carbon monoxide as energy source. Newly designed fluorescence in situ hybridization probes revealed similar filamentous morphology for both species. Raman microspectroscopy was used to quantify the in situ levels of intracellular poly-P. Despite the observed similarities in their physiology (both by genomes and in situ), the two species showed different seasonal dynamics in Danish WWTPs through a 13-years survey, possibly indicating occupation of slightly different niches. The genomic information provides the basis for future research into in situ gene expression and regulation, while the new FISH probes provide a useful tool for further characterization in situ. This study is an important step toward understanding the ecology of Ca. Microthrix in WWTPs, which may eventually lead to optimization of control strategies for its growth in this ecosystem.

Integrated omics analyses reveal differential gene expression and potential for cooperation between denitrifying polyphosphate and glycogen accumulating organisms

Unusually high accumulation of the potent greenhouse gas nitrous oxide (N₂ O) has previously been documented in denitrifying biological phosphorus (P) removal bioprocesses, but the roles of differential denitrification gene expression patterns and ecological interactions between key functional groups in driving these emissions are not well understood. To address these knowledge gaps, we applied genome-resolved metagenomics and metatranscriptomics to a denitrifying bioprocess enriched in as-yet-uncultivated denitrifying polyphosphate accumulating organisms (PAOs) affiliated with Candidatus Accumulibacter. The six transcriptionally most active populations in the community included three co-occurring Accumulibacter strains affiliated with clades IF (a novel clade identified in this study), IA and IC, a competing glycogen accumulating organism (GAO) affiliated with Competibacteraceae (GAO1), a Gammaproteobacteria PR6 and an Anaerolineae CH7. Strongly elevated expression of nitrite reductase genes compared to nitrous oxide reductase genes was observed in the overall community and in Accumulibacter populations, suggesting a strong role for differential gene expression in driving N₂ O accumulation. Surprisingly, while ~90% of the nirS gene transcripts were expressed by the three co-occurring PAO populations, ~93% of the norB gene transcripts were expressed by GAO1 and ~75% of the norZ gene transcripts were mapped to PR6 and several other non-PAO flanking populations. This suggests the potential for cooperation between flanking populations and PAOs in reducing denitrification intermediates. Such cooperation may benefit the community by reducing the accumulation of toxic nitric oxide.

Production of nonulosonic acids in the extracellular polymeric substances of "Candidatus Accumulibacter phosphatis"

Nonulosonic acids (NulOs) are a family of acidic carbohydrates with a nine-carbon backbone, which include different related structures, such as sialic acids. They have mainly been studied for their relevance in animal cells and pathogenic bacteria. Recently, sialic acids have been discovered as an important compound in the extracellular matrix of virtually all microbial life and in "Candidatus Accumulibacter phosphatis", a well-studied polyphosphate-accumulating organism, in particular. Here, bioaggregates highly enriched with these bacteria (approx. 95% based on proteomic data) were used to study the production of NulOs in an enrichment of this microorganism. Fluorescence lectin-binding analysis, enzymatic quantification, and mass spectrometry were used to analyze the different NulOs present, showing a wide distribution and variety of these carbohydrates, such as sialic acids and bacterial NulOs, in the bioaggregates. Phylogenetic analysis confirmed the potential of "Ca. Accumulibacter" to produce different types of NulOs. Proteomic analysis showed the ability of "Ca. Accumulibacter" to reutilize and reincorporate these carbohydrates. This investigation points out the importance of diverse NulOs in non-pathogenic bacteria, which are normally overlooked. Sialic acids and other NulOs should be further investigated for their role in the ecology of "Ca. Accumulibacter" in particular, and biofilms in general. KEY POINTS: •"Ca. Accumulibacter" has the potential to produce a range of nonulosonic acids. •Mass spectrometry and lectin binding can reveal the presence and location of nonulosonic acids. •The role of nonulosonic acid in non-pathogenic bacteria needs to be studied in detail.

Connecting structure to function with the recovery of over 1000 high-quality metagenome-assembled genomes from activated sludge using long-read sequencing

Microorganisms play crucial roles in water recycling, pollution removal and resource recovery in the wastewater industry. The structure of these microbial communities is increasingly understood based on 16S rRNA amplicon sequencing data. However, such data cannot be linked to functional potential in the absence of high-quality metagenome-assembled genomes (MAGs) for nearly all species. Here, we use long-read and short-read sequencing to recover 1083 high-quality MAGs, including 57 closed circular genomes, from 23 Danish full-scale wastewater treatment plants. The MAGs account for ~30% of the community based on relative abundance, and meet the stringent MIMAG high-quality draft requirements including full-length rRNA genes. We use the information provided by these MAGs in combination with >13 years of 16S rRNA amplicon sequencing data, as well as Raman microspectroscopy and fluorescence in situ hybridisation, to uncover abundant undescribed lineages belonging to important functional groups.

Revealing taxon-specific heavy metal-resistance mechanisms in denitrifying phosphorus removal sludge using genome-centric metaproteomics

BACKGROUND

Denitrifying phosphorus removal sludge (DPRS) is widely adopted for nitrogen and phosphorus removal in wastewater treatment but faces threats from heavy metals. However, a lack of understanding of the taxon-specific heavy metal-resistance mechanisms hinders the targeted optimization of DPRS's robustness in nutrient removal.

RESULTS

We obtained 403 high- or medium-quality metagenome-assembled genomes from DPRS treated by elevating cadmium, nickel, and chromium pressure. Then, the proteomic responses of individual taxa under heavy metal pressures were characterized, with an emphasis on functions involving heavy metal resistance and maintenance of nutrient metabolism. When oxygen availability was constrained by high-concentration heavy metals, comammox Nitrospira overproduced highly oxygen-affinitive hemoglobin and electron-transporting cytochrome c-like proteins, underpinning its ability to enhance oxygen acquisition and utilization. In contrast, Nitrosomonas overexpressed ammonia monooxygenase and nitrite reductase to facilitate the partial nitrification and denitrification process for maintaining nitrogen removal. Comparisons between phosphorus-accumulating organisms (PAOs) demonstrated different heavy metal-resistance mechanisms adopted by Dechloromonas and Candidatus Accumulibacter, despite their high genomic similarities. In particular, Dechloromonas outcompeted the canonical PAO Candidatus Accumulibacter in synthesizing polyphosphate, a potential public good for heavy metal detoxification. The superiority of Dechloromonas in energy utilization, radical elimination, and damaged cell component repair also contributed to its dominance under heavy metal pressures. Moreover, the enrichment analysis revealed that functions involved in extracellular polymeric substance formation, siderophore activity, and heavy metal efflux were significantly overexpressed due to the related activities of specific taxa.

CONCLUSIONS

Our study demonstrates that heavy metal-resistance mechanisms within a multipartite community are highly heterogeneous between different taxa. These findings provide a fundamental understanding of how the heterogeneity of individual microorganisms contributes to the metabolic versatility and robustness of microbiomes inhabiting dynamic environments, which is vital for manipulating the adaptation of microbial assemblages under adverse environmental stimuli. Video abstract.

Dissolved organic matter byproducts from combined low-level alkaline and mechanical treatment of sludge: species, and their bioavailability

Low-level alkalinity (pH 9-10) coupled with ultrasonic or mechanical cutting with different energy input for obtaining carbon sources were tested for sludge pretreatment process before anaerobic sludge digestion. The differences between the primary sludge (PS) and waste activated sludge (WAS)-derived dissolved organic matter (DOM) species were evaluated for their bioavailability and affinity (in the form of amino acids) to the bio-nutrient removal (BNR) biomass. Soluble microbial by-product-like substances as the predominant DOM components in the raw PS and WAS increased by 23 and 22%, respectively, after low-level alkaline treatment (pH 9-10) and ultrasonication. In addition, the protein components were degraded further as free amino acids (FAAs). The sludge-derived aspartate, glutamate, followed by arginine were the most commonly used FAAs by the BNR biomass. The pattern of recovering this special sludge-derived carbon source to enhance P removal and recovery in the BNR process is depicted.

The phylogenetic and global distribution of bacterial polyhydroxyalkanoate bioplastic-degrading genes

Polyhydroxyalkanoates (PHAs) are a family of microbially made polyesters commercialized as biodegradable plastics. PHA production rates are predicted to increase as concerns around environmental plastic contamination and limited fossil fuel resources have increased the importance of biodegradable and bio-based plastic alternatives. Microbially produced PHA depolymerases are the key enzymes mediating PHA biodegradation, but only a few PHA depolymerases have been well-characterized and screens employing metagenomic sequence data are lacking. Here, we used 3078 metagenomes to analyse the distribution of PHA depolymerases in microbial communities from diverse aquatic, terrestrial and waste management systems. We significantly expand the recognized diversity of this protein family by screening 1914 Gb of sequence data and identifying 13 869 putative PHA depolymerases in 1295 metagenomes. Our results indicate that PHA depolymerases are unevenly distributed across environments. We predicted the highest frequency of PHA depolymerases in wastewater systems and the lowest in marine and thermal springs. In tandem, we screened 5290 metagenome-assembled genomes to describe the phylogenetic distribution of PHA depolymerases, which is substantially broader compared with current cultured representatives. The Proteobacteria and Bacteroidota are key lineages encoding PHA depolymerases, but PHA depolymerases were predicted from members of the Bdellovibrionota, Methylomirabilota, Actinobacteriota, Firmicutes, Spirochaetota, Desulfobacterota, Myxococcota and Planctomycetota.

Recent advances in understanding the ecophysiology of enhanced biological phosphorus removal

Enhanced biological phosphorus removal (EBPR) is an efficient, cost-effective, and sustainable method for removing excess phosphorus from wastewater. Polyphosphate accumulating organisms (PAOs) exhibit a unique physiology alternating between anaerobic conditions for uptake of carbon substrates and aerobic or anoxic conditions for phosphorus uptake. The implementation of high-throughput sequencing technologies and advanced molecular tools along with biochemical characterization has provided many new perspectives on the EBPR process. These approaches have helped identify a wide range of carbon substrates and electron acceptors utilized by PAOs that in turn influence interactions with microbial community members and determine overall phosphorus removal efficiency. In this review, we systematically discuss the microbial diversity and metabolic response to a range of environmental conditions and process control strategies in EBPR.

Microbiome of Seven Full-Scale Anaerobic Digestion Plants in South Korea: Effect of Feedstock and Operational Parameters

In this study, the microbiomes linked with the operational parameters in seven mesophilic full-scale AD plants mainly treating food waste (four plants) and sewage sludge (three plants) were analyzed. The results obtained indicated lower diversity and evenness of the microbial population in sludge digestion (SD) plants compared to food digestion (FD) plants. Candidatus Accumulibacter dominated (up to 42.1%) in SD plants due to microbial immigration from fed secondary sludge (up to 89%). Its potential activity in SD plants was correlated to H2 production, which was related to the dominance of hydrogenotrophic methanogens (Methanococcus). In FD plants, a balance between the hydrogenotrophic and methylotrophic pathways was found, while Flavobacterium and Levilinea played an important role during acidogenesis. Levilinea also expressed sensitivity to ammonia in FD plants. The substantial differences in hydraulic retention time (HRT), organic loading rate (OLR), and total ammonium nitrogen (TAN) among the studied FD plants did not influence the archaeal methane production pathway. In addition, the bacterial genera responsible for acetate production through syntrophy and homoacetogenesis (Smithella, Treponema) were present in all the plants studied.

Restoring degraded microbiome function with self-assembled communities

The natural microbial functions of many soils are severely degraded. Current state-of-the-art technology to restore these functions is through the isolation, screening, formulation and application of microbial inoculants and synthetic consortia. These approaches have inconsistent success, in part due to the incompatibility between the biofertilizer, crop, climate, existing soil microbiome and physicochemical characteristics of the soils. Here, we review the current state of the art in biofertilization and identify two key deficiencies in current strategies: the difficulty in designing complex multispecies biofertilizers and the bottleneck in scaling the production of complex multispecies biofertilizers. To address the challenge of producing scalable, multispecies biofertilizers, we propose to merge ecological theory with bioprocess engineering to produce 'self-assembled communities' enriched for particular functional guilds and adapted to a target soil and host plant. Using the nitrogen problem as an anchor, we review relevant ecology (microbial, plant and environmental), as well as reactor design strategies and operational parameters for the production of functionally enriched self-assembled communities. The use of self-assembled communities for biofertilization addresses two major hurdles in microbiome engineering: the importance of enriching microbes indigenous to (and targeted for) a specific environment and the recognized potential benefits of microbial consortia over isolates (e.g. functional redundancy). The proposed community enrichment model could also be instrumental for other microbial functions such as phosphorus solubilization, plant growth promotion or disease suppression.

Revealing the Metabolic Flexibility of "Candidatus Accumulibacter phosphatis" through Redox Cofactor Analysis and Metabolic Network Modeling

Environmental fluctuations in the availability of nutrients lead to intricate metabolic strategies. "Candidatus Accumulibacter phosphatis," a polyphosphate-accumulating organism (PAO) responsible for enhanced biological phosphorus removal (EBPR) from wastewater treatment systems, is prevalent in aerobic/anaerobic environments. While the overall metabolic traits of these bacteria are well described, the nonavailability of isolates has led to controversial conclusions on the metabolic pathways used. In this study, we experimentally determined the redox cofactor preferences of different oxidoreductases in the central carbon metabolism of a highly enriched "Ca Accumulibacter phosphatis" culture. Remarkably, we observed that the acetoacetyl coenzyme A reductase engaged in polyhydroxyalkanoate (PHA) synthesis is NADH preferring instead of showing the generally assumed NADPH dependency. This allows rethinking of the ecological role of PHA accumulation as a fermentation product under anaerobic conditions and not just a stress response. Based on previously published metaomics data and the results of enzymatic assays, a reduced central carbon metabolic network was constructed and used for simulating different metabolic operating modes. In particular, scenarios with different acetate-to-glycogen consumption ratios were simulated, which demonstrated optima using different combinations of glycolysis, glyoxylate shunt, or branches of the tricarboxylic acid (TCA) cycle. Thus, optimal metabolic flux strategies will depend on the environment (acetate uptake) and on intracellular storage compound availability (polyphosphate/glycogen). This NADH-related metabolic flexibility is enabled by the NADH-driven PHA synthesis. It allows for maintaining metabolic activity under various environmental substrate conditions, with high carbon conservation and lower energetic costs than for NADPH-dependent PHA synthesis. Such (flexible) metabolic redox coupling can explain the competitiveness of PAOs under oxygen-fluctuating environments.IMPORTANCE Here, we demonstrate how microbial storage metabolism can adjust to a wide range of environmental conditions. Such flexibility generates a selective advantage under fluctuating environmental conditions. It can also explain the different observations reported in PAO literature, including the capacity of "Ca Accumulibacter phosphatis" to act like glycogen-accumulating organisms (GAOs). These observations stem from slightly different experimental conditions, and controversy arises only when one assumes that metabolism can operate only in a single mode. Furthermore, we also show how the study of metabolic strategies is possible when combining omics data with functional cofactor assays and modeling. Genomic information can only provide the potential of a microorganism. The environmental context and other complementary approaches are still needed to study and predict the functional expression of such metabolic potential.

An NADH preferring acetoacetyl-CoA reductase is engaged in poly-3-hydroxybutyrate accumulation in Escherichia coli

Oxygen supply implies higher production cost and reduction of maximum theoretical yields. Thus, generation of fermentation products is more cost-effective. Aiming to find a key piece for the production of (poly)-3-hydroxybutyrate (PHB) as a fermentation product, here we characterize an acetoacetyl-CoA reductase, isolated from a Candidatus Accumulibacter phosphatis-enriched mixed culture, showing a (k_cat^NADH/K_M^NADH)/(k_cat^NADPH/K_M^NADPH)>500. Further kinetic analyses indicate that, at physiological concentrations, this enzyme clearly prefers NADH, presenting the strongest NADH preference so far observed among the acetoacetyl-CoA reductases. Structural and kinetic analyses indicate that residues between E37 and P41 have an important role for the observed NADH preference. Moreover, an operon was assembled combining the phaCA genes from Cupriavidus necator and the gene encoding for this NADH-preferring acetoacetyl-CoA reductase. Escherichia coli cells expressing that assembled operon showed continuous accumulation of PHB under oxygen limiting conditions and PHB titer increased when decreasing the specific oxygen consumption rate. Taken together, these results show that it is possible to generate PHB as a fermentation product in E. coli, opening opportunities for further protein/metabolic engineering strategies envisioning a more efficient anaerobic production of PHB.

Emerging Trends in Biological Treatment of Wastewater From Unconventional Oil and Gas Extraction

Unconventional oil and gas exploration generates an enormous quantity of wastewater, commonly referred to as flowback and produced water (FPW). Limited freshwater resources and stringent disposal regulations have provided impetus for FPW reuse. Organic and inorganic compounds released from the shale/brine formation, microbial activity, and residual chemicals added during hydraulic fracturing bestow a unique as well as temporally varying chemical composition to this wastewater. Studies indicate that many of the compounds found in FPW are amenable to biological degradation, indicating biological treatment may be a viable option for FPW processing and reuse. This review discusses commonly characterized contaminants and current knowledge on their biodegradability, including the enzymes and organisms involved. Further, a perspective on recent novel hybrid biological treatments and application of knowledge gained from omics studies in improving these treatments is explored.

Impact of hydraulic retention time on community assembly and function of photogranules for wastewater treatment

Photogranules are dense, spherical agglomerates of cyanobacteria, microalgae and non-phototrophic microorganisms that have considerable advantages in terms of harvesting and nutrient removal rates for light driven wastewater treatment processes. This ecosystem is poorly understood in terms of the microbial community structure and the response of the community to changing abiotic conditions. To get a better understanding, we investigated the effect of hydraulic retention time (HRT) on photogranule formation and community assembly over a period of 148 days. Three laboratory bioreactors were inoculated with field samples from various locations in the Netherlands and operated in sequencing batch mode. The bioreactors were operated at four different HRTs (2.00, 1.00, 0.67, 0.33 days), while retaining the same solid retention time of 7 days. A microbial community with excellent settling characteristics (95-99% separation efficiency) was established within 2-5 weeks. The observed nutrient uptake rates ranged from 24 to 90 mg_N L^-1 day^-1 and from 3.1 to 5.4 mg_P L^-1 day^-1 depending on the applied HRT. The transition from single-cell suspension culture to floccular agglomeration to granular sludge was monitored by microscopy and 16S/18S sequencing. In particular, two important variables for driving aggregation and granulation, and for the structural integrity of photogranules were identified: 1. Extracellular polymeric substances (EPS) with high protein to polysaccharide ratio and 2. specific microorganisms. The key players were found to be the cyanobacteria Limnothrix and Cephalothrix, the colony forming photosynthetic eukaryotes within Chlamydomonadaceae, and the biofilm producing bacteria Zoogloea and Thauera. Knowing the makeup of the microbial community and the operational conditions influencing granulation and bioreactor function is crucial for successful operation of photogranular systems.

Metabolic Traits of Candidatus Accumulibacter clade IIF Strain SCELSE-1 Using Amino Acids As Carbon Sources for Enhanced Biological Phosphorus Removal

Despite recent evidence from full-scale plants suggesting that Candidatus Accumulibacter may be capable of using amino acids, this metabolic trait has never been confirmed in a bioreactor experiment. Here we show that an enriched culture of Ca. Accumulibacter clade IIF strain SCELSE-1 could metabolize 11 of 20 α-amino acids, with aspartate, glutamate, asparagine, and glutamine resulting in the highest phosphorus removal. The anaerobic uptake of aspartate and glutamate was achieved through a glutamate/aspartate-proton symporter fully powered by the proton motive force (PMF). Under anaerobic conditions aspartate was deaminized and routed into core carbon metabolic pathways to form polyhydroxyalkanoates (PHA). The lack of genes encoding NADH dependent isocitrate dehydrogenase in the Ca. Accumulibacter genome resulted in a kinetic barrier for glutamate to be channelled to the TCA cycle. Glutamate was stored as glutamate polymer. When amino acids (aspartate or glutamate) and acetate were supplied together, Ca. Accumulibacter took up both carbon sources simultaneously, with the uptake rate of each carbon source largely preserved. Overall energy savings (up to 17%) were achieved under mixed carbon scenarios, due to the ability of Ca. Accumulibacter to rearrange its anaerobic carbon metabolism based on the reducing power, PMF and ATP balance.

Compositional variation between high and low prokaryotic diversity coral reef biotopes translates to different predicted metagenomic gene content

In a previous study, we identified host species that housed high and low diversity prokaryotic communities. In the present study, we expand on this and assessed the prokaryotic communities associated with seawater, sediment and 11 host species from 7 different phyla in a Taiwanese coral reef setting. The host taxa sampled included hard, octo- and black corals, molluscs, bryozoans, flatworms, fish and sea urchins. There were highly significant differences in composition among host species and all host species housed distinct communities from those found in seawater and sediment. In a hierarchical clustering analysis, samples from all host species, with the exception of the coral Galaxea astreata, formed significantly supported clusters. In addition to this, the coral G. astreata and the bryozoan Triphyllozoon inornatum on the one hand and the coral Tubastraea coccinea, the hermit crab Calcinus laevimanus and the flatworm Thysanozoon nigropapillosum on the other formed significantly supported clusters. In addition to composition, there were highly pronounced differences in richness and evenness among host species from the most diverse species, the bryozoan T. inornatum at 2518 ± 240 OTUs per 10,000 sequences to the least diverse species, the octocoral Cladiella sp. at 142 ± 14 OTUs per 10,000 sequences. In line with the differences in composition, there were significant differences in predicted metagenomic gene counts among host species. Furthermore, there were pronounced compositional and predicted functional differences between high diversity hosts (Liolophura japonica, G. astreata, T. coccinea, C. laevimanus, T. inornatum) and low diversity hosts (Antipathes sp., Pomacentrus coelestis, Modiolus auriculatus, T. nigropapillosum, Cladiella sp. and Diadema savigny). In particular, we found that all tested low diversity hosts were predicted to be enriched for the phosphotransferase system compared to high diversity hosts.

Cadmium (II) alters the microbial community structure and molecular ecological network in activated sludge system

Cadmium (II) can potentially alter the microbial community structure and molecular ecological network in activated sludge systems. In this study, we used Illumina sequencing combined with an RMT-based network approach to show the response of the microbial community and its network structure to Cd (II) in activated sludge systems. The results demonstrated that 1 mg/L Cd (II) did not have chronic negative effects on chemical oxygen demand (COD) reduction and denitrification processes, but negatively affected the nitrification process and phosphorus removal. In contrast, 10 mg/L Cd (II) adversely affected both COD and nutrient removal, and reduced the microbial diversity and changed the overall microbial community structure. The relative abundances of Nitrosomonadaceae, Nitrospira, Accumulibacter and Acinetobacter, which are involved in nitrogen removal, significantly decreased with increases in the Cd (II) concentration. In addition, molecular ecological network analysis showed that the networks sizes in the presence of higher levels of Cd (II) were smaller than in the control, but the nodes were more closely connected with neighbors. These shifts in bacterial abundance and the bacterial network structure may be responsible for the deterioration of COD and nutrient removal. Overall, this study provides new insights into the effects of Cd (II) on the bacterial community and its interactions in activated sludge systems.

Mining traits for the enrichment and isolation of not-yet-cultured populations

BACKGROUND

The lack of pure cultures limits our understanding into 99% of bacteria. Proper interpretation of the genetic and the transcriptional datasets can reveal clues for the enrichment and even isolation of the not-yet-cultured populations. Unraveling such information requires a proper mining method.

RESULTS

Here, we present a method to infer the hidden traits for the enrichment of not-yet-cultured populations. We demonstrate this method using Candidatus Accumulibacter. Our method constructs a whole picture of the carbon, electron, and energy flows in the not-yet-cultured populations from the genomic datasets. Then, it decodes the coordination across three flows from the transcriptional datasets. Based on it, our method diagnoses the status of the not-yet-cultured populations and provides strategy to optimize the enrichment systems.

CONCLUSION

Our method could shed light to the exploration into the bacterial dark matter in the environments.

Genome-centric metagenomics resolves microbial diversity and prevalent truncated denitrification pathways in a denitrifying PAO-enriched bioprocess

Denitrification is the stepwise microbial reduction of nitrate or nitrite (NO₂^-) to nitrogen gas via the obligate intermediates nitric oxide (NO) and nitrous oxide (N₂O). Substantial N₂O accumulation has been reported in denitrifying enhanced biological phosphorus removal (EBPR) bioreactors enriched in denitrifying polyphosphate accumulating organisms (DPAOs), but little is known about underlying mechanisms for N₂O generation, prevalence of complete versus truncated denitrification pathways, or the impact of NO₂^- feed on DPAO-enriched consortia. To address this knowledge gap, we employed genome-resolved metagenomics to investigate nitrogen transformation potential in a NO₂^- fed denitrifying EBPR bioreactor enriched in Candidatus Accumulibacter and prone to N₂O accumulation. Our analysis yielded 41 near-complete metagenome-assembled genomes (MAGs), including two co-occurring Accumulibacter strains affiliated with clades IA and IC (the first published genome from this clade) and 39 non-PAO flanking bacterial genomes. The dominant Accumulibacter clade IA encoded genes for complete denitrification, while the lower abundance Accumulibacter clade IC harbored all denitrification genes except for a canonical respiratory NO reductase. Analysis of the 39 non-PAO MAGs revealed a high prevalence of taxa harboring an incomplete denitrification pathway. Of the 27 MAGs harboring capacity for at least one step in the denitrification pathway, 10 were putative N₂O producers lacking N₂O reductase, 16 were putative N₂O reducers that lacked at least one upstream denitrification gene, and only one harbored a complete denitrification pathway. We also documented increasing abundance over the course of reactor operation of putative N₂O producers. Our results suggest that the unusually high levels of N₂O production observed in this Accumulibacter-enriched consortium are linked in part to the selection for non-PAO flanking microorganisms with truncated denitrification pathways.