Diversity and genome dynamics of marine cyanophages using metagenomic analyses

Novel virulent and temperate cyanophages predicted to infect Microcoleus associated with anatoxin-producing benthic mats

Cyanophages are crucial for regulating cyanobacterial populations, but their influence on anatoxin-producing Microcoleus mat dynamics remains unexplored. Here, we use metagenomics to explore phage presence in benthic mats from the Wolastoq|Saint John River (New Brunswick, Canada) and the Eel River (California, USA). We recovered multiple viral-like sequences associated with different putative bacterial hosts, including two cyanophage genomes with apparently different replication strategies. A temperate cyanophage was found integrated in the genomes of Microcoleus sp. 3 recovered from the Eel River and is phylogenetically related to Phormidium phages. We also recovered novel virulent cyanophage genomes from Wolastoq and Eel River mats that were dominated by anatoxin-producing Microcoleus species predicted to be the host. Despite the geographical distance, these genomes have similar sizes (circa 239 kbp) and share numerous orthologous genes with high sequence identity. A considerable reduction of the anatoxin-producing Microcoleus species in Wolastoq mats following the emergence of the virulent phage suggests that phage infections have an important role in limiting the abundance of this toxigenic cyanobacterium and releasing anatoxins into the surrounding water. Our results constitute the first report of cyanophages predicted to infect mat-forming Microcoleus species associated with anatoxin production.

Characterization and genomic analysis of a novel Synechococcus phage S-H9-2 belonging to Bristolvirus genus isolated from the Yellow Sea

Cyanophages are known to influence the population dynamics and community structure of cyanobacteria and thus play an important role in biogeochemical cycles in aquatic ecosystems. In this study, a novel Synechococcus phage S-H9-2 infecting Synechococcus sp. WH 8102 was isolated from the coastal water of the Yellow Sea. Synechococcus phage S-H9-2 contains a 187,320 bp genome of double-stranded DNA with a G + C content of 40.3%, 202 potential open reading frames (ORFs), and 15 tRNAs. Phylogenetic analysis and nucleotide-based intergenomic similarity suggest that Synechococcus phage S-H9-2 belongs to the Bristolvirus genus under the family Kyanoviridae. Homologs of the S-H9-2 open reading frame can be found in a variety of marine environments, as shown by the results of mapping the genome sequence of S-H9-2 to the Global Ocean Viromes 2.0 dataset. The presence of auxiliary metabolic genes (AMGs) related to photosynthesis, carbon metabolism, and phosphorus assimilation, as well as phylogenetic relationships based on complete genome sequences, reflect the mechanism of phage-host interaction and host-specific strategies for adaptation to environmental conditions. This study enriches the current genomic database of cyanophage and contributed to our understanding of the virus-host interactions and their adaption to the environment.

A Unique Set of Auxiliary Metabolic Genes Found in an Isolated Cyanophage Sheds New Light on Marine Phage-Host Interactions

Cyanophages, viruses that infect cyanobacteria, are abundant and widely distributed in aquatic ecosystems, playing important roles in regulating the abundance, activity, diversity, and evolution of cyanobacteria. A T4-like cyanophage, S-SCSM1, infecting Synechococcus and Prochlorococcus strains of different ecotypes, was isolated from the South China Sea in this study. For the first time, a mannose-6-phosphate isomerase (MPI) gene was identified in the cultured cyanophage. At least 11 phylogenetic clusters of cyanophage MPIs were retrieved and identified from the marine metagenomic data sets, indicating that cyanophage MPIs in the marine environment are extremely diverse. The existence of 24 genes encoding 2-oxoglutarate (2OG)-Fe(II) oxygenase superfamily proteins in the S-SCSM1 genome emphasizes their potential importance and diverse functions in reprogramming host metabolism during phage infection. Novel cell wall synthesis and modification genes found in the S-SCSM1 genome indicate that diverse phenotypic modifications imposed by phages on cyanobacterial hosts remain to be discovered. Two noncoding RNAs of cis-regulatory elements in the S-SCSM1 genome were predicted to be associated with host exopolysaccharide metabolism and photosynthesis. The isolation and genomic characterization of cyanophage S-SCSM1 provide more information on the genetic diversity of cyanophages and phage-host interactions in the marine environment. IMPORTANCE Cyanophages play important ecological roles in aquatic ecosystems. Genomic and proteomic characterizations of the T4-like cyanophage S-SCSM1 indicate that novel and diverse viral genes and phage-host interactions in the marine environment remain unexplored. The first identified mannose-6-phosphate isomerase (MPI) gene from a cultured cyanophage was found in the S-SCSM1 genome, although MPIs were previously found in viral metagenomes at high frequencies similar to those of the cyanophage photosynthetic gene psbA. The presence of 24 genes encoding 2-oxoglutarate (2OG)-Fe(II) oxygenase superfamily proteins, novel cell wall synthesis and modification genes, a nonbleaching protein A gene, and 2 noncoding RNAs of cis-regulatory elements in the S-SCSM1 genome as well as the presence of a virion-associated regulatory protein indicate the diverse functions that cyanophages have in reprogramming the metabolism and modifying the phenotypes of hosts during infection.

Isolation and Characterization of a Novel Cyanophage Encoding Multiple Auxiliary Metabolic Genes

As significant drivers of cyanobacteria mortality, cyanophages have been known to regulate the population dynamics, metabolic activities, and community structure of this most important marine autotrophic picoplankton and, therefore, influence the global primary production and biogeochemical cycle in aquatic ecosystems. In the present study, a lytic Synechococcus phage, namely S-SZBM1, was isolated and identified. Cyanophage S-SZBM1 has a double-stranded DNA genome of 177,834 bp with a G+C content of 43.31% and contains a total of 218 predicted ORFs and six tRNA genes. Phylogenetic analysis and nucleotide-based intergenomic similarity suggested that cyanophage S-SZBM1 belongs to a new genus under the family Kyanoviridae. A variety of auxiliary metabolic genes (AMGs) that have been proved or speculated to relate to photosynthesis, carbon metabolism, nucleotide synthesis and metabolism, cell protection, and other cell metabolism were identified in cyanophage S-SZBM1 genome and may affect host processes during infection. In addition, 24 of 32 predicted structural proteins were identified by a high-throughput proteome analysis which were potentially involved in the assembly processes of virion. The genomic and proteomic analysis features of cyanophage S-SZBM1 offer a valuable insight into the interactions between cyanophages and their hosts during infection.

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.

Host Cyanobacteria Killing by Novel Lytic Cyanophage YongM: A Protein Profiling Analysis

Cyanobacteria are autotrophic prokaryotes that can proliferate robustly in eutrophic waters through photosynthesis. This can lead to outbreaks of lake “water blooms”, which result in water quality reduction and environmental pollution that seriously affect fisheries and aquaculture. The use of cyanophages to control the growth of cyanobacteria is an important strategy to tackle annual cyanobacterial blooms. YongM is a novel lytic cyanophage with a broad host spectrum and high efficiency in killing its host, cyanobacteria FACHB-596. However, changes in cyanophage protein profile during infestation and killing of the host remains unknown. To characterize the proteins and its regulation networks involved in the killing of host cyanobacteria by YongM and evaluate whether this strain YongM could be used as a chassis for further engineering to be a powerful tool in dealing with cyanobacterial blooms, we herein applied 4D label-free high-throughput quantitative proteomics to analyze differentially expressed proteins (DEPs) involved in cyanobacteria host response infected 1 and 8 h with YongM cyanophage. Metabolic pathways, such as photosynthesis, photosynthesis-antennal protein, oxidative phosphorylation, ribosome, carbon fixation, and glycolysis/glycol-isomerization were significantly altered in the infested host, whereas DEPs were associated with the metabolic processes of photosynthesis, precursor metabolites, energy production, and organic nitrogen compounds. Among these DEPs, key proteins involved in YongM-host interaction may be photosystem I P700 chlorophyll-a apolipoprotein, carbon dioxide concentration mechanism protein, cytochrome B, and some YongM infection lysis-related enzymes. Our results provide comprehensive information of protein profiles during the invasion and killing of host cyanobacteria by its cyanophage, which may shed light on future design and manipulation of artificial cyanophages against water blooms.

Picocyanobacteria Community and Cyanophage Infection Responses to Nutrient Enrichment in a Mesocosms Experiment in Oligotrophic Waters

Prochlorococcus and Synechococcus are pico-sized cyanobacteria that play a fundamental role in oceanic primary production, being particularly important in warm, nutrient-poor waters. Their potential response to nutrient enrichment is expected to be contrasting and to differ from larger phytoplankton species. Here, we used a metagenomic approach to characterize the responses to nutrient enrichment in the community of picocyanobacteria and to analyze the cyanophage response during a mesocosms experiment in the oligotrophic Red Sea. Natural picoplankton community was dominated by Synechococcus clade II, with marginal presence of Prochlorococcus (0.3% bacterial reads). Increased nutrient input triggered a fast Synechococcus bloom, with clade II being the dominant, with no response of Prochlorococcus growth. The largest bloom developed in the mesocosms receiving a single initial input of nutrients, instead of daily additions. The relative abundances of cyanophage sequences in cellular metagenomes increased during the experiment from 12.6% of total virus reads up to 40% in the treatment with the largest Synechococcus bloom. The subsequent collapse of the bloom pointed to a cyanophage infection on Synechococcus that reduced its competitive capacity, and was then followed by a diatom bloom. The cyanophage attack appears to have preferentially affected the most abundant Synechococcus clade II, increasing the evenness within the host population. Our results highlight the relevance of host-phage interactions on determining population dynamics and diversity of Synechococcus populations.

Biogeography and taxonomic overview of terrestrial hot spring thermophilic phages

Bacterial viruses ("phages") play important roles in the regulation and evolution of microbial communities in most ecosystems. Terrestrial hot springs typically contain thermophilic bacterial communities, but the diversity and impacts of its associated viruses ("thermophilic phages") are largely unexplored. Here, we provide a taxonomic overview of phages that have been isolated strictly from terrestrial hot springs around the world. In addition, we placed 17 thermophilic phage genomes in a global phylogenomic context to detect evolutionary patterns. Thermophilic phages have diverse morphologies (e.g., tailed, filamentous), unique virion structures (e.g., extremely long tailed siphoviruses), and span five taxonomic families encompassing strictly thermophilic phage genera. Within the phage proteomic tree, six thermophilic phage-related clades were identified, with evident genomic relatedness between thermophilic phages and archaeal viruses. Moreover, whole proteome analyses showed clustering between phages that infect distinct host phyla, such as Firmicutes and Deinococcus-Thermus. The potential for discovery of novel phage-host systems in terrestrial hot springs remain mostly untapped, thus additional emphasis on thermophilic phages in ecological prospecting is encouraged to gain insights into the microbial population dynamics of these environments.

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.

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.

Novel chaperonins are prevalent in the virioplankton and demonstrate links to viral biology and ecology

Chaperonins are protein-folding machinery found in all cellular life. Chaperonin genes have been documented within a few viruses, yet, surprisingly, analysis of metagenome sequence data indicated that chaperonin-carrying viruses are common and geographically widespread in marine ecosystems. Also unexpected was the discovery of viral chaperonin sequences related to thermosome proteins of archaea, indicating the presence of virioplankton populations infecting marine archaeal hosts. Virioplankton large subunit chaperonin sequences (GroELs) were divergent from bacterial sequences, indicating that viruses have carried this gene over long evolutionary time. Analysis of viral metagenome contigs indicated that: the order of large and small subunit genes was linked to the phylogeny of GroEL; both lytic and temperate phages may carry group I chaperonin genes; and viruses carrying a GroEL gene likely have large double-stranded DNA (dsDNA) genomes (>70 kb). Given these connections, it is likely that chaperonins are critical to the biology and ecology of virioplankton populations that carry these genes. Moreover, these discoveries raise the intriguing possibility that viral chaperonins may more broadly alter the structure and function of viral and cellular proteins in infected host cells.

Metabolic Genes within Cyanophage Genomes: Implications for Diversity and Evolution

Cyanophages, a group of viruses specifically infecting cyanobacteria, are genetically diverse and extensively abundant in water environments. As a result of selective pressure, cyanophages often acquire a range of metabolic genes from host genomes. The host-derived genes make a significant contribution to the ecological success of cyanophages. In this review, we summarize the host-derived metabolic genes, as well as their origin and roles in cyanophage evolution and important host metabolic pathways, such as the light-dependent reactions of photosynthesis, the pentose phosphate pathway, nutrient acquisition and nucleotide biosynthesis. We also discuss the suitability of the host-derived metabolic genes as potential diagnostic markers for the detection of genetic diversity of cyanophages in natural environments.

Complete genome sequence of bacteriophage P8625, the first lytic phage that infects Verrucomicrobia

Bacteriophage P8625 is a lytic bacteriophage that infects the verrucomicrobial strain IMCC8625, a marine bacterium affiliated with Verrucomicrobia subdivision 4. Both the bacteriophage and the host bacterial strain were isolated from surface seawater samples collected off the east coast of Korea. The phage particle has an icosahedral capsid with a diameter of ~47 nm and a long tail of ~75 nm in length, showing the distinctive morphology of the Siphoviridae family. The complete genome sequence of phage P8625 is 32,894 bp long with 51.0 % G + C content. This is the first report of the complete genome sequence of a lytic phage that infects the Verrucomicrobia, for which the name "verrucophage" is proposed.

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.