Enhanced Biological Phosphorus Removal (EBPR) is not well understood at the metabolic level despite being one of the best-studied microbially-mediated industrial processes due to its ecological and economic relevance. Here we present a metagenomic analysis of two lab-scale EBPR sludges dominated by the uncultured bacterium, "Candidatus Accumulibacter phosphatis". This 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. 3Excessive inorganic phosphate (Pi) supply to freshwater negatively affects water quality and ecosystem balance through a process known as eutrophication 1 . Limitations on allowable Pi discharges from municipal and industrial sources via wastewater treatment have proven effective in reducing Pi levels in many waterways 2 . Increasingly stringent Pi limits for effluent wastewater are expected in the future and hence efficient and reliable Pi removal methods are required. Due to the massive quantity of wastewater treated daily (more than 120 billion liters in the US alone 3 ), any improvement in existing methods should have tangible economic and ecological consequences.Enhanced Biological Phosphorus Removal (EBPR) is a treatment process in which microorganisms remove Pi from wastewater by accumulating it inside their cells as polyphosphate. These polyphosphate-accumulating organisms (PAOs) are then allowed to settle in a separate tank (clarifier), leaving the effluent water largely Pi-depleted. EBPR is more economical in the long term 2 and has a lower environmental impact 4 than traditional (chemical) Pi removal 5 , but is prone to unpredictable failures due to loss or reduced activity of microbial populations responsible for Pi removal 6 . This is primarily because the design process is highly empirical due to an incomplete understanding of sludge microbial ecology. Environmental engineers and microbiologists have been studying EBPR since its introduction in municipal wastewater treatment plants over thirty years ago 5 with the goal of making it a more reliable industrial process. Typically, EBPR is studied in lab-scale sequencing batch reactors (SBRs) where the microbial community can be better monitored and perturbed, and PAOs can be enriched to much higher levels than in full scale systems 7 .For th...

Sequencing of 16S rRNA gene tags is a popular method for profiling and comparing microbial communities. The protocols and methods used, however, vary considerably with regard to amplification primers, sequencing primers, sequencing technologies; as well as quality filtering and clustering. How results are affected by these choices, and whether data produced with different protocols can be meaningfully compared, is often unknown. Here we compare results obtained using three different amplification primer sets (targeting V4, V6–V8, and V7–V8) and two sequencing technologies (454 pyrosequencing and Illumina MiSeq) using DNA from a mock community containing a known number of species as well as complex environmental samples whose PCR-independent profiles were estimated using shotgun sequencing. We find that paired-end MiSeq reads produce higher quality data and enabled the use of more aggressive quality control parameters over 454, resulting in a higher retention rate of high quality reads for downstream data analysis. While primer choice considerably influences quantitative abundance estimations, sequencing platform has relatively minor effects when matched primers are used. Beta diversity metrics are surprisingly robust to both primer and sequencing platform biases.

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 laboratoryscale 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.

Termites effectively feed on many types of lignocellulose assisted by their gut microbial symbionts. To better understand the microbial decomposition of biomass with varied chemical profiles, it is important to determine whether termites harbor different microbial symbionts with specialized functionalities geared toward different feeding regimens. In this study, we compared the microbiota in the hindgut paunch of Amitermes wheeleri collected from cow dung and Nasutitermes corniger feeding on sound wood by 16S rRNA pyrotag, comparative metagenomic and metatranscriptomic analyses. We found that Firmicutes and Spirochaetes were the most abundant phyla in A. wheeleri, in contrast to N. corniger where Spirochaetes and Fibrobacteres dominated. Despite this community divergence, a convergence was observed for functions essential to termite biology including hydrolytic enzymes, homoacetogenesis and cell motility and chemotaxis. Overrepresented functions in A. wheeleri relative to N. corniger microbiota included hemicellulose breakdown and fixed-nitrogen utilization. By contrast, glycoside hydrolases attacking celluloses and nitrogen fixation genes were overrepresented in N. corniger microbiota. These observations are consistent with dietary differences in carbohydrate composition and nutrient contents, but may also reflect the phylogenetic difference between the hosts.

The predominance of rRNAs in the transcriptome is a major technical challenge in sequence-based analysis of cDNAs from microbial isolates and communities. Several approaches have been applied to deplete rRNAs from (meta)transcriptomes, but no systematic investigation of potential biases introduced by any of these approaches has been reported. Here we validated the effectiveness and fidelity of the two most commonly used approaches, subtractive hybridization and exonuclease digestion as well as combinations of these treatments, on two synthetic five-microorganism metatranscriptomes using massively parallel sequencing. We found that the effectiveness of rRNA removal was a function of community composition and RNA integrity for these treatments. Subtractive hybridization alone introduced the least bias in relative transcript abundance, whereas exonuclease and in particular combined treatments greatly compromised mRNA abundance fidelity. Illumina sequencing itself also can compromise quantitative data analysis by introducing a G+C bias between runs.Rapid technological advances in ultra-high-throughput sequencing are making de novo sequencing of transcriptomes (RNA-seq) a viable alternative to microarray analysis of microbial isolates and communities 1 . A major technical challenge for de novo transcriptome sequencing is the low relative abundance of mRNAs in total cellular RNA (1-5%; ref.2), the bulk of which is rRNAs and tRNAs 3 . Unlike eukaryotic mRNAs, which can be selectively synthesized into cDNA by virtue of their poly(A) tails 4 , bacterial and archaeal cDNAs are predominantly rRNA -2-sequences 5,6 . Therefore, prokaryotic rRNAs are often removed before sequencing to improve mRNA detection sensitivity. Different methods have been used to eliminate prokaryotic rRNA, including subtractive hybridization with rRNA-specific probes 7,8 , digestion with exonuclease that preferentially acts on rRNA, poly(A) tail addition to discriminate against rRNA 9,10 , reverse transcription with rRNA-specific primers followed by RNase H digestion to degrade rRNA:DNA hybrids 11 , and gel electrophoresis size separation and extraction of non-rRNA bands 12 .Among these methods, subtractive hybridization and exonuclease digestion have become the most popular owing to the availability of commercial kits from Ambion (MICROBExpress Bacterial mRNA Enrichment kit) and Epicentre (mRNA-ONLY Prokaryotic mRNA Isolation kit). The former kit uses a subtractive hybridization with capture oligonucleotides specific to 16S and 23S rRNAs. It has been applied to both bacterial isolates and environmental samples, in one or two rounds 6,[13][14][15][16][17][18] . The Epicentre kit uses exonuclease to preferentially degrade processed RNAs with 5′ monophosphate (the majority of which are believed to be rRNAs) 19,20 . In some instances, these methods have been used in combination to improve rRNA removal [21][22][23] . There is no consensus, however, on the best approach, and existing data are anecdotal. Here we validated the performance of these kits wit...

Summary The capability of ‘Candidatus Accumulibacter’ to use nitrate as an electron acceptor for phosphorus uptake was investigated using two activated sludge communities. The two communities were enriched in Accumulibacter clade IA and clade IIA, respectively. By performing a series of batch experiments, we found that clade IA was able to couple nitrate reduction with phosphorus uptake, but clade IIA could not. These results agree with a previously proposed hypothesis that different populations of Accumulibacter have different nitrate reduction capabilities, and they will help to understand the ecological roles played by these two clades.

Using a combination of bacterial and phage-targeted metagenomics, we analyzed two geographically remote sludge bioreactors enriched in a single bacterial species Candidatus Accumulibacter phosphatis (CAP). We inferred unrestricted global movement of this species and identified aquatic ecosystems as the primary environmental reservoirs facilitating dispersal. Highly related and geographically remote CAP strains differed principally in genomic regions encoding phage defense mechanisms. We found that CAP populations were high density, clonal, and nonrecombining, providing natural targets for “kill-the-winner” phage predation. Community expression analysis demonstrated that phages were consistently active in the bioreactor community. Genomic signatures linking CAP to past phage exposures were observed mostly between local phage and host. We conclude that CAP strains disperse globally but must adapt to phage predation pressure locally.

Extracellular electron transfer (EET) is recognized as a key biochemical process in circumneutral pH Fe(II)-oxidizing bacteria (FeOB). In this study, we searched for candidate EET genes in 73 neutrophilic FeOB genomes, among which 43 genomes are complete or close-to-complete and the rest have estimated genome completeness ranging from 5 to 91%. These neutrophilic FeOB span members of the microaerophilic, anaerobic phototrophic, and anaerobic nitrate-reducing FeOB groups. We found that many microaerophilic and several anaerobic FeOB possess homologs of Cyc2, an outer membrane cytochrome c originally identified in Acidithiobacillus ferrooxidans. The “porin-cytochrome c complex” (PCC) gene clusters homologous to MtoAB/PioAB are present in eight FeOB, accounting for 19% of complete and close-to-complete genomes examined, whereas PCC genes homologous to OmbB-OmaB-OmcB in Geobacter sulfurreducens are absent. Further, we discovered gene clusters that may potentially encode two novel PCC types. First, a cluster (tentatively named “PCC3”) encodes a porin, an extracellular and a periplasmic cytochrome c with remarkably large numbers of heme-binding motifs. Second, a cluster (tentatively named “PCC4”) encodes a porin and three periplasmic multiheme cytochromes c. A conserved inner membrane protein (IMP) encoded in PCC3 and PCC4 gene clusters might be responsible for translocating electrons across the inner membrane. Other bacteria possessing PCC3 and PCC4 are mostly Proteobacteria isolated from environments with a potential niche for Fe(II) oxidation. In addition to cytochrome c, multicopper oxidase (MCO) genes potentially involved in Fe(II) oxidation were also identified. Notably, candidate EET genes were not found in some FeOB, especially the anaerobic ones, probably suggesting EET genes or Fe(II) oxidation mechanisms are different from the searched models. Overall, based on current EET models, the search extends our understanding of bacterial EET and provides candidate genes for future research.