Strong latitudinal and vertical biogeography of Synechococcus diversity in the equatorial Pacific Ocean

Extensive genetic diversity has been observed within Synechococcus including the presence of multiple major clades. However, the biogeography and underlying environmental drivers of these clades remain elusive. Here, we developed a new high-throughput sequencing assay using rpoC1 as marker combined with Illumina sequencing. Using this, we identified the genetic diversity of Synechococcus from ~200 samples in an eastern Pacific Ocean transect between 19˚N and 3˚S. We used a placement method to identify the phylogenetic affiliation of each sequence and detected extensive diversity including multiple previously undescribed clades. We observed clear biogeographical domains, with Clade 2 dominant in the northern part of the transect, Clade CRD peaking at the equator, and Clade 1 dominant deeper in the water column throughout the transect. This biogeography, along with physical and nutrient data, suggests that Clade 2 represents a high temperature, low macronutrient ecotype, CRD a high temperature but low iron ecotype, and at least part of Clade 1 a low-light ecotype. The shift between Clade 2 and CRD occurred at 7˚N, whereas the concentration of macronutrients was low down to 4˚N, before increasing. This biogeography indicates that Synechococcus cells experience iron stress up to 7˚N despite low concentrations of phosphate and nitrate. The overall biogeography closely matched the distribution of Prochlorococcus diversity in this region, suggesting a parallel evolution of ecotypes in these two major lineages of marine Cyanobacteria.