Near-future pH conditions severely impact calcification, metabolism and the nervous system in the pteropod Heliconoides inflatus

Physiological and transcriptomic responses of Aurelia coerulea polyps to acidified seawater conditions

Scyphozoan jellyfish, known for their evolutionary position and ecological significance, are thought to exhibit relatively notable resilience to ocean acidification. However, knowledge regarding the molecular mechanisms underlying the scyphozoan jellyfish response to acidified seawater conditions is currently lacking. In this study, two independent experiments were conducted to determine the physiological and molecular responses of moon jellyfish (Aurelia coerulea) polyps to within- and trans-generational exposure to two reduced pH treatments (pH 7.8 and pH 7.6). The results revealed that the asexual reproduction of A. coerulea polyps significantly declined under acute exposure to pH 7.6 compared with that of polyps at ambient pH conditions. Transcriptomics revealed a notable upregulation of genes involved in immunity and cytoskeleton components. In contrast, genes associated with metabolism were downregulated in response to reduced pH treatments after 6 weeks of within-generational acidified conditions. However, reduced pH treatments had no significant influence on the asexual reproduction of A. coerulea polyps after exposure to acidified conditions over a total of five generations, suggesting that A. coerulea polyps may acclimate to low pH levels. Transcriptomics revealed distinct gene expression profiles between within- and trans-generational exposure groups to two reduced pH treatments. The offspring polyps of A. coerulea subjected to trans-generational acidified conditions exhibited both upregulated and downregulated expression of genes associated with metabolism. These physiological and transcriptomic characteristics of A. coerulea polyps in response to elevated CO2 levels suggest that polyps produced asexually under acidified conditions may be resilient to such conditions in the future.

Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica)

Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.

CO2-Induced Ocean Acidification Alters the Burrowing Behavior of Manila Clam Ruditapes philippinarum by Reversing GABAA Receptor Function

Biological burrowing behavior is an important driver shaping ecosystems that is being threatened by CO₂-induced ocean acidification; however, the effects of ocean acidification on burrowing behavior and its neurological mechanism remain unclear. This study showed that elevated pCO₂ significantly affected the burrowing behaviors of the Manila clam Ruditapes philippinarum, such as increased foot contraction, burrowing time, and intrabottom movement and decreased burrowing depth. Delving deeper into the mechanism, exposure to elevated pCO₂ significantly decreased extracellular pH and increased [HCO₃^-]. Moreover, an indicator GABA_A receptor, a neuroinhibitor for movement, was found to be closely associated with behavioral changes. In situ hybridization confirmed that the GABA_A receptor was widely distributed in ganglia and foot muscles, and elevated pCO₂ significantly increased the mRNA level and GABA concentration. However, the increase in GABA_A receptor and its ligand did not suppress the foot movement, but rather sent "excitatory" signals for foot contraction. The destabilization of acid-base homeostasis was demonstrated to induce an increase in the reversal potential for GABA_A receptor and an alteration in GABA_A receptor function under elevated pCO₂. This study revealed that elevated pCO₂ affects the burrowing behavior of Manila clams by altering GABA_A receptor function from inhibitory to excitatory.

High and diurnally fluctuating carbon dioxide exposure produces lower mercury toxicity in a marine copepod

Coastal waters have experienced fluctuations in partial pressure of carbon dioxide (pCO₂) and mercury (Hg) pollution, yet little is known concerning how natural pCO₂ fluctuations affect Hg biotoxicity. Here, a marine copepod Tigriopus japonicus was interactively exposed to different seawater pCO₂ (ambient 400, steady elevated 1000, and fluctuating elevated 1000 ± 600 μatm) scenarios and Hg (control, 2 μg/L) treatments for 7 d. The results showed that elevated pCO₂ decreased Hg bioaccumulation, and it was even more under fluctuating elevated pCO₂ condition. We found energy depletion and oxidative stress under Hg-treated copepods, while combined exposure initiated compensatory response to alleviate Hg toxicity. Intriguingly, fluctuating acidification presented more immune defense related genes/processes in Hg-treated copepods when compared to steady acidification, probably linking with the greater decrease in Hg bioaccumulation. Collectively, understanding how fluctuating acidification interacts with Hg contaminant will become more crucial in predicting their risks to coastal biota and ecosystems.

RNAi Silencing of the Biomineralization Gene Perlucin Impairs Oyster Ability to Cope with Ocean Acidification

Calcifying marine organisms, including the eastern oyster (Crassostrea virginica), are vulnerable to ocean acidification (OA) because it is more difficult to precipitate calcium carbonate (CaCO₃). Previous investigations of the molecular mechanisms associated with resilience to OA in C. virginica demonstrated significant differences in single nucleotide polymorphism and gene expression profiles among oysters reared under ambient and OA conditions. Converged evidence generated by both of these approaches highlighted the role of genes related to biomineralization, including perlucins. Here, gene silencing via RNA interference (RNAi) was used to evaluate the protective role of a perlucin gene under OA stress. Larvae were exposed to short dicer-substrate small interfering RNA (DsiRNA-perlucin) to silence the target gene or to one of two control treatments (control DsiRNA or seawater) before cultivation under OA (pH ~7.3) or ambient (pH ~8.2) conditions. Two transfection experiments were performed in parallel, one during fertilization and one during early larval development (6 h post-fertilization), before larval viability, size, development, and shell mineralization were monitored. Silenced oysters under acidification stress were the smallest, had shell abnormalities, and had significantly reduced shell mineralization, thereby suggesting that perlucin significantly helps larvae mitigate the effects of OA.

Integrative assessment of biomarker responses in Mytilus galloprovincialis exposed to seawater acidification and copper ions

The interactive effects of ocean acidification (OA) and copper (Cu) ions on the mussel Mytilus galloprovincialis are not well understood. The underlying mechanisms also remain obscure. In this study, individuals of M. galloprovincialis were exposed for 28 days to 25 μg/L and 50 μg/L Cu ions at two pH levels (ambient level - pH 8.1; acidified level - pH 7.6). The mussels were then monitored for 56 days to determine their recovery ability. Physiological parameters (clearance rate and respiration rate), oxidative stress and neurotoxicity biomarkers (activities of superoxide dismutase, lipid peroxidation, catalase, and acetylcholinesterase), as well as the recovery ability of these parameters, were investigated in two typical tissues (i.e., gills and digestive glands). Results showed that (1) OA affected the bioconcentration of Cu in the gills and digestive glands of the mussels; (2) both OA and Cu can lead to physiological disturbance, oxidative stress, cellular damage, energy metabolism disturbance, and neurotoxicity on M. galloprovincialis; (3) gill is more sensitive to OA and Cu than digestive gland; (4) Most of the biochemical and physiological alternations caused by Cu and OA exposures in M. galloprovincialis can be repaired by the recovery experiments; (5) integrated biomarker response (IBR) analysis demonstrated that both OA and Cu ions exposure caused survival stresses to the mussels, with the highest effect shown in the co-exposure treatment. This study highlights the necessity to include OA along with pollutants in future studies to better elucidate the risks of ecological perturbations. The work also sheds light on the recovery of marine animals after short-term environmental stresses when the natural environment has recovered.

Assessment of the juvenile vulnerability of symbiont-bearing giant clams to ocean acidification

Ocean acidification (OA) severely affects marine bivalves, especially their calcification processes. However, very little is known about the fate of symbiont-bearing giant clams in the acidified oceans, which hinders our ability to develop strategies to protect this ecologically and economically important group in coral reef ecosystems. Here, we explored the integrated juvenile responses of fluted giant clam Tridacna squamosa (Lamarck, 1819) to acidified seawater at different levels of biological organization. Our results revealed that OA did not cause a significant reduction in survival and shell growth performance, indicating that T. squamosa juveniles are tolerated to moderate acidification. Yet, significantly reduced net calcification rate demonstrated the calcifying physiology sensitivity to OA, in line with significant declines in symbiont photosynthetic yield and zooxanthellae density which in turn lowered the amount of energy supply for energetically expensive calcification processes. Subsequent transcriptome sequencing and comparative analysis of differentially expressed genes revealed that the regulation of calcification processes, such as transport of calcification substrates, acid-base regulation, synthesis of organic matrix in the calcifying fluid, as well as metabolic depression were the major response to OA. Taken together, the integration of physiological and molecular responses can provide a comprehensive understanding of how the early life history stages of giant clams respond to OA and make an important leap forward in assessing their fate under future ocean conditions.

The Marine Gastropod Crepidula fornicata Remains Resilient to Ocean Acidification Across Two Life History Stages

Rising atmospheric CO₂ reduces seawater pH causing ocean acidification (OA). Understanding how resilient marine organisms respond to OA may help predict how community dynamics will shift as CO₂ continues rising. The common slipper shell snail Crepidula fornicata is a marine gastropod native to eastern North America that has been a successful invader along the western European coastline and elsewhere. It has also been previously shown to be resilient to global change stressors. To examine the mechanisms underlying C. fornicata’s resilience to OA, we conducted two controlled laboratory experiments. First, we examined several phenotypes and genome-wide gene expression of C. fornicata in response to pH treatments (7.5, 7.6, and 8.0) throughout the larval stage and then tested how conditions experienced as larvae influenced juvenile stages (i.e., carry-over effects). Second, we examined genome-wide gene expression patterns of C. fornicata larvae in response to acute (4, 10, 24, and 48 h) pH treatment (7.5 and 8.0). Both C. fornicata larvae and juveniles exhibited resilience to OA and their gene expression responses highlight the role of transcriptome plasticity in this resilience. Larvae did not exhibit reduced growth under OA until they were at least 8 days old. These phenotypic effects were preceded by broad transcriptomic changes, which likely served as an acclimation mechanism for combating reduced pH conditions frequently experienced in littoral zones. Larvae reared in reduced pH conditions also took longer to become competent to metamorphose. In addition, while juvenile sizes at metamorphosis reflected larval rearing pH conditions, no carry-over effects on juvenile growth rates were observed. Transcriptomic analyses suggest increased metabolism under OA, which may indicate compensation in reduced pH environments. Transcriptomic analyses through time suggest that these energetic burdens experienced under OA eventually dissipate, allowing C. fornicata to reduce metabolic demands and acclimate to reduced pH. Carry-over effects from larval OA conditions were observed in juveniles; however, these effects were larger for more severe OA conditions and larvae reared in those conditions also demonstrated less transcriptome elasticity. This study highlights the importance of assessing the effects of OA across life history stages and demonstrates how transcriptomic plasticity may allow highly resilient organisms, like C. fornicata, to acclimate to reduced pH environments.

Projected near-future ocean acidification decreases mercury toxicity in marine copepods

Here, we examined the combinational effect of ocean acidification (OA) and mercury (Hg) in the planktonic copepod Pseudodiaptomus annandalei in cross-factored response to different pCO₂ (400, 800 μatm) and Hg (control, 1.0 and 2.5 μg/L) exposures for three generations (F0-F2), followed by single-generation recovery (F3) under clean condition. Several phenotypic traits and Hg accumulation were analyzed for F0-F3. Furthermore, shotgun-based quantitative proteomics was performed for F0 and F2. Our results showed that OA insignificantly influenced the traits. During F0-F2, combined exposure reduced Hg accumulation as compared with the counterpart Hg treatment, supporting the mitigating effect of OA on Hg toxicity in copepods. Proteomics analysis indicated that the copepods probably increased energy production/storage and stress response to ensure physiological resilience against OA. However, Hg induced many toxic events (e.g., energy depletion and degenerated organomorphogenesis/embryogenesis for F0; cell cycle arrest and detrimental stress-defense for F2), which were translated to the population-level adverse outcome, i.e., compromised growth/reproduction. Particularly, compensatory proteome response was identified (e.g., increased immune defense for F0; energetic compensation and enhanced embryogenesis for F2), accounting for a negative interaction between OA and Hg. Together, this study provides the molecular mechanisms behind the effects of OA and Hg pollution in marine copepods.

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.

The impacts of past, present and future ocean chemistry on predatory planktonic snails

The atlantid heteropods represent the only predatory, aragonite shelled zooplankton. Atlantid shell production is likely to be sensitive to ocean acidification (OA), and yet we know little about their mechanisms of calcification, or their response to changing ocean chemistry. Here, we present the first study into calcification and gene expression effects of short-term OA exposure on juvenile atlantids across three pH scenarios: mid-1960s, ambient and 2050 conditions. Calcification and gene expression indicate a distinct response to each treatment. Shell extension and shell volume were reduced from the mid-1960s to ambient conditions, suggesting that calcification is already limited in today's South Atlantic. However, shell extension increased from ambient to 2050 conditions. Genes involved in protein synthesis were consistently upregulated, whereas genes involved in organismal development were downregulated with decreasing pH. Biomineralization genes were upregulated in the mid-1960s and 2050 conditions, suggesting that any deviation from ambient carbonate chemistry causes stress, resulting in rapid shell growth. We conclude that atlantid calcification is likely to be negatively affected by future OA. However, we also found that plentiful food increased shell extension and shell thickness, and so synergistic factors are likely to impact the resilience of atlantids in an acidifying ocean.

The role of ligand-gated chloride channels in behavioural alterations at elevated CO2 in a cephalopod

Projected future carbon dioxide (CO2) levels in the ocean can alter marine animal behaviours. Disrupted functioning of γ-aminobutyric acid type A (GABAA) receptors (ligand-gated chloride channels) is suggested to underlie CO2-induced behavioural changes in fish. However, the mechanisms underlying behavioural changes in marine invertebrates are poorly understood. We pharmacologically tested the role of GABA-, glutamate-, acetylcholine- and dopamine-gated chloride channels in CO2-induced behavioural changes in a cephalopod, the two-toned pygmy squid (Idiosepius pygmaeus). We exposed squid to ambient (∼450 µatm) or elevated (∼1000 µatm) CO2 for 7 days. Squid were treated with sham, the GABAA receptor antagonist gabazine or the non-specific GABAA receptor antagonist picrotoxin, before measurement of conspecific-directed behaviours and activity levels upon mirror exposure. Elevated CO2 increased conspecific-directed attraction and aggression, as well as activity levels. For some CO2-affected behaviours, both gabazine and picrotoxin had a different effect at elevated compared with ambient CO2, providing robust support for the GABA hypothesis within cephalopods. In another behavioural trait, picrotoxin but not gabazine had a different effect in elevated compared with ambient CO2, providing the first pharmacological evidence, in fish and marine invertebrates, for altered functioning of ligand-gated chloride channels, other than the GABAAR, underlying CO2-induced behavioural changes. For some other behaviours, both gabazine and picrotoxin had a similar effect in elevated and ambient CO2, suggesting altered function of ligand-gated chloride channels was not responsible for these CO2-induced changes. Multiple mechanisms may be involved, which could explain the variability in the CO2 and drug treatment effects across behaviours.

Age dependent physiological tolerances explain population dynamics and distribution in the intertidal zone: A study with porcelain crabs

Population dynamics and their response to environmental stressors have been widely studied in intertidal organisms. However, how these dynamics and responses change with animal age have been largely ignored to date. Traditionally, it is assumed that younger organisms are more sensitive than adults to environmental stressors; under this perspective it could be predicted that fully grown organisms should be able to occupy the harsh upper limit of their intertidal habitat. However, in some intertidal Porcelain crabs the opposite distribution has been observed. Using Petrolisthes laevigatus, we tested the physiological tolerance of crabs of different sizes (i.e. age) and evaluated how this trait shapes population dynamics (distribution and small-scale migrations under different weather conditions). We determined the abundance and size distribution of P. laevigatus at the middle and upper intertidal levels during sunny and rainy days, finding that abundances decreased drastically and size distribution shifted to smaller individuals on rainy days. In the laboratory, survival and behavioural responses of individuals in water at 5, 10, 15 and 33 PSU salinities were evaluated. Young crabs were found in higher proportion in the upper intertidal while fully grown crabs (i.e. adults) mainly occupied the middle intertidal zone. Young crabs had a higher osmoregulatory capacity than adults, as they were better at regulating passive water uptake when challenged with diluted seawater. This was also correlated with a lower lethal salinity LC₅₀ in young crabs compared to adults. Behavioural trials showed that young crabs performed better escaping in both water and air, at intermediate and reduced salinities than adults. Therefore, weather influences small scale migrations from the upper to the lower intertidal zone, and this migration is also age-dependent, with younger crabs being more tolerant to low salinities and therefore allowing them to remain in the upper intertidal zone during raniny days.

Magnitude and Predictability of pH Fluctuations Shape Plastic Responses to Ocean Acidification

AbstractPhenotypic plasticity is expected to facilitate the persistence of natural populations as global change progresses. The attributes of fluctuating environments that favor the evolution of plasticity have received extensive theoretical investigation, yet empirical validation of these findings is still in its infancy. Here, we combine high-resolution environmental data with a laboratory-based experiment to explore the influence of habitat pH fluctuation dynamics on the plasticity of gene expression in two populations of the Mediterranean mussel, Mytilus galloprovincialis. We linked differences in the magnitude and predictability of pH fluctuations in two habitats to population-specific gene expression profiles in ambient and stressful pH treatments. Our results demonstrate population-based differentiation in gene expression plasticity, whereby mussels native to a habitat exhibiting a large magnitude of pH fluctuations with low predictability display reduced phenotypic plasticity between experimentally imposed pH treatments. This work validates recent theoretical findings on evolution in fluctuating environments, suggesting that the predictability of fluctuating selection pressures may play a predominant role in shaping the phenotypic variation observed across natural populations.

Pteropods make thinner shells in the upwelling region of the California Current Ecosystem

Shelled pteropods are widely regarded as bioindicators for ocean acidification, because their fragile aragonite shells are susceptible to increasing ocean acidity. While short-term incubations have demonstrated that pteropod calcification is negatively impacted by ocean acidification, we know little about net calcification in response to varying ocean conditions in natural populations. Here, we examine in situ calcification of Limacina helicina pteropods collected from the California Current Ecosystem, a coastal upwelling system with strong spatial gradients in ocean carbonate chemistry, dissolved oxygen and temperature. Depth-averaged pH ranged from 8.03 in warmer offshore waters to 7.77 in cold CO₂-rich waters nearshore. Based on high-resolution micro-CT technology, we showed that shell thickness declined by ~ 37% along the upwelling gradient from offshore to nearshore water. Dissolution marks covered only ~ 2% of the shell surface area and were not associated with the observed variation in shell thickness. We thus infer that pteropods make thinner shells where upwelling brings more acidified and colder waters to the surface. Probably the thinner shells do not result from enhanced dissolution, but are due to a decline in calcification. Reduced calcification of pteropods is likely to have major ecological and biogeochemical implications for the cycling of calcium carbonate in the oceans.

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.

Transcriptional changes of Pacific oyster Crassostrea gigas reveal essential role of calcium signal pathway in response to CO2-driven acidification

There is increasing evidence that ocean acidification (OA) has a significant impact on marine organisms. However, the ability of most marine organisms to acclimate to OA and the underlying mechanisms are still not well understood. In the present study, whole transcriptome analysis was performed to compare the impacts of short- (7 days, named as short group) and long- (60 days, named as long group) term CO₂ exposure (pH 7.50) on Pacific oyster Crassostrea gigas. The responses of C. gigas to short- and long-term CO₂ exposure shared common mechanisms in metabolism, membrane-associated transportation and binding processes. Long-term CO₂ exposure induced significant expression of genes involved in DNA or RNA binding, indicating the activated transcription after long-term CO₂ exposure. Oysters in the short-term group underwent significant intracellular calcium variation and oxidative stress. In contrast, the intracellular calcium, ROS level in hemocytes and H₂O₂ in serum recovered to normal levels after long-term CO₂ exposure, suggesting the compensation of physiological status and mutual interplay between calcium and oxidative level. The compensation was supported by the up-regulation of a series of calcium binding proteins (CBPs) and calmodulins (CaMs) related signal pathway. The results provided valuable information to understand the molecular mechanism underlying the responses of Pacific oyster to the acidified ocean and might have implications for predicting the possible effects of global climate changes on oyster aquaculture.

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.

Combined stress of ocean acidification and warming influence survival and drives differential gene expression patterns in the Antarctic pteropod, Limacina helicina antarctica

The ecologically important thecosome pteropods in the Limacina spp. complex have recently been the focus of studies examining the impacts global change factors - e.g., ocean acidification (OA) and ocean warming (OW) - on their performance and physiology. This focus is driven by conservation concerns where the health of pteropod populations is threatened by the high susceptibility of their shells to dissolution in low aragonite saturation states associated with OA and how coupling of these stressors may push pteropods past the limits of physiological plasticity. In this manipulation experiment, we describe changes in the transcriptome of the Antarctic pteropod, Limacina helicina antarctica, to these combined stressors. The conditions used in the laboratory treatments met or exceeded those projected for the Southern Ocean by the year 2100. We made two general observations regarding the outcome of the data: (1) Temperature was more influential than pH in terms of changing patterns of gene expression, and (2) these Antarctic pteropods appeared to have a significant degree of transcriptomic plasticity to respond to acute abiotic stress in the laboratory. In general, differential gene expression was observed amongst the treatments; here, for example, transcripts associated with maintaining protein structure and cell proliferation were up-regulated. To disentangle the effects of OA and OW, we used a weighted gene co-expression network analysis to explore patterns of change in the transcriptome. This approach identified gene networks associated with OW that were enriched for transcripts proposed to be involved in increasing membrane fluidity at warmer temperatures. Together these data provide evidence that L.h.antarctica has a limited capacity to acclimate to the combined conditions of OA and OW used in this study. This reduced scope of acclimation argues for continued study of how adaptation to polar aquatic environments may limit the plasticity of present-day populations in responding to future environmental change.

Near-future ocean warming and acidification alter foraging behaviour, locomotion, and metabolic rate in a keystone marine mollusc

Environmentally-induced changes in fitness are mediated by direct effects on physiology and behaviour, which are tightly linked. We investigated how predicted ocean warming (OW) and acidification (OA) affect key ecological behaviours (locomotion speed and foraging success) and metabolic rate of a keystone marine mollusc, the sea hare Stylocheilus striatus, a specialist grazer of the toxic cyanobacterium Lyngbya majuscula. We acclimated sea hares to OW and/or OA across three developmental stages (metamorphic, juvenile, and adult) or as adults only, and compare these to sea hares maintained under current-day conditions. Generally, locomotion speed and time to locate food were reduced ~1.5- to 2-fold when the stressors (OW or OA) were experienced in isolation, but reduced ~3-fold when combined. Decision-making was also severely altered, with correct foraging choice nearly 40% lower under combined stressors. Metabolic rate appeared to acclimate to the stressors in isolation, but was significantly elevated under combined stressors. Overall, sea hares that developed under OW and/or OA exhibited a less severe impact, indicating beneficial phenotypic plasticity. Reduced foraging success coupled with increased metabolic demands may impact fitness in this species and highlight potentially large ecological consequences under unabated OW and OA, namely in regulating toxic cyanobacteria blooms on coral reefs.

CO2-driven ocean acidification weakens mussel shell defense capacity and induces global molecular compensatory responses

Oceanic uptake of atmospheric CO₂ is reducing seawater pH and shifting carbonate chemistry within, a process termed as ocean acidification (OA). Marine mussels are a family of ecologically and economically significant bivalves that are widely distributed along coastal areas worldwide. Studies have demonstrated that OA greatly disrupts mussels' physiological functions. However, the underlying molecular responses (e.g., whether there were any molecular compensation mechanisms) and the extent to which OA affects mussel shell defense capacity remain largely unknown. In this study, the thick shell mussels Mytilus coruscus were exposed to the ambient pH (8.1) or one of two lowered pH levels (7.8 and 7.4) for 40 days. The results suggest that future OA will damage shell structure and weaken shell strength and shell closure strength, ultimately reducing mussel shell defense capacity. In addition, future OA will also disrupt haemolymph pH and Ca²⁺ homeostasis, leading to extracellular acidosis and Ca²⁺ deficiency. Mantle transcriptome analyses indicate that mussels will adopt a series of molecular compensatory responses to mitigate these adverse effects; nevertheless, weakened shell defense capacity will increase mussels' susceptibility to predators, parasites and pathogens, and thereby reduce their fitness. Overall, the findings of this study have significant ecological and economic implications, and will enhance our understanding of the future of the mussel aquaculture industry and coastal ecosystems.

Ocean acidification promotes broad transcriptomic responses in marine metazoans: a literature survey

For nearly a decade, the metazoan-focused research community has explored the impacts of ocean acidification (OA) on marine animals, noting that changes in ocean chemistry can impact calcification, metabolism, acid-base regulation, stress response and behavior in organisms that hold high ecological and economic value. Because OA interacts with several key physiological processes in marine organisms, transcriptomics has become a widely-used method to characterize whole organism responses on a molecular level as well as inform mechanisms that explain changes in phenotypes observed in response to OA. In the past decade, there has been a notable rise in studies that examine transcriptomic responses to OA in marine metazoans, and here we attempt to summarize key findings across these studies. We find that organisms vary dramatically in their transcriptomic responses to pH although common patterns are often observed, including shifts in acid-base ion regulation, metabolic processes, calcification and stress response mechanisms. We also see a rise in transcriptomic studies examining organismal response to OA in a multi-stressor context, often reporting synergistic effects of OA and temperature. In addition, there is an increase in studies that use transcriptomics to examine the evolutionary potential of organisms to adapt to OA conditions in the future through population and transgenerational experiments. Overall, the literature reveals complex organismal responses to OA, in which some organisms will face more dramatic consequences than others. This will have wide-reaching impacts on ocean communities and ecosystems as a whole.

Ocean Acidification and Coastal Marine Invertebrates: Tracking CO2 Effects from Seawater to the Cell

In the last few decades, numerous studies have investigated the impacts of simulated ocean acidification on marine species and communities, particularly those inhabiting dynamic coastal systems. Despite these research efforts, there are many gaps in our understanding, particularly with respect to physiological mechanisms that lead to pathologies. In this review, we trace how carbonate system disturbances propagate from the coastal environment into marine invertebrates and highlight mechanistic links between these disturbances and organism function. We also point toward several processes related to basic invertebrate biology that are severely understudied and prevent an accurate understanding of how carbonate system dynamics influence organismic homeostasis and fitness-related traits. We recommend that significant research effort be directed to studying cellular phenotypes of invertebrates acclimated or adapted to elevated seawater pCO₂ using biochemical and physiological methods.

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.

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.

Condition of pteropod shells near a volcanic CO2 vent region

Natural gradients of pH in the ocean are useful analogues for studying the projected impacts of Ocean Acidification (OA) on marine ecosystems. Here we document the in situ impact of submarine CO₂ volcanic emissions (CO₂ vents) on live shelled-pteropods (planktonic gastropods) species Creseis conica in the Gulf of Naples (Tyrrhenian Sea, Mediterranean). Since the currents inside the Gulf will likely drive those pelagic calcifying organisms into and out of the CO₂ vent zones, we assume that pteropods will be occasionally exposed to the vents during their life cycle. Shell degradation and biomass were investigated in the stations located within and nearby the CO₂ vent emission in relation to the variability of sea water carbonate chemistry. A relative decrease in shell biomass (22%), increase in incidence of shell fractures (38%) and extent of dissolution were observed in Creseis conica collected in the Gulf of Naples compared to those from the Northern Tyrrhenian Sea (control stations). These results suggest that discontinuous but recurrent exposure to highly variable carbonate chemistry could consistently affect the characteristic of the pteropod shells.

Seasonal transcriptomes of the Antarctic pteropod, Limacina helicina antarctica

High latitude seas will be among the first marine systems to be impacted by ocean acidification (OA). Previous research studying the effects of OA on the pteropod, Limacina helicina antarctica, has led this species to be identified as a sentinel organism for OA in polar oceans. Here, we present gene expression data on L. h. antarctica, collected in situ during the seasonal transition from early spring to early summer. Our findings suggest that after over-wintering under seasonal sea ice, pteropods progress toward full maturity in the early summer when food becomes increasingly available. This progression is highlighted by a dramatic shift in gene expression that supports the development of cytoskeletal structures, membrane ion transportation, and metabolically important enzymes associated with glycolysis. In addition, we observed signs of defense of genomic integrity and maturation as evidenced by an up-regulation of genes involved in DNA replication, DNA repair, and gametogenesis. These data contribute to a broader understanding of the life-cycle dynamics for L. h. antarctica and provide key insights into the transcriptomic signals of pteropod maturation and growth during this key seasonal transition.

Comparative analysis of the Pocillopora damicornis genome highlights role of immune system in coral evolution

Comparative analysis of the expanding genomic resources for scleractinian corals may provide insights into the evolution of these organisms, with implications for their continued persistence under global climate change. Here, we sequenced and annotated the genome of Pocillopora damicornis, one of the most abundant and widespread corals in the world. We compared this genome, based on protein-coding gene orthology, with other publicly available coral genomes (Cnidaria, Anthozoa, Scleractinia), as well as genomes from other anthozoan groups (Actiniaria, Corallimorpharia), and two basal metazoan outgroup phlya (Porifera, Ctenophora). We found that 46.6% of P. damicornis genes had orthologs in all other scleractinians, defining a coral ‘core’ genome enriched in basic housekeeping functions. Of these core genes, 3.7% were unique to scleractinians and were enriched in immune functionality, suggesting an important role of the immune system in coral evolution. Genes occurring only in P. damicornis were enriched in cellular signaling and stress response pathways, and we found similar immune-related gene family expansions in each coral species, indicating that immune system diversification may be a prominent feature of scleractinian coral evolution at multiple taxonomic levels. Diversification of the immune gene repertoire may underlie scleractinian adaptations to symbiosis, pathogen interactions, and environmental stress.

Predatory strategies and behaviours in cephalopods are altered by elevated CO2

There is increasing evidence that projected near-future carbon dioxide (CO₂ ) levels can alter predator avoidance behaviour in marine invertebrates, yet little is known about the possible effects on predatory behaviours. Here we tested the effects of elevated CO₂ on the predatory behaviours of two ecologically distinct cephalopod species, the pygmy squid, Idiosepius pygmaeus, and the bigfin reef squid, Sepioteuthis lessoniana. Both species exhibited an increased latency to attack and altered body pattern choice during the attack sequence at elevated CO₂ . I. pygmaeus also exhibited a 20% decrease in predation rate, an increased striking distance, and reduced preference for attacking the posterior end of prey at elevated CO₂ . Elevated CO₂ increased activity levels of S. lessoniana comparable to those previously shown in I. pygmaeus, which could adversely affect their energy budget and increase their potential to be preyed upon. The effects of elevated CO₂ on predatory behaviours, predation strategies and activity levels of cephalopods reported here could have far-reaching consequences in marine ecosystems due to the ecological importance of cephalopods in the marine food web.

Exposure to CO2 influences metabolism, calcification and gene expression of the thecosome pteropod Limacina retroversa

Thecosomatous pteropods, a group of aragonite shell-bearing zooplankton, are becoming an important sentinel organism for understanding the influence of ocean acidification on pelagic organisms. These animals show vulnerability to changing carbonate chemistry conditions, are geographically widespread, and are both biogeochemically and trophically important. The objective of this study was to determine how increasing duration and severity of CO₂ treatment influence the physiology of the thecosome Limacina retroversa, integrating both gene expression and organism-level (respiration and calcification) metrics. We exposed pteropods to over-saturated, near-saturated or under-saturated conditions and sampled individuals at 1, 3, 7, 14 and 21 days of exposure to test for the effect of duration. We found that calcification was affected by borderline and under-saturated conditions by week two, while respiration appeared to be more strongly influenced by an interaction between severity and duration of exposure, showing complex changes by one week of exposure. The organismal metrics were corroborated by specific gene expression responses, with increased expression of biomineralization-associated genes in the medium and high treatments throughout and complex changes in metabolic genes corresponding to both captivity and CO₂ treatment. Genes associated with other physiological processes such as lipid metabolism, neural function and ion pumping had complex responses, influenced by both duration and severity. Beyond these responses, our findings detail the captivity effects for these pelagic organisms, providing information to contextualize the conclusions of previous studies, and emphasizing a need for better culturing protocols.

Ocean acidification alters predator behaviour and reduces predation rate

Ocean acidification poses a range of threats to marine invertebrates; however, the emerging and likely widespread effects of rising carbon dioxide (CO₂) levels on marine invertebrate behaviour are still little understood. Here, we show that ocean acidification alters and impairs key ecological behaviours of the predatory cone snail Conus marmoreus Projected near-future seawater CO₂ levels (975 µatm) increased activity in this coral reef molluscivore more than threefold (from less than 4 to more than 12 mm min^-1) and decreased the time spent buried to less than one-third when compared with the present-day control conditions (390 µatm). Despite increasing activity, elevated CO₂ reduced predation rate during predator-prey interactions with control-treated humpbacked conch, Gibberulus gibberulus gibbosus; 60% of control predators successfully captured and consumed their prey, compared with only 10% of elevated CO₂ predators. The alteration of key ecological behaviours of predatory invertebrates by near-future ocean acidification could have potentially far-reaching implications for predator-prey interactions and trophic dynamics in marine ecosystems. Combined evidence that the behaviours of both species in this predator-prey relationship are altered by elevated CO₂ suggests food web interactions and ecosystem structure will become increasingly difficult to predict as ocean acidification advances over coming decades.

Transcriptomic response of the Antarctic pteropod Limacina helicina antarctica to ocean acidification

Background

Ocean acidification (OA), a change in ocean chemistry due to the absorption of atmospheric CO₂ into surface oceans, challenges biogenic calcification in many marine organisms. Ocean acidification is expected to rapidly progress in polar seas, with regions of the Southern Ocean expected to experience severe OA within decades. Biologically, the consequences of OA challenge calcification processes and impose an energetic cost.

Results

In order to better characterize the response of a polar calcifier to conditions of OA, we assessed differential gene expression in the Antarctic pteropod, Limacina helicina antarctica. Experimental levels of pCO₂ were chosen to create both contemporary pH conditions, and to mimic future pH expected in OA scenarios. Significant changes in the transcriptome were observed when juvenile L. h. antarctica were acclimated for 21 days to low-pH (7.71), mid-pH (7.9) or high-pH (8.13) conditions. Differential gene expression analysis of individuals maintained in the low-pH treatment identified down-regulation of genes involved in cytoskeletal structure, lipid transport, and metabolism. High pH exposure led to increased expression and enrichment for genes involved in shell formation, calcium ion binding, and DNA binding. Significant differential gene expression was observed in four major cellular and physiological processes: shell formation, the cellular stress response, metabolism, and neural function. Across these functional groups, exposure to conditions that mimic ocean acidification led to rapid suppression of gene expression.

Conclusions

Results of this study demonstrated that the transcriptome of the juvenile pteropod, L. h. antarctica, was dynamic and changed in response to different levels of pCO₂. In a global change context, exposure of L. h. antarctica to the low pH, high pCO₂ OA conditions resulted in a suppression of transcripts for genes involved in key physiological processes: calcification, metabolism, and the cellular stress response. The transcriptomic response at both acute and longer-term acclimation time frames indicated that contemporary L. h. antarctica may not have the physiological plasticity necessary for adaptation to OA conditions expected in future decades. Lastly, the differential gene expression results further support the role of shelled pteropods such as L. h. antarctica as sentinel organisms for the impacts of ocean acidification.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4161-0) contains supplementary material, which is available to authorized users.

Responses of neurogenesis and neuroplasticity related genes to elevated CO2 levels in the brain of three teleost species

The continuous increase of anthropogenic CO₂ in the atmosphere resulting in ocean acidification has been reported to affect brain function in some fishes. During adulthood, cell proliferation is fundamental for fish brain growth and for it to adapt in response to external stimuli, such as environmental changes. Here we report the first expression study of genes regulating neurogenesis and neuroplasticity in brains of three-spined stickleback (Gasterosteus aculeatus), cinnamon anemonefish (Amphiprion melanopus) and spiny damselfish (Acanthochromis polyacanthus) exposed to elevated CO₂ The mRNA expression levels of the neurogenic differentiation factor (NeuroD) and doublecortin (DCX) were upregulated in three-spined stickleback exposed to high-CO₂ compared with controls, while no changes were detected in the other species. The mRNA expression levels of the proliferating cell nuclear antigen (PCNA) and the brain-derived neurotrophic factor (BDNF) remained unaffected in the high-CO₂ exposed groups compared to the control in all three species. These results indicate a species-specific regulation of genes involved in neurogenesis in response to elevated ambient CO₂ levels. The higher expression of NeuroD and DCX mRNA transcripts in the brain of high-CO₂-exposed three-spined stickleback, together with the lack of effects on mRNA levels in cinnamon anemonefish and spiny damselfish, indicate differences in coping mechanisms among fish in response to the predicted-future CO₂ level.

Exposure history determines pteropod vulnerability to ocean acidification along the US West Coast

The pteropod Limacina helicina frequently experiences seasonal exposure to corrosive conditions (Ω_ar < 1) along the US West Coast and is recognized as one of the species most susceptible to ocean acidification (OA). Yet, little is known about their capacity to acclimatize to such conditions. We collected pteropods in the California Current Ecosystem (CCE) that differed in the severity of exposure to Ω_ar conditions in the natural environment. Combining field observations, high-CO₂ perturbation experiment results, and retrospective ocean transport simulations, we investigated biological responses based on histories of magnitude and duration of exposure to Ω_ar < 1. Our results suggest that both exposure magnitude and duration affect pteropod responses in the natural environment. However, observed declines in calcification performance and survival probability under high CO₂ experimental conditions do not show acclimatization capacity or physiological tolerance related to history of exposure to corrosive conditions. Pteropods from the coastal CCE appear to be at or near the limit of their physiological capacity, and consequently, are already at extinction risk under projected acceleration of OA over the next 30 years. Our results demonstrate that Ω_ar exposure history largely determines pteropod response to experimental conditions and is essential to the interpretation of biological observations and experimental results.

Soluble adenylyl cyclase mediates mitochondrial pathway of apoptosis and ATP metabolism in oyster Crassostrea gigas exposed to elevated CO2

Ocean acidification (OA) has deleterious impacts on immune response and energy homeostasis status of Mollusca. In the present study, the apoptosis ratio of hemocytes and the adenosine triphosphate (ATP) allocation in gill tissues were determined after Pacific oysters Crassostrea gigas were exposed to elevated CO₂ environment (pH = 7.50) for 16 days.The apoptosis ratio in CO₂ exposure group (35.2%) was significantly higher (p < 0.05) than that in the control group, and the increased apoptosis ratio induced by elevated CO₂ could be significantly inhibited (p < 0.05) by KH7, a specific inhibitor of a bicarbonate sensor soluble adenylyl cyclase (sAC). After CO₂ exposure, sAC in oyster (CgsAC) was found to be clustered with mitochondria in the cytoplasm, and the pro-caspase-3 was cleaved into two small fragments. Moreover, the activities of caspase-3 and caspase-9 also increased post CO₂ exposure and these increases could be inhibited by KH7. However, the activities of caspase-8 did not change significantly compared with that in the control group. After CO₂ exposure, the ATP content in the gill increased significantly (p < 0.05) and such increase could also be inhibited by KH7. The ATP content in purified gill mitochondria decreased significantly (p < 0.05) after CO₂ exposure, which was also inhibited by KH7. These results implied that the elevated CO₂ could activate the mitochondria-CgsAC pathway of apoptosis and ATP metabolism in oyster, and this pathway played essential roles in maintaining the homeostasis and the balance of energy metabolism in response to OA.

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.

Revision of the genus Cuvierina Boas, 1886 based on integrative taxonomic data, including the description of a new species from the Pacific Ocean (Gastropoda, Thecosomata)

Shelled pteropods (Gastropoda, Thecosomata, Euthecosomata) are a group of holoplanktonic gastropods that occur predominantly in the surface layers of the world’s oceans. Accurate species identifications are essential for tracking changes in species assemblages of planktonic gastropods, because different species are expected to have different sensitivities to ocean changes. The genus Cuvierina has a worldwide warm water distribution pattern between ~36°N and ~39°S. Based on an integrative taxonomic approach combining morphometric, genetic, and biogeographic information, the two subgenera of Cuvierina, Cuvierinas. str. and Urceolarica, are rejected. A new species is introduced: Cuvierinatsudaisp. n., which has to date been considered the same species as Cuvierinapacifica. Cuvierinatsudaisp. n. is endemic to the Pacific Ocean and is characterised by a shell height of 7.2-8.0 mm, a moderately cylindrical shell shape, the absence of micro-ornamentation and a triangular aperture. Cuvierinapacifica is restricted to the centre of the oligotrophic southern Pacific gyre, has a shell height of 6.6-8.5 mm, a more cylindrical shell shape, no micro-ornamentation and a less triangular aperture than Cuvierinatsudaisp. n.