Evidence for Carbonate System Mediated Shape Shift in an Intertidal Predatory Gastropod

Phenotypic plasticity represents an important first-line organism response to newly introduced or changing environmental constraints. Knowledge about structural responses to environmental stressors could thus be an essential measure to predict species and ecosystem responses to a world in change. In this study, we combined morphometric analyses with environmental modelling to identify direct shape responses of the predatory gastropod Nucella lapillus to large-scale variability in sea surface temperature and the carbonate system. Our models suggest that the state of the carbonate system and, more specifically, the substrate inhibitor ratio $\left(\left[HC{O}_3^{-}\right]{\left[H^{+}\right]}^{-1}\right)$ (SIR) has a dominant effect on the shell shape of this intertidal muricid. Populations in regions with a lower SIR tend to form narrower shells with a higher spire to body whorl ratio, whereas populations in areas with a higher SIR form wider shells with a much lower spire to body whorl ratio. These results indicate that a widespread phenotypic response of N. lapillus to continuing ocean acidification can be expected, potentially altering the phenotypic response pattern to predator or wave exposure regimes with profound implications for North Atlantic rocky shore communities.

Introduction

Together with genetic predispositions, environmental constraints are the primary causes for shape variability within and among populations. Phenotypic plasticity is the ability of a single genotype to produce a range of different forms in response to its environment (Bradshaw, 1965; West-Eberhard, 1989; Sultan and Stearns, 2005). This trait, if adaptive, is generally perceived advantageous for an organism or population because it facilitates resistance to a range of different stressors allowing similar genotypes to populate heterogeneous habitats (Chevin and Hoffmann, 2017). Therefore, changes in morphology can be strong indicators for increasing or changing environmental stressors, knowledge about which may be important when monitoring the consequences of anthropogenically induced global climate change (Telesca et al., 2018). So far, studies have focused predominately on population and community composition responses (Ashton et al., 2017) or genetic markers to study the effect of global change. However, an understanding of phenotypic consequences was also found to be crucial to forming a complete picture of biological responses to a world in change (Hellberg et al., 2001; Cross et al., 2018; Telesca et al., 2018; Morley et al., 2019; Telesca et al., 2021).

Climate change associated with rising atmospheric carbon dioxide levels, increasing terrestrial and seawater temperatures, ocean acidification (OA) and more frequent extreme weather events is an emerging global threat to most of the Earth’s biosphere (Collins et al., 2013). In the oceans, annually increasing sea surface temperatures and OA pose strong stressors for marine organisms (Byrne, 2011; Fitzer et al., 2015). Species with calcium carbonate (CaCO₃) shells or skeletons are expected to be especially impacted by global climate change (Waldbusser et al., 2013; Waldbusser et al., 2014; Przeslawski et al., 2015; Suckling et al., 2015; Thomsen et al., 2015). The reduction in seawater pH associated with rising atmospheric CO₂ concentration results in potentially harmful effects (e.g., shell dissolution), making an organism more susceptible to predation, parasites, or infectious diseases (Gaylord et al., 2011; Viotti et al., 2019; Barclay et al., 2020). Marine calcifiers are reported to be especially vulnerable at the larval stage to the synergistic effects of OA and ocean warming (Przeslawski et al., 2015), as this can lead to premature death or severe retardation of shell formation (Gazeau et al., 2011; Thomsen et al., 2015; Waldbusser et al., 2015). To overcome these impacts, CaCO₃ bearing organisms are required to spend a larger proportion of their metabolic energy on calcification to maintain the status quo (Melzner et al., 2011), which could have far-reaching implications for their energy balance and competitiveness.

There is a consensus that calcification is linked to CaCO₃ substrate abundance and speciation, i.e., to the CaCO₃ saturation state Ω_CaCO3 in seawater, although its exact role in the acidified ocean is still debated (Bach, 2015; Cyronak et al., 2015; Waldbusser et al., 2016). The Ω_CaCO3 is a function of both calcium ion [Ca²⁺] and bicarbonate $\left[C{O}_3^{2-}\right]$ availability as well as the thermodynamic solubility product K_sp of the respective CaCO₃ polymorphs (for molluscs, most important are aragonite and calcite) at a specific pressure, salinity and temperature (Zeebe and Wolf-Gladrow, 2001) (eq. 1).

$\begin{array}{ll}{\text{Ω}}_{CaC{O}_3}=\frac{\left[C{a}^{2+}\right]\left[C{O}_3^{2-}\right]}{K_{sp}}& (1)\end{array}$

In seawater, both [Ca²⁺] and the dissolved inorganic carbon concentrations (DIC) are tightly coupled to seawater salinity, whereby the speciation of DIC in seawater is pH-dependent. This means that both a reduction in the salinity and the pH can reduce the saturation state of seawater. At Ω_CaCO3 < 1, CaCO₃ dissolution becomes thermodynamically favourable; however, as pointed out by Waldbusser et al. (2016), low Ω_CaCO3 states likely pose more of a kinetic rather than a thermodynamic constraint on marine calcifiers considering the rapid rate of calcification observed especially at the larval stage (Waldbusser et al., 2016). Irrespective of whether thermodynamic or kinetic bottlenecks are responsible, continuously low Ω_CaCO3 states can have serious implications for shell-forming organisms (Doney et al., 2009; Watson et al., 2012; Waldbusser et al., 2014). In recent years, many in-situ and laboratory studies have added to our understanding of the effect of climate change on marine calcifiers (Byrne, 2011; Byrne and Przeslawski, 2013; Kroeker et al., 2013; Parker et al., 2013; Fitzer et al., 2015; Kroeker et al., 2016; Cross et al., 2018). However, relatively short-term experiments (Hofmann et al., 2010), heterogeneous study designs, and partly conflicting results among studies have complicated interpretations and highlight the complexity of possible organism responses to a changing environment. While long-term, multi-generational experiments are key to studying phenotypic and genetic adaptions to environmental change (Hofmann et al., 2010), they require many months to years and thus pose a serious time constraint (Clark et al., 2019a; Clark et al., 2019b). This is especially true if organisms with long lifespans and/or a slow reproductive cycle are to be studied. Thus different approaches to investigating the effect of climate change in a natural setting that also account for both phenotypic and genetic adaptions receive increasingly more attention (Cross et al., 2018; Telesca et al., 2019; Bullard et al., 2021; Telesca et al., 2021). These types of studies focus on organism responses in the temporal domain by comparing archival collections with modern specimens collected from a strongly constrained spatial range (Fisher et al., 2009; Cross et al., 2018; Telesca et al., 2021) or in the spatial domain by studying the effect of environmental constraints on an organism over an extensive geospatial range, made possible by recent advances in ecological modelling (Telesca et al., 2018; Telesca et al., 2019).

The dog whelk, Nucella lapillus (Linnaeus, 1758), is an important predatory gastropod inhabiting the rocky intertidal zone of the North Atlantic (Burrows and Hughes, 1989; Hughes and Burrows, 1994). Nucella lapillus preys predominately on habitat-forming foundation species such as blue mussels or barnacles and thereby exerts a strong top-down control on rocky-shore ecosystems (Trussell et al., 2003; Benedetti-Cecchi and Trussell, 2014; Contolini et al., 2020). Its spatial distribution ranges from the northern tip of Norway to southern Portugal, spanning three climate zones. Nucella lapillus is a direct-developing gastropod without a larval stage (Moore, 1938). Its locomotor activity is limited to movements between prey (usually no more than a few centimetres) on which individuals are reported to feed for hours, if not days (Fretter and Graham, 1962). Nucella lapillus is only occasionally found below the tide marks and does not voluntarily crawl over sandy or muddy areas. The absence of a larval stage and the limited locomotor range of N. lapillus limits its dispersal so that distinct phenotypes exist within different allopatric populations, forming numerous discrete gamodemes (Berry and Crothers, 1974). Its calcareous shell consists of two layers, made of irregular calcite on the outside and crossed lamellar aragonite on the inside, which may be separated by a transitional spherulitic layer (Mayk, 2021). Site-specific exposure to wave action and predators are the most studied natural causes of shape variability among N. lapillus populations (Hughes and Elner, 1979; Vermeij and Currey, 1980; Palmer, 1990; Pascoal et al., 2012). On exposed shores, N. lapillus usually have short squat shells with a larger aperture size to accommodate a larger foot, which improves adherence to substrata. Squatter shells and better substratum adherence reduce the risk of dislodgment and are thus interpreted to be a morphological adaption to increase survival (Kitching et al., 1966; Guerra-Varela et al., 2009). Nucella lapillus on protected shores, where individuals are more regularly exposed to predators, tend to produce stronger and thicker shells (Currey and Hughes, 1982; Appleton and Palmer, 1988) and have increased growth rates (Burrows and Hughes, 1989). Individuals from protected shores also tend to exhibit relatively lighter soft tissue body weight as thicker shells provide comparably less space to accommodate the internal organs (Kitching et al., 1966). Site-specific exposure regimes have also been associated with genetic variability (karyotypes). Nucella lapillus show a pronounced chromosomal polymorphism; first reported in populations from the coast of Roscoff, France (Staiger, 1953). Known chromosome numbers range from 2n = 26 to 2n = 36 between populations (Bantock and Cockayne, 1975; Pascoe and Dixon, 1994). It was suggested that karyotypes correlate with local conditions, as the 2n = 26 karyotype appeared to be more common on exposed shores, whereas the 2n = 36 karyotype appeared to be more common on protected shores (Kirby et al., 1994). However, other studies did not confirm this finding, suggesting that shell shape in N. lapillus is a direct response to its physical and biological environment, expressed within a single generation and without the need for a change in karyotype (Pascoe, 2002; Pascoe, 2006). In any case, N. lapillus shows pronounced phenotypic adaptions to its habitat allowing this intertidal muricid to become a dominant top-down predator. New and continuous stressors, as brought about by OA and ocean warming, could, however, endanger its position, potentially jeopardising the trophic system on large parts of the North Atlantic rocky shore. While a wealth of studies exists in the literature that examines local effects on N. lapillus’ shell morphology, comparatively little is known about large-scale shape trends, despite its wide distribution. However, this would be of particular interest to understand how this important intertidal predatory gastropod adapts to different environments and has the potential to provide new insights and understanding of this species’ capacity to respond to a world in change. This paper aims to close this gap by investigating response patterns of N. lapillus across its entire latitudinal range with a particular focus on identifying significant relationships between environmental gradients and phenotypic responses which may indicate sensitivities of this intertidal muricid to environmental change. To do so, we analysed shape variance among N. lapillus populations, sampled over a 30 degrees latitudinal range, from the northern tip of Norway to the south of Portugal, and used multivariate ecological models to identify significant relationships of shell shape to global change relevant environmental parameters.

Materials and Methods

Shell Collection

Specimens used in this study were sampled from the extensive Crothers collection at the Natural History Museum, London, United Kingdom (NHM) (Table 1). These samples were collected by John Crothers and colleagues during the 1970s-1980s from numerous sites to study shell shape and colour changes with a special focus on wave exposure levels (Crothers, 1971; Crothers, 1979; Crothers, 1980; Crothers, 1981; Crothers, 1982; Crothers, 1983; Crothers, 1985a). From this collection, we sub-sampled a total of 1575 specimens from 19 sites spanning nearly the entire latitudinal range of N. lapillus on the European North Atlantic and Arctic coast (Figure 1). We only included specimens with intact apices that showed no extensive signs of erosion and that had likely reached sexual maturity during their lifetime (Galante-Oliveira et al., 2010).

TABLE 1

Table 1 Sampling details of the Nucella lapillus shells used in this study.

FIGURE 1

Figure 1 Sampling map of Nucella lapillus specimens and environmental gradients of sea surface temperature, seawater salinity, pH and Ω_Ca. The sampling sites ranged from 40°N to almost 70°N, spanning three distinct climate zones. Names of corresponding sampling site IDs can be obtained from Table 1. The colour gradient in the sampling map corresponds to that in the temperature plot, and all gradient colours are in accordance with Thyng et al. (2016).

Environmental Descriptors

To study the effect of environmental constraints on shell morphology, we used data from the global ocean- in-situ reprocessed carbon observations product (Copernicus Marine Service). This comprehensive global data set contains a compilation of in-situ observational gridded data collected from 1950 until 2020 for key parameters such as water temperature, salinity, alkalinity and pH from two up-to-date carbon and biogeochemistry data products, namely the Surface Ocean Carbon Atlas (SOCATv2020) and the Global Ocean Data Analysis Project (GLODAPv2.2020) that underwent rigorous data quality and bias checks (Bakker et al., 2016; Olsen et al., 2020). The data has a 1x1 degree spatial resolution at 33 different water depths. For this study, we used only sea surface data to reflect the conditions in the intertidal zone. Seawater conditions used in the further analysis were inferred from the nearest data grid to each sampling site. To have access to all carbonate system parameters, we calculated missing values using pH and alkalinity as input parameters to the function carb() in the R package seacarb v3.3.0 (Lavigne et al., 2008). The default dissociation constants (K1 and K2) for the respective temperature and salinity ranges were used in the calculation (Waters and Millero, 2013).

Morphometric Analysis

We measured standard dimensions, namely shell height, shell width, aperture height and aperture width, using digital Vernier callipers. Every dimension was measured three times per shell. Aperture size was estimated from height and width measurements using the equation for the area of an ellipse. Continuous shell shapes (outlines) were obtained using elliptic Fourier analysis (EFA) (Giardina and Kuhl, 1977; Kuhl and Giardina, 1982) using an adaptation of the methodology laid out elsewhere (Telesca et al., 2018). For this, shells were individually positioned on an illuminated platform with the aperture facing downward. Images were taken from a fixed distance with a DSLR Camera (Nikon 3000, Sigma 105 mm macro lens) mounted on a photo stand in a dark room. This produced images of the shells as black silhouettes with clear outlines, ideal for automatic outline tracing. We decided to analyse shell shape variability with the shell apertures facing downward in a natural position, rather than in the more common but unnatural aperture upwards position, to study shape changes that are more meaningful for N. lapillus in nature. Shell silhouette images were imported to R (R Core Team, 2020), and outlines were digitised and subsequently turned into closed shape polygons (defined by x-y coordinates) using the R package Momocs (Bonhomme et al., 2014). Shape outlines were smoothed using a running average smoothing algorithm to remove digitisation noise before each shell was individually aligned along its longitudinal axis. This was done by defining a secant line that connected the two furthest apart points (i.e., the shell anterior with the posterior) within the polygon shape. The shape was then rotated until the slope of the secant line was zero. Thereafter, shape polygons were centred and scaled with regard to orientation and size. An EFA with 10 harmonics encompassing 99% of the total harmonic power was computed on the outlines. Principal component analysis (PCA) was used to identify the axes of most variance among individuals (hereafter shape-PCs). For further analyses, we used the first five shape-PCs that encompassed more than 86% of the total shell shape variance among specimens as representative shape descriptors. A graphical representation of the shell shape analysis pipeline can be obtained from the Supplementary Material. To visualise the transition between extremes of the morphospace, we generated deformation grids (Thompson, 1992) and iso deformation lines (Bookstein, 1992) for each shape-PC. Individual relationships of the first five shell shape descriptors with sampling site latitude were analysed using simple linear regression analyses (Figure 2).

FIGURE 2

Figure 2 (A) Latitudinal trends of shell shape descriptors (shape-PC1 to shape-PC5) of the intertidal gastropod Nucella lapillus. Points represent the mean at every sampling site, and error bars denote the standard deviation. Linear regression significances are reported for each shape-PC (Significance levels: n.s. p > 0.05, *p < 0.05, **p < 0.001, ***p < 0.0001). (B) Deformation grids of shape-PC1 to shape-PC5 visualising the transition between extremes of the morphospace represented by each shape-PC (blue: mean - 3σ, red: mean + 3σ). Percentage contributions to total shape variance between individuals are reported for each shape-PC.

Data Analysis and Model Formation

The effect of key environmental parameters on shell shape variability was investigated using generalised additive mixed models (GAMM, mgcv) (Wood, 2011) to accommodate the hierarchical structure and non-linear relationships in our dataset. Initial data exploration was done following Zuur et al. (2010). No significant outliers were detected among shape-PCs. Boxplots of shape-PCs by sampling sites showed heterogeneous variances between shape-PCs which was an expected consequence of the PCA. Since the aim of this work was to study the effect of global change on shell formation, we considered a range of associated environmental parameters of the physical and chemical domain as model predictors. As stated earlier, Ω_CaCO3 states are indicators for spontaneous carbonate precipitation and, as such, have been shown to be good predictors of bio-calcification (Sanders et al., 2021). However, apart from the theoretical framework, there is little evidence supporting its mechanistic role in molluscan shell formation (Bach, 2015). There is growing evidence that hydrogencarbonate concentration $\left[HC{O}_3^{-}\right]$ rather than $\left[C{O}_3^{2-}\right]$ is likely the primary carbon species used in bio-calcification (Herfort et al., 2002; Allemand et al., 2004; Mackinder et al., 2010; Stumpp et al., 2012). Thus, we also included the substrate inhibitor ratio (SIR) along with the other parameters as possible predictors to describe variability in the carbonate system that is likely more mechanistically relevant to biological calcification (Jokiel, 2013; Bach, 2015; Cyronak et al., 2015). The SIR is a ratio of $\left[HC{O}_3^{-}\right]$ concentration by proton availability (H⁺) in mol µmol^-1 (eq. 2) and, as such, is a representative of the inorganic carbon substrate availability in relation to precipitation conditions.

$\begin{array}{ll}SIR=\frac{\left[HC{O}_3^{-}\right]}{\left[H^{+}\right]}& (2)\end{array}$

The predictors included in the initial model were: temperature, salinity, total alkalinity (A_T), pH, $\left[HC{O}_3^{-}\right]$, Ω_Ar, Ω_Ca, SIR, sampling site ID and shell height. We selected the first three shape-PCs that exhibited a significant correlation with latitude (Figure 2) to be included in the GAMM. Given that all three selected shape-PCs are shape-representatives of the same shell, normalised shape-PC1 to shape-PC3 descriptors were included simultaneously within the same GAMM to allow for a holistic description of the shell shape variability. The number of knots per descriptor was manually defined to constrain unnatural degrees of wiggliness, and cubic regression splines were used. The dredge() function in the R package MuMIn (Bartoń, 2020) was used to analyse all possible combinations of predictors, and the Akaike Information Criterion (AIC) was used as the principal tool in the model optimisation process. To overcome potential complications with collinearity in the model caused by correlated independent variables only models with an acceptable VIF factor < 3 of covariates were considered. The sampling site ID was included as a random effect term to account for the dependency of specimens from the same collection site. The best-supported model was then fitted using the restricted maximum likelihood approach (REML) and validated by inspection of standardised residual patterns to verify model assumptions. Variogram inspection suggested no significant spatial autocorrelation. The final model was of the form:

$\begin{array}{ll}\begin{array}{l}shapeP{C}_{ij}\sim N\left({\mu }_{ij};{\sigma }^2\right)\hfill \\ {\mu }_{ij}=f\left(temperatur{e}_i\right)\ast fshap{e}_j+f\left(alkalinit{y}_i\right)\ast fshap{e}_j+f\left(SI{R}_i\right)\ast \hfill \\ fshap{e}_j+f\left(shellheigh{t}_i\right)\ast fshap{e}_j+{\in }_i\hfill \\ {\in }_i\sim N\left(0;{\sigma }_i^2\right)\hfill \end{array}& (3)\end{array}$

where shapePC_ij represents the j^th observation for each shape-PC (j = 1-3) at every sampling site i (i = 1-19). f() * fshape denote the smoothing functions of temperature, alkalinity, SIR and shell height expressed by individual smoother for each level of j and ∈_i is the normally distributed random intercept. Effect size estimates of standardised predictors were calculated from the GAMM following Telesca et al., (2018). The standardisation of model predictors was done by subtracting the predictor’s mean from its standard deviation. 95% confidence intervals were constructed along mean effect sizes to compare the magnitude and significance of environmental predictors on shell shapes. All exploratory data analyses and modelling were performed in R (R Core Team, 2020) and all data and scripts are publicly available (see data availability statement).

Results

Mean plots of environmental parameters revealed clear gradients from north to south (Figure 1). As expected, mean sea surface temperature increased gradually towards the south. Maximum difference in average temperature between sampling sites was 11.49°C. Seawater salinity showed an increasing trend towards the south with pronounced local variability (Figure 1, sampling sites 7-9). Salinity variations between sampling sites were up to 3.14 PSU. Seawater pH decreased between 70°N to 50°N, but sampling sites south of 50°N exhibited higher seawater pH (Figure 1). The maximum seawater pH range encompassed between our sampling sites was 0.13. Minimum seawater Ω_Ca and Ω_Ar values at the sampling sites were 3 and 2 respectively. Seawater Ω_Ca states showed no apparent trend over the whole studied latitudinal range but a significant increase in saturation south of 50°N.

Relative contributions of the first five shape-PCs to the total morphospace range are shown in Figure 2B. Overall, shape variance explained by shape-PC1 (53.17%) was a widening of the shell. Higher shape-PC1 values described wider, more globular shells with almost rectangular shoulders and proportionally larger body whorls. Lower values described a narrowing of the shell with proportionally smaller body whorls and flatter shoulders. Shape-PC2 (14.76%) described a rounding of the body whorl. Higher shape-PC2 values described an overall rounding of the body whorl with a less contoured siphonal canal. Shape-PC3 (8.65%) described an “offsetting” of the body whorl to the spire, which gave the impression of tilting the shell’s longitudinal axis. Shell shape-PC4 (5.94%) and shape-PC5 (4%) only explained minor shape variations. Shape-PC4 described a slight widening and shift of the body whorl, especially around the shoulder, and shape-PC5 described contouring of the shell. A plot of among-population variation described by shape-PC1 and shape-PC2 showed a clear separation among populations, especially along the shape-PC1 axis (Figure 3). Mean plots of selected shell descriptors (shape-PC1 to shape-PC5) also showed a clear effect of locality and significant latitudinal trends. Linear regression analyses confirmed significant shell shape changes in relation to sampling site latitude as expressed by shape-PC1, shape-PC2 and shape-PC3. Shape-PC4 and shape-PC5 did change significantly with latitude (Figure 2).

FIGURE 3

Figure 3 Scatter plot of the first two shape-PCs from PCA analysis performed on elliptic Fourier coefficients revealing the among-population variation of outline and shape features between Nucella lapillus specimens. Confidence intervals are shown for each origin group (ellipses with dashed linetype), and corresponding morphospace extremes are represented as grey shapes.

The best-supported model (Table 2) revealed a significant relationship between the selected shape descriptors and environmental gradients. The model comparison revealed that sea surface temperature, A_T, SIR, shell height, and sampling site ID (as a random effect term) best described the observed changes in shape-PCs. The model identified highly significant non-linear relationships of shape-PC1 with sea surface temperature, A_T, SIR and shell height (Figure 4A). Shape-PC2 showed a non-significant relationship with sea surface temperature but highly significant relationships with A_T and SIR. Shell height had only a minor effect on shape-PC2 (Figure 4B). Shell shape variance described by shape-PC3 was also not associated with sea surface temperature, but showed significant correlations with A_T, SIR and shell height (Figure 4C). Mean effect size estimates (Figure 5) revealed differences in the relative contribution of environmental predictors on observed changes in shape. The SIR and shell height showed a significant effect on shell shape-PC1 to shape-PC3, and in addition, A_T exhibited a significant effect size for shape-PC2. Overall, effect size estimates primarily linked changes in the carbonate system as represented by A_T and especially the SIR to changes in shell shape. Most notably, we identified increasing SIR values resulted in the formation of wider, more globular shells with a more pronounced shoulder. Effect size estimates for sea surface temperature were non-significant for any of the three shape descriptors.

TABLE 2

Table 2 GAMM summary statistics of the best supported model.

FIGURE 4

Figure 4 GAMM shape trend predictions of the best-supported model as expressed by (A) shape-PC1, (B) shape-PC2, and (C) shape-PC3 with selected descriptors of sea surface temperature, A_T, substrate inhibitor ratio (SIR) and normalised shell height. Significance is reported for each predictor (Significance levels: n.s. p > 0.05, *p < 0.05, **p < 0.001, ***p < 0.0001). Colour gradients (red to blue) of significant predictor terms correspond to those in Figure 5.

FIGURE 5

Figure 5 (A) shape-PC1, (B) shape-PC2, and (C) shape-PC3 and the best-supported model’s mean effect size plots, including sea surface temperature, A_T, substrate inhibitor ratio (SIR) and shell height. Error bars represent the 95% CI, and the significance of either predictor is determined if the CI does not cross the zero line (p < 0.05). Asterisks mark significant predictors. Colour gradients (red to blue) correspond to the morphospace extremes visualised to the left of each effect size plot.

Discussion

Using a combined EFA and GAMM approach, we identified significant shape changes in the intertidal N. lapillus that followed latitudinal (Figure 2) and environmental trends (Figures 4, 5). We identified a significant narrowing of shells (i.e., a reduction in the shell aspect ratio) with relatively smaller body whorls from north to south. While we did not identify microclimate effects of individual sites on shell formation, which likely contributed to the observed considerable shape variation between individuals from the same population (Figure 2), our large scale trends are so clear that the variability of mean environmental conditions among sampling sites must have been significant enough to influence the shell formation in the studied specimens. Although the intertidal zone is known to exhibit pronounced temporal and spatial variability in the carbonate system (Krause-Jensen et al., 2015), our results suggest that a shift in the mean has a significant influence on shell formation. Effect size estimates of environmental predictors on shape-PCs revealed that differences in the carbonate system had the strongest influence of all environmental predictors on shell shape. In contrast, sea surface temperature showed a non-significant effect size for any of the three shape-PCs (Figure 5). Despite overall high Ω_Ca and Ω_Ar throughout the sampling sites, shell formation of N. lapillus appears to be remarkably sensitive to variations in the carbonate system. Our models indicate a significant link between shell formation of N. lapillus to variations in the “reactant” $\left[HC{O}_3^{-}\right]$ and “inhibitor” [H⁺] ratio (SIR) (Jokiel, 2011a; Jokiel, 2011b; Jokiel, 2013), which is expected to decrease in the coming years and centuries under increasing OA conditions. This suggests an overall shift in N. lapillus shell shape can be expected as a response or consequence to changing water chemistry constraints.

There is only a handful of studies explicitly investigating shell shape responses of marine gastropods to changes in the carbonate system, and these have produced partly conflicting results. A study of the shell aspect ratio in the gastropod Phorcus saucciatus collected from a CO₂ vent and a control site reported a similar shape response to those here, as specimens from the vent site tend to form narrower and more elongated shells (Viotti et al., 2019). However, seemingly conflicting results were reported from laboratory studies where the combined effect of increasing seawater temperature and decreasing pH resulted in the formation of more globular shells in the intertidal gastropod Littorina littorea and in juvenile N. lapillus (Melatunan et al., 2013; Rühl et al., 2017). This discrepancy among studies is not easily resolved but highlights the complexity of comparing species responses from laboratory and field studies. One of the biggest issues in comparing gastropod shape studies is the lack of a common vocabulary and methodology to describe shape responses. One prominent example here is the interchangeability of the words “globular”, “wide”, “spherical”, “squat”, “round” to describe seemingly similar shape changes for which a visual representation of morphospace extremes (Figures 2B and 5) would provide a better basis of comparison.

The discrepancy between study results also raises the interesting question of whether a universal shape response to changing water chemistry can be expected or if different species will respond in different ways to OA. A crucial point when shell shape responses to a variable carbonate system are to be compared is the careful separation between kinetically, thermodynamically, and stress-induced shape alterations, or in other words, if the phenotypic response is adaptive or non-adaptive. At the simplest level, changes in the seawater carbonate system as brought about by OA result in a reduction of the bio-calcification substrate accompanied by increases in [H⁺] abundance that may constrain crystal formation at the site of calcification through lowering Ω_CaCO3 (Waldbusser et al., 2013; Cyronak et al., 2015). This means that less material or substrate is available to the individual for shell formation and that calcification likely becomes metabolically more expensive as the individual is required to spend more energy on [H⁺] extrusion (Waldbusser et al., 2013; Tresguerres, 2016). While the latitudinal shape trends observed in our study are purely kinetically induced (in accordance with the kinetic–energetic hypothesis) (Waldbusser et al., 2016) and are therefore a direct indicator for variations in substrate and inhibitor abundance, we hypothesise that shape variation in OA laboratory studies is likely a mix of kinetic-, dissolution-, and stress-induced phenotypic responses precluding any direct comparisons.

The literature provides a wealth of information about local effects on N. lapillus shell shape. Shell aspect ratio variations in relation to local wave exposure regimes were first reported by (Cooke and Reed, 1895). The authors showed that N. lapillus on exposed shores exhibited smaller, squatter shells with a larger aspect ratio, whereas individuals from protected shores formed narrower, more elongated shells. This trend was later explained to be a response that improved fitness as smaller squatter shells and a proportionally larger foot of exposed individuals (obvious from relatively larger aperture sizes) conferred an advantage over more elongated shell forms with smaller foot sizes in high wave energy environments (Kitching and Ebling, 1967). However, the exposed morph showed a marked disadvantage on sheltered shores as animals of that shell form could barely retract into the shell, making them more vulnerable to predators and desiccation during low tide (Crothers, 1985b). Crothers (1973) found the shell height to aperture height ratio to be a convenient proxy for wave exposure trends for N. lapillus populations on Pembrokeshire, UK shores and thereby quantified the relationship between wave exposure and relative aperture size (Crothers, 1973). However, unlike the reported wave exposure induced shape changes, mean shape-PCs in this study showed no relationships with aperture size (Figure 6). Likewise, we also did not observe a change in apertural teeth expression across studied populations which would indicate variations in wave exposure (Crothers, 1971) or predatory pressure (Appleton and Palmer, 1988) constraints on shell deposition in the studied N. lapillus populations. This suggests that shape variations, as explained by shape-PC1 to shape-PC3 (over these large spatial scales), were largely decoupled from local wave exposure and predation gradients.

FIGURE 6

Figure 6 Correlation plot of mean aperture size vs mean shape-PCs grouped by sampling site. Error bars denote the standard deviation. The orange dashed line represents the 1:1 correlation line.

While shape variations are an excellent proxy to identify new or changing environmental stressors in a population or species (Marshall et al., 2015), predictions of the long term consequences of plastic adaptions are not trivial - or to use the words of R. J. Berry and J. H. Crothers: “Natural populations of animals rarely behave like suspensions of inanimate matter in a frictionless fluid, responding to and predictable by simple mathematical theory” (Berry and Crothers, 1974). However, newly introduced stressors and associated phenotypic adaptions always represent a trade-off between different optima and can thus challenge an organism’s position in an ecosystem (Ghalambor et al., 2007). Our findings suggest that a trend towards narrower, more elongated shells can be expected under increasing OA conditions with, perhaps, significant implications for the predator-prey relationship on North Atlantic rocky shores. Although shape-PCs were largely decoupled from local wave exposure gradients, a continuing trend towards the formation of narrower and more elongated shells could have conceivable consequences for populations on (extremely) exposed coasts. In particular, since global climate change is associated with increasingly severe storm events IPCC 2014, the larger surface of elongated shells could lead to more frequent dislodgement and, as a consequence, could make N. lapillus populations more susceptible to predatory attacks. The herein observed shell shape response to OA could lead to a disintegration of the delicate balance between shape responses to wave exposure or predatory pressure extremes with the potential consequence of a realignment of rocky-shore community structures. However, our results and predictions of morphological responses to end century OA conditions should be viewed with caution for three reasons: First, the site-specific intra-population variability is many times larger than the average shift in shape observed in this study (Figure 2), suggesting that there is a fair amount of tolerance involved in the shape response to environmental constraints. Second, a sensitivity analysis of shell morphology to different stressors is missing, which would be required to make robust predictions. Third, temporal and spatial constraints of environmental data used in this study limited fine-scale observations. While our models suggest a remarkable sensitivity of shell shape to variations in the carbonate system on a large geographic scale, it is not clear if shell morphology will be more sensitive to water chemistry or site-specific physical constraints under increasing OA conditions.

Irrespective of the direct consequences, our results show that changes in the carbonate system are likely to directly affect N. lapillus shell morphology. Although calcification conditions were, in classical terms, favourable at every sampling site, shell formation and resulting morphology were sensitive to changes in the SIR. Our findings uncovered a, so far unrecognised, large scale pattern in shell morphology of the intertidal N. lapillus, which highlights the coupling of phenotypic adaptions to variations in the carbonate system. However, although our findings suggest an adaptive trend in shell morphology under increasing OA conditions, further research is required to understand the phenotypic sensitivity to local and global constraints, which could present an exciting new research opportunity to understand shell formation constraints in a changing world.

Data Availability Statement

The code and data that support the findings of this study are openly available on Zenodo at https://doi.org/10.5281/zenodo.6520504.

Author Contributions

DM, LP, and EH conceived the original project and designed the study. DM performed the specimen collection, laboratory and modelling work, and analysed the data. DM wrote the first draft of the manuscript and all co-authors contributed substantially to revisions. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by a NERC studentship awarded to DM (NE/L002507/1).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank Andreia Salvador for access to the collections at the Natural History Museum (London).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2022.894182/full#supplementary-material

References

Allemand D., Ferrier-Pagès C., Furla P., Houlbrèque F., Puverel S., Reynaud S., et al. (2004). Biomineralisation in Reef-Building Corals: From Molecular Mechanisms to Environmental Control. C. R. Palevol. 3, 453–467. doi: 10.1016/j.crpv.2004.07.011

CrossRef Full Text | Google Scholar

Appleton R. D., Palmer A. R. (1988). Water-Borne Stimuli Released by Predatory Crabs and Damaged Prey Induce More Predator-Resistant Shells in a Marine Gastropod. Proc. Natl. Acad. Sci. U. S. A. 85, 4387–4391. doi: 10.1073/pnas.85.12.4387

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashton G. V., Morley S. A., Barnes D. K. A., Clark M. S., Peck L. S. (2017). Warming by 1°C Drives Species and Assemblage Level Responses in Antarctica’s Marine Shallows. Curr. Biol. 27, 2698–2705.e3. doi: 10.1016/j.cub.2017.07.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Bach L. T. (2015). Reconsidering the Role of Carbonate Ion Concentration in Calcification by Marine Organisms. Biogeosciences 12, 4939–4951. doi: 10.5194/bg-12-4939-2015

CrossRef Full Text | Google Scholar

Bakker D. C. E., Pfeil B., Landa C. S., Metzl N., O’Brien K. M., Olsen A., et al. (2016). A Multi-Decade Record of High-Quality Fco₂ Data in Version 3 of the Surface Ocean CO₂ Atlas (SOCAT). Earth Syst. Sci. Data 8, 383–413. doi: 10.5194/essd-8-383-2016

CrossRef Full Text | Google Scholar

Bantock C. R., Cockayne W. C. (1975). Chromosomal Polymorphism in Nucella lapillus. Hered. 34. Part 2, 231–245. doi: 10.1038/hdy.1975.26

CrossRef Full Text | Google Scholar

Barclay K. M., Gingras M. K., Packer S. T., Leighton L. R. (2020). The Role of Gastropod Shell Composition and Microstructure in Resisting Dissolution Caused by Ocean Acidification. Mar. Environ. Res. 162, 105105. doi: 10.1016/j.marenvres.2020.105105

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartoń K. (2020). MuMIn: multi-model inference. R package version 1.43.17

Google Scholar

Benedetti-Cecchi L., Trussell G. C. (2014). “Rocky Intertidal Communities,” in Marine Community Ecology and Conservation (Sunderland, MA: Sinauer Associates), 203–225.

Google Scholar

Berry R. J., Crothers J. H. (1974). Visible Variation in the Dog-Whelk, Nucella lapillus. J. Zool. 174, 123–148. doi: 10.1111/j.1469-7998.1974.tb03147.x

CrossRef Full Text | Google Scholar

Bonhomme V., Picq S., Gaucherel C., Claude J. (2014). Momocs: Outline Analysis Using R. J. Stat. Software 56, 1–24. doi: 10.18637/jss.v056.i13

CrossRef Full Text | Google Scholar

Bookstein F. L. (1992). Morphometric Tools for Landmark Data: Geometry and Biology (United Kingdom: Cambridge University Press). doi: 10.1017/CBO9780511573064

CrossRef Full Text | Google Scholar

Bradshaw A. D. (1965). Evolutionary Significance of Phenotypic Plasticity in Plants. Adv. Genet. 13, 115–155. doi: 10.1016/S0065-2660(08)60048-6

CrossRef Full Text | Google Scholar

Bullard E. M., Torres I., Ren T., Graeve O. A., Roy K. (2021). Shell Mineralogy of a Foundational Marine Species, Mytilus californianus, Over Half a Century in a Changing Ocean. Proc. Natl. Acad. Sci. U. S. A. 118. doi: 10.1073/pnas.2004769118

CrossRef Full Text | Google Scholar

Burrows M. T., Hughes R. N. (1989). Natural Foraging of the Dogwhelk, Nucella lapillus (Linnaeus); the Weather and Whether to Feed. J. Moll. Stud. 55, 285–295. doi: 10.1093/mollus/55.2.285

CrossRef Full Text | Google Scholar

Byrne M. (2011). Impact of Ocean Warming and Ocean Acidification on Marine Invertebrate Life History Stages: Vulnerabilities and Potential for Persistence in a Changing Ocean. Oceanogr. Mar. Biol. Annu. Rev. 49, 1–42. doi: 10.1201/b11009-2

CrossRef Full Text | Google Scholar

Byrne M., Przeslawski R. (2013). Multistressor Impacts of Warming and Acidification of the Ocean on Marine Invertebrates’ Life Histories. Integr. Comp. Biol. 53, 582–596. doi: 10.1093/icb/ict049

PubMed Abstract | CrossRef Full Text | Google Scholar

Chevin L.-M., Hoffmann A. A. (2017). Evolution of Phenotypic Plasticity in Extreme Environments. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 372. doi: 10.1098/rstb.2016.0138

CrossRef Full Text | Google Scholar

Clark M. S., Suckling C. C., Cavallo A., Mackenzie C. L., Thorne M. A. S., Davies A. J., et al. (2019a). Molecular Mechanisms Underpinning Transgenerational Plasticity in the Green Sea Urchin Psammech miliaris. Sci. Rep. 9, 952. doi: 10.1038/s41598-018-37255-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Clark M. S., Villota Nieva L., Hoffman J. I., Davies A. J., Trivedi U. H., Turner F., et al. (2019b). Lack of Long-Term Acclimation in Antarctic Encrusting Species Suggests Vulnerability to Warming. Nat. Commun. 10, 3383. doi: 10.1038/s41467-019-11348-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins M., Knutti R., Arblaster J., Dufresne J.-L., Fichefet T., Friedlingstein P., et al. (2013). Long-Term Climate Change: Projections, Commitments and Irreversibility in Climate Change 2013-the Physical Science Basis: Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge: Cambridge University Press). 1029–1136

Google Scholar

Contolini G. M., Reid K., Palkovacs E. P. (2020). Climate Shapes Population Variation in Dogwhelk Predation on Foundational Mussels. Oecologia 192, 553–564. doi: 10.1007/s00442-019-04591-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cooke A. H., Reed F. R. C. (1895). The Cambridge Natural History (London: Macmillan; Company).

Google Scholar

Cross E. L., Harper E. M., Peck L. S. (2018). A 120-Year Record of Resilience to Environmental Change in Brachiopods. Glob. Change Biol. 24, 2262–2271. doi: 10.1111/gcb.14085

CrossRef Full Text | Google Scholar

Crothers J. H. (1971). Further Observations on the Occurrence of “Teeth” in the Dog-Whelk Nucella lapillus. J. Mar. Biol. Assoc. U. K. 51, 623–639. doi: 10.1017/S0025315400015009

CrossRef Full Text | Google Scholar

Crothers J. H. (1973). On Variation in Nucella lapillus (L.): Shell Shape in Populations From Pembrokeshire, South Wales. J. Moll. Stud. 40, 319–327. doi: 10.1093/oxfordjournals.mollus.a065229

CrossRef Full Text | Google Scholar

Crothers J. H. (1979). Variation in the Shell of the Dog-Whelk, Nucella lapillus (L.) From Sullom Voe and Other Parts of the Shetland Islands. Mar. Environ. Res. 2, 311–327. doi: 10.1016/0141-1136(79)90015-1

CrossRef Full Text | Google Scholar

Crothers J. H. (1980). Further Observations on the Growth of the Common Dog-Whelk, Nucella lapillus (L.), in the Laboratory. J. Moll. Stud. 46, 181–185. doi: 10.1093/oxfordjournals.mollus.a065531

CrossRef Full Text | Google Scholar

Crothers J. H. (1981). Shell-Shape Variation in Faroese Dog-Whelks (Nucella lapillus (L.)). Biol. J. Linn. Soc Lond. 15, 327–337. doi: 10.1111/j.1095-8312.1981.tb00767.x

CrossRef Full Text | Google Scholar

Crothers J. H. (1982). Shell Shape Variation in Dog-Whelks (Nucella lapillus (L.)) From the West Coast of Scotland. Biol. J. Linn. Soc Lond. 17, 319–342. doi: 10.1111/j.1095-8312.1982.tb02025.x

CrossRef Full Text | Google Scholar

Crothers J. H. (1983). Variation in Dog-Whelk Shells in Relation to Wave Action and Crab Predation. Biol. J. Linn. Soc Lond. 20, 85–102. doi: 10.1111/j.1095-8312.1983.tb01591.x

CrossRef Full Text | Google Scholar

Crothers J. H. (1985a). Two Different Patterns of Shell-Shape Variation in the Dog-Whelk Nucella lapillus (L.). Biol. J. Linn. Soc Lond. 25 (4), 339–353. doi: 10.1111/j.1095-8312.1985.tb00400.x

CrossRef Full Text | Google Scholar

Crothers J. H. (1985b). Dog-Whelks: An Introduction to the Biology of Nucella lapillus (L.). Field Stud. 6, 291–360.

Google Scholar

Currey J. D., Hughes R. N. (1982). Strength of the Dogwhelk Nucella lapillus and the Winkle Littorina Littorea From Different Habitats. J. Anim. Ecol. 51, 47–56. doi: 10.2307/4309

CrossRef Full Text | Google Scholar

Cyronak T., Schulz K. G., Jokiel P. L. (2015). The Omega Myth: What Really Drives Lower Calcification Rates in an Acidifying Ocean. ICES. J. Mar. Sci. 73, 558–562. doi: 10.1093/icesjms/fsv075

CrossRef Full Text | Google Scholar

Doney S. C., Fabry V. J., Feely R. A., Kleypas J. A. (2009). Ocean Acidification: The Other. CO2. Prob. Ann. Rev. Mar. Sci. 1, 169–192. doi: 10.1146/annurev.marine.010908.163834

CrossRef Full Text | Google Scholar

Fisher J. A. D., Rhile E. C., Liu H., Petraitis P. S. (2009). An Intertidal Snail Shows a Dramatic Size Increase Over the Past Century. Proc. Natl. Acad. Sci. U. S. A. 106, 5209–5212. doi: 10.1073/pnas.0812137106

PubMed Abstract | CrossRef Full Text | Google Scholar

Fitzer S. C., Vittert L., Bowman A., Kamenos N. A., Phoenix V. R., Cusack M.. (2015). Ocean Acidification and Temperature Increase Impact Mussel Shell Shape and Thickness: Problematic for Protection? Ecol. Evol. 5, 4875–4884. doi: 10.1002/ece3.1756

PubMed Abstract | CrossRef Full Text | Google Scholar

Fretter V., Graham A. (1962). British Prosobranch Molluscs: Their Functional Anatomy and Ecology (London: Ray Society).

Google Scholar

Galante-Oliveira S., Oliveira I., Santos J. A., Pereira M. de L., Pacheco M., Barroso C. M., et al. (2010). Factors Affecting RPSI in Imposex Monitoring Studies Using Nucella lapillus (L.) as Bioindicator. J. Environ. Monit. 12, 1055–1063. doi: 10.1039/b921834c

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaylord B., Hill T. M., Sanford E., Lenz E. A., Jacobs L. A., Sato K. N., et al. (2011). Functional Impacts of Ocean Acidification in an Ecologically Critical Foundation Species. J. Exp. Biol. 214, 2586–2594. doi: 10.1242/jeb.055939

PubMed Abstract | CrossRef Full Text | Google Scholar

Gazeau F., Gattuso J.-P., Greaves M., Elderfield H., Peene J., Heip C. H. R., et al. (2011). Effect of Carbonate Chemistry Alteration on the Early Embryonic Development of the Pacific Oyster (Crassostrea gigas). PLoS One 6, e23010. doi: 10.1371/journal.pone.0023010

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghalambor C. K., McKAY J. K., Carroll S. P., Reznick D. N. (2007). Adaptive Versus Non-Adaptive Phenotypic Plasticity and the Potential for Contemporary Adaptation in New Environments. Funct. Ecol. 21, 394–407. doi: 10.1111/j.1365-2435.2007.01283.x

CrossRef Full Text | Google Scholar

Giardina C. R., Kuhl F. P. (1977). Accuracy of Curve Approximation by Harmonically Related Vectors With Elliptical Loci. Comput. Graph. Imag. Proc. 6, 277–285. doi: 10.1016/S0146-664X(77)80029-4

CrossRef Full Text | Google Scholar

Guerra-Varela J., Colson I., Backeljau T., Breugelmans K., Hughes R. N., Rolán-Alvarez E. (2009). The Evolutionary Mechanism Maintaining Shell Shape and Molecular Differentiation Between Two Ecotypes of the Dogwhelk Nucella lapillus. Evol. Ecol. 23, 261–280. doi: 10.1007/s10682-007-9221-5

CrossRef Full Text | Google Scholar

Hellberg M. E., Balch D. P., Roy K. (2001). Climate-Driven Range Expansion and Morphological Evolution in a Marine Gastropod. Science 292, 1707–1710. doi: 10.1126/science.1060102

PubMed Abstract | CrossRef Full Text | Google Scholar

Herfort L., Thake B., Roberts J. (2002). Acquisition and Use of Bicarbonate by Emiliania huxleyi. N. Phytol. 156, 427–436. doi: 10.1046/j.1469-8137.2002.00523.x

CrossRef Full Text | Google Scholar

Hofmann G. E., Barry J. P., Edmunds P. J., Gates R. D., Hutchins D. A., Klinger T., et al. (2010). The Effect of Ocean Acidification on Calcifying Organisms in Marine Ecosystems: An Organism-To-Ecosystem Perspective. Annu. Rev. Ecol. Evol. Syst. 41, 127–147. doi: 10.1146/annurev.ecolsys.110308.120227

CrossRef Full Text | Google Scholar

Hughes R. N., Burrows M. T. (1994). An Interdisciplinary Approach to the Study of Foraging Behaviour in the Predatory Gastropod, Nucella lapillus (L.). Nucella. Lapill. (L.). Ethol. Ecol. Evol. 6, 75–85. doi: 10.1080/08927014.1994.9523010

CrossRef Full Text | Google Scholar

Hughes R. N., Elner R. W. (1979). Tactics of a Predator, Carcinus maenas, and Morphological Responses of the Prey, Nucella lapillus. J. Anim. Ecol. 48, 65–78. doi: 10.2307/4100

CrossRef Full Text | Google Scholar

Jokiel P. L. (2011a). The Reef Coral Two Compartment Proton Flux Model: A New Approach Relating Tissue-Level Physiological Processes to Gross Corallum Morphology. J. Exp. Mar. Bio. Ecol. 409, 1–12. doi: 10.1016/j.jembe.2011.10.008

CrossRef Full Text | Google Scholar

IPCC (2014). Intergovernmental Panel on Climate Change 2013 – the Physical Science Basis: Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press).

Google Scholar

Jokiel P. L. (2011b). Ocean Acidification and Control of Reef Coral Calcification by Boundary Layer Limitation of Proton Flux. Bull. Mar. Sci. 87, 639–657. doi: 10.5343/bms.2010.1107

CrossRef Full Text | Google Scholar

Jokiel P. L. (2013). Coral Reef Calcification: Carbonate, Bicarbonate and Proton Flux Under Conditions of Increasing Ocean Acidification. Proc. Biol. Sci. 280, 20130031. doi: 10.1098/rspb.2013.0031

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirby R. R., Bayne B. L., Berry R. J. (1994). Phenotypic Variation Along a Cline in Allozyme and Karyotype Frequencies, and Its Relationship With Habitat, in the Dog-Whelk Nucella lapillus, L. Biol. J. Linn. Soc Lond. 53, 255–275. doi: 10.1111/j.1095-8312.1994.tb01012.x

CrossRef Full Text | Google Scholar

Kitching J. A., Ebling F. J. (1967). Ecological Studies at Lough Ine Vol. 4. Ed. Cragg J. B. (London: Academic Press), 197–291.

Google Scholar

Kitching J. A., Muntz L., Ebling F. J. (1966). The Ecology of Lough Ine. XV. The Ecological Significance of Shell and Body Forms in Nucella. J. Anim. Ecol. 35, 113–126. doi: 10.2307/2693

CrossRef Full Text | Google Scholar

Krause-Jensen D., Duarte C. M., Hendriks I. E., Meire L., Blicher M. E., Marbà N., et al. (2015). Macroalgae Contribute to Nested Mosaics of pH Variability in a Subarctic Fjord. Biogeosciences 12, 4895–4911. doi: 10.5194/bg-12-4895-2015

CrossRef Full Text | Google Scholar

Kroeker K. J., Kordas R. L., Crim R., Hendriks I. E., Ramajo L., Singh G. S., et al. (2013). Impacts of Ocean Acidification on Marine Organisms: Quantifying Sensitivities and Interaction With Warming. Glob. Change Biol. 19, 1884–1896. doi: 10.1111/gcb.12179

CrossRef Full Text | Google Scholar

Kroeker K. J., Sanford E., Rose J. M., Blanchette C. A., Chan F., Chavez F. P., et al. (2016). Interacting Environmental Mosaics Drive Geographic Variation in Mussel Performance and Predation Vulnerability. Ecol. Lett. 19, 771–779. doi: 10.1111/ele.12613

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuhl F. P., Giardina C. R. (1982). Elliptic Fourier Features of a Closed Contour. Comput. Graph. Imag. Proc. 18, 236–258. doi: 10.1016/0146-664X(82)90034-X

CrossRef Full Text | Google Scholar

Lavigne H., Proye A., Gattuso J. P. (2008) Seacarb 2.0, an R Package to Calculate Parameters of the Seawater Carbonate System. Available at: http://cran.r-project.org/web/packages/seacarb/index.html.

Google Scholar

Mackinder L., Wheeler G., Schroeder D., Riebesell U., Brownlee C. (2010). Molecular Mechanisms Underlying Calcification in Coccolithophores. Geomicrobiol. J. 27, 585–595. doi: 10.1080/01490451003703014

CrossRef Full Text | Google Scholar

Marshall B. J., Thunell R. C., Spero H. J., Henehan M. J., Lorenzoni L., Astor Y., et al. (2015). Morphometric and Stable Isotopic Differentiation of Orbulina Universa Morphotypes From the Cariaco Basin, Venezuela. Mar. Micropaleontol. 120, 46–64. doi: 10.1016/j.marmicro.2015.08.001

CrossRef Full Text | Google Scholar

Mayk D. (2021). Transitional Spherulitic Layer in the Muricid Nucella lapillus. J. Moll. Stud. 87 (1), eyaa035. doi: 10.1093/mollus/eyaa035

CrossRef Full Text | Google Scholar

Melatunan S., Calosi P., Rundle S. D., Widdicombe S., Moody A. J. (2013). Effects of Ocean Acidification and Elevated Temperature on Shell Plasticity and Its Energetic Basis in an Intertidal Gastropod. Mar. Ecol. Prog. Ser. 472, 155–168. doi: 10.3354/meps10046

CrossRef Full Text | Google Scholar

Melzner F., Stange P., Trübenbach K., Thomsen J., Casties I., Panknin U., et al. (2011). Food Supply and Seawater pCO₂ Impact Calcification and Internal Shell Dissolution in the Blue Mussel Mytilus edulis. PLoS One 6, e24223. doi: 10.1371/journal.pone.0024223

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore H. B. (1938). The Biology of. Purpura. Lapill. Part II. Growth. J. Mar. Biol. Assoc. U. K. 23, 57. doi: 10.1017/S0025315400053947

CrossRef Full Text | Google Scholar

Morley S. A., Peck L. S., Sunday J. M., Heiser S., Bates A. E. (2019). Physiological Acclimation and Persistence of Ectothermic Species Under Extreme Heat Events. Glob. Ecol. Biogeogr. 28, 1018–1037. doi: 10.1111/geb.12911

CrossRef Full Text | Google Scholar

Olsen A., Lange N., Key R. M., Tanhua T., Bittig H. C., Kozyr A., et al. (2020). An Updated Version of the Global Interior Ocean Biogeochemical Data Product, Glodapv2.2020. Earth Syst. Sci. Data 12, 3653–3678. doi: 10.5194/essd-12-3653-2020

CrossRef Full Text | Google Scholar

Palmer A. R. (1990). Effect of Crab Effluent and Scent of Damaged Conspecifics on Feeding, Growth, and Shell Morphology of the Atlantic Dogwhelk Nucella lapillus (L.). Hydrobiologia 193, 155–182. doi: 10.1007/BF00028074

CrossRef Full Text | Google Scholar

Parker L. M., Ross P. M., O'Connor W. A., Pörtner H. O., Scanes E., Wright J. M.. (2013). Predicting the Response of Molluscs to the Impact of Ocean Acidification. Biology 2, 651–692. doi: 10.3390/biology2020651

PubMed Abstract | CrossRef Full Text | Google Scholar

Pascoal S., Carvalho G., Creer S., Mendo S., Hughes R. N. (2012). Plastic and Heritable Variation in Shell Thickness of the Intertidal Gastropod Nucella lapillus Associated With Risks of Crab Predation and Wave Action, and Sexual Maturation. PLoS One 7, e52134. doi: 10.1371/journal.pone.0052134

PubMed Abstract | CrossRef Full Text | Google Scholar

Pascoe P. L. (2002). Chromosomal Variation in Nucella lapillus (L.) and Other Muricid Gastropods (Plymouth: University of Plymouth).

Google Scholar

Pascoe P. L. (2006). Chromosomal Polymorphism in the Atlantic Dog-Whelk, Nucella lapillus (Gastropoda: Muricidae): Nomenclature, Variation and Biogeography. Biol. J. Linn. Soc Lond. 87, 195–210. doi: 10.1111/j.1095-8312.2006.00567.x

CrossRef Full Text | Google Scholar

Pascoe P. L., Dixon D. R. (1994). Structural Chromosomal Polymorphism in the Dog-Whelk Nucella lapillus (Mollusca: Neogastropoda). Mar. Biol. 118, 247–253. doi: 10.1007/BF00349791

CrossRef Full Text | Google Scholar

Przeslawski R., Byrne M., Mellin C. (2015). A Review and Meta-Analysis of the Effects of Multiple Abiotic Stressors on Marine Embryos and Larva E. Glob. Change Biol. 21, 2122–2140. doi: 10.1111/gcb.12833

CrossRef Full Text | Google Scholar

R Core Team (2020). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. Available at: https://www.R-project.org/.

Google Scholar

Rühl S., Calosi P., Faulwetter S., Keklikoglou K., Widdicombe S., Queirós A. M. (2017). Long-Term Exposure to Elevated pCO₂ More Than Warming Modifies Early-Life Shell Growth in a Temperate Gastropod. ICES. J. Mar. Sci. 74, 1113–1124. doi: 10.1093/icesjms/fsw242

CrossRef Full Text | Google Scholar

Sanders T., Thomsen J., Müller J. D., Rehder G., Melzner F. (2021). Decoupling Salinity and Carbonate Chemistry: Low Calcium Ion Concentration Rather Than Salinity Limits Calcification in Baltic Sea Mussels. Biogeosciences 18, 2573–2590. doi: 10.5194/bg-18-2573-2021

CrossRef Full Text | Google Scholar

Staiger H. (1953). Der Chromosomendimorphismus Beim Prosobranchier Purpura lapillus in Beziehung Zur Ökologie Der Art. Chromosoma 6, 419–477. doi: 10.1007/BF01259949

CrossRef Full Text | Google Scholar

Stumpp M., Hu M. Y., Melzner F., Gutowska M. A., Dorey N., Himmerkus N., et al. (2012). Acidified Seawater Impacts Sea Urchin Larvae pH Regulatory Systems Relevant for Calcification. Proc. Natl. Acad. Sci. U. S. A. 109, 18192–18197. doi: 10.1073/pnas.1209174109

PubMed Abstract | CrossRef Full Text | Google Scholar

Suckling C. C., Clark M. S., Richard J., Morley S. A., Thorne M. A. S., Harper E. M., et al. (2015). Adult Acclimation to Combined Temperature and pH Stressors Significantly Enhances Reproductive Outcomes Compared to Short-Term Exposures. J. Anim. Ecol. 84, 773–784. doi: 10.1111/1365-2656.12316

PubMed Abstract | CrossRef Full Text | Google Scholar

Sultan S. E., Stearns S. C. (2005). “Chapter 14 - Environmentally Contingent Variation: Phenotypic Plasticity and Norms of Reaction,” in Variation. Eds. Hallgrıímsson B., Hall B. K. (London: Academic Press), 303–332. doi: 10.1016/B978-012088777-4/50016-8

CrossRef Full Text | Google Scholar

Telesca L., Michalek K., Sanders T., Peck L. S., Thyrring J., Harper E. M., et al. (2018). Blue Mussel Shell Shape Plasticity and Natural Environments: A Quantitative Approach. Sci. Rep. 8, 2865. doi: 10.1038/s41598-018-20122-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Telesca L., Peck L. S., Sanders T., Thyrring J., Sejr M. K., Harper E. M., et al. (2019). Biomineralization Plasticity and Environmental Heterogeneity Predict Geographical Resilience Patterns of Foundation Species to Future Change. Glob. Change Biol. 25, 4179–4193. doi: 10.1111/gcb.14758

CrossRef Full Text | Google Scholar

Telesca L., Peck L. S., Backeljau T., Heinig M. F., Harper E. M. (2021). A Century of Coping With Environmental and Ecological Changes via Compensatory Biomineralization in Mussels. Glob. Change Biol. 27, 624–639. doi: 10.1111/gcb.15417

CrossRef Full Text | Google Scholar

Thompson D. W. (1992). On Growth and Form (Cambridge: Cambridge University Press). doi: 10.1017/CBO978110732585210.1017/CBO9781107325852

CrossRef Full Text | Google Scholar

Thomsen J., Haynert K., Wegner K. M., Melzner F. (2015). Impact of Seawater Carbonate Chemistry on the Calcification of Marine Bivalves. Biogeosciences 12, 4209–4220. doi: 10.5194/bg-12-4209-2015

CrossRef Full Text | Google Scholar

Thyng K., Greene C., Hetland R., Zimmerle H., DiMarco S. (2016). True Colors of Oceanography: Guidelines for Effective and Accurate Colormap Selection. Oceanography 29, 9–13. doi: 10.5670/oceanog.2016.66

CrossRef Full Text | Google Scholar

Tresguerres M. (2016). Novel and Potential Physiological Roles of Vacuolar-Type H⁺-ATPase in Marine Organisms. J. Exp. Biol. 219, 2088–2097. doi: 10.1242/jeb.128389

PubMed Abstract | CrossRef Full Text | Google Scholar

Trussell G. C., Ewanchuk P. J., Bertness M. D. (2003). Trait-Mediated Effects in Rocky Intertidal Food Chains: Predator Risk Cues Alter Prey Feeding Rates. Ecology 84, 629–640. doi: 10.1890/0012-9658(2003)084[0629:TMEIRI]2.0.CO;2

CrossRef Full Text | Google Scholar

Vermeij G. J., Currey J. D. (1980). Geographical Variation in the Shell Strength of Thaidid Snail Shells. Biol. Bull. 158, 383–389. doi: 10.2307/1540864

CrossRef Full Text | Google Scholar

Viotti S., Sangil C., Hernández C. A., Hernández J. C. (2019). Effects of Long-Term Exposure to Reduced pH Conditions on the Shell and Survival of an Intertidal Gastropod. Mar. Environ. Res. 152, 104789. doi: 10.1016/j.marenvres.2019.104789

PubMed Abstract | CrossRef Full Text | Google Scholar

Waldbusser G. G., Brunner E. L., Haley B. A., Hales B., Langdon C. J., Prahl F. G., et al. (2013). A Developmental and Energetic Basis Linking Larval Oyster Shell Formation to Acidification Sensitivity. Geophys. Res. Lett. 40, 2171–2176. doi: 10.1002/grl.50449

CrossRef Full Text | Google Scholar

Waldbusser G. G., Hales B., Langdon C. J., Haley B. A., Schrader P., Brunner E. L., et al. (2014). Saturation-State Sensitivity of Marine Bivalve Larvae to Ocean Acidification. Nat. Clim. Change 5, 273–280. doi: 10.1038/nclimate2479

CrossRef Full Text | Google Scholar

Waldbusser G. G., Hales B., Langdon C. J., Haley B. A., Schrader P., Brunner E. L., et al. (2015). Ocean Acidification has Multiple Modes of Action on Bivalve Larvae. PLoS One 10, e0128376. doi: 10.1371/journal.pone.0128376

PubMed Abstract | CrossRef Full Text | Google Scholar

Waldbusser G. G., Hales B., Haley B. A. (2016). Calcium Carbonate Saturation State: On Myths and This or That Stories. ICES. J. Mar. Sci. 73, 563–568. doi: 10.1093/icesjms/fsv174

CrossRef Full Text | Google Scholar

Waters J. F., Millero F. J. (2013). The Free Proton Concentration Scale for Seawater pH. Mar. Chem. 149, 8–22. doi: 10.1016/j.marchem.2012.11.003

CrossRef Full Text | Google Scholar

Watson S.-A., Peck L. S., Tyler P. A., Southgate P. C., Tan K. S., Day R. W., et al. (2012). Marine Invertebrate Skeleton Size Varies With Latitude, Temperature and Carbonate Saturation: Implications for Global Change and Ocean Acidification. Glob. Change Biol. 18, 3026–3038. doi: 10.1111/j.1365-2486.2012.02755.x

CrossRef Full Text | Google Scholar

West-Eberhard M. J. (1989). Phenotypic Plasticity and the Origins of Diversity. Annu. Rev. Ecol. Syst. 20, 249–278. doi: 10.1146/annurev.es.20.110189.001341

CrossRef Full Text | Google Scholar

Wood S. N. (2011). Fast Stable Restricted Maximum Likelihood and Marginal Likelihood Estimation of Semiparametric Generalized Linear Models. J. R. Stat. Soc Ser. B. Stat. Methodol. 73, 3–36. doi: 10.1111/j.1467-9868.2010.00749.x

CrossRef Full Text | Google Scholar

Zeebe R. E., Wolf-Gladrow D. (2001). CO2 in Seawater: Equilibrium, Kinetics, Isotopes (Houston, United States: Gulf Professional Publishing).

Google Scholar

Zuur A. F., Ieno E. N., Elphick C. S. (2010). A Protocol for Data Exploration to Avoid Common Statistical Problems. Meth. Ecol. Evol. 1, 3–14. doi: 10.1111/j.2041-210X.2009.00001.x

CrossRef Full Text | Google Scholar