Evolutionary Mechanisms of Long-Term Genome Diversification Associated With Niche Partitioning in Marine Picocyanobacteria

Uncovering the genomic basis of symbiotic interactions and niche adaptations in freshwater picocyanobacteria

BACKGROUND

Picocyanobacteria from the genera Prochlorococcus, Synechococcus, and Cyanobium are the most widespread photosynthetic organisms in aquatic ecosystems. However, their freshwater populations remain poorly explored, due to uneven and insufficient sampling across diverse inland waterbodies.

RESULTS

In this study, we present 170 high-quality genomes of freshwater picocyanobacteria from non-axenic cultures collected across Central Europe. In addition, we recovered 33 genomes of their potential symbiotic partners affiliated with four genera, Pseudomonas, Mesorhizobium, Acidovorax, and Hydrogenophaga. The genomic basis of symbiotic interactions involved heterotrophs benefiting from picocyanobacteria-derived nutrients while providing detoxification of ROS. The global abundance patterns of picocyanobacteria revealed ecologically significant ecotypes, associated with trophic status, temperature, and pH as key environmental factors. The adaptation of picocyanobacteria in (hyper-)eutrophic waterbodies could be attributed to their colonial lifestyles and CRISPR-Cas systems. The prevailing CRISPR-Cas subtypes in picocyanobacteria were I-G and I-E, which appear to have been acquired through horizontal gene transfer from other bacterial phyla.

CONCLUSIONS

Our findings provide novel insights into the population diversity, ecology, and evolutionary strategies of the most widespread photoautotrophs within freshwater ecosystems. Video Abstract.

A distinct, high-affinity, alkaline phosphatase facilitates occupation of P-depleted environments by marine picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus, the two most abundant phototrophs on Earth, thrive in oligotrophic oceanic regions. While it is well known that specific lineages are exquisitely adapted to prevailing in situ light and temperature regimes, much less is known of the molecular machinery required to facilitate occupancy of these low-nutrient environments. Here, we describe a hitherto unknown alkaline phosphatase, Psip1, that has a substantially higher affinity for phosphomonoesters than other well-known phosphatases like PhoA, PhoX, or PhoD and is restricted to clade III Synechococcus and a subset of high light I-adapted Prochlorococcus strains, suggesting niche specificity. We demonstrate that Psip1 has undergone convergent evolution with PhoX, requiring both iron and calcium for activity and likely possessing identical key residues around the active site, despite generally very low sequence homology. Interrogation of metagenomes and transcriptomes from TARA oceans and an Atlantic Meridional transect shows that psip1 is abundant and highly expressed in picocyanobacterial populations from the Mediterranean Sea and north Atlantic gyre, regions well recognized to be phosphorus (P)-deplete. Together, this identifies psip1 as an important oligotrophy-specific gene for P recycling in these organisms. Furthermore, psip1 is not restricted to picocyanobacteria and is abundant and highly transcribed in some α-proteobacteria and eukaryotic algae, suggesting that such a high-affinity phosphatase is important across the microbial taxonomic world to occupy low-P environments.

Differential acclimation kinetics of the two forms of type IV chromatic acclimaters occurring in marine Synechococcus cyanobacteria

Synechococcus, the second most abundant marine phytoplanktonic organism, displays the widest variety of pigment content of all marine oxyphototrophs, explaining its ability to colonize all spectral niches occurring in the upper lit layer of oceans. Seven Synechococcus pigment types (PTs) have been described so far based on the phycobiliprotein composition and chromophorylation of their light-harvesting complexes, called phycobilisomes. The most elaborate and abundant PT (3d) in the open ocean consists of cells capable of type IV chromatic acclimation (CA4), i.e., to reversibly modify the ratio of the blue light-absorbing phycourobilin (PUB) to the green light-absorbing phycoerythrobilin (PEB) in phycobilisome rods to match the ambient light color. Two genetically distinct types of chromatic acclimaters, so-called PTs 3dA and 3dB, occur at similar global abundance in the ocean, but the precise physiological differences between these two types and the reasons for their complementary niche partitioning in the field remain obscure. Here, photoacclimation experiments in different mixes of blue and green light of representatives of these two PTs demonstrated that they differ by the ratio of blue-to-green light required to trigger the CA4 process. Furthermore, shift experiments between 100% blue and 100% green light, and vice-versa, revealed significant discrepancies between the acclimation pace of the two types of chromatic acclimaters. This study provides novel insights into the finely tuned adaptation mechanisms used by Synechococcus cells to colonize the whole underwater light field.

Uncovering the diversity and distribution of biosynthetic gene clusters of prochlorosins and other putative RiPPs in marine Synechococcus strains

IMPORTANCE

Genome mining studies have revealed the remarkable combinatorial diversity of ribosomally synthesized and post-translationally modified peptides (RiPPs) in marine bacteria, including prochlorosins. However, mining strategies also prove valuable in investigating the genomic landscape of associated genes within biosynthetic gene cluster (BGC) specific to targeted RiPPs of interest. Our study contributes to the enrichment of knowledge regarding prochlorosin diversity. It offers insights into potential mechanisms involved in their biosynthesis and modification, such as hyper-modification, which may give rise to active lantibiotics. Additionally, our study uncovers putative novel promiscuous post-translational enzymes, thereby expanding the chemical space explored within the Synechococcus genus. Moreover, this research extends the applications of mining techniques beyond the discovery of new RiPP-like clusters, allowing for a deeper understanding of genomics and diversity. Furthermore, it holds the potential to reveal previously unknown functions within the intriguing RiPP families, particularly in the case of prochlorosins.

Ecophysiological analysis reveals distinct environmental preferences in closely related Baltic Sea picocyanobacteria

Cluster 5 picocyanobacteria significantly contribute to primary productivity in aquatic ecosystems. Estuarine populations are highly diverse and consist of many co-occurring strains, but their physiology remains largely understudied. In this study, we characterized 17 novel estuarine picocyanobacterial strains. Phylogenetic analysis of the 16S rRNA and pigment genes (cpcB and cpeBA) uncovered multiple estuarine and freshwater-related clusters and pigment types. Assays with five representative strains (three phycocyanin rich and two phycoerythrin rich) under temperature (10-30°C), light (10-190 μmol photons m-2  s-1 ), and salinity (2-14 PSU) gradients revealed distinct growth optima and tolerance, indicating that genetic variability was accompanied by physiological diversity. Adaptability to environmental conditions was associated with differential pigment content and photosynthetic performance. Amplicon sequence variants at a coastal and an offshore station linked population dynamics with phylogenetic clusters, supporting that strains isolated in this study represent key ecotypes within the Baltic Sea picocyanobacterial community. The functional diversity found within strains with the same pigment type suggests that understanding estuarine picocyanobacterial ecology requires analysis beyond the phycocyanin and phycoerythrin divide. This new knowledge of the environmental preferences in estuarine picocyanobacteria is important for understanding and evaluating productivity in current and future ecosystems.

Coexistence of Dominant Marine Phytoplankton Sustained by Nutrient Specialization

Prochlorococcus and Synechococcus are the two dominant picocyanobacteria in the low-nutrient surface waters of the subtropical ocean, but the basis for their coexistence has not been quantitatively demonstrated. Here, we combine in situ microcosm experiments and an ecological model to show that this coexistence can be sustained by specialization in the uptake of distinct nitrogen (N) substrates at low-level concentrations that prevail in subtropical environments. In field incubations, the response of both Prochlorococcus and Synechococcus to nanomolar N amendments demonstrates N limitation of growth in both populations. However, Prochlorococcus showed a higher affinity to ammonium, whereas Synechococcus was more adapted to nitrate uptake. A simple ecological model demonstrates that the differential nutrient preference inferred from field experiments with these genera may sustain their coexistence. It also predicts that as the supply of NO3- decreases, as expected under climate warming, the dominant genera should undergo a nonlinear shift from Synechococcus to Prochlorococcus, a pattern that is supported by subtropical field observations. Our study suggests that the evolution of differential nutrient affinities is an important mechanism for sustaining the coexistence of genera and that climate change is likely to shift the relative abundance of the dominant plankton genera in the largest biomes in the ocean. IMPORTANCE Our manuscript addresses the following fundamental question in microbial ecology: how do different plankton using the same essential nutrients coexist? Prochlorococcus and Synechococcus are the two dominant picocyanobacteria in the low-nutrient surface waters of the subtropical ocean, which support a significant amount of marine primary production. The geographical distributions of these two organisms are largely overlapping, but the basis for their coexistence in these biomes remains unclear. In this study, we combined in situ microcosm experiments and an ecosystem model to show that the coexistence of these two organisms can arise from specialization in the uptake of distinct nitrogen substrates; Prochlorococcus prefers ammonium, whereas Synechococcus prefers nitrate when these nutrients exist at low concentrations. Our framework can be used for simulating and predicting the coexistence in the future ocean and may provide hints toward understanding other similar types of coexistence.

Marine picocyanobacterial PhnD1 shows specificity for various phosphorus sources but likely represents a constitutive inorganic phosphate transporter

Despite being fundamental to multiple biological processes, phosphorus (P) availability in marine environments is often growth-limiting, with generally low surface concentrations. Picocyanobacteria strains encode a putative ABC-type phosphite/phosphate/phosphonate transporter, phnDCE, thought to provide access to an alternative phosphorus pool. This, however, is paradoxical given most picocyanobacterial strains lack known phosphite degradation or carbon-phosphate lyase pathway to utilise alternate phosphorus pools. To understand the function of the PhnDCE transport system and its ecological consequences, we characterised the PhnD1 binding proteins from four distinct marine Synechococcus isolates (CC9311, CC9605, MITS9220, and WH8102). We show the Synechococcus PhnD1 proteins selectively bind phosphorus compounds with a stronger affinity for phosphite than for phosphate or methyl phosphonate. However, based on our comprehensive ligand screening and growth experiments showing Synechococcus strains WH8102 and MITS9220 cannot utilise phosphite or methylphosphonate as a sole phosphorus source, we hypothesise that the picocyanobacterial PhnDCE transporter is a constitutively expressed, medium-affinity phosphate transporter, and the measured affinity of PhnD1 to phosphite or methyl phosphonate is fortuitous. Our MITS9220_PhnD1 structure explains the comparatively lower affinity of picocyanobacterial PhnD1 for phosphate, resulting from a more limited H-bond network. We propose two possible physiological roles for PhnD1. First, it could function in phospholipid recycling, working together with the predicted phospholipase, TesA, and alkaline phosphatase. Second, by having multiple transporters for P (PhnDCE and Pst), picocyanobacteria could balance the need for rapid transport during transient episodes of higher P availability in the environment, with the need for efficient P utilisation in typical phosphate-deplete conditions.

Differential global distribution of marine picocyanobacteria gene clusters reveals distinct niche-related adaptive strategies

The ever-increasing number of available microbial genomes and metagenomes provides new opportunities to investigate the links between niche partitioning and genome evolution in the ocean, especially for the abundant and ubiquitous marine picocyanobacteria Prochlorococcus and Synechococcus. Here, by combining metagenome analyses of the Tara Oceans dataset with comparative genomics, including phyletic patterns and genomic context of individual genes from 256 reference genomes, we show that picocyanobacterial communities thriving in different niches possess distinct gene repertoires. We also identify clusters of adjacent genes that display specific distribution patterns in the field (eCAGs) and are thus potentially involved in the same metabolic pathway and may have a key role in niche adaptation. Several eCAGs are likely involved in the uptake or incorporation of complex organic forms of nutrients, such as guanidine, cyanate, cyanide, pyrimidine, or phosphonates, which might be either directly used by cells, for example for the biosynthesis of proteins or DNA, or degraded to inorganic nitrogen and/or phosphorus forms. We also highlight the enrichment of eCAGs involved in polysaccharide capsule biosynthesis in Synechococcus populations thriving in both nitrogen- and phosphorus-depleted areas vs. low-iron (Fe) regions, suggesting that the complexes they encode may be too energy-consuming for picocyanobacteria thriving in the latter areas. In contrast, Prochlorococcus populations thriving in Fe-depleted areas specifically possess an alternative respiratory terminal oxidase, potentially involved in the reduction of Fe(III) to Fe(II). Altogether, this study provides insights into how phytoplankton communities populate oceanic ecosystems, which is relevant to understanding their capacity to respond to ongoing climate change.

Genomic and Transcriptomic Insights into Salinity Tolerance-Based Niche Differentiation of Synechococcus Clades in Estuarine and Coastal Waters

Cluster 5 Synechococcus is one of the most important primary producers on earth. However, ecotypes of this genus exhibit complex geographical distributions, and the genetic basis of niche partitioning is still not fully understood. Here, we report distinct distributions of subcluster 5.1 (SC5.1) and subcluster 5.2 (SC5.2) Synechococcus in estuarine waters, and we reveal that salinity is the main factor determining their distribution. Clade III (belonging to SC5.1) and CB4 (belonging to SC5.2) are dominant clades in the study region, with different ecological distributions. We further conducted physiological, genomic, and transcriptomic studies of Synechococcus strains YX04-3 and HK05, which are affiliated with clade III and CB4, respectively. Laboratory tests showed that HK05 could grow at low salinity (13 ppt), whereas the growth of YX04-3 was suppressed when salinity decreased to 13 ppt. Genomic and transcriptomic analysis suggested that euryhaline clade CB4 is capable of dealing with a sudden drop of salinity by releasing compatible solutes through mechanosensitive channels that are coded by the mscL gene, decreasing biosynthesis of organic osmolytes, and increasing expression of heat shock proteins and high light-inducible proteins to protect photosystem. Furthermore, CB4 strain HK05 exhibited a higher growth rate when growing at low salinity than at high salinity. This is likely achieved by reducing its biosynthesis of organic osmolyte activity and increasing its photosynthetic activity at low salinity, which allowed it to enhance the assimilation of inorganic carbon and nitrogen. Together, these results provide new insights regarding the ecological distribution of SC5.2 and SC5.1 ecotypes and their underlying molecular mechanisms. IMPORTANCE Synechococcus is a group of unicellular Cyanobacteria that are widely distributed in global aquatic ecosystems. Salinity is a factor that affects the distribution of microorganisms in estuarine and coastal environments. In this study, we studied the distribution pattern of Synechococcus community along the salinity gradient in a subtropical estuary. By using omic methods, we unveiled genetic traits that determine the niche partitioning of euryhaline and strictly marine Synechococcus. We also explored the strategies employed by euryhaline Synechococcus to cope with a sudden drop of salinity, and revealed possible mechanisms for the higher growth rate of euryhaline Synechococcus in low salinity conditions. This study provides new insight into the genetic basis of niche partitioning of Synechococcus clades.

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.

ProX from marine Synechococcus spp. show a sole preference for glycine-betaine with differential affinity between ecotypes

Osmotic stress, caused by high or fluctuating salt concentrations, is a crucial abiotic factor affecting microbial growth in aquatic habitats. Many organisms utilize common responses to osmotic stress, generally requiring active extrusion of toxic inorganic ions and accumulation of compatible solutes to protect cellular machinery. We heterologously expressed and purified predicted osmoprotectant, proline/glycine betaine-binding proteins (ProX) from two phylogenetically distinct Synechococcus spp. MITS9220 and WH8102. Homologues of this protein are conserved only among Prochlorococcus LLIV and Synechococcus clade I, III and CRD1 strains. Our biophysical characterization show Synechococcus ProX exists as a dimer, with specificity solely for glycine betaine but not to other osmoprotectants tested. We discovered that MITS9220_ProX has a 10-fold higher affinity to glycine betaine than WH8102_ProX, which is further elevated (24-fold) in high salt conditions. The stronger affinity and effect of ionic strength on MITS9220_ProX glycine betaine binding but not on WH8102_ProX alludes to a novel regulatory mechanism, providing critical functional insights into the phylogenetic divergence of picocyanobacterial ProX proteins that may be necessary for their ecological success.

α-cyanobacteria possessing form IA RuBisCO globally dominate aquatic habitats

RuBisCO (ribulose 1,5-bisphosphate carboxylase/oxygenase) is one the most abundant enzymes on Earth. Virtually all food webs depend on its activity to supply fixed carbon. In aerobic environments, RuBisCO struggles to distinguish efficiently between CO2 and O2. To compensate, organisms have evolved convergent solutions to concentrate CO2 around the active site. The genetic engineering of such inorganic carbon concentrating mechanisms (CCMs) into plants could help facilitate future global food security for humankind. In bacteria, the carboxysome represents one such CCM component, of which two independent forms exist: α and β. Cyanobacteria are important players in the planet's carbon cycle and the vast majority of the phylum possess a β-carboxysome, including most cyanobacteria used as laboratory models. The exceptions are the exclusively marine Prochlorococcus and Synechococcus that numerically dominate open ocean systems. However, the reason why marine systems favor an α-form is currently unknown. Here, we report the genomes of 58 cyanobacteria, closely related to marine Synechococcus that were isolated from freshwater lakes across the globe. We find all these isolates possess α-carboxysomes accompanied by a form 1A RuBisCO. Moreover, we demonstrate α-cyanobacteria dominate freshwater lakes worldwide. Hence, the paradigm of a separation in carboxysome type across the salinity divide does not hold true, and instead the α-form dominates all aquatic systems. We thus question the relevance of β-cyanobacteria as models for aquatic systems at large and pose a hypothesis for the reason for the success of the α-form in nature.

Elucidating the picocyanobacteria salinity divide through ecogenomics of new freshwater isolates

BACKGROUND

Cyanobacteria are the major prokaryotic primary producers occupying a range of aquatic habitats worldwide that differ in levels of salinity, making them a group of interest to study one of the major unresolved conundrums in aquatic microbiology which is what distinguishes a marine microbe from a freshwater one? We address this question using ecogenomics of a group of picocyanobacteria (cluster 5) that have recently evolved to inhabit geographically disparate salinity niches. Our analysis is made possible by the sequencing of 58 new genomes from freshwater representatives of this group that are presented here, representing a 6-fold increase in the available genomic data.

RESULTS

Overall, freshwater strains had larger genomes (≈2.9 Mb) and %GC content (≈64%) compared to brackish (2.69 Mb and 64%) and marine (2.5 Mb and 58.5%) isolates. Genomic novelties/differences across the salinity divide highlighted acidic proteomes and specific salt adaptation pathways in marine isolates (e.g., osmolytes/compatible solutes - glycine betaine/ggp/gpg/gmg clusters and glycerolipids glpK/glpA), while freshwater strains possessed distinct ion/potassium channels, permeases (aquaporin Z), fatty acid desaturases, and more neutral/basic proteomes. Sulfur, nitrogen, phosphorus, carbon (photosynthesis), or stress tolerance metabolism while showing distinct genomic footprints between habitats, e.g., different types of transporters, did not obviously translate into major functionality differences between environments. Brackish microbes show a mixture of marine (salt adaptation pathways) and freshwater features, highlighting their transitional nature.

CONCLUSIONS

The plethora of freshwater isolates provided here, in terms of trophic status preference and genetic diversity, exemplifies their ability to colonize ecologically diverse waters across the globe. Moreover, a trend towards larger and more flexible/adaptive genomes in freshwater picocyanobacteria may hint at a wider number of ecological niches in this environment compared to the relatively homogeneous marine system.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the “pink berries” of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

Multiple Photolyases Protect the Marine Cyanobacterium Synechococcus from Ultraviolet Radiation

Marine cyanobacteria depend on light for photosynthesis, restricting their growth to the photic zone. The upper part of this layer is exposed to strong UV radiation (UVR), a DNA mutagen that can harm these microorganisms. To thrive in UVR-rich waters, marine cyanobacteria employ photoprotection strategies that are still not well defined. Among these are photolyases, light-activated enzymes that repair DNA dimers generated by UVR. Our analysis of genomes of 81 strains of Synechococcus, Cyanobium, and Prochlorococcus isolated from the world’s oceans shows that they possess up to five genes encoding different members of the photolyase/cryptochrome family, including a photolyase with a novel domain arrangement encoded by either one or two separate genes. We disrupted the putative photolyase-encoding genes in Synechococcus sp. strain RS9916 and discovered that each gene contributes to the overall capacity of this organism to survive UVR. Additionally, each conferred increased survival after UVR exposure when transformed into Escherichia coli lacking its photolyase and SOS response. Our results provide the first evidence that this large set of photolyases endows Synechococcus with UVR resistance that is far superior to that of E. coli, but that, unlike for E. coli, these photolyases provide Synechococcus with the vast majority of its UVR tolerance.

A single sensor controls large variations in zinc quotas in a marine cyanobacterium

Marine cyanobacteria are critical players in global nutrient cycles that crucially depend on trace metals in metalloenzymes, including zinc for CO₂ fixation and phosphorus acquisition. How strains proliferating in the vast oligotrophic ocean gyres thrive at ultra-low zinc concentrations is currently unknown. Using Synechococcus sp. WH8102 as a model we show that its zinc-sensor protein Zur differs from all other known bacterial Zur proteins in overall structure and the location of its sensory zinc site. Uniquely, Synechococcus Zur activates metallothionein gene expression, which supports cellular zinc quotas spanning two orders of magnitude. Thus, a single zinc sensor facilitates growth across pico- to micromolar zinc concentrations with the bonus of banking this precious resource. The resultant ability to grow well at both ultra-low and excess zinc, together with overall lower zinc requirements, likely contribute to the broad ecological distribution of Synechococcus across the global oceans.

The zinc-sensor protein Zur in a marine cyanobacterium is distinct from those in other bacteria in structure and location of its sensory zinc site, and facilitates growth across a range of zinc concentrations via activation of a metallothionein gene.

Comparative Thermophysiology of Marine Synechococcus CRD1 Strains Isolated From Different Thermal Niches in Iron-Depleted Areas

Marine Synechococcus cyanobacteria are ubiquitous in the ocean, a feature likely related to their extensive genetic diversity. Amongst the major lineages, clades I and IV preferentially thrive in temperate and cold, nutrient-rich waters, whilst clades II and III prefer warm, nitrogen or phosphorus-depleted waters. The existence of such cold (I/IV) and warm (II/III) thermotypes is corroborated by physiological characterization of representative strains. A fifth clade, CRD1, was recently shown to dominate the Synechococcus community in iron-depleted areas of the world ocean and to encompass three distinct ecologically significant taxonomic units (ESTUs CRD1A-C) occupying different thermal niches, suggesting that distinct thermotypes could also occur within this clade. Here, using comparative thermophysiology of strains representative of these three CRD1 ESTUs we show that the CRD1A strain MITS9220 is a warm thermotype, the CRD1B strain BIOS-U3-1 a cold temperate thermotype, and the CRD1C strain BIOS-E4-1 a warm temperate stenotherm. Curiously, the CRD1B thermotype lacks traits and/or genomic features typical of cold thermotypes. In contrast, we found specific physiological traits of the CRD1 strains compared to their clade I, II, III, and IV counterparts, including a lower growth rate and photosystem II maximal quantum yield at most temperatures and a higher turnover rate of the D1 protein. Together, our data suggests that the CRD1 clade prioritizes adaptation to low-iron conditions over temperature adaptation, even though the occurrence of several CRD1 thermotypes likely explains why the CRD1 clade as a whole occupies most iron-limited waters.

Novel functional insights into a modified sugar-binding protein from Synechococcus MITS9220

Paradigms of metabolic strategies employed by photoautotrophic marine picocyanobacteria have been challenged in recent years. Based on genomic annotations, picocyanobacteria are predicted to assimilate organic nutrients via ATP-binding cassette importers, a process mediated by substrate-binding proteins. We report the functional characterisation of a modified sugar-binding protein, MsBP, from a marine Synechococcus strain, MITS9220. Ligand screening of MsBP shows a specific affinity for zinc (K_D ~ 1.3 μM) and a preference for phosphate-modified sugars, such as fructose-1,6-biphosphate, in the presence of zinc (K_D ~ 5.8 μM). Our crystal structures of apo MsBP (no zinc or substrate-bound) and Zn-MsBP (with zinc-bound) show that the presence of zinc induces structural differences, leading to a partially-closed substrate-binding cavity. The Zn-MsBP structure also sequesters several sulphate ions from the crystallisation condition, including two in the binding cleft, appropriately placed to mimic the orientation of adducts of a biphosphate hexose. Combined with a previously unseen positively charged binding cleft in our two structures and our binding affinity data, these observations highlight novel molecular variations on the sugar-binding SBP scaffold. Our findings lend further evidence to a proposed sugar acquisition mechanism in picocyanobacteria alluding to a mixotrophic strategy within these ubiquitous photosynthetic bacteria.

Diverse Subclade Differentiation Attributed to the Ubiquity of Prochlorococcus High-Light-Adapted Clade II

Prochlorococcus is the key primary producer in marine ecosystems, and the high-light-adapted clade II (HLII) is the most abundant ecotype. However, the genomic and ecological basis of Prochlorococcus HLII in the marine environment has remained elusive. Here, we show that the ecologically coherent subclade differentiation of HLII corresponds to genomic and ecological characteristics on the basis of analyses of 31 different strains of HLII, including 12 novel isolates. Different subclades of HLII with different core and accessory genes were identified, and their distribution in the marine environment was explored using the TARA Oceans metagenome database. Three major subclade groups were identified, viz., the surface group (HLII-SG), the transition group (HLII-TG), and the deep group (HLII-DG). These subclade groups showed different temperature ranges and optima for distribution. In regression analyses, temperature and nutrient availability were identified as key factors affecting the distribution of HLII subclades. A 35% increase in the relative abundance of HLII-SG by the end of the 21st century was predicted under the Representative Concentration Pathway 8.5 scenario. Our results show that the ubiquity and distribution of Prochlorococcus HLII in the marine environment are associated with the differentiation of diverse subclades. These findings provide insights into the large-scale shifts in the Prochlorococcus community in response to future climate change.

Diversity and evolution of pigment types in marine Synechococcus cyanobacteria

Synechococcus cyanobacteria are ubiquitous and abundant in the marine environment and contribute to an estimated 16% of the ocean net primary productivity. Their light-harvesting complexes, called phycobilisomes (PBS), are composed of a conserved allophycocyanin core, from which radiates six to eight rods with variable phycobiliprotein and chromophore content. This variability allows Synechococcus cells to optimally exploit the wide variety of spectral niches existing in marine ecosystems. Seven distinct pigment types or subtypes have been identified so far in this taxon based on the phycobiliprotein composition and/or the proportion of the different chromophores in PBS rods. Most genes involved in their biosynthesis and regulation are located in a dedicated genomic region called the PBS rod region. Here, we examine the variability of gene content and organization of this genomic region in a large set of sequenced isolates and natural populations of Synechococcus representative of all known pigment types. All regions start with a tRNA-Phe_GAA and some possess mobile elements for DNA integration and site-specific recombination, suggesting that their genomic variability relies in part on a “tycheposon”-like mechanism. Comparison of the phylogenies obtained for PBS and core genes revealed that the evolutionary history of PBS rod genes differs from the core genome and is characterized by the co-existence of different alleles and frequent allelic exchange. We propose a scenario for the evolution of the different pigment types and highlight the importance of incomplete lineage sorting in maintaining a wide diversity of pigment types in different Synechococcus lineages despite multiple speciation events.

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.

Functional characterisation of substrate-binding proteins to address nutrient uptake in marine picocyanobacteria

Marine cyanobacteria are key primary producers, contributing significantly to the microbial food web and biogeochemical cycles by releasing and importing many essential nutrients cycled through the environment. A subgroup of these, the picocyanobacteria (Synechococcus and Prochlorococcus), have colonised almost all marine ecosystems, covering a range of distinct light and temperature conditions, and nutrient profiles. The intra-clade diversities displayed by this monophyletic branch of cyanobacteria is indicative of their success across a broad range of environments. Part of this diversity is due to nutrient acquisition mechanisms, such as the use of high-affinity ATP-binding cassette (ABC) transporters to competitively acquire nutrients, particularly in oligotrophic (nutrient scarce) marine environments. The specificity of nutrient uptake in ABC transporters is primarily determined by the peripheral substrate-binding protein (SBP), a receptor protein that mediates ligand recognition and initiates translocation into the cell. The recent availability of large numbers of sequenced picocyanobacterial genomes indicates both Synechococcus and Prochlorococcus apportion >50% of their transport capacity to ABC transport systems. However, the low degree of sequence homology among the SBP family limits the reliability of functional assignments using sequence annotation and prediction tools. This review highlights the use of known SBP structural representatives for the uptake of key nutrient classes by cyanobacteria to compare with predicted SBP functionalities within sequenced marine picocyanobacteria genomes. This review shows the broad range of conserved biochemical functions of picocyanobacteria and the range of novel and hypothetical ABC transport systems that require further functional characterisation.

Marine Synechococcus picocyanobacteria: Light utilization across latitudes

The most ubiquitous cyanobacteria, Synechococcus, have colonized different marine thermal niches through the evolutionary specialization of lineages adapted to different ranges of temperature seawater. We used the strains of Synechococcus temperature ecotypes to study how light utilization has evolved in the function of temperature. The tropical Synechococcus (clade II) was unable to grow under 16 °C but, at temperatures >25 °C, induced very high growth rates that relied on a strong synthesis of the components of the photosynthetic machinery, leading to a large increase in photosystem cross-section and electron flux. By contrast, the Synechococcus adapted to subpolar habitats (clade I) grew more slowly but was able to cope with temperatures <10 °C. We show that growth at such temperatures was accompanied by a large increase of the photoprotection capacities using the orange carotenoid protein (OCP). Metagenomic analyzes revealed that Synechococcus natural communities show the highest prevalence of the ocp genes in low-temperature niches, whereas most tropical clade II Synechococcus have lost the gene. Moreover, bioinformatic analyzes suggested that the OCP variants of the two cold-adapted Synechococcus clades I and IV have undergone evolutionary convergence through the adaptation of the molecular flexibility. Our study points to an important role of temperature in the evolution of the OCP. We, furthermore, discuss the implications of the different metabolic cost of these physiological strategies on the competitiveness of Synechococcus in a warming ocean. This study can help improve the current hypotheses and models aimed at predicting the changes in ocean carbon fluxes in response to global warming.

Evolution of kaiA, a key circadian gene of cyanobacteria

The circadian system of cyanobacteria is built upon a central oscillator consisting of three genes, kaiA, kaiB, and kaiC. The KaiA protein plays a key role in phosphorylation/dephosphorylation cycles of KaiC, which occur over the 24-h period. We conducted a comprehensive evolutionary analysis of the kaiA genes across cyanobacteria. The results show that, in contrast to the previous reports, kaiA has an ancient origin and is as old as cyanobacteria. The kaiA homologs are present in nearly all analyzed cyanobacteria, except Gloeobacter, and have varying domain architecture. Some Prochlorococcales, which were previously reported to lack the kaiA gene, possess a drastically truncated homolog. The existence of the diverse kaiA homologs suggests significant variation of the circadian mechanism, which was described for the model cyanobacterium, Synechococcus elongatus PCC7942. The major structural modifications in the kaiA genes (duplications, acquisition and loss of domains) have apparently been induced by global environmental changes in the different geological periods.

Temperature Responses of Heterotrophic Bacteria in Co-culture With a Red Sea Synechococcus Strain

Interactions between autotrophic and heterotrophic bacteria are fundamental for marine biogeochemical cycling. How global warming will affect the dynamics of these essential microbial players is not fully understood. The aims of this study were to identify the major groups of heterotrophic bacteria present in a Synechococcus culture originally isolated from the Red Sea and assess their joint responses to experimental warming within the metabolic ecology framework. A co-culture of Synechococcus sp. RS9907 and their associated heterotrophic bacteria, after determining their taxonomic affiliation by 16S rRNA gene sequencing, was acclimated and maintained in the lab at different temperatures (24-34°C). The abundance and cellular properties of Synechococcus and the three dominant heterotrophic bacterial groups (pertaining to the genera Paracoccus, Marinobacter, and Muricauda) were monitored by flow cytometry. The activation energy of Synechococcus, which grew at 0.94-1.38 d-1, was very similar (0.34 ± 0.02 eV) to the value hypothesized by the metabolic theory of ecology (MTE) for autotrophs (0.32 eV), while the values of the three heterotrophic bacteria ranged from 0.16 to 1.15 eV and were negatively correlated with their corresponding specific growth rates (2.38-24.4 d-1). The corresponding carrying capacities did not always follow the inverse relationship with temperature predicted by MTE, nor did we observe a consistent response of bacterial cell size and temperature. Our results show that the responses to future ocean warming of autotrophic and heterotrophic bacteria in microbial consortia might not be well described by theoretical universal rules.

Cyanorak v2.1: a scalable information system dedicated to the visualization and expert curation of marine and brackish picocyanobacteria genomes

Cyanorak v2.1 (http://www.sb-roscoff.fr/cyanorak) is an information system dedicated to visualizing, comparing and curating the genomes of Prochlorococcus, Synechococcus and Cyanobium, the most abundant photosynthetic microorganisms on Earth. The database encompasses sequences from 97 genomes, covering most of the wide genetic diversity known so far within these groups, and which were split into 25,834 clusters of likely orthologous groups (CLOGs). The user interface gives access to genomic characteristics, accession numbers as well as an interactive map showing strain isolation sites. The main entry to the database is through search for a term (gene name, product, etc.), resulting in a list of CLOGs and individual genes. Each CLOG benefits from a rich functional annotation including EggNOG, EC/K numbers, GO terms, TIGR Roles, custom-designed Cyanorak Roles as well as several protein motif predictions. Cyanorak also displays a phyletic profile, indicating the genotype and pigment type for each CLOG, and a genome viewer (Jbrowse) to visualize additional data on each genome such as predicted operons, genomic islands or transcriptomic data, when available. This information system also includes a BLAST search tool, comparative genomic context as well as various data export options. Altogether, Cyanorak v2.1 constitutes an invaluable, scalable tool for comparative genomics of ecologically relevant marine microorganisms.