Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome.

Ocean microbial communities strongly influence the biogeochemistry, food webs, and climate of our planet. Despite recent advances in understanding their taxonomic and genomic compositions, little is known about how their transcriptomes vary globally. Here, we present a dataset of 187 metatranscriptomes and 370 metagenomes from 126 globally distributed sampling stations and establish a resource of 47 million genes to study community-level transcriptomes across depth layers from pole-to-pole. We examine gene expression changes and community turnover as the underlying mechanisms shaping community transcriptomes along these axes of environmental variation and show how their individual contributions differ for multiple biogeochemically relevant processes. Furthermore, we find the relative contribution of gene expression changes to be significantly lower in polar than in non-polar waters and hypothesize that in polar regions, alterations in community activity in response to ocean warming will be driven more strongly by changes in organismal composition than by gene regulatory mechanisms. VIDEO ABSTRACT.

Plankton networks driving carbon export in the oligotrophic ocean.

The biological carbon pump is the process by which CO2 is transformed to organic carbon via photosynthesis, exported through sinking particles, and finally sequestered in the deep ocean. While the intensity of the pump correlates with plankton community composition, the underlying ecosystem structure driving the process remains largely uncharacterized. Here we use environmental and metagenomic data gathered during the Tara Oceans expedition to improve our understanding of carbon export in the oligotrophic ocean. We show that specific plankton communities, from the surface and deep chlorophyll maximum, correlate with carbon export at 150 m and highlight unexpected taxa such as Radiolaria and alveolate parasites, as well as Synechococcus and their phages, as lineages most strongly associated with carbon export in the subtropical, nutrient-depleted, oligotrophic ocean. Additionally, we show that the relative abundance of a few bacterial and viral genes can predict a significant fraction of the variability in carbon export in these regions.

Impact of Mesoscale Eddies on Deep Chlorophyll Maxima.

Deep Chlorophyll Maxima (DCM) are ubiquitous features in stratified oceanic systems. Their establishment and maintenance result from hydrographical stability favoring specific environmental conditions with respect to light and nutrient availability required for phytoplankton growth. This stability can potentially be challenged by mesoscale eddies impacting the water column's vertical structure and thus the environmental parameters that condition the subsistence of DCMs. Here, data from the global BGC-Argo float network are collocated with mesoscale eddies to explore their impact on DCMs. We show that cyclonic eddies, by providing optimal light and nutrient conditions, increase the occurrence of DCMs characterized by Deep Biomass Maxima for phytoplankton. In contrast, DCMs in anticyclonic eddies seem to be driven by photoacclimation as they coincide with Deep Acclimation Maxima without biomass accumulation. These findings suggest that the two types of eddies potentially have different impacts on the role of DCMs in global primary production.

A holistic approach to marine eco-systems biology.

The structure, robustness, and dynamics of ocean plankton ecosystems remain poorly understood due to sampling, analysis, and computational limitations. The Tara Oceans consortium organizes expeditions to help fill this gap at the global level.

The South Atlantic Dipole via multichannel singular spectrum analysis.

This study analyzes coupled atmosphere-ocean variability in the South Atlantic Ocean. To do so, we characterize the spatio-temporal variability of annual mean sea-surface temperature (SST) and sea-level pressure (SLP) using Multichannel Singular Spectrum Analysis (M-SSA). We applied M-SSA to ERA5 reanalysis data (1959-2022) of South Atlantic SST and SLP, both individually and jointly, and identified a nonlinear trend, as well as two climate oscillations. The leading oscillation, with a period of 13 years, consists of a basin-wide southwest-northeast dipole and is observed both in the individual variables and in the coupled analysis. This mode is reminiscent of the already known South Atlantic Dipole, and it is probably related to the Pacific Decadal Oscillation and to El Niño-Southern Oscillation in the Pacific Ocean. The second oscillation has a 5-year period and also displays a dipolar structure. The main difference between the spatial structure of the decadal, 13-year, and the interannual, 5-year mode is that, in the first one, the SST cold tongue region in the southeast Atlantic's Cape Basin is included in the pole closer to the equator. Together, these two oscillatory modes, along with the trend, capture almost 40% of the total interannual variability of the SST and SLP fields, and of their co-variability. These results provide further insights into the spatio-temporal evolution of SST and SLP variability in the South Atlantic, in particular as it relates to the South Atlantic Dipole and its predictability.

Exploring Microdiversity in Novel Kordia sp. (Bacteroidetes) with Proteorhodopsin from the Tropical Indian Ocean via Single Amplified Genomes.

Marine Bacteroidetes constitute a very abundant bacterioplankton group in the oceans that plays a key role in recycling particulate organic matter and includes several photoheterotrophic members containing proteorhodopsin. Relatively few marine Bacteroidetes species have been described and, moreover, they correspond to cultured isolates, which in most cases do not represent the actual abundant or ecologically relevant microorganisms in the natural environment. In this study, we explored the microdiversity of 98 Single Amplified Genomes (SAGs) retrieved from the surface waters of the underexplored North Indian Ocean, whose most closely related isolate is Kordia algicida OT-1. Using Multi Locus Sequencing Analysis (MLSA) we found no microdiversity in the tested conserved phylogenetic markers (16S rRNA and 23S rRNA genes), the fast-evolving Internal Transcribed Spacer and the functional markers proteorhodopsin and the beta-subunit of RNA polymerase. Furthermore, we carried out a Fragment Recruitment Analysis (FRA) with marine metagenomes to learn about the distribution and dynamics of this microorganism in different locations, depths and size fractions. This analysis indicated that this taxon belongs to the rare biosphere, showing its highest abundance after upwelling-induced phytoplankton blooms and sinking to the deep ocean with large organic matter particles. This uncultured Kordia lineage likely represents a novel Kordia species (Kordia sp. CFSAG39SUR) that contains the proteorhodopsin gene and has a widespread spatial and vertical distribution. The combination of SAGs and MLSA makes a valuable approach to infer putative ecological roles of uncultured abundant microorganisms.

EUREC4A: A Field Campaign to Elucidate the Couplings Between Clouds, Convection and Circulation.

Trade-wind cumuli constitute the cloud type with the highest frequency of occurrence on Earth, and it has been shown that their sensitivity to changing environmental conditions will critically influence the magnitude and pace of future global warming. Research over the last decade has pointed out the importance of the interplay between clouds, convection and circulation in controling this sensitivity. Numerical models represent this interplay in diverse ways, which translates into different responses of trade-cumuli to climate perturbations. Climate models predict that the area covered by shallow cumuli at cloud base is very sensitive to changes in environmental conditions, while process models suggest the opposite. To understand and resolve this contradiction, we propose to organize a field campaign aimed at quantifying the physical properties of trade-cumuli (e.g., cloud fraction and water content) as a function of the large-scale environment. Beyond a better understanding of clouds-circulation coupling processes, the campaign will provide a reference data set that may be used as a benchmark for advancing the modelling and the satellite remote sensing of clouds and circulation. It will also be an opportunity for complementary investigations such as evaluating model convective parameterizations or studying the role of ocean mesoscale eddies in air-sea interactions and convective organization.

Open science resources for the discovery and analysis of Tara Oceans data.

The Tara Oceans expedition (2009-2013) sampled contrasting ecosystems of the world oceans, collecting environmental data and plankton, from viruses to metazoans, for later analysis using modern sequencing and state-of-the-art imaging technologies. It surveyed 210 ecosystems in 20 biogeographic provinces, collecting over 35,000 samples of seawater and plankton. The interpretation of such an extensive collection of samples in their ecological context requires means to explore, assess and access raw and validated data sets. To address this challenge, the Tara Oceans Consortium offers open science resources, including the use of open access archives for nucleotides (ENA) and for environmental, biogeochemical, taxonomic and morphological data (PANGAEA), and the development of on line discovery tools and collaborative annotation tools for sequences and images. Here, we present an overview of Tara Oceans Data, and we provide detailed registries (data sets) of all campaigns (from port-to-port), stations and sampling events.

Ocean plankton. Patterns and ecological drivers of ocean viral communities.

Viruses influence ecosystems by modulating microbial population size, diversity, metabolic outputs, and gene flow. Here, we use quantitative double-stranded DNA (dsDNA) viral-fraction metagenomes (viromes) and whole viral community morphological data sets from 43 Tara Oceans expedition samples to assess viral community patterns and structure in the upper ocean. Protein cluster cataloging defined pelagic upper-ocean viral community pan and core gene sets and suggested that this sequence space is well-sampled. Analyses of viral protein clusters, populations, and morphology revealed biogeographic patterns whereby viral communities were passively transported on oceanic currents and locally structured by environmental conditions that affect host community structure. Together, these investigations establish a global ocean dsDNA viromic data set with analyses supporting the seed-bank hypothesis to explain how oceanic viral communities maintain high local diversity.

Ocean plankton. Structure and function of the global ocean microbiome.

Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

Ocean plankton. Environmental characteristics of Agulhas rings affect interocean plankton transport.

Agulhas rings provide the principal route for ocean waters to circulate from the Indo-Pacific to the Atlantic basin. Their influence on global ocean circulation is well known, but their role in plankton transport is largely unexplored. We show that, although the coarse taxonomic structure of plankton communities is continuous across the Agulhas choke point, South Atlantic plankton diversity is altered compared with Indian Ocean source populations. Modeling and in situ sampling of a young Agulhas ring indicate that strong vertical mixing drives complex nitrogen cycling, shaping community metabolism and biogeochemical signatures as the ring and associated plankton transit westward. The peculiar local environment inside Agulhas rings may provide a selective mechanism contributing to the limited dispersal of Indian Ocean plankton populations into the Atlantic.

Ocean plankton. Eukaryotic plankton diversity in the sunlit ocean.

Marine plankton support global biological and geochemical processes. Surveys of their biodiversity have hitherto been geographically restricted and have not accounted for the full range of plankton size. We assessed eukaryotic diversity from 334 size-fractionated photic-zone plankton communities collected across tropical and temperate oceans during the circumglobal Tara Oceans expedition. We analyzed 18S ribosomal DNA sequences across the intermediate plankton-size spectrum from the smallest unicellular eukaryotes (protists, >0.8 micrometers) to small animals of a few millimeters. Eukaryotic ribosomal diversity saturated at ~150,000 operational taxonomic units, about one-third of which could not be assigned to known eukaryotic groups. Diversity emerged at all taxonomic levels, both within the groups comprising the ~11,200 cataloged morphospecies of eukaryotic plankton and among twice as many other deep-branching lineages of unappreciated importance in plankton ecology studies. Most eukaryotic plankton biodiversity belonged to heterotrophic protistan groups, particularly those known to be parasites or symbiotic hosts.

Ocean plankton. Determinants of community structure in the global plankton interactome.

Species interaction networks are shaped by abiotic and biotic factors. Here, as part of the Tara Oceans project, we studied the photic zone interactome using environmental factors and organismal abundance profiles and found that environmental factors are incomplete predictors of community structure. We found associations across plankton functional types and phylogenetic groups to be nonrandomly distributed on the network and driven by both local and global patterns. We identified interactions among grazers, primary producers, viruses, and (mainly parasitic) symbionts and validated network-generated hypotheses using microscopy to confirm symbiotic relationships. We have thus provided a resource to support further research on ocean food webs and integrating biological components into ocean models.