Parallel phylogeography of Prochlorococcus and Synechococcus

Ribosomal protein phylogeography offers quantitative insights into the efficacy of genome-resolved surveys of microbial communities

The increasing availability of microbial genomes is essential to gain insights into microbial ecology and evolution that can propel biotechnological and biomedical advances. Recent advances in genome recovery have significantly expanded the catalogue of microbial genomes from diverse habitats. However, the ability to explain how well a set of genomes account for the diversity in a given environment remains challenging for individual studies or biome-specific databases. Here we present EcoPhylo, a computational workflow to characterize the phylogeography of any gene family through integrated analyses of genomes and metagenomes, and our application of this approach to ribosomal proteins to quantify phylogeny-aware genome recovery rates across three biomes. Our findings show that genome recovery rates vary widely across taxa and biomes, and that single amplified genomes, metagenome-assembled genomes, and isolate genomes have non-uniform yet quantifiable representation of environmental microbes. EcoPhylo reveals highly resolved, reference-free, multi-domain phylogenies in conjunction with distribution patterns of individual clades across environments, providing a means to assess genome recovery in individual studies and benchmark biome-level genome collections.

Prochlorococcus marinus responses to light and oxygen

Prochlorococcus marinus, the smallest picocyanobacterium, comprises multiple clades occupying distinct niches, currently across tropical and sub-tropical oligotrophic ocean regions, including Oxygen Minimum Zones. Ocean warming may open growth-permissive temperatures in new, poleward photic regimes, along with expanded Oxygen Minimum Zones. We used ocean metaproteomic data on current Prochlorococcus marinus niches, to guide testing of Prochlorococcus marinus growth across a matrix of peak irradiances, photoperiods, spectral bands and dissolved oxygen. MED4 from Clade HLI requires greater than 4 h photoperiod, grows at 25 μmol O2 L-1 and above, and exploits high cumulative diel photon doses. MED4, however, relies upon an alternative oxidase to balance electron transport, which may exclude it from growth under our lowest, 2.5 μmol O2 L-1, condition. SS120 from clade LLII/III is restricted to low light under full 250 μmol O2 L-1, shows expanded light exploitation under 25 μmol O2 L-1, but is excluded from growth under 2.5 μmol O2 L-1. Intermediate oxygen suppresses the cost of PSII photoinactivation, and possibly the enzymatic production of H2O2 in SS120, which has limitations on genomic capacity for PSII and DNA repair. MIT9313 from Clade LLIV is restricted to low blue irradiance under 250 μmol O2 L-1, but exploits much higher irradiance under red light, or under lower O2 concentrations, conditions which slow photoinactivation of PSII and production of reactive oxygen species. In warming oceans, range expansions and competition among clades will be governed not only by light levels. Short photoperiods governed by latitude, temperate winters, and depth attenuation of light, will exclude clade HLI (including MED4) from some habitats. In contrast, clade LLII/III (including SS120), and particularly clade LLIV (including MIT9313), may exploit higher light niches nearer the surface, under expanding OMZ conditions, where low O2 relieves the stresses of oxidation stress and PSII photoinhibition.

Interactions between the nitrogen-fixing cyanobacterium Trichodesmium and siderophore-producing cyanobacterium Synechococcus under iron limitation

As diazotrophic cyanobacteria of tremendous biomass, Trichodesmium continuously provide a nitrogen source for carbon-fixing cyanobacteria and drive the generation of primary productivity in marine environments. However, ocean iron deficiencies limit growth and metabolism of Trichodesmium. Recent studies have shown the co-occurrence of Trichodesmium and siderophore-producing Synechococcus in iron-deficient oceans, but whether siderophores secreted by Synechococcus can be used by Trichodesmium to adapt to iron deficiency is not clear. We constructed a mutant Synechococcus strain unable to produce siderophores to explore this issue. Synechococcus filtrates with or without siderophores were added into a Trichodesmium microbial consortium consisting of Trichodesmium erythraeum IMS 101 as the dominant microbe with chronic iron deficiency. By analyzing the physiological phenotype, metagenome, and metatranscriptome, we investigated the interactions between the nitrogen-fixing cyanobacterium Tricodesmium and siderophore-producing cyanobacterium Synechococcus under conditions of iron deficiency. The results indicated that siderophores secreted by Synechococcus are likely to chelate with free iron in the culture medium of the Trichodesmium consortium, reducing the concentration of bioavailable iron and posing greater challenges to the absorption of iron by Trichodesmium. These findings revealed the characteristics of iron-competitive utilization between diazotrophic cyanobacteria and siderophore-producing cyanobacteria, as well as potential interactions, and provide a scientific basis for understanding the regulatory effects of nutrient limitation on marine primary productivity.

Design and use of a new primer pair for the characterization of the cyanobacteria Synechococcus and Prochlorococcus communities targeting petB gene through metabarcoding approaches

During the last years, the application of next-generation sequencing (NGS) technologies to search for specific genetic markers has become a crucial method for the characterization of microbial communities. Illumina MiSeq, likely the most widespread NGS platform for metabarcoding experiments and taxonomic classification, allows processing shorter reads than the classical SANGER sequencing method and therefore requires specific primer pairs that produce shorter amplicons. Specifically, for the analysis of the commonly studied Prochlorococcus and Synechococcus communities, the petB marker gene has recently stood out as able to provide deep coverage to determine the microdiversity of the community. However, current petB primer set produce a 597 bp amplicon that is not suitable for MiSeq chemistry. Here, we designed and tested a petB primer pair that targets both Prochlorococcus and Synechococcus communities producing an appropriate amplicon to be used with state-of-the-art Illumina MiSeq. This new primer set allows the classification of both groups to a low taxonomic level and is therefore suitable for high throughput experiments using MiSeq technologies, therefore constituting a useful, novel tool to facilitate further studies on Prochlorococcus and Synechococcus communities. •This work describes the de novo design of a Prochlorococcus and Synechococcus-specific petB primer pair, allowing the characterization of both populations to a low taxonomic level.•This primer pair is suitable for widespread Illumina MiSeq sequencing technologies.•petB was confirmed as an adequate target for the characterization of both picocyanobacteria.

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.

The Bay of Bengal exposes abundant photosynthetic picoplankton and newfound diversity along salinity-driven gradients

The Bay of Bengal (BoB) is a 2,600,000 km2 expanse in the Indian Ocean upon which many humans rely. However, the primary producers underpinning food chains here remain poorly characterized. We examined phytoplankton abundance and diversity along strong BoB latitudinal and vertical salinity gradients-which have low temperature variation (27-29°C) between the surface and subsurface chlorophyll maximum (SCM). In surface waters, Prochlorococcus averaged 11.7 ± 4.4 × 104  cells ml-1 , predominantly HLII, whereas LLII and 'rare' ecotypes, HLVI and LLVII, dominated in the SCM. Synechococcus averaged 8.4 ± 2.3 × 104  cells ml-1 in the surface, declined rapidly with depth, and population structure of dominant Clade II differed between surface and SCM; Clade X was notable at both depths. Across all sites, Ostreococcus Clade OII dominated SCM eukaryotes whereas communities differentiated strongly moving from Arabian Sea-influenced high salinity (southerly; prasinophytes) to freshwater-influenced low salinity (northerly; stramenopiles, specifically, diatoms, pelagophytes, and dictyochophytes, plus the prasinophyte Micromonas) surface waters. Eukaryotic phytoplankton peaked in the south (1.9 × 104  cells ml-1 , surface) where a novel Ostreococcus was revealed, named here Ostreococcus bengalensis. We expose dominance of a single picoeukaryote and hitherto 'rare' picocyanobacteria at depth in this complex ecosystem where studies suggest picoplankton are replacing larger phytoplankton due to climate change.

Ecological mechanisms and current systems shape the modular structure of the global oceans' prokaryotic seascape

Major biogeographic features of the microbial seascape in the oceans have been established and their underlying ecological mechanisms in the (sub)tropical oceans and the Pacific Ocean identified. However, we still lack a unifying understanding of how prokaryotic communities and biogeographic patterns are affected by large-scale current systems in distinct ocean basins and how they are globally shaped in line with ecological mechanisms. Here we show that prokaryotic communities in the epipelagic Pacific and Atlantic Ocean, in the southern Indian Ocean, and the Mediterranean Sea are composed of modules of co-occurring taxa with similar environmental preferences. The relative partitioning of these modules varies along latitudinal and longitudinal gradients and are related to different hydrographic and biotic conditions. Homogeneous selection and dispersal limitation were identified as the major ecological mechanisms shaping these communities and their free-living (FL) and particle-associated (PA) fractions. Large-scale current systems govern the dispersal of prokaryotic modules leading to the highest diversity near subtropical fronts.

Global scale phylogeography of functional traits and microdiversity in Prochlorococcus

Prochlorococcus is the most numerically abundant photosynthetic organism in the surface ocean. The Prochlorococcus high-light and warm-water adapted ecotype (HLII) is comprised of extensive microdiversity, but specific functional differences between microdiverse sub-clades remain elusive. Here we characterized both functional and phylogenetic diversity within the HLII ecotype using Bio-GO-SHIP metagenomes. We found widespread variation in gene frequency connected to local environmental conditions. Metagenome-assembled marker genes and genomes revealed a globally distributed novel HLII haplotype defined by adaptation to chronically low P conditions (HLII-P). Environmental correlation analysis revealed different factors were driving gene abundances verses phylogenetic differences. An analysis of cultured HLII genomes and metagenome-assembled genomes revealed a subclade within HLII, which corresponded to the novel HLII-P haplotype. This work represents the first global assessment of the HLII ecotype's phylogeography and corresponding functional differences. These findings together expand our understanding of how microdiversity structures functional differences and reveals the importance of nutrients as drivers of microdiversity in Prochlorococcus.

Investigating the eco-evolutionary response of microbiomes to environmental change

Microorganisms are the primary engines of biogeochemical processes and foundational to the provisioning of ecosystem services to human society. Free-living microbial communities (microbiomes) and their functioning are now known to be highly sensitive to environmental change. Given microorganisms' capacity for rapid evolution, evolutionary processes could play a role in this response. Currently, however, few models of biogeochemical processes explicitly consider how microbial evolution will affect biogeochemical responses to environmental change. Here, we propose a conceptual framework for explicitly integrating evolution into microbiome-functioning relationships. We consider how microbiomes respond simultaneously to environmental change via four interrelated processes that affect overall microbiome functioning (physiological acclimation, demography, dispersal and evolution). Recent evidence in both the laboratory and the field suggests that ecological and evolutionary dynamics occur simultaneously within microbiomes; however, the implications for biogeochemistry under environmental change will depend on the timescales over which these processes contribute to a microbiome's response. Over the long term, evolution may play an increasingly important role for microbially driven biogeochemical responses to environmental change, particularly to conditions without recent historical precedent.

Nested PCR Approach for petB Gene Metabarcoding of Marine Synechococcus Populations

The molecular diversity of marine picocyanobacterial populations, an important component of phytoplankton communities, is better characterized using high-resolution marker genes than the 16S rRNA gene as they have greater sequence divergence to differentiate between closely related picocyanobacteria groups. Although specific ribosomal primers have been developed, another general disadvantage of bacterial ribosome-based diversity analyses is the variable number of rRNA gene copies. To overcome these issues, the single-copy petB gene, encoding the cytochrome b6 subunit of the cytochrome b6f complex, has been used as a high-resolution marker gene to characterize Synechococcus diversity. We have designed new primers targeting the petB gene and proposed a nested PCR method (termed Ong_2022) for metabarcoding of marine Synechococcus populations obtained by flow cytometry cell sorting. We evaluated the specificity and sensitivity of Ong_2022 against the standard amplification protocol (termed Mazard_2012) using filtered seawater samples. The Ong_2022 approach was also tested on flow cytometry-sorted Synechococcus populations. Samples (filtered and sorted) were obtained in the Southwest Pacific Ocean, from subtropical (ST) and subantarctic (SA) water masses. The two PCR approaches using filtered samples recovered the same dominant subclades, Ia, Ib, IVa, and IVb, with small differences in relative abundance across the distinct samples. For example, subclade IVa was dominant in ST samples with the Mazard_2012 approach, while the same samples processed with Ong_2022 showed similar contributions of subclades IVa and Ib to the total community. The Ong_2022 approach generally captured a higher genetic diversity of Synechococcus subcluster 5.1 than the Mazard_2012 approach while having a lower proportion of incorrectly assigned amplicon sequence variants (ASVs). All flow cytometry-sorted Synechococcus samples could be amplified only by our nested approach. The taxonomic diversity obtained with our primers on both sample types was in agreement with the clade distribution observed by previous studies that applied other marker genes or PCR-free metagenomic approaches under similar environmental conditions. IMPORTANCE The petB gene has been proposed as a high-resolution marker gene to access the diversity of marine Synechococcus populations. A systematic metabarcoding approach based on the petB gene would improve the characterization/assessment of the Synechococcus community structure in marine planktonic ecosystems. We have designed and tested specific primers to be applied in a nested PCR protocol (Ong_2022) for metabarcoding the petB gene. The Ong_2022 protocol can be applied to samples with low DNA content, such as those obtained by flow cytometry cell sorting, allowing the simultaneous assessment of the genetic diversity of Synechococcus populations and cellular properties and activities (e.g., nutrient cell ratios or carbon uptake rates). Our approach will allow future studies using flow cytometry to investigate the link between ecological traits and taxonomic diversity of marine Synechococcus.

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.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml^-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml^-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

A Large Silicon Pool in Small Picophytoplankton

Marine picophytoplankton (<2 μm) play a key role in supporting food web and energy flow in the ocean, and are major contributors to the global marine carbon (C) cycle. In recent years, picophytoplankton have been found to have significant silica (Si) accumulation, a finding which provides a new sight into the interaction of marine C and Si cycles and questions the overwhelming role of large diatoms (>2 μm) in the Si cycle. As picophytoplankton have high cell abundance and wide distribution in the open ocean, exploring their influences on the C and Si cycles as well as other element cycles are becoming new scientific hotspots. However, there are still few studies on the physiology and ecology of picophytoplankton, especially their potential roles in the biogeochemical Si cycle at present. Thus, it is necessary to accurately evaluate and quantify the contributions of picophytoplankton to the C and Si cycles, and to further understand their C and Si sinking mechanisms. In this review, we expect to have a novel understanding of picophytoplankton Si pool and regulation mechanism by conducting targeted studies on these scientific issues. This also provides a premise foundation and theoretical framework for further study of the role of small cells in the global ocean Si cycle and the coupling of C and Si cycles.

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.

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.

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.

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.

Differential Response of Bacterial Microdiversity to Simulated Global Change

Global change experiments often observe shifts in bacterial community composition based on 16S rRNA gene sequences. However, this genetic region can mask a large amount of genetic and phenotypic variation among bacterial strains sharing even identical 16S regions. As such, it remains largely unknown whether variation at the sub-16S level, sometimes termed microdiversity, responds to environmental perturbations and whether such changes are relevant to ecosystem processes. Here, we investigated microdiversity within Curtobacterium, the dominant bacterium found in the leaf litter layer of soil, to simulated drought and nitrogen addition in a field experiment. We first developed and validated Curtobacterium-specific primers of the groEL gene to assess microdiversity within this lineage. We then tracked the response of this microdiversity to simulated global change in two adjacent plant communities, grassland and coastal sage scrub (CSS). Curtobacterium microdiversity responded to drought but not nitrogen addition, indicating variation within the genus of drought tolerance but not nitrogen response. Further, the response of microdiversity to drought depended on the ecosystem, suggesting that litter substrate selects for a distinct composition of microdiversity that is constrained in its response, perhaps related to tradeoffs in resource acquisition traits. Supporting this interpretation, a metagenomic analysis revealed that the composition of Curtobacterium-encoded CAZymes varied distinctly across the two ecosystems. Identifying the degree to which relevant traits are phylogenetically conserved may help to predict when the aggregated response of a 16S-defined taxon masks differential responses of finer-scale bacterial diversity to global change. Importance Microbial communities play an integral role in global biogeochemical cycling, but our understanding of how global change will affect microbial community structure and functioning remains limited. Microbiome analyses typically aggregate large amounts of genetic diversity which may obscure finer variation in traits. This study found that fine-scale diversity (or microdiversity) within the bacterial genus, Curtobacterium, was affected by simulated global changes. However, the degree to which this was true depended on the type of global change as the composition of Curtobacterium microdiversity was affected by drought, but not by nitrogen addition. Further, these changes were associated with variation in carbon degradation traits. Future work might improve predictions of microbial community responses to global change by considering microdiversity.

Significance of gene variants for the functional biogeography of the near-surface Atlantic Ocean microbiome

Microbial communities are major drivers of global elemental cycles in the oceans due to their high abundance and enormous taxonomic and functional diversity. Recent studies assessed microbial taxonomic and functional biogeography in global oceans but microbial functional biogeography remains poorly studied. Here we show that in the near-surface Atlantic and Southern Ocean between 62°S and 47°N microbial communities exhibit distinct taxonomic and functional adaptations to regional environmental conditions. Richness and diversity showed maxima around 40° latitude and intermediate temperatures, especially in functional genes (KEGG-orthologues, KOs) and gene profiles. A cluster analysis yielded three clusters of KOs but five clusters of genes differing in the abundance of genes involved in nutrient and energy acquisition. Gene profiles showed much higher distance-decay rates than KO and taxonomic profiles. Biotic factors were identified as highly influential in explaining the observed patterns in the functional profiles, whereas temperature and biogeographic province mainly explained the observed taxonomic patterns. Our results thus indicate fine-tuned genetic adaptions of microbial communities to regional biotic and environmental conditions in the Atlantic and Southern Ocean.

The taxonomic and functional diversity of marine microbial communities are shaped by both environmental and biotic factors. Here, the authors investigate the functional biogeography of epipelagic prokaryotic communities along a 13,000-km transect in the Southern and Atlantic Oceans, showing finely tuned genetic adaptations to regional conditions.

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.

Microdiversity characterizes prevalent phylogenetic clades in the glacier-fed stream microbiome

Glacier-fed streams (GFSs) are extreme and rapidly vanishing ecosystems, and yet they harbor diverse microbial communities. Although our understanding of the GFS microbiome has recently increased, we do not know which microbial clades are ecologically successful in these ecosystems, nor do we understand potentially underlying mechanisms. Ecologically successful clades should be more prevalent across GFSs compared to other clades, which should be reflected as clade-wise distinctly low phylogenetic turnover. However, methods to assess such patterns are currently missing. Here we developed and applied a novel analytical framework, “phyloscore analysis”, to identify clades with lower spatial phylogenetic turnover than other clades in the sediment microbiome across twenty GFSs in New Zealand. These clades constituted up to 44% and 64% of community α-diversity and abundance, respectively. Furthermore, both their α-diversity and abundance increased as sediment chlorophyll a decreased, corroborating their ecological success in GFS habitats largely devoid of primary production. These clades also contained elevated levels of putative microdiversity than others, which could potentially explain their high prevalence in GFSs. This hitherto unknown microdiversity may be threatened as glaciers shrink, urging towards further genomic and functional exploration of the GFS microbiome.

A drift-barrier model drives the genomic landscape of a structured bacterial population

Bacterial populations differentiate over time and space to form distinct genetic units. The mechanisms governing this diversification are presumed to result from the ecological context of living units to adapt to specific niches. Recently, a model assuming the acquisition of advantageous genes among populations rather than whole genome sweeps has emerged to explain population differentiation. However, the characteristics of these exchanged, or flexible, genes and whether their evolution is driven by adaptive or neutral processes remain controversial. By analysing the flexible genome of single-amplified genomes of co-occurring populations of the marine Prochlorococcus HLII ecotype, we highlight that genomic compartments - rather than population units - are characterized by different evolutionary trajectories. The dynamics of gene fluxes vary across genomic compartments and therefore the effectiveness of selection depends on the fluctuation of the effective population size along the genome. Taken together, these results support the drift-barrier model of bacterial evolution.

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.

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

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus are the most abundant photosynthetic organisms on Earth, an ecological success thought to be linked to the differential partitioning of distinct ecotypes into specific ecological niches. However, the underlying processes that governed the diversification of these microorganisms and the appearance of niche-related phenotypic traits are just starting to be elucidated. Here, by comparing 81 genomes, including 34 new Synechococcus, we explored the evolutionary processes that shaped the genomic diversity of picocyanobacteria. Time-calibration of a core-protein tree showed that gene gain/loss occurred at an unexpectedly low rate between the different lineages, with for instance 5.6 genes gained per million years (My) for the major Synechococcus lineage (sub-cluster 5.1), among which only 0.71/My have been fixed in the long term. Gene content comparisons revealed a number of candidates involved in nutrient adaptation, a large proportion of which are located in genomic islands shared between either closely or more distantly related strains, as identified using an original network construction approach. Interestingly, strains representative of the different ecotypes co-occurring in phosphorus-depleted waters (Synechococcus clades III, WPC1, and sub-cluster 5.3) were shown to display different adaptation strategies to this limitation. In contrast, we found few genes potentially involved in adaptation to temperature when comparing cold and warm thermotypes. Indeed, comparison of core protein sequences highlighted variants specific to cold thermotypes, notably involved in carotenoid biosynthesis and the oxidative stress response, revealing that long-term adaptation to thermal niches relies on amino acid substitutions rather than on gene content variation. Altogether, this study not only deciphers the respective roles of gene gains/losses and sequence variation but also uncovers numerous gene candidates likely involved in niche partitioning of two key members of the marine phytoplankton.

Linking regional shifts in microbial genome adaptation with surface ocean biogeochemistry

Linking 'omics measurements with biogeochemical cycles is a widespread challenge in microbial community ecology. Here, we propose applying genomic adaptation as 'biosensors' for microbial investments to overcome nutrient stress. We then integrate this genomic information with a trait-based model to predict regional shifts in the elemental composition of marine plankton communities. We evaluated this approach using metagenomic and particulate organic matter samples from the Atlantic, Indian and Pacific Oceans. We find that our genome-based trait model significantly improves our prediction of particulate C : P (carbon : phosphorus) across ocean regions. Furthermore, we detect previously unrecognized ocean areas of iron, nitrogen and phosphorus stress. In many ecosystems, it can be very challenging to quantify microbial stress. Thus, a carefully calibrated genomic approach could become a widespread tool for understanding microbial responses to environmental changes and the biogeochemical outcomes. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

Unveiling membrane thermoregulation strategies in marine picocyanobacteria

The wide latitudinal distribution of marine Synechococcus cyanobacteria partly relies on the differentiation of lineages adapted to distinct thermal environments. Membranes are highly thermosensitive cell components, and the ability to modulate their fluidity can be critical for the fitness of an ecotype in a particular thermal niche. We compared the thermophysiology of Synechococcus strains representative of major temperature ecotypes in the field. We measured growth, photosynthetic capacities and membrane lipidome variations. We carried out a metagenomic analysis of stations of the Tara Oceans expedition to describe the latitudinal distribution of the lipid desaturase genes in the oceans. All strains maintained efficient photosynthetic capacities over their different temperature growth ranges. Subpolar and cold temperate strains showed enhanced capacities for lipid monodesaturation at low temperature thanks to an additional, poorly regiospecific Δ9-desaturase. By contrast, tropical and warm temperate strains displayed moderate monodesaturation capacities but high proportions of double unsaturations in response to cold, thanks to regiospecific Δ12-desaturases. The desaturase genes displayed specific distributions directly related to latitudinal variations in ocean surface temperature. This study highlights the critical importance of membrane fluidity modulation by desaturases in the adaptive strategies of Synechococcus cyanobacteria during the colonization of novel thermal niches.

Microdiversity ensures the maintenance of functional microbial communities under changing environmental conditions

Microdiversity can lead to different ecotypes within the same species. These are assumed to provide stability in time and space to those species. However, the role of microdiversity in the stability of whole microbial communities remains underexplored. Understanding the drivers of microbial community stability is necessary to predict community response to future disturbances. Here, we analyzed 16S rRNA gene amplicons from eight different temperate bog lakes at the 97% OTU and amplicon sequence variant (ASV) levels and found ecotypes within the same OTU with different distribution patterns in space and time. We observed that these ecotypes are adapted to different values of environmental factors such as water temperature and oxygen concentration. Our results showed that the existence of several ASVs within a OTU favored its persistence across changing environmental conditions. We propose that microdiversity aids the stability of microbial communities in the face of fluctuations in environmental factors.

Biogeochemical controls of surface ocean phosphate

Surface ocean phosphate is commonly below the standard analytical detection limits, leading to an incomplete picture of the global variation and biogeochemical role of phosphate. A global compilation of phosphate measured using high-sensitivity methods revealed several previously unrecognized low-phosphate areas and clear regional differences. Both observational climatologies and Earth system models (ESMs) systematically overestimated surface phosphate. Furthermore, ESMs misrepresented the relationships between phosphate, phytoplankton biomass, and primary productivity. Atmospheric iron input and nitrogen fixation are known important controls on surface phosphate, but model simulations showed that differences in the iron-to-macronutrient ratio in the vertical nutrient supply and surface lateral transport are additional drivers of phosphate concentrations. Our study demonstrates the importance of accurately quantifying nutrients for understanding the regulation of ocean ecosystems and biogeochemistry now and under future climate conditions.

Differential expression of the glucose transporter gene glcH in response to glucose and light in marine picocyanobacteria

Background

Our team discovered that Prochlorococcus can take up glucose, in a process that changes the transcriptional pattern of several genes involved in glucose metabolization. We have also shown that glcH encodes a very high affinity glucose transporter, and that glucose is taken up by natural Prochlorococcus populations. We demonstrated that the kinetic parameters of glucose uptake show significant diversity in different Prochlorococcus and Synechococcus strains. Here, we tested whether the transcriptional response of glcH to several glucose concentrations and light conditions was also different depending on the studied strain.

Methods

Cultures were grown in the light, supplemented with five different glucose concentrations or subjected to darkness, and cells harvested after 24 h of treatment. qRT-PCR was used to determine glcH expression in four Prochlorococcus and two Synechococcus strains.

Results

In all studied strains glcH was expressed in the absence of glucose, and it increased upon glucose addition to cultures. The changes differed depending on the strain, both in the magnitude and in the way cells responded to the tested glucose concentrations. Unlike the other strains, Synechococcus BL107 showed the maximum glucose uptake at 5 nM glucose. Darkness induced a strong decrease in glcH expression, especially remarkable in Prochlorococcus MIT9313.

Discussion

Our results suggest that marine picocyanobacteria are actively monitoring the availability of glucose, to upregulate glcH expression in order to exploit the presence of sugars in the environment. The diverse responses observed in different strains suggest that the transcriptional regulation of glucose uptake has been adjusted by evolutive selection. Darkness promotes a strong decrease in glcH expression in all studied strains, which fits with previous results on glucose uptake in Prochlorococcus. Overall, this work reinforces the importance of mixotrophy for marine picocyanobacteria.