Raphaël Morard

Eukaryotic plankton diversity in the sunlit ocean

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

The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote Small Sub-Unit rRNA sequences with curated taxonomy

The interrogation of genetic markers in environmental meta-barcoding studies is currently seriously hindered by the lack of taxonomically curated reference data sets for the targeted genes. The Protist Ribosomal Reference database (PR2, http://ssu-rrna.org/) provides a unique access to eukaryotic small sub-unit (SSU) ribosomal RNA and DNA sequences, with curated taxonomy. The database mainly consists of nuclear-encoded protistan sequences. However, metazoans, land plants, macrosporic fungi and eukaryotic organelles (mitochondrion, plastid and others) are also included because they are useful for the analysis of high-troughput sequencing data sets. Introns and putative chimeric sequences have been also carefully checked. Taxonomic assignation of sequences consists of eight unique taxonomic fields. In total, 136 866 sequences are nuclear encoded, 45 708 (36 501 mitochondrial and 9657 chloroplastic) are from organelles, the remaining being putative chimeric sequences. The website allows the users to download sequences from the entire and partial databases (including representative sequences after clustering at a given level of similarity). Different web tools also allow searches by sequence similarity. The presence of both rRNA and rDNA sequences, taking into account introns (crucial for eukaryotic sequences), a normalized eight terms ranked-taxonomy and updates of new GenBank releases were made possible by a long-term collaboration between experts in taxonomy and computer scientists.

The cryptic and the apparent reversed: lack of genetic differentiation within the morphologically diverse plexus of the planktonic foraminifer Globigerinoides sacculifer

Abstract Previous genetic studies of extant planktonic foraminifera have provided evidence that the traditional, strictly morphological definition of species in these organisms underestimates their biodiversity. Here, we report the first case where this pattern is reversed. The modern (sub)tropical species plexus Globigerinoides sacculifer is characterized by large morphological variability, which has led to the proliferation of taxonomic names attributed to morphological end-members within the plexus. In order to clarify the taxonomic status of its morphotypes and to investigate the genetic connectivity among its currently partly disjunct (sub)tropical populations, we carried out a global survey of two ribosomal RNA regions (SSU and ITS-1) in all recent morphotypes of the plexus collected throughout (sub)tropical surface waters of the global ocean. Unexpectedly, we find an extremely reduced genetic variation within the plexus and no correlation between genetic and morphological divergence, suggesting taxonomical overinterpretation. The genetic homogeneity within the morphospecies is unexpected, considering its partly disjunct range in the (sub)tropical Atlantic and Indo-Pacific and its old age (early Miocene). A sequence variant in the rapidly evolving ITS-1 region indicates the existence of an exclusively Atlantic haplotype, which suggests an episode of relatively recent (last glacial) isolation, followed by subsequent resumption of unidirectional gene flow from the Indo-Pacific into the Atlantic. This is the first example in planktonic foraminifera where the morphological variability in a morphospecies exceeds its rDNA genetic variability. Such evidence for inconsistent scaling of morphological and genetic diversity in planktonic foraminifera could complicate the interpretation of evolutionary patterns in their fossil record.

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.

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.

Recent Invasion of the Symbiont-Bearing Foraminifera Pararotalia into the Eastern Mediterranean Facilitated by the Ongoing Warming Trend.

The eastern Mediterranean is a hotspot of biological invasions. Numerous species of Indo-pacific origin have colonized the Mediterranean in recent times, including tropical symbiont-bearing foraminifera. Among these is the species Pararotalia calcariformata. Unlike other invasive foraminifera, this species was discovered only two decades ago and is restricted to the eastern Mediterranean coast. Combining ecological, genetic and physiological observations, we attempt to explain the recent invasion of this species in the Mediterranean Sea. Using morphological and genetic data, we confirm the species attribution to P. calcariformata McCulloch 1977 and identify its symbionts as a consortium of diatom species dominated by Minutocellus polymorphus. We document photosynthetic activity of its endosymbionts using Pulse Amplitude Modulated Fluorometry and test the effects of elevated temperatures on growth rates of asexual offspring. The culturing of asexual offspring for 120 days shows a 30-day period of rapid growth followed by a period of slower growth. A subsequent 48-day temperature sensitivity experiment indicates a similar developmental pathway and high growth rate at 28°C, whereas an almost complete inhibition of growth was observed at 20°C and 35°C. This indicates that the offspring of this species may have lower tolerance to cold temperatures than what would be expected for species native to the Mediterranean. We expand this hypothesis by applying a Species Distribution Model (SDM) based on modern occurrences in the Mediterranean using three environmental variables: irradiance, turbidity and yearly minimum temperature. The model reproduces the observed restricted distribution and indicates that the range of the species will drastically expand westwards under future global change scenarios. We conclude that P. calcariformata established a population in the Levant because of the recent warming in the region. In line with observations from other groups of organisms, our results indicate that continued warming of the eastern Mediterranean will facilitate the invasion of more tropical marine taxa into the Mediterranean, disturbing local biodiversity and ecosystem structure.

PFR2: 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 PFR2, 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 PFR2 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.

Methodology for single-cell genetic analysis of planktonic foraminifera for studies of protist diversity and evolution

Single-cell genetic analysis is an essential method to investigate the biodiversity and evolutionary ecology of marine protists. In protist groups that do not reproduce under laboratory conditions, this approach provides the only means to directly associate molecular sequences with cell morphology. The resulting unambiguous taxonomic identification of the DNA sequences is a prerequisite for barcoding and analyses of environmental metagenomic data. Extensive single-cell genetic studies have been carried out on planktonic foraminifera over the past 20 years to elucidate their phylogeny, cryptic diversity, biogeography and the relationship between genetic and morphological variability. In the course of these investigations, it has become evident that genetic analysis at the individual specimen level is confronted by innumerable challenges ranging from the negligible amount of DNA present in the single cell to the substantial amount of DNA contamination introduced by endosymbionts or food particles. Consequently, a range of methods has been developed and applied throughout the years for the genetic analysis of planktonic foraminifera in order to enhance DNA amplification success rates. Yet, the description of these methods in the literature rarely occurred with equivalent levels of detail and the different approaches have never been compared in terms of their efficiency and reproducibility. Here, aiming at a standardization of methods, we provide a comprehensive review of all methods that have been employed for the single-cell genetic analysis of planktonic foraminifera. We compile data on success rates of DNA amplification and use these to evaluate the effects of key parameters associated with the methods of sample collection, storage and extraction of single-cell DNA. We show that the chosen methods influence the success rates of single-cell genetic studies, but the differences between them are not sufficient to hinder comparisons between studies carried out by different methods. The review thus not only provides a comprehensive reference with guidelines for future genetic studies on foraminifera, but it also establishes an important benchmark for investigations using existing single-cell datasets. The methods are widely applicable and the review may help to establish similar standard principles for their utilization in other protist groups.

Caught in the act: anatomy of an ongoing benthic–planktonic transition in a marine protist

The transition from benthos to plankton requires multiple adaptations, yet so far it remains unclear how these are acquired in the course of the transition. To investigate this process, we analyzed the genetic diversity and distribution patterns of a group of foraminifera of the genus Bolivina with a tychopelagic mode of life (same species occurring both in benthos and plankton). We assembled a global sequence data set for this group from single-cell DNA extractions and occurrences in metabarcodes from pelagic environmental samples. The pelagic sequences all cluster within a single monophyletic clade within Bolivina. This clade harbors three distinct genetic lineages, which are associated with incipient morphological differentiation. All lineages occur in the plankton and benthos, but only one lineage exhibits no limit to offshore dispersal and has been shown to grow in the plankton. These observations indicate that the emergence of buoyancy regulation within the clade preceded the evolution of pelagic feeding and that the evolution of both traits was not channeled into a full transition into the plankton. We infer that in foraminifera, colonization of the planktonic niche may occur by sequential cooptation of independently acquired traits, with holoplanktonic species being recruited from tychopelagic ancestors

Tracing shifts of oceanic fronts using the cryptic diversity of the planktonic foraminifera Globorotalia inflata

The use of planktonic foraminifera in paleoceanographic studies relies on the assumption that morphospecies represent biological species with ecological preferences that are stable through time and space. However, genetic surveys unveiled a considerable level of diversity in most morphospecies of planktonic foraminifera. This diversity is significant for paleoceanographic applications because cryptic species were shown to display distinct ecological preferences that could potentially help refine paleoceanographic proxies. Subtle morphological differences between cryptic species of planktonic foraminifera have been reported, but so far their applicability within paleoceanographic studies remains largely unexplored. Here we show how information on genetic diversity can be transferred to paleoceanography using Globorotalia inflata as a case study. The two cryptic species of G. inflata are separated by the Brazil-Malvinas Confluence (BMC), a major oceanographic feature in the South Atlantic. Based on this observation, we developed a morphological model of cryptic species detection in core top material. The application of the cryptic species detection model to Holocene samples implies latitudinal oscillations in the position of the confluence that are largely consistent with reconstructions obtained from stable isotope data. We show that the occurrence of cryptic species in G. inflata, can be detected in the fossil record and used to trace the migration of the BMC. Since a similar degree of morphological separation as in G. inflata has been reported from other species of planktonic foraminifera, the approach presented in this study can potentially yield a wealth of new paleoceanographical proxies.