Time-calibrated molecular phylogeny of pteropods

Evolution and biomineralization of pteropod shells

Shelled pteropods, known as sea butterflies, are a group of small gastropods that spend their entire lives swimming and drifting in the open ocean. They build thin shells of aragonite, a metastable polymorph of calcium carbonate. Pteropod shells have been shown to experience dissolution and reduced thickness with a decrease in pH and therefore represent valuable bioindicators to monitor the impacts of ocean acidification. Over the past decades, several studies have highlighted the striking diversity of shell microstructures in pteropods, with exceptional mechanical properties, but their evolution and future in acidified waters remains uncertain. Here, we revisit the body-of-work on pteropod biomineralization, focusing on shell microstructures and their evolution. The evolutionary history of pteropods was recently resolved, and thus it is timely to examine their shell microstructures in such context. We analyse new images of shells from fossils and recent species providing a comprehensive overview of their structural diversity. Pteropod shells are made of the crossed lamellar and prismatic microstructures common in molluscs, but also of curved nanofibers which are proposed to form a helical three-dimensional structure. Our analyses suggest that the curved fibres emerged before the split between coiled and uncoiled pteropods and that they form incomplete to multiple helical turns. The curved fibres are seen as an important trait in the adaptation to a planktonic lifestyle, giving maximum strength and flexibility to the pteropod thin and lightweight shells. Finally, we also elucidate on the candidate biomineralization genes underpinning the shell diversity in these important indicators of ocean health.

Oceanic dispersal barriers in a holoplanktonic gastropod

Pteropods, a group of holoplanktonic gastropods, are regarded as bioindicators of the effects of ocean acidification on open ocean ecosystems, because their thin aragonitic shells are susceptible to dissolution. While there have been recent efforts to address their capacity for physiological acclimation, it is also important to gain predictive understanding of their ability to adapt to future ocean conditions. However, little is known about the levels of genetic variation and large-scale population structuring of pteropods, key characteristics enabling local adaptation. We examined the spatial distribution of genetic diversity in the mitochondrial cytochrome c oxidase I (COI) and nuclear 28S gene fragments, as well as shell shape variation, across a latitudinal transect in the Atlantic Ocean (35°N-36°S) for the pteropod Limacina bulimoides. We observed high levels of genetic variability (COI π = 0.034, 28S π = 0.0021) and strong spatial structuring (COI ΦST  = 0.230, 28S ΦST  = 0.255) across this transect. Based on the congruence of mitochondrial and nuclear differentiation, as well as differences in shell shape, we identified a primary dispersal barrier in the southern Atlantic subtropical gyre (15-18°S). This barrier is maintained despite the presence of expatriates, a gyral current system, and in the absence of any distinct oceanographic gradients in this region, suggesting that reproductive isolation between these populations must be strong. A secondary dispersal barrier supported only by 28S pairwise ΦST comparisons was identified in the equatorial upwelling region (between 15°N and 4°S), which is concordant with barriers observed in other zooplankton species. Both oceanic dispersal barriers were congruent with regions of low abundance reported for a similar basin-scale transect that was sampled 2 years later. Our finding supports the hypothesis that low abundance indicates areas of suboptimal habitat that result in barriers to gene flow in widely distributed zooplankton species. Such species may in fact consist of several populations or (sub)species that are adapted to local environmental conditions, limiting their potential for adaptive responses to ocean changes. Future analyses of genome-wide diversity in pteropods could provide further insight into the strength, formation and maintenance of oceanic dispersal barriers.

The origin and diversification of pteropods precede past perturbations in the Earth's carbon cycle

Pteropods are a group of planktonic gastropods that are widely regarded as biological indicators for assessing the impacts of ocean acidification. Their aragonitic shells are highly sensitive to acute changes in ocean chemistry. However, to gain insight into their potential to adapt to current climate change, we need to accurately reconstruct their evolutionary history and assess their responses to past changes in the Earth's carbon cycle. Here, we resolve the phylogeny and timing of pteropod evolution with a phylogenomic dataset (2,654 genes) incorporating new data for 21 pteropod species and revised fossil evidence. In agreement with traditional taxonomy, we recovered molecular support for a division between "sea butterflies" (Thecosomata; mucus-web feeders) and "sea angels" (Gymnosomata; active predators). Molecular dating demonstrated that these two lineages diverged in the early Cretaceous, and that all main pteropod clades, including shelled, partially-shelled, and unshelled groups, diverged in the mid- to late Cretaceous. Hence, these clades originated prior to and subsequently survived major global change events, including the Paleocene-Eocene Thermal Maximum (PETM), the closest analog to modern-day ocean acidification and warming. Our findings indicate that planktonic aragonitic calcifiers have shown resilience to perturbations in the Earth's carbon cycle over evolutionary timescales.

Fossil-calibrated molecular phylogeny of atlantid heteropods (Gastropoda, Pterotracheoidea)

BACKGROUND

The aragonite shelled, planktonic gastropod family Atlantidae (shelled heteropods) is likely to be one of the first groups to be impacted by imminent ocean changes, including ocean warming and ocean acidification. With a fossil record spanning at least 100 Ma, atlantids have experienced and survived global-scale ocean changes and extinction events in the past. However, the diversification patterns and tempo of evolution in this family are largely unknown.

RESULTS

Based on a concatenated maximum likelihood phylogeny of three genes (cytochrome c oxidase subunit 1 mitochondrial DNA, 28S and 18S ribosomal rRNA) we show that the three extant genera of the family Atlantidae, Atlanta, Protatlanta and Oxygyrus, form monophyletic groups. The genus Atlanta is split into two groups, one exhibiting smaller, well ornamented shells, and the other having larger, less ornamented shells. The fossil record, in combination with a fossil-calibrated phylogeny, suggests that large scale atlantid extinction was accompanied by considerable and rapid diversification over the last 25 Ma, potentially driven by vicariance events.

CONCLUSIONS

Now confronted with a rapidly changing modern ocean, the ability of atlantids to survive past global change crises gives some optimism that they may be able to persist through the Anthropocene.

Novel genomic resources for shelled pteropods: a draft genome and target capture probes for Limacina bulimoides, tested for cross-species relevance

BACKGROUND

Pteropods are planktonic gastropods that are considered as bio-indicators to monitor impacts of ocean acidification on marine ecosystems. In order to gain insight into their adaptive potential to future environmental changes, it is critical to use adequate molecular tools to delimit species and population boundaries and to assess their genetic connectivity. We developed a set of target capture probes to investigate genetic variation across their large-sized genome using a population genomics approach. Target capture is less limited by DNA amount and quality than other genome-reduced representation protocols, and has the potential for application on closely related species based on probes designed from one species.

RESULTS

We generated the first draft genome of a pteropod, Limacina bulimoides, resulting in a fragmented assembly of 2.9 Gbp. Using this assembly and a transcriptome as a reference, we designed a set of 2899 genome-wide target capture probes for L. bulimoides. The set of probes includes 2812 single copy nuclear targets, the 28S rDNA sequence, ten mitochondrial genes, 35 candidate biomineralisation genes, and 41 non-coding regions. The capture reaction performed with these probes was highly efficient with 97% of the targets recovered on the focal species. A total of 137,938 single nucleotide polymorphism markers were obtained from the captured sequences across a test panel of nine individuals. The probes set was also tested on four related species: L. trochiformis, L. lesueurii, L. helicina, and Heliconoides inflatus, showing an exponential decrease in capture efficiency with increased genetic distance from the focal species. Sixty-two targets were sufficiently conserved to be recovered consistently across all five species.

CONCLUSION

The target capture protocol used in this study was effective in capturing genome-wide variation in the focal species L. bulimoides, suitable for population genomic analyses, while providing insights into conserved genomic regions in related species. The present study provides new genomic resources for pteropods and supports the use of target capture-based protocols to efficiently characterise genomic variation in small non-model organisms with large genomes.

Resolving species boundaries in the Atlanta brunnea species group (Gastropoda, Pterotracheoidea)

Atlantid heteropods are a family of holoplanktonic marine gastropods that occur primarily in tropical and subtropical latitudes. Atlantids bear a delicate aragonitic shell (<14 mm) and live in the upper ocean, where ocean acidification and ocean warming have a pronounced effect. Therefore, atlantids are likely to be sensitive to these ocean changes. However, we lack sufficiently detailed information on atlantid taxonomy and biogeography, which is needed to gain a deeper understanding of the consequences of a changing ocean. To date, atlantid taxonomy has mainly relied on morphometrics and shell ornamentation, but recent molecular work has highlighted hidden diversity. This study uses an integrated approach in a global analysis of biogeography, variation in shell morphology and molecular phylogenies based on three genes (CO1, 28S and 18S) to resolve the species boundaries within the Atlanta brunnea group. Results identify a new species, Atlanta vanderspoeli, from the Equatorial and South Pacific Ocean, and suggest that individuals of A. brunnea living in the Atlantic Ocean are an incipient species. Our results provide an important advance in atlantid taxonomy and will enable identification of these species in future studies of living and fossil plankton.

Assessing species boundaries in the open sea: an integrative taxonomic approach to the pteropod genus Diacavolinia

To track changes in pelagic biodiversity in response to climate change, it is essential to accurately define species boundaries. Shelled pteropods are a group of holoplanktonic gastropods that have been proposed as bio-indicators because of their vulnerability to ocean acidification. A particularly suitable, yet challenging group for integrative taxonomy is the pteropod genus Diacavolinia, which has a circumglobal distribution and is the most species-rich pteropod genus, with 24 described species. We assessed species boundaries in this genus, with inferences based on geometric morphometric analyses of shell-shape variation, genetic (cytochrome c oxidase subunit I, 28S rDNA sequences) and geographic data. We found support for a total of 13 species worldwide, with observations of 706 museum and 263 freshly collected specimens across a global collection of material, including holo‐ and paratype specimens for 14 species. In the Atlantic Ocean, two species are well supported, in contrast to the eight currently described, and in the Indo‐Pacific we found a maximum of 11 species, partially merging 13 of the described species. Distributions of these revised species are congruent with well-known biogeographic provinces. Combining varied datasets in an integrative framework may be suitable for many diverse taxa and is an important first step to predicting species-specific responses to global change.

The old and the new plankton: ecological replacement of associations of mollusc plankton and giant filter feeders after the Cretaceous?

Owing to their great diversity and abundance, ammonites and belemnites represented key elements in Mesozoic food webs. Because of their extreme ontogenetic size increase by up to three orders of magnitude, their position in the food webs likely changed during ontogeny. Here, we reconstruct the number of eggs laid by large adult females of these cephalopods and discuss developmental shifts in their ecologic roles. Based on similarities in conch morphology, size, habitat and abundance, we suggest that similar niches occupied in the Cretaceous by juvenile ammonites and belemnites were vacated during the extinction and later partially filled by holoplanktonic gastropods. As primary consumers, these extinct cephalopod groups were important constituents of the plankton and a principal food source for planktivorous organisms. As victims or, respectively, profiteers of this case of ecological replacement, filter feeding chondrichthyans and cetaceans likely filled the niches formerly occupied by large pachycormid fishes during the Jurassic and Cretaceous.

Time-calibrated molecular phylogeny of pteropods

Pteropods are a widespread group of holoplanktonic gastropod molluscs and are uniquely suitable for study of long-term evolutionary processes in the open ocean because they are the only living metazoan plankton with a good fossil record. Pteropods have been proposed as bioindicators to monitor the impacts of ocean acidification and in consequence have attracted considerable research interest, however, a robust evolutionary framework for the group is still lacking. Here we reconstruct their phylogenetic relationships and examine the evolutionary history of pteropods based on combined analyses of Cytochrome Oxidase I, 28S, and 18S ribosomal rRNA sequences and a molecular clock calibrated using fossils and the estimated timing of the formation of the Isthmus of Panama. Euthecosomes with uncoiled shells were monophyletic with Creseis as the earliest diverging lineage, estimated at 41-38 million years ago (mya). The coiled euthecosomes (Limacina, Heliconoides, Thielea) were not monophyletic contrary to the accepted morphology-based taxonomy; however, due to their high rate heterogeneity no firm conclusions can be drawn. We found strong support for monophyly of most euthecosome genera, but Clio appeared as a polyphyletic group, and Diacavolinia grouped within Cavolinia, making the latter genus paraphyletic. The highest evolutionary rates were observed in Heliconoides inflatus and Limacina bulimoides for both 28S and 18S partitions. Using a fossil-calibrated phylogeny that sets the first occurrence of coiled euthecosomes at 79-66 mya, we estimate that uncoiled euthecosomes evolved 51-42 mya and that most extant uncoiled genera originated 40-15 mya. These findings are congruent with a molecular clock analysis using the Isthmus of Panama formation as an independent calibration. Although not all phylogenetic relationships could be resolved based on three molecular markers, this study provides a useful resource to study pteropod diversity and provides general insight into the processes that generate and maintain their diversity in the open ocean.