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.

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.

Machine learning techniques to characterize functional traits of plankton from image data.

Plankton imaging systems supported by automated classification and analysis have improved ecologists' ability to observe aquatic ecosystems. Today, we are on the cusp of reliably tracking plankton populations with a suite of lab-based and in situ tools, collecting imaging data at unprecedentedly fine spatial and temporal scales. But these data have potential well beyond examining the abundances of different taxa; the individual images themselves contain a wealth of information on functional traits. Here, we outline traits that could be measured from image data, suggest machine learning and computer vision approaches to extract functional trait information from the images, and discuss promising avenues for novel studies. The approaches we discuss are data agnostic and are broadly applicable to imagery of other aquatic or terrestrial organisms.

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.

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.

New insights into global biogeography, population structure and natural selection from the genome of the epipelagic copepod Oithona.

In the epipelagic ocean, the genus Oithona is considered as one of the most abundant and widespread copepods and plays an important role in the trophic food web. Despite its ecological importance, little is known about Oithona and cyclopoid copepods genomics. Therefore, we sequenced, assembled and annotated the genome of Oithona nana. The comparative genomic analysis integrating available copepod genomes highlighted the expansions of genes related to stress response, cell differentiation and development, including genes coding Lin12-Notch-repeat (LNR) domain proteins. The Oithona biogeography based on 28S sequences and metagenomic reads from the Tara Oceans expedition showed the presence of O. nana mostly in the Mediterranean Sea (MS) and confirmed the amphitropical distribution of Oithona similis. The population genomics analyses of O. nana in the Northern MS, integrating the Tara Oceans metagenomic data and the O. nana genome, led to the identification of genetic structure between populations from the MS basins. Furthermore, 20 loci were found to be under positive selection including four missense and eight synonymous variants, harbouring soft or hard selective sweep patterns. One of the missense variants was localized in the LNR domain of the coding region of a male-specific gene. The variation in the B-allele frequency with respect to the MS circulation pattern showed the presence of genomic clines between O. nana and another undefined Oithona species possibly imported through Atlantic waters. This study provides new approaches and results in zooplankton population genomics through the integration of metagenomic and oceanographic data.

Global census of the significance of giant mesopelagic protists to the marine carbon and silicon cycles.

Thriving in both epipelagic and mesopelagic layers, Rhizaria are biomineralizing protists, mixotrophs or flux-feeders, often reaching gigantic sizes. In situ imaging showed their contribution to oceanic carbon stock, but left their contribution to element cycling unquantified. Here, we compile a global dataset of 167,551 Underwater Vision Profiler 5 Rhizaria images, and apply machine learning models to predict their organic carbon and biogenic silica biomasses in the uppermost 1000 m. We estimate that Rhizaria represent up to 1.7% of mesozooplankton carbon biomass in the top 500 m. Rhizaria biomass, dominated by Phaeodaria, is more than twice as high in the mesopelagic than in the epipelagic layer. Globally, the carbon demand of mesopelagic, flux-feeding Phaeodaria reaches 0.46 Pg C y-1, representing 3.8 to 9.2% of gravitational carbon export. Furthermore, we show that Rhizaria are a unique source of biogenic silica production in the mesopelagic layer, where no other silicifiers are present. Our global census further highlights the importance of Rhizaria for ocean biogeochemistry.

Monitoring and modelling marine zooplankton in a changing climate.

Zooplankton are major consumers of phytoplankton primary production in marine ecosystems. As such, they represent a critical link for energy and matter transfer between phytoplankton and bacterioplankton to higher trophic levels and play an important role in global biogeochemical cycles. In this Review, we discuss key responses of zooplankton to ocean warming, including shifts in phenology, range, and body size, and assess the implications to the biological carbon pump and interactions with higher trophic levels. Our synthesis highlights key knowledge gaps and geographic gaps in monitoring coverage that need to be urgently addressed. We also discuss an integrated sampling approach that combines traditional and novel techniques to improve zooplankton observation for the benefit of monitoring zooplankton populations and modelling future scenarios under global changes.

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.

Machine Learning for the Study of Plankton and Marine Snow from Images.

Quantitative imaging instruments produce a large number of images of plankton and marine snow, acquired in a controlled manner, from which the visual characteristics of individual objects and their in situ concentrations can be computed. To exploit this wealth of information, machine learning is necessary to automate tasks such as taxonomic classification. Through a review of the literature, we highlight the progress of those machine classifiers and what they can and still cannot be trusted for. Several examples showcase how the combination of quantitative imaging with machine learning has brought insights on pelagic ecology. They also highlight what is still missing and how images could be exploited further through trait-based approaches. In the future, we suggest deeper interactions with the computer sciences community, the adoption of data standards, and the more systematic sharing of databases to build a global community of pelagic image providers and users.

Marine snow morphology illuminates the evolution of phytoplankton blooms and determines their subsequent vertical export.

The organic carbon produced in the ocean's surface by phytoplankton is either passed through the food web or exported to the ocean interior as marine snow. The rate and efficiency of such vertical export strongly depend on the size, structure and shape of individual particles, but apart from size, other morphological properties are still not quantitatively monitored. With the growing number of in situ imaging technologies, there is now a great possibility to analyze the morphology of individual marine snow. Thus, automated methods for their classification are urgently needed. Consequently, here we present a simple, objective categorization method of marine snow into a few ecologically meaningful functional morphotypes using field data from successive phases of the Arctic phytoplankton bloom. The proposed approach is a promising tool for future studies aiming to integrate the diversity, composition and morphology of marine snow into our understanding of the biological carbon pump.

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.

Biogeography and diversity of Collodaria (Radiolaria) in the global ocean.

Collodaria are heterotrophic marine protists that exist either as large colonies composed of hundreds of cells or as large solitary cells. All described species so far harbour intracellular microalgae as photosymbionts. Although recent environmental diversity surveys based on molecular methods demonstrated their consistently high contribution to planktonic communities and their worldwide occurrence, our understanding of their diversity and biogeography is still very limited. Here we estimated the 18S ribosomal DNA (rDNA) gene copies per collodarian cell for solitary (5770±1960 small subunit (SSU) rDNA copies) and colonial specimens (37 474±17 799 SSU rDNA copies, for each individual cell within a colony) using single-specimen quantitative PCR. We then investigated the environmental diversity of Collodaria within the photic zone through the metabarcoding survey from the Tara Oceans expedition and found that the two collodarian families Collosphaeridae and Sphaerozoidae contributed the most to the collodarian diversity and encompassed mostly cosmopolitan taxa. Although the biogeographical patterns were homogeneous within each biogeochemical biome considered, we observed that coastal biomes were consistently less diverse than oceanic biomes and were dominated by the Sphaerozoidae while the Collosphaeridae were dominant in the open oceans. The significant relationships with six environmental variables suggest that collodarian diversity is influenced by the trophic status of oceanic provinces and increased towards more oligotrophic regions.

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.

Temporal Variability of Zooplankton (2000-2013) in the Levantine Sea: Significant Changes Associated to the 2005-2010 EMT-like Event?

In this study, we investigated, for the first time, the potential impact of environmental changes on zooplankton abundance over a fourteen year period (2000-2013) at an offshore station in the Eastern Mediterranean Sea (the Levantine basin, offshore Lebanon). Samples were collected monthly and analyzed using the semi-automated system ZooScan. Salinity, temperature and phytoplankton abundance (nano and microphytoplankton) were also measured. Results show no significant temporal trend in sea surface temperature over the years. Between 2005-2010, salinity in the upper layer (0-80 m) of the Levantine basin increased (~0.3°C). During this 5 year period, total zooplankton abundance significantly increased. These modifications were concomitant to the activation of Aegean Sea as a source of dense water formation as part of the "Eastern Mediterranean Transient-like" event. The results of the present study suggested that zooplankton benefited from enhanced phytoplankton production during the mixing years of the event. Changes in the phenology of some taxa were observed accordingly with a predominantly advanced peak of zooplankton abundance. In conclusion, long-term changes in zooplankton abundance were related to the Levantine thermohaline circulation rather than sea surface warming. Sampling must be maintained to assess the impact of long-term climate change on zooplankton communities.

Nitrogen Fuelling of the Pelagic Food Web of the Tropical Atlantic.

We estimated the relative contribution of atmosphere (ic Nitrogen (N) input (wet and dry deposition and N fixation) to the epipelagic food web by measuring N isotopes of different functional groups of epipelagic zooplankton along 23°W (17°N-4°S) and 18°N (20-24°W) in the Eastern Tropical Atlantic. Results were related to water column observations of nutrient distribution and vertical diffusive flux as well as colony abundance of Trichodesmium obtained with an Underwater Vision Profiler (UVP5). The thickness and depth of the nitracline and phosphocline proved to be significant predictors of zooplankton stable N isotope values. Atmospheric N input was highest (61% of total N) in the strongly stratified and oligotrophic region between 3 and 7°N, which featured very high depth-integrated Trichodesmium abundance (up to 9.4×10(4) colonies m(-2)), strong thermohaline stratification and low zooplankton δ15N (~2‰). Relative atmospheric N input was lowest south of the equatorial upwelling between 3 and 5°S (27%). Values in the Guinea Dome region and north of Cape Verde ranged between 45 and 50%, respectively. The microstructure-derived estimate of the vertical diffusive N flux in the equatorial region was about one order of magnitude higher than in any other area (approximately 8 mmol m(-2) d(1)). At the same time, this region received considerable atmospheric N input (35% of total). In general, zooplankton δ15N and Trichodesmium abundance were closely correlated, indicating that N fixation is the major source of atmospheric N input. Although Trichodesmium is not the only N fixing organism, its abundance can be used with high confidence to estimate the relative atmospheric N input in the tropical Atlantic (r2 = 0.95). Estimates of absolute N fixation rates are two- to tenfold higher than incubation-derived rates reported for the same regions. Our approach integrates over large spatial and temporal scales and also quantifies fixed N released as dissolved inorganic and organic N. In a global analysis, it may thus help to close the gap in oceanic N budgets.

Comprehensive model of annual plankton succession based on the whole-plankton time series approach.

Ecological succession provides a widely accepted description of seasonal changes in phytoplankton and mesozooplankton assemblages in the natural environment, but concurrent changes in smaller (i.e. microbes) and larger (i.e. macroplankton) organisms are not included in the model because plankton ranging from bacteria to jellies are seldom sampled and analyzed simultaneously. Here we studied, for the first time in the aquatic literature, the succession of marine plankton in the whole-plankton assemblage that spanned 5 orders of magnitude in size from microbes to macroplankton predators (not including fish or fish larvae, for which no consistent data were available). Samples were collected in the northwestern Mediterranean Sea (Bay of Villefranche) weekly during 10 months. Simultaneously collected samples were analyzed by flow cytometry, inverse microscopy, FlowCam, and ZooScan. The whole-plankton assemblage underwent sharp reorganizations that corresponded to bottom-up events of vertical mixing in the water-column, and its development was top-down controlled by large gelatinous filter feeders and predators. Based on the results provided by our novel whole-plankton assemblage approach, we propose a new comprehensive conceptual model of the annual plankton succession (i.e. whole plankton model) characterized by both stepwise stacking of four broad trophic communities from early spring through summer, which is a new concept, and progressive replacement of ecological plankton categories within the different trophic communities, as recognised traditionally.