Temporal variation of Synechococcus clades at a coastal Pacific Ocean monitoring site

Temporal and Spatial Dynamics of Synechococcus Clade II and Other Microbes in the Eutrophic Subtropical San Diego Bay

The diversity of the marine cyanobacterium Synechococcus can be broadly separated into clades, with clade II typically present in warm oligotrophic water, and clades I and IV found in cooler coastal water. We found amplicon sequence variants (ASVs) belonging to clade II in the nutrient-replete waters of San Diego Bay (SDB). Using the 16S rRNA gene, 18S rRNA gene and internal transcribed spacer region sequencing, we analysed multiple locations in SDB monthly for over a year, with additional samples dating back to 2015. Synechococcus community composition differed from the nearby coast into SDB in terms of dominant clade and ASVs. Specific clade II ASVs became relatively more abundant towards the back of the bay and showed seasonality, with higher relative abundance in the warm months. Select ASVs group phylogenetically and show similar seasonal and spatial distribution patterns, indicating these ASVs have adapted to SDB. Isolates matching clade II ASVs from SDB show pigment composition that is better adapted to the green light available in SDB, further supporting our findings. Other microbial taxa also show SDB enrichment, providing evidence that SDB is a chemostat-like environment where circulation, temperature, light and other environmental conditions create a zone for microbial evolution and diversification.

Dissolved organic phosphorus bond-class utilization by Synechococcus

Dissolved organic phosphorus (DOP) contains compounds with phosphoester, phosphoanhydride, and phosphorus-carbon bonds. While DOP holds significant nutritional value for marine microorganisms, the bioavailability of each bond-class to the widespread cyanobacterium Synechococcus remains largely unknown. This study evaluates bond-class specific DOP utilization by Synechococcus strains from open and coastal oceans. Both strains exhibited comparable growth rates when provided phosphate, a phosphoanhydride [3-polyphosphate and 45-polyphosphate], or a DOP compound with both phosphoanhydride and phosphoester bonds (adenosine 5'-triphosphate). Growth rates on phosphoesters [glucose-6-phosphate, adenosine 5'-monophosphate, bis(4-methylumbelliferyl) phosphate] were variable, and neither strain grew on selected phosphorus-carbon compounds. Both strains hydrolyzed 3-polyphosphate, then adenosine 5'-triphosphate, and lastly adenosine 5'-monophosphate, exhibiting preferential enzymatic hydrolysis of phosphoanhydride bonds. The strains' exoproteomes contained phosphorus hydrolases, which combined with enhanced cell-free hydrolysis of 3-polyphosphate and adenosine 5'-triphosphate under phosphate deficiency, suggests active mineralization of phosphoanhydride bonds by these exoproteins. Synechococcus alkaline phosphatases presented broad substrate specificities, including activity toward the phosphoanhydride 3-polyphosphate, with varying affinities between strains. Collectively, these findings underscore the potentially significant role of compounds with phosphoanhydride bonds in Synechococcus phosphorus nutrition and highlight varied growth and enzymatic responses to molecular diversity within DOP bond-classes, thereby expanding our understanding of microbially mediated DOP cycling in marine ecosystems.

Substrate specificity and ecological significance of PstS homologs in phosphorus uptake in marine Synechococcus sp. WH8102

Phosphorus, a vital macronutrient, often limits primary productivity in marine environments. Marine Synechococcus strains, including WH8102, rely on high-affinity phosphate-binding proteins (PstS) to scavenge inorganic phosphate in oligotrophic oceans. However, WH8102 possesses three distinct PstS homologs whose substrate specificity and ecological roles are unclear. The three PstS homologs were heterologously expressed and purified to investigate their substrate specificity and binding kinetics. Our study revealed that all three PstS homologs exhibited a high degree of specificity for phosphate but differed in phosphate binding affinities. Notably, PstS1b displayed nearly 10-fold higher binding affinity (KD = 0.44 µM) compared to PstS1a (KD = 3.3 μM) and PstS2 (KD = 4.3 μM). Structural modeling suggested a single amino acid variation in the binding pocket of PstS1b (threonine instead of serine in PstS1a and PstS2) likely contributed to its higher Pi affinity. Genome context data, together with the protein biophysical data, suggest distinct ecological roles for the three PstS homologs. We propose that PstS1b may be involved in scavenging inorganic phosphorus in oligotrophic conditions and that PstS1a may be involved in transporting recycled phosphate derived from organic phosphate cleavage. The role of PstS2 is less clear, but it may be involved in phosphate uptake when environmental phosphate concentrations are transiently higher. The conservation of three distinct PstS homologs in Synechococcus clade III strains likely reflects distinct adaptations for P acquisition under varying oligotrophic conditions.IMPORTANCEPhosphorus is an essential macronutrient that plays a key role in marine primary productivity and biogeochemistry. However, intense competition for bioavailable phosphorus in the marine environment limits growth and productivity of ecologically important cyanobacteria. In oligotrophic oceans, marine Synechococcus strains, like WH8102, utilize high-affinity phosphate-binding proteins (PstS) to scavenge inorganic phosphate. However, WH8102 possesses three distinct PstS homologs, with unclear substrate specificity and ecological roles, creating a knowledge gap in understanding phosphorus acquisition mechanisms in picocyanobacteria. Through genomic, functional, biophysical, and structural analysis, our study unravels the ecological functions of these homologs. Our findings enhance our understanding of cyanobacterial nutritional uptake strategies and shed light on the crucial role of these conserved nutrient uptake systems in adaptation to specific niches, which ultimately underpins the success of marine Synechococcus across a diverse array of marine ecosystems.

Seasonal and spatial variations of Synechococcus in abundance, pigment types, and genetic diversity in a temperate semi-enclosed bay

Synechococcus is abundant and globally widespread in various marine environments. Seasonal and spatial variations in Synechococcus abundance, pigment types, and genetic diversity were investigated based on flow cytometric analysis and high-throughput sequencing of cpcBA operon (encoding phycocyanin) and rpoC1 gene (encoding RNA polymerase) in a temperate semi-enclosed bay. Synechococcus abundance exhibited seasonal variations with the highest value in summer and the lowest value in winter, which was consistent with temperature variation. Three pigment types of Synechococcus type 1, type 2, and type 3 were distinguished based on cpcBA operon, which displayed obvious variations spatially between the inner and the outer bay. Freshwater discharge and water turbidity played important roles in regulating Synechococcus pigment types. Synechococcus assemblages were phylogenetically diverse (12 different lineages) based on rpoC1 gene and dominated by three core lineages S5.1-I, S5.1-IX, and S5.2-CB5 in different seasons. Our study demonstrated that Synechococcus abundance, pigment types, and genetic diversity displayed variations seasonally and spatially by different techniques, which were mainly driven by temperature, salinity, nutrients, and turbidity. The combination of more technical means provides more information for studying Synechococcus distribution. In this study, three pigment types of Synechococcus were discriminated simultaneously by dual lasers flow cytometer for the first time.

Seasonal Patterns of Picocyanobacterial Community Structure in the Kuroshio Current

The nutrient-scarce, warm, and high-salinity Kuroshio current has a profound impact on both the marine ecology of the northwestern Pacific Ocean and the global climate. This study aims to reveal the seasonal dynamics of picoplankton in the subtropical Kuroshio current. Our results showed that one of the picocyanobacteria, Synechococcus, mainly distributed in the surface water layer regardless of seasonal changes, and the cell abundance ranged from 104 to 105 cells mL-1. In contrast, the maximum concentration of the other picocyanobacteria, Prochlorococcus, was maintained at more than 105 cells mL-1 throughout the year. In the summer and the autumn, Prochlorococcus were mainly concentrated at the water layer near the bottom of the euphotic zone. They were evenly distributed in the euphotic zone in the spring and winter. The stirring effect caused by the monsoon determined their distribution in the water column. In addition, the results of 16S rRNA gene diversity analysis showed that the seasonal changes in the relative abundance of Synechococcus and Prochlorococcus in the surface water of each station accounted for 20 to 40% of the total reads. The clade II of Synechococcus and the High-light II of Prochlorococcus were the dominant strains in the waters all year round. Regarding other picoplankton, Proteobacteria and Actinobacteria occupied 45% and 10% of the total picoplankton in the four seasons. These data should be helpful for elucidating the impacts of global climate changes on marine ecology and biogeochemical cycles in the Western Boundary Currents in the future.

Network pharmacology‒based analysis of marine cyanobacteria derived bioactive compounds for application to Alzheimer's disease

In recent years, the Alzheimer's disease (AD) epidemic has become one of the largest global healthcare crises. Besides, the available systemic therapies for AD are still inadequate. Due to the insufficient therapeutic options, new treatment strategies are urgently needed to achieve a satisfactory therapeutic effect. Marine bio-resources have been accepted as one of the most economically viable and sustainable sources with potential applications for drug discovery and development. In this study, a marine cyanobacteria-Synechococcus sp. XM-24 was selected as the object of research, to systematically investigate its therapeutic potential mechanisms for AD. The major active compounds derived from the Synechococcus sp. biomass were identified via pyrolysis-gas chromatography-mass spectrometry (GC-MS), and 22 compounds were identified in this strain. The most abundant chemical compounds was (E)-octadec-11-enoic acid, with the peak area of 30.6%. Follow by tridecanoic acid, 12-methyl- and hexadecanoic acid, with a peak area of 23.26% and 18.23%, respectively. GC-MS analysis also identified indolizine, isoquinoline, 3,4-dihydro- and Phthalazine, 1-methyl-, as well as alkene and alkane from the strain. After the chemical toxicity test, 10 compounds were finally collected to do the further analysis. Then, network pharmacology and molecular docking were adopted to systematically study the potential anti-AD mechanism of these compounds. Based on the analysis, the 10 Synechococcus-derived active compounds could interact with 128 related anti-AD targets. Among them, epidermal growth factor receptor (EGFR), vascular endothelial growth factor A (VEGFA) and mitogen-activated protein kinase 3 (MAPK3) were the major targets. Furthermore, the compounds N-capric acid isopropyl ester, (E)-octadec-11-enoic acid, and 2H-Pyran-2,4(3H)-dione, dihydro-6-methyl- obtained higher degrees in the compounds-intersection targets network analysis, indicating these compounds may play more important role in the process of anti-AD. In addition, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these active compounds exert the anti-AD effects mainly through PI3K-Akt signaling pathway, neuroactive ligand-receptor interaction and ras signaling pathway. Our study identified Synechococcus-derived bioactive compounds have the potential for application to AD by targeting multiple targets and related pathways, which will provide a foundation for future research on applications of marine cyanobacteria in the functional drug industry.

Spatio-temporal patterns of Synechococcus oligotypes in Moroccan lagoonal environments

Synechococcus are unicellular cyanobacteria susceptible to environmental fluctuations and can be used as bioindicators of eutrophication in marine ecosystems. We examined their distribution in two Moroccan lagoons, Marchica on the Mediterranean coast and Oualidia on the Atlantic, in the summers of 2014 and 2015 using 16S rRNA amplicon oligotyping. Synechococcus representatives recruited a higher number of reads from the 16S rRNA in Marchica in comparison to Oualidia. We identified 31 Synechococcus oligotypes that clustered into 10 clades with different distribution patterns. The Synechococcus community was mainly represented by oligotype 1 (clade III) in Marchica. Cooccurring clades IV and I had an important relative abundance in Marchica in the summer of 2014, which is unusual, as these clades are widespread in cold waters. Moreover, Clades VII and subcluster "5.3" formed a sizeable percentage of the Synechococcus community in Marchica. Notably, we found low Synechococcus sequence counts in the Atlantic Lagoon. These results showed that the relative abundance of Synechococcus reads is not constant over space and time and that rare members of the Synechococcus community did not follow a consistent pattern. Further studies are required to decipher Synechococcus dynamics and the impact of environmental parameters on their spatial and temporal distributions.

Functional flexibility of cyanobacterial light harvesting phycobilisomes enable acclimation to the complex light regime of mixing marine water columns

The light environment in a mixing water column is arguably the most erratic condition under which photosynthesis functions. Shifts in light intensity, by an order of magnitude, can occur over the time scale of hours. In marine Synechococcus, light is harvested by massive, membrane attached, phycobilisome chromophore-protein complexes (PBS). We examined the ability of a phycobilisome-containing marine Synechococcus strain (WH8102) to acclimate to illumination perturbations on this scale. Although changes in pigment composition occurred gradually over the course of days, we did observe significant and reversible changes in the pigment's fluorescence emission spectra on a time scale of hours. Upon transition to ten-fold higher intensities, we observed a decrease in the energy transferred to Photosystem II. At the same time, the spectral composition of PBS fluorescence emission shifted. Unlike fluorescence quenching mechanisms, this phenomenon resulted in increased fluorescence intensities. These data suggest a mechanism by which marine Synechococcus WH8102 detaches hexamers from the phycobilisome structure. The fluorescence yield of these uncoupled hexamers is high. The detachment process does not require protein synthesis as opposed to reattachment. Hence, the most likely process would be the degradation and resynthesis of labile PBS linker proteins. Experiments with additional species yielded similar results, suggesting that this novel mechanism might be broadly used among PBS-containing organisms.

Quantification of Marine Picocyanobacteria on Water Column Particles and in Sediments Using Real-Time PCR Reveals Their Role in Carbon Export

Picocyanobacteria are the most abundant primary producers in the ocean and play a fundamental role in marine carbon cycling. Quantification of picocyanobacteria on sinking particles and in sediments is essential to understanding their contribution to the biological carbon pump. We designed a primer set targeting the 16S-23S rRNA internal transcribed spacer (ITS) sequence of cyanobacteria and established a quantitative PCR (qPCR) method for quantifying the ITS sequence abundance. High-throughput sequencing confirmed that this primer set can cover broad diversities of marine picocyanobacteria and avoid amplification of other marine cyanobacteria such as Trichodesmium and Crocosphaera. Amplification efficiencies were slightly different when seven marine Synechococcus and Prochlorococcus strains were assayed. The qPCR results were comparable with flow cytometry for water samples. Using this method, we found that, in the dark ocean, picocyanobacterial ITS sequence abundances were 10 to 100 copies/mL in the size fraction of 0.2 to 3 μm, which were 1 to 3 orders of magnitude more abundant than on the >3-μm particles. We also found that picocyanobacterial ITS abundance in sediment ranged from 10⁵ to 10⁷ copies/g along two nearshore-to-offshore transects in the northern South China Sea. These results further explain the important role of picocyanobacteria in carbon export. Collectively, we provide a qPCR method quantifying the total abundance of marine picocyanobacteria on water column particles and in sediments. Moreover, this newly designed primer set can be also applied to investigate the community of picocyanobacteria via high-throughput sequencing. IMPORTANCE Picocyanobacteria are the most abundant primary producers in the ocean. However, quantification of picocyanobacteria on the sinking particles and in sediments remains challenging using flow cytometry or epifluorescence microscopy. Here, we developed a real-time PCR method to quantify picocyanobacteria using a newly designed primer set specifically targeting the 16S-23S rRNA ITS sequence of cyanobacteria. We showed that in the dark ocean, picocyanobacteria are 1 to 3 orders of magnitude more abundant in small particles (0.2 to 3 μm) than in larger particles (>3 μm). This result supports the important role of direct sinking free-living picocyanobacteria cells in the carbon export to deep ocean. We also found that the picocyanobacterial ITS sequence abundance were 10⁵ to 10⁷ copies per gram in sediments, suggesting significant accumulation of sinking picocyanobacteria in the benthic ecosystem. This qPCR method can be used to quantify the contribution of picocyanobacteria to the biological carbon pump.

Distributional Pattern of Bacteria, Protists, and Diatoms in Ocean according to Water Depth in the Northern South China Sea

Ocean microbiomes provide insightful details about the condition of water and the global impact of marine ecosystems. A fine-scale analysis of ocean microbes may shed light on the dynamics and function of the ocean microbiome community. In this study, we evaluated the changes in the community and function of marine bacteria, protists, and diatoms corresponding to different ocean depths using next-generation sequencing methods. We found that diatoms displayed a potential water-depth pattern in species richness (alpha diversity) and community composition (beta diversity). However, for bacteria and protists, there was no significant relationship between water depth and species richness. This may be related to the biological characteristics of diatoms. The photosynthesis of diatoms and their distribution may be associated with the fluctuating light regime in the underwater climate. Moreover, salinity displayed negative effects on the abundance of some diatom and bacterial groups, which indicates that salinity may be one of the factors restricting ocean microorganism diversity. In addition, compared to the global ocean microbiome composition, function, and antibiotic resistance genes, a water depth pattern due to the fine-scale region was not observed in this study. IMPORTANCE Fine-scale analysis of ocean microbes provides insights into the dynamics and functions of the ocean microbiome community. Here, using amplicon and metagenome sequencing methods, we found that diatoms in the northern South China Sea displayed a potential water-depth pattern in species richness and community composition, which may be related to their biological characteristics. The potential effects of the differences in geographic sites mainly occurred in the diatom and bacterial communities. Moreover, given the correlation between the environmental factors and relative abundance of antibiotic resistance genes (ARGs), the study of ocean ARG distribution patterns should integrate the potential effects of environmental factors.

Global Phylogeography of Marine Synechococcus in Coastal Areas Reveals Strong Community Shifts

Marine Synechococcus comprise a numerically and ecologically prominent phytoplankton group, playing a major role in both carbon cycling and trophic networks in all oceanic regions except in the polar oceans. Despite their high abundance in coastal areas, our knowledge of Synechococcus communities in these environments is based on only a few local studies. Here, we use the global metagenome data set of the Ocean Sampling Day (June 21^st, 2014) to get a snapshot of the taxonomic composition of coastal Synechococcus communities worldwide, by recruitment on a reference database of 141 picocyanobacterial genomes, representative of the whole Prochlorococcus, Synechococcus, and Cyanobium diversity. This allowed us to unravel drastic community shifts over small to medium scale gradients of environmental factors, in particular along European coasts. The combined analysis of the phylogeography of natural populations and the thermophysiological characterization of eight strains, representative of the four major Synechococcus lineages (clades I to IV), also brought novel insights about the differential niche partitioning of clades I and IV, which most often co-dominate the Synechococcus community in cold and temperate coastal areas. Altogether, this study reveals several important characteristics and specificities of the coastal communities of Synechococcus worldwide. IMPORTANCE Synechococcus is the second most abundant phytoplanktonic organism on Earth, and its wide genetic diversity allowed it to colonize all the oceans except for polar waters, with different clades colonizing distinct oceanic niches. In recent years, the use of global metagenomics data sets has greatly improved our knowledge of "who is where" by describing the distribution of Synechococcus clades or ecotypes in the open ocean. However, little is known about the global distribution of Synechococcus ecotypes in coastal areas, where Synechococcus is often the dominant phytoplanktonic organism. Here, we leverage the global Ocean Sampling Day metagenomics data set to describe Synechococcus community composition in coastal areas worldwide, revealing striking community shifts, in particular along the coasts of Europe. As temperature appears as an important driver of the community composition, we also characterize the thermal preferenda of 8 Synechococcus strains, bringing new insights into the adaptation to temperature of the dominant Synechococcus clades.

Antibiotic Resistomes and Microbiomes in the Surface Water along the Code River in Indonesia Reflect Drainage Basin Anthropogenic Activities

Water and sanitation are important factors in the emergence of antimicrobial resistance in low- and middle-income countries. Drug residues, metals, and various wastes foster the spread of antibiotic resistance genes (ARGs) with the help of mobile genetic elements (MGEs), and therefore, rivers receiving contaminants and effluents from multiple sources are of special interest. We followed both the microbiome and resistome of the Code River in Indonesia from its pristine origin at the Merapi volcano through rural and then city areas to the coast of the Indian Ocean. We used a SmartChip quantitative PCR with 382 primer pairs for profiling the resistome and MGEs and 16S rRNA gene amplicon sequencing to analyze the bacterial communities. The community structure explained the resistome composition in rural areas, while the city sampling sites had lower bacterial diversity and more ARGs, which correlated with MGEs, suggesting increased mobility potential in response to pressures from human activities. Importantly, the vast majority of ARGs and MGEs were no longer detectable in marine waters at the ocean entrance. Our work provides information on the impact of different influents on river health as well as sheds light on how land use contributes to the river resistome and microbiome.

Seasonal dynamics in picocyanobacterial abundance and clade composition at coastal and offshore stations in the Baltic Sea

Picocyanobacteria (< 2 µm in diameter) are significant contributors to total phytoplankton biomass. Due to the high diversity within this group, their seasonal dynamics and relationship with environmental parameters, especially in brackish waters, are largely unknown. In this study, the abundance and community composition of phycoerythrin rich picocyanobacteria (PE-SYN) and phycocyanin rich picocyanobacteria (PC-SYN) were monitored at a coastal (K-station) and at an offshore station (LMO; ~ 10 km from land) in the Baltic Sea over three years (2018–2020). Cell abundances of picocyanobacteria correlated positively to temperature and negatively to nitrate (NO₃) concentration. While PE-SYN abundance correlated to the presence of nitrogen fixers, PC-SYN abundance was linked to stratification/shallow waters. The picocyanobacterial targeted amplicon sequencing revealed an unprecedented diversity of 2169 picocyanobacterial amplicons sequence variants (ASVs). A unique assemblage of distinct picocyanobacterial clades across seasons was identified. Clade A/B dominated the picocyanobacterial community, except during summer when low NO_3, high phosphate (PO₄) concentrations and warm temperatures promoted S5.2 dominance. This study, providing multiyear data, links picocyanobacterial populations to environmental parameters. The difference in the response of the two functional groups and clades underscore the need for further high-resolution studies to understand their role in the ecosystem.

Effect of Rising Temperature and Carbon Dioxide on the Growth, Photophysiology, and Elemental Ratios of Marine Synechococcus: A Multistressor Approach

Marine picocyanobacteria belonging to the genus Synechococcus are one of the most abundant photosynthetic organisms on Earth. They are often exposed to large fluctuations in temperature and CO2 concentrations in the ocean, which are expected to further change in the coming decades due to ocean acidification and warming resulting from rising atmospheric CO2 levels. To decipher the effect of changing temperature and CO2 levels on Synechococcus, six Synechococcus strains previously isolated from various coastal and open ocean sites were exposed to a matrix of three different temperatures (22 °C, 24 °C and 26 °C) and CO2 levels (400 ppm, 600 ppm and 800 ppm). Thereafter, the specific growth rates, photophysiological parameters (σPSII and Fv/Fm), C/N (mol/mol) ratios and the nitrogen stable isotopic composition (δ15N (‰)) of the strains were measured. Temperature was found to be a stronger driver of the changes in specific growth rates and photophysiology in the Synechococcus strains. Carbon-concentrating mechanisms (CCM) operational in these strains that shield the photosynthetic machinery from directly sensing ambient changes in CO2 possibly played a major role in causing minimal changes in the specific growth rates under the varying CO2 levels.

Seasonal and Spatial Variations in Synechococcus Abundance and Diversity Throughout the Gullmar Fjord, Swedish Skagerrak

The picophytoplankton Synechococcus is a globally abundant autotroph that contributes significantly to primary production in the oceans and coastal areas. These cyanobacteria constitute a diverse genus of organisms that have developed independent niche spaces throughout aquatic environments. Here, we use the 16S V3–V4 rRNA gene region and flow cytometry to explore the diversity of Synechococcus within the picophytoplankton community in the Gullmar Fjord, on the west coast of Sweden. We conducted a station-based 1-year time series and two transect studies of the fjord. Our analysis revealed that within the large number of Synechococcus amplicon sequence variants (ASVs; 239 in total), prevalent ASVs phylogenetically clustered with clade representatives in both marine subcluster 5.1 and 5.2. The near-surface composition of ASVs shifted from spring to summer, when a 5.1 subcluster dominated community developed along with elevated Synechococcus abundances up to 9.3 × 104 cells ml–1. This seasonal dominance by subcluster 5.1 was observed over the length of the fjord (25 km), where shifts in community composition were associated with increasing depth. Unexpectedly, the community shift was not associated with changes in salinity. Synechococcus abundance dynamics also differed from that of the photosynthetic picoeukaryote community. These results highlight how seasonal variations in environmental conditions influence the dynamics of Synechococcus clades in a high latitude threshold fjord.

Picocyanobacterial Synechococcus in marine ecosystem: Insights from genetic diversity, global distribution, and potential function

Marine Synechococcus, a main group of picocyanobacteria, has been ubiquitously observed across the global oceans. Synechococcus exhibits high phylogenetical and phenotypical diversity, and horizontal gene transfer makes its genetic evolution much more intricate. With the development of measurement technologies and analysis methods, the genomic information and niche partition of each Synechococcus lineage tend to be precisely described, but the global analysis is still lacking. Therefore, it is necessary to summarize existing studies and integrate published data to gain a comprehensive understanding of Synechococcus on genetic variation, niche division, and potential functions. In this review, the maximum likelihood trees are constructed based on existing sequence data, including both phylogenetic and pigmentary gene markers. The global distribution characteristics of abundance, lineages, and pigment types are concluded through pooled analysis of more than 700 samples obtained from approximately 50 scientific research cruises. The potential functions of Synechococcus are explored in element cycles and biological interactions. Future work on Synechococcus is suggested to focus on not only elucidating the nature of Synechococcus biodiversity but also demonstrating its interactions with the ecosystem by combining bioinformatics and macroscopic isotope-labeled environmental parameters.

Phylogenetic Structure of Synechococcus Assemblages and Its Environmental Determinants in the Bay and Strait Areas of a Continental Sea

Marine Synechococcus, a significant contributor to primary production, shows high phylogenetic diversity. However, studies on its phylogenetic composition in the Bohai Sea, the largest continental sea in China, are lacking. We sequenced rpoC1 (encodes the RNA polymerase β’ subunit protein) in samples from the Laizhou Bay (LZB) and Bohai Strait (BS) in June and November using high-throughput sequencing to reveal the phylogenetic composition of Synechococcus assemblages in the bay and strait areas of the Bohai Sea. In total, 12 lineages representing Synechococcus subclusters S5.1, S5.2, and S5.3 were identified. Spatially, clade I was obligately dominant in BS. In contrast, the Synechococcus assemblage in LZB was more diverse, with clades VI and III being highly abundant. In addition, we detected strong variation in Synechococcus structure between June and November in the Bohai Sea. Clades II, III, XX, and miyav were only detected in November. Vertically, variation in Synechococcus assemblage was not apparent among the water layers probably due to the shallow water depth with intense water mixing. Results of redundancy analysis (RDA) and random forest (RF) analysis together highlighted the key role of silicate in the Synechococcus assemblage. Our results suggested that the Bohai Sea provides various niches for different Synechococcus clades, resulting in a special phylogenetic composition of the Synechococcus assemblage, compared with that in the adjacent shelf sea and other continental seas in the world.

Spatio-Temporal Variation of Synechococcus Assemblages at DNA and cDNA Levels in the Tropical Estuarine and Coastal Waters

Synechococcus is a major contributor to global marine primary production. Here, its spatio-temporal variations in abundance and phylogenetic structure were studied at three stations of the South China Sea at both DNA and cDNA levels. Synechococcus cell abundance was lowest in March, but highest in October at two coastal stations. Its abundance was higher at the estuarine station, which reached a peak value of 1.36 × 105 cells/ml in April, owing to the nitrogen nutrients discharged from the Sanya River. Gene and gene transcript abundances of four Synechococcus lineages, clades II, III, VIII, and S5.3, were studied by quantitative PCR, which showed that clade II was the most abundant lineage at both DNA and cDNA levels. High-throughput sequencing revealed that, at the DNA level, Synechococcus assemblage was dominated by clade SY4 (a novel clade defined in this study), S5.2, and clade II in the coastal waters and was dominated by freshwater/S5.2 Synechococcus, reaching a value up to 88.61% in June, in estuarine waters. Changes in salinity and nutrient concentration caused by seasonal monsoonal forcing and river discharge were the key determinants of the spatio-temporal variation in Synechococcus assemblages at the DNA level. In comparison, high dissimilation among samples at the same stations and in the same seasons leads to the imperceptible spatio-temporal variation pattern of Synechococcus assemblages at the cDNA level. Furthermore, co-occurrence networks disclosed that Synechococcus community had closer and more complex internal interactions at the cDNA level. These discrepancies highlighted the necessity to study Synechococcus assemblages at both DNA and cDNA levels.

The "Dark Side" of Picocyanobacteria: Life as We Do Not Know It (Yet)

Picocyanobacteria of the genus Synechococcus (together with Cyanobium and Prochlorococcus) have captured the attention of microbial ecologists since their description in the 1970s. These pico-sized microorganisms are ubiquitous in aquatic environments and are known to be some of the most ancient and adaptable primary producers. Yet, it was only recently, and thanks to developments in molecular biology and in the understanding of gene sequences and genomes, that we could shed light on the depth of the connection between their evolution and the history of life on the planet. Here, we briefly review the current understanding of these small prokaryotic cells, from their physiological features to their role and dynamics in different aquatic environments, focussing particularly on the still poorly understood ability of picocyanobacteria to adapt to dark conditions. While the recent discovery of Synechococcus strains able to survive in the deep Black Sea highlights how adaptable picocyanobacteria can be, it also raises more questions-showing how much we still do not know about microbial life. Using available information from brackish Black Sea strains able to perform and survive in dark (anoxic) conditions, we illustrate how adaptation to narrow ecological niches interacts with gene evolution and metabolic capacity.

Phylogenetic and Phenogenetic Diversity of Synechococcus along a Yellow Sea Section Reveal Its Environmental Dependent Distribution and Co-Occurrence Microbial Pattern

Synechococcus is a dominant genus of the coastal phytoplankton with an effective contribution to primary productivity. Here, the phylogenetic and phenogenetic composition of Synechococcus in the coastal Yellow Sea was addressed by sequencing marker gene methods. Meanwhile, its co-occurrence pattern with bacterial and eukaryotic microbes was further investigated based on the construction of networks. The result revealed that Synechococcus abundance ranged from 9.8 × 102 cells mL−1 to 1.6 × 105 cells mL−1, which was significantly correlated to sampling depth and nutrient contents of nitrite, ammonia, and dissolved silicon. A total of eight Synechococcus phylogenetic lineages were detected, of which clade III was dominant in most of the samples. Meanwhile, clade I increased along the water column and even reached a maximum value of 76.13% at 20 m of station B. Phenogenetically, Synechococcus PT3 was always the predominant pigment type across the whole study zone. Only salinity was significantly correlated to the phenogenetic constitution. The networks revealed that Synechococcus co-occurred with 159 prokaryotes, as well as 102 eukaryotes including such possible grazers as Gymnodinium clades and Alveolata. Potential function prediction further showed that microbes co-occurring with Synechococcus were associated with diverse element cycles, but the exact mechanism needed further experimentation to verify. This research promotes exploring regularity in the genomic composition and niche position of Synechococcus in the coastal ecosystem and is significant to further discuss its potential participation in materials circulation and bottom-up effects in microbial food webs.

Seasonal environmental variability drives microdiversity within a coastal Synechococcus population

Marine microbes often show a high degree of physiological or ecological diversity below the species level. This microdiversity raises questions about the processes that drive diversification and permit coexistence of diverse yet closely related marine microbes, especially given the theoretical efficiency of competitive exclusion. Here, we provide insight with an 8‐year time series of diversity within Synechococcus, a widespread and important marine picophytoplankter. The population of Synechococcus on the Northeast U.S. Shelf is comprised of six main types, each of which displays a distinct and consistent seasonal pattern. With compositional data analysis, we show that these patterns can be reproduced with a simple model that couples differential responses to temperature and light with the seasonal cycle of the physical environment. These observations support the hypothesis that temporal variability in environmental factors can maintain microdiversity in marine microbial populations. We also identify how seasonal diversity patterns directly determine overarching Synechococcus population abundance features.

Dynamics and Distribution of Marine Synechococcus Abundance and Genotypes during Seasonal Hypoxia in a Coastal Marine Ranch

Marine Synechococcus are an ecologically important picocyanobacterial group widely distributed in various oceanic environments. Little is known about the dynamics and distribution of Synechococcus abundance and genotypes during seasonal hypoxia in coastal zones. In this study, an investigation was conducted in a coastal marine ranch along two transects in Muping, Yantai, where hypoxic events (defined here as the dissolved oxygen concentration <3 mg L−1) occurred in the summer of 2015. The hypoxia occurred in the bottom waters from late July and persisted until late August. It was confined at nearshore stations of the two transects, one running across a coastal ranch and the other one outside. During this survey, cell abundance of Synechococcus was determined with flow cytometry, showing great variations ranging from 1 × 104 to 3.0 × 105 cells mL−1, and a bloom of Synechococcus occurred when stratification disappeared and hypoxia faded out outside the ranch. Regression analysis indicated that dissolved oxygen, pH, and inorganic nutrients were the most important abiotic factors in explaining the variation in Synechococcus cell abundance. Diverse genotypes (mostly belonged to the sub-clusters 5.1 and 5.2) were detected using clone library sequencing and terminal restriction fragment length polymorphism analysis of the 16S–23S rRNA internal transcribed spacer region. The richness of genotypes was significantly related to salinity, temperature, silicate, and pH, but not dissolved oxygen. Two environmental factors, temperature and salinity, collectively explained 17% of the variation in Synechococcus genotype assemblage. With the changes in population composition in diverse genotypes, the Synechococcus assemblages survived in the coastal hypoxia event and thrived when hypoxia faded out.

Recurrent microbial community types driven by nearshore and seasonal processes in coastal Southern California

A multitude of concurrent biological and physical processes contribute to microbial community turnover, especially in highly dynamic coastal environments. Characterizing what factors contribute most to shifts in microbial community structure and the specific organisms that correlate with changes in the products of photosynthesis improves our understanding of nearshore microbial ecosystem functions. We conducted high frequency sampling in nearshore Southern California in order to capture sub-weekly microbial community dynamics. Microbial communities were characterized by flow cytometry and 16S rRNA gene sequencing, and placed in the context of physicochemical parameters. Within our time-series, season and nutrient availability corresponded to changes in dominant microbial community members. Concurrent aseasonal drivers with overlapping scales of variability were also apparent when we used network analysis to assess the microbial community as subsets of the whole. Our analyses revealed the microbial community as a mosaic, with overlapping groups of taxa that varied on different timescales and correlated with unique abiotic and biotic factors. Specifically, a subnetwork associated with chlorophyll a exhibited rapid turnover, indicating that ecologically important subsets of the microbial community can change on timescales different than and in response to factors other than those that govern turnover of most members of the assemblage.

Spatial and temporal variations in Synechococcus microdiversity in the Southern California coastal ecosystem

The Synechococcus cyanobacterial population at the Scripps Institution of Oceanography pier in La Jolla, CA, shows large increases in abundance, typically in the spring and summer followed, by rapid declines within weeks. Here we used amplicon sequencing of the ribosomal RNA internal transcribed spacer region to examine the microdiversity within this cyanobacterial genus during these blooms as well as further offshore in the Southern California coastal ecosystem (CCE). These analyses revealed numerous Synechococcus amplicon sequence variants (ASVs) and that clade and ASV composition can change over the course of blooms. We also found that a large bloom in August 2016 was highly anomalous both in its overall Synechococcus abundance and in terms of the presence of normally oligotrophic Synechococcus clade II. The dominant ASVs at the pier were found further offshore and in the California Current, but we did observe more oligotrophic ASVs and clades along with depth variation in Synechococcus diversity. We also observed that the dominant sequence variant switched during the peak of multiple Synechococcus blooms, with this switch occurring in multiple clades, but we present initial evidence that this apparent ASV switch is a physiological response rather than a change in the dominant population.

Spatial Distribution Patterns of Bacterioplankton in the Oxygen Minimum Zone of the Tropical Mexican Pacific

Microbial communities within oxygen minimum zones (OMZs) are crucial drivers of marine biogeochemical cycles; however, we still lack an understanding of how these communities are distributed across an OMZ. We explored vertical (from 5 to 500 m depth) and horizontal (coast to open ocean) distribution of bacterioplankton and its relationships with the main oceanographic conditions in three transects of the tropical Mexican Pacific OMZ. The distribution of the microbial diversity and the main clades changed along the transition from oxygen-rich surface water to the OMZ core, demonstrating the sensitivity of key bacterial groups to deoxygenation. The euphotic zone was dominated by Synechococcales, followed by Flavobacteriales, Verrucomicrobiales, Rhodobacterales, SAR86, and Cellvibrionales, whereas the OMZ core was dominated by SAR11, followed by SAR406, SAR324, SAR202, UBA10353 marine group, Thiomicrospirales and Nitrospinales. The marked environmental gradients along the water column also supported a high potential for niche partitioning among OMZ microorganisms. Additionally, in the OMZ core, bacterial assemblages from the same water mass were more similar to each other than those from another water mass. There were also important differences between coastal and open-ocean communities: Flavobacteriales, Verrucomicrobiales, Rhodobacterales, SAR86, and Cellvibrionales were more abundant in coastal areas, while Synechococcales, SAR406, SAR324, SAR202, UBA10353 marine group, and Thiomicrospirales were more abundant in the open ocean. Our results suggest a biogeographic structure of the bacterioplankton in this OMZ region, with limited community mixing across water masses, except in upwelling events, and little dispersion of the community by currents in the euphotic zone.

Persistent El Niño driven shifts in marine cyanobacteria populations

In the California Current Ecosystem, El Niño acts as a natural phenomenon that is partially representative of climate change impacts on marine bacteria at timescales relevant to microbial communities. Between 2014–2016, the North Pacific warm anomaly (a.k.a., the “blob”) and an El Niño event resulted in prolonged ocean warming in the Southern California Bight (SCB). To determine whether this “marine heatwave” resulted in shifts in microbial populations, we sequenced the rpoC1 gene from the biogeochemically important picocyanobacteria Prochlorococcus and Synechococcus at 434 time points from 2009–2018 in the MICRO time series at Newport Beach, CA. Across the time series, we observed an increase in the abundance of Prochlorococcus relative to Synechococcus as well as elevated frequencies of ecotypes commonly associated with low-nutrient and high-temperature conditions. The relationships between environmental and ecotype trends appeared to operate on differing temporal scales. In contrast to ecotype trends, most microdiverse populations were static and possibly reflect local habitat conditions. The only exceptions were microdiversity from Prochlorococcous HLI and Synechococcus Clade II that shifted in response to the 2015 El Niño event. Overall, Prochlorococcus and Synechococcus populations did not return to their pre-heatwave composition by the end of this study. This research demonstrates that extended warming in the SCB can result in persistent changes in key microbial populations.

Genomic Comparison and Spatial Distribution of Different Synechococcus Phylotypes in the Black Sea

Picocyanobacteria of the genus Synechococcus are major contributors to global primary production and nutrient cycles due to their oxygenic photoautotrophy, their abundance, and the extensive distribution made possible by their wide-ranging biochemical capabilities. The recent recovery and isolation of strains from the deep euxinic waters of the Black Sea encouraged us to expand our analysis of their adaptability also beyond the photic zone of aquatic environments. To this end, we quantified the total abundance and distribution of Synechococcus along the whole vertical profile of the Black Sea by flow cytometry, and analyzed the data obtained in light of key environmental factors. Furthermore, we designed phylotype-specific primers using the genomes of two new epipelagic coastal strains – first described here – and of two previously described mesopelagic strains, analyzed their presence/abundance by qPCR, and tested this parameter also in metagenomes from two stations at different depths. Together, whole genome sequencing, metagenomics and qPCR techniques provide us with a higher resolution of Synechococcus dynamics in the Black Sea. Both phylotypes analyzed are abundant and successful in epipelagic coastal waters; but while the newly described epipelagic strains are specifically adapted to this environment, the strains previously isolated in mesopelagic waters are able, in low numbers, to withstand the aphotic and oxygen depleted conditions of deep layers. This heterogeneity allows different Synechococcus phylotypes to occupy different niches and underscores the importance of a more detailed characterization of the abundance, distribution, and dynamics of individual populations of these picocyanobacteria.

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.

Comparative membrane proteomics reveal contrasting adaptation strategies for coastal and oceanic marine Synechococcus cyanobacteria

Marine cyanobacteria genus Synechococcus are among the most abundant and widespread primary producers in the open ocean. Synechococcus strains belonging to different clades have adapted distinct strategies for growth and survival across a range of marine conditions. Clades I and IV are prevalent in colder, mesotrophic, coastal waters, while clades II and III prefer warm, oligotrophic open oceans. To gain insight into the cellular resources these unicellular organisms invest in adaptation strategies we performed shotgun membrane proteomics of four Synechococcus spp. strains namely CC9311 (clade I), CC9605 (clade II), WH8102 (clade III) and CC9902 (clade IV). Comparative membrane proteomes analysis demonstrated that CC9902 and WH8102 showed high resource allocation for phosphate uptake, accounting for 44% and 38% of overall transporter protein expression of the species. WH8102 showed high expression of the iron uptake ATP-binding cassette binding protein FutA, suggesting that a high binding affinity for iron is possibly a key adaptation strategy for some strains in oligotrophic ocean environments. One protein annotated as a phosphatase 2c (Sync_2505 and Syncc9902_0387) was highly expressed in the coastal mesotrophic strains CC9311 and CC9902, constituting 14%-16% of total membrane protein, indicating a vital, but undefined function, for strains living in temperate mesotrophic environments.

Multi-year dynamics of fine-scale marine cyanobacterial populations are more strongly explained by phage interactions than abiotic, bottom-up factors

Currently defined ecotypes in marine cyanobacteria Prochlorococcus and Synechococcus likely contain subpopulations that themselves are ecologically distinct. We developed and applied high-throughput sequencing for the 16S-23S rRNA internally transcribed spacer (ITS) to examine ecotype and fine-scale genotypic community dynamics for monthly surface water samples spanning 5 years at the San Pedro Ocean Time-series site. Ecotype-level structure displayed regular seasonal patterns including succession, consistent with strong forcing by seasonally varying abiotic parameters (e.g. temperature, nutrients, light). We identified tens to thousands of amplicon sequence variants (ASVs) within ecotypes, many of which exhibited distinct patterns over time, suggesting ecologically distinct populations within ecotypes. Community structure within some ecotypes exhibited regular, seasonal patterns, but not for others, indicating other more irregular processes such as phage interactions are important. Network analysis including T4-like phage genotypic data revealed distinct viral variants correlated with different groups of cyanobacterial ASVs including time-lagged predator-prey relationships. Variation partitioning analysis indicated that phage community structure more strongly explains cyanobacterial community structure at the ASV level than the abiotic environmental factors. These results support a hierarchical model whereby abiotic environmental factors more strongly shape niche partitioning at the broader ecotype level while phage interactions are more important in shaping community structure of fine-scale variants within ecotypes.

Primer Design for an Accurate View of Picocyanobacterial Community Structure by Using High-Throughput Sequencing

High-throughput sequencing (HTS) of the 16S rRNA gene has been used successfully to describe the structure and dynamics of microbial communities. Picocyanobacteria are important members of bacterioplankton communities, and, so far, they have predominantly been targeted using universal bacterial primers, providing a limited resolution of the picocyanobacterial community structure and dynamics. To increase such resolution, the study of a particular target group is best approached with the use of specific primers. Here, we aimed to design and evaluate specific primers for aquatic picocyanobacterial genera to be used with high-throughput sequencing. Since the various regions of the 16S rRNA gene have different degrees of conservation in different bacterial groups, we therefore first determined which hypervariable region of the 16S rRNA gene provides the highest taxonomic and phylogenetic resolution for the genera Synechococcus, Prochlorococcus, and Cyanobium An in silico analysis showed that the V5, V6, and V7 hypervariable regions appear to be the most informative for this group. We then designed primers flanking these hypervariable regions and tested them in natural marine and freshwater communities. We successfully detected that most (97%) of the obtained reads could be assigned to picocyanobacterial genera. We defined operational taxonomic units as exact sequence variants (zero-radius operational taxonomic units [zOTUs]), which allowed us to detect higher genetic diversity and infer ecologically relevant information about picocyanobacterial community composition and dynamics in different aquatic systems. Our results open the door to future studies investigating picocyanobacterial diversity in aquatic systems.IMPORTANCE The molecular diversity of the aquatic picocyanobacterial community cannot be accurately described using only the available universal 16S rRNA gene primers that target the whole bacterial and archaeal community. We show that the hypervariable regions V5, V6, and V7 of the 16S rRNA gene are better suited to study the diversity, community structure, and dynamics of picocyanobacterial communities at a fine scale using Illumina MiSeq sequencing. Due to its variability, it allows reconstructing phylogenies featuring topologies comparable to those generated when using the complete 16S rRNA gene sequence. Further, we successfully designed a new set of primers flanking the V5 to V7 region whose specificity for picocyanobacterial genera was tested in silico and validated in several freshwater and marine aquatic communities. This work represents a step forward for understanding the diversity and ecology of aquatic picocyanobacteria and sets the path for future studies on picocyanobacterial diversity.

Latitudinal and Vertical Variation of Synechococcus Assemblage Composition Along 170° W Transect From the South Pacific to the Arctic Ocean

Synechococcus is one of the most widely distributed and abundant picocyanobacteria in the global oceans. Although latitudinal variation of Synechococcus assemblage in marine surface waters has been observed, few studies compared Synechococcus assemblage composition in surface and subsurface waters at the basin scale. Here, we report marine Synechococcus diversity in the surface and deep chlorophyll maximum (DCM) layers along 170° W from the South Pacific to the Arctic Ocean in summer. Along the transect, spatial niche partitioning of Synechococcus lineages in the surface waters was clearly observed. Species richness of surface Synechococcus assemblage was positively correlated with water temperature. Clade CRD1 was dominant in the areas (15° S-10° N and 35-40° N) associated with upwelling, and there were 3 different subclades with distinct distribution. CRD1-A was restricted in the North Equatorial Current (5-10° N), CRD1-B dominated in the equatorial upwelling region (15° S-0.17° N), and CRD1-C was only distributed in the North Pacific Current (35-40° N). Similarities between the Synechococcus assemblages in the surface and DCM layers were high at the upwelling regions and areas where the mixed layer was deep, while low in the Subtropical Gyres with strong stratification. Clade I, CRD1-B, and CRD1-C were major Synechococcus lineages in the DCM layer. In particular, clade I, which is composed of 7 subclades with distinct thermal niches, was widely distributed in the DCM layer. Overall, our results provide new insights into not only the latitudinal distribution of Synechococcus assemblages, but also their vertical variation in the central Pacific.

UCYN-A3, a newly characterized open ocean sublineage of the symbiotic N2 -fixing cyanobacterium Candidatus Atelocyanobacterium thalassa

The symbiotic unicellular cyanobacterium Candidatus Atelocyanobacterium thalassa (UCYN-A) is one of the most abundant and widespread nitrogen (N₂ )-fixing cyanobacteria in the ocean. Although it remains uncultivated, multiple sublineages have been detected based on partial nitrogenase (nifH) gene sequences, including the four most commonly detected sublineages UCYN-A1, UCYN-A2, UCYN-A3 and UCYN-A4. However, very little is known about UCYN-A3 beyond the nifH sequences from nifH gene diversity surveys. In this study, single cell sorting, DNA sequencing, qPCR and CARD-FISH assays revealed discrepancies involving the identification of sublineages, which led to new information on the diversity of the UCYN-A symbiosis. 16S rRNA and nifH gene sequencing on single sorted cells allowed us to identify the 16S rRNA gene of the uncharacterized UCYN-A3 sublineage. We designed new CARD-FISH probes that allowed us to distinguish and observe UCYN-A2 in a coastal location (SIO Pier; San Diego) and UCYN-A3 in an open ocean location (Station ALOHA; Hawaii). Moreover, we reconstructed about 13% of the UCYN-A3 genome from Tara Oceans metagenomic data. Finally, our findings unveil the UCYN-A3 symbiosis in open ocean waters suggesting that the different UCYN-A sublineages are distributed along different size fractions of the plankton defined by the cell-size ranges of their prymnesiophyte hosts.

Parallel phylogeography of Prochlorococcus and Synechococcus

The globally abundant marine Cyanobacteria Prochlorococcus and Synechococcus share many physiological traits but presumably have different evolutionary histories and associated phylogeography. In Prochlorococcus, there is a clear phylogenetic hierarchy of ecotypes, whereas multiple Synechococcus clades have overlapping physiologies and environmental distributions. However, microbial traits are associated with different phylogenetic depths. Using this principle, we reclassified diversity at different phylogenetic levels and compared the phylogeography. We sequenced the genetic diversity of Prochlorococcus and Synechococcus from 339 samples across the tropical Pacific Ocean and North Atlantic Ocean using a highly variable phylogenetic marker gene (rpoC1). We observed clear parallel niche distributions of ecotypes leading to high Pianka's Index values driven by distinct shifts at two transition points. The first transition point at 6°N distinguished ecotypes adapted to warm waters but separated by macronutrient content. At 39°N, ecotypes adapted to warm, low macronutrient vs. colder, high macronutrient waters shifted. Finally, we detected parallel vertical and regional single-nucleotide polymorphism microdiversity within clades from both Prochlorococcus and Synechococcus, suggesting uniquely adapted populations at very specific depths, as well as between the Atlantic and Pacific Oceans. Overall, this study demonstrates that Prochlorococcus and Synechococcus have shared phylogenetic organization of traits and associated phylogeography.

Temporal dynamics of eukaryotic microbial diversity at a coastal Pacific site

High-throughput sequencing of ocean biomes has revealed vast eukaryotic microbial diversity, a significant proportion of which remains uncharacterized. Here we use a temporal approach to understanding eukaryotic diversity at the Scripps Pier, La Jolla, California, USA, via high-throughput amplicon sequencing of the 18S rRNA gene, the abundances of both Synechococcus and Synechococcus grazers, and traditional oceanographic parameters. We also exploit our ability to track operational taxonomic units (OTUs) temporally to evaluate the ability of 18S sequence-based OTU assignments to meaningfully reflect ecological dynamics. The eukaryotic community is highly dynamic in terms of both species richness and composition, although proportional representation of higher-order taxa remains fairly consistent over time. Synechococcus abundance fluctuates throughout the year. OTUs unique to dates of Synechococcus blooms and crashes or enriched in Synechococcus addition incubation experiments suggest that the prasinophyte Tetraselmis sp. and Gymnodinium-like dinoflagellates are likely Synechococcus grazers under certain conditions, and may play an important role in their population fluctuations.

Seasonal variability of picophytoplankton under contrasting environments in northern Tunisian coasts, southwestern Mediterranean Sea

We investigated at the single cell level during 16months (June 2012 to September 2013) the temporal distribution of picophytoplankton (picoeukaryotes, Synechococcus and Prochlorococcus) communities in two contrasted ecosystems: the Bay of Bizerte characterised by an oligotrophic regime typical of the Mediterranean Sea and the Bizerte Lagoon that exhibits a mesotrophic/eutrophic state. We aimed at depicting seasonal variations and quantifying the relationships between the environmental factors and the structure and abundance of picophytoplankton communities. Results showed that picophytoplankton groups were able to grow under a wide range of environmental conditions varying seasonally, although their abundances and contributions to the total chlorophyll biomass significantly varied and showed importance in the Bay of Bizerte. Synechococcus was the most abundant group reaching 225∗10³cells·cm^-3 in the Bay and 278∗10³cells·cm^-3 in the lagoon. This group was present all over the year in both ecosystems. Structural equation model results pointed out a different configuration regarding the picophytoplankton environmental drivers. The complexity of the configuration, i.e. number of significant links within the system, decreased under enhanced eutrophication conditions. The less exposure to anthropogenic stress, i.e. in the Bay of Bizerte, highlight a larger role of nutrient and hydrological conditions on the seasonal variations of picophytoplankton, whereas a negative effect of eutrophication on picophytoplankton communities was unveiled in the Bizerte Lagoon. We stress that such influence may be exacerbated under expected scenarios of Mediterranean warming conditions and nutrient release in coastal ecosystems.

Photosynthetic functions of Synechococcus in the ocean microbiomes of diverse salinity and seasons

Synechococcus is an important photosynthetic picoplankton in the temperate to tropical oceans. As a photosynthetic bacterium, Synechococcus has an efficient mechanism to adapt to the changes in salinity and light intensity. The analysis of the distributions and functions of such microorganisms in the ever changing river mouth environment, where freshwater and seawater mix, should help better understand their roles in the ecosystem. Toward this objective, we have collected and sequenced the ocean microbiome in the river mouth of Kwangyang Bay, Korea, as a function of salinity and temperature. In conjunction with comparative genomics approaches using the sequenced genomes of a wide phylogeny of Synechococcus, the ocean microbiome was analyzed in terms of their composition and clade-specific functions. The results showed significant differences in the compositions of Synechococcus sampled in different seasons. The photosynthetic functions in such enhanced Synechococcus strains were also observed in the microbiomes in summer, which is significantly different from those in other seasons.

The microbiome associated with two Synechococcus ribotypes at different levels of ecological interaction

Planktonic cyanobacteria belonging to the genus Synechococcus are ubiquitously distributed in marine and fresh waters, substantially contributing to total carbon fixation on a global scale. While their ecological relevance is acknowledged, increasing resolution in molecular techniques allows disentangling cyanobacteria's role at the micro-scale, where complex microbial interactions may drive the overall community assembly. The interplay between phylogenetically different Synechococcus clades and their associated bacterial communities can affect their ecological fate and susceptibility to protistan predation. In this study, we experimentally promoted different levels of ecological interaction by mixing two Synechococcus ribotypes (MW101C3 and LL) and their associated bacteria, with and without a nanoflagellate grazer (Poterioochromonas sp.) in laboratory cultures. The beta-diversity of the Synechococcus-associated microbiome in laboratory cultures indicated that the presence of the LL ribotype was the main factor determining community composition changes (41% of total variance), and prevailed over the effect of protistan predation (18% of total variance). Our outcomes also showed that species coexistence and predation may promote microbial diversity, thus highlighting the underrated ecological relevance of such micro-scale factors.

Microbiome Dynamics in a Shrimp Grow-out Pond with Possible Outbreak of Acute Hepatopancreatic Necrosis Disease

Acute hepatopancreatic necrosis disease (AHPND) (formerly, early mortality syndrome) is a high-mortality-rate shrimp disease prevalent in shrimp farming areas. Although AHPND is known to be caused by pathogenic Vibrio parahaemolyticus hosting the plasmid-related PirAB^vp toxin gene, the effects of disturbances in microbiome have not yet been studied. We took 62 samples from a grow-out pond during an AHPND developing period from Days 23 to 37 after stocking white postlarvae shrimp and sequenced the 16S rRNA genes with Illumina sequencing technology. The microbiomes of pond seawater and shrimp stomachs underwent varied dynamic succession during the period. Despite copies of PirAB^vp, principal co-ordinates analysis revealed two distinctive stages of change in stomach microbiomes associated with AHPND. AHPND markedly changed the bacterial diversity in the stomachs; it decreased the Shannon index by 53.6% within approximately 7 days, shifted the microbiome with Vibrio and Candidatus Bacilloplasma as predominant populations, and altered the species-to-species connectivity and complexity of the interaction network. The AHPND-causing Vibrio species were predicted to develop a co-occurrence pattern with several resident and transit members within Candidatus Bacilloplasma and Cyanobacteria. This study’s insights into microbiome dynamics during AHPND infection can be valuable for minimising this disease in shrimp farming ponds.

Seasonal Succession and Spatial Patterns of Synechococcus Microdiversity in a Salt Marsh Estuary Revealed through 16S rRNA Gene Oligotyping

Synechococcus are ubiquitous and cosmopolitan cyanobacteria that play important roles in global productivity and biogeochemical cycles. This study investigated the fine scale microdiversity, seasonal patterns, and spatial distributions of Synechococcus in estuarine waters of Little Sippewissett salt marsh (LSM) on Cape Cod, MA. The proportion of Synechococcus reads was higher in the summer than winter, and higher in coastal waters than within the estuary. Variations in the V4-V6 region of the bacterial 16S rRNA gene revealed 12 unique Synechococcus oligotypes. Two distinct communities emerged in early and late summer, each comprising a different set of statistically co-occurring Synechococcus oligotypes from different clades. The early summer community included clades I and IV, which correlated with lower temperature and higher dissolved oxygen levels. The late summer community included clades CB5, I, IV, and VI, which correlated with higher temperatures and higher salinity levels. Four rare oligotypes occurred in the late summer community, and their relative abundances more strongly correlated with high salinity than did other co-occurring oligotypes. The analysis revealed that multiple, closely related oligotypes comprised certain abundant clades (e.g., clade 1 in the early summer and clade CB5 in the late summer), but the correlations between these oligotypes varied from pair to pair, suggesting they had slightly different niches despite being closely related at the clade level. Lack of tidal water exchange between sampling stations gave rise to a unique oligotype not abundant at other locations in the estuary, suggesting physical isolation plays a role in generating additional microdiversity within the community. Together, these results contribute to our understanding of the environmental and ecological factors that influence patterns of Synechococcus microbial community composition over space and time in salt marsh estuarine waters.

Unscrambling Cyanobacteria Community Dynamics Related to Environmental Factors

Future climate scenarios in the Baltic Sea project an increase of cyanobacterial bloom frequency and duration, attributed to eutrophication and climate change. Some cyanobacteria can be toxic and their impact on ecosystem services is relevant for a sustainable sea. Yet, there is limited understanding of the mechanisms regulating cyanobacterial diversity and biogeography. Here we unravel successional patterns and changes in cyanobacterial community structure using a 2-year monthly time- series during the productive season in a 100 km coastal-offshore transect using microscopy and high-throughput sequencing of 16S rRNA gene fragments. A total of 565 cyanobacterial OTUs were found, of which 231 where filamentous/colonial and 334 picocyanobacterial. Spatial differences in community structure between coastal and offshore waters were minor. An “epidemic population structure” (dominance of asingle cluster) was found for Aphanizomenon/Dolichospermum within the filamentous/colonial cyanobacterial community. In summer, this clusters imultaneously occurred with opportunistic clusters/OTUs, e.g., Nodularia spumigena and Pseudanabaena. Picocyanobacteria, Synechococcus/Cyanobium, formeda consistent but highly diverse group. Overall, the potential drivers structuring summer cyanobacterial communities were temperature and salinity. However, the different responses to environmental factors among and within genera suggest high niche specificity for individual OTUs. The recruitment and occurrence of potentially toxic filamentous/colonial clusters was likely related to disturbance such as mixing events and short-term shifts in salinity, and not solely dependent on increasing temperature and nitrogen-limiting conditions. Nutrients did not explain further the changes in cyanobacterial community composition. Novel occurrence patterns were identified as a strong seasonal succession revealing a tight coupling between the emergence of opportunistic picocynobacteria and the bloom offilamentous/colonialclusters. These findings highlight that if environmental conditions can partially explain the presence of opportunistic picocyanobacteria, microbial and trophic interactions with filamentous/colonial cyanobacteria should also be considered as potential shaping factors for single-celled communities. Regional climate change scenarios in the Baltic Sea predict environmental shifts leading to higher temperature and lower salinity; conditions identified here as favorable for opportunistic filamentous/colonial cyanobacteria. Altogether, the diversity and complexity of cyanobacterial communities reported here is far greater than previously known, emphasizing the importance of microbial interactions between filamentous and picocyanobacteria in the context of environmental disturbances.

Response of Spring Diatoms to CO2 Availability in the Western North Pacific as Determined by Next-Generation Sequencing

Next-generation sequencing (NGS) technologies have enabled us to determine phytoplankton community compositions at high resolution. However, few studies have adopted this approach to assess the responses of natural phytoplankton communities to environmental change. Here, we report the impact of different CO₂ levels on spring diatoms in the Oyashio region of the western North Pacific as estimated by NGS of the diatom-specific rbcL gene (DNA), which encodes the large subunit of RubisCO. We also examined the abundance and composition of rbcL transcripts (cDNA) in diatoms to assess their physiological responses to changing CO₂ levels. A short-term (3-day) incubation experiment was carried out on-deck using surface Oyashio waters under different pCO₂ levels (180, 350, 750, and 1000 μatm) in May 2011. During the incubation, the transcript abundance of the diatom-specific rbcL gene decreased with an increase in seawater pCO₂ levels. These results suggest that CO₂ fixation capacity of diatoms decreased rapidly under elevated CO₂ levels. In the high CO₂ treatments (750 and 1000 μatm), diversity of diatom-specific rbcL gene and its transcripts decreased relative to the control treatment (350 μatm), as well as contributions of Chaetocerataceae, Thalassiosiraceae, and Fragilariaceae to the total population, but the contributions of Bacillariaceae increased. In the low CO₂ treatment, contributions of Bacillariaceae also increased together with other eukaryotes. These suggest that changes in CO₂ levels can alter the community composition of spring diatoms in the Oyashio region. Overall, the NGS technology provided us a deeper understanding of the response of diatoms to changes in CO₂ levels in terms of their community composition, diversity, and photosynthetic physiology.

Diversity of Synechococcus at the Martha's Vineyard Coastal Observatory: Insights from Culture Isolations, Clone Libraries, and Flow Cytometry

The cyanobacterium Synechococcus is a ubiquitous, important phytoplankter across the world's oceans. A high degree of genetic diversity exists within the marine group, which likely contributes to its global success. Over 20 clades with different distribution patterns have been identified. However, we do not fully understand the environmental factors that control clade distributions. These factors are likely to change seasonally, especially in dynamic coastal systems. To investigate how coastal Synechococcus assemblages change temporally, we assessed the diversity of Synechococcus at the Martha's Vineyard Coastal Observatory (MVCO) over three annual cycles with culture-dependent and independent approaches. We further investigated the abundance of both phycoerythrin (PE)-containing and phycocyanin (PC)-only Synechococcus with a flow cytometric setup that distinguishes PC-only Synechococcus from picoeukaryotes. We found that the Synechococcus assemblage at MVCO is diverse (13 different clades identified), but dominated by clade I representatives. Many clades were only isolated during late summer and fall, suggesting more favorable conditions for isolation at this time. PC-only strains from four different clades were isolated, but these cells were only detected by flow cytometry in a few samples over the time series, suggesting they are rare at this site. Within clade I, we identified four distinct subclades. The relative abundances of each subclade varied over the seasonal cycle, and the high Synechococcus cell concentration at MVCO may be maintained by the diversity found within this clade. This study highlights the need to understand how temporal aspects of the environment affect Synechococcus community structure and cell abundance.

Comparative genomics of Synechococcus and proposal of the new genus Parasynechococcus

Synechococcus is among the most important contributors to global primary productivity. The genomes of several strains of this taxon have been previously sequenced in an effort to understand the physiology and ecology of these highly diverse microorganisms. Here we present a comparative study of Synechococcus genomes. For that end, we developed GenTaxo, a program written in Perl to perform genomic taxonomy based on average nucleotide identity, average amino acid identity and dinucleotide signatures, which revealed that the analyzed strains are drastically distinct regarding their genomic content. Phylogenomic reconstruction indicated a division of Synechococcus in two clades (i.e. Synechococcus and the new genus Parasynechococcus), corroborating evidences that this is in fact a polyphyletic group. By clustering protein encoding genes into homologue groups we were able to trace the Pangenome and core genome of both marine and freshwater Synechococcus and determine the genotypic traits that differentiate these lineages.

Comparison of the seasonal variations of Synechococcus assemblage structures in estuarine waters and coastal waters of Hong Kong

Seasonal variation in the phylogenetic composition of Synechococcus assemblages in estuarine and coastal waters of Hong Kong was examined through pyrosequencing of the rpoC1 gene. Sixteen samples were collected in 2009 from two stations representing estuarine and ocean-influenced coastal waters, respectively. Synechococcus abundance in coastal waters gradually increased from 3.6 × 10(3) cells ml(-1) in March, reaching a peak value of 5.7 × 10(5) cells ml(-1) in July, and then gradually decreased to 9.3 × 10(3) cells ml(-1) in December. The changes in Synechococcus abundance in estuarine waters followed a pattern similar to that in coastal waters, whereas its composition shifted from being dominated by phycoerythrin-rich (PE-type) strains in winter to phycocyanin-only (PC-type) strains in summer owing to the increase in freshwater discharge from the Pearl River and higher water temperature. The high abundance of PC-type Synechococcus was composed of subcluster 5.2 marine Synechococcus, freshwater Synechococcus (F-PC), and Cyanobium. The Synechococcus assemblage in the coastal waters, on the other hand, was dominated by marine PE-type Synechococcus, with subcluster 5.1 clades II and VI as the major lineages from April to September, when the summer monsoon prevailed. Besides these two clades, clade III cooccurred with clade V at relatively high abundance in summer. During winter, the Synechococcus assemblage compositions at the two sites were similar and were dominated by subcluster 5.1 clades II and IX and an undescribed clade (represented by Synechococcus sp. strain miyav). Clade IX Synechococcus was a relatively ubiquitous PE-type Synechococcus found at both sites, and our study demonstrates that some strains of the clade have the ability to deal with large variation of salinity in subtropical estuarine environments. Our study suggests that changes in seawater temperature and salinity caused by the seasonal variation of monsoonal forcing are two major determinants of the community composition and abundance of Synechococcus assemblages in Hong Kong waters.

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.