Marine Cyanophages Research Articles

Cyanobacteria-cyanophage interactions between freshwater and marine ecosystems based on large-scale cyanophage genomic analysis

The disparities in harmful algal blooms dynamics are largely attributed to variations in cyanobacteria populations within aquatic ecosystems. However, cyanobacteria-cyanophage interactions and their role in shaping cyanobacterial populations has been previously underappreciated. To address this knowledge gap, we isolated and sequenced 42 cyanophages from diverse water sources in China, with the majority (n = 35) originating from freshwater sources. We designated these sequences as the “Novel Cyanophage Genome sequence Collection” (NCGC). NCGC displayed notable genetic variations, with 95 % (40/42) of the sequences representing previously unidentified taxonomic ranks. By integrating NCGC with public data of cyanophages and cyanobacteria, we found evidence for more frequent historical cyanobacteria-cyanophage interactions in freshwater ecosystems. This was evidenced by a higher prevalence of prophage integrase-related genes in freshwater cyanophages (37.97 %) than marine cyanophages (7.42 %). In addition, freshwater cyanophages could infect a broader range of cyanobacteria orders (n = 4) than marine ones (n = 0). Correspondingly, freshwater cyanobacteria harbored more defense systems per million base pairs in their genomes, indicating more frequent phage infections. Evolutionary and cyanophage epidemiological studies suggest that interactions between cyanobacteria and cyanophages in freshwater and marine ecosystems are interconnected, and that brackish water can act as a transitional zone for freshwater and marine cyanophages. In conclusion, our research significantly expands the genetic information database of cyanophage, offering a wider selection of cyanophages to control harmful cyanobacterial blooms. Additionally, we represent a pioneering large-scale and comprehensive analysis of cyanobacteria and cyanophage sequencing data, and it provides theoretical guidance for the application of cyanophages in different environments.

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

ABSTRACTCyanobacteria 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.

Genetic engineering of marine cyanophages reveals integration but not lysogeny in T7-like cyanophages

Marine cyanobacteria of the genera Synechococcus and Prochlorococcus are the most abundant photosynthetic organisms on earth, spanning vast regions of the oceans and contributing significantly to global primary production. Their viruses (cyanophages) greatly influence cyanobacterial ecology and evolution. Although many cyanophage genomes have been sequenced, insight into the functional role of cyanophage genes is limited by the lack of a cyanophage genetic engineering system. Here, we describe a simple, generalizable method for genetic engineering of cyanophages from multiple families, that we named REEP for REcombination, Enrichment and PCR screening. This method enables direct investigation of key cyanophage genes, and its simplicity makes it adaptable to other ecologically relevant host-virus systems. T7-like cyanophages often carry integrase genes and attachment sites, yet exhibit lytic infection dynamics. Here, using REEP, we investigated their ability to integrate and maintain a lysogenic life cycle. We found that these cyanophages integrate into the host genome and that the integrase and attachment site are required for integration. However, stable lysogens did not form. The frequency of integration was found to be low in both lab cultures and the oceans. These findings suggest that T7-like cyanophage integration is transient and is not part of a classical lysogenic cycle.

Novel PhoH-encoding vibriophages with lytic activity against environmental Vibriostrains.

Cholera is a devastating diarrheal disease that accounts for more than 10% of children's lives worldwide, but its treatment is hampered by a rise in antibiotic resistance. One promising alternative to antibiotic therapy is the use of bacteriophages to treat antibiotic-resistant cholera infections, and control Vibrio cholera in clinical cases and in the environment, respectively. Here, we report four novel, closely related environmental myoviruses, VP4, VP6, VP18, and VP24, which we isolated from two environmental toxigenic Vibrio cholerae strains from river Kuja and Usenge beach in Kenya. High-throughput sequencing followed by bioinformatics analysis indicated that the genomes of the four bacteriophages have closely related sequences, with sizes of 148,180bp, 148,181bp, 148,179bp, and 148,179bp, and a G + C content of 36.4%. The four genomes carry the phoH gene, which is overrepresented in marine cyanophages. The isolated phages displayed a lytic activity against 15 environmental, as well as one clinical, Vibrio cholerae strains. Thus, these novel lytic vibriophages represent potential biocontrol candidates for water decontamination against pathogenic Vibrio cholerae and ought to be considered for future studies of phage therapy.

An uncultured marine cyanophage encodes an active phycobilisome proteolysis adaptor protein NblA.

Phycobilisomes (PBS) are large water-soluble membrane-associated complexes in cyanobacteria and some chloroplasts that serve as light-harvesting antennae for the photosynthetic apparatus. When deplete of nitrogen or sulphur, cyanobacteria readily degrade their phycobilisomes allowing the cell to replenish these vanishing nutrients. The key regulator in the degradation process is NblA, a small protein (∼6 kDa), which recruits proteases to the PBS. It was discovered previously that not only do cyanobacteria possess nblA genes but also that they are encoded by genomes of some freshwater cyanophages. A recent study, using assemblies from oceanic metagenomes, revealed genomes of a novel uncultured marine cyanophage lineage, representatives of which contain genes coding for the PBS degradation protein. Here, we examined the functionality of nblA-like genes from these marine cyanophages by testing them in a freshwater model cyanobacterial nblA knockout. One of the viral NblA variants could complement the non-bleaching phenotype and restore PBS degradation. Our findings reveal a functional NblA from a novel marine cyanophage lineage. Furthermore, we shed new light on the distribution of nblA genes in cyanobacteria and cyanophages.

Rates of Molecular Evolution in a Marine Synechococcus Phage Lineage.

Cyanophages are characterized by vast genomic diversity and the formation of stable ecotypes over time. The evolution of phage diversity includes vertical processes, such as mutation, and horizontal processes, such as recombination and gene transfer. Here, we study the contribution of vertical and horizontal processes to short-term evolution of marine cyanophages. Analyzing time series data of Synechococcus-infecting Myoviridae ecotypes spanning up to 17 years, we found a high contribution of recombination relative to mutation (r/m) in all ecotypes. Additionally, we found a molecular clock of substitution and recombination in one ecotype, RIM8. The estimated RIM8 evolutionary rates are 2.2 genome-wide substitutions per year (1.275 × 10−5 substitutions/site/year) and 29 genome-wide nucleotide alterations due to recombination per year. We found 26 variable protein families, of which only two families have a predicted functional annotation, suggesting that they are auxiliary metabolic genes with bacterial homologs. A comparison of our rate estimates to other phage evolutionary rate estimates in the literature reveals a negative correlation of phage substitution rates with their genome size. A comparison to evolutionary rates in bacterial organisms further shows that phages have high rates of mutation and recombination compared to their bacterial hosts. We conclude that the increased recombination rate in phages likely contributes to their vast genomic diversity.

A novel uncultured marine cyanophage lineage with lysogenic potential linked to a putative marine Synechococcus 'relic' prophage.

Marine cyanobacteria are important contributors to primary production in the ocean and their viruses (cyanophages) affect the ocean microbial communities. Despite reports of lysogeny in marine cyanobacteria, a genome sequence of such temperate cyanophages remains unknown although genomic analysis indicate potential for lysogeny in certain marine cyanophages. Using assemblies from Red Sea and Tara Oceans metagenomes, we recovered genomes of a novel uncultured marine cyanophage lineage, which contain, in addition to common cyanophage genes, a phycobilisome degradation protein NblA, an integrase and a split DNA polymerase. The DNA polymerase forms a monophyletic clade with a DNA polymerase from a genomic island in Synechococcus WH8016. The island contains a relic prophage that does not resemble any previously reported cyanophage but shares several genes with the newly identified cyanophages reported here. Metagenomic recruitment indicates that the novel cyanophages are widespread, albeit at low abundance. Here, we describe a novel potentially lysogenic cyanophage family, their abundance and distribution in the marine environment.

Phagen-codierte metabolische Hilfsgene krempeln den Wirtsstoffwechsel um

Auxiliary metabolic genes (AMGs) are commonly found in the genomes of marine cyanophages and on metagenomic scaffolds originating from phages. Many of them encode proteins involved in photosynthesis and are thought to increase phage fitness during infection. Here, we summarize the biochemical properties of AMG-encoded proteins involved in light harvesting pigment biosynthesis and assembly.

Transcriptome Analysis of a Bloom-Forming Cyanobacterium Microcystis aeruginosa during Ma-LMM01 Phage Infection.

Microcystis aeruginosa forms massive blooms in eutrophic freshwaters, where it is constantly exposed to lytic cyanophages. Unlike other marine cyanobacteria, M. aeruginosa possess remarkably abundant and diverse potential antiviral defense genes. Interestingly, T4-like cyanophage Ma-LMM01, which is the sole cultured lytic cyanophage infecting M. aeruginosa, lacks the host-derived genes involved in maintaining host photosynthesis and directing host metabolism that are abundant in other marine cyanophages. Based on genomic comparisons with closely related cyanobacteria and their phages, Ma-LMM01 is predicted to employ a novel infection program that differs from that of other marine cyanophages. Here, we used RNA-seq technology and in silico analysis to examine transcriptional dynamics during Ma-LMM01 infection to reveal host transcriptional responses to phage infection, and to elucidate the infection program used by Ma-LMM01 to avoid the highly abundant host defense systems. Phage-derived reads increased only slightly at 1 h post-infection, but significantly increased from 16% of total cellular reads at 3 h post-infection to 33% of all reads by 6 h post-infection. Strikingly, almost none of the host genes (0.17%) showed a significant change in expression during infection. However, like other lytic dsDNA phages, including marine cyanophages, phage gene dynamics revealed three expression classes: early (host-takeover), middle (replication), and late (virion morphogenesis). The early genes were concentrated in a single ∼5.8-kb window spanning 10 open reading frames (gp054–gp063) on the phage genome. None of the early genes showed homology to the early genes of other T4-like phages, including known marine cyanophages. Bacterial RNA polymerase (σ70) recognition sequences were also found in the upstream region of middle and late genes, whereas phage-specific motifs were not found. Our findings suggest that unlike other known T4-like phages, Ma-LMM01 achieves three sequential gene expression patterns with no change in host promoter activity. This type of infection that does not cause significant change in host transcriptional levels may be advantageous in allowing Ma-LMM01 to escape host defense systems while maintaining host photosynthesis.

Genomic differentiation among wild cyanophages despite widespread horizontal gene transfer.

BackgroundGenetic recombination is a driving force in genome evolution. Among viruses it has a dual role. For genomes with higher fitness, it maintains genome integrity in the face of high mutation rates. Conversely, for genomes with lower fitness, it provides immediate access to sequence space that cannot be reached by mutation alone. Understanding how recombination impacts the cohesion and dissolution of individual whole genomes within viral sequence space is poorly understood across double-stranded DNA bacteriophages (a.k.a phages) due to the challenges of obtaining appropriately scaled genomic datasets.ResultsHere we explore the role of recombination in both maintaining and differentiating whole genomes of 142 wild double-stranded DNA marine cyanophages. Phylogenomic analysis across the 51 core genes revealed ten lineages, six of which were well represented. These phylogenomic lineages represent discrete genotypic populations based on comparisons of intra- and inter- lineage shared gene content, genome-wide average nucleotide identity, as well as detected gaps in the distribution of pairwise differences between genomes. McDonald-Kreitman selection tests identified putative niche-differentiating genes under positive selection that differed across the six well-represented genotypic populations and that may have driven initial divergence. Concurrent with patterns of recombination of discrete populations, recombination analyses of both genic and intergenic regions largely revealed decreased genetic exchange across individual genomes between relative to within populations.ConclusionsThese findings suggest that discrete double-stranded DNA marine cyanophage populations occur in nature and are maintained by patterns of recombination akin to those observed in bacteria, archaea and in sexual eukaryotes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3286-x) contains supplementary material, which is available to authorized users.

Genomic diversification of marine cyanophages into stable ecotypes.

Understanding the structure and origin of natural bacteriophage genomic diversity is important in elucidating how bacteriophages influence the mortality rates and composition of their host communities. Here, we examine the genetic structure and genomic diversification of naturally occurring bacteriophages by analyzing the full genomic sequences of over 100 isolates of Synechococcus-infecting cyanophages collected over 15 years from coastal waters of Southern New England, USA. Our analysis revealed well-supported cyanophage genomic clusters (genome-wide average nucleotide identity (ANI) >93%) and subclusters (genome-wide ANI >98%) that remained consistent for a decade or longer. Furthermore, by combining the genomic data with genetic analysis of an additional 800 isolates and environmental amplicon sequence data both genomic clusters and subclusters were found to exhibit clear temporal and/or spatial patterns of abundance, suggesting that these units represent distinct viral ecotypes. The processes responsible for diversification of cyanophages into genomic clusters and subclusters were similar across genetic scales and included allelic exchange as well as gene gain and loss. Isolates belonging to different subclusters were found to differ in genes that encoded auxiliary metabolic functions, restriction modification enzymes, and virion structural proteins, although the specific traits and selection pressures responsible for the maintenance of distinct ecotypes remain unknown.

A Novel Strategy for Exploitation of Host RNase E Activity by a Marine Cyanophage.

Previous studies have shown that infection of Prochlorococcus MED4 by the cyanophage P-SSP7 leads to increased transcript levels of host endoribonuclease (RNase) E. However, it has remained enigmatic whether this is part of a host defense mechanism to degrade phage messenger RNA (mRNA) or whether this single-strand RNA-specific RNase is utilized by the phage. Here we describe a hitherto unknown means through which this cyanophage increases expression of RNase E during phage infection and concomitantly protects its own RNA from degradation. We identified two functionally different RNase E mRNA variants, one of which is significantly induced during phage infection. This transcript lacks the 5' UTR, is considerably more stable than the other transcript, and is likely responsible for increased RNase E protein levels during infection. Furthermore, selective enrichment and in vivo analysis of double-stranded RNA (dsRNA) during infection revealed that phage antisense RNAs (asRNAs) sequester complementary mRNAs to form dsRNAs, such that the phage protein-coding transcriptome is nearly completely covered by asRNAs. In contrast, the host protein-coding transcriptome is only partially covered by asRNAs. These data suggest that P-SSP7 orchestrates degradation of host RNA by increasing RNase E expression while masking its own transcriptome from RNase E degradation in dsRNA complexes. We propose that this combination of strategies contributes significantly to phage progeny production.

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.

Viruses Inhibit CO2 Fixation in the Most Abundant Phototrophs on Earth

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most numerous photosynthetic organisms on our planet [1, 2]. With a global population size of 3.6× 10(27) [3], they are responsible for approximately 10% of global primary production [3, 4]. Viruses that infect Prochlorococcus and Synechococcus (cyanophages) can be readily isolated from ocean waters [5-7] and frequently outnumber their cyanobacterial hosts [8]. Ultimately, cyanophage-induced lysis of infected cells results in the release of fixed carbon into the dissolved organic matter pool [9]. What is less well known is the functioning of photosynthesis during the relatively long latent periods of many cyanophages [10, 11]. Remarkably, the genomes of many cyanophage isolates contain genes involved in photosynthetic electron transport (PET) [12-18] as well as central carbon metabolism [14, 15, 19, 20], suggesting that cyanophages may play an active role in photosynthesis. However, cyanophage-encoded gene products are hypothesized to maintain or even supplement PET for energy generation while sacrificing wasteful CO2 fixation during infection [17, 18, 20]. Yet this paradigm has not been rigorously tested. Here, we measured the ability of viral-infected Synechococcus cells to fix CO2 as well as maintain PET. We compared two cyanophage isolates that share different complements of PET and central carbon metabolism genes. We demonstrate cyanophage-dependent inhibition of CO2 fixation early in the infection cycle. In contrast, PET is maintained throughout infection. Our data suggest a generalized strategy among marine cyanophages to redirect photosynthesis to support phage development, which has important implications for estimates of global primary production.

Synthesis of 2'-Fluoro RNA by Syn5 RNA polymerase.

The substitution of 2′-fluoro for 2′-hydroxyl moieties in RNA substantially improves the stability of RNA. RNA stability is a major issue in RNA research and applications involving RNA. We report that the RNA polymerase from the marine cyanophage Syn5 has an intrinsic low discrimination against the incorporation of 2′-fluoro dNMPs during transcription elongation. The presence of both magnesium and manganese ions at high concentrations further reduce this discrimination without decreasing the efficiency of incorporation. We have constructed a Syn5 RNA polymerase in which tyrosine 564 is replaced with phenylalanine (Y564F) that further decreases the discrimination against 2′-fluoro-dNTPs during RNA synthesis. Sequence elements in DNA templates that affect the yield of RNA and incorporation of 2′-fluoro-dNMPs by Syn5 RNA polymerase have been identified.

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.

An RNA polymerase from the ocean improves in vitro RNA synthesis (560.1)

Phages in the ocean harbor the greatest diversity of protein on the earth. Since a large portion of these proteins must play roles in nucleic acids metabolism, one would expect numerous novel mechanisms underlying the fundamental processes including transcription, DNA replication, and recombination and thus abundant molecular reagents therein. We recently characterized the single subunit RNA polymerase from marine cyanophage Syn5 and identified its promoter sequence. This marine enzyme catalyzes RNA synthesis over a wide range of temperature and salinity. Its processivity is greater than 30,000 nt without any significant intermediate products. Most important, Syn5 RNA polymerase produces precise run‐off transcripts with homogeneous 3’ termini. The enzyme predominantly used for in vitro run‐off RNA synthesis‐T7 RNA polymerase synthesizes, in addition to run‐off products of precise length, transcripts with an additional non‐base‐paired nucleotide at the 3’‐terminus (N+1 product). Therefore, the enzyme is advantageous for the production of RNAs that require precise 3’‐termini, such as tRNAs and RNA fragments that are used for subsequent assembly.

Two Novel Proteins of Cyanophage Syn5 Compose Its Unusual Horn Structure

The marine cyanophage Syn5 can be propagated to a high titer in the laboratory on marine photosynthetic Synechococcus sp. strain WH8109. The purified particles carry a novel slender horn structure projecting from the vertex opposite the tail vertex. The genome of Syn5 includes a number of genes coding for novel proteins. Using immune-electron microscopy with gold-labeled antibodies, we show that two of these novel proteins, products of genes 53 and 54, are part of the horn structure. A third novel protein, the product of gene 58, is assembled onto the icosahedral capsid lattice. Characterization of radioactively labeled precursor procapsids by sucrose gradient centrifugation shows that there appear to be three classes of particles-procapsids, scaffold-deficient procapsids, and expanded capsids. These lack fully assembled horn appendages. The horn presumably assembles onto the virion just before or after DNA packaging. Antibodies raised to the recombinant novel Syn5 proteins did not interfere with phage infectivity, suggesting that the functions of these proteins are not directly involved in phage attachment or infection of the host WH8109. The horn structure may represent some adaption to the marine environment, whose function will require additional investigation.

Syn5 RNA polymerase synthesizes precise run-off RNA products

The enzyme predominantly used for in vitro run-off RNA synthesis is bacteriophage T7 RNA polymerase. T7 RNA polymerase synthesizes, in addition to run-off products of precise length, transcripts with an additional non-base-paired nucleotide at the 3′-terminus (N + 1 product). This contaminating product is extremely difficult to remove. We recently characterized the single-subunit RNA polymerase from marine cyanophage Syn5 and identified its promoter sequence. This marine enzyme catalyses RNA synthesis over a wider range of temperature and salinity than does T7 RNA polymerase. Its processivity is >30 000 nt without significant intermediate products. The requirement for the initiating nucleotide at the promoter is less stringent for Syn5 RNA polymerase as compared to T7 RNA polymerase. A major difference is the precise run-off transcripts with homogeneous 3′-termini synthesized by Syn5 RNA polymerase. Therefore, the enzyme is advantageous for the production of RNAs that require precise 3′-termini, such as tRNAs and RNA fragments that are used for subsequent assembly.

Prevalence of psbA-containing cyanobacterial podoviruses in the ocean

Podoviruses that infect marine picocyanobacteria are abundant and could play a significant role on regulating host populations due to their specific phage-host relationship. Genome sequencing of cyanophages has unveiled that many marine cyanophages encode certain photosynthetic genes like psbA. It appears that psbA is only present in certain groups of cyanopodovirus isolates. In order to better understand the prevalence of psbA in cyanobacterial podoviruses, we searched the marine metagenomic database (GOS, BATS, HOT and MarineVirome). Our study suggests that 89% of recruited cyanopodovirus scaffolds from the GOS database contained the psbA gene, supporting the ecological relevance of the photosynthesis gene for surface oceanic cyanophages. Diversification between Clade A and B are consistent with recent finding of two major groups of cyanopodoviruses. All the data also shows that Clade B cyanopodoviruses dominate the surface ocean water, while Clade A cyanopodoviruses become more important in the coastal and estuarine environments.