Genomes reveal nutrient stress patterns: Within the surface ocean, nitrogen, iron, and phosphorous can all be limiting nutrients for phytoplankton depending on location. Ustick et al. used the prevalence of Prochlorococcus genes involved in nutrient acquisition to develop maps of inferred nutrient stress across the global ocean (see the Perspective by Coleman). They found broad patterns of limitation consistent with an Earth system model and nutrient addition experiments. Leveraging metagenomic data in this manner is an appealing approach that will help to expand our understanding of the biogeochemistry in the vast open ocean. Science , this issue p. 287 ; see also p. 239

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.

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.

The distribution of major clades of Prochlorococcus tracks light, temperature and other environmental variables; yet, the drivers of genomic diversity within these ecotypes and the net effect on biodiversity of the larger community are poorly understood. We examined high light (HL) adapted Prochlorococcus communities across spatial and temporal environmental gradients in the Pacific Ocean to determine the ecological drivers of population structure and diversity across taxonomic ranks. We show that the Prochlorococcus community has the highest diversity at low latitudes, but seasonality driven by temperature, day length and nutrients adds complexity. At finer taxonomic resolution, some ‘sub-ecotype’ clades have unique, cohesive responses to environmental variables and distinct biogeographies, suggesting that presently defined ecotypes can be further partitioned into ecologically meaningful units. Intriguingly, biogeographies of the HL-I sub-ecotypes are driven by unique combinations of environmental traits, rather than through trait hierarchy, while the HL-II sub-ecotypes appear ecologically similar, thus demonstrating differences among these dominant HL ecotypes. Examining biodiversity across taxonomic ranks reveals high-resolution dynamics of Prochlorococcus evolution and ecology that are masked at phylogenetically coarse resolution. Spatial and seasonal trends of Prochlorococcus communities suggest that the future ocean may be comprised of different populations, with implications for ecosystem structure and function.

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.