Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Host-Associated and Free-Living Phage Communities Differ Profoundly in Phylogenetic Composition

  • J. Gregory Caporaso,

    Affiliation Department of Chemistry and Biochemistry, University of Colorado at Boulder, Boulder, Colorado, United States of America

  • Rob Knight,

    Affiliations Department of Chemistry and Biochemistry, University of Colorado at Boulder, Boulder, Colorado, United States of America, Howard Hughes Medical Institute, Boulder, Colorado, United States of America

  • Scott T. Kelley

    gregcaporaso@gmail.com

    Affiliation Department of Biology, San Diego State University, San Diego, California, United States of America

Abstract

Phylogenetic profiling has been widely used for comparing bacterial communities, but has so far been impossible to apply to viruses because of the lack of a single marker gene analogous to 16S rRNA. Here we developed a reference tree approach for matching viral sequences and applied it to the largest viral datasets available. The resulting technique, Shotgun UniFrac, was used to compare host-associated and non-host-associated phage communities (130 total metagenomes), and revealed a profound split similar to that found with bacterial communities. This new informatics approach complements analysis of bacterial communities and promises to provide new insights into viral community dynamics, such as top-down versus bottom-up control of bacterial communities by viruses in a range of systems.

Introduction

The phylogenetic composition of bacterial communities is primarily determined by whether they are found in host-associated or free-living environments [1]. Much less is known about the phylogenetic composition of viral communities, which may comprise most of the genetic diversity on Earth. If viral communities follow this pattern, microbial and viral community composition should be correlated, adding to recent evidence that phage predation can exert top-down control on microbial communities [2], [3].

The lack of a single marker gene in viral genomes complicates phylogenetic profiling of viral communities, a powerful technique for studying microbial communities, and previous studies have focused on profiling viral gene functions [4]. To complement these data with phylogenetic profiles of phage community composition, we developed Shotgun UniFrac (Figure 1). Shotgun UniFrac matches metagenomic reads against full phage genomes from the Phage Proteomic Tree [5] using BLAST. OTUs are assigned to reads by best hit, discarding reads with no significant hit, and UniFrac is applied using QIIME [6] and the Phage Proteomic Tree.

thumbnail
Figure 1. Schematic of the Shotgun UniFrac analysis pipeline.

https://doi.org/10.1371/journal.pone.0016900.g001

Results

We applied Shotgun UniFrac to 130 phage metagenomes from diverse environments. As observed with microbial communities, the primary factor separating metagenomes was whether they were derived from a free-living or host-associated environment. Host-associated environments vary more than a variety of free-living communities (considering only matches to the subset of viruses in the reference tree), and phage communities from the same host species tended to cluster (Figure 2a).

thumbnail
Figure 2. Principal Coordinates plot of weighted Shotgun UniFrac distances between viral communities where each point represents a metagenome colored by (a) host type and (b) data source.

https://doi.org/10.1371/journal.pone.0016900.g002

Our analysis also included 26 human feces phage metagenomes from 12 individuals with between 1 and 4 metagenomes per individual (recently presented in [7]). To include a metagenome in this analysis, we required a minimum of 200 reads assignable to a viral genome. We observed clustering of metagenomes by individual, although some aberrant clustering occurred (Figure 3a). This is likely due to the limited number of phage genomes currently available, which limits the resolution of Shotgun UniFrac (see Discussion). Confirming the observations of [7], [8] we found between-individual Shotgun UniFrac distances to be significantly greater than within-individual distances (Figure 3b; p = 3×10−23, one-tailed t-test; p<0.001, Monte Carlo t-test with 1000 iterations), suggesting stability in distal gut phage community membership over time.

thumbnail
Figure 3. (a) UPGMA clustering of individuals by weighted Shotgun UniFrac distances between metagenomes.

Cases where metagenomes from a single individual cluster monophyletically are highlighted in red. Cases where only a single metagenome for an individual was included are highlighted in blue. 1000 jackknife iterations were performed at a depth of 200 sequences per metagenome, and jackknife support values are provided for each node. The Reyes et al. analysis from which these samples were derived studied gut microbial communities from human twins and their mothers. The labels for each sample indicate the individual where: Fn corresponds to family number n; M corresponds to mother; and T1 and T2 refer to twin 1 and twin 2, respectively. (b) Histograms of within individual (grey) and between individual (pink) Shotgun UniFrac distances.

https://doi.org/10.1371/journal.pone.0016900.g003

Discussion

Taken together, our results suggest that phage communities mirror microbial communities, and that comparison of phage communities by phylogenetic identity of viral types, even with relatively few sequenced phage genomes available to assign sequences, can be a powerful complement to functional profiles of the communities. Collecting viral metagenomes, microbial metagenomes, and 16S reads from the same samples and comparing these data with techniques such as Procrustes analysis [9] will provide insight into fundamental parameters of microbial ecosystems, such as whether control occurs in a top-down or bottom-up manner.

Currently the limiting factor in applying Shotgun UniFrac to phage data is the availability of phage genomes, because sequences not matching known genomes are excluded from the analysis. For some metagenome types less than 1% of the viral metagenomic sequences could be classified (Table 1, Table S1) resulting in relatively few sequences per metagenome for comparing communities. The UniFrac results presented in Figures 2 and 3 are based on exactly 200 sequences per metagenome. Data sets of this size are useful for comparing microbial communities [10] and phage communities (Figure 2), but increasing the database of sequenced phage genomes and their phylogenies will further enhance the resolution of these techniques. Better resolution will aid understanding the complex dynamics and large compositional shifts seen in the human infant microbiome and virome [11], [12] that might be due to predator-prey cycling leading to chaos. Understanding such disruptions might be key to developing an understanding of probiotics and a wide range of time-variable diseases, such as Crohn's disease.

Materials and Methods

Viral community metagenomic data was compiled from CAMERA [13], MG-RAST [14], and study authors [7] (Table S2, Table S3). There was no community clustering by data source (Figure 2b). Sequences were assigned to source viral genomes using Shotgun UniFrac, an extension of the reference-based OTU picking strategy presented by [15], using the open source QIIME and PyCogent [16] toolkits. Shotgun UniFrac was applied against full phage genomes from the Phage Proteomic Tree, and the associated reference tree was used for phylogenetic beta diversity analysis. Sequences were assigned to a viral genome if they achieved an E-value of less than 0.001, resulting in the viral OTU table (Table S4). The viral OTU table was then sub-sampled to 200 sequences per metagenome (Table S5) to control for depth of coverage. The UniFrac diversity metric was applied to the sub-sampled viral OTU table using the Phage Proteomic Tree. The version of the Phage Proteomic Tree used here contains 651 tips built from fully sequenced phage genomes as described in [5]. Community clustering and within- versus between-individual Shotgun UniFrac distances were calculated using Weighted UniFrac. Shotgun UniFrac analysis, Principal Coordinates Analysis, distance calculations and plotting were all performed using QIIME, and Shotgun UniFrac is accessible in QIIME v1.2.0-dev using the pick_reference_otus_through_otu_table.py workflow.

The number of input metagenomes by type were: Reclaimed water at discharge point (n = 1); Reclaimed water at point-of-use (n = 2); Freshwater stromatolite (n = 2); Hot Spring, Yellowstone National Park (n = 2); Potable water (n = 1); Saltern (medium salinity) (n = 5); Ocean (db:MG-RAST) (n = 4); Saltern (high salinity) (n = 3); Northern Islands Coral (n = 4); Marine stromatolite (n = 1); Ocean (db:CAMERA) (n = 4); Freshwater (n = 4); Human feces (n = 80); Saltern (low salinity) (n = 3); Healthy human lung (n = 2); Mosquito-associated (n = 3); Cystic fibrosis human lung (n = 3); Porites compressa (coral, wild and experimentally treated) (n = 6). Four overlapping metagenomes (Ocean (db:MG-RAST) and Ocean (db:CAMERA)), were used as controls to ensure that the source database did not affect the clustering results which is possible, for example, if one required preprocessing that the other did not.

Supporting Information

Table S1.

OTU assignment statistics by metagenome.

https://doi.org/10.1371/journal.pone.0016900.s001

(XLS)

Table S2.

Description of metagenome types and sources.

https://doi.org/10.1371/journal.pone.0016900.s002

(XLS)

Table S3.

Full QIIME metadata mapping file.

https://doi.org/10.1371/journal.pone.0016900.s003

(XLS)

Table S4.

Full viral OTU table (i.e., metagenome × viral OTU abundance matrix). These data were used in jackknifed weighted Shotgun UniFrac calculations (Figure 3a).

https://doi.org/10.1371/journal.pone.0016900.s004

(XLS)

Table S5.

Viral OTU table sub-sampled to 200 sequences per metagenome. These data were used in weighted UniFrac calculations (Figure 2 and Figure 3b).

https://doi.org/10.1371/journal.pone.0016900.s005

(XLS)

Acknowledgments

We thank Robert A. Edwards for access to the latest version of the Phage Proteomic Tree and useful discussion about the manuscript; Jeffrey I. Gordon, Alejandro Reyes, Forrest Rohwer and Ruth Ley for comments on drafts; and Ingrid Caporaso for suggesting the name “Shotgun UniFrac.”

Author Contributions

Conceived and designed the experiments: STK. Performed the experiments: JGC. Analyzed the data: JGC. Contributed reagents/materials/analysis tools: JGC STK RK. Wrote the manuscript: JGC STK RK.

References

  1. 1. Ley RE, Lozupone CA, Hamady M, Knight R, Gordon JI (2008) Worlds within worlds: evolution of the vertebrate gut microbiota. Nat Rev Microbiol 6: 776–788.
  2. 2. Allen B, Willner D, Oechel WC, Lipson D (2009) Top-down control of microbial activity and biomass in an Arctic soil ecosystem. Environ Microbiol 12: 642–648.
  3. 3. Sandaa RA, Gomez-Consarnau L, Pinhassi J, Riemann L, Malits A, et al. (2009) Viral control of bacterial biodiversity–evidence from a nutrient-enriched marine mesocosm experiment. Environ Microbiol 11: 2585–2597.
  4. 4. Dinsdale EA, Edwards RA, Hall D, Angly F, Breitbart M, et al. (2008) Functional metagenomic profiling of nine biomes. Nature 452: 629–632.
  5. 5. Rohwer F, Edwards R (2002) The Phage Proteomic Tree: a genome-based taxonomy for phage. J Bacteriol 184: 4529–4535.
  6. 6. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, et al. (2010) QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7: 335–336.
  7. 7. Reyes A, Haynes M, Hanson N, Angly FE, Heath AC, et al. (2010) Viruses in the faecal microbiota of monozygotic twins and their mothers. Nature 466: 334–338.
  8. 8. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, et al. (2009) A core gut microbiome in obese and lean twins. Nature 457: 480–484.
  9. 9. Gower JC (1975) Generalized Procrustes Analysis. Psychometrika 40:
  10. 10. Kuczynski J, Liu Z, Lozupone C, McDonald D, Fierer N, et al. (2010) Microbial community resemblance methods differ in their ability to detect biologically relevant patterns. Nat Methods 7: 813–819.
  11. 11. Breitbart M, Haynes M, Kelley S, Angly F, Edwards RA, et al. (2008) Viral diversity and dynamics in an infant gut. Res Microbiol 159: 367–373.
  12. 12. Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO (2007) Development of the human infant intestinal microbiota. PLoS Biol 5: e177.
  13. 13. Seshadri R, Kravitz SA, Smarr L, Gilna P, Frazier M (2007) CAMERA: a community resource for metagenomics. PLoS Biol 5: e75.
  14. 14. Meyer F, Paarmann D, D'Souza M, Olson R, Glass EM, et al. (2008) The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9: 386.
  15. 15. Hamady M, Lozupone C, Knight R (2009) Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J.
  16. 16. Knight R, Maxwell P, Birmingham A, Carnes J, Caporaso JG, et al. (2007) PyCogent: a toolkit for making sense from sequence. Genome Biol 8: R171.