Genetic and morphological divergence in the warm-water planktonic foraminifera genus Globigerinoides
Authors:
Raphaël Morard aff001; Angelina Füllberg aff001; Geert-Jan A. Brummer aff002; Mattia Greco aff001; Lukas Jonkers aff001; André Wizemann aff004; Agnes K. M. Weiner aff001; Kate Darling aff006; Michael Siccha aff001; Ronan Ledevin aff008; Hiroshi Kitazato aff009; Thibault de Garidel-Thoron aff010; Colomban de Vargas aff011; Michal Kucera aff001
Authors place of work:
MARUM Center for Marine Environmental Sciences, University of Bremen, Leobener Strasse, Bremen, Germany
aff001; NIOZ Royal Netherlands Institute for Sea Research, Department of Ocean Systems, and Utrecht University, Den Burg, and Utrecht University, The Netherlands
aff002; Vrije Universiteit Amsterdam, Department of Earth Sciences, Faculty of Science, Amsterdam, The Netherlands
aff003; Leibniz Centre for Tropical Marine Research, Bremen, Germany
aff004; Department of Biological Sciences, Smith College, Northampton, Massachusetts, United States of America
aff005; School of GeoSciences, University of Edinburgh, Edinburgh, Scotland, United Kingdom
aff006; School of Geography and Sustainable Development, University of St Andrews, St Andrews, Scotland, United Kingdom
aff007; UMR5199 PACEA, Université de Bordeaux, Allée Geoffroy Saint Hilaire, Pessac, France
aff008; Japan Agency for Marine Earth Science and Technology (JAMSTEC), Yokosuka, Kanagawa, Japan
aff009; Aix-Marseille Université, CNRS, IRD, Collège de France, INRA, CEREGE, Aix-en-Provence, France
aff010; Sorbonne Université, CNRS, Station Biologique de Roscoff, UMR 7144, ECOMAP, Roscoff, France
aff011; Research Federation for the Study of Global Ocean Systems Ecology and Evolution, FR2022/Tara GOSEE, Paris, France
aff012
Published in the journal:
PLoS ONE 14(12)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0225246
Summary
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.
Keywords:
Phylogenetics – Phylogenetic analysis – Taxonomy – Phylogeography – Plankton – Morphogenesis – Paleogenetics – Cryptic speciation
Introduction
Species of the genus Globigerinoides are the dominant constituent of tropical-subtropical planktonic foraminifera assemblages throughout the Neogene and represent a cornerstone for paleoceanography. The extant members of the genus feature one of the most iconic species of planktonic foraminifera that was formally described from the Atlantic by d’Orbigny in 1839 as Globigerina rubra, after the reddish coloration of its test. The species definition was later widened to include colorless specimens as variants with the same morphology because shell color was not considered taxonomically relevant at the species level [1,2]. It was further broadened by Parker [3] to include the morphologically similar Globigerinoides elongatus (d’Orbigny) and Globigerinoides pyramidalis (van den Broeck) that were originally distinguished using characteristics such as the compression of the last chamber and a higher trochospire. Parker [3] considered that the three species formed a morphological continuum with G. ruber and this broad definition was endorsed by Kennett and Srinivasan in 1983 [4], who interpreted G. elongatus, G. pyramidalis and also G. cyclostomus (Galloway and Wissler) as ecophenotypic variants of G. ruber. This broad species definition has remained stable since, but most researchers continued to distinguish the two “chromotypes” as G. ruber “white” and G. ruber “pink”, because of differences in biogeography, seasonality and isotopic composition [5]. Their distinction is particularly highlighted by the extinction of G. ruber “pink” in the Indian and Pacific Oceans 120,000 years ago, while persisting in the Atlantic to the present day [6].
The lumping of G. elongatus, G. pyramidalis, G. cyclostomus with G. ruber was questioned by Robbins and Healy-Williams [7], who identified stable isotopic differences among morphological variants. This motivated Wang [8] to further test for isotopic differences between morphological variants of G. ruber “white”. Wang [8] informally re-created the split between G. ruber and G. elongatus, that had already been identified by d’Orbigny and referred to the original G. ruber as G. ruber sensu stricto (s.s.) and lumped the specimens matching the description of G. elongatus, G. pyramidalis and G. cyclostomus into G. ruber sensu lato (s.l.). Wang [8] showed subtle but statistically significant differences of 0.21 ± 0.21‰ for δ18O and −0.28±0.29‰ for δ13C between the two informal taxonomic units in the South China Sea and suggested that G. ruber s.s. lived in the upper 30 meters of the water column and G. ruber s.l. lived below 30 meters. Wang used this feature to reconstruct the variation of the thermal structure of the water column during the last glacial cycle. The work of Wang [8] triggered a series of studies during the last two decades that examined chemical/compositional, morphological and ecological differences between G. ruber s.s. and G. ruber s.l. [8–21] to assess their usefulness for paleoceanography.
In parallel to the investigation of the ecology of G. ruber s.l. and s.s., sequencing of the small sub-unit of the ribosomal RNA gene (SSU rDNA) shed new light on the diversity within the genus Globigerinoides. The earliest molecular phylogenies by Darling et al. [22,23] demonstrated that the two chromotypes of G. ruber are genetically distinct, in line with the well-established biogeographical and ecological differences [5]. Later, Darling and Wade [24] described further genetic diversity within G. ruber “white” and Kuroyanagi et al. [25] suggested that the genetic discontinuity observed within G. ruber “white” mirrored the sensu stricto/sensu lato division of Wang [8]. These observations were confirmed by Aurahs et al. [26] who identified four genotypes in G. ruber “white” (Ia, Ib, IIa and IIb) and in a second study [27], these authors analysed images of the barcoded specimens to show that genotypes Ia and Ib matched the diagnosis of G. ruber s.s., whilst genotype IIa matched the diagnosis of G. ruber s.l. As a result, they proposed to reinstate Globigerinoides elongatus as a valid name for genotype IIa. Irrespective of the complicated taxonomy, all genetic studies consistently identified G. elongatus (or G. ruber s.l.) as a sister to the morphologically distinct species G. conglobatus. The contrast between the genetic divergence and morphological similarity of G. elongatus and G. ruber implies a disconnection between genetic and morphological evolution in the genus. Thus, next to the need to clarify and stabilize its nomenclature, the complex diversification pattern in the genus also calls for a comprehensive study of the pattern of speciation and morphological diversification leading to the present-day diversity in Globigerinoides.
To this end, we assembled a global dataset of single-cell SSU rDNA sequences covering all morphospecies of the genus, applied an objective molecular nomenclature system [28] to parse the genetic variability and used the shell morphology of the barcoded specimens to map the genetic units onto a morphological taxonomic framework. To explore patterns of morphological evolution within the genus, we used CT scanning to quantify the ontogenetic trajectory of the five morphospecies [29,30]. This allowed us to investigate whether the diagnostic differences in adult morphology between closely related species in the genus could be the result of heterochrony, with slight alteration in the developmental sequence leading to large differences in adult shape and size. Finally, we use our collection of globally distributed samples to analyze the ecology of the morphological and cryptic species in the genus and discuss the potential drivers of their evolution.
Material
Living planktonic foraminifera of the morphospecies Globoturborotalita rubescens, Globigerinoides ruber, Globigerinoides conglobatus, Globigerinoides elongatus and Globigerinoides tenellus were sampled between 1993 and 2015 during 23 research cruises and 6 near shore sampling campaigns (Fig 1) in all oceans. No sampling permit was needed for planktonic foraminifera. G. rubescens was included in the analysis to serve as outgroup in phylogenetic analyses. The specimens were sampled using different open-closing plankton net systems, simple plankton nets or ship pump systems between 0 and 700 m water depth and mesh sizes from 63 to 200 μm. The specimens were separated from other plankton, cleaned with brushes and either transferred onto cardboard slides and air-dried or directly transferred into DNA extraction buffer and stored at -20°C or -80°C. The specimens stored on cardboard slides were transferred into DNA extraction buffer later in the laboratory.
Methods
Molecular analyses
DNA extraction was performed using either the DOC protocol, the GITC* protocol or the Urea Protocol [31]. A fragment located at the 3’end of the of the SSU rDNA between the primers S14F1 or S14p and 1528R [32] was amplified and the PCR products obtained were purified and sequenced directly with Sanger sequencing by several service providers (LGC Genomics Berlin, University of Edinburgh Gene Pool, AGOWA and Station Biologique de Roscoff). In addition, we randomly selected eight specimens for cloning in order to quantify potential intragenomic variability and used the TOPO TA cloning kit (Invitrogen) according to manufacturer instructions. Between 2 and 13 clones were sequenced per individual. All chromatograms were carefully checked to ensure sequence quality and were deposited on NCBI under the accession numbers MN383323 to MN384218. The methodologies used for sampling, DNA extraction, amplification and cloning of single planktonic foraminifera cells are described in Weiner et al. [31].
Public databases
We completed our dataset with sequences already made available by earlier studies. First, we retrieved all 359 SSU rDNA sequences of the six morphospecies that were stored in the PFR2 database v 1.0 [33]. We then manually queried the NCBI portal (last accession: 15.11.2018) and retrieved seven additional sequences of G. ruber (Accession numbers KY397454-KY397460).
Detailed information on handling procedures, sequences and associated metadata of the newly generated data and those retrieved from public databases are provided in S1 Table.
Molecular nomenclature
The genetic diversity within the six morphospecies was classified into a three-tier hierarchical scheme of Molecular Taxonomic Units following the system described in Morard et al [28]. The system uses the amplified ~1000 bp long sequence fragment located at the 3’end of the SSU rDNA between stems 32 and 50 as molecular marker [32], which is the barcode selected for benthic foraminifera [34] that covers six variable regions, three of which are foraminifera-specific. To exclude potential sequencing errors when constructing the nomenclature, we retained only sequences for which the individual sequence pattern was observed at least three times across our dataset. All distinct sequences in the resulting trimmed dataset were considered as basetypes. Basetypes co-occurring within one or several individuals (because of intra-individual variability among tandem copies of the gene) were assembled into basegroups, and constitute the lowest level of the nomenclature (MOTUs lvl-3). The variability observed between the basetypes represents at least the intragenomic (intra-individual) variability and the variability observed among different basegroups is considered to represent at least the level of population variability. If a unique basetype is observed within a single specimen, which is the majority of cases in our dataset (see Results), the resulting basegroup contains a single basetype. The levels 1 and 2 of the nomenclature (following Morard et al. [28]) were constructed using a combination of two automated delimitation methods, the Automated Barcode Gap Discovery method (ABGD; [35]) and the Poisson Tree Process (PTP; [36]). The sequences were aligned with MAFFT v.7 [37] and a phylogenetic inference was calculated with 1000 non-parametric bootstrapping pseudo replicates based on a BioNJ starting tree using PhyML [38]. The best substitution models were selected using the Smart Model Selection [39] under Akaike Information Criterion and the model GTR+I+G was selected. The resulting trees were submitted to the PTP server (http://species.h-its.org/) under default settings. The same alignment that served to generate the tree was submitted to the online ABGD server (http://wwwabi.snv.jussieu.fr/public/abgd/abgdweb.html) using the Kimura K80 distance and default options. We retained the initial (coarsest delimitation) and recursive partition (finest delimitation) provided with the lowest prior intraspecific divergence. We defined the MOTU lvl-2 as the finest delimitation proposed by either ABGD or PTP and the MOTU lvl-1 as the coarsest. The proposed delimitations are retained as working hypotheses provided that two clones belonging to the same basegroup were not attributed to different partitions (oversplit) and that sequences belonging to different morphospecies were not grouped in the same partition (lumping). The delimitation proposed by ABGD and PTP as well as the retained delimitation are reported in Fig 2. As multiple, but partly overlapping, nomenclatural schemes were proposed by successive studies [21–24, 40, 41], we reported the correspondence between these schemes and their equivalent in our system (Fig 3 and S2 Table).
A significant part of the sequences had insufficient quality and/or coverage to be included in the assessment of the diversity within the Globigerinoides plexus, but carried enough information to be attributed to at least one MOTU level of our nomenclatural system. The Sanger sequences not meeting the quality criteria were compared to the basetype sequences and received the finest taxonomic attribution possible based on the availability of diagnostic sites in the region they covered (See S1 Table). Biogeography and temporal occurrences of the genotypes and basegroups are shown in Fig 4.
Sample coverage and environmental parameters
We calculated rarefaction curves at MOTUs lvl-2 and lvl-3 (Fig 5) and complemented the approach with a first order Jackknifing to evaluate the coverage of our dataset (Table 1). Because G. tenellus and G. conglobatus were under-sampled (60 sequences in a dataset of 1251 sequences), we calculated the rarefaction curves to include all species and selectively only for G. ruber and G. elongatus separately (Fig 5). Likewise, the Jackknifing was applied to G. ruber and G. elongatus combined at the MOTUs lvl-2 and on each species separately at the lvl-3 (Table 1). We then applied the analyses to the global dataset and separately on three main biogeographic regions: North Atlantic Ocean, Indian Ocean and Pacific Ocean.
The dataset constituted for this study is the result of the efforts by multiple research teams and re-exploitation of public data, therefore it was difficult to recover and harmonize the environmental parameters measured during each sampling campaign. In order to analyze the ecological preferences of the sampled genotypes and basegroups, we chose to use geographic coordinates and collection date to extract the monthly average values of the following environmental parameters from public databases: Sea Surface Temperature (SST), Mixed Layer Salinity (MLS), Chlorophyll concentration (CHL), Particulate Organic Carbon (POC) and Productivity (PROD). The SST, CHL and POC parameters were extracted from the MODIS-Aqua (NASA, Greenbelt, MD, USA) database [42–44], the MLS was extracted from the Isopycnal/Mixed-layer Ocean Climatology (MIMOC) database [45] and PROD was calculated following the Vertically Generalized Production Model from Behrenfeld and Falkowski [46]. In this way, we could gather a homogeneous environmental dataset although it is less precise than in-situ measurements. We display the environmental parameter values at the morphospecies, genotype and basegroup levels in Fig 6 and tested if the distribution of values of sister taxa at each taxonomic level was the same (null hypothesis) with a simple non-parametric Wilcoxon-Mann-Whitney U-test using the Bonferroni correction (Table 2). All statistical analyses were performed in PAST 3.21 [47].
Phylogeny and molecular clock
To reconstruct the evolutionary history of the genus Globigerinoides, we applied a molecular clock estimation using the same alignment as for the maximum likelihood tree inference (Fig 2). We used the divergence between G. rubescens and the genus Globigerinoides (23.8 Ma [48]), the First Appearance Datum (FAD) of G. conglobatus (8–8.6 Ma) and G. tenellus (2.5 Ma), which are known from the fossil record [49], as minimum ages to constrain the phylogeny. We used a relaxed clock model implemented in BEAST v.1.8.4 [50]. Model parameters were set using BEAUti v1.8.4. The distribution of the fixed node age prior was considered normal and the speciation rate was assumed constant under the Yule-Process. The GTR (Generalised Time Reversible) model was selected as substitution model and an UPGMA (Unweighted Pair Group method with arithmetic mean) tree was calculated as starting tree. Markov-Chain-Monte Carlo (MCMC) analyses were conducted for 10,000,000 generations, with a burn-in of 1000 generations and saving each 1000th generation. The maximum clade credibility tree with median node heights was calculated in TREEAnnotator from the BEAST package, with a burn-in of 100 trees and a posterior probability limit of 0. The resulting tree was then visualized in FigTree v. 1.3.1 [51] and is shown in Fig 7.
3D morphology
We produced CT-scans of G. rubescens, G. ruber albus n.subsp., G. conglobatus, G. elongatus and G. tenellus to assess the ontogenetic development of each species. To ensure that the specimens had completed their life cycle, which usually is not the case for the living specimens collected in the water column, we used specimens recently deposited on the seafloor from a core top sample retrieved south of Barbados at station GeoB3935 (12°36.8 N, 59°23.2 W; bottom depth 1554 meters) [52]. We chose this sample because of the exceptional preservation of the tests, which were free of fine-grained sediment. Moreover, its provenance is close to the sampling localitions where Globigerinoides spp. were previously analysed for their ontogeny [53]. From this sample, we selected one specimen per morphospecies, choosing specimens with well-developed characteristic features. We choose specimens that were of large size and had a thick test, which indicates maturity and facilitates CT-scanning at good resolution. Indeed, four of the five species have a diminutive final chamber indicative of reproduction by gametogenesis (the reproductive terminal stage sensu Brummer et al. [53]), while G. ruber is normalform. We realize that planktonic foraminifera are morphologically variable, not only in their adult shape but also throughout their ontogeny [53], so the decision to analyze only a single specimen per morphospecies was made in order to achieve a first rough assessment of the main differences of ontogenetic trajectories among the morphospecies. Such trajectories are known to differ between species but are stable within species, with much variability correlated with proloculus size [29,30,53]. The selected specimens were individually mounted on a stub and scanned at a cubic resolution of 1.2 μm with a General Electrics V/Tome/x micro-scanner (PACEA, Bordeaux University). Each scan was performed at 80 kV and 180 μA without filter as the shell had a low X-ray absorption rate. The smaller specimen of G. rubescens was analyzed with a cubic resolution of 0.68 μm with a Zeiss Versa 500 at 80kV, 7W and with a filter LE1. Semi-automated segmentation was used to reconstruct three-dimensional (3D) virtual surfaces of the calcite volume (external morphology) of each specimen (Fig 8), and the inner volume of individual chambers (Fig 9) were produced by manual segmentation with the ITK-SNAP v 3.6 software [54] to reconstruct the ontogenetic trajectory of each morphospecies. We automatically extracted the volume, centroid position and major axis of individual chambers using a custom script in MATLAB R2017b to calculate growth parameters of the trochospire, following the model of Raup [55]. We calculated the whorl expansion rate W, the relative distance between the generating curve and the axis of coiling D, the translation rate T and the shape of the generating curve S. The calculated growth parameters of each species are displayed in Fig 10 and the numerical values are provided in S3 Table.
Results
Genetic diversity within Globigerinoides
Our dataset on the molecular diversity within the genus Globigerinoides and its sister species G. rubescens includes 1251 Sanger sequences, of which 893 are new. All 1251 sequences cover the same rDNA barcode region and originated from a total of 1159 individuals collected at 179 sampling stations (Fig 1). Among the 1251 sequences, 147 met the quality criteria to derive molecular taxonomy and served to define a total of 17 basetypes (unique, replicable sequence motifs). We observed three basetypes that co-occurred within two single individuals of G. rubescens that were consequently grouped into a single basegroup. Additionally, we identified the co-occurrence of two basetypes within three clones from a single individual of G. ruber, published by Kuroyanagi et al. [25]. Since this is the only observation of intragenomic variability within the SSU rRNA gene in G. ruber, we consider it likely that it resulted from contamination or PCR/sequencing error and we thus reject this single observation as evidence for intragenomic variability in the species. As a result, we retained 15 basegroups (Fig 2), 14 of these consisting of a single basetype, which provided a basis for the construction of a molecular nomenclature of the group. The automated taxa partitions proposed by ABGD and PTP did not violate any of the conditions of the taxonomic system (lumping of sequences belonging to different morphotaxa or splitting of basetypes belonging to the same basegroup) and were thus retained. Partitions by ABGD reflected the morphological species concept of the group. The PTP analysis identified three partitions within G. ruber albus n.subsp. and two within G. conglobatus, which were retained as distinct genotypes. No partitions were identified within the morphospecies G. rubescens, G. elongatus, G. tenellus and G. ruber ruber indicating that these morphospecies consist of only a single genotype.
Molecular and morphological revision of existing taxonomic concepts
The first DNA sequences of members of the genus Globigerinoides were made available in the earliest publications on the genetic diversity of planktonic foraminifera [23,56–58], but nomenclatural schemes to describe the cryptic diversity in the genus were presented only a decade later in parallel and independently by Darling and Wade [24] and Kuroyanagi et al. [25], who both identified five cryptic species within G. ruber (Fig 3). The complexity of naming cryptic species further increased in the following year when Ujiié and Lipps [40] produced a distinct nomenclatural scheme with only four cryptic species within G. ruber, whereas Aurahs et al. [26] further developed the scheme initially proposed by Darling and Wade [24], but chose to interpret all subtle sequence differences across their dataset and produced a scheme with 14 different cryptic species. Two years later, Aurahs et al. [27] reduced the diversity to only eight cryptic species by considering only the most repeatable sequence pattern in their dataset. Furthermore, they split the genetic diversity between G. ruber s. s. (G. ruber Ia, Ib, Ib2 and pink) and G. ruber s. l. (G. ruber IIa, IIa1, IIa2 and IIb) with the genotype IIa being considered as G. elongatus. André et al. [59] proposed a revision of the available nomenclature based on automated methods for species delimitation to define cryptic diversity in planktonic foraminifera, which reduced the diversity to five cryptic species only.
In our study, we quadrupled the size of the dataset compared to previous studies and placed all the formally described cryptic species into a new framework, and we identified only one new basegroup in G. ruber albus n.subsp. (Ib2). We have also generated a SSU rDNA sequence from a specimen identified on collection and by later observations as G. tenellus (S1 Fig) that was identical to the sequences obtained from specimens of G. ruber Type IIb of Aurahs et al. [27]. This allowed us to recognize this type as G. tenellus and thus return the species to Globigerinoides, as a sister to G. elongatus. The extended and strictly curated dataset allowed for identifying one new basegroup in G. conglobatus as well as in G. tenellus and reducing the number of basetypes to three within G. elongatus.
Since the genetic distance separating the “pink” and “white” chromospecies of G. ruber is greater than the distance separating G. elongatus and G. tenellus (Fig 2), we feel compelled to express the genetic and phenotypic distinction between the two lineages in formal taxonomy. However, the phenotypic distinction reflects the color of the shell, not shell morphology, and this character fades with age, rendering it impossible to distinguish the lineages in fossil material older than ~750 kyr [6]. Therefore, we propose to use the subspecies names G. ruber albus n.subsp. and G. ruber ruber, facilitating continuity by allowing the use of the nominotype “G. ruber” at the species level in situations where the chromospecies cannot be differentiated. Moreover, the physical holotype designated here combines the shell morphology with the SSU rDNA drawn from the same individual, the first for a planktonic foraminifer. The physical specimens have been deposited at the Naturalis Biodiversity Center, Leiden, the Netherlands.
Systematics
Phylum Foraminifera d’Orbigny, 1826
Class Globothalamea Pawlowski, Holzmann & Tyszka, 2013
Order Rotaliida Delage and Hérouard, 1896
Superfamily Globigerinoidea Carpenter, Parker & Jones, 1862
Family Globigerinidae Carpenter, Parker & Jones, 1862
Genus Globigerinoides Cushman, 1927, amended by Spezzaferri et al., 2015
Type species Globigerina rubra d’Orbigny, 1839
Species Globigerinoides ruber (d’Orbigny, 1839)
Subspecies Globigerinoides ruber albus n. subsp.
Type material: Holotype: Voucher C319 collected at 7.409°S, 165.274°E on 12.03.2013 between 0–20 meters water depth (Museum number: RGM.1332320). Paratypes: Voucher C208 collected at 6.414°N, 143.024°E on 18.03.2013 between 80–100 meters water depth (Museum Number: RGM.1332321), Voucher C281 collected at 22.719°S, 170.918°E on 08.03.2013 between 60–80 meters water depth (Museum Number: RGM.1332322) and Vouchers C329 collected at 7.409°S, 165.274°E on 12.03.2013 between 0–20 meters water depth (Museum number: RGM.133233). Light microscopy images of the type specimens are provided in S2 Fig.
Diagnosis: Differs from G. ruber ruber by the absence of reddish color of the shell, by the presence of a distinct sequence motive in the SSU rDNA gene, by its seasonality and depth habitat in the modern Atlantic and its presence in the Indopacific throughout the last 120 ka. The two subspecies cannot be distinguished prior to 750 ka due to the fading of the color with time and both are then captured as G. ruber well into the Neogene.
Description. The new subspecies largely overlaps with G. ruber ruber in test morphology, but differs in the color of the test, which develops during the neanic stage [53]. The morphology of the species and its changes during the ontogeny have been described in detail by Brummer et al. [53] and is formalized accordingly below. The holotype has been selected such that the test shows all key features of the species, but lacks color and because it yielded a SSU rDNA sequence of genetic type G. ruber albus n.subsp. Ia (Voucher C319).
Prolocular stage. Proloculus small, 12.5 ± 1.5 μm (10–16 μm), wall imperforate, smooth and non-spinose; aperture interiomarginal, circular with thickened rim, in multi-chambered tests larger than deuteroconch and truncated by flat wall shared with deuteroconch.
Juvenile stage. Starting with deuteroconch, test lobate, umbilico-convex, umbilicus open, wide, narrowing after completion of initial whorl; chambers hemispherical, 7–12 (9.3 ± 1.2) added in ± 1.5 whorls of near planispire, with 5–6 in initial whorl, totaling 8–13 (9.7 ± 1.2) chambers in tests 54–76 (65.3 ± 5.8) μm in diameter. Aperture interiomarginal-marginal, a small, low arch with marked rim. Spines sparse, thin, flexible; microspines present; pores sparse, exclusively along sutures on spiral side; wall texture spinose, non-cancellate. No preferential shell coiling direction; algal symbionts acquired.
Neanic stage. Test rapidly changing towards adult morphology, becoming sphaeroidal with umbilicus closing; chambers globose, 3–4 in half to complete whorl of low trochospire, decreasing to 3 in last whorl, totaling 12–16 (14 ± 1.3) chambers in tests 120–190 (140 ± 25) μm in diameter. Aperture widening to a wide, high arch and migrating to the umbilicus. Spines and pores becoming numerous and evenly distributed; spines becoming thicker and more rigid; spine bases, inter-spine ridges and pore pits develop; wall becoming coarsely perforate and cancellate.
Adult stage. Test sphaeroidal to elongate with reddish color, chambers globose in a low-medium trochospire, at least 1, usually 2 to 3, up to 4 chambers are added, totaling 14–18 chambers in test >180, up to 510 μm in diameter, until reproduction (gametogenesis). Secondary aperture(s) develop. Wall texture cancellate-spinose and macroperforate.
Terminal stage. Usually one, occasionally two normalform and/or diminutive (kummerform) chambers are added, rarely one or two bullate chambers capping the secondary apertures. Spines progressively shed, wall coarsely perforate, smooth to coarsely cancellate. Loss of algal symbionts, loss of buoyancy. Terminal shells 230–560 μm in diameter with 15–19 chambers in 3–4 whorls of low to medium trochospire.
Distribution and ecological preferences of Globigerinoides MOTUs
Although our study benefits from a globally distributed sampling, we unfortunately lack sampling points in the Southern Atlantic. The rarefaction curves, however, confirm that the genotype diversity within Globigerinoides likely has been entirely captured by our global dataset as well as in the individual ocean basins when considering all morphospecies and the better sampled G. ruber ruber, G. ruber albus n.subsp. and G. elongatus respectively (Fig 5). We are confident that all existing genotypes and the majority of basegroups have been detected, so that we are able to interpret their biogeographic patterns (Fig 4A). We observe that the genotypes G. ruber albus n.subsp. Ia and Ib are cosmopolitan whilst the genotype G. ruber albus n.subsp. Ic was not found in the North Atlantic. A similar pattern could hold for the basegroup G. ruber albus n.subsp. Ib2 as well, as it has not been found in the North Atlantic. This may be a sampling bias because its genotype has been encountered only at two stations in the Caribbean Sea. Also, G. elongatus basegroups Ia1 and Ia3 have a cosmopolitan distribution whilst basegroup Ia2 was not found in the North Atlantic. The unique basetype detected in G. ruber ruber Ia1 was only found in the North Atlantic in our dataset. Unfortunately, the biogeography of the MOTUs of G. conglobatus and G. tenellus remains unknown due to the low number of observations.
While saturation is also reached at the basegroup level in the global dataset for the three morphospecies, it is not reached for the Indian and North Atlantic oceans, indicating that our sampling was not sufficient in these two basins. Jackknifing analysis indicates that it is likely that two basegroups of G. ruber have not been sampled in the North Atlantic, while it is possible that one basegroup of G. elongatus may still be discovered in the Pacific Ocean. However, this seems unlikely for G. elongatus because the diversity in the Indian Ocean would thus be higher (three observed and four estimated genotypes) than in the global dataset (three observed and estimated genotypes). These results may be the consequence of our unevenly distributed sampling and the fact that the detection of basegroups depends on the fragment of SSU rDNA covered, which depends of the primer used in each study. Therefore, it is impossible to say whether we failed to capture the diversity in G. ruber or G. elongatus in every basin, also given the lack of data from the South Atlantic, or if these results reflect an existing bias in our sample set.
We observe a significant difference in the sea surface temperature and mixed layer salinity at which G. ruber albus n.subsp., G. ruber ruber and G. elongatus were collected (Fig 6, Table 2). However, the apparent preference of G. ruber ruber for higher salinity may be artificial because most of our sampling for this species originates from the Caribbean and Mediterranean Seas (characterized by higher salinity) and the central Atlantic has not been sampled yet precluding a robust assessment of the true preferences of this taxa. Our sampling suggests differences between the basegroups G. elongatus Ia1 and Ia2, which occupy the lower and upper end of the thermal range of the morphological species. We also find G. ruber albus n.subsp., G. ruber albus n.subsp. Ib and G. elongatus Ia2 in more productive waters compared to G. ruber ruber, G. ruber albus n.subsp. Ia and G. elongatus Ia1, but do not observe differences with respect to chlorophyll content or particulate organic carbon. Our dataset does not reveal any seasonality in the occurrence of either the genotypes or basegroups (Fig 4B), but we stress that the sample set may not be suited to reveal such patterns.
Phylogeny of Globigerinoides
The topology and timing of diversification between members of the genus Globigerinoides (Figs 2 and 7) is largely congruent with the phylogeny proposed by Aurahs et al. [27]. The deepest split in the molecular clock phylogeny (Fig 7) separates G. ruber from G. conglobatus, G. elongatus and G. tenellus and is dated at 17.59 Ma but with a large credible interval on the age of the split (23.25 to 12.58 Ma). The Maximum-likelihood inference (Fig 2) does not support the monophyly of this clade and it is not possible to conclude from the molecular perspective alone if G. conglobatus is more closely related to G. elongatus and G. tenellus or to the G. ruber clade. The next diversification event in each lineage occurred in the late Miocene, when G. conglobatus diverged from the ancestor of G. elongatus and G. tenellus (ca. 8.29 Ma) and G. ruber albus n.subsp. and G. ruber ruber separated (ca. 6.74 Ma). Further diversification occurred between the late Pliocene and early Quaternary, when G. elongatus and G. tenellus separated concomitantly with the deepest split among the constitutive genotypes of G. ruber albus n.subsp. and G. conglobatus. A further divergence occurred in the course of the Quaternary between the genotypes Ib and Ic of G. ruber albus n.subsp., but all the remaining six divergences at the level of basetypes emerged into the Pleistocene, estimated between ~9 and 224 ka.
3D ontogenetic morphology
The largest shell diameter of the analyzed specimens ranges from 250 μm in G. rubescens and G. tenellus, to 700 μm for G. conglobatus (Fig 8), and the CT scans revealed that the specimens consist of 15 to 18 chambers (Fig 9A). The number of chambers is not fixed within a species and specimens with smaller proloculus seem to have more chambers [53]. For example, the chamber number can vary from 15 to 19 chambers in G. ruber albus, and the onset of the ontogenetic stage is not tied to the development of a particular chamber [53]. In this study, we use the chamber number as a descriptive term for convenience to explore only our results, and do not mean to imply a fixed boundary between the ontogenetic stages. In all five morphospecies the proloculus is consistently larger than the deuteroconch. Proloculus diameters differ among species, ranging from 9 μm in G. elongatus to 17 μm in G. conglobatus (Fig 10) and the ontogenetic development is accompanied by marked differences in the pattern of chamber addition among the species (Figs 9 and 10).
The ontogenetic trajectory of G. rubescens is the most stable. It begins with a steady logarithmic increase of chamber size from chambers two to thirteen, then levels off towards chamber 16 and ends with a diminutive final chamber 17 after 3.5 whorls (Figs 9 and 10). While the chamber shape (S) remains the same throughout its ontogeny, the whorl expansion rate (W) first drops steeply to chamber 5, then decreases slowly until chamber 13, slightly increases until chamber 15 and then decreases again to the final chamber. Inversely, the translation rate (T) increases slowly until chamber 13, then drops until chamber 15, to rise sharply over the two last chambers, while the relative distance between the coiling axis and the chamber centroid (D) decreases steadily throughout the ontogeny except in chambers 14–15.
Ontogenetic trajectories of G. tenellus and G. elongatus are initially similar and only diverge in the last stages. The analyzed specimen of G. tenellus produced slightly larger chambers but terminated its growth with two chambers less than G. elongatus (Fig 10). The chambers of G. elongatus gradually flatten between chambers 14–18, resulting in a decreasing S, whilst in G. tenellus they become rounder between chambers 14–16, which results in the divergent final shape that distinguishes between the sister species. As for G. rubescens, the final chamber of the scanned specimens of G. tenellus and G. elongatus is smaller than the penultimate chamber, which is indicative of the terminal reproductive stage.
Largest shells are typically found in G. conglobatus and G. ruber, but shell size is clearly not associated with the growth of more chambers: G. ruber has only 15 chambers in our dataset. The ontogenetic trajectory of G. ruber differs from all other species in its whorl number, which increases more steeply from chamber 9 onwards (Fig 10) in line with the higher angular increment between successive chambers. However, its expansion rate W is close to all other species except for G. conglobatus (Fig 10). G. ruber and G. conglobatus show a higher rate of size increase in consecutive chambers that the other species, such that for G. ruber the three last chambers occupy 94% of the total chamber volume (Fig 9B). The rates D and T are mirrored in their unevenness due to the abrupt decrease of the radius during the ontogeny (see S3 Table), the elevation of the trochospire and the tighter coiling axis. Finally, G. conglobatus has the largest test but its most distinctive feature is the increase of the sphericity between chambers 1 to 10 that is followed by compression between chambers 11 to 18. Its whorl expansion rate (W) is the highest throughout its ontogeny, but the formation of its high trochospire occurs over the last two chambers with an increase of T and a decrease of D.
Our Raupian analysis of the 3D ontogenetic trajectory of the five species could be used to determine changes in the position in the growth sequence when the juvenile, near-planispiral, many-chambered stage ends (onset of neanic stage sensu Brummer et al. [53]) and when the diagnostic, reproductive morphology is established (onset of adult stage sensu Brummer et al. [53]). The distinction of the ontogenetic stages in the CT reconstructions is based mainly on the parameters of chamber addition, but in several cases, the observed transitions could also be correlated with the emergence of further indicative traits, such as supplementary apertures. The analysis of the ontogenetic trajectories reveals that the allocation of chamber number and chamber volume to the ontogenetic trajectory remained similar between G. rubescens and G. ruber (Fig 9), but the other species show distinct differences in allocation. G. conglobatus differs most from the other species, exhibiting distinct juvenile-neanic stage with radially elongated chambers. G. elongatus shows a morphologically normal juvenile stage with 10 chambers and becomes trochospiral late in its ontogeny. Both species develop compressed chambers but the compression starts during the neanic stage at chamber 11 for G. conglobatus and at the onset of adult stage at chambers 14–15 for G. elongatus. By comparison, G. tenellus is much smaller, does not develop chamber compression and has fewer chambers (16).
Discussion
Strict dataset curation of the genetic dataset associated with the application of our nomenclature system confirms recent metabarcoding results which indicate that the biological diversity in planktonic foraminifera is limited [60,61]. We identified only eight genotypes and 14 basegroups within the five sequenced morphospecies of Globigerinoides, which likely covers the entire genotypic diversity in the genus. At the basegroup level, Globigerinoides conglobatus and Globigerinoides tenellus remain undersampled, but for the Globigerinoides ruber plexus and Globigerinoides elongatus, the sampling effort is sufficient to analyze the distribution of genetic diversity at all hierarchical levels (Fig 5, Table 2).
Our data confirm earlier work [26] in their conclusions that G. ruber ruber occurs only in the Atlantic, is the only type with test color and constitutes a single basegroup. We observe no other basegroup or genotype restricted to the Atlantic within the genus (Fig 4), but instead note the apparent absence of the basegroups G. elongatus Ia2 as well as G. ruber albus n.subsp. Ic1 and potentially Ib2 from the North Atlantic. Despite the fact that our first order Jackknifing (Table 1) and rarefaction analyses (Fig 4) suggest that the diversity in the North Atlantic may not have been captured entirely for G. ruber albus n.subsp. at the basegroup level, it does seem to be the case for G. elongatus at the basegroup level and for G. ruber albus n.subsp. at the genotype level. Therefore, the observed distribution pattern likely highlights an isolation of the tropical Atlantic from the Indian and Pacific Oceans.
Because of the equatorial position of the continents, the subtropical-tropical waters of the world oceans are only connected to a limited degree. At present, transport of tropical/subtropical marine plankton is largely unidirectional, from the Pacific to the Indian Ocean via the Indonesian throughflow, and from the Indian Ocean into the Atlantic via the Agulhas leakage. During glacial times, these connections likely became even more restricted [62]. Indeed, the disappearance of G. ruber ruber from the Indian and Pacific Oceans 120 kyrs ago [6] and its persistence in the Atlantic indicate a reduced ability to re-invade the Indian Ocean from the Atlantic. Dispersal from the Indian and Pacific Oceans into the Atlantic via Agulhas leakage is evidenced by the existence of a number of cosmopolitan basetypes (G. elongatus Ia1/Ia3 and G. ruber albus n.subsp. Ia1/Ia2/Ib2). In this scenario, the absence of G. elongatus Ia2 and G. ruber albus n.subsp. Ib2/Ic1 in the North Atlantic cannot be the result of dispersal limitation. Instead, the apparent accumulation of recently diverged endemic basegroups in the Pacific rather than the Atlantic (Figs 4 and 5) is reminiscent of the pattern observed in the hyperdiverse Globigerinella [63], where it has been ascribed to incumbency (expansion of a species into a new environment being prevented by an incumbent species with similar ecological preferences [64]). In our case, it might be that the Atlantic residents G. ruber ruber Ia1, G. ruber albus n.subsp. Ia1/Ia2/Ib1 and G. elongatus Ia1/Ia3 impede the establishment of invading genotypes recently diverged in the Indian and Pacific Oceans. The lack of diversity in the Atlantic endemic G. ruber ruber, compared to the cosmopolitan sister clade (Figs 2 and 4) suggests that no diversification occurs in the North Atlantic. Therefore, the Indian and Pacific Oceans seem to act as the primary source for biodiversity and the North Atlantic as a sink within the Globigerinoides genus.
Notwithstanding the pattern of limited connectivity between the Atlantic and the Indian and Pacific Oceans, the majority of the MOTUs has a cosmopolitan distribution within the (sub)tropical habitat of Globigerinoides, with co-occurrences at all taxonomic levels at the same stations (Figs 1 and 4), consistent with their apparently similar ecological niches (Fig 6). Although we did not sample G. ruber ruber in the South Atlantic, the distribution of the better-covered taxa is associated with higher SST in G. ruber albus n.subsp. compared to G. elongatus (Fig 6, Table 2). We acknowledge that our sampling of G. ruber ruber, with more sampling stations in the Caribbean and Mediterranean Seas compared to the central Atlantic, may have produced a biased view on the ecological preferences of this morphospecies. However, we are confident that our dataset of G. ruber albus and G. elongatus does not suffer from this limitation (Fig 4). The difference in thermal niches between G. ruber albus n.subsp. and G. elongatus has been a matter of debate since the seminal work of Wang [8]. Several studies replicated the observation of the preference of G. elongatus for colder waters compared to G. ruber albus n.subsp. akin to our observations [9,14–16,18,19,65,66], but observations of the absence of such differences have also been made. Indeed, a global synthesis of seasonally and depth-resolved sediment trap and plankton net observations [11] showed no statistically significant difference between G. ruber albus n.subsp. and G. elongatus in Mg/Ca composition of the shell. Studies conducted in the Gulf of Mexico [10,21] and in the central North Atlantic [67] showed similar absence of oxygen isotopic offsets between the morphospecies and argued that the difference in habitat, seasonal and calcifying depth is not systematic. Downcore analyses of Mg/Ca ratios from the southwest Pacific [15,20] showed that the difference between the two morphospecies was not stable though time and varied between 0 and 2°C in temperature space. This is consistent with the findings of Numberger et al. [18] in Mediterranean sediments, who noted oxygen isotopic offsets between the species, but the value and direction of the offset changed during the last 400 kyrs. Altogether, the niches of the two morphospecies may differ, but temperature sensitivity alone is unlikely to be the sole factor explaining the niche difference.
The conflicting observations on the degree of overlap between the ecological niches of G. ruber albus n.subsp. and G. elongatus raise the question of whether the degree of the overlap could be driven by ongoing diversification at the genotype and basegroup levels. In our analysis, we observe little to no ecological differences between the genotypes and basetypes of G. ruber albus n.subsp. and G. elongatus, except for (small) differences in temperature, salinity and productivity niches between G. elongatus basegroups Ia1 and Ia2 (Fig 6 and Table 2). Therefore, the regionally and temporally varying overlap between the ecological niches of the two morphospecies is unlikely to be the result of ecological differentiation among the constituent MOTUs. There is no evidence for the existence of ecological or biogeographic differentiation between the genotypes of G. ruber albus n.subsp. nor G. elongatus such as those that were discovered in morphospecies like Orbulina universa [68–70], Globorotalia inflata [71,72], Globorotalia truncatulinoides [73–75], Globigerina bulloides [76–79], Neogloboquadrina pachyderma [80–83] and Pulleniatina obliquiloculata [84,85]. An explanation invoking a vertical niche separation as observed in Hastigerina pelagica [86] is unlikely, because G. ruber albus n.subsp. and G. elongatus are both symbiont-bearing taxa limited to the photic zone and a consistent separation with depth or season would result in a constant isotopic offset, which contrasts general observations (see above).
Although abiotic factors, such as temperature, are important drivers of plankton community structure [87,88], recent studies have shown that biotic interactions may be even more important drivers of plankton diversification. Analyses of plankton metacommunity structure showed that abiotic factors alone explained only 18% of the variability in the distribution of environmental OTUs [89], leaving biotic interactions as the main driver of ecological and biological diversification in the open ocean. Photosymbiosis is the biotic interaction that has been most studied in foraminifera [90] and is of interest to paleoceanographers, not only because it ties photosymbiotic species to photic depths, but also because it impacts the incorporation of stable carbon isotopes and trace elements in the calcareous shell [91–93]. Photophysiology [92,94–98] investigations have documented the dynamic relationship between the foraminifera and their photosymbionts, but the diversity of these interactions, including other interactions such as parasitism or commensalism, has not yet been systematically resolved. Indeed, Shaked and de Vargas [99] found 21 phylotypes of the dinoflagellate Symbiodinium hosted by four morphospecies of tropical planktonic foraminifera, including G. ruber and G. conglobatus, and suggested that this number most likely represents the lower bound of the true symbiotic diversity, leaving ample space for differentiation due to preference for different symbiont strains.
Planktonic foraminifera, like many protists living in the oligotrophic ocean, are capable of mixotrophy (capable of autotrophy by symbiosis and heterotrophy) and the type of mixotrophy influences the biogeography and seasonality of the mixotrophs hosting the symbionts [100]. We hypothesize that the position in the trophic network occupied by planktonic foraminifera may control when and where they calcify their shell. The control of temperature on planktonic foraminifera individual species abundance and occurrence could be indirect and the physico-chemical condition of the water column that the planktonic foraminifera record may reflect their relationships with other organisms rather than a mere thermal response. In this scenario, temperature alone would not explain evolution in planktonic foraminifera [101] and vital effects impacting the incorporation of carbon isotopes could have varied through time as a function of varying symbiotic association and mixotrophy level [93]. Indeed, a prominent role of biotic factors in the diversification of Globigerinoides species is consistent with the lack of physical niche differentiation at the level of genotypes and basegroups. The large number of apparently recently diverging basegroups could result from a high turnover driven by biotic interactions which rarely leads to persistent separation of lineages, resulting in a continuous diversification in the genus throughout the late Neogene and Quaternary (Fig 7), without a clear partitioning of the ecological space along abiotic factors.
Diversification at the cryptic level in the genus likely reflects biotic interactions, but it remains to be explained why and how the morphological evolution and genetic divergence are disconnected at the morphospecies level. For instance, G. ruber ruber and G. ruber albus n.subsp. diverged around ~6.7 Ma and remained morphologically identical, whereas G. elongatus and G. conglobatus diverged around 8.3 Ma (Fig 7) but are morphologically distinct from juvenile to adult. Similarly, G. tenellus and G. elongatus, which are morphologically dissimilar diverged around 2.4 Ma and this event could be concomitant with the divergence time of the constitutive genotype of G. conglobatus and G. ruber albus n.subsp (Fig 7). Because of a similarity in shape, G. tenellus was previously considered a sister species of G. rubescens. The apparent similarity motivated us to analyze the ontogeny of this species as well. Our strategy was to recover the potential phylogenetic information contained in the ontogenetic development of the five extant morphospecies of Globigerinoides and to use Globoturborotalita rubescens as an outgroup. Because of the time-consuming nature of 3D analysis, we limited our approach to a single representative specimen per species to obtain the main differences in the ontogenetic development between species. We acknowledge that intra-species variability in the ontogenetic development exists [53] and that our study design prevents assessing the magnitude of this variability. Nevertheless, the observed contrasting patterns of growth allocation to ontogenetic stages are substantial and associated with systematic changes in chamber shape and growth pattern (Fig 9), in a manner that can be best described in the light of heterochrony [102]. Heterochrony is defined as evolutionary change in the rate and timing of ontogenetic development. Although heterochrony is a concept developed to understand the connection between evolution and development in multicellular organisms, we apply it in a broad sense to planktonic foraminifera because the sequential growth of their tests preserves the sequence of shapes during individual growth. Also, we stress that heterochrony as a concept does not explain the mechanistic cause for evolutionary change, but provides a framework in which the emergence of the divergent adult shapes can be described through changes in the ontogenetic trajectory [102].
In this heterochronic framework, we observe that the G. rubescens specimen displays the most stable development with relatively little change in the shape of its chambers during ontogeny compared to the other species (Figs 9 and 10). Considering this morphospecies as outgroup (given its phylogenetic position; Figs 2 and 7), we explore the divergence of adult morphologies of the individual species in terms of Raupian alterations in the ontogenetic trajectory and the successive emergence of new characters. Compared to G. rubescens, the morphological innovations in Globigerinoides are the emergence of elongate chambers, compressed chambers and supplementary apertures. Chamber elongation is restricted to the juvenile stage of the G. conglobatus specimen and it is followed by compression in the neanic-adult stages of large G. conglobatus. Chamber compression also occurs in the adult stage of G. elongatus and its absence in small G. tenellus hints at heterochrony by dwarfing. Supplementary apertures are lacking in the small ancestral G. rubescens but are typically found in in the sister clade and their reduction to the last 1–2 chambers in G. tenellus is consistent with heterochrony by dwarfing. In G. tenellus a single secondary aperture is typically present in the final chamber, whereas all other species of the genus develop at least in the final chambers two supplementary apertures per chamber.
The ratio S describing the evolution of the roundness of the chambers is more stable during the ontogeny of G. ruber compared to the four other species (Fig 10). The analyzed specimen is large (400 μm) for the few (15) chambers it has, and lacks chamber compression in comparison to G. elongatus and G. conglobatus, indicating that G. ruber may have a neotenic ontogenetic trajectory. Neoteny is characterized by a conservation of juvenile features during the adult stage, reduced compression of the last chamber in the case of G. ruber, without a change of size. It is associated with a steeper increase of chamber size at a higher angular increment towards the end of the growth. This scenario would be consistent with the hypothesis that G. ruber evolved from G. obliquus (which has more compressed chambers) as proposed by Aurahs [27]. In contrast, the ontogenetic trajectory of G. conglobatus appears hypermorphic, which is characterized by larger final size. Finally, G. elongatus and G. tenellus seem to follow similar ontogenetic paths and to differ in the last three chambers, with the compression of the chambers of G. elongatus and the increase of the roundness of G. tenellus chambers. Also, G. tenellus has larger chambers through its ontogeny and its final size is smaller than G. elongatus, suggesting progenesis. Progenesis is defined as a loss of an adult feature, the final compressed chamber akin to what we hypothesize for G. ruber, but in this case associated with a reduction in size due to a premature interruption of the growth. In terms of size, G. tenellus is one of the few known examples of dwarfing in planktonic foraminifera, but unlike the fossil species Globorotalia exilis, Globorotalia miocenica and Morozovelloides crassatus the dwarfing in G. tenellus does not (yet) seem to be associated with a reduction of abundance preceding extinction [103].
Evolution through heterochrony could provide an explanation for the erroneous taxonomic placement of G. elongatus as a sister to G. ruber that led to the informal delimitation G. ruber s.l. and s.s. by Wang [8]. Indeed, we hypothesize that G. elongatus may not attain the size and shape of G. conglobatus because it has smaller chambers, which are less compressed, and could consequently converge towards the size and shape of G. ruber. Similarly, G. tenellus may create a morphological convergence with G. rubescens despite having markedly different pre-adult ontogenetic trajectories (Fig 9). The presence of supplementary apertures in G. tenellus is thus an apomorphy of Globigerinoides. Based on our observations, we proposed several interpretations of the molecular phylogeny topology that would be in agreement with the morphology, taking into account the heterochronic development within Globigerinoides genus (Fig 11).
Similar to previous studies [29,30] our results show that CT-scanning offers a promising avenue for ontogenetic analysis and resolve phylogenetic relationships among extinct species of planktonic foraminifera [104]. We recognize that we cannot draw firm conclusions from our analysis because of the limited amount of specimen analyzed, and stress the need for replicate analysis to confirm our results. Even though ontogenetic analysis may not explain what triggered the divergence and convergence of juvenile and adult morphologies, it could provide a viable explanation for the apparent disconnection between morphological and genetic divergence. Heterochrony is a process through which large changes in adult morphology could be achieved at genetically low cost [102], creating an impression of large change not matched by the degree of genetic kinship.
Supporting information
S1 Fig [a]
Light microscopy images of the specimen CA1261 identified as and from which sequence match the type IIb of Aurahs et al [,].
S2 Fig [tif]
Light microscopy images of the holotype of (C319) and paratypes (C208, C281, C329) of . n.subsp.
S1 Table [xlsx]
Metadata and taxonomy of the Sanger sequences used in the study.
S2 Table [xlsx]
Taxonomic equivalence between the existing taxonomic nomenclatures proposed in the literature and our updated molecular taxonomy.
S3 Table [xlsx]
Volume, Cartesian coordinates and parameters of the Raup’s model measured on individual chambers of the five selected morphological species (Figs and ).
Zdroje
1. Cushman JA. An outline of a re-classification of the foraminifera. Contrib from Cushman Found Foraminifer Res. 1927;3: 1–105.
2. Banner FT, Blow FW. Some primary types of species belonging to the superfamily Globigerinaceae. Contrib from Cushman Found Foraminifer Res. 1960;11: 1–41.
3. Parker FL. Planktonic Foraminiferal Species in Pacific Sediments. Micropaleontology. 1962;8: 219. doi: 10.2307/1484745
4. Kennett JP, Srinivasan MS. Neogene Planktonic Foraminifera: A Phylogenetic Atlas. Hutchinson. Troudsburg; 1983.
5. Deuser WG, Ross EH. Seasonally abundant planktonic foraminifera of the Sargasso Sea; succession, deep-water fluxes, isotopic compositions, and paleoceanographic implications. J Foraminifer Res. 1989;19: 268–293. doi: 10.2113/gsjfr.19.4.268
6. Thompson PR, Bé H, Duplessy J-C, Shackleton NJ. Disappearance of pink-pigmented Globigerinoides ruber at 120,000 yr BP in the Indian and Pacific Oceans. Nature. 1979;280: 554–558.
7. Robbins LL, Healy-Williams N. Toward a classification of planktonic foraminifera based on biochemical, geochemical, and morphological criteria. J Foraminifer Res. 2009;21: 159–167. doi: 10.2113/gsjfr.21.2.159
8. Wang L. Isotopic signals in two morphotypes of Globigerinoides ruber (white) from the South China Sea: implications for monsoon climate change during the last glacial cycle. Palaeogeogr Palaeoclimatol Palaeoecol. 2000;161: 381–394. doi: 10.1016/S0031-0182(00)00094-8
9. Steinke S, Chiu H-Y, Yu P-S, Shen C-C, Löwemark L, Mii H-S, et al. Mg/Ca ratios of two Globigerinoides ruber (white) morphotypes: Implications for reconstructing past tropical/subtropical surface water conditions. Geochemistry, Geophys Geosystems. 2005;6: 1–12. doi: 10.1029/2005GC000926
10. Richey JN, Thirumalai K, Khider D, Reynolds CE, Partin JW, Quinn TM. Considerations for Globigerinoides ruber (White and Pink) Paleoceanography: Comprehensive Insights From a Long‐Running Sediment Trap. Paleoceanogr Paleoclimatology. 2019;34: 353–373. doi: 10.1029/2018PA003417
11. Gray WR, Weldeab S, Lea DW, Rosenthal Y, Gruber N, Donner B, et al. The effects of temperature, salinity, and the carbonate system on Mg/Ca in Globigerinoides ruber (white): A global sediment trap calibration. Earth Planet Sci Lett. 2018;482: 607–620. doi: 10.1016/j.epsl.2017.11.026
12. Jentzen A, Schönfeld J, Schiebel R. Assessment of the Effect of Increasing Temperature On the Ecology and Assemblage Structure of Modern Planktic Foraminifers in the Caribbean and Surrounding Seas. J Foraminifer Res. 2018;48: 251–272. doi: 10.2113/gsjfr.48.3.251
13. Mojtahid M, Manceau R, Schiebel R, Hennekam R, de Lange GJ. Thirteen thousand years of southeastern Mediterranean climate variability inferred from an integrative planktic foraminiferal-based approach. Paleoceanography. 2015;30: 402–422. doi: 10.1002/2014PA002705
14. Antonarakou A, Kontakiotis G, Mortyn PG, Drinia H, Sprovieri M, Besiou E, et al. Biotic and geochemical (δ18O, δ13C, Mg/Ca, Ba/Ca) responses of Globigerinoides ruber morphotypes to upper water column variations during the last deglaciation, Gulf of Mexico. Geochim Cosmochim Acta. 2015;170: 69–93. doi: 10.1016/j.gca.2015.08.003
15. Regoli F, de Garidel-Thoron T, Tachikawa K, Jian Z, Ye L, Droxler AW, et al. Progressive shoaling of the equatorial Pacific thermocline over the last eight glacial periods. Paleoceanography. 2015;30: 439–455. doi: 10.1002/2014PA002696
16. Kawahata H. Stable isotopic composition of two morphotypes of Globigerinoides ruber (white) in the subtropical gyre in the North Pacific. Paleontol Res. 2005;9: 27–35.
17. Löwemark L, Hong W-L, Yui T-F, Hung G-W. A test of different factors influencing the isotopic signal of planktonic foraminifera in surface sediments from the northern South China Sea. Mar Micropaleontol. 2005;55: 49–62. doi: 10.1016/j.marmicro.2005.02.004
18. Numberger L, Hemleben C, Hoffmann R, Mackensen A, Schulz H, Wunderlich J-M, et al. Habitats, abundance patterns and isotopic signals of morphotypes of the planktonic foraminifer Globigerinoides ruber (d’Orbigny) in the eastern Mediterranean Sea since the Marine Isotopic Stage 12. Mar Micropaleontol. 2009;73: 90–104. doi: 10.1016/j.marmicro.2009.07.004
19. Carter A, Clemens S, Kubota Y, Holbourn A, Martin A. Differing oxygen isotopic signals of two Globigerinoides ruber (white) morphotypes in the East China Sea: Implications for paleoenvironmental reconstructions. Mar Micropaleontol. 2017;131: 1–9. doi: 10.1016/j.marmicro.2017.01.001
20. Schmitt A, Elliot M, Thirumalai K, La C, Bassinot F, Petersen J, et al. Single foraminifera Mg/Ca analyses of past glacial-interglacial temperatures derived from G. ruber sensu stricto and sensu lato morphotypes. Chem Geol. 2018; doi: 10.1016/j.chemgeo.2018.11.007
21. Thirumalai K, Richey JN, Quinn TM, Poore RZ. Globigerinoides ruber morphotypes in the Gulf of Mexico: A test of null hypothesis. Sci Rep. 2014;4: 1–7. doi: 10.1038/srep06018 25109442
22. Darling KF, Wade CM, Kroon D, Brown AJL. Planktic foraminiferal molecular evolution and their polyphyletic origins from benthic taxa. Mar Micropaleontol. 1997;30: 251–266.
23. Darling KF, Wade CM, Kroon D, Brown AJL, Bijma J. The diversity and distribution of modern planktic foraminiferal small subunit ribosomal RNA genotypes and their potential as tracers of present and past circulation. Paleoceanography. 1999;14: 3–12. doi: 10.1029/1998PA900002
24. Darling KF, Wade CM. The genetic diversity of planktic foraminifera and the global distribution of ribosomal RNA genotypes. Mar Micropaleontol. 2008;67: 216–238. doi: 10.1016/j.marmicro.2008.01.009
25. Kuroyanagi A, Tsuchiya M, Kawahata H, Kitazato H. The occurrence of two genotypes of the planktonic foraminifer Globigerinoides ruber (white) and paleo-environmental implications. Mar Micropaleontol. 2008;68: 236–243. doi: 10.1016/j.marmicro.2008.04.004
26. Aurahs R, Grimm GW, Hemleben V, Hemleben C, Kucera M. Geographical distribution of cryptic genetic types in the planktonic foraminifer Globigerinoides ruber. Mol Ecol. 2009;18: 1692–1706. doi: 10.1111/j.1365-294X.2009.04136.x 19302352
27. Aurahs R, Treis Y, Darling K, Kucera M. A revised taxonomic and phylogenetic concept for the planktonic foraminifer species Globigerinoides ruber based on molecular and morphometric evidence. Mar Micropaleontol. 2011;79: 1–14. doi: 10.1016/j.marmicro.2010.12.001
28. Morard R, Escarguel G, Weiner AKM, André A, Douady CJ, Wade CM, et al. Nomenclature for the Nameless: A Proposal for an Integrative Molecular Taxonomy of Cryptic Diversity Exemplified by Planktonic Foraminifera. Syst Biol. 2016;65: 925–940. doi: 10.1093/sysbio/syw031 27073250
29. Caromel AGM, Schmidt DN, Rayfield EJ. Ontogenetic constraints on foraminiferal test construction. Evol Dev. 2017;19: 157–168. doi: 10.1111/ede.12224 28463472
30. Caromel AGM, Schmidt DN, Fletcher I, Rayfield EJ. Morphological Change During The Ontogeny Of The Planktic Foraminifera. J Micropalaeontology. 2015; 2014–017.
31. Weiner AKM, Morard R, Weinkauf MF, Darling KF, André A, Quillévéré F, et al. Methodology for single-cell genetic analysis of planktonic foraminifera for studies of protist diversity and evolution. Front Mar Sci. 2016;3: 1–15. doi: 10.3389/fmars.2016.00255
32. Wuyts J, Van de Peer Y, De Wachter R. Distribution of substitution rates and location of insertion sites in the tertiary structure of ribosomal RNA. Nucleic Acids Res. 2001;29: 5017–5028. doi: 10.1093/nar/29.24.5017 11812832
33. Morard R, Darling KF, Mahé F, Audic S, Ujiié Y, Weiner AKM, et al. PFR2: a curated database of planktonic foraminifera 18S ribosomal DNA as a resource for studies of plankton ecology, biogeography and evolution. Mol Ecol Resour. 2015;15: 1472–1485. doi: 10.1111/1755-0998.12410 25828689
34. Pawlowski J, Holzmann M. A plea for DNA barcoding of foraminifera. J Foraminifer Res. 2014;44: 62–67. doi: 10.2113/gsjfr.44.1.62
35. Puillandre N, Lambert A, Brouillet S, Achaz G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol Ecol. 2012;21: 1864–1877. doi: 10.1111/j.1365-294X.2011.05239.x 21883587
36. Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29: 2869–2876. doi: 10.1093/bioinformatics/btt499 23990417
37. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol Biol Evol. 2013;30: 772–780. doi: 10.1093/molbev/mst010 23329690
38. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst Biol. 2010;59: 307–321. doi: 10.1093/sysbio/syq010 20525638
39. Lefort V, Longueville J, Gascuel O. SMS: Smart Model Selection in PhyML. Mol Biol Evol. 2017;34: 2422–2424. doi: 10.1093/molbev/msx149 28472384
40. Ujiié Y, Lipps JH. Cryptic diversity in planktonic foraminifera in the northwest Pacific ocean. J Foraminifer Res. 2009;39: 145–154.
41. André A, Quillévéré F, Morard R, Ujiié Y, Escarguel G, De Vargas C, et al. SSU rDNA divergence in planktonic foraminifera: Molecular taxonomy and biogeographic implications. PLoS One. 2014;9. doi: 10.1371/journal.pone.0104641 25119900
42. NASA Goddard Space Flight Center, Ocean Ecology Laboratory, Ocean Biology Processing Group. Moderate-resolution Imaging Spectroradiometer (MODIS) Aqua Sea Surface Temperature Data; 2014 Reprocessing. NASA OB.DAAC, Greenbelt, MD, USA.
43. NASA Goddard Space Flight Center, Ocean Ecology Laboratory, Ocean Biology Processing Group. Moderate-resolution Imaging Spectroradiometer (MODIS) Aqua Chlorophyll Data; 2018 Reprocessing. NASA OB.DAAC, Greenbelt, MD, USA.
44. NASA Goddard Space Flight Center, Ocean Ecology Laboratory, Ocean Biology Processing Group. Moderate-resolution Imaging Spectroradiometer (MODIS) Aqua Particulate Organic Carbon Data; 2018 Reprocessing. NASA OB.DAAC, Greenbelt, MD, USA.
45. Schmidtko S, Johnson GC, Lyman JM. MIMOC: A global monthly isopycnal upper-ocean climatology with mixed layers. J Geophys Res Ocean. 2013;118: 1658–1672. doi: 10.1002/jgrc.20122
46. Behrenfeld MJ, Falkowski PG. A consumer’s guide to phytoplankton primary productivity models. Limnol Oceanogr. 1997;42: 1479–1491. doi: 10.4319/lo.1997.42.7.1479
47. Hammer Ø, Harper D a. T, Ryan PD. Paleontological statistics software package for education and data analysis. Palaeontol Electron. 2001;4: 9–18.
48. Aze T, Ezard THG, Purvis A, Coxall HK, Stewart DRM, Wade BS, et al. A phylogeny of Cenozoic macroperforate planktonic foraminifera from fossil data. Biol Rev. 2011;86: 900–927. doi: 10.1111/j.1469-185X.2011.00178.x 21492379
49. Kucera M, Schönfeld J. The origin of modern oceanic foraminiferal faunas and Neogene climate change. In: The Micropalaeontological Society SP, editor. Deep-Time Perspectives on Climate Change: Marrying the Signal from Computer Models and Biological Proxies. London: The Geological Society; 2007. pp. 409–425.
50. Drummond AJ, Suchard M a. Bayesian random local clocks, or one rate to rule them all. BMC Biol. 2010;8: 114. doi: 10.1186/1741-7007-8-114 20807414
51. Rambaut A. Tree Figure Drawing Tool Version 1.3.1. University of Edinburgh; 2009. http://tree.bio.ed.ac.uk/software/figtree/
52. Fisher G, Arz H, Baschek B, Bassek D, Costa E, Dierssen H, et al. Report and preliminary results of METEOR-Cruise M34/4, Recife-Bridgetown, 19.3–15.4.1996. Berichte, Fachbereich Geowissenschaften, Universität Bremen, No 80, 105 pp. 1996;
53. Brummer G-J a., Hemleben C, Spindler M. Ontogeny of extant spinose planktonic foraminifera (Globigerinidae): A concept exemplified by Globigerinoides sacculifer (Brady) andG. Ruber (d’Orbigny). Mar Micropaleontol. 1987;12: 357–381. doi: 10.1016/0377-8398(87)90028-4
54. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, et al. User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. Neuroimage. 2006;31: 1116–1128. doi: 10.1016/j.neuroimage.2006.01.015 16545965
55. Raup DM. Geometric Analysis of Shell Coiling: Coiling in Ammonoids. J Paleontol. 1967;41: 43–65.
56. Darling KF, Kroon D, Wade CM, Leigh J. Molecular Phylogeny of the planktic foraminifera. J Foraminifer Res. 1996;26: 324–330.
57. Pawlowski J, Bolivar I, Fahrni JF, de Vargas C, Gouy M, Zaninetti L. Extreme differences in rates of molecular evolution of foraminifera revealed by comparison of ribosomal DNA sequences and the fossil record. Mol Biol Evol. 1997;14: 498–505. doi: 10.1093/oxfordjournals.molbev.a025786 9159927
58. de Vargas C, Zaninetti L, Hilbrecht H, Pawlowski J. Phylogeny and rates of molecular evolution of planktonic foraminifera: SSU rDNA sequences compared to the fossil record. J Mol Evol. 1997;45: 285–294. doi: 10.1007/pl00006232 9302323
59. André A, Quillévéré F, Morard R, Ujiié Y, Escarguel G, De Vargas C, et al. SSU rDNA divergence in planktonic foraminifera: molecular taxonomy and biogeographic implications. PLoS One. 2014;9: e104641. doi: 10.1371/journal.pone.0104641 25119900
60. Morard R, Mahé F, Romac S, Poulain J, Kucera M, De Vargas C. Surface ocean metabarcoding confirms limited diversity in planktonic foraminifera but reveals unknown hyper-abundant lineages. Sci Rep. 2018;8:2539: 1–10. doi: 10.1038/s41598-018-20833-z 29416071
61. Morard R, Vollmar NM, Greco M, Kucera M. Unassigned diversity of planktonic foraminifera from environmental sequencing revealed as known but neglected species. PLoS One. 2019;14: e0213936. doi: 10.1371/journal.pone.0213936 30897140
62. Peeters FJC, Acheson R, Brummer G-J a, De Ruijter WPM, Schneider RR, Ganssen GM, et al. Vigorous exchange between the Indian and Atlantic oceans at the end of the past five glacial periods. Nature. 2004;430: 661–5. doi: 10.1038/nature02785 15295596
63. Weiner AKM, Weinkauf MFG, Kurasawa A, Darling KF, Kucera M, Grimm GW. Phylogeography of the tropical planktonic foraminifera lineage Globigerinella reveals isolation inconsistent with passive dispersal by ocean currents. PLoS One. 2014;9: e92148. doi: 10.1371/journal.pone.0092148 24663038
64. Algar AC, Mahler DL, Glor RE, Losos JB. Niche incumbency, dispersal limitation and climate shape geographical distributions in a species-rich island adaptive radiation. Glob Ecol Biogeogr. 2013;22: 391–402. doi: 10.1111/geb.12003
65. Kuroyanagi A, Kawahata H. Vertical distribution of living planktonic foraminifera in the seas around Japan. Mar Micropaleontol. 2004;53: 173–196. doi: 10.1016/j.marmicro.2004.06.001
66. Capotondi L, Girone A, Lirer F, Bergami C, Verducci M, Vallefuoco M, et al. Central Mediterranean Mid-Pleistocene paleoclimatic variability and its association with global climate. Palaeogeogr Palaeoclimatol Palaeoecol. 2016;442: 72–83. doi: 10.1016/j.palaeo.2015.11.009
67. Bonfardeci A., Caruso A., Bartolini A., Bassinot F., Blanc-Valleron M.-M., 2018. Distribution and ecology of the Globigerinoides ruber—Globigerinoides elongatus morphotypes in the Azores region during the late Pleistocene-Holocene. Palaeogeogr. Palaeoclimatol. Palaeoecol. 491, 92–111. doi: 10.1016/j.palaeo.2017.11.052
68. de Vargas C, Norris R, Zaninetti L, Gibb SW, Pawlowski J. Molecular evidence of cryptic speciation in planktonic foraminifers and their relation to oceanic provinces. Proc Natl Acad Sci U S A. 1999;96: 2864–2868. doi: 10.1073/pnas.96.6.2864 10077602
69. Morard R, Quillévéré F, Escarguel G, Ujiie Y, de Garidel-Thoron T, Norris RD, et al. Morphological recognition of cryptic species in the planktonic foraminifer Orbulina universa. Mar Micropaleontol. 2009;71: 148–165. doi: 10.1016/j.marmicro.2009.03.001
70. Marshall BJ, Thunell RC, Spero HJ, Henehan MJ, Lorenzoni L, Astor Y. Morphometric and stable isotopic differentiation of Orbulina universa morphotypes from the Cariaco Basin, Venezuela. Mar Micropaleontol. 2015;120: 46–64. doi: 10.1016/j.marmicro.2015.08.001
71. Morard R, Quillévéré F, Douady CJ, de Vargas C, de Garidel-Thoron T, Escarguel G. Worldwide genotyping in the planktonic foraminifer Globoconella inflata: Implications for life history and paleoceanography. PLoS One. 2011;6: e26665. doi: 10.1371/journal.pone.0026665 22028935
72. Morard R, Reinelt M, Chiessi CM, Groeneveld J, Kucera M. Tracing shifts of oceanic fronts using the cryptic diversity of the planktonic foraminifera Globorotalia inflata. Paleoceanography. 2016;31.
73. de Vargas C, Renaud S, Hilbrecht H, Pawlowski J. Pleistocene adaptive radiation in Globorotalia truncatulinoides: genetic, morphologic, and environmental evidence. Paleobiology. 2001;27: 104–125.
74. Quillévéré F, Morard R, Escarguel G, Douady CJ, Ujiié Y, de Garidel-Thoron T, et al. Global scale same-specimen morpho-genetic analysis of Truncorotalia truncatulinoides: A perspective on the morphological species concept in planktonic foraminifera. Palaeogeogr Palaeoclimatol Palaeoecol. 2013;391. doi: 10.1016/j.palaeo.2011.03.013
75. Ujiié Y, Asami T. Temperature is not responsible for left-right reversal in pelagic unicellular zooplanktons. J Zool. 2014;293: 16–24. doi: 10.1111/jzo.12095
76. Darling KF, Wade CM, Stewart I a, Kroon D, Dingle R, Brown a J. Molecular evidence for genetic mixing of Arctic and Antarctic subpolar populations of planktonic foraminifers. Nature. 2000;405: 43–7. doi: 10.1038/35011002 10811211
77. Kucera M, Darling KF. Cryptic species of planktonic foraminifera: their effect on palaeoceanographic reconstructions. Philos Trans A Math Phys Eng Sci. 2002;360: 695–718. doi: 10.1098/rsta.2001.0962 12804300
78. Morard R, Quillévéré F, Escarguel G, de Garidel-Thoron T, de Vargas C, Kucera M. Ecological modeling of the temperature dependence of cryptic species of planktonic Foraminifera in the Southern Hemisphere. Palaeogeogr Palaeoclimatol Palaeoecol. 2013;391: 13–33. doi: 10.1016/j.palaeo.2013.05.011
79. Sadekov AY, Darling KF, Ishimura T, Wade CM, Kimoto K, Singh AD, et al. Geochemical imprints of genotypic variants of Globigerina bulloides in the Arabian Sea. Paleoceanography. 2016;
80. Darling KF, Kucera M, Pudsey CJ, Wade CM. Molecular evidence links cryptic diversification in polar planktonic protists to Quaternary climate dynamics. Proc Natl Acad Sci U S A. 2004;101: 7657–7662. doi: 10.1073/pnas.0402401101 15136732
81. Darling KF, Kucera M, Kroon D, Wade CM. A resolution for the coiling direction paradox in Neogloboquadrina pachyderma. Paleoceanography. 2006;21: PA2011. doi: 10.1029/2005PA001189
82. Darling KF, Kucera M, Wade CM. Global molecular phylogeography reveals persistent Arctic circumpolar isolation in a marine planktonic protist. Proc Natl Acad Sci U S A. 2007;104: 5002–5007. doi: 10.1073/pnas.0700520104 17360336
83. André A, Quillévéré F, Schiebel R, Morard R, Howa H, Meilland J, et al. Disconnection between genetic and morphological diversity in the planktonic foraminifer Neogloboquadrina pachyderma from the Indian sector of the Southern Ocean. Mar Micropaleontol. 2018;144: 14–24. doi: 10.1016/j.marmicro.2018.10.001
84. Ujiié Y, Asami T, de Garidel-Thoron T, Liu H, Ishitani Y, de Vargas C. Longitudinal differentiation among pelagic populations in a planktic foraminifer. Ecol Evol. 2012;2: 1725–1737. doi: 10.1002/ece3.286 22957176
85. Ujiié Y, Ishitani Y. Evolution of a Planktonic Foraminifer during Environmental Changes in the Tropical Oceans. PLoS One. 2016;11: 1–16. doi: 10.1371/journal.pone.0148847 26886349
86. Weiner A, Aurahs R, Kurasawa A, Kitazato H, Kucera M. Vertical niche partitioning between cryptic sibling species of a cosmopolitan marine planktonic protist. Mol Ecol. 2012;21: 4063–4073. doi: 10.1111/j.1365-294X.2012.05686.x 22738662
87. Sunagawa S, Coelho LP, Chaffron S, Kultima JR, Labadie K, Salazar G, et al. Structure and function of the global ocean microbiome. 2015;348: 1–10.
88. Rutherford S D’Hondt S, Prell W. Environmental controls on the geographic distribution of zooplankton diversity. Nature. 1999;400: 749–753. doi: 10.1038/23449
89. Lima-Mendez G, Faust K, Henry N, Decelle J, Colin S, Carcillo F, et al. Determinants of community structure in the global plankton interactome. Science. 2015;348: 1262073–1262073. doi: 10.1126/science.1262073 25999517
90. Decelle J, Colin S, Foster RA. Photosymbiosis in Marine Planktonic Protists. In: Ohtsuka S, Suzaki T, Horiguchi T, Suzuki N, Not F, editors. Marine Protists. Tokyo: Springer Japan; 2015. pp. 465–500.
91. Ezard THG, Edgar KM, Hull PM. Environmental and biological controls on size-specific δ 13 C and δ 18 O in recent planktonic foraminifera. Paleoceanography. 2015;30: 151–173. doi: 10.1002/2014PA002735
92. Takagi H, Moriya K, Ishimura T, Suzuki A, Kawahata H, Hirano H. Exploring photosymbiotic ecology of planktic foraminifers from chamber-by-chamber isotopic history of individual foraminifers. Paleobiology. 2015;41: 108–121. doi: 10.1017/pab.2014.7
93. Edgar KM, Hull PM, Ezard THG. Evolutionary history biases inferences of ecology and environment from δ13C but not δ18O values. Nat Commun. Springer US; 2017;8: 1–9. doi: 10.1038/s41467-017-01154-7 29062052
94. Takagi H., Kimoto K., Fujiki T., Moriya K., 2018. Effect of nutritional condition on photosymbiotic consortium of cultured Globigerinoides sacculifer (Rhizaria, Foraminifera). Symbiosis 76, 25–39. doi: 10.1007/s13199-017-0530-3 30147222
95. Takagi H, Kimoto K, Fujiki T, Kurasawa A, Moriya K, Hirano H. Ontogenetic dynamics of photosymbiosis in cultured planktic foraminifers revealed by fast repetition rate fluorometry. Mar Micropaleontol. 2016;122: 44–52. doi: 10.1016/j.marmicro.2015.10.003
96. Fujiki T, Takagi H, Kimoto K, Kurasawa a., Yuasa T, Mino Y. Assessment of algal photosynthesis in planktic foraminifers by fast repetition rate fluorometry. J Plankton Res. 2014;36: 1403–1407. doi: 10.1093/plankt/fbu083
97. LeKieffre C, Spero HJ, Russell AD, Fehrenbacher JS, Geslin E, Meibom A. Assimilation, translocation, and utilization of carbon between photosynthetic symbiotic dinoflagellates and their planktic foraminifera host. Mar Biol. Springer Berlin Heidelberg; 2018;165: 1–15. doi: 10.1007/s00227-018-3362-7
98. Takagi H, Kimoto K, Fujiki T, Saito H, Schmidt C, Kucera M, et al. Characterizing photosymbiosis in modern planktonic foraminifera. Biogeosciences Discuss. 2019; 1–32. doi: 10.5194/bg-2019-145
99. Shaked Y, de Vargas C. Pelagic photosymbiosis: rDNA assessment of diversity and evolution of dinoflagellate symbionts and planktonic foraminiferal hosts. Mar Ecol Prog Ser. 2006;325: 59–71.
100. Leles SG, Mitra A, Flynn KJ, Stoecker DK, Hansen PJ, Calbet A, et al. Oceanic protists with different forms of acquired phototrophy display contrasting biogeographies and abundance. Proc R Soc B Biol Sci. 2017;284. doi: 10.1098/rspb.2017.0664 28768886
101. Brombacher A, Wilson PA, Bailey I, Ezard THG. Temperature is a poor proxy for synergistic climate forcing of plankton evolution. Proc R Soc B Biol Sci. 2018;285: 1–8. doi: 10.1098/rspb.2018.0665 30051846
102. Klingenberg CP. Heterochrony and allometry: The analysis of evolutionary change in ontogeny. Biol Rev. 1998;73: 79–123. doi: 10.1017/s000632319800512x 9569772
103. Wade BS, Olsson RK. Investigation of pre-extinction dwarfing in Cenozoic planktonic foraminifera. Palaeogeogr Palaeoclimatol Palaeoecol. 2009;284: 39–46. doi: 10.1016/j.palaeo.2009.08.026
104. Fordham BG, Aze T, Haller C, Zehady AK, Pearson PN, Ogg JG, et al. Future-proofing the Cenozoic macroperforate planktonic foraminifera phylogeny of Aze & others (2011). Jonkers L, editor. PLoS One. 2018;13: e0204625. doi: 10.1371/journal.pone.0204625 30379910
105. Spring S. NOAA Atlas NESDIS 73 WORLD OCEAN ATLAS 2013 Volume 1: Temperature. 2013;1.
106. Siccha M, Kucera M. ForCenS, a curated database of planktonic foraminifera census counts in marine surface sediment samples. Sci Data. 2017;4:170109: 1–12. doi: 10.1038/sdata.2017.109 28829434
107. Schlitzer R. Ocean Data View. 2018. Available: http://odv.awi.de
108. R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria; 2014. http://www.r-project.org/
109. Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer-Verlag; 2009.
Článek vyšel v časopise
PLOS One
2019 Číslo 12
- S diagnostikou Parkinsonovy nemoci může nově pomoci AI nástroj pro hodnocení mrkacího reflexu
- Je libo čepici místo mozkového implantátu?
- Metamizol jako analgetikum první volby: kdy, pro koho, jak a proč?
- Pomůže v budoucnu s triáží na pohotovostech umělá inteligence?
- AI může chirurgům poskytnout cenná data i zpětnou vazbu v reálném čase
Nejčtenější v tomto čísle
- Methylsulfonylmethane increases osteogenesis and regulates the mineralization of the matrix by transglutaminase 2 in SHED cells
- Oregano powder reduces Streptococcus and increases SCFA concentration in a mixed bacterial culture assay
- The characteristic of patulous eustachian tube patients diagnosed by the JOS diagnostic criteria
- Parametric CAD modeling for open source scientific hardware: Comparing OpenSCAD and FreeCAD Python scripts