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

Single-cell transcriptomics supports presence of cryptic species and reveals low levels of population genetic diversity in two testate amoebae morphospecies with large population sizes

The enormous population sizes and wide biogeographical distribution of many microbial eukaryotes set the expectation of high levels of intraspecific genetic variation. However, studies investigating protist populations remain scarce, mostly due to limited 'omics data. Instead, most genetics studies of microeukaryotes have thus far relied on single loci, which can be misleading and do not easily allow for detection of recombination, a hallmark of sexual reproduction. Here, we analyze >40 genes from 72 single-cell transcriptomes from two morphospecies-Hyalosphenia papilio and Hyalosphenia elegans-of testate amoebae (Arcellinida, Amoebozoa) to assess genetic diversity in samples collected over four years from New England bogs. We confirm the existence of cryptic species based on our multilocus dataset, which provides evidence of recombination within and high levels of divergence between the cryptic species. At the same time, total levels of genetic diversity within cryptic species are low, suggesting that these abundant organisms have small effective population sizes, perhaps due to extinction and repopulation events coupled with efficient modes of dispersal. This study is one of the first to investigate population genetics in uncultivable heterotrophic protists using transcriptomics data and contributes towards understanding cryptic species of nonmodel microeukaryotes.

Macroevolutionary patterns in intragenomic rDNA variability among planktonic foraminifera

Ribosomal intragenomic variability in prokaryotes and eukaryotes is a genomic feature commonly studied for its inflationary impact on molecular diversity assessments. However, the evolutionary mechanisms and distribution of this phenomenon within a microbial group are rarely explored. Here, we investigate the intragenomic variability in 33 species of planktonic foraminifera, calcifying marine protists, by inspecting 2,403 partial SSU sequences obtained from single-cell clone libraries. Our analyses show that polymorphisms are common among planktonic foraminifera species, but the number of polymorphic sites significantly differs among clades. With our molecular simulations, we could assess that most of these mutations are located in paired regions that do not affect the secondary structure of the SSU fragment. Finally, by mapping the number of polymorphic sites on the phylogeny of the clades, we were able to discuss the evolution and potential sources of intragenomic variability in planktonic foraminifera, linking this trait to the distinctive nuclear and genomic dynamics of this microbial group.

Newly designed foraminifera primers identify habitat-specific lineages through metabarcoding analyses

Foraminifera include diverse shell-building lineages found in a wide array of aquatic habitats from the deep-sea to intertidal zones to brackish and freshwater ecosystems. Recent estimates of morphological and molecular foraminifera diversity have increased the knowledge of foraminiferal diversity, which is critical as these lineages are used as bioindicators of past and present environmental perturbation. However, a comparative analysis of foraminiferal biodiversity between their major habitats (freshwater, brackish, intertidal, and marine) is underexplored, particularly using molecular tools. Here, we present metabarcoding survey of foraminiferal diversity across different ecosystems using newly designed foraminifera-specific primers that target the hypervariable regions of the foraminifera SSU-rRNA gene (~250-300bp long). We tested these primer sets on four foraminifera species and then across several environments: the intertidal zone, coastal ecosystems, and freshwater vernal pools. We retrieved 655 operational taxonomic units (OTUs); the majority are undetermined taxa that have no closely-matching sequences in the database. Furthermore, we identified 163 OTUs with distinct habitat preferences. Most of the observed OTUs belonged to lineages of single-chambered foraminifera, including poorly explored freshwater foraminifera which encompass a clade of Reticulomyxa-like forms. Our pilot study provides the community with an additional set of newly designed and taxon-specific primers to elucidate foraminiferal diversity across different habitats.

Reproduction dynamics of planktonic microbial eukaryotes in the open ocean

Understanding the biology of reproduction of an organismal lineage is important for retracing key evolutionary processes, yet gaining detailed insights often poses major challenges. Planktonic Foraminifera are globally distributed marine microbial eukaryotes and important contributors to the global carbon cycle. They cannot routinely be cultured under laboratory conditions across generations, and thus details of their life cycle remain incomplete. The production of flagellated gametes has long been taken as an indication of exclusively sexual reproduction, but recent research suggests the existence of an additional asexual generation in the life cycle. To gain a better understanding of the reproductive biology of planktonic Foraminifera, we applied a dynamic, individual-based modelling approach with parameters based on laboratory and field observations to test if sexual reproduction is sufficient for maintaining viable populations. We show that temporal synchronization and potentially spatial concentration of gamete release seems inevitable for maintenance of the population under sexual reproduction. We hypothesize that sexual reproduction is likely beneficial during the adaptation to new environments, while population sustenance in stable environments can be ensured through asexual reproduction.

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.

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.

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.

Integrated species-phenon trees: visualizing infraspecific diversity within lineages

The unprecedented detail with which contemporary molecular phylogenetics are visualizing infraspecific relationships within living species and species complexes cannot as yet be reliably extended into deep time. Yet paleontological systematics has routinely dealt in (mainly) morphotaxa envisaged in various ways to have been components of past species lineages. Bridging these perspectives can only enrich both. We present a visualization tool that digitally depicts infraspecific diversity within species through deep time. Our integrated species-phenon tree merges ancestor-descendant trees for fossil morphotaxa (phena) into reconstructed phylogenies of lineages (species) by expanding the latter into "species boxes" and placing the phenon trees inside. A key programming strategy to overcome the lack of a simple overall parent-child hierarchy in the integrated tree has been the progressive population of a species-phenon relationship map which then provides the graphical footprint for the overarching species boxes. Our initial case has been limited to planktonic foraminfera via Aze & others' important macroevolutionary dataset. The tool could potentially be appropriated for other organisms, to detail other kinds of infraspecific granularity within lineages, or more generally to visualize two nested but loosely coupled trees.

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 progress? A phylogeographic study among populations of Hemitrichia serpula (Myxomycetes)

Myxomycetes (plasmodial slime molds, Amoebozoa) are often perceived as widely distributed, confounding to the “everything is everywhere” hypothesis. To test if gene flow within these spore-dispersed protists is restricted by geographical barriers, we chose the widespread but morphologically unmistakable species Hemitrichia serpula for a phylogeographic study. Partial sequences from nuclear ribosomal RNA genes (SSU) revealed 40 ribotypes among 135 specimens, belonging to three major clades. Each clade is dominated by specimens from a certain region and by one of two morphological varieties which can be differentiated by SEM micrographs. Partial sequences of the protein elongation factor 1 alpha (EF1A) showed each clade to possess a unique combination of SSU and EF1A genotypes. This pattern is best explained assuming the existence of several putative biospecies dominating in a particular geographical region. However, occasional mismatches between molecular data and morphological characters, but as well heterogeneous SSU and heterozygous EF1A sequences, point to ongoing speciation. Environmental niche models suggest that the putative biospecies are rather restricted by geographical barriers than by macroecological conditions. Like other protists, myxomycetes seem to follow the moderate endemicity hypothesis and are in active speciation, which is most likely shaped by limited gene flow and reproductive isolation.

Environmental Predictors of Diversity in Recent Planktonic Foraminifera as Recorded in Marine Sediments

Global diversity patterns are thought to result from a combination of environmental and historical factors. This study tests the set of ecological and evolutionary hypotheses proposed to explain the global variation in present-day coretop diversity in the macroperforate planktonic foraminifera, a clade with an exceptional fossil record. Within this group, marine surface sediment assemblages are thought to represent an accurate, although centennial to millennial time-averaged, representation of recent diversity patterns. Environmental variables chosen to capture ocean temperature, structure, productivity and seasonality were used to model a range of diversity measures across the world's oceans. Spatial autoregressive models showed that the same broad suite of environmental variables were important in shaping each of the four largely independent diversity measures (rarefied species richness, Simpson's evenness, functional richness and mean evolutionary age). Sea-surface temperature explains the largest portion of diversity in all four diversity measures, but not in the way predicted by the metabolic theory of ecology. Vertical structure could be linked to increased diversity through the strength of stratification, but not through the depth of the mixed layer. There is limited evidence that seasonal turnover explains diversity patterns. There is evidence for functional redundancy in the low-latitude sites. The evolutionary mechanism of deep-time stability finds mixed support whilst there is relatively little evidence for an out-of-the-tropics model. These results suggest the diversity patterns of planktonic foraminifera cannot be explained by any one environmental variable or proposed mechanism, but instead reflect multiple processes acting in concert.

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.

Fossil and genetic evidence for the polyphyletic nature of the planktonic foraminifera "Globigerinoides", and description of the new genus Trilobatus

Planktonic foraminifera are one of the most abundant and diverse protists in the oceans. Their utility as paleo proxies requires rigorous taxonomy and comparison with living and genetically related counterparts. We merge genetic and fossil evidence of "Globigerinoides", characterized by supplementary apertures on spiral side, in a new approach to trace their "total evidence phylogeny" since their first appearance in the latest Paleogene. Combined fossil and molecular genetic data indicate that this genus, as traditionally understood, is polyphyletic. Both datasets indicate the existence of two distinct lineages that evolved independently. One group includes "Globigerinoides" trilobus and its descendants, the extant "Globigerinoides" sacculifer, Orbulina universa and Sphaeroidinella dehiscens. The second group includes the Globigerinoides ruber clade with the extant G. conglobatus and G. elongatus and ancestors. In molecular phylogenies, the trilobus group is not the sister taxon of the ruber group. The ruber group clusters consistently together with the modern Globoturborotalita rubescens as a sister taxon. The re-analysis of the fossil record indicates that the first "Globigerinoides" in the late Oligocene are ancestral to the trilobus group, whereas the ruber group first appeared at the base of the Miocene with representatives distinct from the trilobus group. Therefore, polyphyly of the genus "Globigerinoides" as currently defined can only be avoided either by broadening the genus concept to include G. rubescens and a large number of fossil species without supplementary apertures, or if the trilobus group is assigned to a separate genus. Since the former is not feasible due to the lack of a clear diagnosis for such a broad genus, we erect a new genus Trilobatus for the trilobus group (type species Globigerina triloba Reuss) and amend Globoturborotalita and Globigerinoides to clarify morphology and wall textures of these genera. In the new concept, Trilobatus n. gen. is paraphyletic and gave rise to the Praeorbulina/Orbulina and Sphaeroidinellopsis/Sphaeroidinella lineages.

The Interdependency of the Morphological Variations of the Planktonic Foraminiferal Species Globigerina bulloides in Surface Sediments on the Environmental Parameters of the Southwestern Indian Ocean

18 surface sediment samples collected from a north-south transect along the Indian Ocean have been analyzed for planktonic Foraminifera content. Among the other planktonic foraminiferal faunas, Globigerina bulloides was present substantially in all samples. Census data of G. bulloides were measured for different parameters (average size, mean proloculus size, coiling direction, and number of chambers) and a Q-mode cluster analysis was applied on these data. Samples were segregated into two homogeneous clusters, each reflecting particular environmental conditions. Two clusters are as follows: (1) Cluster A, comprised of 6 samples and characterized by the highest range of foraminiferal and ecological parameters, except sea surface temperature and salinity which shows the lowest range, and (2) Cluster B, comprised of 12 samples and characterized by the lowest range of foraminiferal parameters and ecological parameters, except sea surface temperature and salinity which shows the highest range. The study suggests that the ecological parameters are the governing factors for the morphological characteristics of planktonic foraminiferal species G. bulloides.