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

From Gene to Structure: Unraveling Genomic Dark Matter in Ca. Accumulibacter

"Candidatus Accumulibacter" is a unique and pivotal genus of polyphosphate-accumulating organisms prevalent in wastewater treatment plants and plays mainstay roles in the global phosphorus cycle. However, the efforts to fully understand their genetic and metabolic characteristics are largely hindered by major limitations in existing sequence-based annotation methods. Here, we reported an integrated approach combining pangenome analysis, protein structure prediction and clustering, and meta-omic characterization, to uncover genetic and metabolic traits previously unexplored for Ca. Accumulibacter. The identification of a previously overlooked pyrophosphate-fructose 6-phosphate 1-phosphotransferase gene (pfp) suggested that all Ca. Accumulibacter encoded a complete Embden-Meyerhof-Parnas pathway. A homologue of the phosphate-specific transport system accessory protein (PhoU) was suggested to be an inorganic phosphate transport (Pit) accessory protein (Pap) conferring effective and efficient phosphate transport. Additional lineage members were found to encode complete denitrification pathways. A pipeline was built, generating a pan-Ca. Accumulibacter annotation reference database, covering >200,000 proteins and their encoding genes. Benchmarking on 27 Ca. Accumulibacter genomes showed major improvement in the average annotation coverage from 51% to 82%. This pipeline is readily applicable to diverse cultured and uncultured bacteria to establish high-coverage annotation reference databases, facilitating the exploration of genomic dark matter in the bacterial domain.

Distinct microdiversity of phosphate accumulating organisms (PAOs) between side-stream and conventional enhanced biological phosphorus removal (EBPR) systems with performance implications

Polyphosphate Accumulating Organisms (PAOs) microdiversity is a key factor to elucidate the mechanisms involved in the side-stream enhanced biological phosphorus removal (S2EBPR) systems, which has been shown to improve the process stability over conventional EBPR. However, fast, effective and cost-efficient methods to resolve PAO microdiversity in real-world activate sludge samples is still in absence. In this study, we applied oligotyping analysis following the regular 16S rRNA gene amplicon sequencing standard operation pipeline (SOP) to resolve subgenus-level PAO oligotypes, which cannot be achieved using traditional 16S rRNA sequencing SOP. The identified oligotype profiles of PAO-containing genera Ca. Accumulibacter, Tetrasphaera and Comamonas showed distinguished community-level differences across 12 water resource recovery facilities (WRRFs), which would not be revealed at the genus level. The WRRF-level differences were observed larger than the temporal differences in the same WRRF, indicating intrinsic sub-genus level microdiversity fingerprint between EBPR/S2EBPR systems. The identified oligotypes can be associated with known PAO clades phylogenetically, suggesting that oligotyping can suffice as a fast and cost-efficient approach for PAO microdiversity profiling. In addition, network analysis can be used to identify coexistence patterns between oligotypes with respect to EBPR/S2EBPR configurations and performance, enabling more detailed analysis between EBPR system performance and PAOs microdiversity. Correlation analyses between oligotype profiles and key EBPR performance parameters revealed potential different biological functional traits among these PAO species with P-removal performance implications.

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.

Refinement of the "Candidatus Accumulibacter" genus based on metagenomic analysis of biological nutrient removal (BNR) pilot-scale plants operated with reduced aeration

Members of the "Candidatus Accumulibacter" genus are widely studied as key polyphosphate-accumulating organisms (PAOs) in biological nutrient removal (BNR) facilities performing enhanced biological phosphorus removal (EBPR). This diverse lineage includes 18 "Ca. Accumulibacter" species, which have been proposed based on the phylogenetic divergence of the polyphosphate kinase 1 (ppk1) gene and genome-scale comparisons of metagenome-assembled genomes (MAGs). Phylogenetic classification based on the 16S rRNA genetic marker has been difficult to attain because most "Ca. Accumulibacter" MAGs are incomplete and often do not include the rRNA operon. Here, we investigate the "Ca. Accumulibacter" diversity in pilot-scale treatment trains performing BNR under low dissolved oxygen (DO) conditions using genome-resolved metagenomics. Using long-read sequencing, we recovered medium- and high-quality MAGs for 5 of the 18 "Ca. Accumulibacter" species, all with rRNA operons assembled, which allowed a reassessment of the 16S rRNA-based phylogeny of this genus and an analysis of phylogeny based on the 23S rRNA gene. In addition, we recovered a cluster of MAGs that based on 16S rRNA, 23S rRNA, ppk1, and genome-scale phylogenetic analyses do not belong to any of the currently recognized "Ca. Accumulibacter" species for which we propose the new species designation "Ca. Accumulibacter jenkinsii" sp. nov. Relative abundance evaluations of the genus across all pilot plant operations revealed that regardless of the operational mode, "Ca. A. necessarius" and "Ca. A. propinquus" accounted for more than 40% of the "Ca. Accumulibacter" community, whereas the newly proposed "Ca. A. jenkinsii" accounted for about 5% of the "Ca. Accumulibacter" community.IMPORTANCEOne of the main drivers of energy use and operational costs in activated sludge processes is the amount of oxygen provided to enable biological phosphorus and nitrogen removal. Wastewater treatment facilities are increasingly considering reduced aeration to decrease energy consumption, and whereas successful BNR has been demonstrated in systems with minimal aeration, an adequate understanding of the microbial communities that facilitate nutrient removal under these conditions is still lacking. In this study, we used genome-resolved metagenomics to evaluate the diversity of the "Candidatus Accumulibacter" genus in pilot-scale plants operating with minimal aeration. We identified the "Ca. Accumulibacter" species enriched under these conditions, including one novel species for which we propose "Ca. Accumulibacter jenkinsii" sp. nov. as its designation. Furthermore, the MAGs obtained for five additional "Ca. Accumulibacter" species further refine the phylogeny of the "Ca. Accumulibacter" genus and provide new insight into its diversity within unconventional biological nutrient removal systems.

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.

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.

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.

CRISPR-Cas phage defense systems and prophages in Candidatus Accumulibacter

Candidatus Accumulibacter plays a major role in enhanced biological phosphorus removal (EBPR) from wastewater. Although bacteriophages have been shown to represent fatal threats to Ca. Accumulibacter organisms and thus interfere with the stability of the EBPR process, little is known about the ability of different Ca. Accumulibacter strains to resist phage infections. We conducted a systematic analysis of the occurrence and characteristics of clustered regularly interspaced short palindromic repeats and associated proteins (CRISPR-Cas) systems and prophages in Ca. Accumulibacter lineage members (43 in total, including 10 newly recovered genomes). Results indicate that 28 Ca. Accumulibacter genomes encode CRISPR-Cas systems. They were likely acquired via horizontal gene transfer, conveying a distinct adaptivity to phage predation to different Ca. Accumulibacter members. Major differences in the number of spacers show the unique phage resistance of these members. A comparison of the spacers in closely related Ca. Accumulibacter members from distinct geographical locations indicates that habitat isolation may have resulted in the acquisition of resistance to different phages by different Ca. Accumulibacter. Long-term operation of three laboratory-scale EBPR bioreactors revealed high relative abundances of Ca. Accumulibacter with CRISPSR-Cas systems. Their specific resistance to phages in these reactors was indicated by spacer analysis. Metatranscriptomic analyses showed the activation of the CRISPR-Cas system under both anaerobic and aerobic conditions. Additionally, 133 prophage regions were identified in 43 Ca. Accumulibacter genomes. Twenty-seven of them (in 19 genomes) were potentially active. Major differences in the occurrence of CRISPR-Cas systems and prophages in Ca. Accumulibacter will lead to distinct responses to phage predation. This study represents the first systematic analysis of CRISPR-Cas systems and prophages in the Ca. Accumulibacter lineage, providing new perspectives on the potential impacts of phages on Ca. Accumulibacter and EBPR systems.

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.

Integrating Genome-Resolved Metagenomics with Trait-Based Process Modeling to Determine Biokinetics of Distinct Nitrifying Communities within Activated Sludge

Conventional bioprocess models for wastewater treatment are based on aggregated bulk biomass concentrations and do not incorporate microbial physiological diversity. Such a broad aggregation of microbial functional groups can fail to predict ecosystem dynamics when high levels of physiological diversity exist within trophic guilds. For instance, functional diversity among nitrite-oxidizing bacteria (NOB) can obfuscate engineering strategies for their out-selection in activated sludge (AS), which is desirable to promote energy-efficient nitrogen removal. Here, we hypothesized that different NOB populations within AS can have different physiological traits that drive process performance, which we tested by estimating biokinetic growth parameters using a combination of highly replicated respirometry, genome-resolved metagenomics, and process modeling. A lab-scale AS reactor subjected to a selective pressure for over 90 days experienced resilience of NOB activity. We recovered three coexisting Nitrospira population genomes belonging to two sublineages, which exhibited distinct growth strategies and underwent a compositional shift following the selective pressure. A trait-based process model calibrated at the NOB genus level better predicted nitrite accumulation than a conventional process model calibrated at the NOB guild level. This work demonstrates that trait-based modeling can be leveraged to improve our prediction, control, and design of functionally diverse microbiomes driving key environmental biotechnologies.

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.

Conceptual strategies for characterizing interactions in microbial communities

Understanding the sets of inter- and intraspecies interactions in microbial communities is a fundamental goal of microbial ecology. However, the study and quantification of microbial interactions pose several challenges owing to their complexity, dynamic nature, and the sheer number of unique interactions within a typical community. To overcome such challenges, microbial ecologists must rely on various approaches to distill the system of study to a functional and conceptualizable level, allowing for a practical understanding of microbial interactions in both simplified and complex systems. This review broadly addresses the role of several conceptual approaches available for the microbial ecologist’s arsenal, examines specific tools used to accomplish such approaches, and describes how the assumptions, expectations, and philosophies underlying these tools change across scales of complexity.