Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses.

Ocean microbes drive biogeochemical cycling on a global scale. However, this cycling is constrained by viruses that affect community composition, metabolic activity, and evolutionary trajectories. Owing to challenges with the sampling and cultivation of viruses, genome-level viral diversity remains poorly described and grossly understudied, with less than 1% of observed surface-ocean viruses known. Here we assemble complete genomes and large genomic fragments from both surface- and deep-ocean viruses sampled during the Tara Oceans and Malaspina research expeditions, and analyse the resulting 'global ocean virome' dataset to present a global map of abundant, double-stranded DNA viruses complete with genomic and ecological contexts. A total of 15,222 epipelagic and mesopelagic viral populations were identified, comprising 867 viral clusters (defined as approximately genus-level groups). This roughly triples the number of known ocean viral populations and doubles the number of candidate bacterial and archaeal virus genera, providing a near-complete sampling of epipelagic communities at both the population and viral-cluster level. We found that 38 of the 867 viral clusters were locally or globally abundant, together accounting for nearly half of the viral populations in any global ocean virome sample. While two-thirds of these clusters represent newly described viruses lacking any cultivated representative, most could be computationally linked to dominant, ecologically relevant microbial hosts. Moreover, we identified 243 viral-encoded auxiliary metabolic genes, of which only 95 were previously known. Deeper analyses of four of these auxiliary metabolic genes (dsrC, soxYZ, P-II (also known as glnB) and amoC) revealed that abundant viruses may directly manipulate sulfur and nitrogen cycling throughout the epipelagic ocean. This viral catalog and functional analyses provide a necessary foundation for the meaningful integration of viruses into ecosystem models where they act as key players in nutrient cycling and trophic networks.

Plankton networks driving carbon export in the oligotrophic ocean.

The biological carbon pump is the process by which CO2 is transformed to organic carbon via photosynthesis, exported through sinking particles, and finally sequestered in the deep ocean. While the intensity of the pump correlates with plankton community composition, the underlying ecosystem structure driving the process remains largely uncharacterized. Here we use environmental and metagenomic data gathered during the Tara Oceans expedition to improve our understanding of carbon export in the oligotrophic ocean. We show that specific plankton communities, from the surface and deep chlorophyll maximum, correlate with carbon export at 150 m and highlight unexpected taxa such as Radiolaria and alveolate parasites, as well as Synechococcus and their phages, as lineages most strongly associated with carbon export in the subtropical, nutrient-depleted, oligotrophic ocean. Additionally, we show that the relative abundance of a few bacterial and viral genes can predict a significant fraction of the variability in carbon export in these regions.

Illuminating structural proteins in viral "dark matter" with metaproteomics.

Viruses are ecologically important, yet environmental virology is limited by dominance of unannotated genomic sequences representing taxonomic and functional "viral dark matter." Although recent analytical advances are rapidly improving taxonomic annotations, identifying functional dark matter remains problematic. Here, we apply paired metaproteomics and dsDNA-targeted metagenomics to identify 1,875 virion-associated proteins from the ocean. Over one-half of these proteins were newly functionally annotated and represent abundant and widespread viral metagenome-derived protein clusters (PCs). One primarily unannotated PC dominated the dataset, but structural modeling and genomic context identified this PC as a previously unidentified capsid protein from multiple uncultivated tailed virus families. Furthermore, four of the five most abundant PCs in the metaproteome represent capsid proteins containing the HK97-like protein fold previously found in many viruses that infect all three domains of life. The dominance of these proteins within our dataset, as well as their global distribution throughout the world's oceans and seas, supports prior hypotheses that this HK97-like protein fold is the most abundant biological structure on Earth. Together, these culture-independent analyses improve virion-associated protein annotations, facilitate the investigation of proteins within natural viral communities, and offer a high-throughput means of illuminating functional viral dark matter.

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

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.

An inexpensive, accurate, and precise wet-mount method for enumerating aquatic viruses.

Viruses affect biogeochemical cycling, microbial mortality, gene flow, and metabolic functions in diverse environments through infection and lysis of microorganisms. Fundamental to quantitatively investigating these roles is the determination of viral abundance in both field and laboratory samples. One current, widely used method to accomplish this with aquatic samples is the "filter mount" method, in which samples are filtered onto costly 0.02-µm-pore-size ceramic filters for enumeration of viruses by epifluorescence microscopy. Here we describe a cost-effective (ca. 500-fold-lower materials cost) alternative virus enumeration method in which fluorescently stained samples are wet mounted directly onto slides, after optional chemical flocculation of viruses in samples with viral concentrations of <5×10(7) viruses ml(-1). The concentration of viruses in the sample is then determined from the ratio of viruses to a known concentration of added microsphere beads via epifluorescence microscopy. Virus concentrations obtained by using this wet-mount method, with and without chemical flocculation, were significantly correlated with, and had precision equivalent to, those obtained by the filter mount method across concentrations ranging from 2.17×10(6) to 1.37×10(8) viruses ml(-1) when tested by using cultivated viral isolates and natural samples from marine and freshwater environments. In summary, the wet-mount method is significantly less expensive than the filter mount method and is appropriate for rapid, precise, and accurate enumeration of aquatic viruses over a wide range of viral concentrations (≥1×10(6) viruses ml(-1)) encountered in field and laboratory samples.

Viral community analysis in a marine oxygen minimum zone indicates increased potential for viral manipulation of microbial physiological state.

Microbial communities in oxygen minimum zones (OMZs) are known to have significant impacts on global biogeochemical cycles, but viral influence on microbial processes in these regions are much less studied. Here we provide baseline ecological patterns using microscopy and viral metagenomics from the Eastern Tropical North Pacific (ETNP) OMZ region that enhance our understanding of viruses in these climate-critical systems. While extracellular viral abundance decreased below the oxycline, viral diversity and lytic infection frequency remained high within the OMZ, demonstrating that viral influences on microbial communities were still substantial without the detectable presence of oxygen. Viral community composition was strongly related to oxygen concentration, with viral populations in low-oxygen portions of the water column being distinct from their surface layer counterparts. However, this divergence was not accompanied by the expected differences in viral-encoded auxiliary metabolic genes (AMGs) relating to nitrogen and sulfur metabolisms that are known to be performed by microbial communities in these low-oxygen and anoxic regions. Instead, several abundant AMGs were identified in the oxycline and OMZ that may modulate host responses to low-oxygen stress. We hypothesize that this is due to selection for viral-encoded genes that influence host survivability rather than modulating host metabolic reactions within the ETNP OMZ. Together, this study shows that viruses are not only diverse throughout the water column in the ETNP, including the OMZ, but their infection of microorganisms has the potential to alter host physiological state within these biogeochemically important regions of the ocean.

Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems.

Biogeographical studies have traditionally focused on readily visible organisms, but recent technological advances are enabling analyses of the large-scale distribution of microscopic organisms, whose biogeographical patterns have long been debated. Here we assessed the global structure of plankton geography and its relation to the biological, chemical, and physical context of the ocean (the 'seascape') by analyzing metagenomes of plankton communities sampled across oceans during the Tara Oceans expedition, in light of environmental data and ocean current transport. Using a consistent approach across organismal sizes that provides unprecedented resolution to measure changes in genomic composition between communities, we report a pan-ocean, size-dependent plankton biogeography overlying regional heterogeneity. We found robust evidence for a basin-scale impact of transport by ocean currents on plankton biogeography, and on a characteristic timescale of community dynamics going beyond simple seasonality or life history transitions of plankton.

Oceans are brimming with life invisible to our eyes, a myriad of species of bacteria, viruses and other microscopic organisms essential for the health of the planet. These 'marine plankton' are unable to swim against currents and should therefore be constantly on the move, yet previous studies have suggested that distinct species of plankton may in fact inhabit different oceanic regions. However, proving this theory has been challenging; collecting plankton is logistically difficult, and it is often impossible to distinguish between species simply by examining them under a microscope. However, within the last decade, a research schooner called Tara has travelled the globe to gather thousands of plankton samples. At the same time, advances in genomics have made it possible to identify species based only on fragments of their DNA sequence. To understand the hidden geography of plankton communities in Earth's oceans, Richter et al. pored over DNA from the Tara Oceans expedition. This revealed that, despite being unable to resist the flow of water, various planktonic species which live close to the surface manage to occupy distinct, stable provinces shaped by currents. Different sizes of plankton are distributed in different sized provinces, with the smallest organisms tending to inhabit the smallest areas. Comparing DNA similarities and speeds of currents at the ocean surface revealed how these might stretch and mix plankton communities. Plankton play a critical role in the health of the ocean and the chemical cycles of planet Earth. These results could allow deeper investigation by marine modellers, ecologists, and evolutionary biologists. Meanwhile, work is already underway to investigate how climate change might impact this hidden geography.

Microbial Ecology of Oxygen Minimum Zones Amidst Ocean Deoxygenation.

Oxygen minimum zones (OMZs) have substantial effects on the global ecology and biogeochemical processes of marine microbes. However, the diversity and activity of OMZ microbes and their trophic interactions are only starting to be documented, especially in regard to the potential roles of viruses and protists. OMZs have expanded over the past 60 years and are predicted to expand due to anthropogenic climate change, furthering the need to understand these regions. This review summarizes the current knowledge of OMZ formation, the biotic and abiotic factors involved in OMZ expansion, and the microbial ecology of OMZs, emphasizing the importance of bacteria, archaea, viruses, and protists. We describe the recognized roles of OMZ microbes in carbon, nitrogen, and sulfur cycling, the potential of viruses in altering host metabolisms involved in these cycles, and the control of microbial populations by grazers and viruses. Further, we highlight the microbial community composition and roles of these organisms in oxic and anoxic depths within the water column and how these differences potentially inform how microbial communities will respond to deoxygenation. Additionally, the current literature on the alteration of microbial communities by other key climate change parameters such as temperature and pH are considered regarding how OMZ microbes might respond to these pressures. Finally, we discuss what knowledge gaps are present in understanding OMZ microbial communities and propose directions that will begin to close these gaps.

Morphological characterization of viruses in the stratified water column of alkaline, hypersaline Mono Lake.

Concentrations of viruses and prokaryotes in the alkaline, moderately hypersaline, seasonally stratified Mono Lake are among the highest reported for a natural aquatic environment. We used electron microscopy to test whether viral morphological characteristics differed among the epilimnion, metalimnion, and the anoxic hypolimnion of the lake and to determine how the properties of viruses in Mono Lake compare to other aquatic environments. Viral capsid size distributions were more similar in the metalimnion and hypolimnion of Mono Lake, while viral tail lengths were more similar in the epilimnion and metalimnion. The percentage of tailed viruses decreased with depth and the relative percentages of tailed phage families changed with depth. The presence of large (>125 nm capsid), untailed viruses in the metalimnion and hypolimnion suggests that eukaryotic viruses are produced in these suboxic and anoxic, hypersaline environments. Capsid diameters of viruses were larger on average in Mono Lake compared to other aquatic environments, and no lemon-shaped or filamentous viruses were found, in contrast to other high-salinity or high-altitude lakes and seas. Our data suggest that the physically and chemically distinct layers of Mono Lake harbor different viral assemblages, and that these assemblages are distinct from other aquatic environments that have been studied. Furthermore, we found that filtration of a sample through a 0.22-µm pore-size filter significantly altered the distribution of viral capsid diameters and tail lengths, resulting in a relative depletion of viruses having larger capsids and longer tails. This observation highlights the potential for bias in molecular surveys of viral diversity, which typically rely on filtration through 0.2- or 0.22-µm pore-size membrane filters to remove bacteria during sample preparation.

Host-linked soil viral ecology along a permafrost thaw gradient.

Climate change threatens to release abundant carbon that is sequestered at high latitudes, but the constraints on microbial metabolisms that mediate the release of methane and carbon dioxide are poorly understood1-7. The role of viruses, which are known to affect microbial dynamics, metabolism and biogeochemistry in the oceans8-10, remains largely unexplored in soil. Here, we aimed to investigate how viruses influence microbial ecology and carbon metabolism in peatland soils along a permafrost thaw gradient in Sweden. We recovered 1,907 viral populations (genomes and large genome fragments) from 197 bulk soil and size-fractionated metagenomes, 58% of which were detected in metatranscriptomes and presumed to be active. In silico predictions linked 35% of the viruses to microbial host populations, highlighting likely viral predators of key carbon-cycling microorganisms, including methanogens and methanotrophs. Lineage-specific virus/host ratios varied, suggesting that viral infection dynamics may differentially impact microbial responses to a changing climate. Virus-encoded glycoside hydrolases, including an endomannanase with confirmed functional activity, indicated that viruses influence complex carbon degradation and that viral abundances were significant predictors of methane dynamics. These findings suggest that viruses may impact ecosystem function in climate-critical, terrestrial habitats and identify multiple potential viral contributions to soil carbon cycling.

Putative archaeal viruses from the mesopelagic ocean.

Oceanic viruses that infect bacteria, or phages, are known to modulate host diversity, metabolisms, and biogeochemical cycling, while the viruses that infect marine Archaea remain understudied despite the critical ecosystem roles played by their hosts. Here we introduce "MArVD", for Metagenomic Archaeal Virus Detector, an annotation tool designed to identify putative archaeal virus contigs in metagenomic datasets. MArVD is made publicly available through the online iVirus analytical platform. Benchmarking analysis of MArVD showed it to be >99% accurate and 100% sensitive in identifying the 127 known archaeal viruses among the 12,499 viruses in the VirSorter curated dataset. Application of MArVD to 10 viral metagenomes from two depth profiles in the Eastern Tropical North Pacific (ETNP) oxygen minimum zone revealed 43 new putative archaeal virus genomes and large genome fragments ranging in size from 10 to 31 kb. Network-based classifications, which were consistent with marker gene phylogenies where available, suggested that these putative archaeal virus contigs represented six novel candidate genera. Ecological analyses, via fragment recruitment and ordination, revealed that the diversity and relative abundances of these putative archaeal viruses were correlated with oxygen concentration and temperature along two OMZ-spanning depth profiles, presumably due to structuring of the host Archaea community. Peak viral diversity and abundances were found in surface waters, where Thermoplasmata 16S rRNA genes are prevalent, suggesting these archaea as hosts in the surface habitats. Together these findings provide a baseline for identifying archaeal viruses in sequence datasets, and an initial picture of the ecology of such viruses in non-extreme environments.

Viral to metazoan marine plankton nucleotide sequences from the Tara Oceans expedition.

A unique collection of oceanic samples was gathered by the Tara Oceans expeditions (2009-2013), targeting plankton organisms ranging from viruses to metazoans, and providing rich environmental context measurements. Thanks to recent advances in the field of genomics, extensive sequencing has been performed for a deep genomic analysis of this huge collection of samples. A strategy based on different approaches, such as metabarcoding, metagenomics, single-cell genomics and metatranscriptomics, has been chosen for analysis of size-fractionated plankton communities. Here, we provide detailed procedures applied for genomic data generation, from nucleic acids extraction to sequence production, and we describe registries of genomics datasets available at the European Nucleotide Archive (ENA, www.ebi.ac.uk/ena). The association of these metadata to the experimental procedures applied for their generation will help the scientific community to access these data and facilitate their analysis. This paper complements other efforts to provide a full description of experiments and open science resources generated from the Tara Oceans project, further extending their value for the study of the world's planktonic ecosystems.

Seasonal time bombs: dominant temperate viruses affect Southern Ocean microbial dynamics.

Rapid warming in the highly productive western Antarctic Peninsula (WAP) region of the Southern Ocean has affected multiple trophic levels, yet viral influences on microbial processes and ecosystem function remain understudied in the Southern Ocean. Here we use cultivation-independent quantitative ecological and metagenomic assays, combined with new comparative bioinformatic techniques, to investigate double-stranded DNA viruses during the WAP spring-summer transition. This study demonstrates that (i) temperate viruses dominate this region, switching from lysogeny to lytic replication as bacterial production increases, and (ii) Southern Ocean viral assemblages are genetically distinct from lower-latitude assemblages, primarily driven by this temperate viral dominance. This new information suggests fundamentally different virus-host interactions in polar environments, where intense seasonal changes in bacterial production select for temperate viruses because of increased fitness imparted by the ability to switch replication strategies in response to resource availability. Further, temperate viral dominance may provide mechanisms (for example, bacterial mortality resulting from prophage induction) that help explain observed temporal delays between, and lower ratios of, bacterial and primary production in polar versus lower-latitude marine ecosystems. Together these results suggest that temperate virus-host interactions are critical to predicting changes in microbial dynamics brought on by warming in polar marine systems.

Rising to the challenge: accelerated pace of discovery transforms marine virology.

Marine viruses have important roles in microbial mortality, gene transfer, metabolic reprogramming and biogeochemical cycling. In this Review, we discuss recent technological advances in marine virology including the use of near-quantitative, reproducible metagenomics for large-scale investigation of viral communities and the emergence of gene-based viral ecology. We also describe the reprogramming of microbially driven processes by viral metabolic genes, the identification of novel viruses using cultivation-dependent and cultivation-independent tools, and the potential for modelling studies to provide a framework for studying virus-host interactions. These transformative advances have set a rapid pace in exploring and predicting how marine viruses manipulate and respond to their environment.

Depth-stratified functional and taxonomic niche specialization in the 'core' and 'flexible' Pacific Ocean Virome.

Microbes drive myriad ecosystem processes, and their viruses modulate microbial-driven processes through mortality, horizontal gene transfer, and metabolic reprogramming by viral-encoded auxiliary metabolic genes (AMGs). However, our knowledge of viral roles in the oceans is primarily limited to surface waters. Here we assess the depth distribution of protein clusters (PCs) in the first large-scale quantitative viral metagenomic data set that spans much of the pelagic depth continuum (the Pacific Ocean Virome; POV). This established 'core' (180 PCs; one-third new to science) and 'flexible' (423K PCs) community gene sets, including niche-defining genes in the latter (385 and 170 PCs are exclusive and core to the photic and aphotic zones, respectively). Taxonomic annotation suggested that tailed phages are ubiquitous, but not abundant (<5% of PCs) and revealed depth-related taxonomic patterns. Functional annotation, coupled with extensive analyses to document non-viral DNA contamination, uncovered 32 new AMGs (9 core, 20 photic and 3 aphotic) that introduce ways in which viruses manipulate infected host metabolism, and parallel depth-stratified host adaptations (for example, photic zone genes for iron-sulphur cluster modulation for phage production, and aphotic zone genes for high-pressure deep-sea survival). Finally, significant vertical flux of photic zone viruses to the deep sea was detected, which is critical for interpreting depth-related patterns in nature. Beyond the ecological advances outlined here, this catalog of viral core, flexible and niche-defining genes provides a resource for future investigation into the organization, function and evolution of microbial molecular networks to mechanistically understand and model viral roles in the biosphere.

Ocean plankton. Patterns and ecological drivers of ocean viral communities.

Viruses influence ecosystems by modulating microbial population size, diversity, metabolic outputs, and gene flow. Here, we use quantitative double-stranded DNA (dsDNA) viral-fraction metagenomes (viromes) and whole viral community morphological data sets from 43 Tara Oceans expedition samples to assess viral community patterns and structure in the upper ocean. Protein cluster cataloging defined pelagic upper-ocean viral community pan and core gene sets and suggested that this sequence space is well-sampled. Analyses of viral protein clusters, populations, and morphology revealed biogeographic patterns whereby viral communities were passively transported on oceanic currents and locally structured by environmental conditions that affect host community structure. Together, these investigations establish a global ocean dsDNA viromic data set with analyses supporting the seed-bank hypothesis to explain how oceanic viral communities maintain high local diversity.

Ocean plankton. Environmental characteristics of Agulhas rings affect interocean plankton transport.

Agulhas rings provide the principal route for ocean waters to circulate from the Indo-Pacific to the Atlantic basin. Their influence on global ocean circulation is well known, but their role in plankton transport is largely unexplored. We show that, although the coarse taxonomic structure of plankton communities is continuous across the Agulhas choke point, South Atlantic plankton diversity is altered compared with Indian Ocean source populations. Modeling and in situ sampling of a young Agulhas ring indicate that strong vertical mixing drives complex nitrogen cycling, shaping community metabolism and biogeochemical signatures as the ring and associated plankton transit westward. The peculiar local environment inside Agulhas rings may provide a selective mechanism contributing to the limited dispersal of Indian Ocean plankton populations into the Atlantic.