The diversity of cyanomyovirus populations along a North-South Atlantic Ocean transect

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.

How might bacteriophages shape biological invasions?

Invasions by eukaryotes dependent on environmentally acquired bacterial mutualists are often limited by the ability of bacterial partners to survive and establish free-living populations. Focusing on the model legume-rhizobium mutualism, we apply invasion biology hypotheses to explain how bacteriophages can impact the competitiveness of introduced bacterial mutualists. Predicting how phage-bacteria interactions affect invading eukaryotic hosts requires knowing the eco-evolutionary constraints of introduced and native microbial communities, as well as their differences in abundance and diversity. By synthesizing research from invasion biology, as well as bacterial, viral, and community ecology, we create a conceptual framework for understanding and predicting how phages can affect biological invasions through their effects on bacterial mutualists.

Novel Freshwater Cyanophages Provide New Insights into Evolutionary Relationships between Freshwater and Marine Cyanophages

Cyanobacteria and cyanophages are present widely in both freshwater and marine environments. However, freshwater cyanophages remain unknown largely due to the small numbers of cyanophage isolates despite their ecological and environmental significance. In this study, we present the characterization of two novel lytic freshwater cyanophages isolated from a tropical inland lake in Singapore, namely, cyanopodovirus S-SRP01 and cyanomyovirus S-SRM01, infecting two different strains of Synechococcus spp. Functional annotation of S-SRP01 and S-SRM01 genomes revealed a high degree of homology with marine cyanophages. Phylogenetic trees of concatenated genes and whole-genome alignment provided further evidence that S-SRP01 is close evolutionarily to marine cyanopodoviruses, while S-SRM01 is evolutionarily close to marine cyanomyoviruses. Few genetic similarities between freshwater and marine cyanophages have been identified in previous studies. The isolation of S-SRP01 and S-SRM01 expand current knowledge on freshwater cyanophages infecting Synechococcus spp. Their high degree of gene sharing provides new insights into the evolutionary relationships between freshwater and marine cyanophages. This relatedness is further supported by the discovery of similar phenomenon from other freshwater viral metagenomes. IMPORTANCE This study expands the current knowledge on freshwater cyanophage isolates and cyanophage genetic diversity, indicating that freshwater and marine cyanophages infecting Synechococcus spp. may share close genetic similarity and evolutionary relationships.

Spatial patterns in phage-Rhizobium coevolutionary interactions across regions of common bean domestication

Bacteriophages play significant roles in the composition, diversity, and evolution of bacterial communities. Despite their importance, it remains unclear how phage diversity and phage-host interactions are spatially structured. Local adaptation may play a key role. Nitrogen-fixing symbiotic bacteria, known as rhizobia, have been shown to locally adapt to domesticated common bean at its Mesoamerican and Andean sites of origin. This may affect phage-rhizobium interactions. However, knowledge about the diversity and coevolution of phages with their respective Rhizobium populations is lacking. Here, through the study of four phage-Rhizobium communities in Mexico and Argentina, we show that both phage and host diversity is spatially structured. Cross-infection experiments demonstrated that phage infection rates were higher overall in sympatric rhizobia than in allopatric rhizobia except for one Argentinean community, indicating phage local adaptation and host maladaptation. Phage-host interactions were shaped by the genetic identity and geographic origin of both the phage and the host. The phages ranged from specialists to generalists, revealing a nested network of interactions. Our results suggest a key role of local adaptation to resident host bacterial communities in shaping the phage genetic and phenotypic composition, following a similar spatial pattern of diversity and coevolution to that in the host.

Rumen Virus Populations: Technological Advances Enhancing Current Understanding

The rumen contains a multi-kingdom, commensal microbiome, including protozoa, bacteria, archaea, fungi and viruses, which enables ruminant herbivores to ferment and utilize plant feedstuffs that would be otherwise indigestible. Within the rumen, virus populations are diverse and highly abundant, often out-numbering the microbial populations that they both predate on and co-exist with. To date the research effort devoted to understanding rumen-associated viral populations has been considerably less than that given to the other microbial populations, yet their contribution to maintaining microbial population balance, intra-ruminal microbial lysis, fiber breakdown, nutrient cycling and genetic transfer may be highly significant. This review follows the technological advances which have contributed to our current understanding of rumen viruses and drawing on knowledge from other environmental and animal-associated microbiomes, describes the known and potential roles and impacts viruses have on rumen function and speculates on the future directions of rumen viral research.

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.

Genetic Diversity and Cooccurrence Patterns of Marine Cyanopodoviruses and Picocyanobacteria

Picocyanobacteria Prochlorococcus and Synechococcus are abundant in the global oceans and subject to active viral infection. In this study, the genetic diversity of picocyanobacteria and the genetic diversity of cyanopodoviruses were synchronously investigated along water columns in the equatorial Indian Ocean and over a seasonal time course in the coastal Sanya Bay, South China Sea. Using the 16S-23S rRNA internal transcribed spacer (ITS)-based clone library and quantitative PCR (qPCR) analyses, the picocyanobacterial community composition and abundance were determined. Sanya Bay was dominated by clade II Synechococcus during all the seasons, and a typical population shift from high-light-adapted Prochlorococcus to low-light-adapted Prochlorococcus was found along the vertical profiles. Strikingly, the DNA polymerase gene sequences of cyanopodoviruses revealed a much greater genetic diversity than we expected. Nearly one-third of the phylogenetic groups were newly described here. No apparent seasonal pattern was observed for the Sanya Bay picocyanobacterial or cyanopodoviral communities. Different dominant cyanopodovirus lineages were identified for the coastal area, upper euphotic zone, and middle-to-lower euphotic zone of the open ocean. Diversity indices of both picocyanobacteria and cyanopodoviruses were highest in the middle euphotic zone and both were lower in the upper euphotic zone, reflecting a host-virus interaction. Cyanopodoviral communities differed significantly between the upper euphotic zone and the middle-to-lower euphotic zone, showing a vertical pattern similar to that of picocyanobacteria. However, in the surface waters of the open ocean, cyanopodoviruses exhibited no apparent biogeographic pattern, differing from picocyanobacteria. This study demonstrates correlated distribution patterns of picocyanobacteria and cyanopodoviruses, as well as the complex biogeography of cyanopodoviruses.IMPORTANCE Picocyanobacteria are highly diverse and abundant in the ocean and display remarkable global biogeography and a vertical distribution pattern. However, how the diversity and distribution of picocyanobacteria affect those of the viruses that infect them remains largely unknown. Here we synchronously analyzed the community structures of cyanopodoviruses and picocyanobacteria at spatial and temporal scales. Both spatial and temporal variations of cyanopodoviral communities can be linked to those of picocyanobacteria. The coastal area, upper euphotic zone, and middle-to-lower euphotic zone of the open ocean have distinct cyanopodoviral communities, showing horizontal and vertical variation patterns closely related to those of picocyanobacteria. These findings emphasize the driving force of host community in shaping the biogeographic structure of viruses. Our work provides important information for future assessments of the ecological roles of viruses and hosts for each other.

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.

Improved sensitivity of lateral flow assay using paper-based sample concentration technique

Lateral flow assays (LFAs) hold great promise for point-of-care testing, especially in resource-poor settings. However, the poor sensitivity of LFAs limits their widespread applications. To address this, we developed a novel device by integrating dialysis-based concentration method into LFAs. The device successfully achieved 10-fold signal enhancement in Human Immunodeficiency Virus (HIV) nucleic acid detection with a detection limit of 0.1 nM and 4-fold signal enhancement in myoglobin (MYO) detection with a detection limit of 1.56 ng/mL in less than 25 min. This simple, low-cost and portable integrated device holds great potential for highly sensitive detection of various target analytes for medical diagnostics, food safety analysis and environmental monitoring.

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.

Population dynamics of cyanomyovirus in a tropical eutrophic reservoir

Samples from three stations in Kranji Reservoir, Singapore (n = 21) were collected and analyzed for cyanomyovirus abundance and diversity. A total of 73 different g20 (viral capsid assembly protein genes) amino acid sequences were obtained from this study. A phylogenetic analysis revealed that the 73 segments were distributed in six major clusters (α to ζ), with four unique subclusters, which were identified as KRM-I, KRM-II, KRM-III, and KRM-IV. The cyanophage community in Kranji Reservoir exhibited a large degree of diversity; the clones obtained in this study showed similarities to those from many different environments, including oceans, lakes, bays, and paddy floodwater, as well as clones from paddy field soils. However, the sequences in this study were generally found to be more closely related to the g20 sequences of freshwaters and brackish waters than those from marine environments. The rarefaction curves and Chao 1 indices from this study showed that the diversity of the cyanomyovirus community was greater during the Inter-monsoon periods than the Southwest and Northeast Monsoons. A few seasonal changes in the taxa were observed: (i) Cluster ζ was absent during the Southwest Monsoon, and (ii) most of the samples fell into Group 3 in the PCA score plot during the Northeast Monsoon, and the fraction of Cluster ɛ increased significantly.

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.

Marine cyanophages demonstrate biogeographic patterns throughout the global ocean

Myoviruses and podoviruses that infect cyanobacteria are the two major groups of marine cyanophages, but little is known of how their phylogenetic lineages are distributed in different habitats. In this study, we analyzed the phylogenetic relationships of cyanopodoviruses and cyanomyoviruses based on the existing genomes. The 28 cyanomyoviruses were classified into four clusters (I to IV), and 19 of the 20 cyanopodoviruses were classified into two clusters, MPP-A and MPP-B, with four subclusters within cluster MPP-B. These genomes were used to recruit cyanophage-like fragments from microbial and viral metagenomes to estimate the relative abundances of these cyanophage lineages. Our results showed that cyanopodoviruses and cyanomyoviruses are both abundant in various marine environments and that clusters MPP-B, II and III appear to be the most dominant lineages. Cyanopodoviruses and cluster I and IV cyanomyoviruses exhibited habitat-related variability in their relative levels of abundance, while cluster II and III cyanomyoviruses appeared to be consistently dominant in various habitats. Multivariate analyses showed that reads that mapped to Synechococcus phages and Prochlorococcus phages had distinct distribution patterns that were significantly correlated to those of Synechococcus and Prochlorococcus, respectively. The Mantel test also revealed a strong correlation between the community compositions of cyanophages and picocyanobacteria. Given that cyanomyoviruses tend to have a broad host range and some can cross-infect Synechococcus and Prochlorococcus, while cyanopodoviruses are commonly host specific, the observation that their community compositions both correlated significantly with that of picocyanobacteria was unexpected. Although cyanomyoviruses and cyanopodoviruses differ in host specificity, their biogeographic distributions are likely both constrained by the picocyanobacterial community.

Phylogenetic distribution of the capsid assembly protein gene (g20) of cyanophages in paddy floodwaters in Northeast China

Numerous studies have revealed the high diversity of cyanophages in marine and freshwater environments, but little is currently known about the diversity of cyanophages in paddy fields, particularly in Northeast (NE) China. To elucidate the genetic diversity of cyanophages in paddy floodwaters in NE China, viral capsid assembly protein gene (g20) sequences from five floodwater samples were amplified with the primers CPS1 and CPS8. Denaturing gradient gel electrophoresis (DGGE) was applied to distinguish different g20 clones. In total, 54 clones differing in g20 nucleotide sequences were obtained in this study. Phylogenetic analysis showed that the distribution of g20 sequences in this study was different from that in Japanese paddy fields, and all the sequences were grouped into Clusters α, β, γ and ε. Within Clusters α and β, three new small clusters (PFW-VII∼-IX) were identified. UniFrac analysis of g20 clone assemblages demonstrated that the community compositions of cyanophage varied among marine, lake and paddy field environments. In paddy floodwater, community compositions of cyanophage were also different between NE China and Japan.

The genome, proteome and phylogenetic analysis of Sinorhizobium meliloti phage ΦM12, the founder of a new group of T4-superfamily phages

Phage ΦM12 is an important transducing phage of the nitrogen-fixing rhizobial bacterium Sinorhizobium meliloti. Here we report the genome, phylogenetic analysis, and proteome of ΦM12, the first report of the genome and proteome of a rhizobium-infecting T4-superfamily phage. The structural genes of ΦM12 are most similar to T4-superfamily phages of cyanobacteria. ΦM12 is the first reported T4-superfamily phage to lack genes encoding class I ribonucleotide reductase (RNR) and exonuclease dexA, and to possess a class II coenzyme B12-dependent RNR. ΦM12's novel collection of genes establishes it as the founder of a new group of T4-superfamily phages, fusing features of cyanophages and phages of enteric bacteria.

Prevalence of viral photosynthetic and capsid protein genes from cyanophages in two large and deep perialpine lakes

Cyanophages are important components of aquatic ecosystems, but their genetic diversity has been little investigated in freshwaters. A yearlong survey was conducted in surface waters of the two largest natural perialpine lakes in France (Lake Annecy and Lake Bourget) to investigate part of this cyanophage diversity through the analysis of both structural (e.g., g20) and functional (e.g., psbA) genes. We found that these cyanophage signature genes were prevalent throughout the year but that the community compositions of g20 cyanomyoviruses were significantly different between the two lakes. In contrast, psbA-containing cyanophages seemed to be more similar between the two ecosystems. We also found that a large proportion of g20 sequences grouped with cyanomyophage isolates. psbA sequences, belonging to phages of Synechococcus spp., were characterized by distinct triplet motifs (with a novel viral triplet motif, EFE). Thus, our results show that cyanophages (i) are a diverse viral community in alpine lakes and (ii) are clearly distinct from some other freshwater and marine environments, suggesting the influence of unique biogeographic factors.

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.

Seasonality and monthly dynamics of marine myovirus communities

Marine myoviruses (i.e. bacteriophages with a contractile tail sheath) are numerically abundant and genetically diverse. We developed a terminal restriction fragment length polymorphism assay (TRFLP) for g23, the conserved gene encoding the major capsid protein, to investigate T4-like myovirus communities at USC's Microbial Observatory at the San Pedro Ocean Time-series (SPOT), where we previously reported bacterial seasonality. Between 71 and 154 operational taxonomic units (OTUs) were observed monthly over 3 years. Roughly 25% of OTUs were detected in 31 or more months. T4-like myoviral community structure varied seasonally with some OTUs peaking repeatedly in spring-summer and others in fall-winter, while moderately abundant OTUs persisted year-round. Recurring community structure was demonstrated using discriminant function analysis (DFA, selecting taxa that best predict months) and average Bray-Curtis similarity. DFA showed communities from adjacent months or 12 months apart were positively auto-correlated, while communities 3-7 months apart were negatively auto-correlated. Bray-Curtis similarity was highest between adjacent months - with a local maximum at 12-month and local minima at 6- and 18- to 20-month lags. The T4-like virus community at SPOT exhibited seasonality, yet the somewhat unexpected persistence of moderately abundant OTUs and predictability of the community add new twists to existing conceptual models of marine viruses.