|Home | About | Journals | Submit | Contact Us | Français|
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.
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.
The phylogenetic composition of bacterial communities is primarily determined by whether they are found in host-associated or free-living environments . 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 , .
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 . 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  using BLAST. OTUs are assigned to reads by best hit, discarding reads with no significant hit, and UniFrac is applied using QIIME  and the Phage Proteomic Tree.
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).
Our analysis also included 26 human feces phage metagenomes from 12 individuals with between 1 and 4 metagenomes per individual (recently presented in ). 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 ,  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.
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  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 and33 are based on exactly 200 sequences per metagenome. Data sets of this size are useful for comparing microbial communities  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 ,  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.
Viral community metagenomic data was compiled from CAMERA , MG-RAST , and study authors  (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 , using the open source QIIME and PyCogent  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 . 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.
OTU assignment statistics by metagenome.
Description of metagenome types and sources.
Full QIIME metadata mapping file.
Full viral OTU table (i.e., metagenome × viral OTU abundance matrix). These data were used in jackknifed weighted Shotgun UniFrac calculations (Figure 3a).
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.”
Competing Interests: The authors have declared that no competing interests exist.
Funding: This work was funded with grants from Crohn's and Colitis Foundation of America (http://www.ccfa.org), National Institutes of Health grant number HG004872 (http://www.nih.gov/), the Bill and Melinda Gates Foundation (http://www.gatesfoundation.org/), and the Howard Hughes Medical Institute (http://www.hhmi.org/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.