Abstract Background Marine ecosystem function is largely determined by matter and energy transformations mediated by microbial community interaction networks. Viral infection modulates network properties through mortality, gene transfer and metabolic reprogramming. Results Here we explore the nature and extent of viral metabolic reprogramming throughout the Pacific Ocean depth continuum. We describe 35 marine viral gene families with potential to reprogram metabolic flux through central metabolic pathways recovered from Pacific Ocean waters. Four of these families have been previously reported but 31 are novel. These known and new carbon pathway auxiliary metabolic genes were recovered from a total of 22 viral metagenomes in which viral auxiliary metabolic genes were differentiated from low-level cellular DNA inputs based on small subunit ribosomal RNA gene content, taxonomy, fragment recruitment and genomic context information. Auxiliary metabolic gene distribution patterns reveal that marine viruses target overlapping, but relatively distinct pathways in sunlit and dark ocean waters to redirect host carbon flux towards energy production and viral genome replication under low nutrient, niche-differentiated conditions throughout the depth continuum. Conclusions Given half of ocean microbes are infected by viruses at any given time, these findings of broad viral metabolic reprogramming suggest the need for renewed consideration of viruses in global ocean carbon models.

The ecological importance of viruses is now widely recognized, yet our limited knowledge of viral sequence space and virus-host interactions precludes accurate prediction of their roles and impacts. In this study, we mined publicly available bacterial and archaeal genomic data sets to identify 12,498 high-confidence viral genomes linked to their microbial hosts. These data augment public data sets 10-fold, provide first viral sequences for 13 new bacterial phyla including ecologically abundant phyla, and help taxonomically identify 7-38% of 'unknown' sequence space in viromes. Genome- and network-based classification was largely consistent with accepted viral taxonomy and suggested that (i) 264 new viral genera were identified (doubling known genera) and (ii) cross-taxon genomic recombination is limited. Further analyses provided empirical data on extrachromosomal prophages and coinfection prevalences, as well as evaluation of in silico virus-host linkage predictions. Together these findings illustrate the value of mining viral signal from microbial genomes.

Background

Viral metagenomics (viromics) is increasingly used to obtain uncultivated viral genomes, evaluate community diversity, and assess ecological hypotheses. While viromic experimental methods are relatively mature and widely accepted by the research community, robust bioinformatics standards remain to be established. Here we used in silico mock viral communities to evaluate the viromic sequence-to-ecological-inference pipeline, including (i) read pre-processing and metagenome assembly, (ii) thresholds applied to estimate viral relative abundances based on read mapping to assembled contigs, and (iii) normalization methods applied to the matrix of viral relative abundances for alpha and beta diversity estimates.

Results

Tools specifically designed for metagenomes, specifically metaSPAdes, MEGAHIT, and IDBA-UD, were the most effective at assembling viromes. Read pre-processing, such as partitioning, had virtually no impact on assembly output, but may be useful when hardware is limited. Viral populations with 2-5 × coverage typically assembled well, whereas lesser coverage led to fragmented assembly. Strain heterogeneity within populations hampered assembly, especially when strains were closely related (average nucleotide identity, or ANI ≥97%) and when the most abundant strain represented <50% of the population. Viral community composition assessments based on read recruitment were generally accurate when the following thresholds for detection were applied: (i) ≥10 kb contig lengths to define populations, (ii) coverage defined from reads mapping at ≥90% identity, and (iii) ≥75% of contig length with ≥1 × coverage. Finally, although data are limited to the most abundant viruses in a community, alpha and beta diversity patterns were robustly estimated (±10%) when comparing samples of similar sequencing depth, but more divergent (up to 80%) when sequencing depth was uneven across the dataset. In the latter cases, the use of normalization methods specifically developed for metagenomes provided the best estimates.

Conclusions

These simulations provide benchmarks for selecting analysis cut-offs and establish that an optimized sample-to-ecological-inference viromics pipeline is robust for making ecological inferences from natural viral communities. Continued development to better accessing RNA, rare, and/or diverse viral populations and improved reference viral genome availability will alleviate many of viromics remaining limitations.

Microbes affect nutrient and energy transformations throughout the world's ecosystems, yet they do so under viral constraints. In complex communities, viral metagenome (virome) sequencing is transforming our ability to quantify viral diversity and impacts. Although some bottlenecks, for example, few reference genomes and nonquantitative viromics, have been overcome, the void of centralized data sets and specialized tools now prevents viromics from being broadly applied to answer fundamental ecological questions. Here we present iVirus, a community resource that leverages the CyVerse cyberinfrastructure to provide access to viromic tools and data sets. The iVirus Data Commons contains both raw and processed data from 1866 samples and 73 projects derived from global ocean expeditions, as well as existing and legacy public repositories. Through the CyVerse Discovery Environment, users can interrogate these data sets using existing analytical tools (software applications known as 'Apps') for assembly, open reading frame prediction and annotation, as well as several new Apps specifically developed for analyzing viromes. Because Apps are web based and powered by CyVerse supercomputing resources, they enable scalable analyses for a broad user base. Finally, a use-case scenario documents how to apply these advances toward new data. This growing iVirus resource should help researchers utilize viromics as yet another tool to elucidate viral roles in nature.

Viruses influence the ecology and evolutionary trajectory of microbial communities. Yet our understanding of their roles in ecosystems is limited by the paucity of model systems available for hypothesis generation and testing. Further, virology is limited by the lack of a broadly accepted conceptual framework to classify viral diversity into evolutionary and ecologically cohesive units. Here, we introduce genomes, structural proteomes, and quantitative host range data for eight Pseudoalteromonas phages isolated from Helgoland (North Sea, Germany) and use these data to advance a genome-based viral operational taxonomic unit (OTU) definition. These viruses represent five new genera and inform 498 unaffiliated or unannotated protein clusters (PCs) from global virus metagenomes. In a comparison of previously sequenced Pseudoalteromonas phage isolates (n = 7) and predicted prophages (n = 31), the eight phages are unique. They share a genus with only one other isolate, Pseudoalteromonas podophage RIO-1 (East Sea, South Korea) and two Pseudoalteromonas prophages. Mass-spectrometry of purified viral particles identified 12-20 structural proteins per phage. When combined with 3-D structural predictions, these data led to the functional characterization of five previously unidentified major capsid proteins. Protein functional predictions revealed mechanisms for hijacking host metabolism and resources. Further, they uncovered a hybrid sipho-myovirus that encodes genes for Mu-like infection rarely described in ocean systems. Finally, we used these data to evaluate a recently introduced definition for virus populations that requires members of the same population to have >95% average nucleotide identity across at least 80% of their genes. Using physiological traits and genomics, we proposed a conceptual model for a viral OTU definition that captures evolutionarily cohesive and ecologically distinct units. In this trait-based framework, sensitive hosts are considered viral niches, while host ranges and infection efficiencies are tracked as viral traits. Quantitative host range assays revealed conserved traits within virus OTUs that break down between OTUs, suggesting the defined units capture niche and fitness differentiation. Together these analyses provide a foundation for model system-based hypothesis testing that will improve our understanding of marine copiotrophs, as well as phage-host interactions on the ocean particles and aggregates where Pseudoalteromonas thrive.

The Au’au Channel between the islands of Maui and Lanai, Hawaii comprises critical breeding habitat for humpback whales(Megaptera novaeangliae) of the Central North Pacific stock. However, like many regions where marine mega-fauna gather,these waters are also the focus of a flourishing local eco-tourism and whale watching industry. Our aim was to establishcurrent trends in habitat preference in female-calf humpback whale pairs within this region, focusing specifically on thebusy, eastern portions of the channel. We used an equally-spaced zigzag transect survey design, compiled our results in aGIS model to identify spatial trends and calculated Neu’s Indices to quantify levels of habitat use. Our study revealed thatwhile mysticete female-calf pairs on breeding grounds typically favor shallow, inshore waters, female-calf pairs in the Au’auChannel avoided shallow waters (,20 m) and regions within 2 km of the shoreline. Preferred regions for female-calf pairscomprised water depths between 40–60 m, regions of rugged bottom topography and regions that lay between 4 and6 km from a small boat harbor (Lahaina Harbor) that fell within the study area. In contrast to other humpback whalebreeding grounds, there was only minimal evidence of typical patterns of stratification or segregation according to groupcomposition. A review of habitat use by maternal females across Hawaiian waters indicates that maternal habitat choicevaries between localities within the Hawaiian Islands, suggesting that maternal females alter their use of habitat accordingto locally varying pressures. This ability to respond to varying environments may be the key that allows wildlife species topersist in regions where human activity and critical habitat overlap.

Antibiotic resistance genes (ARGs) are pervasive in gut microbiota, but it remains unclear how often ARGs are transferred, particularly to pathogens. Traditionally, ARG spread is attributed to horizontal transfer mediated either by DNA transformation, bacterial conjugation or generalized transduction. However, recent viral metagenome (virome) analyses suggest that ARGs are frequently carried by phages, which is inconsistent with the traditional view that phage genomes rarely encode ARGs. Here we used exploratory and conservative bioinformatic strategies found in the literature to detect ARGs in phage genomes, and experimentally assessed a subset of ARG predicted using exploratory thresholds. ARG abundances in 1181 phage genomes were vastly overestimated using exploratory thresholds (421 predicted vs 2 known), due to low similarities and matches to protein unrelated to antibiotic resistance. Consistent with this, four ARGs predicted using exploratory thresholds were experimentally evaluated and failed to confer antibiotic resistance in Escherichia coli. Reanalysis of available human- or mouse-associated viromes for ARGs and their genomic context suggested that bona fide ARG attributed to phages in viromes were previously overestimated. These findings provide guidance for documentation of ARG in viromes, and reassert that ARGs are rarely encoded in phages.

Taxonomic classification of archaeal and bacterial viruses is challenging, yet also fundamental for developing a predictive understanding of microbial ecosystems. Recent identification of hundreds of thousands of new viral genomes and genome fragments, whose hosts remain unknown, requires a paradigm shift away from traditional classification approaches and towards the use of genomes for taxonomy. Here we revisited the use of genomes and their protein content as a means for developing a viral taxonomy for bacterial and archaeal viruses. A network-based analytic was evaluated and benchmarked against authority-accepted taxonomic assignments and found to be largely concordant. Exceptions were manually examined and found to represent areas of viral genome 'sequence space' that are under-sampled or prone to excessive genetic exchange. While both cases are poorly resolved by genome-based taxonomic approaches, the former will improve as viral sequence space is better sampled and the latter are uncommon. Finally, given the largely robust taxonomic capabilities of this approach, we sought to enable researchers to easily and systematically classify new viruses. Thus, we established a tool, vConTACT, as an app at iVirus, where it operates as a fast, highly scalable, user-friendly app within the free and powerful CyVerse cyberinfrastructure.

Virophages are small viruses that co-infect eukaryotic cells alongside giant viruses (Mimiviridae) and hijack their machinery to replicate. While two types of virophages have been isolated, their genomic diversity and ecology remain largely unknown. Here we use time series metagenomics to identify and study the dynamics of 25 uncultivated virophage populations, 17 of which represented by complete or near-complete genomes, in two North American freshwater lakes. Taxonomic analysis suggests that these freshwater virophages represent at least three new candidate genera. Ecologically, virophage populations are repeatedly detected over years and evolutionary stable, yet their distinct abundance profiles and gene content suggest that virophage genera occupy different ecological niches. Co-occurrence analyses reveal 11 virophages strongly associated with uncultivated Mimiviridae, and three associated with eukaryotes among the Dinophyceae, Rhizaria, Alveolata, and Cryptophyceae groups. Together, these findings significantly augment virophage databases, help refine virophage taxonomy, and establish baseline ecological hypotheses and tools to study virophages in nature.Virophages are recently-identified small viruses that infect larger viruses, yet their diversity and ecological roles are poorly understood. Here, Roux and colleagues present time series metagenomics data revealing new virophage genera and their putative ecological interactions in two freshwater lakes.