|
|
|
|
Published online before print
December 12, 2007, 10.1101/gr.6835308 Genome Res. 18:293-297, 2008 ©2008 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/08 $5.00 OPEN ACCESS ARTICLE
Letter A bacterial metapopulation adapts locally to phage predation despite global dispersal1 Microbial Ecology Program, Department of Energy Joint Genome Institute, Walnut Creek, California 94598, USA; 2 Department of Civil and Environmental Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA; 3 Department of Plant Pathology, University of Wisconsin-Madison, Madison, Wisconsin 53706 USA; 4 Department of Biology, San Diego State University, San Diego, California 92182, USA; 5 Genome Biology Program, Department of Energy Joint Genome Institute, Walnut Creek, California 94598, USA; 6 Advanced Wastewater Management Centre, University of Queensland, St Lucia 4072, Queensland, Australia; 7 University of South Florida, St. Petersburg, Florida 33701 USA
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.
Ecological theory is largely grounded on the study of macroscopic communities (Begon et al. 2006
Molecular methods developed over the past decade are addressing these limitations and allowing microbial ecology to mature as a discipline and, in the process, are challenging long-held assumptions about microbial populations. For example, multilocus sequence typing (MLST) (Maiden et al. 1998
Metagenomics, the application of shotgun sequencing to environmental samples, holds the promise of providing the least biased (culture-independent) and most comprehensive (genome-wide) resolution of sympatric populations (Whitaker and Banfield 2006
We began by searching for evidence of geographic isolation of the CAP populations by analyzing the phylogenetic distribution of 48 single-copy genes (Supplemental Table S1). Using single-copy genes ensures that any given strain is not represented by more than one sequence and therefore minimizes the possibility of misinterpreting paralogs as orthologs (Venter et al. 2004
Contrary to recent findings in hot springs using high-resolution molecular methods (Papke et al. 2003
To date, CAP has only been detected in activated sludges (Hesselmann et al. 1999
To identify CAP habitats, we surveyed a range of environmental samples using Accumulibacter-specific PCR targeting the 16S rRNA and ppk genes that were subsequently confirmed by sequencing. Accumulibacter species were detected in both fresh and estuarine waters and associated sediments but were rarely observed in soil samples (Supplemental Table S2). We therefore suggest that CAP populations are distributed in the environment as sparse high-density point sources (EBPR sludges) linked by dispersal via widespread diffuse reservoirs (aquatic environments), conforming to the ecological definition of a metapopulation as a collection of contained populations connected by a small amount of gene flow (Hanski 1999
The presence of multiple strains in each sludge sample allowed us to investigate CAP for evidence of homologous recombination between strains (Supplemental Fig. S1). Unlike recent studies in which microbial populations were found to be highly recombining (Tyson et al. 2004 While most of the CAP strains were represented by unassembled reads or short contigs in the metagenomic data (indicating low abundance), one strain dominated each sludge, producing large contigs with high read depths, allowing assessment of within-strain heterogeneity. The dominant strain populations were found to be extremely homogeneous in the US and OZ sludges with an average of one confirmed single nucleotide polymorphism (SNP) per 163.2 and 65.6 kb, respectively (Supplemental Table S3). This indicates that both dominant CAP strains are virtually clonal.
The near clonality of the dominant strains, and their inability to recombine, means that the bulk of the biomass in each laboratory-scale EBPR sludge is composed of genetically identical cells. Such populations are natural targets for phage predation, via the so-called "kill-the-winner" phenomenon (Thingstad and Lignell 1997
Another phage defense mechanism are clustered regularly interspaced short palindromic repeat (CRISPR) elements (Jansen et al. 2002
CRISPR elements and EPS gene clusters were among the most notable differences between closely related strains of Streptococcus thermophilus, which is used in coculture with Lactobacillus species for industrial yogurt and cheese production (Bolotin et al. 2004 To test the hypothesis that phages are playing a major role in structuring CAP populations in EBPR, we sampled the phage virion metagenome of the US sludge 7 mo after sampling the bacterial metagenome. Eleven US CRISPR spacers, eight of which belonged to the dominant CAP strain, had matches to phage genome fragments (Supplemental Table S5), with some phages being targeted by multiple spacers (Fig. 2C) and some spacers targeting multiple related phages. This provides a direct link between the uncultivated bacterial host and phage virions and confirms that the CAP population had previously been infected by these phages. Two CRISPR spacers found only in the dominant OZ CAP strain had matches to the US phage community, supporting geographic dispersal of the host and/or phage. To confirm that phages are active in the sludge ecosystem, we monitored the US sludge at three time points spanning 3 mo using expression arrays targeting both phage and bacterial genes obtained from the metagenomic data sets. We found that large numbers of genes originating from the phage virion metagenome and some genes in the bacterial metagenome of putative prophage origin were highly expressed (Table 1 and Supplemental Table S6). These included many hypothetical proteins but also proteins associated with phage tail assembly, a phage-specific endonuclease and terminase (Supplemental Table S6), suggesting that phages are continuously active in the sludge. Since the microarray was based on phage virion genes sampled almost 2 yr prior to the expression analysis, some phages must persist for long periods in the sludge. These data imply that the bacterial community is under persistent local predation pressure by phages and live in a volatile but relatively stable coexistence.
In summary, we have shown that (1) CAP is globally dispersed, (2) highly related and geographically remote CAP strains differ principally in genomic regions encoding phage defense mechanisms, (3) high-density, clonal, nonrecombining CAP populations in EBPR bioreactors are natural targets of "kill-the-winner" phage predation, (4) phages are consistently active in EBPR bioreactor communities, and (5) signatures of past phage infections in CAP are observed mostly between local phage and host. We therefore conclude that CAP strains disperse globally but must adapt to local persistent phage predation pressure. The present study illustrates the value of combining high-throughput sequence and gene expression data from the bacterial and viral fractions of an ecosystem to elucidate population structure, biogeography, and host–parasite interdependencies.
Metagenomic sequencing Sludge samples for the US and Australian (OZ) bacterial metagenomes were obtained on July 3 and August 18, 2004, respectively. Sequencing, assembly, and gene prediction of the bacterial metagenomes are described elsewhere (Garcia Martin et al. 2006
Bioinformatic analyses To refine the resolution of the single-copy gene analysis, we PCR-amplified the ppk gene from the sludge biomass and environmental samples. The amplification product was cloned into Escherichia coli, and 96 clones were picked from each library for sequencing. See Supplemental Research Data for further details.
Quantification of SNP frequency was done using CONSED program (Gordon et al. 1998
Microarrays
This work was performed under the auspices of the U.S. Department of Energys Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48, Lawrence Berkeley National Laboratory under contract No. DE-AC02-05CH11231, and Los Alamos National Laboratory under contract No. DE-AC02-06NA25396.
8 Corresponding author.
E-mail phugenholtz{at}lbl.gov; fax (925) 296-5720. [Supplemental material is available online at www.genome.org.] Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.6835308
Altschul, S.F., Madden, T.L., Schaffer, A.A., Zhang, J., Zhang, Z., Miller, W., and Lipman, D.J. 1997. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 25: 3389–3402. doi: 10.1093/nar/25.17.3389. Barrangou, R., Fremaux, C., Deveau, H., Richards, M., Boyaval, P., Moineau, S., Romero, D.A., and Horvath, P. 2007. CRISPR provides acquired resistance against viruses in prokaryotes. Science 315: 1709–1712. Bateman, A., Coin, L., Durbin, R., Finn, R.D., Hollich, V., Griffiths-Jones, S., Khanna, A., Marshall, M., Moxon, S., Sonnhammer, E.L., et al. 2004. The Pfam protein families database. Nucleic Acids Res. 32: D138–D141. doi: 10.1093/nar/gkh121. Begon, M., Townsend, C.R., and Harper, J.L. 2006. Ecology: From individuals to ecosystems. Blackwell, Malden, MA. Bolotin, A., Quinquis, B., Renault, P., Sorokin, A., Ehrlich, S.D., Kulakauskas, S., Lapidus, A., Goltsman, E., Mazur, M., Pusch, G.D., et al. 2004. Complete sequence and comparative genome analysis of the dairy bacterium Streptococcus thermophilus. Nat. Biotechnol. 22: 1554–1558.[CrossRef][Medline] Buckling, A. and Rainey, P.B. 2002. The role of parasites in sympatric and allopatric host diversification. Nature 420: 496–499.[CrossRef][Medline] Crocetti, G.R., Hugenholtz, P., Bond, P.L., Schuler, A., Keller, J., Jenkins, D., and Blackall, L.L. 2000. Identification of polyphosphate-accumulating organisms and design of 16S rRNA-directed probes for their detection and quantitation. Appl. Environ. Microbiol. 66: 1175–1182. Edgar, R.C. 2007. PILER-CR: Fast and accurate identification of CRISPR repeats. BMC Bioinformatics 8: 18. doi: 10.1186/1471-2105-8-18.[CrossRef][Medline] Feil, E.J., Enright, M.C., and Spratt, B.G. 2000. Estimating the relative contributions of mutation and recombination to clonal diversification: A comparison between Neisseria meningitidis and Streptococcus pneumoniae. Res. Microbiol. 151: 465–469.[Medline] Garcia Martin, H., Ivanova, N., Kunin, V., Warnecke, F., Barry, K.W., McHardy, A.C., Yeates, C., He, S., Salamov, A.A., Szeto, E., et al. 2006. Metagenomic analysis of two enhanced biological phosphorus removal (EBPR) sludge communities. Nat. Biotechnol. 24: 1263–1269.[CrossRef][Medline] Gordon, D., Abajian, C., and Green, P. 1998. Consed: A graphical tool for sequence finishing. Genome Res. 8: 195–202. Hanski, I. 1999. Metapopulation ecology. Oxford University Press, Oxford, UK. Hesselmann, R.P., Werlen, C., Hahn, D., van der Meer, J.R., and Zehnder, A.J. 1999. Enrichment, phylogenetic analysis and detection of a bacterium that performs enhanced biological phosphate removal in activated sludge. Syst. Appl. Microbiol. 22: 454–465.[Medline] Jansen, R., Embden, J.D., Gaastra, W., and Schouls, L.M. 2002. Identification of genes that are associated with DNA repeats in prokaryotes. Mol. Microbiol. 43: 1565–1575.[CrossRef][Medline] Jessup, C.M., Kassen, R., Forde, S.E., Kerr, B., Buckling, A., Rainey, P.B., and Bohannan, B.J.M. 2004. Big questions, small worlds: Microbial model systems in ecology. Trends Ecol. Evol. 19: 189–197.[CrossRef][Medline] Maiden, M.C., Bygraves, J.A., Feil, E., Morelli, G., Russell, J.E., Urwin, R., Zhang, Q., Zhou, J., Zurth, K., Caugant, D.A., et al. 1998. Multilocus sequence typing: A portable approach to the identification of clones within populations of pathogenic microorganisms. Proc. Natl. Acad. Sci. 95: 3140–3145. Markowitz, V.M., Korzeniewski, F., Palaniappan, K., Szeto, E., Werner, G., Padki, A., Zhao, X., Dubchak, I., Hugenholtz, P., Anderson, I., et al. 2006. The integrated microbial genomes (IMG) system. Nucleic Acids Res. 34: D344–D348. doi: 10.1093/nar/gkj024. McMahon, K.D., Dojka, M.A., Pace, N.R., Jenkins, D., and Keasling, J.D. 2002. Polyphosphate kinase from activated sludge performing enhanced biological phosphorus removal. Appl. Environ. Microbiol. 68: 4971–4978. Nesbø, C.L., Dlutek, M., and Doolittle, W.F. 2006. Recombination in Thermotoga: Implications for species concepts and biogeography. Genetics 172: 759–769. Pace, N.R. 1997. A molecular view of microbial diversity and the biosphere. Science 276: 734–740. Papke, R.T., Ramsing, N.B., Bateson, M.M., and Ward, D.M. 2003. Geographical isolation in hot spring cyanobacteria. Environ. Microbiol. 5: 650–659.[CrossRef][Medline] Pernthaler, J. 2005. Predation on prokaryotes in the water column and its ecological implications. Nat. Rev. Microbiol. 3: 537–546.[CrossRef][Medline] Shah, N., Teplitsky, M.V., Minovitsky, S., Pennacchio, L.A., Hugenholtz, P., Hamann, B., and Dubchak, I. 2005. SNP-VISTA: An interactive SNP visualization tool. BMC Bioinformatics 6: 292. doi: 10.1186/1471-2105-6-292.[CrossRef][Medline] Sutherland, I. 2001. Biofilm exopolysaccharides: A strong and sticky framework. Microbiol. 147: 3–9. Tchobanoglous, G., Burton, F.L., and Stensel, H.D. 2003. Wastewater engineering: Treatment and reuse. Metcalf & Eddy–McGraw-Hill, Boston. Thingstad, T. and Lignell, R. 1997. Theoretical models for the control of bacterial growth rate, abundance, diversity and carbon demand. Aquat. Microb. Ecol. 13: 19–27.[CrossRef] Thompson, J.D., Higgins, D.G., and Gibson, T.J. 1994. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22: 4673–4680. doi: 10.1093/nar/22.22.4673. Tyson, G.W., Chapman, J., Hugenholtz, P., Allen, E.E., Ram, R.J., Richardson, P.M., Solovyev, V.V., Rubin, E.M., Rokhsar, D.S., and Banfield, J.F. 2004. Community structure and metabolism through reconstruction of microbial genomes from the environment. Nature 428: 37–43.[CrossRef][Medline] Venter, J.C., Remington, K., Heidelberg, J.F., Halpern, A.L., Rusch, D., Eisen, J.A., Wu, D., Paulsen, I., Nelson, K.E., Nelson, W., et al. 2004. Environmental genome shotgun sequencing of the Sargasso Sea. Science 304: 66–74. Whitaker, R.J. and Banfield, J.F. 2006. Population genomics in natural microbial communities. Trends Ecol. Evol. 21: 508–516.[CrossRef][Medline] Whitaker, R.J., Grogan, D.W., and Taylor, J.W. 2003. Geographic barriers isolate endemic populations of hyperthermophilic archaea. Science 301: 976–978. Wong, M.T., Mino, T., Seviour, R.J., Onuki, M., and Liu, W.T. 2005. In situ identification and characterization of the microbial community structure of full-scale enhanced biological phosphorous removal plants in Japan. Water Res. 39: 2901–2914.[Medline] Zilles, J.L., Peccia, J., Kim, M.W., Hung, C.H., and Noguera, D.R. 2002. Involvement of Rhodocyclus-related organisms in phosphorus removal in full-scale wastewater treatment plants. Appl. Environ. Microbiol. 68: 2763–2769.
Received June 22, 2007; accepted in revised format October 14, 2007.
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||