Ecology of uncultured Prochlorococcus clades revealed through single-cell genomics and biogeographic analysis

Prochlorococcus is the numerically dominant photosynthetic organism throughout much of the world''s oceans, yet little is known about the ecology and genetic diversity of populations inhabiting tropical waters. To help close this gap, we examined natural Prochlorococcus communities in the tropical Pacific Ocean using a single-cell whole-genome amplification and sequencing. Analysis of the gene content of just 10 single cells from these waters added 394 new genes to the Prochlorococcus pan-genome—that is, genes never before seen in a Prochlorococcus cell. Analysis of marker genes, including the ribosomal internal transcribed sequence, from dozens of individual cells revealed several representatives from two uncultivated clades of Prochlorococcus previously identified as HNLC1 and HNLC2. While the HNLC clades can dominate Prochlorococcus communities under certain conditions, their overall geographic distribution was highly restricted compared with other clades of Prochlorococcus. In the Atlantic and Pacific oceans, these clades were only found in warm waters with low Fe and high inorganic P levels. Genomic analysis suggests that at least one of these clades thrives in low Fe environments by scavenging organic-bound Fe, a process previously unknown in Prochlorococcus. Furthermore, the capacity to utilize organic-bound Fe appears to have been acquired horizontally and may be exchanged among other clades of Prochlorococcus. Finally, one of the single Prochlorococcus cells sequenced contained a partial genome of what appears to be a prophage integrated into the genome.     

A novel clade of Prochlorococcus found in high nutrient low chlorophyll waters in the South and Equatorial Pacific Ocean

A novel high-light (HL)-adapted Prochlorococcus clade was discovered in high nutrient and low chlorophyll (HNLC) waters in the South Pacific Ocean by phylogenetic analyses of 16S ribosomal RNA (rRNA) and 16S–23S internal transcribed spacer (ITS) sequences. This clade, named HNLC fell within the HL-adapted Prochlorococcus clade with sequences above 99% similarity to one another, and was divided into two subclades, HNLC1 and HNLC2. The distribution of the whole HNLC clade in a northwest to southeast transect in the South Pacific (HNLC-to-gyre) and two 8°N to 8°S transects in the Equatorial Pacific was determined by quantitative PCR using specific primers targeting ITS regions. HNLC was the dominant HL Prochlorococcus clade (2–9% of bacterial 16S rRNA genes) at the three westernmost stations in the South Pacific but decreased to less than 0.1% at the other stations being replaced by the eMIT9312 ecotype in the hyperoligotrophic gyre. The highest contributions of HNLC Prochlorococcus in both Equatorial Pacific transects along the latitudinal lines of 170°W and 155°W were observed at the southernmost stations, reaching 16 and 6% of bacterial 16S rRNA genes, respectively, whereas eMIT9312 dominated near the Equator. Spearman Rank Order correlation analysis indicated that although both the HNLC clade and eMIT9312 were correlated with temperature, they showed different correlations with regard to nutrients. HNLC only showed significant correlations to ammonium uptake and regeneration rates, whereas eMIT9312 was negatively correlated with inorganic nutrients.     

Genomic mosaicism underlies the adaptation of marine Synechococcus ecotypes to distinct oceanic iron niches

Phytoplankton are limited by iron (Fe) in ~40% of the world's oceans including high-nutrient low-chlorophyll (HNLC) regions. While low-Fe adaptation has been well-studied in large eukaryotic diatoms, less is known for small, prokaryotic marine picocyanobacteria. This study reveals key physiological and genomic differences underlying Fe adaptation in marine picocyanobacteria. HNLC ecotype CRD1 strains have greater physiological tolerance to low Fe congruent with their expanded repertoire of Fe transporter, storage and regulatory genes compared to other ecotypes. From metagenomic analysis, genes encoding ferritin, flavodoxin, Fe transporters and siderophore uptake genes were more abundant in low-Fe waters, mirroring paradigms of low-Fe adaptation in diatoms. Distinct Fe-related gene repertories of HNLC ecotypes CRD1 and CRD2 also highlight how coexisting ecotypes have evolved independent approaches to life in low-Fe habitats. Synechococcus and Prochlorococcus HNLC ecotypes likewise exhibit independent, genome-wide reductions of predicted Fe-requiring genes. HNLC ecotype CRD1 interestingly was most similar to coastal ecotype I in Fe physiology and Fe-related gene content, suggesting populations from these different biomes experience similar Fe-selective conditions. This work supports an improved perspective that phytoplankton are shaped by more nuanced Fe niches in the oceans than previously implied from mostly binary comparisons of low- versus high-Fe habitats and populations.     

Ocean time-series reveals recurring seasonal patterns of virioplankton dynamics in the northwestern Sargasso Sea

There are an estimated 10³⁰ virioplankton in the world oceans, the majority of which are phages (viruses that infect bacteria). Marine phages encompass enormous genetic diversity, affect biogeochemical cycling of elements, and partially control aspects of prokaryotic production and diversity. Despite their importance, there is a paucity of data describing virioplankton distributions over time and depth in oceanic systems. A decade of high-resolution time-series data collected from the upper 300 m in the northwestern Sargasso Sea revealed recurring temporal and vertical patterns of virioplankton abundance in unprecedented detail. An annual virioplankton maximum developed between 60 and 100 m during periods of summer stratification and eroded during winter convective mixing. The timing and vertical positioning of this seasonal pattern was related to variability in water column stability and the dynamics of specific picophytoplankton and heterotrophic bacterioplankton lineages. Between 60 and 100 m, virioplankton abundance was negatively correlated to the dominant heterotrophic bacterioplankton lineage SAR11, as well as the less abundant picophytoplankton, Synechococcus. In contrast, virioplankton abundance was positively correlated to the dominant picophytoplankton lineage Prochlorococcus, and the less abundant alpha-proteobacteria, Rhodobacteraceae. Seasonally, virioplankton abundances were highly synchronous with Prochlorococcus distributions and the virioplankton to Prochlorococcus ratio remained remarkably constant during periods of water column stratification. The data suggest that a significant fraction of viruses in the mid-euphotic zone of the subtropical gyres may be cyanophages and patterns in their abundance are largely determined by Prochlorococcus dynamics in response to water column stability. This high-resolution, decadal survey of virioplankton abundance provides insight into the possible controls of virioplankton dynamics in the open ocean.     

Resolution of Prochlorococcus and Synechococcus Ecotypes by Using 16S-23S Ribosomal DNA Internal Transcribed Spacer Sequences

Cultured isolates of the marine cyanobacteria Prochlorococcus and Synechococcus vary widely in their pigment compositions and growth responses to light and nutrients, yet show greater than 96% identity in their 16S ribosomal DNA (rDNA) sequences. In order to better define the genetic variation that accompanies their physiological diversity, sequences for the 16S-23S rDNA internal transcribed spacer (ITS) region were determined in 32 Prochlorococcus isolates and 25 Synechococcus isolates from around the globe. Each strain examined yielded one ITS sequence that contained two tRNA genes. Dramatic variations in the length and G+C content of the spacer were observed among the strains, particularly among Prochlorococcus strains. Secondary-structure models of the ITS were predicted in order to facilitate alignment of the sequences for phylogenetic analyses. The previously observed division of Prochlorococcus into two ecotypes (called high and low-B/A after their differences in chlorophyll content) were supported, as was the subdivision of the high-B/A ecotype into four genetically distinct clades. ITS-based phylogenies partitioned marine cluster A Synechococcus into six clades, three of which can be associated with a particular phenotype (motility, chromatic adaptation, and lack of phycourobilin). The pattern of sequence divergence within and between clades is suggestive of a mode of evolution driven by adaptive sweeps and implies that each clade represents an ecologically distinct population. Furthermore, many of the clades consist of strains isolated from disparate regions of the world's oceans, implying that they are geographically widely distributed. These results provide further evidence that natural populations of Prochlorococcus and Synechococcus consist of multiple coexisting ecotypes, genetically closely related but physiologically distinct, which may vary in relative abundance with changing environmental conditions.     

Co-occurring Synechococcus ecotypes occupy four major oceanic regimes defined by temperature,macronutrients and iron

Marine picocyanobacteria, comprised of the genera Synechococcus and Prochlorococcus, are the most abundant and widespread primary producers in the ocean. More than 20 genetically distinct clades of marine Synechococcus have been identified, but their physiology and biogeography are not as thoroughly characterized as those of Prochlorococcus. Using clade-specific qPCR primers, we measured the abundance of 10 Synechococcus clades at 92 locations in surface waters of the Atlantic and Pacific Oceans. We found that Synechococcus partition the ocean into four distinct regimes distinguished by temperature, macronutrients and iron availability. Clades I and IV were prevalent in colder, mesotrophic waters; clades II, III and X dominated in the warm, oligotrophic open ocean; clades CRD1 and CRD2 were restricted to sites with low iron availability; and clades XV and XVI were only found in transitional waters at the edges of the other biomes. Overall, clade II was the most ubiquitous clade investigated and was the dominant clade in the largest biome, the oligotrophic open ocean. Co-occurring clades that occupy the same regime belong to distinct evolutionary lineages within Synechococcus, indicating that multiple ecotypes have evolved independently to occupy similar niches and represent examples of parallel evolution. We speculate that parallel evolution of ecotypes may be a common feature of diverse marine microbial communities that contributes to functional redundancy and the potential for resiliency.     

Spatial distribution of small phytoplankton composition in the Chukchi Sea

This study was conducted as part of the second Russian American long-term census of the Arctic in 2009 with sampling across the territorial of the Russian Federation in the Chukchi Sea since recent information regarding the standing stocks of pico- and nano-plankton is very limited. Using flow cytometry, small size fractions of the phytoplankton (<20 μm) obtained during the cruise were analyzed for the contributions of pico- and nano-eukaryotic phytoplankton and prokaryotic phytoplankton for the first time in the Chukchi Sea. The salinity and temperature in the southern part were significantly higher than those in the northern part of the Chukchi Sea. The abundance of Prochlorococcus and Synechococcus represented about 30 % of total small phytoplankton cells with their cell abundances reaching 295 and 590 cells ml^?1, respectively, although a significant portion (about 70 %) of the small phytoplankton community (<20 μm) was pico- and nano-eukaryotic plankton in the Chukchi Sea. Among different environmental factors, we found that temperature and NH₄ concentrations were significantly positively correlated with the abundance of total phytoplankton, pico- and nano-eukaryotes, Prochlorococcus, and Synechococcus. Projected higher water temperature and increase in NH₄ concentration condition in the Arctic Ocean as well as the Chukchi Sea could have fostered more small phytoplankton communities especially Prochlorococcus and Synechococcus.     

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

Viruses that infect the marine cyanobacterium Prochlorococcus have the potential to impact the growth, productivity, diversity and abundance of their hosts. In this study, changes in the microdiversity of cyanomyoviruses were investigated in 10 environmental samples taken along a North–South Atlantic Ocean transect using a myoviral-specific PCR-sequencing approach. Phylogenetic analyses of 630 viral g20 clones from this study, with 786 published g20 sequences, revealed that myoviral populations in the Atlantic Ocean had higher diversity than previously reported, with several novel putative g20 clades. Some of these clades were detected throughout the Atlantic Ocean. Multivariate statistical analyses did not reveal any significant correlations between myoviral diversity and environmental parameters, although myoviral diversity appeared to be lowest in samples collected from the north and south of the transect where Prochlorococcus diversity was also lowest. The results were correlated to the abundance and diversity of the co-occurring Prochlorococcus and Synechococcus populations, but revealed no significant correlations to either of the two potential host genera. This study provides evidence that cyanophages have extremely high and variable diversity and are distributed over large areas of the Atlantic Ocean.     

Ecological Genomics of Marine Picocyanobacteria

Summary: Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus numerically dominate the picophytoplankton of the world ocean, making a key contribution to global primary production. Prochlorococcus was isolated around 20 years ago and is probably the most abundant photosynthetic organism on Earth. The genus comprises specific ecotypes which are phylogenetically distinct and differ markedly in their photophysiology, allowing growth over a broad range of light and nutrient conditions within the 45°N to 40°S latitudinal belt that they occupy. Synechococcus and Prochlorococcus are closely related, together forming a discrete picophytoplankton clade, but are distinguishable by their possession of dissimilar light-harvesting apparatuses and differences in cell size and elemental composition. Synechococcus strains have a ubiquitous oceanic distribution compared to that of Prochlorococcus strains and are characterized by phylogenetically discrete lineages with a wide range of pigmentation. In this review, we put our current knowledge of marine picocyanobacterial genomics into an environmental context and present previously unpublished genomic information arising from extensive genomic comparisons in order to provide insights into the adaptations of these marine microbes to their environment and how they are reflected at the genomic level.     

Codon usage patterns and adaptive evolution of marine unicellular cyanobacteria Synechococcus and Prochlorococcus

Marine unicellular cyanobacteria, represented by Synechococcus and Prochlorococcus, dominate the total phytoplankton biomass and production in oligotrophic ocean. In this study, we employed comparative genomics approaches to extensively investigate synonymous codon usage bias and evolutionary rates in a large number of closely related species of marine unicellular cyanobacteria. Although these two groups of marine cyanobacteria have a close phylogenetic relationship, we find that they are highly divergent not only in codon usage patterns but also in the driving forces behind the diversification. It is revealed that in Prochlorococcus, mutation and genome compositional constraints are the main forces contributing to codon usage bias, whereas in Synechococcus, translational selection. In addition, nucleotide substitution rate analysis indicates that they are not evolving at a constant rate after the divergence and that the average d_N/d_S values of core genes in Synechococcus are significantly higher than those in Prochlorococcus. Our evolutionary genomic analysis provides the first insight into codon usage, evolutionary genetic mechanisms and environmental adaptation of Synechococcus and Prochlorococcus after divergence.     

East Tibetan Lakes Harbour Novel Clusters of Picocyanobacteria as Inferred from the 16S–23S rRNA Internal Transcribed Spacer Sequences

Planktonic picocyanobacteria abundance and diversity were investigated in nine lakes on the East Tibetan Plateau spanning a salinity gradient of 0.4–22.6 g l⁻¹. The investigation was conducted using epifluorescence microscopy (EFM) and terminal restriction fragment polymorphism analysis of 16S–23S rRNA internal transcribed spacer (ITS) PCR amplicons followed by sequence analyses of large ITS clone libraries of seven selected samples. EFM showed that picocyanobacteria comprised 7–19% of the total prokaryotic cells found in surface water. Most of the clones were classified into six clusters and grouped within the “picocyanobacterial clade”, which consists exclusively of freshwater Synechococcus. Four new phylogenetic clusters and one new subcluster of Synechococcus spp. were found, none of which are members of any known picocyanobacterial clusters. The new clusters and subcluster were the most abundant picocyanobacteria (about 96% of the sequences) in the samples collected. Sequence analyses indicated that members of the four new Synechococcus groups were only found in freshwater lakes (<1.0 g l⁻¹ of total dissolved solid), while members of the new subcluster were found in all the investigated Tibetan lakes, over a large salinity gradient of 0.4–22.6 g l ⁻¹. This suggests that there is ecologically significant microdiversity within the observed Synechococcus group as defined by ITS sequences. Collectively our study demonstrated abundant and potentially novel Synechococcus in East Tibetan lakes that are likely the result of evolutionary adaptations to regional conditions.     

Potential for phosphite and phosphonate utilization by Prochlorococcus

Phosphonates (Pn) are diverse organic phosphorus (P) compounds containing C–P bonds and comprise up to 25% of the high-molecular weight dissolved organic P pool in the open ocean. Pn bioavailability was suggested to influence markedly bacterial primary production in low-P areas. Using metagenomic data from the Global Ocean Sampling expedition, we show that the main potential microbial contributor in Pn utilization in oceanic surface water is the globally important marine primary producer Prochlorococcus. Moreover, a number of Prochlorococcus strains contain two distinct putative Pn uptake operons coding for ABC-type Pn transporters. On the basis of microcalorimetric measurements, we find that each of the two different putative Pn-binding protein (PhnD) homologs transcribed from these operons possesses different Pn- as well as inorganic phosphite-binding specificities. Our results suggest that Prochlorococcus adapt to low-P environments by increasing the number of Pn transporters with different specificities towards phosphite and different Pns.     

Cyanophycin Production in a Phycoerythrin-Containing Marine Synechococcus Strain of Unusual Phylogenetic Affinity

Thirty-two strains of phycoerythrin-containing marine picocyanobacteria were screened for the capacity to produce cyanophycin, a nitrogen storage compound synthesized by some, but not all, cyanobacteria. We found that one of these strains, Synechococcus sp. strain G2.1 from the Arabian Sea, was able to synthesize cyanophycin. The cyanophycin extracted from the cells was composed of roughly equimolar amounts of arginine and aspartate (29 and 35 mol%, respectively), as well as a small amount of glutamate (15 mol%). Phylogenetic analysis, based on partial 16S ribosomal DNA (rDNA) sequence data, showed that Synechococcus sp. strain G2.1 formed a well-supported clade with several strains of filamentous cyanobacteria. It was not closely related to several other well-studied marine picocyanobacteria, including Synechococcus strains PCC7002, WH7805, and WH8018 and Prochlorococcus sp. strain MIT9312. This is the first report of cyanophycin production in a phycoerythrin-containing strain of marine or halotolerant Synechococcus, and its discovery highlights the diversity of this ecologically important functional group.     

Presence of toxin-antitoxin systems in picocyanobacteria and their ecological implications

Picocyanobacteria (mainly Synechococcus and Prochlorococcus) contribute significantly to ocean’s primary production. Toxin-Antitoxin (TA) systems present in bacteria and archaea are known to regulate cell growth in response to environmental stresses. However, little is known about the presence of TA systems in picocyanobacteria. This study investigated complete genomes of Synechococcus and Prochlorococcus to understand the prevalence of TA systems in picocyanobacteria. Using the TAfinder software, Type II TA systems were predicted in 27 of 33 (81%) Synechococcus strains, but none of 38 Prochlorococcus strains contain TA genes. Synechococcus strains with larger genomes tend to contain more putative type II TA systems. The number of TA pairs varies from 0 to 42 in Synechococcus strains isolated from various environments. A linear correlation between the genome size and the number of putative TA systems in both coastal and freshwater Synechococcus was established. In general, open ocean Synechococcus contain no or few TA systems, while coastal and freshwater Synechococcus contain more TA systems. The type II TA systems inhibit microbial translation via ribonucleases and allow cells to enter the “dormant” stage in adverse environments. Inheritance of TA genes in freshwater and coastal Synechococcus could confer a recoverable persister mechanism important to survive in variable environments.Subject terms: Microbial ecology, Biooceanography     

Modeling Selective Pressures on Phytoplankton in the Global Ocean

Our view of marine microbes is transforming, as culture-independent methods facilitate rapid characterization of microbial diversity. It is difficult to assimilate this information into our understanding of marine microbe ecology and evolution, because their distributions, traits, and genomes are shaped by forces that are complex and dynamic. Here we incorporate diverse forces—physical, biogeochemical, ecological, and mutational—into a global ocean model to study selective pressures on a simple trait in a widely distributed lineage of picophytoplankton: the nitrogen use abilities of Synechococcus and Prochlorococcus cyanobacteria. Some Prochlorococcus ecotypes have lost the ability to use nitrate, whereas their close relatives, marine Synechococcus, typically retain it. We impose mutations for the loss of nitrogen use abilities in modeled picophytoplankton, and ask: in which parts of the ocean are mutants most disadvantaged by losing the ability to use nitrate, and in which parts are they least disadvantaged? Our model predicts that this selective disadvantage is smallest for picophytoplankton that live in tropical regions where Prochlorococcus are abundant in the real ocean. Conversely, the selective disadvantage of losing the ability to use nitrate is larger for modeled picophytoplankton that live at higher latitudes, where Synechococcus are abundant. In regions where we expect Prochlorococcus and Synechococcus populations to cycle seasonally in the real ocean, we find that model ecotypes with seasonal population dynamics similar to Prochlorococcus are less disadvantaged by losing the ability to use nitrate than model ecotypes with seasonal population dynamics similar to Synechococcus. The model predictions for the selective advantage associated with nitrate use are broadly consistent with the distribution of this ability among marine picocyanobacteria, and at finer scales, can provide insights into interactions between temporally varying ocean processes and selective pressures that may be difficult or impossible to study by other means. More generally, and perhaps more importantly, this study introduces an approach for testing hypotheses about the processes that underlie genetic variation among marine microbes, embedded in the dynamic physical, chemical, and biological forces that generate and shape this diversity.     

Clade-Specific 16S Ribosomal DNA Oligonucleotides Reveal the Predominance of a Single Marine Synechococcus Clade throughout a Stratified Water Column in the Red Sea

Phylogenetic relationships among members of the marine Synechococcus genus were determined following sequencing of the 16S ribosomal DNA (rDNA) from 31 novel cultured isolates from the Red Sea and several other oceanic environments. This revealed a large genetic diversity within the marine Synechococcus cluster consistent with earlier work but also identified three novel clades not previously recognized. Phylogenetic analyses showed one clade, containing halotolerant isolates lacking phycoerythrin (PE) and including strains capable, or not, of utilizing nitrate as the sole N source, which clustered within the MC-A (Synechococcus subcluster 5.1) lineage. Two copies of the 16S rRNA gene are present in marine Synechococcus genomes, and cloning and sequencing of these copies from Synechococcus sp. strain WH 7803 and genomic information from Synechococcus sp. strain WH 8102 reveal these to be identical. Based on the 16S rDNA sequence information, clade-specific oligonucleotides for the marine Synechococcus genus were designed and their specificity was optimized. Using dot blot hybridization technology, these probes were used to determine the in situ community structure of marine Synechococcus populations in the Red Sea at the time of a Synechococcus maximum during April 1999. A predominance of genotypes representative of a single clade was found, and these genotypes were common among strains isolated into culture. Conversely, strains lacking PE, which were also relatively easily isolated into culture, represented only a minor component of the Synechococcus population. Genotypes corresponding to well-studied laboratory strains also appeared to be poorly represented in this stratified water column in the Red Sea.     

COMPARATIVE MOLECULAR EVOLUTION OF NEWLY DISCOVERED PICOCYANOBACTERIAL STRAINS REVEALS A PHYLOGENETICALLY INFORMATIVE VARIABLE REGION OF β‐PHYCOERYTHRIN1

The genetic diversity and phylogenetic position of 10 strains of picocyanobacteria from the Arabian Sea were examined using partial sequences from three loci: 16S rDNA, RNA polymerase rpoC1, and two elements of the phycoerythrin (PE) locus, cpeA and cpeB which encode for the α and β subunit of PE. Nine of the strains showed nearly identical spectral phenotypes based on the in vivo excitation spectrum for PE fluorescence emission and appear to be strains synthesizing a phycourobilin (PUB)–lacking PE. These strains include one, Synechococcus sp. G2.1, already known to be closely related to filamentous cyanobacteria and not to the commonly studied 5.1 subcluster of marine Synechococcus. The 10th strain was a PE‐lacking strain that was of interest because it was isolated from open‐ocean conditions where picocyanobacteria with this phenotype are relatively uncommon. Phylogenetic analysis of the concatenated 16S rDNA and rpoC1 data sets showed that none of the previously described strains were members of the 5.1 subcluster of marine Synechococcus, nor were they closely related to strain G2.1. Instead, they form a well‐supported and previously undescribed clade of cyanobacteria that is sister to Cyanobium. Thus, these strains represent the first PE‐containing Cyanobium from oceanic waters, and the lineage they define includes a strain with a PE‐lacking phenotype from the same environment. Analysis of the PE sequence data showed the PE apoprotein has evolved independently in the G2.1 lineage and the Cyanobium‐like lineage represented by the study strains. It also revealed a hypervariable region of the β‐subunit not described previously; variation in this region shows a pattern among a wide range of PE‐containing organisms congruent with the phylogenetic relationships inferred from other genes. This suggests that the PUB‐lacking spectral phenotype is more likely to have evolved in distantly related phylogenetic lineages by either divergent or convergent evolution than by lateral gene transfer. Both the conserved PE gene sequences and the inferred amino acid sequences for the hypervariable region show considerable divergence among Prochlorococcus PEs, red algal PEs, PUB‐containing PEs from the marine Synechococcus 5.1 subcluster, PEs from the Cyanobium‐like strains, and PEs from other cyanobacteria (including strain G2.1). Thus, it appears that the hypervariable region of the PE gene can be used as a taxon‐specific marker.     

Comparative Genomic and Phylogenomic Analyses Reveal a Conserved Core Genome Shared by Estuarine and Oceanic Cyanopodoviruses

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.