Modeling ecological drivers in marine viral communities using comparative metagenomics and network analyses

B. Hurwitz,A. Westveld,Jennifer R. Brum,M. Sullivan

DOI: https://doi.org/10.1073/pnas.1319778111

2014-07-07

Abstract:Significance Microorganisms and their viruses are increasingly recognized as drivers of myriad ecosystem processes. However, our knowledge of their roles is limited by the inability of culture-dependent and culture-independent (e.g., metagenomics) methods to be fully implemented at scales relevant to the diversity found in nature. Here we combine advances in bioinformatics (shared k-mer analyses) and social networking (regression modeling) to develop an annotation- and assembly-free visualization and analytical strategy for comparative metagenomics that uses all the data in a unified statistical framework. Application to 32 Pacific Ocean viromes, the first large-scale quantitative viral metagenomic dataset, tested existing and generated further hypotheses about ecological drivers of viral community structure. Highly computationally scalable, this new approach enables diverse sequence-based large-scale comparative studies. Long-standing questions in marine viral ecology are centered on understanding how viral assemblages change along gradients in space and time. However, investigating these fundamental ecological questions has been challenging due to incomplete representation of naturally occurring viral diversity in single gene- or morphology-based studies and an inability to identify up to 90% of reads in viral metagenomes (viromes). Although protein clustering techniques provide a significant advance by helping organize this unknown metagenomic sequence space, they typically use only ∼75% of the data and rely on assembly methods not yet tuned for naturally occurring sequence variation. Here, we introduce an annotation- and assembly-free strategy for comparative metagenomics that combines shared k-mer and social network analyses (regression modeling). This robust statistical framework enables visualization of complex sample networks and determination of ecological factors driving community structure. Application to 32 viromes from the Pacific Ocean Virome dataset identified clusters of samples broadly delineated by photic zone and revealed that geographic region, depth, and proximity to shore were significant predictors of community structure. Within subsets of this dataset, depth, season, and oxygen concentration were significant drivers of viral community structure at a single open ocean station, whereas variability along onshore–offshore transects was driven by oxygen concentration in an area with an oxygen minimum zone and not depth or proximity to shore, as might be expected. Together these results demonstrate that this highly scalable approach using complete metagenomic network-based comparisons can both test and generate hypotheses for ecological investigation of viral and microbial communities in nature.