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.

Variation in ocean C:N:P of particulate organic matter (POM) has led to competing hypotheses for the underlying drivers. Each hypothesis predicts C:N:P equally well due to regional co-variance in environmental conditions and biodiversity. The Indian Ocean offers a unique positive temperature and nutrient supply relationship to test these hypotheses. Here we show how elemental concentrations and ratios vary over daily and regional scales. POM concentrations were lowest in the southern gyre, elevated across the equator, and peaked in the Bay of Bengal. Elemental ratios were highest in the gyre, but approached Redfield proportions northwards. As Prochlorococcus dominated the phytoplankton community, biodiversity changes could not explain the elemental variation. Instead, our data supports the nutrient supply hypothesis. Finally, gyre dissolved iron concentrations suggest extensive iron stress, leading to depressed ratios compared to other gyres. We propose a model whereby differences in iron supply and N2-fixation influence C:N:P levels across ocean gyres.

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’.

While the majority of studies use the environment to describe microbial populations, the high diversity of microbes can conversely be used as a resource to understand subtle environmental variability. Here, we used a high‐resolution spatial and temporal analysis of Prochlorococcus sp. in the Eastern Indian Ocean to determine whether ecotypes and microdiverse taxa can be used to identify fine‐scale biogeochemical regimes in this under‐studied region. A total of 246 DNA samples were collected every 4–6 h in April 2016 on GO‐SHIP cruise I09N, which transected gyre, equatorial, and monsoonal ecosystems between Western Australia and the Bay of Bengal. Using amplicon sequencing of the highly variable rpoC1 marker, we found that the region was largely dominated by the Prochlorococcus HL‐II clade. Conserved single nucleotide polymorphisms (SNPs) were used to identify four microdiverse haplotypes, or SNP‐delineated taxa, within the HL‐II clade of Prochlorococcus. The haplotypes showed regional patterns of relative gene count abundance that were significantly correlated with environmental conditions. Additionally, we used nonlinear least squares models to fit the sine wave function to our data and demonstrate that the haplotypes show distinct patterns in relative diel frequency, providing evidence that these microdiverse populations are ecologically and evolutionarily distinct. Overall, we show how the integration of a genomics data set into a biogeochemical framework can reveal a more nuanced understanding of a complex ocean basin.

No abstract

Climate-driven depletion of ocean oxygen strongly impacts the global cycles of carbon and nutrients as well as the survival of many animal species. One of the main uncertainties in predicting changes to marine oxygen levels is the regulation of the biological respiration demand associated with the biological pump. Derived from the Redfield ratio, the molar ratio of oxygen to organic carbon consumed during respiration (i.e., the respiration quotient, r−O2:C) is consistently assumed constant but rarely, if ever, measured. Using a prognostic Earth system model, we show that a 0.1 increase in the respiration quotient from 1.0 leads to a 2.3% decline in global oxygen, a large expansion of low-oxygen zones, additional water column denitrification of 38 Tg N/y, and the loss of fixed nitrogen and carbon production in the ocean. We then present direct chemical measurements of r−O2:C using a Pacific Ocean meridional transect crossing all major surface biome types. The observed r−O2:C has a positive correlation with temperature, and regional mean values differ significantly from Redfield proportions. Finally, an independent global inverse model analysis constrained with nutrients, oxygen, and carbon concentrations supports a positive temperature dependence of r−O2:C in exported organic matter. We provide evidence against the common assumption of a static biological link between the respiration of organic carbon and the consumption of oxygen. Furthermore, the model simulations suggest that a changing respiration quotient will impact multiple biogeochemical cycles and that future warming can lead to more intense deoxygenation than previously anticipated.

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.

Ocean temperatures will increase significantly over the next 100 years due to global climate change. As temperatures increase beyond current ranges, it is unclear how adaptation will impact the distribution and ecological role of marine microorganisms. To address this major unknown, we imposed a stressful high-temperature regime for 500 generations on a strain from the abundant marine Roseobacter clade. High-temperature-adapted isolates significantly improved their fitness but also increased biofilm formation at the air-liquid interface. Furthermore, this altered lifestyle was coupled with genomic changes linked to biofilm formation in individual isolates, and was also dominant in evolved populations. We hypothesize that the increasing biofilm formation was driven by lower oxygen availability at elevated temperature, and we observe a relative fitness increase at lower oxygen. The response is uniquely different from that of Escherichia coli adapted to high temperature (only 3% of mutated genes were shared in both studies). Thus, future increased temperatures could have a direct effect on organismal physiology and an indirect effect via a decrease in ocean oxygen solubility, leading to an alteration in microbial lifestyle.