The globally abundant marine Cyanobacteria Prochlorococcus and Synechococcus share many physiological traits but presumably have different evolutionary histories and associated phylogeography. In Prochlorococcus, there is a clear phylogenetic hierarchy of ecotypes, whereas multiple Synechococcus clades have overlapping physiologies and environmental distributions. However, microbial traits are associated with different phylogenetic depths. Using this principle, we reclassified diversity at different phylogenetic levels and compared the phylogeography. We sequenced the genetic diversity of Prochlorococcus and Synechococcus from 339 samples across the tropical Pacific Ocean and North Atlantic Ocean using a highly variable phylogenetic marker gene (rpoC1). We observed clear parallel niche distributions of ecotypes leading to high Pianka’s Index values driven by distinct shifts at two transition points. The first transition point at 6°N distinguished ecotypes adapted to warm waters but separated by macronutrient content. At 39°N, ecotypes adapted to warm, low macronutrient vs. colder, high macronutrient waters shifted. Finally, we detected parallel vertical and regional single-nucleotide polymorphism microdiversity within clades from both Prochlorococcus and Synechococcus, suggesting uniquely adapted populations at very specific depths, as well as between the Atlantic and Pacific Oceans. Overall, this study demonstrates that Prochlorococcus and Synechococcus have shared phylogenetic organization of traits and associated phylogeography.
Niche theory suggests that the realized niche occupied by an organism in the field is a subset of the fundamental niche space of the organism, absent additional biotic and abiotic factors. Though often assumed, this discrepancy is rarely tested for specific organisms, and could act as a source of error in model predictions of biogeographical shifts resulting from temperature change which assume niche theory constraints. Here, we quantify the difference between fundamental and realized temperature niches for four dominant ecotypes of
- Award ID(s):
- 1937715
- NSF-PAR ID:
- 10449328
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Limnology and Oceanography
- Volume:
- 66
- Issue:
- 7
- ISSN:
- 0024-3590
- Page Range / eLocation ID:
- p. 2632-2647
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract -
Abstract 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.
-
Summary Currently defined ecotypes in marine cyanobacteria
Prochlorococcus andSynechococcus likely contain subpopulations that themselves are ecologically distinct. We developed and applied high‐throughput sequencing for the 16S‐23S rRNA internally transcribed spacer (ITS) to examine ecotype and fine‐scale genotypic community dynamics for monthly surface water samples spanning 5 years at the San Pedro Ocean Time‐series site. Ecotype‐level structure displayed regular seasonal patterns including succession, consistent with strong forcing by seasonally varying abiotic parameters (e.g. temperature, nutrients, light). We identified tens to thousands of amplicon sequence variants (ASVs) within ecotypes, many of which exhibited distinct patterns over time, suggesting ecologically distinct populations within ecotypes. Community structure within some ecotypes exhibited regular, seasonal patterns, but not for others, indicating other more irregular processes such as phage interactions are important. Network analysis including T4‐like phage genotypic data revealed distinct viral variants correlated with different groups of cyanobacterial ASVs including time‐lagged predator–prey relationships. Variation partitioning analysis indicated that phage community structure more strongly explains cyanobacterial community structure at the ASV level than the abiotic environmental factors. These results support a hierarchical model whereby abiotic environmental factors more strongly shape niche partitioning at the broader ecotype level while phage interactions are more important in shaping community structure of fine‐scale variants within ecotypes. -
Abstract 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.
-
Multiple abrupt warming events (“hyperthermals”) punctuated the Early Eocene and were associated with deep-sea temperature increases of 2 to 4 °C, seafloor carbonate dissolution, and negative carbon isotope (δ13C) excursions. Whether hyperthermals were associated with changes in the global ocean overturning circulation is important for understanding their driving mechanisms and feedbacks and for gaining insight into the circulation’s sensitivity to climatic warming. Here, we present high-resolution benthic foraminiferal stable isotope records (δ13C and δ18O) throughout the Early Eocene Climate Optimum (~53.26 to 49.14 Ma) from the deep equatorial and North Atlantic. Combined with existing records from the South Atlantic and Pacific, these indicate consistently amplified δ13C excursion sizes during hyperthermals in the deep equatorial Atlantic. We compare these observations with results from an intermediate complexity Earth system model to demonstrate that this spatial pattern of δ13C excursion size is a predictable consequence of global warming-induced changes in ocean overturning circulation. In our model, transient warming drives the weakening of Southern Ocean-sourced overturning circulation, strengthens Atlantic meridional water mass aging gradients, and amplifies the magnitude of negative δ13C excursions in the equatorial to North Atlantic. Based on model-data consistency, we conclude that Eocene hyperthermals coincided with repeated weakening of the global overturning circulation. Not accounting for ocean circulation impacts on δ13C excursions will lead to incorrect estimates of the magnitude of carbon release driving hyperthermals. Our finding of weakening overturning in response to past transient climatic warming is consistent with predictions of declining Atlantic Ocean overturning strength in our warm future.