Comparative genomics of two 'Candidatus Accumulibacter' clades performing biological phosphorus removal.

Members of the genus Candidatus Accumulibacter are important in many wastewater treatment systems performing enhanced biological phosphorus removal (EBPR). The Accumulibacter lineage can be subdivided phylogenetically into multiple clades, and previous work showed that these clades are ecologically distinct. The complete genome of Candidatus Accumulibacter phosphatis strain UW-1, a member of Clade IIA, was previously sequenced. Here, we report a draft genome sequence of Candidatus Accumulibacter spp. strain UW-2, a member of Clade IA, assembled following shotgun metagenomic sequencing of laboratory-scale bioreactor sludge. We estimate the genome to be 80-90% complete. Although the two clades share 16S rRNA sequence identity of >98.0%, we observed a remarkable lack of synteny between the two genomes. We identified 2317 genes shared between the two genomes, with an average nucleotide identity (ANI) of 78.3%, and accounting for 49% of genes in the UW-1 genome. Unlike UW-1, the UW-2 genome seemed to lack genes for nitrogen fixation and carbon fixation. Despite these differences, metabolic genes essential for denitrification and EBPR, including carbon storage polymer and polyphosphate metabolism, were conserved in both genomes. The ANI from genes associated with EBPR was statistically higher than that from genes not associated with EBPR, indicating a high selective pressure in EBPR systems. Further, we identified genomic islands of foreign origins including a near-complete lysogenic phage in the Clade IA genome. Interestingly, Clade IA appeared to be more phage susceptible based on it containing only a single Clustered Regularly Interspaced Short Palindromic Repeats locus as compared with the two found in Clade IIA. Overall, the comparative analysis provided a genetic basis to understand physiological differences and ecological niches of Accumulibacter populations, and highlights the importance of diversity in maintaining system functional resilience.

'Candidatus Accumulibacter' gene expression in response to dynamic EBPR conditions.

Enhanced biological phosphorus removal (EBPR) activated sludge communities enriched in 'Candidatus Accumulibacter' relatives are widely used in wastewater treatment, but much remains to be learned about molecular-level controls on the EBPR process. The expression of genes found in the carbon and polyphosphate metabolic pathways in Accumulibacter was investigated using reverse transcription quantitative PCR. During a normal anaerobic/aerobic EBPR cycle, gene expression exhibited a dynamic change in response to external acetate, oxygen, phosphate concentrations and probably internal chemical pools. Anaerobic acetate addition induced expression of genes associated with the methylmalonyl-CoA pathway enabling the split mode of the tricarboxylic acid (TCA) cycle. Components of the full TCA cycle were induced after the switch to aerobic conditions. The induction of a key gene in the glyoxylate shunt pathway was observed under both anaerobic and aerobic conditions, with a higher induction by aeration. Polyphosphate kinase 1 from Accumulibacter was expressed, but did not appear to be regulated by phosphate limitation. To understand how Accumulibacter responds to disturbed electron donor and acceptor conditions, we perturbed the process by adding acetate aerobically. When high concentrations of oxygen were present simultaneously with acetate, phosphate-release was almost completely inhibited, and polyphosphate kinase 1 transcript abundance decreased. Genes associated with the methylmalonyl-CoA pathway were repressed and genes associated with the aerobic TCA cycle exhibited higher expression under this perturbation, suggesting that more acetyl-CoA was metabolized through the TCA cycle. These findings suggest that several genes involved in EBPR are tightly regulated at the transcriptional level.

Bacterial community and "Candidatus Accumulibacter" population dynamics in laboratory-scale enhanced biological phosphorus removal reactors.

"Candidatus Accumulibacter" and total bacterial community dynamics were studied in two lab-scale enhanced biological phosphorus removal (EBPR) reactors by using a community fingerprint technique, automated ribosomal intergenic spacer analysis (ARISA). We first evaluated the quantitative capability of ARISA compared to quantitative real-time PCR (qPCR). ARISA and qPCR provided comparable relative quantification of the two dominant "Ca. Accumulibacter" clades (IA and IIA) detected in our reactors. The quantification of total "Ca. Accumulibacter" 16S rRNA genes relative to that from the total bacterial community was highly correlated, with ARISA systematically underestimating "Ca. Accumulibacter" abundance, probably due to the different normalization techniques applied. During 6 months of normal (undisturbed) operation, the distribution of the two clades within the total "Ca. Accumulibacter" population was quite stable in one reactor while comparatively dynamic in the other reactor. However, the variance in the clade distribution did not appear to affect reactor performance. Instead, good EBPR activity was positively associated with the abundance of total "Ca. Accumulibacter." Therefore, we concluded that the different clades in the system provided functional redundancy. We disturbed the reactor operation by adding nitrate together with acetate feeding in the anaerobic phase to reach initial reactor concentrations of 10 mg/liter NO(3)-N for 35 days. The reactor performance deteriorated with a concomitant decrease in the total "Ca. Accumulibacter" population, suggesting that a population shift was the cause of performance upset after a long exposure to nitrate in the anaerobic phase.

Metatranscriptomic array analysis of 'Candidatus Accumulibacter phosphatis'-enriched enhanced biological phosphorus removal sludge.

Here we report the first metatranscriptomic analysis of gene expression and regulation of 'Candidatus Accumulibacter'-enriched lab-scale sludge during enhanced biological phosphorus removal (EBPR). Medium density oligonucleotide microarrays were generated with probes targeting most predicted genes hypothesized to be important for the EBPR phenotype. RNA samples were collected at the early stage of anaerobic and aerobic phases (15 min after acetate addition and switching to aeration respectively). We detected the expression of a number of genes involved in the carbon and phosphate metabolisms, as proposed by EBPR models (e.g. polyhydroxyalkanoate synthesis, a split TCA cycle through methylmalonyl-CoA pathway, and polyphosphate formation), as well as novel genes discovered through metagenomic analysis. The comparison between the early stage anaerobic and aerobic gene expression profiles showed that expression levels of most genes were not significantly different between the two stages. The majority of upregulated genes in the aerobic sample are predicted to encode functions such as transcription, translation and protein translocation, reflecting the rapid growth phase of Accumulibacter shortly after being switched to aerobic conditions. Components of the TCA cycle and machinery involved in ATP synthesis were also upregulated during the early aerobic phase. These findings support the predictions of EBPR metabolic models that the oxidation of intracellularly stored carbon polymers through the TCA cycle provides ATP for cell growth when oxygen becomes available. Nitrous oxide reductase was among the very few Accumulibacter genes upregulated in the anaerobic sample, suggesting that its expression is likely induced by the deprivation of oxygen.

Progress toward understanding the distribution of Accumulibacter among full-scale enhanced biological phosphorus removal systems.

This study investigated the role of Accumulibacter-related bacterial populations and factors influencing their distribution in enhanced biological phosphorus removal (EBPR) systems in the USA. For this purpose, five full-scale wastewater treatment facilities performing EBPR were surveyed. The facilities had different configurations but were all treating primarily domestic wastewater. Two facilities had history of poor EBPR performance. Batch-scale acetate uptake and inorganic phosphate (P(i)) release and uptake experiments were conducted to evaluate the EBPR activity of each sludge. Typical P(i) and acetate profiles were observed, and EBPR activity was found to be positively correlated to polyphosphate (polyP)-accumulating organism (PAO) abundance, as determined by staining intracellular polyP. The abundance of Accumulibacter-related organisms was investigated using fluorescent in situ hybridization. Accumulibacter-related organisms were present in all full-scale EBPR facilities, at levels ranging from 9 to 24% of total cells. More than 80% of Accumulibacter-related organisms were estimated to have high polyP content, confirming their involvement in EBPR in these five facilities. However, Accumulibacter-related PAOs were only a fraction (40-69%) of the total PAO population. The variation of Accumulibacter-related PAO abundance among these EBPR systems suggests that multiple interacting factors such as wastewater characteristics and operational conditions are structuring PAO communities.

"Candidatus Accumulibacter" population structure in enhanced biological phosphorus removal sludges as revealed by polyphosphate kinase genes.

We investigated the fine-scale population structure of the "Candidatus Accumulibacter" lineage in enhanced biological phosphorus removal (EBPR) systems using the polyphosphate kinase 1 gene (ppk1) as a genetic marker. We retrieved fragments of "Candidatus Accumulibacter" 16S rRNA and ppk1 genes from one laboratory-scale and several full-scale EBPR systems. Phylogenies reconstructed using 16S rRNA genes and ppk1 were largely congruent, with ppk1 granting higher phylogenetic resolution and clearer tree topology and thus serving as a better genetic marker than 16S rRNA for revealing population structure within the "Candidatus Accumulibacter" lineage. Sequences from at least five clades of "Candidatus Accumulibacter" were recovered by ppk1-targeted PCR, and subsequently, specific primer sets were designed to target the ppk1 gene for each clade. Quantitative real-time PCR (qPCR) assays using "Candidatus Accumulibacter"-specific 16S rRNA and "Candidatus Accumulibacter" clade-specific ppk1 primers were developed and conducted on three laboratory-scale and nine full-scale EBPR samples and two full-scale non-EBPR samples to determine the abundance of the total "Candidatus Accumulibacter" lineage and the relative distributions and abundances of the five "Candidatus Accumulibacter" clades. The qPCR-based estimation of the total "Candidatus Accumulibacter" fraction as a proportion of the bacterial community as measured using 16S rRNA genes was not significantly different from the estimation measured using ppk1, demonstrating the power of ppk1 as a genetic marker for detection of all currently defined "Candidatus Accumulibacter" clades. The relative distributions of "Candidatus Accumulibacter" clades varied among different EBPR systems and also temporally within a system. Our results suggest that the "Candidatus Accumulibacter" lineage is more diverse than previously realized and that different clades within the lineage are ecologically distinct.

Metagenomic analysis of two enhanced biological phosphorus removal (EBPR) sludge communities.

Enhanced biological phosphorus removal (EBPR) is one of the best-studied microbially mediated industrial processes because of its ecological and economic relevance. Despite this, it is not well understood at the metabolic level. Here we present a metagenomic analysis of two lab-scale EBPR sludges dominated by the uncultured bacterium, "Candidatus Accumulibacter phosphatis." The analysis sheds light on several controversies in EBPR metabolic models and provides hypotheses explaining the dominance of A. phosphatis in this habitat, its lifestyle outside EBPR and probable cultivation requirements. Comparison of the same species from different EBPR sludges highlights recent evolutionary dynamics in the A. phosphatis genome that could be linked to mechanisms for environmental adaptation. In spite of an apparent lack of phylogenetic overlap in the flanking communities of the two sludges studied, common functional themes were found, at least one of them complementary to the inferred metabolism of the dominant organism. The present study provides a much needed blueprint for a systems-level understanding of EBPR and illustrates that metagenomics enables detailed, often novel, insights into even well-studied biological systems.