Linking regional shifts in microbial genome adaptation with surface ocean biogeochemistry

Proteome trait regulation of marine Synechococcus elemental stoichiometry under global change

Recent studies have demonstrated regional differences in marine ecosystem C:N:P with implications for carbon and nutrient cycles. Due to strong co-variance, temperature and nutrient stress explain variability in C:N:P equally well. A reductionistic approach can link changes in individual environmental drivers with changes in biochemical traits and cell C:N:P. Thus, we quantified effects of temperature and nutrient stress on Synechococcus chemistry using laboratory chemostats, chemical analyses, and data-independent acquisition mass spectrometry proteomics. Nutrient supply accounted for most C:N:Pcell variability and induced tradeoffs between nutrient acquisition and ribosomal proteins. High temperature prompted heat-shock, whereas thermal effects via the "translation-compensation hypothesis" were only seen under P-stress. A Nonparametric Bayesian Local Clustering algorithm suggested that changes in lipopolysaccharides, peptidoglycans, and C-rich compatible solutes may also contribute to C:N:P regulation. Physiological responses match field-based trends in ecosystem stoichiometry and suggest a hierarchical environmental regulation of current and future ocean C:N:P.

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.

Seasonal shifts in Fe-acquisition strategies in Southern Ocean microbial communities revealed by metagenomics and autonomous sampling

Iron (Fe) governs the cycling of organic carbon in large parts of the Southern Ocean. The strategies of diverse microbes to acquire the different chemical forms of Fe under seasonally changing organic carbon regimes remain, however, poorly understood. Here, we report high-resolution seasonal metagenomic observations from the region off Kerguelen Island (Indian Sector of the Southern Ocean) where natural Fe-fertilization induces consecutive spring and summer phytoplankton blooms. Our data illustrate pronounced, but distinct seasonal patterns in the abundance of genes implicated in the transport of different forms of Fe and organic substrates, of siderophore biosynthesis and carbohydrate-active enzymes. The seasonal dynamics suggest a temporal decoupling in the prokaryotic requirements of Fe and organic carbon during the spring phytoplankton bloom and a concerted access to these resources after the summer bloom. Taxonomic assignments revealed differences in the prokaryotic groups harbouring genes of a given Fe-related category and pronounced seasonal successions were observed. Using MAGs we could decipher the respective Fe- and organic substrate-related genes of individual taxa assigned to abundant groups. The ecological strategies related to Fe-acquisition provide insights on how this element could shape microbial community composition with potential implications on organic matter transformations in the Southern Ocean.

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.

Basin-scale biogeography of Prochlorococcus and SAR11 ecotype replication

Establishing links between microbial diversity and environmental processes requires resolving the high degree of functional variation among closely related lineages or ecotypes. Here, we implement and validate an improved metagenomic approach that estimates the spatial biogeography and environmental regulation of ecotype-specific replication patterns (RObs) across ocean regions. A total of 719 metagenomes were analyzed from meridional Bio-GO-SHIP sections in the Atlantic and Indian Ocean. Accounting for sequencing bias and anchoring replication estimates in genome structure were critical for identifying physiologically relevant biological signals. For example, ecotypes within the dominant marine cyanobacteria Prochlorococcus exhibited distinct diel cycles in RObs that peaked between 19:00-22:00. Additionally, both Prochlorococcus ecotypes and ecotypes within the highly abundant heterotroph Pelagibacter (SAR11) demonstrated systematic biogeographies in RObs that differed from spatial patterns in relative abundance. Finally, RObs was significantly regulated by nutrient stress and temperature, and explained by differences in the genomic potential for nutrient transport, energy production, cell wall structure, and replication. Our results suggest that our new approach to estimating replication is reflective of gross population growth. Moreover, this work reveals that the interaction between adaptation and environmental change drives systematic variability in replication patterns across ocean basins that is ecotype-specific, adding an activity-based dimension to our understanding of microbial niche space.

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. Metagenomically 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 metagenomically 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 .

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.

Functional responses of key marine bacteria to environmental change – toward genetic counselling for coastal waters

Coastal ecosystems deteriorate globally due to human-induced stress factors, like nutrient loading and pollution. Bacteria are critical to marine ecosystems, e.g., by regulating nutrient cycles, synthesizing vitamins, or degrading pollutants, thereby providing essential ecosystem services ultimately affecting economic activities. Yet, until now bacteria are overlooked both as mediators and indicators of ecosystem health, mainly due to methodological limitations in assessing bacterial ecosystem functions. However, these limitations are largely overcome by the advances in molecular biology and bioinformatics methods for characterizing the genetics that underlie functional traits of key bacterial populations – “key” in providing important ecosystem services, being abundant, or by possessing high metabolic rates. It is therefore timely to analyze and define the functional responses of bacteria to human-induced effects on coastal ecosystem health. We posit that categorizing the responses of key marine bacterial populations to changes in environmental conditions through modern microbial oceanography methods will allow establishing the nascent field of genetic counselling for our coastal waters. This requires systematic field studies of linkages between functional traits of key bacterial populations and their ecosystem functions in coastal seas, complemented with systematic experimental analyses of the responses to different stressors. Research and training in environmental management along with dissemination of results and dialogue with societal actors are equally important to ensure the role of bacteria is understood as fundamentally important for coastal ecosystems. Using the responses of microorganisms as a tool to develop genetic counselling for coastal ecosystems can ultimately allow for integrating bacteria as indicators of environmental change.

Global Ocean Particulate Organic Phosphorus, Carbon, Oxygen for Respiration, and Nitrogen (GO-POPCORN)

Concentrations and elemental stoichiometry of suspended particulate organic carbon, nitrogen, phosphorus, and oxygen demand for respiration (C:N:P:−O2) play a vital role in characterizing and quantifying marine elemental cycles. Here, we present Version 2 of the Global Ocean Particulate Organic Phosphorus, Carbon, Oxygen for Respiration, and Nitrogen (GO-POPCORN) dataset. Version 1 is a previously published dataset of particulate organic matter from 70 different studies between 1971 and 2010, while Version 2 is comprised of data collected from recent cruises between 2011 and 2020. The combined GO-POPCORN dataset contains 2673 paired surface POC/N/P measurements from 70°S to 73°N across all major ocean basins at high spatial resolution. Version 2 also includes 965 measurements of oxygen demand for organic carbon respiration. This new dataset can help validate and calibrate the next generation of global ocean biogeochemical models with flexible elemental stoichiometry. We expect that incorporating variable C:N:P:-O2 into models will help improve our estimates of key ocean biogeochemical fluxes such as carbon export, nitrogen fixation, and organic matter remineralization.

Global patterns and predictors of C:N:P in marine ecosystems

Oceanic nutrient cycles are coupled, yet carbon-nitrogen-phosphorus (C:N:P) stoichiometry in marine ecosystems is variable through space and time, with no clear consensus on the controls on variability. Here, we analyze hydrographic, plankton genomic diversity, and particulate organic matter data from 1970 stations sampled during a global ocean observation program (Bio-GO-SHIP) to investigate the biogeography of surface ocean particulate organic matter stoichiometry. We find latitudinal variability in C:N:P stoichiometry, with surface temperature and macronutrient availability as strong predictors of stoichiometry at high latitudes. Genomic observations indicated community nutrient stress and suggested that nutrient supply rate and nitrogen-versus-phosphorus stress are predictive of hemispheric and regional variations in stoichiometry. Our data-derived statistical model suggests that C:P and N:P ratios will increase at high latitudes in the future, however, changes at low latitudes are uncertain. Our findings suggest systematic regulation of elemental stoichiometry among ocean ecosystems, but that future changes remain highly uncertain.

Ecological modelling approaches for predicting emergent properties in microbial communities

Recent studies have brought forward the critical role of emergent properties in shaping microbial communities and the ecosystems of which they are a part. Emergent properties-patterns or functions that cannot be deduced linearly from the properties of the constituent parts-underlie important ecological characteristics such as resilience, niche expansion and spatial self-organization. While it is clear that emergent properties are a consequence of interactions within the community, their non-linear nature makes mathematical modelling imperative for establishing the quantitative link between community structure and function. As the need for conservation and rational modulation of microbial ecosystems is increasingly apparent, so is the consideration of the benefits and limitations of the approaches to model emergent properties. Here we review ecosystem modelling approaches from the viewpoint of emergent properties. We consider the scope, advantages and limitations of Lotka-Volterra, consumer-resource, trait-based, individual-based and genome-scale metabolic models. Future efforts in this research area would benefit from capitalizing on the complementarity between these approaches towards enabling rational modulation of complex microbial ecosystems.

Siderophores as an iron source for picocyanobacteria in deep chlorophyll maximum layers of the oligotrophic ocean

Prochlorococcus and Synechococcus are the most abundant photosynthesizing organisms in the oceans. Gene content variation among picocyanobacterial populations in separate ocean basins often mirrors the selective pressures imposed by the region's distinct biogeochemistry. By pairing genomic datasets with trace metal concentrations from across the global ocean, we show that the genomic capacity for siderophore-mediated iron uptake is widespread in Synechococcus and low-light adapted Prochlorococcus populations from deep chlorophyll maximum layers of iron-depleted regions of the oligotrophic Pacific and S. Atlantic oceans: Prochlorococcus siderophore consumers were absent in the N. Atlantic ocean (higher new iron flux) but constituted up to half of all Prochlorococcus genomes from metagenomes in the N. Pacific (lower new iron flux). Picocyanobacterial siderophore consumers, like many other bacteria with this trait, also lack siderophore biosynthesis genes indicating that they scavenge exogenous siderophores from seawater. Statistical modeling suggests that the capacity for siderophore uptake is endemic to remote ocean regions where atmospheric iron fluxes are the smallest, especially at deep chlorophyll maximum and primary nitrite maximum layers. We argue that abundant siderophore consumers at these two common oceanographic features could be a symptom of wider community iron stress, consistent with prior hypotheses. Our results provide a clear example of iron as a selective force driving the evolution of marine picocyanobacteria.

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.

Evolution of Phytoplankton as Estimated from Genetic Diversity

Phytoplankton are photosynthetic, single-celled organisms producing almost half of all oxygen on Earth and play a central role as prey for higher organisms, making them irreplaceable in the marine food web. As Global Change proceeds, imposing rapidly intensifying selection pressures, phytoplankton are forced to undergo evolution, local extinction, or redistribution, with potentially cascading effects throughout the marine ecosystem. Recent results from the field of population genetics display high levels of standing genetic diversity in natural phytoplankton populations, providing ample ‘evolutionary options’ and implying high adaptive potential to changing conditions. This potential for adaptive evolution is realized in several studies of experimental evolution, even though most of these studies investigate the evolution of only single strains. This, however, shows that phytoplankton not only evolve from standing genetic diversity, but also rely on de novo mutations. Recent global sampling campaigns show that the immense intraspecific diversity of phytoplankton in the marine ecosystem has been significantly underestimated, meaning we are only studying a minor portion of the relevant variability in the context of Global Change and evolution. An increased understanding of genomic diversity is primarily hampered by the low number of ecologically representative reference genomes of eukaryotic phytoplankton and the functional annotation of these. However, emerging technologies relying on metagenome and transcriptome data may offer a more realistic understanding of phytoplankton diversity.

Environmental control of marine phytoplankton stoichiometry in the North Atlantic Ocean

The stoichiometric coupling of carbon to limiting nutrients in marine phytoplankton regulates the magnitude of biological carbon sequestration in the ocean. While clear links between plankton C:N ratios and environmental drivers have been identified, the nature and direction of these links, as well as their underlying physiological and ecological controls, remain uncertain. We show, with a well-constrained mechanistic model of plankton ecophysiology, that while nitrogen availability and temperature emerge as the main drivers of phytoplankton C:N stoichiometry in the North Atlantic, the biological mechanisms involved vary depending on the spatiotemporal scale and region considered. We find that phytoplankton C:N stoichiometry is overall controlled by nitrogen availability below 40° N, predominantly driven by ecoevolutionary shifts in the functional composition of the phytoplankton communities, while phytoplankton stoichiometric plasticity in response to dropping temperatures and increased grazing pressure dominates at higher latitudes. Our findings highlight the potential of "organisms-to-ecosystems" modeling approaches based on mechanistic models of plankton biology accounting for physiology, ecology, and trait evolution to explore and explain complex observational data and ultimately improve the predictions of global ocean models.

High spatial resolution global ocean metagenomes from Bio-GO-SHIP repeat hydrography transects

Detailed descriptions of microbial communities have lagged far behind physical and chemical measurements in the marine environment. Here, we present 971 globally distributed surface ocean metagenomes collected at high spatio-temporal resolution. Our low-cost metagenomic sequencing protocol produced 3.65 terabases of data, where the median number of base pairs per sample was 3.41 billion. The median distance between sampling stations was 26 km. The metagenomic libraries described here were collected as a part of a biological initiative for the Global Ocean Ship-based Hydrographic Investigations Program, or "Bio-GO-SHIP." One of the primary aims of GO-SHIP is to produce high spatial and vertical resolution measurements of key state variables to directly quantify climate change impacts on ocean environments. By similarly collecting marine metagenomes at high spatiotemporal resolution, we expect that this dataset will help answer questions about the link between microbial communities and biogeochemical fluxes in a changing ocean.

Metagenomic analysis reveals global-scale patterns of ocean nutrient limitation

Nutrient supply regulates the activity of phytoplankton, but the global biogeography of nutrient limitation and co-limitation is poorly understood. Prochlorococcus adapt to local environments by gene gains and losses, and we used genomic changes as an indicator of adaptation to nutrient stress. We collected metagenomes from all major ocean regions as part of the Global Ocean Ship-based Hydrographic Investigations Program (Bio-GO-SHIP) and quantified shifts in genes involved in nitrogen, phosphorus, and iron assimilation. We found regional transitions in stress type and severity as well as widespread co-stress. Prochlorococcus stress genes, bottle experiments, and Earth system model predictions were correlated. We propose that the biogeography of multinutrient stress is stoichiometrically linked by controls on nitrogen fixation. Our omics-based description of phytoplankton resource use provides a nuanced and highly resolved description of nutrient stress in the global ocean.

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.