Podoviruses are among the major viral groups that infect marine picocyanobacteria Prochlorococcus and Synechococcus. Here, we reported the genome sequences of five Synechococcus podoviruses isolated from the estuarine environment, and performed comparative genomic and phylogenomic analyses based on a total of 20 cyanopodovirus genomes. The genomes of all the known marine cyanopodoviruses are highly syntenic. A pan-genome of 349 clustered orthologous groups was determined, among which 15 were core genes. These core genes make up nearly half of each genome in length, reflecting the high level of genome conservation among this cyanophage type. The whole genome phylogenies based on concatenated core genes and gene content were highly consistent and confirmed the separation of two discrete marine cyanopodovirus clusters MPP-A and MPP-B. The genomes within cluster MPP-B grouped into subclusters mainly corresponding to Prochlorococcus or Synechococcus host types. Auxiliary metabolic genes tend to occur in a specific phylogenetic group of these cyanopodoviruses. All the MPP-B phages analyzed here encode the photosynthesis gene psbA, which are absent in all the MPP-A genomes thus far. Interestingly, all the MPP-B and two MPP-A Synechococcus podoviruses encode the thymidylate synthase gene thyX, while at the same genome locus all the MPP-B Prochlorococcus podoviruses encode the transaldolase gene talC. Both genes are hypothesized to have the potential to facilitate the biosynthesis of deoxynucleotide for phage replication. Inheritance of specific functional genes could be important to the evolution and ecological fitness of certain cyanophage genotypes. Our analyses demonstrate that cyanopodoviruses of estuarine and oceanic origins share a conserved core genome and suggest that accessory genes may be related to environmental adaptation.

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.

The phylogenetic diversity of picocyanobacteria in seven alkaline lakes on the Tibetan Plateau was analyzed using the molecular marker 16S-23S rRNA internal transcribed spacer sequence. A total of 1,077 environmental sequences retrieved from the seven lakes were grouped into seven picocyanobacterial clusters, with two clusters newly described here. Each of the lakes was dominated by only one or two clusters, while different lakes could have disparate communities, suggesting low alpha diversity but high beta diversity of picocyanobacteria in these high-altitude freshwater and saline lakes. Several globally distributed clusters were found in these Tibetan lakes, such as subalpine cluster I and the Cyanobium gracile cluster. Although other clusters likely exhibit geographic restriction to the plateau temporally, reflecting endemicity, they can indeed be distributed widely on the plateau. Lakes with similar salinities may have similar genetic populations despite a large geographic distance. Canonical correspondence analysis identified salinity as the only environmental factor that may in part explain the diversity variations among lakes. Mantel tests suggested that the community similarities among lakes are independent of geographic distance. A portion of the picocyanobacterial clusters appear to be restricted to a narrow salinity range, while others are likely adapted to a broad range. A seasonal survey of Lake Namucuo across 3 years did not show season-related variations in diversity, and depth-related population partitioning was observed along a vertical profile of the lake. Our study emphasizes the high dispersive potential of picocyanobacteria and suggests that the regional distribution may result from adaptation to specified environments.