The global genetic diversity of planktonic foraminifera reveals the structure of cryptic speciation in plankton.

The nature and extent of diversity in the plankton has fascinated scientists for over a century. Initially, the discovery of many new species in the remarkably uniform and unstructured pelagic environment appeared to challenge the concept of ecological niches. Later, it became obvious that only a fraction of plankton diversity had been formally described, because plankton assemblages are dominated by understudied eukaryotic lineages with small size that lack clearly distinguishable morphological features. The high diversity of the plankton has been confirmed by comprehensive metabarcoding surveys, but interpretation of the underlying molecular taxonomies is hindered by insufficient integration of genetic diversity with morphological taxonomy and ecological observations. Here we use planktonic foraminifera as a study model and reveal the full extent of their genetic diversity and investigate geographical and ecological patterns in their distribution. To this end, we assembled a global data set of ~7600 ribosomal DNA sequences obtained from morphologically characterised individual foraminifera, established a robust molecular taxonomic framework for the observed diversity, and used it to query a global metabarcoding data set covering ~1700 samples with ~2.48 billion reads. This allowed us to extract and assign 1 million reads, enabling characterisation of the structure of the genetic diversity of the group across ~1100 oceanic stations worldwide. Our sampling revealed the existence of, at most, 94 distinct molecular operational taxonomic units (MOTUs) at a level of divergence indicative of biological species. The genetic diversity only doubles the number of formally described species identified by morphological features. Furthermore, we observed that the allocation of genetic diversity to morphospecies is uneven. Only 16 morphospecies disguise evolutionarily significant genetic diversity, and the proportion of morphospecies that show genetic diversity increases poleward. Finally, we observe that MOTUs have a narrower geographic distribution than morphospecies and that in some cases the MOTUs belonging to the same morphospecies (cryptic species) have different environmental preferences. Overall, our analysis reveals that even in the light of global genetic sampling, planktonic foraminifera diversity is modest and finite. However, the extent and structure of the cryptic diversity reveals that genetic diversification is decoupled from morphological diversification, hinting at different mechanisms acting at different levels of divergence.

The FORCIS database: A global census of planktonic Foraminifera from ocean waters.

Planktonic Foraminifera are unique paleo-environmental indicators through their excellent fossil record in ocean sediments. Their distribution and diversity are affected by different environmental factors including anthropogenically forced ocean and climate change. Until now, historical changes in their distribution have not been fully assessed at the global scale. Here we present the FORCIS (Foraminifera Response to Climatic Stress) database on foraminiferal species diversity and distribution in the global ocean from 1910 until 2018 including published and unpublished data. The FORCIS database includes data collected using plankton tows, continuous plankton recorder, sediment traps and plankton pump, and contains ~22,000, ~157,000, ~9,000, ~400 subsamples, respectively (one single plankton aliquot collected within a depth range, time interval, size fraction range, at a single location) from each category. Our database provides a perspective of the distribution patterns of planktonic Foraminifera in the global ocean on large spatial (regional to basin scale, and at the vertical scale), and temporal (seasonal to interdecadal) scales over the past century.

Contamination of planktonic food webs in the Mediterranean Sea: Setting the frame for the MERITE-HIPPOCAMPE oceanographic cruise (spring 2019).

This paper looks at experiential feedback and the technical and scientific challenges tied to the MERITE-HIPPOCAMPE cruise that took place in the Mediterranean Sea in spring 2019. This cruise proposes an innovative approach to investigate the accumulation and transfer of inorganic and organic contaminants within the planktonic food webs. We present detailed information on how the cruise worked, including 1) the cruise track and sampling stations, 2) the overall strategy, based mainly on the collection of plankton, suspended particles and water at the deep chlorophyll maximum, and the separation of these particles and planktonic organisms into various size fractions, as well as the collection of atmospheric deposition, 3) the operations performed and material used at each station, and 4) the sequence of operations and main parameters analysed. The paper also provides the main environmental conditions that were prevailing during the campaign. Lastly, we present the types of articles produced based on work completed by the cruise that are part of this special issue.

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.

Genetic and morphological divergence in the warm-water planktonic foraminifera genus Globigerinoides.

The planktonic foraminifera genus Globigerinoides provides a prime example of a species-rich genus in which genetic and morphological divergence are uncorrelated. To shed light on the evolutionary processes that lead to the present-day diversity of Globigerinoides, we investigated the genetic, ecological and morphological divergence of its constituent species. We assembled a global collection of single-cell barcode sequences and show that the genus consists of eight distinct genetic types organized in five extant morphospecies. Based on morphological evidence, we reassign the species Globoturborotalita tenella to Globigerinoides and amend Globigerinoides ruber by formally proposing two new subspecies, G. ruber albus n.subsp. and G. ruber ruber in order to express their subspecies level distinction and to replace the informal G. ruber "white" and G. ruber "pink", respectively. The genetic types within G. ruber and Globigerinoides elongatus show a combination of endemism and coexistence, with little evidence for ecological differentiation. CT-scanning and ontogeny analysis reveal that the diagnostic differences in adult morphologies could be explained by alterations of the ontogenetic trajectories towards final (reproductive) size. This indicates that heterochrony may have caused the observed decoupling between genetic and morphological diversification within the genus. We find little evidence for environmental forcing of either the genetic or the morphological diversification, which allude to biotic interactions such as symbiosis, as the driver of speciation in Globigerinoides.

Nomenclature for the Nameless: A Proposal for an Integrative Molecular Taxonomy of Cryptic Diversity Exemplified by Planktonic Foraminifera.

Investigations of biodiversity, biogeography, and ecological processes rely on the identification of "species" as biologically significant, natural units of evolution. In this context, morphotaxonomy only provides an adequate level of resolution if reproductive isolation matches morphological divergence. In many groups of organisms, morphologically defined species often disguise considerable genetic diversity, which may be indicative of the existence of cryptic species. The diversity hidden by morphological species can be disentangled through genetic surveys, which also provide access to data on the ecological distribution of genetically circumscribed units. These units can be identified by unique DNA sequence motifs and allow studies of evolutionary and ecological processes at different levels of divergence. However, the nomenclature of genetically circumscribed units within morphological species is not regulated and lacks stability. This represents a major obstacle to efforts to synthesize and communicate data on genetic diversity for multiple stakeholders. We have been confronted with such an obstacle in our work on planktonic foraminifera, where the stakeholder community is particularly diverse, involving geochemists, paleoceanographers, paleontologists, and biologists, and the lack of stable nomenclature beyond the level of formal morphospecies prevents effective transfer of knowledge. To circumvent this problem, we have designed a stable, reproducible, and flexible nomenclature system for genetically circumscribed units, analogous to the principles of a formal nomenclature system. Our system is based on the definition of unique DNA sequence motifs collocated within an individual, their typification (in analogy with holotypes), utilization of their hierarchical phylogenetic structure to define levels of divergence below that of the morphospecies, and a set of nomenclature rules assuring stability. The resulting molecular operational taxonomic units remain outside the domain of current nomenclature codes, but are linked to formal morphospecies as regulated by the codes. Subsequently, we show how this system can be applied to classify genetically defined units using the SSU rDNA marker in planktonic foraminifera and we highlight its potential use for other groups of organisms where similarly high levels of connectivity between molecular and formal taxonomies can be achieved.

PFR²: a curated database of planktonic foraminifera 18S ribosomal DNA as a resource for studies of plankton ecology, biogeography and evolution.

Planktonic foraminifera (Rhizaria) are ubiquitous marine pelagic protists producing calcareous shells with conspicuous morphology. They play an important role in the marine carbon cycle, and their exceptional fossil record serves as the basis for biochronostratigraphy and past climate reconstructions. A major worldwide sampling effort over the last two decades has resulted in the establishment of multiple large collections of cryopreserved individual planktonic foraminifera samples. Thousands of 18S rDNA partial sequences have been generated, representing all major known morphological taxa across their worldwide oceanic range. This comprehensive data coverage provides an opportunity to assess patterns of molecular ecology and evolution in a holistic way for an entire group of planktonic protists. We combined all available published and unpublished genetic data to build PFR(2), the Planktonic foraminifera Ribosomal Reference database. The first version of the database includes 3322 reference 18S rDNA sequences belonging to 32 of the 47 known morphospecies of extant planktonic foraminifera, collected from 460 oceanic stations. All sequences have been rigorously taxonomically curated using a six-rank annotation system fully resolved to the morphological species level and linked to a series of metadata. The PFR(2) website, available at http://pfr2.sb-roscoff.fr, allows downloading the entire database or specific sections, as well as the identification of new planktonic foraminiferal sequences. Its novel, fully documented curation process integrates advances in morphological and molecular taxonomy. It allows for an increase in its taxonomic resolution and assures that integrity is maintained by including a complete contingency tracking of annotations and assuring that the annotations remain internally consistent.

SSU rDNA divergence in planktonic foraminifera: molecular taxonomy and biogeographic implications.

The use of planktonic foraminifera in paleoceanography requires taxonomic consistency and precise assessment of the species biogeography. Yet, ribosomal small subunit (SSUr) DNA analyses have revealed that most of the modern morpho-species of planktonic foraminifera are composed of a complex of several distinct genetic types that may correspond to cryptic or pseudo-cryptic species. These genetic types are usually delimitated using partial sequences located at the 3'end of the SSUrDNA, but typically based on empirical delimitation. Here, we first use patristic genetic distances calculated within and among genetic types of the most common morpho-species to show that intra-type and inter-type genetic distances within morpho-species may significantly overlap, suggesting that genetic types have been sometimes inconsistently defined. We further apply two quantitative and independent methods, ABGD (Automatic Barcode Gap Detection) and GMYC (General Mixed Yule Coalescent) to a dataset of published and newly obtained partial SSU rDNA for a more objective assessment of the species status of these genetic types. Results of these complementary approaches are highly congruent and lead to a molecular taxonomy that ranks 49 genetic types of planktonic foraminifera as genuine (pseudo)cryptic species. Our results advocate for a standardized sequencing procedure allowing homogenous delimitations of (pseudo)cryptic species. On the ground of this revised taxonomic framework, we finally provide an integrative taxonomy synthesizing geographic, ecological and morphological differentiations that can occur among the genuine (pseudo)cryptic species. Due to molecular, environmental or morphological data scarcities, many aspects of our proposed integrative taxonomy are not yet fully resolved. On the other hand, our study opens up the potential for a correct interpretation of environmental sequence datasets.

Longitudinal differentiation among pelagic populations in a planktic foraminifer.

Evolutionary processes in marine plankton have been assumed to be dependent on the oceanic circulation system, which transports plankton between populations in marine surface waters. Gene flow facilitated by oceanic currents along longitudinal gradients may efficiently impede genetic differentiation of pelagic populations in the absence of confounding marine environmental effects. However, how responsible oceanic currents are for the geographic distribution and dispersal of plankton is poorly understood. We examined the phylogeography of the planktic foraminifer Pulleniatina obliquiloculata in the Indo-Pacific Warm Pool (IPWP) by using partial small subunit ribosomal DNA (SSU rDNA) sequences. We found longitudinal clines in the frequencies of three distinct genetic types in the IPWP area. These frequencies were correlated with environmental factors that are characteristic of three water masses in the IPWP. Noteworthy, populations inhabiting longitudinally distant water masses at the Pacific and Indian sides of the IPWP were genetically different, despite transportation of individuals via oceanic currents. These results demonstrate that populations of pelagic plankton have diverged genetically among different water masses within a single climate zone. Changes of the oceanic circulation system could have impacted the geographic patterns of dispersal and divergence of pelagic plankton.

Worldwide genotyping in the planktonic foraminifer Globoconella inflata: implications for life history and paleoceanography.

The planktonic foraminiferal morpho-species Globoconella inflata is widely used as a stratigraphic and paleoceanographic index. While G. inflata was until now regarded as a single species, we show that it rather constitutes a complex of two pseudo-cryptic species. Our study is based on SSU and ITS rDNA sequence analyses and genotyping of 497 individuals collected at 49 oceanic stations covering the worldwide range of the morpho-species. Phylogenetic analyses unveil the presence of two divergent genotypes. Type I inhabits transitional and subtropical waters of both hemispheres, while Type II is restricted to the Antarctic subpolar waters. The two genetic species exhibit a strictly allopatric distribution on each side of the Antarctic Subpolar Front. On the other hand, sediment data show that G. inflata was restricted to transitional and subtropical environments since the early Pliocene, and expanded its geographic range to southern subpolar waters ∼700 kyrs ago, during marine isotopic stage 17. This datum may correspond to a peripatric speciation event that led to the partition of an ancestral genotype into two distinct evolutionary units. Biometric measurements performed on individual G. inflata from plankton tows north and south of the Antarctic Subpolar Front indicate that Types I and II display slight but significant differences in shell morphology. These morphological differences may allow recognition of the G. inflata pseudo-cryptic species back into the fossil record, which in turn may contribute to monitor past movements of the Antarctic Subpolar Front during the middle and late Pleistocene.

Stable sea surface temperatures in the western Pacific warm pool over the past 1.75 million years.

About 850,000 years ago, the period of the glacial cycles changed from 41,000 to 100,000 years. This mid-Pleistocene climate transition has been attributed to global cooling, possibly caused by a decrease in atmospheric carbon dioxide concentrations. However, evidence for such cooling is currently restricted to the cool upwelling regions in the eastern equatorial oceans, although the tropical warm pools on the western side of the ocean basins are particularly sensitive to changes in radiative forcing. Here we present high-resolution records of sea surface temperatures spanning the past 1.75 million years, obtained from oxygen isotopes and Mg/Ca ratios in planktonic foraminifera from the western Pacific warm pool. In contrast with the eastern equatorial regions, sea surface temperatures in the western Pacific warm pool are relatively stable throughout the Pleistocene epoch, implying little long-term change in the tropical net radiation budget. Our results challenge the hypothesis of a gradual decrease in atmospheric carbon dioxide concentrations as a dominant trigger of the longer glacial cycles since 850,000 years ago. Instead, we infer that the temperature contrast across the equatorial Pacific Ocean increased, which might have had a significant influence on the mid-Pleistocene climate transition.

Evidence for large methane releases to the atmosphere from deep-sea gas-hydrate dissociation during the last glacial episode.

Past atmospheric methane-concentration oscillations recorded in polar ice cores vary together with rapid global climatic changes during the last glacial episode. In the "clathrate gun hypothesis," massive releases of deep-sea methane from marine gas-hydrate dissociation led to these well known, global, abrupt warmings in the past. If evidence for such releases in the water column exists, however, the mechanism and eventual transfer to the atmosphere has not yet been documented clearly. Here we describe a high-resolution marine-sediment record of stable carbon isotopic changes from the Papua Gulf, off Papua New Guinea, which exhibits two extremely depleted excursions (down to -9 per thousand ) at approximately 39,000 and approximately 55,000 years. Morphological, isotopic, and trace metal evidence dismisses authigenic calcite as the main source of depleted carbon. Massive methane release associated with deep-sea gas-hydrate dissociation is the most likely cause for such large depletions of delta(13)C. The absence of a delta(13)C gradient in the water column during these events implies that the methane rose through the entire water column, reaching the sea-air interface and thus the atmosphere. Foraminiferal delta(18)O composition suggests that the rise of the methane in the water column created an upwelling flow. These inferred emission events suggest that during the last glacial episode, this process was likely widespread, including tropical regions. Thus, the release of methane from the ocean floor into the atmosphere cannot be dismissed as a strong positive feedback in climate dynamics processes.