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

Global biogeography of microbes driving ocean ecological status under climate change

Microbial communities play a crucial role in ocean ecology and global biogeochemical processes. However, understanding the intricate interactions among diversity, taxonomical composition, functional traits, and how these factors respond to climate change remains a significant challenge. Here, we propose seven distinct ecological statuses by systematically considering the diversity, structure, and biogeochemical potential of the ocean microbiome to delineate their biogeography. Anthropogenic climate change is expected to alter the ecological status of the surface ocean by influencing environmental conditions, particularly nutrient and oxygen contents. Our predictive model, which utilizes machine learning, indicates that the ecological status of approximately 32.44% of the surface ocean may undergo changes from the present to the end of this century, assuming no policy interventions. These changes mainly include poleward shifts in the main taxa, increases in photosynthetic carbon fixation and decreases in nutrient metabolism. However, this proportion can decrease significantly with effective control of greenhouse gas emissions. Our study underscores the urgent necessity for implementing policies to mitigate climate change, particularly from an ecological perspective.

Decoding populations in the ocean microbiome

Understanding the characteristics and structure of populations is fundamental to comprehending ecosystem processes and evolutionary adaptations. While the study of animal and plant populations has spanned a few centuries, microbial populations have been under scientific scrutiny for a considerably shorter period. In the ocean, analyzing the genetic composition of microbial populations and their adaptations to multiple niches can yield important insights into ecosystem function and the microbiome's response to global change. However, microbial populations have remained elusive to the scientific community due to the challenges associated with isolating microorganisms in the laboratory. Today, advancements in large-scale metagenomics and metatranscriptomics facilitate the investigation of populations from many uncultured microbial species directly from their habitats. The knowledge acquired thus far reveals substantial genetic diversity among various microbial species, showcasing distinct patterns of population differentiation and adaptations, and highlighting the significant role of selection in structuring populations. In the coming years, population genomics is expected to significantly increase our understanding of the architecture and functioning of the ocean microbiome, providing insights into its vulnerability or resilience in the face of ongoing global change. Video Abstract.

Genome-scale community modelling reveals conserved metabolic cross-feedings in epipelagic bacterioplankton communities

Marine microorganisms form complex communities of interacting organisms that influence central ecosystem functions in the ocean such as primary production and nutrient cycling. Identifying the mechanisms controlling their assembly and activities is a major challenge in microbial ecology. Here, we integrated Tara Oceans meta-omics data to predict genome-scale community interactions within prokaryotic assemblages in the euphotic ocean. A global genome-resolved co-activity network revealed a significant number of inter-lineage associations across diverse phylogenetic distances. Identified co-active communities include species displaying smaller genomes but encoding a higher potential for quorum sensing, biofilm formation, and secondary metabolism. Community metabolic modelling reveals a higher potential for interaction within co-active communities and points towards conserved metabolic cross-feedings, in particular of specific amino acids and group B vitamins. Our integrated ecological and metabolic modelling approach suggests that genome streamlining and metabolic auxotrophies may act as joint mechanisms shaping bacterioplankton community assembly in the global ocean surface.

Marine picoplankton metagenomes and MAGs from eleven vertical profiles obtained by the Malaspina Expedition

The Ocean microbiome has a crucial role in Earth's biogeochemical cycles. During the last decade, global cruises such as Tara Oceans and the Malaspina Expedition have expanded our understanding of the diversity and genetic repertoire of marine microbes. Nevertheless, there are still knowledge gaps regarding their diversity patterns throughout depth gradients ranging from the surface to the deep ocean. Here we present a dataset of 76 microbial metagenomes (MProfile) of the picoplankton size fraction (0.2-3.0 µm) collected in 11 vertical profiles covering contrasting ocean regions sampled during the Malaspina Expedition circumnavigation (7 depths, from surface to 4,000 m deep). The MProfile dataset produced 1.66 Tbp of raw DNA sequences from which we derived: 17.4 million genes clustered at 95% sequence similarity (M-GeneDB-VP), 2,672 metagenome-assembled genomes (MAGs) of Archaea and Bacteria (Malaspina-VP-MAGs), and over 100,000 viral genomic sequences. This dataset will be a valuable resource for exploring the functional and taxonomic connectivity between the photic and bathypelagic tropical and sub-tropical ocean, while increasing our general knowledge of the Ocean microbiome.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59

IMPORTANCE

These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

A Pipeline for Constructing Reference Genomes for Large Cohort-Specific Metagenome Compression

Metagenomic data compression is very important as metagenomic projects are facing the challenges of larger data volumes per sample and more samples nowadays. Reference-based compression is a promising method to obtain a high compression ratio. However, existing microbial reference genome databases are not suitable to be directly used as references for compression due to their large size and redundancy, and different metagenomic cohorts often have various microbial compositions. We present a novel pipeline that generated simplified and tailored reference genomes for large metagenomic cohorts, enabling the reference-based compression of metagenomic data. We constructed customized reference genomes, ranging from 2.4 to 3.9 GB, for 29 real metagenomic datasets and evaluated their compression performance. Reference-based compression achieved an impressive compression ratio of over 20 for human whole-genome data and up to 33.8 for all samples, demonstrating a remarkable 4.5 times improvement than the standard Gzip compression. Our method provides new insights into reference-based metagenomic data compression and has a broad application potential for faster and cheaper data transfer, storage, and analysis.

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.

Identification of microbial metabolic functional guilds from large genomic datasets

Heterotrophic microbes play an important role in the Earth System as key drivers of major biogeochemical cycles. Specifically, the consumption rate of organic matter is set by the interaction between diverse microbial communities and the chemical and physical environment in which they reside. Modeling these dynamics requires reducing the complexity of microbial communities and linking directly with biogeochemical functions. Microbial metabolic functional guilds provide one approach for reducing microbial complexity and incorporating microbial biogeochemical functions into models. However, we lack a way to identify these guilds. In this study, we present a method for defining metabolic functional guilds from annotated genomes, which are derived from both uncultured and cultured organisms. This method utilizes an Aspect Bernoulli (AB) model and was tested on three large genomic datasets with 1,733-3,840 genomes each. Ecologically relevant microbial metabolic functional guilds were identified including guilds related to DMSP degradation, dissimilatory nitrate reduction to ammonia, and motile copiotrophy. This method presents a way to generate hypotheses about functions co-occurring within individual microbes without relying on cultured representatives. Applying the concept of metabolic functional guilds to environmental samples will provide new insight into the role that heterotrophic microbial communities play in setting rates of carbon cycling.

The Landscape of Global Ocean Microbiome: From Bacterioplankton to Biofilms

The development of metagenomics has opened up a new era in the study of marine microbiota, which play important roles in biogeochemical cycles. In recent years, the global ocean sampling expeditions have spurred this research field toward a deeper understanding of the microbial diversities and functions spanning various lifestyles, planktonic (free-living) or sessile (biofilm-associated). In this review, we deliver a comprehensive summary of marine microbiome datasets generated in global ocean expeditions conducted over the last 20 years, including the Sorcerer II GOS Expedition, the Tara Oceans project, the bioGEOTRACES project, the Micro B3 project, the Bio-GO-SHIP project, and the Marine Biofilms. These datasets have revealed unprecedented insights into the microscopic life in our oceans and led to the publication of world-leading research. We also note the progress of metatranscriptomics and metaproteomics, which are confined to local marine microbiota. Furthermore, approaches to transforming the global ocean microbiome datasets are highlighted, and the state-of-the-art techniques that can be combined with data analyses, which can present fresh perspectives on marine molecular ecology and microbiology, are proposed.

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

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.