Seasonality and monthly dynamics of marine myovirus communities

Phylogenetic proximity drives temporal succession of marine giant viruses in a five-year metagenomic time-series

Nucleocytoplasmic Large DNA Viruses (NCLDVs, also called giant viruses) are widespread in marine systems and infect a broad range of microbial eukaryotes (protists). Recent biogeographic work has provided global snapshots of NCLDV diversity and community composition across the world's oceans, yet little information exists about the guiding 'rules' underpinning their community dynamics over time. We leveraged a five-year monthly metagenomic time-series to quantify the community composition of NCLDVs off the coast of Southern California and characterize these populations' temporal dynamics. NCLDVs were dominated by Algavirales (Phycodnaviruses, 59%) and Imitervirales (Mimiviruses, 36%). We identified clusters of NCLDVs with distinct classes of seasonal and non-seasonal temporal dynamics. Overall, NCLDV population abundances were often highly dynamic with a strong seasonal signal. The Imitervirales group had highest relative abundance in the more oligotrophic late summer and fall, while Algavirales did so in winter. Generally, closely related strains had similar temporal dynamics, suggesting that evolutionary history is a key driver of the temporal niche of marine NCLDVs. However, a few closely-related strains had drastically different seasonal dynamics, suggesting that while phylogenetic proximity often indicates ecological similarity, occasionally phenology can shift rapidly, possibly due to host-switching. Finally, we identified distinct functional content and possible host interactions of two major NCLDV orders-including connections of Imitervirales with primary producers like the diatom Chaetoceros and widespread marine grazers like Paraphysomonas and Spirotrichea ciliates. Together, our results reveal key insights on season-specific effect of phylogenetically distinct giant virus communities on marine protist metabolism, biogeochemical fluxes and carbon cycling.

The structure and assembly mechanisms of T4-like cyanophages community in the South China Sea

Marine ecosystems contain an immense diversity of phages, many of which infect cyanobacteria (cyanophage) that are largely responsible for primary productivity. To characterize the genetic diversity and biogeographic distribution of the marine T4-like cyanophage community in the northern South China Sea, the T4-like cyanophage portal protein gene (g20) was amplified. Phylogenetic analysis revealed that marine T4-like cyanophages were highly diverse, with g20 operational taxonomic units being affiliated with five defined clades (Clusters I-V). Cluster II had a wide geographic distribution, Cluster IV was the most abundant in the open sea, and Cluster I was dominant in coastal shelf environments. Our results showed T4-like cyanophages (based on g20) community was generally shaped via heterogeneous selection. Highly variable environmental factors (such as salinity and temperature) can heterogeneously select different cyanophage communities. Nevertheless, the dominant drivers of the T4-like cyanophage community based on the g20 and g23 (T4-like phage major capsid protein gene) were different, probably due to different coverages by the primer sets. Furthermore, the community assembly processes of T4-like cyanophages were affected by host traits (abundance and distribution), viral traits (latent period, burst size, and host range), and environmental properties (temperature and salinity).IMPORTANCECyanophages are abundant and ubiquitous in the oceans, altering population structures and evolution of cyanobacteria, which account for a large portion of global carbon fixation, through host mortality, horizontal gene transfer, and the modulation of host metabolism. However, little is known about the biogeography and ecological drivers that shape the cyanophage community. Here, we use g20 and g23 genes to examine the biogeographic patterns and the assembly mechanisms of T4-like cyanophage community in the northern part of the South China Sea. The different coverages of primer sets might lead to the different dominant drivers of T4-like cyanophage community based on g20 and g23 genes. Our results showed that characteristics of viral traits (latent period, burst size, and host range) and host traits (abundance and distribution) were found to either limit or enhance the biogeographic distribution of T4-like cyanophages. Overall, both virus and host properties are critical to consider when determining rules of community assembly for viruses.

Spatial and temporal metagenomics of river compartments reveals viral community dynamics in an urban impacted stream

Although river ecosystems comprise less than 1% of Earth's total non-glaciated area, they are critical modulators of microbially and virally orchestrated global biogeochemical cycles. However, most studies either use data that is not spatially resolved or is collected at timepoints that do not reflect the short life cycles of microorganisms. As a result, the relevance of microbiome interactions and the impacts they have over time on biogeochemical cycles are poorly understood. To assess how viral and microbial communities change over time, we sampled surface water and pore water compartments of the wastewater-impacted River Erpe in Germany every 3 hours over a 48-hour period resulting in 32 metagenomes paired to geochemical and metabolite measurements. We reconstructed 6,500 viral and 1,033 microbial genomes and found distinct communities associated with each river compartment. We show that 17% of our vMAGs clustered to viruses from other ecosystems like wastewater treatment plants and rivers. Our results also indicated that 70% of the viral community was persistent in surface waters, whereas only 13% were persistent in the pore waters taken from the hyporheic zone. Finally, we predicted linkages between 73 viral genomes and 38 microbial genomes. These putatively linked hosts included members of the Competibacteraceae, which we suggest are potential contributors to carbon and nitrogen cycling. Together, these findings demonstrate that microbial and viral communities in surface waters of this urban river can exist as stable communities along a flowing river; and raise important considerations for ecosystem models attempting to constrain dynamics of river biogeochemical cycles.

Prevalence of Viral Frequency-Dependent Infection in Coastal Marine Prokaryotes Revealed Using Monthly Time Series Virome Analysis

Viruses infecting marine prokaryotes have a large impact on the diversity and dynamics of their hosts. Model systems suggest that viral infection is frequency dependent and constrained by the virus-host encounter rate. However, it is unclear whether frequency-dependent infection is pervasive among the abundant prokaryotic populations with different temporal dynamics. To address this question, we performed a comparison of prokaryotic and viral communities using 16S rRNA amplicon and virome sequencing based on samples collected monthly for 2 years at a Japanese coastal site, Osaka Bay. Concurrent seasonal shifts observed in prokaryotic and viral community dynamics indicated that the abundance of viruses correlated with that of their predicted host phyla (or classes). Cooccurrence network analysis between abundant prokaryotes and viruses revealed 6,423 cooccurring pairs, suggesting a tight coupling of host and viral abundances and their "one-to-many" correspondence. Although stable dominant species, such as SAR11, showed few cooccurring viruses, a fast succession of their viruses suggests that viruses infecting these populations changed continuously. Our results suggest that frequency-dependent viral infection prevails in coastal marine prokaryotes regardless of host taxa and temporal dynamics. IMPORTANCE There is little room for doubt that viral infection is prevalent among abundant marine prokaryotes regardless of their taxa or growth strategy. However, comprehensive evaluations of viral infections in natural prokaryotic communities are still technically difficult. In this study, we examined viral infection in abundant prokaryotes by monitoring the monthly dynamics of prokaryotic and viral communities at a eutrophic coastal site, Osaka Bay. We compared the community dynamics of viruses with those of their putative hosts based on genome-based in silico host prediction. We observed frequent cooccurrence among the predicted virus-host pairs, suggesting that viral infection is prevalent in abundant prokaryotes regardless of their taxa or temporal dynamics. This likely indicates that frequent lysis of the abundant prokaryotes via viral infection has a considerable contribution to the biogeochemical cycling and maintenance of prokaryotic community diversity.

Diverse Marine T4-like Cyanophage Communities Are Primarily Comprised of Low-Abundance Species Including Species with Distinct Seasonal, Persistent, Occasional, or Sporadic Dynamics

Cyanophages exert important top-down controls on their cyanobacteria hosts; however, concurrent analysis of both phage and host populations is needed to better assess phage-host interaction models. We analyzed picocyanobacteria Prochlorococcus and Synechococcus and T4-like cyanophage communities in Pacific Ocean surface waters using five years of monthly viral and cellular fraction metagenomes. Cyanophage communities contained thousands of mostly low-abundance (<2% relative abundance) species with varying temporal dynamics, categorized as seasonally recurring or non-seasonal and occurring persistently, occasionally, or sporadically (detected in ≥85%, 15-85%, or <15% of samples, respectively). Viromes contained mostly seasonal and persistent phages (~40% each), while cellular fraction metagenomes had mostly sporadic species (~50%), reflecting that these sample sets capture different steps of the infection cycle-virions from prior infections or within currently infected cells, respectively. Two groups of seasonal phages correlated to Synechococcus or Prochlorococcus were abundant in spring/summer or fall/winter, respectively. Cyanophages likely have a strong influence on the host community structure, as their communities explained up to 32% of host community variation. These results support how both seasonally recurrent and apparent stochastic processes, likely determined by host availability and different host-range strategies among phages, are critical to phage-host interactions and dynamics, consistent with both the Kill-the-Winner and the Bank models.

A Tale of Two Seasons: Distinct Seasonal Viral Communities in a Thermokarst Lake

Thermokarst lakes are important features of subarctic landscapes and are a substantial source of greenhouse gases, although the extent of gas produced varies seasonally. Microbial communities are responsible for the production of methane and CO₂ but the "top down" forces that influence microbial dynamics (i.e., grazers and viruses) and how they vary temporally within these lakes are still poorly understood. The aim of this study was to examine viral diversity over time to elucidate the seasonal structure of the viral communities in thermokarst lakes. We produced virus-enriched metagenomes from a subarctic peatland thermokarst lake in the summer and winter over three years. The vast majority of vOTUs assigned to viral families belonged to Caudovirales (Caudoviricetes), notably the morphological groups myovirus, siphovirus and podovirus. We identified two distinct communities: a dynamic, seasonal community in the oxygenated surface layer during the summer and a stable community found in the anoxic water layer at the bottom of the lake in summer and throughout much of the water column in winter. Comparison with other permafrost and northern lake metagenomes highlighted the distinct composition of viral communities in this permafrost thaw lake ecosystem.

Effects of phytoplankton, viral communities, and warming on free-living and particle-associated marine prokaryotic community structure

Free-living and particle-associated marine prokaryotes have physiological, genomic, and phylogenetic differences, yet factors influencing their temporal dynamics remain poorly constrained. In this study, we quantify the entire microbial community composition monthly over several years, including viruses, prokaryotes, phytoplankton, and total protists, from the San-Pedro Ocean Time-series using ribosomal RNA sequencing and viral metagenomics. Canonical analyses show that in addition to physicochemical factors, the double-stranded DNA viral community is the strongest factor predicting free-living prokaryotes, explaining 28% of variability, whereas the phytoplankton (via chloroplast 16S rRNA) community is strongest with particle-associated prokaryotes, explaining 31% of variability. Unexpectedly, protist community explains little variability. Our findings suggest that biotic interactions are significant determinants of the temporal dynamics of prokaryotes, and the relative importance of specific interactions varies depending on lifestyles. Also, warming influenced the prokaryotic community, which largely remained oligotrophic summer-like throughout 2014-15, with cyanobacterial populations shifting from cold-water ecotypes to warm-water ecotypes.

Isolation, characterization, and comparative genomic analysis of vB_PlaM_Pd22F, a new bacteriophage of the family Myoviridae

The use of phage and phage-based products for the prevention and treatment of bee disease is one of the promising natural alternatives to chemical or antibiotic treatments in beekeeping. A novel lysogenic bacteriophage, phage Pd22F (vB_PlaM_Pd22F), was isolated from Paenibacillus dendritiformis by the prophage induction method. This phage, which is capable of infecting Paenibacillus larvae and P. dendritiformis strains, was characterized by microbiological and comparative genomic analysis. Transmission electron microscopy images showed that phage Pd22F had the morphology of a myovirus. Whole-genome sequencing results showed that vB_Pla M_Pd22F has an 86,388-bp linear dsDNA genome with a GC content of 50.68%. This genome has 124 coding sequences (CDSs), 53% of which encode functionally unknown proteins and 57 of which encode proteins that show similarity to known proteins. In addition, one tRNA gene was found. The phage Pd22F genome does not contain any antimicrobial resistance genes. The similarity between the genome sequence of phage Pd22F and the whole genome sequences of other Paenibacillus phages available in the NCBI Virus Database was found to be below 50% (42%), indicating that phage Pd22F differs greatly from previously characterized phages at the DNA level. The results of comparative genomics and phylogenetic analysis revealed that Pd22F is a new phage belonging to the family Myoviridae, order Caudovirales. This is the first report of genomic and morphological characterization of a Paenibacillus dendritiformis prophage.

Viral diversity is linked to bacterial community composition in alpine stream biofilms

Biofilms play pivotal roles in fluvial ecosystems, yet virtually nothing is known about viruses in these communities. Leveraging an optimized sample-to-sequence pipeline, we studied the spatiotemporal turnover of dsDNA viruses associated with stream biofilms and found an astounding diversity to be structured by seasons and along the longitudinal gradient in the stream. While some vOTUs were region- or season-specific, we also identified a large group of permanent biofilm phages, taxonomically dominated by Myoviridae. Comparison of the observed viral distribution with predictions based on neutral community assembly indicated that chance and dispersal may be important for structuring stream biofilm viral communities. Deviation from neutral model predictions suggests that certain phages distribute efficiently across distant locations within the stream network. This dispersal capacity appears to be linked to EPS depolymerases that enable phages to efficiently overcome the biofilm barrier. Other phages, particularly vOTUs classified as Siphoviridae, appear locally overrepresented and to rely on a lysogenic life cycle, potentially to exploit the spatial distribution of bacterial populations in stream biofilms. Overall, biofilm viral and bacterial community turnover were significantly coupled. Yet, viral communities were linked to the presence of the most abundant bacterial community members. With this work, we provide a foundational ecological perspective on factors that structure viral diversity in stream biofilms and identify potentially important viral traits related to the biofilm mode of life.

Year-round dynamics of amplicon sequence variant communities differ among eukaryotes, Imitervirales, and prokaryotes in a coastal ecosystem

Coastal microbial communities are affected by seasonal environmental change, biotic interactions, and fluctuating nutrient availability. We investigated the seasonal dynamics of communities of eukaryotes, a major group of double-stranded DNA viruses that infect eukaryotes (order Imitervirales; phylum Nucleocytoviricota), and prokaryotes in the Uranouchi Inlet, Kochi, Japan. We performed metabarcoding using ribosomal RNA genes and viral polB genes as markers in 43 seawater samples collected over 20 months. Eukaryotes, prokaryotes, and Imitervirales communities characterized by the compositions of amplicon sequence variants (ASVs) showed synchronic seasonal cycles. However, the community dynamics showed intriguing differences in several aspects, such as the recovery rate after a year. We also showed that the differences in community dynamics were at least partially explained by differences in recurrence/persistence levels of individual ASVs among eukaryotes, prokaryotes, and Imitervirales. Prokaryotic ASVs were the most persistent, followed by eukaryotic ASVs and Imitervirales ASVs, which were the least persistent. We argue that the differences in the specificity of interactions (virus-eukaryote vs. prokaryote-eukaryote) as well as the niche breadth of community members were at the origin of the distinct community dynamics among eukaryotes, their viruses, and prokaryotes.

Role of Phylogenetic Structure in the Dynamics of Coastal Viral Assemblages

Marine microbes, including viruses, are an essential part of the marine ecosystem, forming the base of the food web and driving biogeochemical cycles. Within this system, the composition of viral assemblages changes markedly with time, and some of these changes are repeatable through time; however, the extent to which these dynamics are reflected within versus among evolutionarily related groups of viruses is largely unexplored. To examine these dynamics, changes in the composition of two groups of ecologically important viruses and communities of their potential hosts were sampled every 2 weeks for 13 months at a coastal site in British Columbia, Canada. We sequenced two marker genes for viruses-the gene encoding the major capsid protein of T4-like phages and their relatives (gp23) and the RNA-dependent RNA polymerase (RdRp) gene of marnavirus-like RNA viruses-as well as marker genes for their bacterial and eukaryotic host communities, the genes encoding 16S rRNA and 18S rRNA. There were strong lagged correlations between viral diversity and community similarity of putative hosts, implying that the viruses influenced the composition of the host communities. The results showed that for both viral assemblages, the dominant clusters of phylogenetically related viruses shifted over time, and this was correlated with environmental changes. Viral clusters contained many ephemeral taxa and few persistent taxa, but within a viral assemblage, the ephemeral and persistent taxa were closely related, implying ecological dynamics within these clusters. Furthermore, these dynamics occurred in both the RNA and DNA viral assemblages surveyed, implying that this structure is common in natural viral assemblages.IMPORTANCE Viruses are major agents of microbial mortality in marine systems, yet little is known about changes in the composition of viral assemblages in relation to those of the microbial communities that they infect. Here, we sampled coastal seawater every 2 weeks for 1 year and used high-throughput sequencing of marker genes to follow changes in the composition of two groups of ecologically important viruses, as well as the communities of bacteria and protists that serve as their respective hosts. Different subsets of genetically related viruses dominated at different times. These results demonstrate that although the genetic composition of viral assemblages is highly dynamic temporally, for the most part the shuffling of genotypes occurs within a few clusters of phylogenetically related viruses. Thus, it appears that even in temperate coastal waters with large seasonal changes, the highly dynamic shuffling of viral genotypes occurs largely within a few subsets of related individuals.

Interannual dynamics, diversity and evolution of the virome in Sclerotinia sclerotiorum from a single crop field

Mycovirus diversity is generally analyzed from isolates of fungal culture isolates at a single point in time as a snapshot. The stability of mycovirus composition within the same geographical location over time remains unclear. Not knowing how the population fluctuates in the field can be a source of unpredictability in the successful application of virocontrol. To better understand the changes over time, we monitored the interannual dynamics and abundance of mycoviruses infecting Sclerotinia sclerotiorum at a rapeseed-growing field for three years. We found that the virome in S. sclerotiorum harbors unique mycovirus compositions each year. In total, sixty-eight mycoviruses were identified, among which twenty-four were detected in all three successive years. These twenty-four mycoviruses can be classified as the members of the core virome in this S. sclerotiorum population, which show persistence and relatively high transmissibility under field conditions. Nearly two-thirds of the mycoviruses have positive-sense, single-stranded RNA genomes and were found consistently across all three years. Moreover, twenty-eight mycoviruses are newly described, including four novel, multi-segmented narnaviruses, and four unique bunyaviruses. Overall, the newly discovered mycoviruses in this study belong to as many as twenty families, into which eight were first identified in S. sclerotiorum, demonstrating evolutionarily diverse viromes. Our findings not only shed light on the annual variation of mycovirus diversity but also provide important virus evolutionary clues.

Seasonal impact of grazing, viral mortality, resource availability and light on the group-specific growth rates of coastal Mediterranean bacterioplankton

Estimation of prokaryotic growth rates is critical to understand the ecological role and contribution of different microbes to marine biogeochemical cycles. However, there is a general lack of knowledge on what factors control the growth rates of different prokaryotic groups and how these vary between sites and along seasons at a given site. We carried out several manipulation experiments during the four astronomical seasons in the coastal NW Mediterranean in order to evaluate the impact of grazing, viral mortality, resource competition and light on the growth and loss rates of prokaryotes. Gross and net growth rates of different bacterioplankton groups targeted by group-specific CARD-FISH probes and infrared microscopy (for aerobic anoxygenic phototrophs, AAP), were calculated from changes in cell abundances. Maximal group-specific growth rates were achieved when both predation pressure and nutrient limitation were experimentally minimized, while only a minimal effect of viral pressure on growth rates was observed; nevertheless, the response to predation removal was more remarkable in winter, when the bacterial community was not subjected to nutrient limitation. Although all groups showed increases in their growth rates when resource competition as well as grazers and viral pressure were reduced, Alteromonadaceae consistently presented the highest rates in all seasons. The response to light availability was generally weaker than that to the other factors, but it was variable between seasons. In summer and spring, the growth rates of AAP were stimulated by light whereas the growth of the SAR11 clade (likely containing proteorhodopsin) was enhanced by light in all seasons. Overall, our results set thresholds on bacterioplankton group-specific growth and mortality rates and contribute to estimate the seasonally changing contribution of various bacterioplankton groups to the function of microbial communities. Our results also indicate that the least abundant groups display the highest growth rates, contributing to the recycling of organic matter to a much greater extent than what their abundances alone would predict.

Seasonal and diel patterns of abundance and activity of viruses in the Red Sea

Virus-microbe interactions have been studied in great molecular details for many years in cultured model systems, yielding a plethora of knowledge on how viruses use and manipulate host machinery. Since the advent of molecular techniques and high-throughput sequencing, methods such as cooccurrence, nucleotide composition, and other statistical frameworks have been widely used to infer virus-microbe interactions, overcoming the limitations of culturing methods. However, their accuracy and relevance is still debatable as cooccurrence does not necessarily mean interaction. Here we introduce an ecological perspective of marine viral communities and potential interaction with their hosts, using analyses that make no prior assumptions on specific virus-host pairs. By size fractionating water samples into free viruses and microbes (i.e., also viruses inside or attached to their hosts) and looking at how viral group abundance changes over time along both fractions, we show that the viral community is undergoing a change in rank abundance across seasons, suggesting a seasonal succession of viruses in the Red Sea. We use abundance patterns in the different size fractions to classify viral clusters, indicating potential diverse interactions with their hosts and potential differences in life history traits between major viral groups. Finally, we show hourly resolved variations of intracellular abundance of similar viral groups, which might indicate differences in their infection cycles or metabolic capacities.

T4-like myovirus community shaped by dispersal and deterministic processes in the South China Sea

As the most abundant and genetically diverse biological entities, viruses significantly influence ecological, biogeographical and evolutionary processes in the ocean. However, the biogeography of marine viruses and the drivers shaping viral community are unclear. Here, the biogeographic patterns of T4-like viruses and the relative impacts of deterministic (environmental selection) and dispersal (spatial distance) processes were investigated in the northern South China Sea. The dominant viral operational taxonomic units were affiliated with previously defined Marine, Estuary, Lake and Paddy Groups. A clear viral biogeographic pattern was observed along the environmental gradient from the estuary to open sea. Marine Groups I and IV had a wide geographical distribution, whereas Marine Groups II, III and V were abundant in lower-salinity continental or eutrophic environments. A significant distance-decay pattern was noted for the T4-like viral community, especially for those infecting cyanobacteria. Both deterministic and dispersal processes influenced viral community assembly, although environmental selection (e.g. temperature, salinity, bacterial abundance and community, etc.) had a greater impact than spatial distance. Network analysis confirmed the strong association between viral and bacterial community composition, and suggested a diverse ecological relationship (e.g. lysis, co-infection or mutualistic) between and within viruses and their potential bacterial hosts.

Double-stranded DNA virioplankton dynamics and reproductive strategies in the oligotrophic open ocean water column

Microbial communities are critical to ecosystem dynamics and biogeochemical cycling in the open oceans. Viruses are essential elements of these communities, influencing the productivity, diversity, and evolution of cellular hosts. To further explore the natural history and ecology of open-ocean viruses, we surveyed the spatiotemporal dynamics of double-stranded DNA (dsDNA) viruses in both virioplankton and bacterioplankton size fractions in the North Pacific Subtropical Gyre, one of the largest biomes on the planet. Assembly and clustering of viral genomes revealed a peak in virioplankton diversity at the base of the euphotic zone, where virus populations and host species richness both reached their maxima. Simultaneous characterization of both extracellular and intracellular viruses suggested depth-specific reproductive strategies. In particular, analyses indicated elevated lytic interactions in the mixed layer, more temporally variable temperate phage interactions at the base of the euphotic zone, and increased lysogeny in the mesopelagic ocean. Furthermore, the depth variability of auxiliary metabolic genes suggested habitat-specific strategies for viral influence on light-energy, nitrogen, and phosphorus acquisition during host infection. Most virus populations were temporally persistent over several years in this environment at the 95% nucleic acid identity level. In total, our analyses revealed variable distributional patterns and diverse reproductive and metabolic strategies of virus populations in the open-ocean water column.

Long-term stability and Red Queen-like strain dynamics in marine viruses

Viruses that infect microorganisms dominate marine microbial communities numerically, with impacts ranging from host evolution to global biogeochemical cycles^1,2. However, virus community dynamics, necessary for conceptual and mechanistic model development, remains difficult to assess. Here, we describe the long-term stability of a viral community by analysing the metagenomes of near-surface 0.02-0.2 μm samples from the San Pedro Ocean Time-series³ that were sampled monthly over 5 years. Of 19,907 assembled viral contigs (>5 kb, mean 15 kb), 97% were found in each sample (by >98% ID metagenomic read recruitment) to have relative abundances that ranged over seven orders of magnitude, with limited temporal reordering of rank abundances along with little change in richness. Seasonal variations in viral community composition were superimposed on the overall stability; maximum community similarity occurred at 12-month intervals. Despite the stability of viral genotypic clusters that had 98% sequence identity, viral sequences showed transient variations in single-nucleotide polymorphisms (SNPs) and constant turnover of minor population variants, each rising and falling over a few months, reminiscent of Red Queen dynamics⁴. The rise and fall of variants within populations, interpreted through the perspective of known virus-host interactions⁵, is consistent with the hypothesis that fluctuating selection acts on a microdiverse cloud of strains, and this succession is associated with ever-shifting virus-host defences and counterdefences. This results in long-term virus-host coexistence that is facilitated by perpetually changing minor variants.

Multi-year dynamics of fine-scale marine cyanobacterial populations are more strongly explained by phage interactions than abiotic, bottom-up factors

Currently defined ecotypes in marine cyanobacteria Prochlorococcus and Synechococcus likely contain subpopulations that themselves are ecologically distinct. We developed and applied high-throughput sequencing for the 16S-23S rRNA internally transcribed spacer (ITS) to examine ecotype and fine-scale genotypic community dynamics for monthly surface water samples spanning 5 years at the San Pedro Ocean Time-series site. Ecotype-level structure displayed regular seasonal patterns including succession, consistent with strong forcing by seasonally varying abiotic parameters (e.g. temperature, nutrients, light). We identified tens to thousands of amplicon sequence variants (ASVs) within ecotypes, many of which exhibited distinct patterns over time, suggesting ecologically distinct populations within ecotypes. Community structure within some ecotypes exhibited regular, seasonal patterns, but not for others, indicating other more irregular processes such as phage interactions are important. Network analysis including T4-like phage genotypic data revealed distinct viral variants correlated with different groups of cyanobacterial ASVs including time-lagged predator-prey relationships. Variation partitioning analysis indicated that phage community structure more strongly explains cyanobacterial community structure at the ASV level than the abiotic environmental factors. These results support a hierarchical model whereby abiotic environmental factors more strongly shape niche partitioning at the broader ecotype level while phage interactions are more important in shaping community structure of fine-scale variants within ecotypes.

Dynamic marine viral infections and major contribution to photosynthetic processes shown by spatiotemporal picoplankton metatranscriptomes

Viruses provide top-down control on microbial communities, yet their direct study in natural environments was hindered by culture limitations. The advance of bioinformatics enables cultivation-independent study of viruses. Many studies assemble new viral genomes and study viral diversity using marker genes from free viruses. Here we use cellular metatranscriptomics to study active community-wide viral infections. Recruitment to viral contigs allows tracking infection dynamics over time and space. Our assemblies represent viral populations, but appear biased towards low diversity viral taxa. Tracking relatives of published T4-like cyanophages and pelagiphages reveals high genomic continuity. We determine potential hosts by matching dynamics of infection with abundance of particular microbial taxa. Finally, we quantify the relative contribution of cyanobacteria and viruses to photosystem-II psbA (reaction center) expression in our study sites. We show sometimes >50% of all cyanobacterial+viral psbA expression is of viral origin, highlighting the contribution of viruses to photosynthesis and oxygen production.

Here, Sieradzki et al. use metatranscriptomics to study active community-wide viral infections at three coastal California sites throughout a year, identify potential viral hosts, and show that viruses can contribute a substantial amount to photosystem-II psbA expression.

The Environment and Cyanophage Diversity: Insights From Environmental Sequencing of DNA Polymerase

Globally distributed and abundant cyanophages in the family Myoviridae have dsDNA genomes with variable gene content, including host-derived auxiliary metabolic genes (AMGs) that potentially can facilitate viral replication. However, it is not well understood how this variation in gene content interacts with environmental variables to shape cyanomyovirus communities. This project correlated the genetic repertoire of cyanomyoviruses with their phyologeny, and investigated cyanomyovirus ecotype distribution as a function of environmental conditions across locations and seasons. Reference cyanomyovirus genomes were compared for their overlap in gene content to infer phyologenetic distances, and these distances were compared to distances calculated based on DNA polymerase (gp43) gene sequences. In turn, gp43 partial gene sequences amplified from natural cyanophage communities were used to describe cyanomyovirus community composition and to assess the relationship between environmental variables. The results showed the following: (1) DNA polymerase gene phylogeny generally correlated with the similarity in gene content among reference cyanomyoviruses, and thus can be used to describe environmental cyanomyovirus communities; (2) spatial and seasonal patterns in cyanomyovirus communities were related to environmental variables; (3) salinity and temperature, combined with nutrient concentration were predictors of cyanomyovirus richness, diversity and community composition. This study shows that environmental variables shape viral communities by drawing on a diverse seed bank of viral genotypes. From these results it is evident that that viral ecotypes with their corresponding genetic repertoires underlie selection pressures. However, the mechanisms involved in selecting for specific viral genotypes remain to be fully understood.

Family A DNA Polymerase Phylogeny Uncovers Diversity and Replication Gene Organization in the Virioplankton

Shotgun metagenomics, which allows for broad sampling of viral diversity, has uncovered genes that are widely distributed among virioplankton populations and show linkages to important biological features of unknown viruses. Over 25% of known dsDNA phage carry the DNA polymerase I (polA) gene, making it one of the most widely distributed phage genes. Because of its pivotal role in DNA replication, this enzyme is linked to phage lifecycle characteristics. Previous research has suggested that a single amino acid substitution might be predictive of viral lifestyle. In this study Chesapeake Bay virioplankton were sampled by shotgun metagenomic sequencing (using long and short read technologies). More polA sequences were predicted from this single viral metagenome (virome) than from 86 globally distributed virome libraries (ca. 2,100, and 1,200, respectively). The PolA peptides predicted from the Chesapeake Bay virome clustered with 69% of PolA peptides from global viromes; thus, remarkably the Chesapeake Bay virome captured the majority of known PolA peptide diversity in viruses. This deeply sequenced virome also expanded the diversity of PolA sequences, increasing the number of PolA clusters by 44%. Contigs containing polA sequences were also used to examine relationships between phylogenetic clades of PolA and other genes within unknown viral populations. Phylogenic analysis revealed five distinct groups of phages distinguished by the amino acids at their 762 (Escherichia coli IAI39 numbering) positions and replication genes. DNA polymerase I sequences from Tyr762 and Phe762 groups were most often neighbored by ring-shaped superfamily IV helicases and ribonucleotide reductases (RNRs). The Leu762 groups had non-ring shaped helicases from superfamily II and were further distinguished by an additional helicase gene from superfamily I and the lack of any identifiable RNR genes. Moreover, we found that the inclusion of ribonucleotide reductase associated with PolA helped to further differentiate phage diversity, chiefly within lytic podovirus populations. Altogether, these data show that DNA Polymerase I is a useful marker for observing the diversity and composition of the virioplankton and may be a driving factor in the divergence of phage replication components.

Selective Maintenance of Multiple CRISPR Arrays Across Prokaryotes

Prokaryotes are under nearly constant attack by viral pathogens. To protect against this threat of infection, bacteria and archaea have evolved a wide array of defense mechanisms, singly and in combination. While immune diversity in a single organism likely reduces the chance of pathogen evolutionary escape, it remains puzzling why many prokaryotes also have multiple, seemingly redundant, copies of the same type of immune system. Here, we focus on the highly flexible CRISPR adaptive immune system, which is present in multiple copies in a surprising 28% of the prokaryotic genomes in RefSeq. We use a comparative genomics approach looking across all prokaryotes to demonstrate that on average, organisms are under selection to maintain more than one CRISPR array. Given this surprising conclusion, we consider several hypotheses concerning the source of selection and include a theoretical analysis of the possibility that a trade-off between memory span and learning speed could select for both "long-term memory" and "short-term memory" CRISPR arrays.

Phage puppet masters of the marine microbial realm

Viruses numerically dominate our oceans; however, we have only just begun to document the diversity, host range and infection dynamics of marine viruses, as well as the subsequent effects of infection on both host cell metabolism and oceanic biogeochemistry. Bacteriophages (that is, phages: viruses that infect bacteria) are highly abundant and are known to play critical roles in bacterial mortality, biogeochemical cycling and horizontal gene transfer. This Review Article summarizes current knowledge of marine viral ecology and highlights the importance of phage particles to the dissolved organic matter pool, as well as the complex interactions between phages and their bacterial hosts. We emphasize the newly recognized roles of phages as puppet masters of their bacterial hosts, where phages are capable of altering the metabolism of infected bacteria through the expression of auxiliary metabolic genes and the redirection of host gene expression patterns. Finally, we propose the 'royal family model' as a hypothesis to describe successional patterns of bacteria and phages over time in marine systems, where despite high richness and significant seasonal differences, only a small number of phages appear to continually dominate a given marine ecosystem. Although further testing is required, this model provides a framework for assessing the specificity and ecological consequences of phage-host dynamics.

High Diversity of Myocyanophage in Various Aquatic Environments Revealed by High-Throughput Sequencing of Major Capsid Protein Gene With a New Set of Primers

Myocyanophages, a group of viruses infecting cyanobacteria, are abundant and play important roles in elemental cycling. Here we investigated the particle-associated viral communities retained on 0.2 μm filters and in sediment samples (representing ancient cyanophage communities) from four ocean and three lake locations, using high-throughput sequencing and a newly designed primer pair targeting a gene fragment (∼145-bp in length) encoding the cyanophage gp23 major capsid protein (MCP). Diverse viral communities were detected in all samples. The fragments of 142-, 145-, and 148-bp in length were most abundant in the amplicons, and most sequences (>92%) belonged to cyanophages. Additionally, different sequencing depths resulted in different diversity estimates of the viral community. Operational taxonomic units obtained from deep sequencing of the MCP gene covered the majority of those obtained from shallow sequencing, suggesting that deep sequencing exhibited a more complete picture of cyanophage community than shallow sequencing. Our results also revealed a wide geographic distribution of marine myocyanophages, i.e., higher dissimilarities of the myocyanophage communities corresponded with the larger distances between the sampling sites. Collectively, this study suggests that the newly designed primer pair can be effectively used to study the community and diversity of myocyanophage from different environments, and the high-throughput sequencing represents a good method to understand viral diversity.

Assessing the diversity of the g23 gene of T4-like bacteriophages from Lake Baikal with high-throughput sequencing

Based on second generation sequencing (MiSeq platform, Illumina), we determined the genetic diversity of T4-like bacteriophages of the family Myoviridae by analysing fragments of the major capsid protein gene g23 in the plankton of Lake Baikal. The sampling depth in our study was significantly higher than in those obtained by the Sanger method before. We obtained 33 701 sequences of the g23 gene fragments, 141 operational taxonomic units (OTUs) of which were identified. 86 OTUs (60.9%) had the closest relatives from lakes Bourget and Annecy, and 28 OTUs (19.8%) had the highest identity with the Baikal g23 clones, which had been previously identified in the northern and southern basins of the lake by the Sanger method. The remaining OTUs were similar to the clones from other ecosystems. We showed a high genetic diversity of T4-type bacteriophages and a genetic difference with the phage communities from other ecosystems.

Temporal dynamics of uncultured viruses: a new dimension in viral diversity

Recent work has vastly expanded the known viral genomic sequence space, but the seasonal dynamics of viral populations at the genome level remain unexplored. Here we followed the viral community in a freshwater lake for 1 year using genome-resolved viral metagenomics, combined with detailed analyses of the viral community structure, associated bacterial populations and environmental variables. We reconstructed 8950 complete and partial viral genomes, the majority of which were not persistent in the lake throughout the year, but instead continuously succeeded each other. Temporal analysis of 732 viral genus-level clusters demonstrated that one-fifth were undetectable at specific periods of the year. Based on host predictions for a subset of reconstructed viral genomes, we for the first time reveal three distinct patterns of host-pathogen dynamics, where the viruses may peak before, during or after the peak in their host's abundance, providing new possibilities for modelling of their interactions. Time series metagenomics opens up a new dimension in viral profiling, which is essential to understand the full scale of viral diversity and evolution, and the ecological roles of these important factors in the global ecosystem.

Benchmarking viromics: an in silico evaluation of metagenome-enabled estimates of viral community composition and diversity

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.

Marine archaeal dynamics and interactions with the microbial community over 5 years from surface to seafloor

Marine archaea are critical contributors to global carbon and nitrogen redox cycles, but their temporal variability and microbial associations across the water column are poorly known. We evaluated seasonal variability of free living (0.2-1 μm size fraction) Thaumarchaea Marine Group I (MGI) and Euryarchaea Marine Group II (MGII) communities and their associations with the microbial community from surface to seafloor (890 m) over 5 years by 16S rRNA V4-V5 gene sequencing. MGI and MGII communities demonstrated distinct compositions at different depths, and seasonality at all depths. Microbial association networks at 150 m, 500 m and 890 m, revealed diverse assemblages of MGI (presumed ammonia oxidizers) and Nitrospina taxa (presumed dominant nitrite oxidizers, completing the nitrification process), suggesting distinct MGI-Nitrospina OTUs are responsible for nitrification at different depths and seasons, and depth- related and seasonal variability in nitrification could be affected by alternating MGI-Nitrospina assemblages. MGII taxa also showed distinct correlations to possibly heterotrophic bacteria, most commonly to members of Marine Group A, Chloroflexi, Marine Group B, and SAR86. Thus, both MGI and MGII likely have dynamic associations with bacteria based on similarities in activity or other interactions that select for distinct microbial assemblages over time. The importance of MGII taxa as members of the heterotrophic community previously reported for photic zone appears to apply throughout the water column.

Ecological dynamics and co-occurrence among marine phytoplankton, bacteria and myoviruses shows microdiversity matters

Numerous ecological processes, such as bacteriophage infection and phytoplankton-bacterial interactions, often occur via strain-specific mechanisms. Therefore, studying the causes of microbial dynamics should benefit from highly resolving taxonomic characterizations. We sampled daily to weekly over 5 months following a phytoplankton bloom off Southern California and examined the extent of microdiversity, that is, significant variation within 99% sequence similarity clusters, operational taxonomic units (OTUs), of bacteria, archaea, phytoplankton chloroplasts (all via 16S or intergenic spacer (ITS) sequences) and T4-like-myoviruses (via g23 major capsid protein gene sequence). The extent of microdiversity varied between genes (ITS most, g23 least) and only temporally common taxa were highly microdiverse. Overall, 60% of taxa exhibited microdiversity; 59% of these had subtypes that changed significantly as a proportion of the parent taxon, indicating ecologically distinct taxa. Pairwise correlations between prokaryotes and myoviruses or phytoplankton (for example, highly microdiverse Chrysochromulina sp.) improved when using single-base variants. Correlations between myoviruses and SAR11 increased in number (172 vs 9, Spearman>0.65) and became stronger (0.61 vs 0.58, t-test: P<0.001) when using SAR11 ITS single-base variants vs OTUs. Whole-community correlation between SAR11 and myoviruses was much improved when using ITS single-base variants vs OTUs, with Mantel rho=0.49 vs 0.27; these results are consistent with strain-specific interactions. Mantel correlations suggested >1 μm (attached/large) prokaryotes are a major myovirus source. Consideration of microdiversity improved observation of apparent host and virus networks, and provided insights into the ecological and evolutionary factors influencing the success of lineages, with important implications to ecosystem resilience and microbial function.

A comprehensive method for amplicon-based and metagenomic characterization of viruses, bacteria, and eukaryotes in freshwater samples

BACKGROUND

Studies of environmental microbiota typically target only specific groups of microorganisms, with most focusing on bacteria through taxonomic classification of 16S rRNA gene sequences. For a more holistic understanding of a microbiome, a strategy to characterize the viral, bacterial, and eukaryotic components is necessary.

RESULTS

We developed a method for metagenomic and amplicon-based analysis of freshwater samples involving the concentration and size-based separation of eukaryotic, bacterial, and viral fractions. Next-generation sequencing and culture-independent approaches were used to describe and quantify microbial communities in watersheds with different land use in British Columbia. Deep amplicon sequencing was used to investigate the distribution of certain viruses (g23 and RdRp), bacteria (16S rRNA and cpn60), and eukaryotes (18S rRNA and ITS). Metagenomic sequencing was used to further characterize the gene content of the bacterial and viral fractions at both taxonomic and functional levels.

CONCLUSION

This study provides a systematic approach to separate and characterize eukaryotic-, bacterial-, and viral-sized particles. Methodologies described in this research have been applied in temporal and spatial studies to study the impact of land use on watershed microbiomes in British Columbia.

Biogeographic Variation in Host Range Phenotypes and Taxonomic Composition of Marine Cyanophage Isolates

Despite the important role of phages in marine systems, little is understood about how their diversity is distributed in space. Biogeographic patterns of marine phages may be difficult to detect due to their vast genetic diversity, which may not be accurately represented by conserved marker genes. To investigate the spatial biogeographic structure of marine phages, we isolated over 400 cyanophages on Synechococcus host strain WH7803 at three coastal locations in the United States (Rhode Island, Washington, and southern California). Approximately 90% of the cyanophage isolates were myoviruses, while the other 10% were podoviruses. The diversity of isolates was further characterized in two ways: (i) taxonomically, using conserved marker genes and (ii) phenotypically, by testing isolates for their ability to infect a suite of hosts, or their "host range." Because host range is a highly variable trait even among closely related isolates, we hypothesized that host range phenotypes of cyanophage isolates would vary more strongly among locations than would taxonomic composition. Instead, we found evidence for strong biogeographic variation both in taxonomic composition and host range phenotypes, with little taxonomic overlap among the three coastal regions. For both taxonomic composition and host range phenotypes, cyanophage communities from California and Rhode Island were the most dissimilar, while Washington communities exhibited similarity to each of the other two locations. These results suggest that selection imposed by spatial variation in host dynamics influence the biogeographic distribution of cyanophages.

Biogeography of bacteriophages at four hydrothermal vent sites in the Antarctic based on g23 sequence diversity

In this study, which was carried out within the ChEsSO consortium project (Chemosynthetically driven ecosystems south of the Polar Front), we sampled two hydrothermal vent sites on the East Scotia Ridge, Scotia Sea, one in the Kemp Caldera, South Sandwich Arc and one in the Bransfield Strait, north-west of the Antarctic Peninsula, which exhibit strong differences in their chemical characteristics. We compared a subset of their bacteriophage population by Sanger- and 454-sequencing of g23, which codes for the major capsid protein of T4likeviruses. We found that the sites differ vastly in their bacteriophage diversity, which reflects the differences in the chemical conditions and therefore putatively the differences in microbial hosts living at these sites. Comparing phage diversity in the vent samples to other aquatic samples, the vent samples formed a distinct separate cluster, which also included the non-vent control samples that were taken several hundred meters above the vent chimneys. This indicates that the influence of the vents on the microbial population and therefore also the bacteriophage population extends much further than anticipated.

Dynamics of Viral Abundance and Diversity in a Sphagnum-Dominated Peatland: Temporal Fluctuations Prevail Over Habitat

Viruses impact microbial activity and carbon cycling in various environments, but their diversity and ecological importance in Sphagnum-peatlands are unknown. Abundances of viral particles and prokaryotes were monitored bi-monthly at a fen and a bog at two different layers of the peat surface. Viral particle abundance ranged from 1.7 x 10⁶ to 5.6 x 10⁸ particles mL^-1, and did not differ between fen and bog but showed seasonal fluctuations. These fluctuations were positively correlated with prokaryote abundance and dissolved organic carbon, and negatively correlated with water-table height and dissolved oxygen. Using shotgun metagenomics we observed a shift in viral diversity between winter/spring and summer/autumn, indicating a seasonal succession of viral communities, mainly driven by weather-related environmental changes. Based on the seasonal asynchrony between viral and microbial diversity, we hypothesize a seasonal shift in the active microbial communities associated with a shift from lysogenic to lytic lifestyles. Our results suggest that temporal variations of environmental conditions rather than current habitat differences control the dynamics of virus-host interactions in Sphagnum-dominated peatlands.

Biogeography of Viruses in the Sea

Viral ecology is a rapidly progressing area of research, as molecular methods have improved significantly for targeted research on specific populations and whole communities. To interpret and synthesize global viral diversity and distribution, it is feasible to assess whether macroecology concepts can apply to marine viruses. We review how viral and host life history and physical properties can influence viral distribution in light of biogeography and metacommunity ecology paradigms. We highlight analytical approaches that can be applied to emerging global data sets and meta-analyses to identify individual taxa with global influence and drivers of emergent properties that influence microbial community structure by drawing on examples across the spectrum of viral taxa, from RNA to ssDNA and dsDNA viruses.

Using signature genes as tools to assess environmental viral ecology and diversity

Viruses (including bacteriophages) are the most abundant biological entities on the planet. As such, they are thought to have a major impact on all aspects of microbial community structure and function. Despite this critical role in ecosystem processes, the study of virus/phage diversity has lagged far behind parallel studies of the bacterial and eukaryotic kingdoms, largely due to the absence of any universal phylogenetic marker. Here we review the development and use of signature genes to investigate viral diversity, as a viable strategy for data sets of specific virus groups. Genes that have been used include those encoding structural proteins, such as portal protein, major capsid protein, and tail sheath protein, auxiliary metabolism genes, such as psbA, psbB,and phoH, and several polymerase genes. These marker genes have been used in combination with PCR-based fingerprinting and/or sequencing strategies to investigate spatial, temporal, and seasonal variations and diversity in a wide range of habitats.

Deep sequencing of the viral phoH gene reveals temporal variation, depth-specific composition, and persistent dominance of the same viral phoH genes in the Sargasso Sea

Deep sequencing of the viral phoH gene, a host-derived auxiliary metabolic gene, was used to track viral diversity throughout the water column at the Bermuda Atlantic Time-series Study (BATS) site in the summer (September) and winter (March) of three years. Viral phoH sequences reveal differences in the viral communities throughout a depth profile and between seasons in the same year. Variation was also detected between the same seasons in subsequent years, though these differences were not as great as the summer/winter distinctions. Over 3,600 phoH operational taxonomic units (OTUs; 97% sequence identity) were identified. Despite high richness, most phoH sequences belong to a few large, common OTUs whereas the majority of the OTUs are small and rare. While many OTUs make sporadic appearances at just a few times or depths, a small number of OTUs dominate the community throughout the seasons, depths, and years.

Marine microbial community dynamics and their ecological interpretation

Recent advances in studying the dynamics of marine microbial communities have shown that the composition of these communities follows predictable patterns and involves complex network interactions, which shed light on the underlying processes regulating these globally important organisms. Such 'holistic' (or organism- and system-based) studies of these communities complement popular reductionist, often culture-based, approaches for understanding organism function one gene or protein at a time. In this Review, we summarize our current understanding of marine microbial community dynamics at various scales, from hours to decades. We also explain how the data illustrate community resilience and seasonality, and reveal interactions among microorganisms.

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.

High temporal and spatial diversity in marine RNA viruses implies that they have an important role in mortality and structuring plankton communities

Viruses in the order Picornavirales infect eukaryotes, and are widely distributed in coastal waters. Amplicon deep-sequencing of the RNA dependent RNA polymerase (RdRp) revealed diverse and highly uneven communities of picorna-like viruses in the coastal waters of British Columbia (BC), Canada. Almost 300 000 pyrosequence reads revealed 145 operational taxonomic units (OTUs) based on 95% sequence similarity at the amino-acid level. Each sample had between 24 and 71 OTUs and there was little overlap among samples. Phylogenetic analysis revealed that some clades of OTUs were only found at one site; whereas, other clades included OTUs from all sites. Since most of these OTUs are likely from viruses that infect eukaryotic phytoplankton, and viral isolates infecting phytoplankton are strain-specific; each OTU probably arose from the lysis of a specific phytoplankton taxon. Moreover, the patchiness in OTU distribution, and the high turnover of viruses in the mixed layer, implies continuous infection and lysis by RNA viruses of a diverse array of eukaryotic phytoplankton taxa. Hence, these viruses are likely important elements structuring the phytoplankton community, and play a significant role in nutrient cycling and energy transfer.

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.

Top-down controls on bacterial community structure: microbial network analysis of bacteria, T4-like viruses and protists

Characterizing ecological relationships between viruses, bacteria and protists in the ocean are critical to understanding ecosystem function, yet these relationships are infrequently investigated together. We evaluated these relationships through microbial association network analysis of samples collected approximately monthly from March 2008 to January 2011 in the surface ocean (0-5 m) at the San Pedro Ocean Time series station. Bacterial, T4-like myoviral and protistan communities were described by Automated Ribosomal Intergenic Spacer Analysis and terminal restriction fragment length polymorphism of the gene encoding the major capsid protein (g23) and 18S ribosomal DNA, respectively. Concurrent shifts in community structure suggested similar timing of responses to environmental and biological parameters. We linked T4-like myoviral, bacterial and protistan operational taxonomic units by local similarity correlations, which were then visualized as association networks. Network links (correlations) potentially represent synergistic and antagonistic relationships such as viral lysis, grazing, competition or other interactions. We found that virus-bacteria relationships were more cross-linked than protist-bacteria relationships, suggestive of increased taxonomic specificity in virus-bacteria relationships. We also found that 80% of bacterial-protist and 74% of bacterial-viral correlations were positive, with the latter suggesting that at monthly and seasonal timescales, viruses may be following their hosts more often than controlling host abundance.

New approaches indicate constant viral diversity despite shifts in assemblage structure in an Australian hypersaline lake

It is widely stated that viruses represent the most significant source of biodiversity on Earth, yet characterizing the diversity of viral assemblages in natural systems remains difficult. Viral diversity studies are challenging because viruses lack universally present, phylogenetically informative genes. Here, we developed an approach to estimate viral diversity using a series of functional and novel conserved genes. This approach provides direct estimates of viral assemblage diversity while retaining resolution at the level of individual viral populations in a natural system. We characterized viral assemblages in eight samples from hypersaline Lake Tyrrell (LT), Victoria, Australia, using 39,636 viral contigs. We defined viral operational taxonomic units (OTUs) in two ways. First, we used genes with three different functional predictions that were abundantly represented in the data set. Second, we clustered proteins of unknown function based on sequence similarity, and we chose genes represented by three clusters with numerous members to define OTUs. In combination, diversity metrics indicated between 412 and 735 sampled populations, and the number of populations remained relatively constant across samples. We determined the relative representation of each viral OTU in each sample and found that viral assemblage structures correlate with salinity and solution chemistry. LT viral assemblages were near-replicates from the same site sampled a few days apart but differed significantly on other spatial and temporal scales. The OTU definition approach proposed here paves the way for metagenomics-based analyses of viral assemblages using ecological models previously applied to bacteria and archaea.

Strong seasonality and interannual recurrence in marine myovirus communities

The temporal community dynamics and persistence of different viral types in the marine environment are still mostly obscure. Polymorphism of the major capsid protein gene, g23, was used to investigate the community composition dynamics of T4-like myoviruses in a North Atlantic fjord for a period of 2 years. A total of 160 unique operational taxonomic units (OTUs) were identified by terminal restriction fragment length polymorphism (TRFLP) of the gene g23. Three major community profiles were identified (winter-spring, summer, and autumn), which resulted in a clear seasonal succession pattern. These seasonal transitions were recurrent over the 2 years and significantly correlated with progression of seawater temperature, Synechococcus abundance, and turbidity. The appearance of the autumn viral communities was concomitant with the occurrence of prominent Synechococcus blooms. As a whole, we found a highly dynamic T4-like viral community with strong seasonality and recurrence patterns. These communities were unexpectedly dominated by a group of persistently abundant viruses.

Temporal variability and coherence of euphotic zone bacterial communities over a decade in the Southern California Bight

Time-series are critical to understanding long-term natural variability in the oceans. Bacterial communities in the euphotic zone were investigated for over a decade at the San Pedro Ocean Time-series station (SPOT) off southern California. Community composition was assessed by Automated Ribosomal Intergenic Spacer Analysis (ARISA) and coupled with measurements of oceanographic parameters for the surface ocean (0-5 m) and deep chlorophyll maximum (DCM, average depth ≈ 30 m). SAR11 and cyanobacterial ecotypes comprised typically more than one-third of the measured community; diversity within both was temporally variable, although a few operational taxonomic units (OTUs) were consistently more abundant. Persistent OTUs, mostly Alphaproteobacteria (SAR11 clade), Actinobacteria and Flavobacteria, tended to be abundant, in contrast to many rarer yet intermittent and ephemeral OTUs. Association networks revealed potential niches for key OTUs from SAR11, cyanobacteria, SAR86 and other common clades on the basis of robust correlations. Resilience was evident by the average communities drifting only slightly as years passed. Average Bray-Curtis similarity between any pair of dates was ≈ 40%, with a slight decrease over the decade and obvious near-surface seasonality; communities 8-10 years apart were slightly more different than those 1-4 years apart with the highest rate of change at 0-5 m between communities <4 years apart. The surface exhibited more pronounced seasonality than the DCM. Inter-depth Bray-Curtis similarities repeatedly decreased as the water column stratified each summer. Environmental factors were better predictors of shifts in community composition than months or elapsed time alone; yet, the best predictor was community composition at the other depth (that is, 0-5 m versus DCM).

Short-term observations of marine bacterial and viral communities: patterns, connections and resilience

Observation of short-term temporal variation in bacterial and viral communities is important for understanding patterns of aquatic microbial diversity. We collected surface seawater once daily for 38 consecutive days with seven more samples interspersed over 40 more days at one location ∼2 km from Santa Catalina Island, California. Bacterial communities were analyzed by automated ribosomal intergenic spacer analysis (ARISA) and viral communities were analyzed by terminal restriction fragment length polymorphism (TRFLP) of the conserved T4-like myoviral gene encoding the major capsid protein (g23). Common bacterial and viral taxa were consistently dominant, and relatively few displayed dramatic increases/decreases or 'boom/bust' patterns that might be expected from dynamic predator-prey interactions. Association network analysis showed most significant covariations (associations) occurred among bacterial taxa or among viral taxa and there were several modular (highly-interconnected) associations (P≤0.005). Associations observed between bacteria and viruses (P≤0.005) occurred with a median time lag of 2 days. Regression of all pairwise Bray-Curtis similarities between samples indicated a rate of bacterial community change that slows from 2.1%-0.18% per day over a week to 2 months; the rate stays around 0.4% per day for viruses. Our interpretation is that, over the scale of days, individual bacterial and viral OTUs can be dynamic and patterned; resulting in statistical associations regarded as potential ecological interactions. However, over the scale of weeks, average bacterial community variation is slower, suggesting that there is strong community-level ecological resilience, that is, a tendency to converge towards a 'mean' microbial community set by longer-term controlling factors.

Abundant SAR11 viruses in the ocean

Several reports proposed that the extraordinary dominance of the SAR11 bacterial clade in ocean ecosystems could be a consequence of unusual mechanisms of resistance to bacteriophage infection, including 'cryptic escape' through reduced cell size and/or K-strategist defence specialism. Alternatively, the evolution of high surface-to-volume ratios coupled with minimal genomes containing high-affinity transporters enables unusually efficient metabolism for oxidizing dissolved organic matter in the world's oceans that could support vast population sizes despite phage susceptibility. These ideas are important for understanding plankton ecology because they emphasize the potentially important role of top-down mechanisms in predation, thus determining the size of SAR11 populations and their concomitant role in biogeochemical cycling. Here we report the isolation of diverse SAR11 viruses belonging to two virus families in culture, for which we propose the name 'pelagiphage', after their host. Notably, the pelagiphage genomes were highly represented in marine viral metagenomes, demonstrating their importance in nature. One of the new phages, HTVC010P, represents a new podovirus subfamily more abundant than any seen previously, in all data sets tested, and may represent one of the most abundant virus subfamilies in the biosphere. This discovery disproves the theory that SAR11 cells are immune to viral predation and is consistent with the interpretation that the success of this highly abundant microbial clade is the result of successfully evolved adaptation to resource competition.

Marine cyanophages exhibit local and regional biogeography

Biogeographic patterns have been demonstrated for a wide range of microorganisms. Nevertheless, the biogeography of marine viruses has been slower to emerge. Here we investigate biogeographic patterns of marine cyanophages that infect Synechococcus sp. WH7803 across multiple spatial and temporal scales. We compared cyanophage myoviral communities from nine coastal sites in Southern New England (SNE), USA, one site in Long Island NY, and four sites from Bermuda's inshore waters by assaying cyanophage isolates using the myoviral g43 DNA polymerase gene. Cyanophage community composition varied temporally at each of the sites. Further, 6 years of sampling at one Narragansett Bay site revealed annual seasonal variations in community composition, driven by the seasonal reoccurrence of specific viral taxa. Although the four Bermuda communities were similar to one another, they were significantly different than the North American coastal communities, with almost no overlap of taxa between the two regions. Among the SNE sites, cyanophage community composition also varied significantly and was correlated with the body of water sampled (e.g. Narragansett Bay, Cape Cod Bay, Vineyard Sound), although here, the same viral taxa were found at multiple sites. This study demonstrates that marine cyanophages display striking seasonal and spatial biogeographic patterns.