Longitudinal differentiation among pelagic populations in a planktic foraminifer

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.

Coiling directions in the planktonic foraminifer Pulleniatina: A complex eco-evolutionary dynamic spanning millions of years

Planktonic foraminifera are heterotrophic sexually reproducing marine protists with an exceptionally complete fossil record that provides unique insights into long-term patterns and processes of evolution. Populations often exhibit strong biases towards either right (dextral) or left (sinistral) shells. Deep-sea sediment cores spanning millions of years reveal that some species show large and often rapid fluctuations in their dominant coiling direction through time. This is useful for biostratigraphic correlation but further work is required to understand the population dynamical processes that drive these fluctuations. Here we address the case of coiling fluctuations in the planktonic foraminifer genus Pulleniatina based on new high-resolution counts from two recently recovered sediment cores from either side of the Indonesian through-flow in the tropical west Pacific and Indian Oceans (International Ocean Discovery Program Sites U1486 and U1483). We use single-specimen stable isotope analyses to show that dextral and sinistral shells from the same sediment samples can show significant differences in both carbon and oxygen isotopes, implying a degree of ecological separation between populations. In one case we detect a significant difference in size between dextral and sinistral specimens. We suggest that major fluctuations in coiling ratio are caused by cryptic populations replacing one another in competitive sweeps, a mode of evolution that is more often associated with asexual organisms than with the classical 'biological species concept'.

Intraspecific size variation in planktonic foraminifera cannot be consistently predicted by the environment

The size structure of plankton communities is an important determinant of their functions in marine ecosystems. However, few studies have quantified how organism size varies within species across biogeographical scales. Here, we investigate how planktonic foraminifera, a ubiquitous zooplankton group, vary in size across the tropical and subtropical oceans of the world. Using a recently digitized museum collection, we measured shell area of 3,799 individuals of nine extant species in 53 seafloor sediments. We first analyzed potential size biases in the collection. Then, for each site, we obtained corresponding local values of mean annual sea-surface temperature (SST), net primary productivity (NPP), and relative abundance of each species. Given former studies, we expected species to reach largest shell sizes under optimal environmental conditions. In contrast, we observe that species differ in how much their size variation is explained by SST, NPP, and/or relative abundance. While some species have predictable size variation given these variables (Trilobatus sacculifer, Globigerinoides conglobatus, Globigerinella siphonifera, Pulleniatina obliquiloculata, Globorotalia truncatulinoides), other species show no relationships between size and the studied covariates (Globigerinoides ruber, Neogloboquadrina dutertrei, Globorotalia menardii, Globoconella inflata). By incorporating intraspecific variation and sampling broader geographical ranges compared to previous studies, we conclude that shell size variation in planktonic foraminifera species cannot be consistently predicted by the environment. Our results caution against the general use of size as a proxy for planktonic foraminifera environmental optima. More generally, our work highlights the utility of natural history collections and the importance of studying intraspecific variation when interpreting macroecological patterns.

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.

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.

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).

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.

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.

Spatially varying selection shapes life history clines among populations of Drosophila melanogaster from sub-Saharan Africa

Clines in life history traits, presumably driven by spatially varying selection, are widespread. Major latitudinal clines have been observed, for example, in Drosophila melanogaster, an ancestrally tropical insect from Africa that has colonized temperate habitats on multiple continents. Yet, how geographic factors other than latitude, such as altitude or longitude, affect life history in this species remains poorly understood. Moreover, most previous work has been performed on derived European, American and Australian populations, but whether life history also varies predictably with geography in the ancestral Afro-tropical range has not been investigated systematically. Here, we have examined life history variation among populations of D. melanogaster from sub-Saharan Africa. Viability and reproductive diapause did not vary with geography, but body size increased with altitude, latitude and longitude. Early fecundity covaried positively with altitude and latitude, whereas lifespan showed the opposite trend. Examination of genetic variance-covariance matrices revealed geographic differentiation also in trade-off structure, and QST -FST analysis showed that life history differentiation among populations is likely shaped by selection. Together, our results suggest that geographic and/or climatic factors drive adaptive phenotypic differentiation among ancestral African populations and confirm the widely held notion that latitude and altitude represent parallel gradients.

Global biogeography and evolution of Cuvierina pteropods

BACKGROUND

Shelled pteropods are planktonic gastropods that are potentially good indicators of the effects of ocean acidification. They also have high potential for the study of zooplankton evolution because they are metazoan plankton with a good fossil record. We investigated phenotypic and genetic variation in pteropods belonging to the genus Cuvierina in relation to their biogeographic distribution across the world's oceans. We aimed to assess species boundaries and to reconstruct their evolutionary history.

RESULTS

We distinguished six morphotypes based on geometric morphometric analyses of shells from 926 museum and 113 fresh specimens. These morphotypes have distinct geographic distributions across the Atlantic, Pacific and Indian oceans, and belong to three major genetic clades based on COI and 28S DNA sequence data. Using a fossil-calibrated phylogeny, we estimated that these clades separated in the Late Oligocene and Early to Middle Miocene. We found evidence for ecological differentiation among all morphotypes based on ecological niche modelling with sea surface temperature, salinity and phytoplankton biomass as primary determinants. Across all analyses, we found highly congruent patterns of differentiation suggesting species level divergences between morphotypes. However, we also found distinct morphotypes (e.g. in the Atlantic Ocean) that were ecologically, but not genetically differentiated.

CONCLUSIONS

Given the distinct ecological and phenotypic specializations found among both described and undescribed Cuvierina taxa, they may not respond equally to future ocean changes and may not be equally sensitive to ocean acidification. Our findings support the view that ecological differentiation may be an important driving force in the speciation of zooplankton.

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.

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.