Geographical distribution of cryptic genetic types in the planktonic foraminifer Globigerinoides ruber

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.

Latitudinal Differentiation among Modern Planktonic Foraminiferal Populations of Central Mediterranean: Species–Specific Distribution Patterns and Size Variability

Studies of the spatial distribution and size of modern planktonic foraminifera are still lacking in the Mediterranean Sea. In this study, 17 core-top sediments collected from a north-south transect along the central Mediterranean have been analyzed for planktonic foraminiferal content, in terms of their distributional pattern and intraspecific size variability. Among the analyzed planktonic foraminiferal species, Globigerina bulloides and Globigerinoides ruber (w) were the most abundant, presenting an antagonistic behavior and an overall decreasing trend in their average size values from Adriatic to Ionian sub-basins. Intraspecific differences have been also documented for G. ruber (w), with the dominant sensu stricto morphotype to present generally higher frequencies and more constant shell sizes than sensu lato. The greater size variability of the latter is possibly related to its adaptation in particular hydrographic conditions based on its depth habitat preference and ecological characteristics to reach the (sub)optimum growth conditions. The rest of the species occur in minor percentages and show on average 11% increase with decreasing latitude characterized by distinct species-specific size variations along the transect. Our results show that the relationship between planktonic foraminifera shell size and abundance or sea surface temperature are either absent or weaker than previously reported for other regions and that in central Mediterranean assemblages’ size may be mainly related to nutrient availability. Besides the environmental parameters (sea surface temperature, primary productivity, water depth, stratification), the possible hidden cryptic diversity, still lingers to be consistently determined, could give a better understanding of the geographic and morphological differentiation within the Mediterranean planktonic populations.

Redescription of Gnorimosphaeroma oregonense (Dana, 1853) (Crustacea, Isopoda, Sphaeromatidae), designation of neotype, and 16S-rDNA molecular phylogeny of the north-eastern Pacific species

Gnorimosphaeroma oregonense (Dana, 1852) is revised, a male neotype is designated, photographed, and illustrated; the species occurs from Vancouver British Columbia to the central California coast. 16S-rDNA sequences (~650 bp) for all available ethanol preserved species of Gnorimosphaeroma were used to hypothesize their relationships. Our analyses revealed a sister taxon relationship between the fully marine G. oregonense and the brackish to freshwater species, G. noblei. The oyster associated and introduced G. rayi is sister to a previously not recognized or identified, but genetically distinct, Gnorimosphaeroma sp. collected at two sites in San Francisco Bay. Gnorimosphaeroma sp. is probably also a western Pacific species based on its genetic relationship to G. rayi. Photographic comparisons are offered for G. oregonense (marine), G. noblei (freshwater), G. rayi (introduced), G. sp. (presumably introduced), and G. insulare (San Nicolas Island). Records of the holdings at the Los Angeles County Museum of Natural History are summarized. Without material available north of Vancouver through Alaska, the range of G. oregonense could not be genetically verified. This review includes a diagnosis and description of the genus Gnorimosphaeroma Menzies, 1954.

Key Environmental Factors Controlling Planktonic Foraminiferal and Pteropod Community’s Response to Late Quaternary Hydroclimate Changes in the South Aegean Sea (Eastern Mediterranean)

A multidisciplinary study was conducted in order to investigate the environmental factors affecting the planktonic foraminiferal and pteropod communities of the south Aegean Sea. Aspects of the Late Quaternary paleoceanographic evolution were revealed by means of quantitative analyses of planktonic foraminiferal and pteropod assemblages (including multivariate statistical approach; principal component analysis (PCA)), the oxygen (δ18O) and carbon (δ13C) isotopic composition of planktonic foraminifera and related paleoceanographic (planktonic paleoclimatic curve (PPC), productivity (E-index), stratification (S-index), seasonality) indices, extracted by the gravity core KIM-2A derived from the submarine area between Kimolos and Sifnos islands. Focusing on the last ~21 calibrated thousands of years before present (ka BP), cold and eutrophicated conditions were identified during the Late Glacial period (21.1–15.7 ka BP) and were followed by warmer and wetter conditions during the deglaciation phase. The beginning of the Holocene was marked by a climatic amelioration and increased seasonality. The more pronounced environmental changes were identified during the deposition of the sapropel sublayers S1a (9.4–7.7 ka BP) and S1b (6.9–6.4 ka BP), with extremely warm and stratified conditions. Pteropod fauna during the sapropel deposition were recorded for the first time in the south Aegean Sea, suggesting arid conditions towards the end of S1a. Besides sea surface temperature (SST), which shows the highest explanatory power for the distribution of the analyzed fauna, water column stratification, primary productivity, and seasonality also control their communities during the Late Quaternary.

High dispersal capacity and biogeographic breaks shape the genetic diversity of a globally distributed reef-dwelling calcifier

Understanding the role of dispersal and adaptation in the evolutionary history of marine species is essential for predicting their response to changing conditions. We analyzed patterns of genetic differentiation in the key tropical calcifying species of large benthic foraminifera Amphistegina lobifera to reveal the evolutionary processes responsible for its biogeographic distribution. We collected specimens from 16 sites encompassing the entire range of the species and analyzed hypervariable fragments of the 18S SSU rDNA marker. We identified six hierarchically organized genotypes with mutually exclusive distribution organized along a longitudinal gradient. The distribution is consistent with diversification occurring in the Indo-West Pacific (IWP) followed by dispersal toward the periphery. This pattern can be explained by: (a) high dispersal capacity of the species, (b) habitat heterogeneity driving more recent differentiation in the IWP, and (c) ecological-scale processes such as niche incumbency reinforcing patterns of genotype mutual exclusion. The dispersal potential of this species drives the ongoing range expansion into the Mediterranean Sea, indicating that A. lobifera is able to expand its distribution by tracking increases in temperature. The genetic structure reveals recent diversification and high rate of extinction in the evolutionary history of the clade suggesting a high turnover rate of the diversity at the cryptic level. This diversification dynamic combined with high dispersal potential, allowed the species to maintain a widespread distribution over periods of geological and climatic upheaval. These characteristics are likely to allow the species to modify its geographic range in response to ongoing global warming without requiring genetic differentiation.

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.

Speciation in protists: Spatial and ecological divergence processes cause rapid species diversification in a freshwater chrysophyte

Although eukaryotic microorganisms are extremely numerous, diverse and essential to global ecosystem functioning, they are largely understudied by evolutionary biologists compared to multicellular macroscopic organisms. In particular, very little is known about the speciation mechanisms which may give rise to the diversity of microscopic eukaryotes. It was postulated that the enormous population sizes and ubiquitous distribution of these organisms could lead to a lack of population differentiation and therefore very low speciation rates. However, such assumptions have traditionally been based on morphospecies, which may not accurately reflect the true diversity, missing cryptic taxa. In this study, we aim to articulate the major diversification mechanisms leading to the contemporary molecular diversity by using a colonial freshwater flagellate, Synura sphagnicola, as an example. Phylogenetic analysis of five sequenced loci showed that S. sphagnicola differentiated into two morphologically distinct lineages approximately 15.4 million years ago, which further diverged into several evolutionarily recent haplotypes during the late Pleistocene. The most recent haplotypes are ecologically and biogeographically much more differentiated than the old lineages, presumably because of their persistent differentiation after the allopatric speciation events. Our study shows that in microbial eukaryotes, species diversification via the colonization of new geographical regions or ecological resources occurs much more readily than was previously thought. Consequently, divergence times of microorganisms in some lineages may be equivalent to the estimated times of speciation in plants and animals.

Repeated species radiations in the recent evolution of the key marine phytoplankton lineage Gephyrocapsa

Phytoplankton account for nearly half of global primary productivity and strongly affect the global carbon cycle, yet little is known about the forces that drive the evolution of these keystone microscopic organisms. Here we combine morphometric data from the fossil record of the ubiquitous coccolithophore genus Gephyrocapsa with genomic analyses of extant species to assess the genetic processes underlying Pleistocene palaeontological patterns. We demonstrate that all modern diversity in Gephyrocapsa (including Emiliania huxleyi) originated in a rapid species radiation during the last 0.6 Ma, coincident with the latest of the three pulses of Gephyrocapsa diversification and extinction documented in the fossil record. Our evolutionary genetic analyses indicate that new species in this genus have formed in sympatry or parapatry, with occasional hybridisation between species. This sheds light on the mode of speciation during evolutionary radiation of marine phytoplankton and provides a model of how new plankton species form.

The phytoplankton Gephyrocapsa have gone through repeated macroevolutionary shifts in size. Here, Bendif et al. combine fossil and genomic data to show the latest shift was coincident with a species radiation and suggest that previous shifts have also resulted from cycles of radiation and extinction.

Metabarcoding Insights Into the Trophic Behavior and Identity of Intertidal Benthic Foraminifera

Foraminifera are ubiquitous marine protists with an important role in the benthic carbon cycle. However, morphological observations often fail to resolve their exact taxonomic placement and there is a lack of field studies on their particular trophic preferences. Here, we propose the application of metabarcoding as a tool for the elucidation of the in situ feeding behavior of benthic foraminifera, while also allowing the correct taxonomic assignment of the feeder, using the V9 region of the 18S (small subunit; SSU) rRNA gene. Living foraminiferal specimens were collected from two intertidal mudflats of the Wadden Sea and DNA was extracted from foraminiferal individuals and from the surrounding sediments. Molecular analysis allowed us to confirm that our foraminiferal specimens belong to three genetic types: Ammonia sp. T6, Elphidium sp. S5 and Haynesina sp. S16. Foraminiferal intracellular eukaryote communities reflected to an extent those of the surrounding sediments but at different relative abundances. Unlike sediment eukaryote communities, which were largely determined by the sampling site, foraminiferal intracellular eukaryote communities were driven by foraminiferal species, followed by sediment depth. Our data suggests that Ammonia sp. T6 can predate on metazoan classes, whereas Elphidium sp. S5 and Haynesina sp. S16 are more likely to ingest diatoms. These observations, alongside the use of metabarcoding in similar ecological studies, significantly contribute to our overall understanding of the ecological roles of these protists in intertidal benthic environments and their position and function in the benthic food webs.

Advanced approach to analyzing calcareous protists for present and past pelagic ecology: Comprehensive analysis of 3D-morphology, stable isotopes, and genes of planktic foraminifers

Marine protists play an important role in oceanic ecosystems and biogeochemical cycles. However, the difficulties in culturing pelagic protists indicate that their ecology and behavior remain poorly understood; phylogeographic studies based on single-cell genetic analyses have often shown that they are highly divergent at the biological species level, with variable geographic distributions. This indicates that their ecology could be complex. On the other hand, the biomineral (calcareous) shells of planktic foraminifers are widely used in geochemical analyses to estimate marine paleoenvironmental characteristics (i.e., temperature), because the shell chemical composition reflects ambient seawater conditions. Among the pelagic protists, planktic foraminifers are ideal study candidates to develop a combined approach of genetic, morphological, and geochemical methods, thus reflecting environmental and ecological characteristics. The present study precisely tested whether the DNA extraction process physically and chemically affects the shells of the planktic foraminifer Globigerinoides ruber. We used a nondestructive method for analyzing physical changes (micro-focus X-ray computed tomography (MXCT) scanning) to compare specimens at the pre- and post-DNA extraction stages. Our results demonstrate that DNA extraction has no significant effect on shell density and thickness. We measured stable carbon and oxygen isotopes on the shell of each individual in a negative control or one of two DNA-extracted groups and detected no significant differences in isotopic values among the three groups. Moreover, we evaluated isotopic variations at the biological species level with regard to their ecological characteristics such as depth habitat, life stages, and symbionts. Thus, our examination of the physiochemical effects on biomineral shells through DNA extraction shows that morphological and isotopic analyses of foraminifers can be combined with genetic analysis. These analytical methods are applicable to other shell-forming protists and microorganisms. In this study, we developed a powerful analytical tool for use in ecological and environmental studies of modern and past oceans.

Peijnenburg K. Biogeography and genetic diversity of the atlantid heteropods

The atlantid heteropods are regularly encountered, but rarely studied marine planktonic gastropods. Relying on a small (<14 mm), delicate aragonite shell and living in the upper ocean means that, in common with pteropods, atlantids are likely to be affected by imminent ocean changes. Variable shell morphology and widespread distributions indicate that the family is more diverse than the 23 currently known species. Uncovering this diversity is fundamental to determining the distribution of atlantids and to understanding their environmental tolerances. Here we present phylogenetic analyses of all described species of the family Atlantidae using 437 new and 52 previously published cytochrome c oxidase subunit 1 mitochondrial DNA (mtCO1) sequences. Specimens and published sequences were gathered from 32 Atlantic Ocean stations, 14 Indian Ocean stations and 21 Pacific Ocean stations between 35°N and 43°S. DNA barcoding and Automatic Barcode Gap Discovery (ABGD) proved to be valuable tools for the identification of described atlantid species, and also revealed ten additional distinct clades, suggesting that the diversity within this family has been underestimated. Only two of these clades displayed obvious morphological characteristics, demonstrating that much of the newly discovered diversity is hidden from morphology-based identification techniques. Investigation of six large atlantid collections demonstrated that 61% of previously described (morpho) species have a circumglobal distribution. Of the remaining 39%, two species were restricted to the Atlantic Ocean, five occurred in the Indian and Pacific oceans, one species was only found in the northeast Pacific Ocean, and one occurred only in the Southern Subtropical Convergence Zone. Molecular analysis showed that seven of the species with wide distributions were comprised of two or more clades that occupied distinct oceanographic regions. These distributions may suggest narrower environmental tolerances than the described morphospecies. Results provide an updated biogeography and mtCO1 reference dataset of the Atlantidae that may be used to identify atlantid species and provide a first step in understanding their evolutionary history and accurate distribution, encouraging the inclusion of this family in future plankton research.

Assessing SSU rDNA Barcodes in Foraminifera: A Case Study using Bolivina quadrata

The Small Subunit Ribosomal RNA gene (SSU rDNA) is a widely used tool to reconstruct phylogenetic relationships among foraminiferal species. Recently, the highly variable regions of this gene have been proposed as DNA barcodes to identify foraminiferal species. However, the resolution of these barcodes has not been well established, yet. In this study, we evaluate four SSU rDNA hypervariable regions (37/f, 41/f, 43/e, and 45/e) as DNA barcodes to distinguish among species of the genus Bolivina, with particular emphasis on Bolivina quadrata for which ten new sequences (KY468817-KY468826) were obtained during this study. Our analyses show that a single SSU rDNA hypervariable sequence is insufficient to resolve all Bolivina species and that some regions (37/f and 41/f) are more useful than others (43/e and 45/e) to distinguish among closely related species. In addition, polymorphism analyses reveal a high degree of variability. In the context of barcoding studies, these results emphasize the need to assess the range of intraspecific variability of DNA barcodes prior to their application to identify foraminiferal species in environmental samples; our results also highlight the possibility that a longer SSU rDNA region might be required to distinguish among species belonging to the same taxonomic group (i.e. genus).

Clarifying the Cryptic Host Specificity of Blastocystis spp. Isolates from Alouatta palliata and A. pigra Howler Monkeys

Although the presence of cryptic host specificity has been documented in Blastocystis, differences in infection rates and high genetic polymorphism within and between populations of some subtypes (ST) have impeded the clarification of the generalist or specialist specificity of this parasite. We assessed the genetic variability and host specificity of Blastocystis spp. in wild howler monkeys from two rainforest areas in the southeastern region of Mexico. Fecal samples of 225 Alouatta palliata (59) and A. pigra (166) monkeys, belonging to 16 sylvatic sites, were analyzed for infection with Blastocystis ST using a region of the small subunit rDNA (SSUrDNA) gene as a marker. Phylogenetic and genetic diversity analyses were performed according to the geographic areas where the monkeys were found. Blastocystis ST2 was the most abundant (91.9%), followed by ST1 and ST8 with 4.6% and 3.5%, respectively; no association between Blastocystis ST and Alouatta species was observed. SSUrDNA sequences in GenBank from human and non-human primates (NHP) were used as ST references and included in population analyses. The haplotype network trees exhibited different distributions: ST1 showed a generalist profile since several haplotypes from different animals were homogeneously distributed with few mutational changes. For ST2, a major dispersion center grouped the Mexican samples, and high mutational differences were observed between NHP. Furthermore, nucleotide and haplotype diversity values, as well as migration and genetic differentiation indexes, showed contrasting values for ST1 and ST2. These data suggest that ST1 populations are only minimally differentiated, while ST2 populations in humans are highly differentiated from those of NHP. The host generalist and specialist specificities exhibited by ST1 and ST2 Blastocystis populations indicate distinct adaptation processes. Because ST1 exhibits a generalist profile, this haplotype can be considered a metapopulation; in contrast, ST2 exists as a set of local populations with preferences for either humans or NHP.

Evolution of a Planktonic Foraminifer during Environmental Changes in the Tropical Oceans

Ecological adaptation to environmental changes is a strong driver of evolution, enabling speciation of pelagic plankton in the open ocean without the presence of effective physical barriers to gene flow. The tropical ocean environment, which plays an important role in shaping marine biodiversity, has drastically and frequently changed since the Pliocene. Nevertheless, the evolutionary history of tropical pelagic plankton has been poorly understood, as phylogeographic investigations are still in the developing state and paleontological approaches are insufficient to obtain a sequential record from the deep-sea sediments. The planktonic foraminifer Pulleniatina obliquiloculata is widely distributed in the tropical area throughout the world’s oceans, and its phylogeography is well established. It is thus one of the best candidates to examine how past environmental changes may have shifted the spatial distribution and affected the diversification of tropical pelagic plankton. Such an examination requires the divergence history of the planktonic foraminifers, yet the gene marker (partial small subunit (SSU) rDNA) previously used for phylogeographic studies was not powerful enough to achieve a high accuracy in estimating the divergence times. The present study focuses on improving the precision of divergence time estimates for the splits between sibling species (genetic types) of planktonic foraminifers by increasing the number of genes as well as the number of nucleotide bases used for molecular clock estimates. We have amplified the entire coding regions of two ribosomal RNA genes (SSU rDNA and large subunit (LSU) rDNA) of three genetic types of P. obliquiloculata and two closely related species for the first time and applied them to the Bayesian relaxed clock method. The comparison of the credible intervals of the four datasets consisting either of sequences of the partial SSU rDNA, the complete SSU rDNA, LSU rDNA, or a combination of both genes (SSU+LSU) clearly demonstrated that the two-gene dataset improved the accuracy of divergence time estimates. The P. obliquiloculata lineage diverged twice, first at the end of the Pliocene (3.1 Ma) and again in the middle Pleistocene (1.4 Ma). Both timings coincided with the environmental changes, which indirectly involved geographic separation of populations. The habitat of P. obliquiloculata was expanded toward the higher latitudinal zones during the stable warm periods and subsequently placed on the steep environmental gradients following the global cooling. Different environmental conditions in the stable warm tropics and unstable higher latitudes may have triggered ecological divergence among the populations, leading to adaptive differentiation and eventually speciation. A comprehensive analysis of divergence time estimates combined with phylogeography enabled us to reveal the evolutionary history of the pelagic plankton and to find the potential paleoenvironmental events, which could have changed their biogeography and ecology.

Molecular Phylogeny and Ecology of Textularia agglutinans d'Orbigny from the Mediterranean Coast of Israel: A Case of a Successful New Incumbent

Textularia agglutinans d’Orbigny is a non-symbiont bearing and comparatively large benthic foraminiferal species with a widespread distribution across all oceans. In recent years, its populations have considerably expanded along the Israeli Mediterranean coast of the eastern Levantine basin. Despite its exceptionally widespread occurrence, no molecular data have yet been obtained. This study provides the first ribosomal DNA sequences of T. agglutinans complemented with morphological and ecological characterization, which are based on material collected during environmental monitoring of the hard bottom habitats along the Israeli Mediterranean coast, and from the Gulf of Elat (northern Red Sea). Our phylogenetic analyses reveal that all specimens from both provinces belong to the same genetic population, regardless their morphological variability. These results indicate that modern population of T. agglutinans found on the Mediterranean coast of Israel is probably Lessepsian. Our study also reveals that T. agglutinans has an epiphytic life mode, which probably enabled its successful colonization of the hard bottom habitats, at the Mediterranean coast of Israel, which consist of a diverse community of macroalgae. Our study further indicates that the species does not tolerate high SST (> 35°C), which will probably prevent its future expansion in the easternmost Mediterranean in light of the expected rise in temperatures.

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.

Uncovering sibling species in Radiolaria: evidence for ecological partitioning in a marine planktonic protist

Phylogeography of unicellular plankton, as representative pelagic organisms, is fundamental to understanding their evolution in the ocean. Historically, these microplankton were believed to have cosmopolitan distributions achieved through passive transport and little potential for speciation because of a lack of geographic barriers in the oceans. Recent phylogeographic studies of these microplankton, however, have often revealed high diversity and fine-scale geographic distributions. These apparent contradictions may result from poor knowledge of the spatial distributions of pelagic microplankton in the water column. More information about both geographic and vertical distributions of pelagic populations could reveal the dispersal pathways, gene flow, and resulting diversifications in the open ocean. Here we demonstrate that two genetic types of the radiolarian morphospecies Spongotrochus glacialis with morphological differences are vertically segregated into the upper and lower surface waters within the pycnocline of the North Pacific Subtropical Water. This vertically separated distribution of two sister species is associated with distinct ecological partitioning. These two species could survive on different food resources from their respective environments: one in oligotrophic surface waters by using nutrients from symbionts, and the other at greater depths by depending on both heterotrophic and symbiotic nutrition. Moreover, molecular divergence-time estimates suggest that the two species diverged during the period of oligotrophic surface-water development in the Pacific Ocean. Our findings suggest that genetic isolation in the vertical dimension occurs through ecological partitioning even in the absence of physical barriers in the pelagic oceans.

Phylogeography of the tropical planktonic foraminifera lineage globigerinella reveals isolation inconsistent with passive dispersal by ocean currents

Morphologically defined species of marine plankton often harbor a considerable level of cryptic diversity. Since many morphospecies show cosmopolitan distribution, an understanding of biogeographic and evolutionary processes at the level of genetic diversity requires global sampling. We use a database of 387 single-specimen sequences of the SSU rDNA of the planktonic foraminifera Globigerinella as a model to assess the biogeographic and phylogenetic distributions of cryptic diversity in marine microplankton on a global scale. Our data confirm the existence of multiple, well isolated genetic lineages. An analysis of their abundance and distribution indicates that our sampling is likely to approximate the actual total diversity. Unexpectedly, we observe an uneven allocation of cryptic diversity among the phylogenetic lineages. We show that this pattern is neither an artifact of sampling intensity nor a function of lineage age. Instead, we argue that it reflects an ongoing speciation process in one of the three major lineages. Surprisingly, four of the six genetic types in the hyperdiverse lineage are biogeographically restricted to the Indopacific. Their mutual co-occurrence and their hierarchical phylogenetic structure provide no evidence for an origin through sudden habitat fragmentation and their limitation to the Indopacific challenges the view of a global gene flow within the warm-water provinces. This phenomenon shows that passive dispersal is not sufficient to describe the distribution of plankton diversity. Rather, these organisms show differentiated distribution patterns shaped by species interactions and reflecting phylogenetic contingency with unique histories of diversification rates.

High cryptic diversity across the global range of the migratory planktonic copepods Pleuromamma piseki and P. gracilis

Although holoplankton are ocean drifters and exhibit high dispersal potential, a number of studies on single species are finding highly divergent genetic clades. These cryptic species complexes are important to discover and describe, as identification of common marine species is fundamental to understanding ecosystem dynamics. Here we investigate the global diversity within Pleuromamma piseki and P. gracilis, two dominant members of the migratory zooplankton assemblage in subtropical and tropical waters worldwide. Using DNA sequence data from the mitochondrial gene cytochrome c oxidase subunit II (mtCOII) from 522 specimens collected across the Pacific, Atlantic and Indian Oceans, we discover twelve well-resolved genetically distinct clades in this species complex (Bayesian posterior probabilities >0.7; 6.3–17% genetic divergence between clades). The morphologically described species P. piseki and P. gracilis did not form monophyletic groups, rather they were distributed throughout the phylogeny and sometimes co-occurred within well-resolved clades: this result suggests that morphological characters currently used for taxonomic identification of P. gracilis and P. piseki may be inaccurate as indicators of species’ boundaries. Cryptic clades within the species complex ranged from being common to rare, and from cosmopolitan to highly restricted in distribution across the global ocean. These novel lineages appear to be ecologically divergent, with distinct biogeographic distributions across varied pelagic habitats. We hypothesize that these mtDNA lineages are distinct species and suggest that resolving their systematic status is important, given the ecological significance of the genus Pleuromamma in subtropical-tropical waters worldwide.

Using environmental niche models to test the 'everything is everywhere' hypothesis for Badhamia

It is often discussed whether the biogeography of free-living protists is better explained by the 'everything is everywhere'(EiE) hypothesis, which postulates that only ecology drives their distribution, or by the alternative hypothesis of 'moderate endemicity' in which geographic barriers can limit their dispersal. To formally test this, it would be necessary not only to find organisms restricted to a geographical area but also to check for their presence in any other place with a similar ecology. We propose the use of environmental niche models to generate and test null EiE distributions. Here we have analysed the distribution of 18S rDNA variants (ribotypes) of the myxomycete Badhamia melanospora (belonging to the protozoan phylum Amoebozoa) using 125 specimens from 91 localities. Two geographically structured groups of ribotypes congruent with slight morphological differences in the spores can be distinguished. One group comprises all populations from Argentina and Chile, and the other is formed by populations from North America together with human-introduced populations from other parts of the world. Environmental climatic niche models constructed separately for the two groups have significant differences, but show several overlapping areas. However, only specimens from one group were found in an intensively surveyed area in South America where both niche models overlap. It can be concluded that everything is not everywhere for B. melanospora. This taxon constitutes a complex formed by at least two cryptic species that probably diverged allopatrically in North and South America.

High evolutionary potential of marine zooplankton

Open ocean zooplankton often have been viewed as slowly evolving species that have limited capacity to respond adaptively to changing ocean conditions. Hence, attention has focused on the ecological responses of zooplankton to current global change, including range shifts and changing phenology. Here, we argue that zooplankton also are well poised for evolutionary responses to global change. We present theoretical arguments that suggest plankton species may respond rapidly to selection on mildly beneficial mutations due to exceptionally large population size, and consider the circumstantial evidence that supports our inference that selection may be particularly important for these species. We also review all primary population genetic studies of open ocean zooplankton and show that genetic isolation can be achieved at the scale of gyre systems in open ocean habitats (100s to 1000s of km). Furthermore, population genetic structure often varies across planktonic taxa, and appears to be linked to the particular ecological requirements of the organism. In combination, these characteristics should facilitate adaptive evolution to distinct oceanographic habitats in the plankton. We conclude that marine zooplankton may be capable of rapid evolutionary as well as ecological responses to changing ocean conditions, and discuss the implications of this view. We further suggest two priority areas for future research to test our hypothesis of high evolutionary potential in open ocean zooplankton, which will require (1) assessing how pervasive selection is in driving population divergence and (2) rigorously quantifying the spatial and temporal scales of population differentiation in the open ocean. Recent attention has focused on the ecological responses of open ocean zooplankton to current global change, including range shifts and changing phenology. Here, we argue that marine zooplankton also are well poised for evolutionary responses to global change.

Assessing the role of cladogenesis in macroevolution by integrating fossil and molecular evidence

Assessing the extent to which population subdivision during cladogenesis is necessary for long-term phenotypic evolution is of fundamental importance in a broad range of biological disciplines. Differentiating cladogenesis from anagenesis, defined as evolution within a species, has generally been hampered by dating precision, insufficient fossil data, and difficulties in establishing a direct link between morphological changes detectable in the fossil record and biological species. Here we quantify the relative frequencies of cladogenesis and anagenesis for macroperforate planktic Foraminifera, which arguably have the most complete fossil record currently available, to address this question. Analyzing this record in light of molecular evidence, while taking into account the precision of fossil dating techniques, we estimate that the fraction of speciation events attributable to anagenesis is <19% during the Cenozoic era (last 65 Myr) and <10% during the Neogene period (last 23 Myr). Our central conclusion--that cladogenesis is the predominant mode by which new planktic Foraminifera taxa become established at macroevolutionary time scales--differs markedly from the conclusion reached in a recent study based solely on fossil data. These disparate findings demonstrate that interpretations of macroevolutionary dynamics in the fossil record can be fundamentally altered in light of genetic evidence.

A legacy of contrasting spatial genetic structure on either side of the Atlantic-Mediterranean transition zone in a marine protist

The mechanisms that underpin the varied spatial genetic structures exhibited by free-living marine microorganisms remain controversial, with most studies emphasizing a high dispersal capability that should redistribute genetic diversity in contrast to most macroorganisms whose populations often retain a genetic signature of demographic response to historic climate fluctuations. We quantified the European phylogeographic structure of the marine flagellate Oxyrrhis marina and found a marked difference in spatial genetic structure, population demography, and genetic diversity between the northwest Atlantic and Mediterranean Sea that reflects the persistent separation of these regions as well as context-dependent population responses to contrasting environments. We found similar geographic variation in the level of genetic diversity in the sister species Oxyrrhis maritima. Because the capacity for wide dispersal is not always realized, historic genetic footprints of range expansion and contraction persist in contemporary populations of marine microbes, as they do in larger species. Indeed, the well-described genetic effects of climatic variation on macroorganisms provide clear, testable hypotheses about the processes that drive genetic divergence in marine microbes and thus about the response to future environmental change.

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.

Vertical niche partitioning between cryptic sibling species of a cosmopolitan marine planktonic protist

A large portion of the surface-ocean biomass is represented by microscopic unicellular plankton. These organisms are functionally and morphologically diverse, but it remains unclear how their diversity is generated. Species of marine microplankton are widely distributed because of passive transport and lack of barriers in the ocean. How does speciation occur in a system with a seemingly unlimited dispersal potential? Recent studies using planktonic foraminifera as a model showed that even among the cryptic genetic diversity within morphological species, many genetic types are cosmopolitan, lending limited support for speciation by geographical isolation. Here we show that the current two-dimensional view on the biogeography and potential speciation mechanisms in the microplankton may be misleading. By depth-stratified sampling, we present evidence that sibling genetic types in a cosmopolitan species of marine microplankton, the planktonic foraminifer Hastigerina pelagica, are consistently separated by depth throughout their global range. Such strong separation between genetically closely related and morphologically inseparable genetic types indicates that niche partitioning in marine heterotrophic microplankton can be maintained in the vertical dimension on a global scale. These observations indicate that speciation along depth (depth-parapatric speciation) can occur in vertically structured microplankton populations, facilitating diversification without the need for spatial isolation.

Ecological partitioning and diversity in tropical planktonic foraminifera

BACKGROUND

Ecological processes are increasingly being viewed as an important mode of diversification in the marine environment, where the high dispersal potential of pelagic organisms, and a lack of absolute barriers to gene flow may limit the occurrence of allopatric speciation through vicariance. Here we focus on the potential role of ecological partitioning in the diversification of a widely distributed group of marine protists, the planktonic foraminifera. Sampling was conducted in the tropical Arabian Sea, during the southwest (summer) monsoon, when pronounced environmental conditions result in a strong disparity in temperature, salinity and productivity between distinct northern and southern water masses.

RESULTS

We uncovered extensive genetic diversity within the Arabian Sea planktonic foraminifera, identifying 13 morphospecies, represented by 20 distinct SSU rRNA genetic types. Several morphospecies/genetic types displayed non-random biogeographical distributions, partitioning between the northern and southern water masses, giving a strong indication of independent ecological adaptations.

CONCLUSIONS

We propose sea-surface primary productivity as the main factor driving the geographical segregation of Arabian Sea planktonic foraminifera, during the SW monsoon, with variations in symbiotic associations possibly playing a role in the specific ecological adaptations observed. Our findings suggest that ecological partitioning could be contributing to the high levels of 'cryptic' genetic diversity observed within the planktonic foraminifera, and support the view that ecological processes may play a key role in the diversification of marine pelagic organisms.

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.

Global diversity of aloricate Oligotrichea (Protista, Ciliophora, Spirotricha) in marine and brackish sea water

Oligotrichids and choreotrichids are ciliate taxa contributing to the multi-step microbial food web and episodically dominating the marine microzooplankton. The global diversity and distribution of aloricate Oligotrichea are unknown. Here, the geographic ranges of the 141 accepted species and their synonyms in marine and brackish sea water are analyzed, using hundreds of taxonomical and ecological studies; the quality of the records is simultaneously evaluated. The aloricate Oligotrichea match the moderate endemicity model, i.e., the majority (94) of morphospecies has a wide, occasionally cosmopolitan distribution, while 47 morphospecies show biogeographic patterns: they are restricted to single geographic regions and probably include 12 endemic morphospecies. These endemics are found in the Antarctic, North Pacific, and Black Sea, whereas the "flagship" species Strombidinopsis cercionis is confined to the Caribbean Sea. Concerning genera, again several geographic patterns are recognizable. The species richness is distinctly lower in the southern hemisphere than in the northern, ranging from nine morphospecies in the South Pacific to 95 in the North Atlantic; however, this pattern is probably caused by undersampling. Since the loss of species might affect higher trophical levels substantially, the aloricate Oligotrichea should not any longer be ignored in conservation issues. The ecophysiological diversity is considerably larger than the morphological, and even tops the richness of SSrRNA and ITS haplotypes, indicating that probably more than 83-89% of the diversity in aloricate Oligotrichea are unknown. The huge challenge to discover all these species can only be managed by combining the expertises of morphological taxonomists, molecular biologists, ecologists, and physiologists.

A Clustering Optimization Strategy for Molecular Taxonomy Applied to Planktonic Foraminifera SSU rDNA

Identifying species is challenging in the case of organisms for which primarily molecular data are available. Even if morphological features are available, molecular taxonomy is often necessary to revise taxonomic concepts and to analyze environmental DNA sequences. However, clustering approaches to delineate molecular operational taxonomic units often rely on arbitrary parameter choices. Also, distance calculation is difficult for highly alignment-ambiguous sequences. Here, we applied a recently described clustering optimization method to highly divergent planktonic foraminifera SSU rDNA sequences. We determined the distance function and the clustering setting that result in the highest agreement with morphological reference data. Alignment-free distance calculation, when adapted to the use with partly non-homologous sequences caused by distinct primer pairs, outperformed multiple sequence alignment. Clustering optimization offers new perspectives for the barcoding of species diversity and for environmental sequencing. It bridges the gap between traditional and modern taxonomic disciplines by specifically addressing the issue of how to optimally account for both genetic divergence and given species concepts.

Limits to gene flow in a cosmopolitan marine planktonic diatom

The role of geographic isolation in marine microbial speciation is hotly debated because of the high dispersal potential and large population sizes of planktonic microorganisms and the apparent lack of strong dispersal barriers in the open sea. Here, we show that gene flow between distant populations of the globally distributed, bloom-forming diatom species Pseudo-nitzschia pungens (clade I) is limited and follows a strong isolation by distance pattern. Furthermore, phylogenetic analysis implies that under appropriate geographic and environmental circumstances, like the pronounced climatic changes in the Pleistocene, population structuring may lead to speciation and hence may play an important role in diversification of marine planktonic microorganisms. A better understanding of the factors that control population structuring is thus essential to reveal the role of allopatric speciation in marine microorganisms.

Using the Multiple Analysis Approach to Reconstruct Phylogenetic Relationships among Planktonic Foraminifera from Highly Divergent and Length-polymorphic SSU rDNA Sequences

The high sequence divergence within the small subunit ribosomal RNA gene (SSU rDNA) of foraminifera makes it difficult to establish the homology of individual nucleotides across taxa. Alignment-based approaches so far relied on time-consuming manual alignments and discarded up to 50% of the sequenced nucleotides prior to phylogenetic inference. Here, we investigate the potential of the multiple analysis approach to infer a molecular phylogeny of all modern planktonic foraminiferal taxa by using a matrix of 146 new and 153 previously published SSU rDNA sequences. Our multiple analysis approach is based on eleven different automated alignments, analysed separately under the maximum likelihood criterion. The high degree of congruence between the phylogenies derived from our novel approach, traditional manually homologized culled alignments and the fossil record indicates that poorly resolved nucleotide homology does not represent the most significant obstacle when exploring the phylogenetic structure of the SSU rDNA in planktonic foraminifera. We show that approaches designed to extract phylogenetically valuable signals from complete sequences show more promise to resolve the backbone of the planktonic foraminifer tree than attempts to establish strictly homologous base calls in a manual alignment.