Marine DNA Viral Macro- and Microdiversity from Pole to Pole.

Microbes drive most ecosystems and are modulated by viruses that impact their lifespan, gene flow, and metabolic outputs. However, ecosystem-level impacts of viral community diversity remain difficult to assess due to classification issues and few reference genomes. Here, we establish an ∼12-fold expanded global ocean DNA virome dataset of 195,728 viral populations, now including the Arctic Ocean, and validate that these populations form discrete genotypic clusters. Meta-community analyses revealed five ecological zones throughout the global ocean, including two distinct Arctic regions. Across the zones, local and global patterns and drivers in viral community diversity were established for both macrodiversity (inter-population diversity) and microdiversity (intra-population genetic variation). These patterns sometimes, but not always, paralleled those from macro-organisms and revealed temperate and tropical surface waters and the Arctic as biodiversity hotspots and mechanistic hypotheses to explain them. Such further understanding of ocean viruses is critical for broader inclusion in ecosystem models.

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.

Global Trends in Marine Plankton Diversity across Kingdoms of Life.

The ocean is home to myriad small planktonic organisms that underpin the functioning of marine ecosystems. However, their spatial patterns of diversity and the underlying drivers remain poorly known, precluding projections of their responses to global changes. Here we investigate the latitudinal gradients and global predictors of plankton diversity across archaea, bacteria, eukaryotes, and major virus clades using both molecular and imaging data from Tara Oceans. We show a decline of diversity for most planktonic groups toward the poles, mainly driven by decreasing ocean temperatures. Projections into the future suggest that severe warming of the surface ocean by the end of the 21st century could lead to tropicalization of the diversity of most planktonic groups in temperate and polar regions. These changes may have multiple consequences for marine ecosystem functioning and services and are expected to be particularly significant in key areas for carbon sequestration, fisheries, and marine conservation. 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.

Meta-Omics Reveals Genetic Flexibility of Diatom Nitrogen Transporters in Response to Environmental Changes.

Diatoms (Bacillariophyta), one of the most abundant and diverse groups of marine phytoplankton, respond rapidly to the supply of new nutrients, often out-competing other phytoplankton. Herein, we integrated analyses of the evolution, distribution, and expression modulation of two gene families involved in diatom nitrogen uptake (DiAMT1 and DiNRT2), in order to infer the main drivers of divergence in a key functional trait of phytoplankton. Our results suggest that major steps in the evolution of the two gene families reflected key events triggering diatom radiation and diversification. Their expression is modulated in the contemporary ocean by seawater temperature, nitrate, and iron concentrations. Moreover, the differences in diversity and expression of these gene families throughout the water column hint at a possible link with bacterial activity. This study represents a proof-of-concept of how a holistic approach may shed light on the functional biology of organisms in their natural environment.

Into the bloom: Molecular response of pelagic tunicates to fluctuating food availability.

The planktonic tunicates appendicularians and thaliaceans are highly efficient filter feeders on a wide range of prey size including bacteria and have shorter generation times than any other marine grazers. These traits allow some tunicate species to reach high population densities and ensure their success in a favourable environment. However, there are still few studies focusing on which genes and gene pathways are associated with responses of pelagic tunicates to environmental variability. Herein, we present the effect of food availability increase on tunicate community and gene expression at the Marquesas Islands (South-East Pacific Ocean). By using data from the Tara Oceans expedition, we show that changes in phytoplankton density and composition trigger the success of a dominant larvacean species (an undescribed appendicularian). Transcriptional signature to the autotroph bloom suggests key functions in specific physiological processes, i.e., energy metabolism, muscle contraction, membrane trafficking, and proteostasis. The relative abundance of reverse transcription-related Pfams was lower at bloom conditions, suggesting a link with adaptive genetic diversity in tunicates in natural ecosystems. Downstream of the bloom, pelagic tunicates were outcompeted by copepods. Our work represents the first metaomics study of the biological effects of phytoplankton bloom on a key zooplankton taxon.

Delineating ecologically significant taxonomic units from global patterns of marine picocyanobacteria.

Prochlorococcus and Synechococcus are the two most abundant and widespread phytoplankton in the global ocean. To better understand the factors controlling their biogeography, a reference database of the high-resolution taxonomic marker petB, encoding cytochrome b6, was used to recruit reads out of 109 metagenomes from the Tara Oceans expedition. An unsuspected novel genetic diversity was unveiled within both genera, even for the most abundant and well-characterized clades, and 136 divergent petB sequences were successfully assembled from metagenomic reads, significantly enriching the reference database. We then defined Ecologically Significant Taxonomic Units (ESTUs)-that is, organisms belonging to the same clade and occupying a common oceanic niche. Three major ESTU assemblages were identified along the cruise transect for Prochlorococcus and eight for Synechococcus Although Prochlorococcus HLIIIA and HLIVA ESTUs codominated in iron-depleted areas of the Pacific Ocean, CRD1 and the yet-to-be cultured EnvB were the prevalent Synechococcus clades in this area, with three different CRD1 and EnvB ESTUs occupying distinct ecological niches with regard to iron availability and temperature. Sharp community shifts were also observed over short geographic distances-for example, around the Marquesas Islands or between southern Indian and Atlantic Oceans-pointing to a tight correlation between ESTU assemblages and specific physico-chemical parameters. Together, this study demonstrates that there is a previously overlooked, ecologically meaningful, fine-scale diversity within some currently defined picocyanobacterial ecotypes, bringing novel insights into the ecology, diversity, and biology of the two most abundant phototrophs on Earth.

Insights into global diatom distribution and diversity in the world's ocean.

Diatoms (Bacillariophyta) constitute one of the most diverse and ecologically important groups of phytoplankton. They are considered to be particularly important in nutrient-rich coastal ecosystems and at high latitudes, but considerably less so in the oligotrophic open ocean. The Tara Oceans circumnavigation collected samples from a wide range of oceanic regions using a standardized sampling procedure. Here, a total of ∼12 million diatom V9-18S ribosomal DNA (rDNA) ribotypes, derived from 293 size-fractionated plankton communities collected at 46 sampling sites across the global ocean euphotic zone, have been analyzed to explore diatom global diversity and community composition. We provide a new estimate of diversity of marine planktonic diatoms at 4,748 operational taxonomic units (OTUs). Based on the total assigned ribotypes, Chaetoceros was the most abundant and diverse genus, followed by Fragilariopsis, Thalassiosira, and Corethron We found only a few cosmopolitan ribotypes displaying an even distribution across stations and high abundance, many of which could not be assigned with confidence to any known genus. Three distinct communities from South Pacific, Mediterranean, and Southern Ocean waters were identified that share a substantial percentage of ribotypes within them. Sudden drops in diversity were observed at Cape Agulhas, which separates the Indian and Atlantic Oceans, and across the Drake Passage between the Atlantic and Southern Oceans, indicating the importance of these ocean circulation choke points in constraining diatom distribution and diversity. We also observed high diatom diversity in the open ocean, suggesting that diatoms may be more relevant in these oceanic systems than generally considered.

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 dynamics of sexual phase in the marine diatom Pseudo-nitzschia multistriata (Bacillariophyceae).

Sexual reproduction represents a fundamental phase in the life cycle of diatoms, linked to both the production of genotypic diversity and the formation of large-sized initial cells. Only cells below a certain size threshold can be sexualized, but various environmental factors can modulate the success of sexual reproduction. We investigated the role of cell density and physiological conditions of parental strains in affecting the success and timing of sexual reproduction in the marine heterothallic diatom Pseudo-nitzschia multistriata. We also studied the dynamics of the sexual phase in still conditions allowing cell sedimentation and in gently mixed conditions that keep cells in suspension. Our results showed that successful sexual reproduction can only be achieved when crossing parental strains in the exponential growth phase. Evidence was provided for the fact that sexual reproduction is a density-dependent event and requires a threshold cell concentration to start, although this might vary considerably amongst strains. Moreover, the onset of the sexual phase was coupled to a marked reduction in growth of the vegetative parental cells. The crosses carried out in physically mixed conditions produced a significantly reduced number of sexual stages as compared to crosses in still conditions, showing that mixing impairs sexualization. The results of our experiments suggest that the signaling that triggers the sexual phase is favored when cells can accumulate, reducing the distance between them and facilitating contacts and/or the perception of chemical cues. Information on the progression of the sexual phase in laboratory conditions help understanding the conditions at which sex occurs in the natural environment.

Deviation from neutral species abundance distributions unveils geographical differences in the structure of diatom communities.

In recent years, the application of metagenomics techniques has advanced our understanding of plankton communities and their global distribution. Despite this progress, the relationship between the abundance distribution of diatom species and varying marine environmental conditions remains poorly understood. This study, leveraging data from the Tara Oceans expedition, tests the hypothesis that diatoms in sampled stations display a consistent species abundance distribution structure, as though they were sampled from a single ocean-wide metacommunity. Using a neutral sampling theory, we thus develop a framework to estimate the structure and diversity of diatom communities at each sampling station given the shape of the species abundance distribution of the metacommunity and the information of a reference station. Our analysis reveals a substantial temperature gradient in the discrepancies between predicted and observed biodiversity across the sampled stations. These findings challenge the hypothesis of a single neutral metacommunity, indicating that environmental differences substantially influence both the composition and structure of diatom communities.

Ocean-wide comparisons of mesopelagic planktonic community structures.

For decades, marine plankton have been investigated for their capacity to modulate biogeochemical cycles and provide fishery resources. Between the sunlit (epipelagic) layer and the deep dark waters, lies a vast and heterogeneous part of the ocean: the mesopelagic zone. How plankton composition is shaped by environment has been well-explored in the epipelagic but much less in the mesopelagic ocean. Here, we conducted comparative analyses of trans-kingdom community assemblages thriving in the mesopelagic oxygen minimum zone (OMZ), mesopelagic oxic, and their epipelagic counterparts. We identified nine distinct types of intermediate water masses that correlate with variation in mesopelagic community composition. Furthermore, oxygen, NO3- and particle flux together appeared as the main drivers governing these communities. Novel taxonomic signatures emerged from OMZ while a global co-occurrence network analysis showed that about 70% of the abundance of mesopelagic plankton groups is organized into three community modules. One module gathers prokaryotes, pico-eukaryotes and Nucleo-Cytoplasmic Large DNA Viruses (NCLDV) from oxic regions, and the two other modules are enriched in OMZ prokaryotes and OMZ pico-eukaryotes, respectively. We hypothesize that OMZ conditions led to a diversification of ecological niches, and thus communities, due to selective pressure from limited resources. Our study further clarifies the interplay between environmental factors in the mesopelagic oxic and OMZ, and the compositional features of communities.

Trade-off between sex and growth in diatoms: Molecular mechanisms and demographic implications.

Diatoms are fast-growing and winning competitors in aquatic environments, possibly due to optimized growth performance. However, their life cycles are complex, heteromorphic, and not fully understood. Here, we report on the fine control of cell growth and physiology during the sexual phase of the marine diatom Pseudo-nitzschia multistriata. We found that mating, under nutrient replete conditions, induces a prolonged growth arrest in parental cells. Transcriptomic analyses revealed down-regulation of genes related to major metabolic functions from the early phases of mating. Single-cell photophysiology also pinpointed an inhibition of photosynthesis and storage lipids accumulated in the arrested population, especially in gametes and zygotes. Numerical simulations revealed that growth arrest affects the balance between parental cells and their siblings, possibly favoring the new generation. Thus, in addition to resources availability, life cycle traits contribute to shaping the species ecological niches and must be considered to describe and understand the structure of plankton communities.

Cryptic and abundant marine viruses at the evolutionary origins of Earth's RNA virome.

Whereas DNA viruses are known to be abundant, diverse, and commonly key ecosystem players, RNA viruses are insufficiently studied outside disease settings. In this study, we analyzed ≈28 terabases of Global Ocean RNA sequences to expand Earth's RNA virus catalogs and their taxonomy, investigate their evolutionary origins, and assess their marine biogeography from pole to pole. Using new approaches to optimize discovery and classification, we identified RNA viruses that necessitate substantive revisions of taxonomy (doubling phyla and adding >50% new classes) and evolutionary understanding. "Species"-rank abundance determination revealed that viruses of the new phyla "Taraviricota," a missing link in early RNA virus evolution, and "Arctiviricota" are widespread and dominant in the oceans. These efforts provide foundational knowledge critical to integrating RNA viruses into ecological and epidemiological models.

Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems.

Biogeographical studies have traditionally focused on readily visible organisms, but recent technological advances are enabling analyses of the large-scale distribution of microscopic organisms, whose biogeographical patterns have long been debated. Here we assessed the global structure of plankton geography and its relation to the biological, chemical, and physical context of the ocean (the 'seascape') by analyzing metagenomes of plankton communities sampled across oceans during the Tara Oceans expedition, in light of environmental data and ocean current transport. Using a consistent approach across organismal sizes that provides unprecedented resolution to measure changes in genomic composition between communities, we report a pan-ocean, size-dependent plankton biogeography overlying regional heterogeneity. We found robust evidence for a basin-scale impact of transport by ocean currents on plankton biogeography, and on a characteristic timescale of community dynamics going beyond simple seasonality or life history transitions of plankton.

Oceans are brimming with life invisible to our eyes, a myriad of species of bacteria, viruses and other microscopic organisms essential for the health of the planet. These 'marine plankton' are unable to swim against currents and should therefore be constantly on the move, yet previous studies have suggested that distinct species of plankton may in fact inhabit different oceanic regions. However, proving this theory has been challenging; collecting plankton is logistically difficult, and it is often impossible to distinguish between species simply by examining them under a microscope. However, within the last decade, a research schooner called Tara has travelled the globe to gather thousands of plankton samples. At the same time, advances in genomics have made it possible to identify species based only on fragments of their DNA sequence. To understand the hidden geography of plankton communities in Earth's oceans, Richter et al. pored over DNA from the Tara Oceans expedition. This revealed that, despite being unable to resist the flow of water, various planktonic species which live close to the surface manage to occupy distinct, stable provinces shaped by currents. Different sizes of plankton are distributed in different sized provinces, with the smallest organisms tending to inhabit the smallest areas. Comparing DNA similarities and speeds of currents at the ocean surface revealed how these might stretch and mix plankton communities. Plankton play a critical role in the health of the ocean and the chemical cycles of planet Earth. These results could allow deeper investigation by marine modellers, ecologists, and evolutionary biologists. Meanwhile, work is already underway to investigate how climate change might impact this hidden geography.

Global drivers of eukaryotic plankton biogeography in the sunlit ocean.

Eukaryotic plankton are a core component of marine ecosystems with exceptional taxonomic and ecological diversity, yet how their ecology interacts with the environment to drive global distribution patterns is poorly understood. In this work, we use Tara Oceans metabarcoding data, which cover all major ocean basins, combined with a probabilistic model of taxon co-occurrence to compare the biogeography of 70 major groups of eukaryotic plankton. We uncover two main axes of biogeographic variation. First, more-diverse groups display clearer biogeographic patterns. Second, large-bodied consumers are structured by oceanic basins, mostly through the main current systems, whereas small-bodied phototrophs are structured by latitude and follow local environmental conditions. Our study highlights notable differences in biogeographies across plankton groups and investigates their determinants at the global scale.

Environmental vulnerability of the global ocean epipelagic plankton community interactome.

Marine plankton form complex communities of interacting organisms at the base of the food web, which sustain oceanic biogeochemical cycles and help regulate climate. Although global surveys are starting to reveal ecological drivers underlying planktonic community structure and predicted climate change responses, it is unclear how community-scale species interactions will be affected by climate change. Here, we leveraged Tara Oceans sampling to infer a global ocean cross-domain plankton co-occurrence network-the community interactome-and used niche modeling to assess its vulnerabilities to environmental change. Globally, this revealed a plankton interactome self-organized latitudinally into marine biomes (Trades, Westerlies, Polar) and more connected poleward. Integrated niche modeling revealed biome-specific community interactome responses to environmental change and forecasted the most affected lineages for each community. These results provide baseline approaches to assess community structure and organismal interactions under climate scenarios while identifying plausible plankton bioindicators for ocean monitoring of climate change.

Tara Oceans: towards global ocean ecosystems biology.

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

The Water Mass Transformation Framework for Ocean Physics and Biogeochemistry.

The water mass transformation (WMT) framework weaves together circulation, thermodynamics, and biogeochemistry into a description of the ocean that complements traditional Eulerian and Lagrangian methods. In so doing, a WMT analysis renders novel insights and predictive capabilities for studies of ocean physics and biogeochemistry. In this review, we describe fundamentals of the WMT framework and illustrate its practical analysis capabilities. We show how it provides a robust methodology to characterize and quantify the impact of physical processes on buoyancy and other thermodynamic fields. We also detail how to extend WMT to insightful analysis of biogeochemical cycles.

Discovering millions of plankton genomic markers from the Atlantic Ocean and the Mediterranean Sea.

Comparison of the molecular diversity in all plankton populations present in geographically distant water columns may allow for a holistic view of the connectivity, isolation and adaptation of organisms in the marine environment. In this context, a large-scale detection and analysis of genomic variants directly in metagenomic data appeared as a powerful strategy for the identification of genetic structures and genes under natural selection in plankton. Here, we used discosnp++, a reference-free variant caller, to produce genetic variants from large-scale metagenomic data and assessed its accuracy on the copepod Oithona nana in terms of variant calling, allele frequency estimation and population genomic statistics by comparing it to the state-of-the-art method. discosnp ++ produces variants leading to similar conclusions regarding the genetic structure and identification of loci under natural selection. discosnp++ was then applied to 120 metagenomic samples from four size fractions, including prokaryotes, protists and zooplankton sampled from 39 tara Oceans sampling stations located in the Atlantic Ocean and the Mediterranean Sea to produce a new set of marine genomic markers containing more than 19 million of variants. This new genomic resource can be used by the community to relocate these markers on their plankton genomes or transcriptomes of interest. This resource will be updated with new marine expeditions and the increase of metagenomic data (availability: http://bioinformatique.rennes.inria.fr/taravariants/).