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.

Morphological diversity increases with decreasing resources along a zooplankton time series.

Biodiversity is studied notably because of its reciprocal relationship with ecosystem functions such as production. Diversity is traditionally described from a taxonomic, genetic or functional point of view but the diversity in organism morphology is seldom explicitly considered, except for body size. We describe morphological diversity of marine zooplankton seasonally and over 12 years using quantitative imaging of weekly plankton samples, in the northwestern Mediterranean Sea. We extract 45 morphological features on greater than 800 000 individuals, which we summarize into four main morphological traits (size, transparency, circularity and shape complexity). In this morphological space, we define objective morphological groups and, from those, compute morphological diversity indices (richness, evenness and divergence) using metrics originally defined for functional diversity. On both time scales, morphological diversity increased when nutritive resources and plankton concentrations were low, thus matching the theoretical reciprocal relationship. Over the long term at least, this diversity increase was not fully attributable to taxonomic diversity changes. The decline in the most common plankton forms and the increase in morphological variance and in extreme morphologies suggest a mechanism akin to specialization under low production, with likely consequences for trophic structure and carbon flux.

Genomic differentiation of three pico-phytoplankton species in the Mediterranean Sea.

For more than a decade, high-throughput sequencing has transformed the study of marine planktonic communities and has highlighted the extent of protist diversity in these ecosystems. Nevertheless, little is known relative to their genomic diversity at the species-scale as well as their major speciation mechanisms. An increasing number of data obtained from global scale sampling campaigns is becoming publicly available, and we postulate that metagenomic data could contribute to deciphering the processes shaping protist genomic differentiation in the marine realm. As a proof of concept, we developed a findable, accessible, interoperable and reusable (FAIR) pipeline and focused on the Mediterranean Sea to study three a priori abundant protist species: Bathycoccus prasinos, Pelagomonas calceolata and Phaeocystis cordata. We compared the genomic differentiation of each species in light of geographic, environmental and oceanographic distances. We highlighted that isolation-by-environment shapes the genomic differentiation of B. prasinos, whereas P. cordata is impacted by geographic distance (i.e. isolation-by-distance). At present time, the use of metagenomics to accurately estimate the genomic differentiation of protists remains challenging since coverages are lower compared to traditional population surveys. However, our approach sheds light on ecological and evolutionary processes occurring within natural marine populations and paves the way for future protist population metagenomic studies.

Length, width, shape regularity, and chain structure: time series analysis of phytoplankton morphology from imagery.

Functional traits are increasingly used to assess changes in phytoplankton community structure and to link individual characteristics to ecosystem functioning. However, they are usually inferred from taxonomic identification or manually measured for each organism, both time consuming approaches. Instead, we focus on high throughput imaging to describe the main temporal variations of morphological changes of phytoplankton in Narragansett Bay, a coastal time-series station. We analyzed a 2-yr dataset of morphological features automatically extracted from continuous imaging of individual phytoplankton images (~ 105 million images collected by an Imaging FlowCytobot). We identified synthetic morphological traits using multivariate analysis and revealed that morphological variations were mainly due to changes in length, width, shape regularity, and chain structure. Morphological changes were especially important in winter with successive peaks of larger cells with increasing complexity and chains more clearly connected. Small nanophytoplankton were present year-round and constituted the base of the community, especially apparent during the transitions between diatom blooms. High inter-annual variability was also observed. On a weekly timescale, increases in light were associated with more clearly connected chains while more complex shapes occurred at lower nitrogen concentrations. On an hourly timescale, temperature was the determinant variable constraining cell morphology, with a general negative influence on length and a positive one on width, shape regularity, and chain structure. These first insights into the phytoplankton morphology of Narragansett Bay highlight the possible morphological traits driving the phytoplankton succession in response to light, temperature, and nutrient changes.

In situ imaging reveals the biomass of giant protists in the global ocean.

Planktonic organisms play crucial roles in oceanic food webs and global biogeochemical cycles. Most of our knowledge about the ecological impact of large zooplankton stems from research on abundant and robust crustaceans, and in particular copepods. A number of the other organisms that comprise planktonic communities are fragile, and therefore hard to sample and quantify, meaning that their abundances and effects on oceanic ecosystems are poorly understood. Here, using data from a worldwide in situ imaging survey of plankton larger than 600 µm, we show that a substantial part of the biomass of this size fraction consists of giant protists belonging to the Rhizaria, a super-group of mostly fragile unicellular marine organisms that includes the taxa Phaeodaria and Radiolaria (for example, orders Collodaria and Acantharia). Globally, we estimate that rhizarians in the top 200 m of world oceans represent a standing stock of 0.089 Pg carbon, equivalent to 5.2% of the total oceanic biota carbon reservoir. In the vast oligotrophic intertropical open oceans, rhizarian biomass is estimated to be equivalent to that of all other mesozooplankton (plankton in the size range 0.2-20 mm). The photosymbiotic association of many rhizarians with microalgae may be an important factor in explaining their distribution. The previously overlooked importance of these giant protists across the widest ecosystem on the planet changes our understanding of marine planktonic ecosystems.

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.

The Underwater Vision Profiler 6: an imaging sensor of particle size spectra and plankton, for autonomous and cabled platforms.

Autonomous and cabled platforms are revolutionizing our understanding of ocean systems by providing 4D monitoring of the water column, thus going beyond the reach of ship-based surveys and increasing the depth of remotely sensed observations. However, very few commercially available sensors for such platforms are capable of monitoring large particulate matter (100-2000 µm) and plankton despite their important roles in the biological carbon pump and as trophic links from phytoplankton to fish. Here, we provide details of a new, commercially available scientific camera-based particle counter, specifically designed to be deployed on autonomous and cabled platforms: the Underwater Vision Profiler 6 (UVP6). Indeed, the UVP6 camera-and-lighting and processing system, while small in size and requiring low power, provides data of quality comparable to that of previous much larger UVPs deployed from ships. We detail the UVP6 camera settings, its performance when acquiring data on aquatic particles and plankton, their quality control, analysis of its recordings, and streaming from in situ acquisition to users. In addition, we explain how the UVP6 has already been integrated into platforms such as BGC-Argo floats, gliders and long-term mooring systems (autonomous platforms). Finally, we use results from actual deployments to illustrate how UVP6 data can contribute to addressing longstanding questions in marine science, and also suggest new avenues that can be explored using UVP6-equipped autonomous platforms.

Light color acclimation is a key process in the global ocean distribution of Synechococcus cyanobacteria.

Marine Synechococcus cyanobacteria are major contributors to global oceanic primary production and exhibit a unique diversity of photosynthetic pigments, allowing them to exploit a wide range of light niches. However, the relationship between pigment content and niche partitioning has remained largely undetermined due to the lack of a single-genetic marker resolving all pigment types (PTs). Here, we developed and employed a robust method based on three distinct marker genes (cpcBA, mpeBA, and mpeW) to estimate the relative abundance of all known Synechococcus PTs from metagenomes. Analysis of the Tara Oceans dataset allowed us to reveal the global distribution of Synechococcus PTs and to define their environmental niches. Green-light specialists (PT 3a) dominated in warm, green equatorial waters, whereas blue-light specialists (PT 3c) were particularly abundant in oligotrophic areas. Type IV chromatic acclimaters (CA4-A/B), which are able to dynamically modify their light absorption properties to maximally absorb green or blue light, were unexpectedly the most abundant PT in our dataset and predominated at depth and high latitudes. We also identified populations in which CA4 might be nonfunctional due to the lack of specific CA4 genes, notably in warm high-nutrient low-chlorophyll areas. Major ecotypes within clades I-IV and CRD1 were preferentially associated with a particular PT, while others exhibited a wide range of PTs. Altogether, this study provides important insights into the ecology of Synechococcus and highlights the complex interactions between vertical phylogeny, pigmentation, and environmental parameters that shape Synechococcus community structure and evolution.

Statistical distributions of trace metal concentrations in the northwestern Mediterranean atmospheric aerosol.

The concentrations of 11 crustal and anthropogenic trace metals (Li, Al, V, Mn, Fe, Co, Ni, Cu, Zn, Cd, Pb) were measured from 2006 to 2008 in the atmospheric aerosol at a northwestern Mediterranean coast (station of Cap Ferrat, situated on the southeastern coast of France). Statistical models (lognormal, Weibull, and gamma) that best represented the trace metal distribution for this environment are described. The lognormal model was selected for the distributions of (in decreasing strength of the fit) Al, Co, Li, Zn, Mn, Cu, Pb, and Cd, i.e., metals that are introduced into the atmospheric aerosol by pulses inducing temporal variability in their concentrations. The gamma model was associated with Fe, i.e., metals that exhibit less inter-annual variability than the former trace metals. The third mode (Weibull) represented the distribution of the concentrations of V and Ni. The statistical approach presented in this study contributed to better define and constrain the distribution of the 11 trace metals of the atmospheric aerosol from the northwestern Mediterranean coast. In a close future, knowledge of these statistical distributions will allow using convolution models to separate their natural and anthropogenic contributions, therefore increasing our ability to study anthropogenic emissions of trace metals and their impact on the environment.

A large population of king crabs in Palmer Deep on the west Antarctic Peninsula shelf and potential invasive impacts.

Lithodid crabs (and other skeleton-crushing predators) may have been excluded from cold Antarctic continental shelf waters for more than 14 Myr. The west Antarctic Peninsula shelf is warming rapidly and has been hypothesized to be soon invaded by lithodids. A remotely operated vehicle survey in Palmer Deep, a basin 120 km onto the Antarctic shelf, revealed a large, reproductive population of lithodids, providing the first evidence that king crabs have crossed the Antarctic shelf. DNA sequencing and morphology indicate the lithodid is Neolithodes yaldwyni Ahyong & Dawson, previously reported only from Ross Sea waters. We estimate a N. yaldwyni population density of 10 600 km(-2) and a population size of 1.55 × 10(6) in Palmer Deep, a density similar to lithodid populations of commercial interest around Alaska and South Georgia. The lithodid occurred at depths of more than 850 m and temperatures of more than 1.4°C in Palmer Deep, and was not found in extensive surveys of the colder shelf at depths of 430-725 m. Where N. yaldwyni occurred, crab traces were abundant, megafaunal diversity reduced and echinoderms absent, suggesting that the crabs have major ecological impacts. Antarctic Peninsula shelf waters are warming at approximately 0.01°C yr(-1); if N. yaldwyni is currently limited by cold temperatures, it could spread up onto the shelf (400-600 m depths) within 1-2 decades. The Palmer Deep N. yaldwyni population provides an important model for the potential invasive impacts of crushing predators on vulnerable Antarctic shelf ecosystems.

Diversity and ecological footprint of Global Ocean RNA viruses.

DNA viruses are increasingly recognized as influencing marine microbes and microbe-mediated biogeochemical cycling. However, little is known about global marine RNA virus diversity, ecology, and ecosystem roles. In this study, we uncover patterns and predictors of marine RNA virus community- and "species"-level diversity and contextualize their ecological impacts from pole to pole. Our analyses revealed four ecological zones, latitudinal and depth diversity patterns, and environmental correlates for RNA viruses. Our findings only partially parallel those of cosampled plankton and show unexpectedly high polar ecological interactions. The influence of RNA viruses on ecosystems appears to be large, as predicted hosts are ecologically important. Moreover, the occurrence of auxiliary metabolic genes indicates that RNA viruses cause reprogramming of diverse host metabolisms, including photosynthesis and carbon cycling, and that RNA virus abundances predict ocean carbon export.

Patterns of eukaryotic diversity from the surface to the deep-ocean sediment.

Remote deep-ocean sediment (DOS) ecosystems are among the least explored biomes on Earth. Genomic assessments of their biodiversity have failed to separate indigenous benthic organisms from sinking plankton. Here, we compare global-scale eukaryotic DNA metabarcoding datasets (18S-V9) from abyssal and lower bathyal surficial sediments and euphotic and aphotic ocean pelagic layers to distinguish plankton from benthic diversity in sediment material. Based on 1685 samples collected throughout the world ocean, we show that DOS diversity is at least threefold that in pelagic realms, with nearly two-thirds represented by abundant yet unknown eukaryotes. These benthic communities are spatially structured by ocean basins and particulate organic carbon (POC) flux from the upper ocean. Plankton DNA reaching the DOS originates from abundant species, with maximal deposition at high latitudes. Its seafloor DNA signature predicts variations in POC export from the surface and reveals previously overlooked taxa that may drive the biological carbon pump.

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.

Eukaryotic virus composition can predict the efficiency of carbon export in the global ocean.

The biological carbon pump, in which carbon fixed by photosynthesis is exported to the deep ocean through sinking, is a major process in Earth's carbon cycle. The proportion of primary production that is exported is termed the carbon export efficiency (CEE). Based on in-lab or regional scale observations, viruses were previously suggested to affect the CEE (i.e., viral "shunt" and "shuttle"). In this study, we tested associations between viral community composition and CEE measured at a global scale. A regression model based on relative abundance of viral marker genes explained 67% of the variation in CEE. Viruses with high importance in the model were predicted to infect ecologically important hosts. These results are consistent with the view that the viral shunt and shuttle functions at a large scale and further imply that viruses likely act in this process in a way dependent on their hosts and ecosystem dynamics.

Macroscale patterns of oceanic zooplankton composition and size structure.

Ocean plankton comprise organisms from viruses to fish larvae that are fundamental to ecosystem functioning and the provision of marine services such as fisheries and CO2 sequestration. The latter services are partly governed by variations in plankton community composition and the expression of traits such as body size at community-level. While community assembly has been thoroughly studied for the smaller end of the plankton size spectrum, the larger end comprises ectotherms that are often studied at the species, or group-level, rather than as communities. The body size of marine ectotherms decreases with temperature, but controls on community-level traits remain elusive, hindering the predictability of marine services provision. Here, we leverage Tara Oceans datasets to determine how zooplankton community composition and size structure varies with latitude, temperature and productivity-related covariates in the global surface ocean. Zooplankton abundance and median size decreased towards warmer and less productive environments, as a result of changes in copepod composition. However, some clades displayed the opposite relationships, which may be ascribed to alternative feeding strategies. Given that climate models predict increasingly warmed and stratified oceans, our findings suggest that zooplankton communities will shift towards smaller organisms which might weaken their contribution to the biological carbon pump.

Probabilistic modeling to estimate jellyfish ecophysiological properties and size distributions.

While Ocean modeling has made significant advances over the last decade, its complex biological component is still oversimplified. In particular, modeling organisms in the ocean system must integrate parameters to fit both physiological and ecological behaviors that are together very difficult to determine. Such difficulty occurs for modeling Pelagia noctiluca. This jellyfish has a high abundance in the Mediterranean Sea and could contribute to several biogeochemical processes. However, gelatinous zooplanktons remain poorly represented in biogeochemical models because uncertainties about their ecophysiology limit our understanding of their potential role and impact. To overcome this issue, we propose, for the first time, the use of the Statistical Model Checking Engine (SMCE), a probability-based computational framework that considers a set of parameters as a whole. Contrary to standard parameter inference techniques, SMCE identifies sets of parameters that fit both laboratory-culturing observations and in situ patterns while considering uncertainties. Doing so, we estimated the best parameter sets of the ecophysiological model that represents the jellyfish growth and degrowth in laboratory conditions as well as its size. Behind this application, SMCE remains a computational framework that supports the projection of a model with uncertainties in broader contexts such as biogeochemical processes to drive future studies.

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.

A global ocean atlas of eukaryotic genes.

While our knowledge about the roles of microbes and viruses in the ocean has increased tremendously due to recent advances in genomics and metagenomics, research on marine microbial eukaryotes and zooplankton has benefited much less from these new technologies because of their larger genomes, their enormous diversity, and largely unexplored physiologies. Here, we use a metatranscriptomics approach to capture expressed genes in open ocean Tara Oceans stations across four organismal size fractions. The individual sequence reads cluster into 116 million unigenes representing the largest reference collection of eukaryotic transcripts from any single biome. The catalog is used to unveil functions expressed by eukaryotic marine plankton, and to assess their functional biogeography. Almost half of the sequences have no similarity with known proteins, and a great number belong to new gene families with a restricted distribution in the ocean. Overall, the resource provides the foundations for exploring the roles of marine eukaryotes in ocean ecology and biogeochemistry.