Exoskeletal predator defenses of juvenile California spiny lobsters (Panulirus interruptus) are affected by fluctuating ocean acidification-like conditions

Introduction

Marine invertebrates possess a diverse array of weapons and defense structures that use biomineralization to confer sufficient hardness and stiffness for optimal function. Structures spanning spines, plates, claws, mandibles, antennae, raptorial appendages, shells, and exoskeletons all have hierarchical structural design elements fortified with minerals, such as calcium carbonates and calcium phosphates, that help produce the material properties for a specific function. Yet, biomineralization, particularly calcification, in marine organisms is demonstrably sensitive to ocean acidification-like reduced pH/high pCO₂ conditions (Kroeker et al., 2010). Most invertebrates with calcified structures experience reductions in mineral content (Chan and Connolly, 2013; Kroeker et al., 2013; Davis et al., 2021) under reduced pH/high pCO₂ exposure due to lowered carbonate ion concentration and carbonate mineral saturation state. In contrast, some crustacean species experience increases in calcification (McDonald et al., 2009; Long et al., 2013; Taylor et al., 2015; deVries et al., 2016; Glandon et al., 2018), which is potentially accomplished by drawing upon bicarbonate accumulated for acid-base regulation (Truchot, 1979; Pane and Barry, 2007; Spicer et al., 2007; Knapp et al., 2015; Whiteley et al., 2018). The degree of mineralization is strongly linked with material properties, whereby greater mineral content corresponds to greater stiffness and brittleness (Currey, 1984). Thus, disruptions to the biomineralization process under ocean acidification (OA) conditions, resulting in either decreases or increases to calcification, have the potential to alter the mechanical integrity of defense structures and ultimately the potential outcomes of predator-prey interactions in future ocean conditions.

OA research often targets biomineralization responses, but rarely are the effects on material properties and mechanical functioning, which are consequential for survival during predator-prey interactions, studied concurrently. Reduced calcification is generally purported to have negative consequences on material properties and mechanical efficiency. This is exemplified by reduced calcification corresponding with lower stem stiffness in a calcified green alga and decreased shell elasticity in a Mediterranean top snail, both studied at CO₂ seep sites (Newcomb et al., 2015; Duquette et al., 2017), and decreased abalone shell fracture toughness concurrent with reductions in shell mineralization and thickness under reduced pH conditions (Auzoux-Bordenave et al., 2020; Avignon et al., 2020). Increased calcification, however, can also create weaker barnacle shell wall plates (McDonald et al., 2009) due to increased brittleness. In mantis shrimp, exoskeletal hardness and stiffness remained the same despite increases to Mg content (deVries et al., 2016). Altered mineralization, either through increases or decreases, has the potential to disrupt the material properties and mechanical efficiency of organism structures, but the thresholds for these effects are likely specific to structures and species (Siegel et al., 2022). As such, there is still much to learn about the relationship between changes to material composition and the degree to which the material is compromised. Unraveling this relationship will help inform how single-species responses can be inferred to impact species interactions, such as predator-prey encounters.

A better understanding of potential impacts is especially needed for crustaceans, which have generally garnered less OA research attention than other marine calcifiers (Whiteley, 2011; Siegel et al., 2022), but have an exoskeleton that is complex, multifunctional, and highly variable within a single individual. The decapod crustacean exoskeleton is heterogenous and comprised of mineral, usually in the form of calcite, magnesian calcite, amorphous calcium carbonate, and sometimes calcium phosphate, that is embedded among lamellae of sugars and protein, which stack to form the distinct layers of the epicuticle, exocuticle, and endocuticle (Travis, 1963; Huner et al., 1979; Vigh and Dendinger, 1982; Addadi et al., 2003; Dillaman et al., 2005; Boßelmann et al., 2007; Kunkel et al., 2012). Within the exocuticle and endocuticle, layer thickness, lamellae stacking density, and mineral composition can all be adjusted to produce a spectrum of material properties that support a range of exoskeleton functions (Raabe et al., 2005; Chen et al., 2008). Crustaceans can construct exoskeletal regions with specificity to meet functional demands (Melnick et al., 1996; Dutil et al., 2000; Chen et al., 2008) and likely have complex responses to OA-like conditions (Siegel et al., 2022).

Exoskeleton versatility is exemplified in the California spiny lobster, Panulirus interruptus, which begin living in coastal kelp forests as juveniles and use their exoskeleton in a suite of defenses against predators such as octopus, eels, and large fish, particularly the California sheephead, Semicossyphus pulcher (Lindberg, 1955; Loflen and Hovel, 2010; Higgins et al., 2018). One of the first strategies that spiny lobsters can employ to avoid harm is using the rostral horns to lock onto the roof of a rocky shelter to avoid expulsion (Lindberg, 1955). Another option is to use their long, spiny antennae to push against predators, keeping the bulk of their body a safe distance away (Barshaw et al., 2003; Briones-Fourzán et al., 2006). As a last line of defense, the carapace is coated in flattened spines that likely enhance pierce-resistance (Tarsitano et al., 2006). Each of these defense behaviors are dependent on the exoskeleton being sufficiently stiff, hard, and strong under their different loading regimes, and their distinct composition and structure affords them the proper mechanical properties to achieve their function (Lowder, 2019). These defense structures, however, may be differentially vulnerable to the impacts of changing carbonate chemistry and therefore yield greater insight into the variable effects of OA on the structure-property relationships of calcified materials.

Understanding the impacts on spiny lobsters is not only complicated by the variability in their exoskeleton, but also by the variability in their environment. Spiny lobsters, like a number of other coastal organisms, live in dynamic environments that may shape their response to changing ocean conditions. Juvenile P. interruptus spend two to three years in shallow eelgrass or kelp environments (Engle, 1979), where they must be able to successfully exercise their full range of predator defense strategies while contending with diurnal changes in their environment (Frieder et al., 2012). In coastal environments, carbonate chemistry fluctuates diurnally, tidally, and seasonally (Johnson et al., 2013), producing variable conditions to which organisms in these habitats are subjected. In environments dominated by a large biomass of photosynthetic organisms, such as kelp forests or seagrass meadows, daily ranges of pH may be as great as the decrease in pH expected over the next century, 0.3-0.4 units (Frieder et al., 2012; Challener et al., 2015; James et al., 2019)driven in part by the uptake of CO₂via photosynthesis during the day and organismal respiration at night. These ranges can be even greater than 1 unit (Semesi et al., 2009) depending on currents, tides, depth, and community composition (Hofmann et al., 2011; Page et al., 2016; Silbiger and Sorte, 2018). In addition, in coastal regions where seasonal upwelling is prominent, such as the U.S. West Coast, further decreases in near-shore pH may occur owing to intrusion of deep seawater with high CO₂ and low pH (Takeshita et al., 2015; Chan et al., 2017; Kekuewa, 2020).

For organisms that inhabit dynamic environments, the exposure to a range of conditions has potentially led to a degree of adaptation that may be beneficial as ocean carbonate chemistry continues to change. A growing, but still limited, number of studies have imposed fluctuations around low and high pH setpoints and discovered varied responses. For example, fluctuating pH conditions may ameliorate some or all of the negative growth and sensory impacts observed at stable reduced pH in pink salmon (Ou et al., 2015) and coral reef fish, where it is hypothesized to help prevent physiological changes in the brain (Jarrold et al., 2017). Mytilus edulis mussels take advantage of higher pH during the day to increase calcification rates, lessening the overall negative impact to calcification (Wahl et al., 2018), while Mytilus galloprovincialis veligers experiencing fluctuating reduced pH did not experience the same reductions in shell length that those in stable reduced pH did (Frieder et al., 2014). Yet, in some organisms, fluctuating pH conditions either have no discernable effect from stable reduced pH conditions (Clark and Gobler, 2016), or they appear to increase stress, producing more negative responses than stable reduced pH conditions alone (Alenius and Munguia, 2012; Mangan et al., 2017; Onitsuka et al., 2017; Johnson et al., 2019). Thus, to establish realistic impacts on organisms like the California spiny lobster, it is important to incorporate the natural fluctuating conditions of their habitat into OA experiments (Andersson and Mackenzie, 2012; McElhany and Busch, 2013; Andersson et al., 2015).

In this study, we exposed juvenile California spiny lobsters, Panulirus interruptus, to stable and fluctuating high pCO₂/reduced pH conditions that mimic their current and future habitat and examined the integrity of their predator defenses. Specifically, we examined the rostral horns, antennae, and carapace for elemental composition (Ca and Mg content), structure (cuticle and layer thickness), and material properties (hardness and stiffness) (Figure 1). We hypothesized that lobsters under low, stable pH conditions would have relatively higher exoskeleton mineralization than those in ambient conditions, resulting in altered mechanical properties, while these effects would be lessened in the fluctuating reduced pH conditions due to a respite from constant reduced pH conditions. Alternatively, the fluctuating reduced pH conditions may exceed potential pH thresholds, eliciting the same or exacerbated responses. We also hypothesized that the magnitude of response would vary across defensive structures due to localized control and functional specificity.

Figure 1 Exoskeletal defense structures of Panulirus interruptus. (A) antenna, (B) horn tip, (C) carapace spine, (D) abdominal segment in image columns I and II. The analyses performed for each structure are listed in column III. Arrows indicate areas where nanoindentation was performed normal to the surface. The horn tip in B-II is affixed to an aluminum block for nanoindentation.

Animal collection

Sixty-four juvenile lobsters were collected from offshore La Jolla, San Diego, CA over the course of six days in October 2016. Commercial lobster traps were modified with smaller mesh and set at depths of 8.5 to 15.2 m. Lobsters were transported to the experimental aquarium at Scripps Institution of Oceanography (SIO) within 2.5 hours of collection and placed in flow-through seawater pumped from the Scripps Pier, approximately 2 km northeast from the site of collection.

Experimental design

The experimental ocean acidification system consisted of four 81 L header tanks that each received filtered seawater pumped from the SIO pier (3-4 m depth, 300 m offshore) at ambient pH (7.97 ± 0.03), pCO₂ (464 ± 39 μatm), and salinity (33.4 ± 0.2 PSU) (mean ± s.d. during the experimental period). A mix of chilled and ambient temperature seawater was fed into the header tanks, which were maintained at a stable 16.5 ± 0.6°C with the aid of heaters. Two 1585 GPH aquarium powerheads were placed in each header tank to ensure adequate mixing.

One header tank was kept at ambient pH to provide the control treatment while the three others were adjusted for high pCO₂/low pH treatments. One treatment was adjusted to a stable pH of 7.7, congruent with current predictions for decreases in global ocean surface pH values of 0.3 pH units by the year 2100 (IPCC, 2014). pH levels as low as 7.77 have been observed at 10 m in the habitat of spiny lobsters during upwelling seasons (Kekuewa, 2020). In the kelp forest less than 1000 m from the site of lobster collection, pH was found to average 8.05 units at 7 m depth with minima around 7.8 during a springtime period of observation (Frieder et al., 2012). The two other treatments were set to fluctuate around this reduced pH, but to different degrees. The reduced/low fluctuating pH treatment was set to 7.7 ± 0.05 (total range: 7.65-7.75) pH units based on published measurements of southern California kelp forest diurnal pH variations, which depend upon tide, season, and weather conditions, and which are approximately 0.1 pH units on average (Frieder et al., 2012; Kapsenberg and Hofmann, 2016). The reduced/high fluctuating pH treatment was set to 7.7 ± 0.10 (total range: 7.6-7.8) pH units, which mimics future diurnal fluctuations, calculated in CO₂Sys for decreases in buffering capacity associated with higher anthropogenic pCO₂ values and thus larger pH swings assuming the same autotrophic forcing of DIC (Pierrot et al., 2006). For both of the fluctuating treatments, high pH values were held stable for approximately 10 hours to match daylight (and photosynthesis) hours at the time of the experiment. Nighttime low pH values were held stable for 12 hours. pH was set at a midpoint value for one hour between each of these highs and lows as a transition period.

Experimental pH conditions were achieved by bubbling 100% CO₂ into the three treatment header tanks and were controlled by an Apex Lite aquarium controller with Apex Neptune pH and temperature probes (0.01 pH accuracy; temperature accuracy 0.1°C, Neptune Systems, Morgan Hill, CA, USA). pH and temperature of the header tanks were logged every 10 minutes throughout the duration of the experiment. The Apex pH probes were calibrated prior to the experimental period, but because stability and not accuracy was desired, they were not calibrated again.

Each header tank supplied flow-through seawater separately to 16 experimental tanks (8.4 L) that were randomly intermixed in the experimental area. Each experimental tank contained a single lobster, for a combined total of 64 individuals. Experimental water parameters were gradually altered over the course of three days to minimize potential acute stress. Once the targeted water parameters were reached, the experiment was run for 102 days.

Water chemistry

Daily readings of pH and temperature were taken from each lobster tank using a portable probe (HQ40d, probe PHC201, precision 0.01 pH, 0.1 temperature, Hach, Loveland, CO, USA) as a method of high-frequency measurements in between discrete carbonate chemistry water sample measurements (Figure S1). This probe was calibrated with NBS buffer solutions (Fisher Scientific, Fair Lawn, NJ, USA) approximately every three weeks to control its slow drift. Daily readings were taken at consistent times during the daytime to capture the high value in the fluctuating treatments. Calculated pH from discrete water samples (see below) was averaged for each sampling period and used to correct the portable probe values to pH_total. Each water sampling round corresponded to one calibration period of the portable probe.

Nighttime pH was checked using two methods to confirm that stable treatments were constant and that target fluctuations were achieved and consistent. In addition to four probes constantly monitoring and logging header tank pH, two additional Neptune Apex electrodes were rotated between random lobster tanks once per day and recorded pH every 10 minutes over the course of each 24-hour period. Nighttime pH in the fluctuating treatments was recorded to be 0.11 ± 0.02 (n = 54) and 0.21 ± 0.03 (n = 46) pH units lower than daytime pH, consistent with our targeted values. Periodic measurements of nighttime pH with the handheld probe (n=8) were in agreement with the readings of the Apex probes.

Discrete samples were collected for water chemistry analysis four times from three randomized lobster tanks per treatment over the experimental period (48 samples total). Samples were collected in 500 mL Pyrex glass bottles and poisoned with 120 μL of a saturated solution of HgCl₂ in accordance with standard operating procedures (Dickson et al., 2007). Each water sample was analyzed for total alkalinity (A_T) by potentiometric acid titration (0.1 N HCl) using an 876 Dosimat (Metrohm, Riverview, Florida, USA) and pH electrode (Ecotrode Plus, Metrohm, Riverview, Florida, USA) and for dissolved inorganic carbon (DIC) using an AIRICA system (Marianda, Kiel, Germany) with integrated infrared CO₂ analyzer (Li-COR 7000, Li-COR, Lincoln, Nebraska, USA). Certified reference material (CRM) provided by the Dickson laboratory at SIO were used to verify the accuracy and precision, both to ± 1–4 μmol kg⁻¹, before and after measuring samples as well as between every five samples for DIC and ten for A_T. Other parameters of the carbonate system were calculated using CO₂Sys 25.06 (Pierrot et al., 2006) (Table 1). For calculations, dissociation constants of K₁, K₂ were from Mehrbach et al., 1973, refit by Dickson and Millero, 1987. The HSO₄^- constant was from Dickson, 1990 and the [B]T value was from Uppström, 1974.

Table 1 Water parameters during the experimental period as mean ± sd.

Survival, molting, and growth

Lobsters were fed equal amounts of squid five days per week in excess, with remaining food removed daily. Tanks were checked daily for deaths and molting events; exuviae were promptly removed to prevent fouling and consumption. Growth was recorded at the start and end of the experiment by weighing damp lobsters in air after removal from tanks and allowing excess water to run off. Carapace length was measured from the notch at the base of the rostral horns to the posterior edge of the carapace at the dorsal midline. Initial size ranged from 43.5 to 58.5 mm carapace length (CL) and 75.4 to 198.2 g. Of the 64 lobsters, 20 were male and 44 were female, and 30% showed at least some physical indications of sexual maturity (e.g., softened sternal plates and enlarged testes identified during dissection). Carapace length and sex were used to assign individuals to each treatment so that these factors were evenly distributed.

Cuticle structure and elemental composition

At the end of the experiment, lobsters were anesthetized by brief placement in a -20°C freezer and then sacrificed by piercing the exoskeleton and tissue between the rostral horns with a ceramic knife. Samples of cuticle were dissected by cutting ~ 0.75 x 0.75 cm² around a consistent carapace spine located below the cephalothorax line, a 0.5 x 1.5 cm² section of the antennae 1 cm above the posterior base on the dorsal side, and around the base of the rostral horn. Each cuticle sample was rinsed with DI water and allowed to air dry. Samples were then freeze-fractured with liquid nitrogen so that an orthogonal cross-section was produced and critical-point dried (AutoSamdri 815 Series A, Tousimis, Rockville, MD, USA) before being mounted on a 90-degree SEM tip and sputter-coated with iridium.

Cross-sections of these cuticle samples were examined with ultra-high-resolution scanning electron microscopy under high vacuum (XL30 SFEG with Sirion column and Apreo LoVac, FEI, Hillsboro, OR, USA with Oxford X-MAX 80 EDS detector, Concord, MA, USA) at 10 or 20 kV. One to two samples each of the carapace spine and antenna from individual lobsters were imaged and measured for the total cuticle thickness (epicuticle, exocuticle, and endocuticle), as well as thickness of the individual exo- and endocuticle layers, which were differentiated by the abrupt change in lamellae stacking density (Figure S3). The rostral horn has an unusual structure, consisting of only two layers: an inner region, that we will refer to as the core, and an outer layer that we will refer to as the outer ring (Lowder, 2019). Layer thickness (and the total thickness) for which the sample did not produce a clean cross section, most commonly observed as a flaky and uneven endocuticle, were not included in analyses (n=0-2 per layer per treatment).

Elemental composition of cuticle samples was quantified using both energy-dispersive x-ray spectroscopy (EDX) and inductively coupled plasma mass spectrometry (ICP-MS). The spatial resolution of EDX enabled elemental sampling of the surface of exocuticle and endocuticle layers separately, whereas ICP-MS provided precise quantification of elements present in the whole volume of cuticle sample. EDX was measured with two machines (XL30 SFEG with Sirion column and Apreo LoVac, FEI, Hillsboro, OR, USA with Oxford X-MAX 80 EDS detector, Concord, MA, USA) at 20 kV acceleration voltage. Spectra were taken on the cross-sectional surface of the exocuticle and the endocuticle layers of the carapace spine and antenna base and the core and outer ring of the horn tip. Eleven elements were consistently present in varying amounts: C, N, O, Ca, Mg, Na, Al, P, Si, S, and Cl. The amount of Ca and Mg were expressed as weight percent (% wt) relative to the other nine elements.

Inductively-coupled plasma mass spectrometry (ICP-MS) was performed on cuticle samples at the Scripps Isotope Geochemistry Laboratory (SIGL) for a precise quantification of elements. The carapace spine was air-dried and then trimmed so only the spine remained with no setae, and the abdominal segment was cut as a 1 x 1 cm² from the center of the second abdominal segment and air-dried. Samples were weighed and placed in Teflon vials for digestion with 0.5 ml of concentrated Teflon-distilled (TD) nitric acid (HNO₃) on a hotplate at 100°C for >24 h. Samples were dried down and diluted by a factor of 4000 with 2% TD HNO₃ before being transferred to pre-cleaned centrifuge tubes for analysis. Samples were doped with an indium solution to monitor instrumental drift. Measurements were done using a ThermoScientific iCAPq c ICP-MS (Thermo Fisher Scientific GmbH, Bremen, Germany) in standard mode. Masses of Mg and Ca were sequentially measured for 30 ratios, resulting in internal precision of <2% (2 s.d.). Elements were corrected for total mole fraction. Total procedural blanks represented <0.3% of the measurement for Mg and Ca. Raw data were corrected offline for instrument background and drift. Samples were bracketed by internal standards of crab carapace (n=3), which allowed for calculation of absolute values. The standards yielded external precision of 2% and 3% for Mg and Ca, respectively. We targeted the following isotopes: ¹⁰B, ²⁵Mg, ²⁶Mg, ²⁷Al, ²⁹Si, ³¹P, ⁴³Ca, ⁴⁸Ca, ⁵⁴Fe, ⁵⁷Fe, ⁶⁵Cu, ⁶⁶Zn, ⁸⁶Sr, ¹¹¹Cd, ¹³⁷Ba, ²⁰⁸Pb, and ²³⁸U (Figures S4, S5). A weighted average based on the natural abundance of isotopes was calculated for Ca and Mg and used in analyses.

Material properties

The carapace spine and rostral horn tip were tested for hardness and stiffness using a nanoindentation materials testing machine (Nano Hardness Tester, Nanovea, Irvine, CA, USA) equipped with a Berkovich tip. Fresh samples (<12 hours, except for two samples tested within 24 hours) were kept hydrated in seawater until testing. Samples were secured to an aluminum block with cyanoacrylate glue such that the outer surface was facing up. Indentations were performed by applying a load of 40 mN to the outer surface of the sample at loading and unloading rates of 80 mN/min with a 30 sec hold for creep. At least three indents were taken and averaged per sample.

Antenna stiffness

One antenna was removed from each lobster at the base (first segment above the antennae knuckle). The proximal 5 mm of the antennal base was embedded in epoxy (EpoxiCure, Buehler, Lake Bluff, IL, USA) and allowed to cure overnight. Antennae were covered with seawater-soaked paper towels to remain hydrated until testing. The embedded ends of the antennae were clamped in a universal materials testing machine (E1000, 250N load cell, Instron, Norwood, MA, USA) such that the antennae extended laterally. A force was applied at a rate of 20 mm/min to two locations on the antennae, one at 1/3 of the total antenna length (“proximal”) and one at 2/3 of the total antenna length (“distal”), until a displacement of 10% of antenna length was reached. The stiffer proximal region of the antenna has a thicker ovular cross section than the distal location, and it lies directly over the carapace when the lobster raises its antennae in a defensive position, so it is believed to be the primary location where force is applied by predators (Lowder, pers. obs.). The tapering of the antennae imbues more flexibility towards the distal end, perhaps to avoid breakage. Flexural stiffness, EI, was calculated as:

Statistics

All analyses were carried out in R version 3.5.1 (R Core Team, 2019) and associated packages (Wickham, 2016a, 2016b, 2017).

We tested survival, growth, cuticle measurements (elemental composition, thickness, material properties), and antennae flexural stiffness data for normality and homogeneity of variances with the Shapiro-Wilk and Bartlett’s tests (after the removal of outliers, see below). If data met these conditions, we compared metrics across treatments using one-way ANOVAs and followed with Tukey’s Honest Significant Difference test if appropriate. Otherwise, the Kruskal-Wallis with Dunn’s test or Welch’s corrected ANOVA with Games-Howell post-hoc test was used for instances of non-normality and heteroscedasticity, respectively. Carapace length ANOVA residuals were checked for test assumptions visually and with the Shapiro-Wilks test, then log-transformed to meet assumptions. α = 0.05 with the exception of elemental data, where we used a Bonferroni correction for repeated measures.

Some cuticle samples had one or two material properties tests that resulted in unusually high values for crustacean cuticle. In many of these instances, the indent was misshapen (stretched in one direction), indicating that the force was not applied perpendicularly to the surface in the manner that is required for accurate measurements. As such, material property outliers of each treatment were calculated and removed via the Tukey “fence” or “boxplot rules” method of Q₃ + 1.5 x IQR, where Q₃ is the third quartile and IQR is the interquartile range from the first to the third quartile (Tukey, 1977). The lower “fence” was also calculated, but was always below zero. This method identified outliers ≥ 13.4 GPA in stiffness for the horn and ≥13.1 GPa for the carapace spine. Hardness outliers were ≥ 0.49 GPa for the horn and ≥ 0.52 GPa for the carapace spine.

Results

Survival and molting

All lobsters survived the experimental period and appeared healthy at its conclusion. Molting only occurred in nine out of the sixty-four animals (n=4 in ambient/stable; n=1 in reduced/stable pH; n=3 in reduced/low fluctuating; n=1 in reduced/high fluctuating treatment) and were not included in subsequent analyses. Growth results are included in the supplementary material. The remaining lobsters were molt-staged based on exoskeleton hardness and evidence of apolysis (Travis, 1954; Lyle and MacDonald, 1983) and were determined to be in intermolt, excluding two lobsters that showed signs of apolysis (one from the reduced/stable pH treatment and one from the reduced/high fluctuating pH treatment).

Cuticle structure

There were no significant differences in total cuticle thickness of the carapace (ANOVA: F_3,47= 1.283, p=0.291) and antennae (ANOVA: F_3,51= 2.006, p=0.125) among treatments (Figures 2A, B;Lowder et al., 2022; Table 2). Additionally, there were no significant differences in thickness of either the exocuticle or endocuticle layers of either structure (Figures 2C, D; ANOVA, carapace exocuticle: F_3,50= 0.522, p=0.669; endocuticle: F_3,47= 1.874, p=0.147; antenna exocuticle: ANOVA: F_3,50= 1.675, p=0.184; endocuticle: ANOVA: F_3,50= 0.562, p=0.642).

Figure 2 Total cuticle thickness of the carapace (A) and antennae (B) and cuticle layer thickness of the carapace (C) and antennae (D). Exocuticle layers are right boxes of each treatment (stripes in C, D) and endocuticle layers are left boxes of each treatment (solids in C, D). There were no differences in total thickness or layer thickness of the carapace or antennae among treatments. Box boundaries are the first and third quartiles, whiskers are 1.5 times the interquartile range, the center line is the median, black and white diamonds are the mean, and outliers are solid dots. Values on the bottom of panels denote sample size.

Table 2 Mineralization (µmol/mg from ICP-MS measurements; % wt from EDX measurements), ultrastructure (µm from SEM measurements), and mechanical properties (GPa from nanoindentation tests) for the four studied structures across treatments.

Cuticle elemental composition

Elemental data from the XL30 SFEG EDX were inconsistent with that of the other EDX instrument and with previous data on lobster cuticle (Lowder, 2019), and were thus excluded from analysis. Analyses carried out with the full set of data did not result in different conclusions. The remaining data (n=5-9 samples per treatment) revealed no significant differences in % wt Ca in either of the carapace layers (Figure 3A; exocuticle: Kruskal-Wallis, χ² = 2.154, df=3, p = 0.541; endocuticle: ANOVA, F_3,25 = 1.081, p = 0.375), horn layers (Figure 3C; outer ring: Kruskal-Wallis, χ² = 0.453, df=3, p = 0.929; core: ANOVA, F_3,21 = 0.04, p = 0.989), and antenna layers (Figure 3B; exocuticle: ANOVA, F_3,25 = 1.873, p = 0.160; endocuticle: ANOVA, F_3,25 = 0.658, p = 0.585; Table 2).

Figure 3 Calcium and (A-C) Mg (D-F) % wt of three exoskeletal structures from EDX measurements. There were no significant differences in either element across treatment conditions. Each structure’s exocuticle, or the outer layer in the case of the horn, is represented as the right boxes of each treatment (striped shades), and the endocuticle layers or core of the horn are the left boxes of each treatment (solid shades). Box boundaries are the first and third quartiles, whiskers are 1.5 times the interquartile range, the center line is the median, and black and white diamonds are the mean. Values on the bottom of panels (D-F) denote sample size for each structure.

There were also no significant differences in % wt Mg in either of the carapace layers (Figure 3D; exocuticle: Welch’s ANOVA, F_3,11.80 = 2.389, p = 0.121; endocuticle: Kruskal-Wallis, χ² = 1.266, df=3, p = 0.737), horn layers (Figure 3F); outer ring: Kruskal-Wallis, χ² = 0.754, df=3, p = 0.861; core: Kruskal-Wallis, χ² =  2.559, df=3, p = 0.465), and antenna layers (Figure 3E; exocuticle: Kruskal-Wallis, χ² = 6.061, df=3, p = 0.109; endocuticle: Kruskal-Wallis, χ² = 0.968, df=3, p = 0.809; Table 2).

Elemental analysis of the carapace using ICP-MS contrasted with the results obtained from the EDX analysis; the former measured elements over whole volume cuticle whereas the latter measured elements over a surface area whose features can influence electron scattering and thus element quantification. There were significant differences in calcium and magnesium content in the carapace across treatments (calcium: ANOVA, F_3,51 = 4.313, p = 0.009; magnesium: Kruskal-Wallis, χ² = 15.184, df=3, p = 0.002; Figures 4A, C; Table 2). Calcium was 6% lower in the reduced/high fluctuating treatment than the reduced/stable treatment (Tukey HSD, p = 0.043) and 7% lower than the ambient/stable conditions (Tukey HSD, p = 0.010).

Figure 4 Calcium (A, B) and magnesium (C, D) concentrations measured with ICP-MS. In the carapace of lobsters in the reduced/high fluctuating treatment (7.67 ± 0.10), there was significantly less calcium (A) present than in lobsters in ambient/stable conditions. Lobsters in both fluctuating treatments had significantly less magnesium than those in ambient/stable treatments. There were no differences found in the cuticle of the abdominal segment, which had lower and less variable amounts of calcium but similar amounts of Mg to the carapace. Note different axes for panels (A, B) α=0.025 after a Bonferroni correction. Box boundaries are the first and third quartiles, whiskers are 1.5 times the interquartile range, the center line is the median, black and white diamonds are the mean, and outliers are solid dots. Values on the bottom of panels denote sample size.

Carapace magnesium was significantly lower, by 13%, in the reduced/high fluctuating treatment in comparison to both ambient/stable and the reduced/stable treatments (Tukey HSD, p ≤ 0.006; Figure 4C; Table 2). The reduced/low fluctuating treatment was also significantly decreased by 8% compared to the ambient/stable treatment (Dunn’s test, p = 0.011). Magnesium tended to be lower in the reduced/stable treatment than in ambient/stable (Dunn’s test, p = 0.047), and the two fluctuating treatments showed slightly different concentrations (Dunn’s test, p = 0.042), although these trends were not statistically significant, as the Bonferroni correction for repeated measures was applied.

In the abdominal segment, there were no significant differences in the calcium (ANOVA: F_3,51= 2.21, p=0.098; Figure 4B), and magnesium (ANOVA: F_3,51= 1.36, p=0.266; Figure 4D) concentrations between treatments (Table 2).

Cuticle material properties

The carapace spine and rostral horn responded differently to treatment conditions. The carapace spine did not show a significant difference in hardness (ANOVA: F_3,37= 2.708, p=0.06; Figure 5A), although the carapaces that were exposed to reduced fluctuating pH treatments tended to show lower hardness than those in stable pH treatment, averaging 28% lower. Stiffness was not significantly different across treatments (ANOVA: F_3,37 =  1.376, p =0.27; Figure 5C; Table 2).

Figure 5 Material properties of the carapace spine (A, C) and rostral horn (B, D). The carapace spine did not show any differences in material properties among treatments. The horn significantly differed in hardness. Lobsters exposed to reduced/high fluctuating pH conditions had horns that were less hard than those exposed to stable pH conditions (both ambient and reduced pH). Box boundaries are the first and third quartiles, whiskers are 1.5 times the interquartile range, the center line is the median, black and white diamonds are the mean, and outliers are solid dots. Values on the bottom of panels denote sample size.

The horns differed significantly in hardness among treatments (ANOVA: F_3,36 = 7.307, p<0.001) (Figure 5B; Table 2). Horns from lobsters exposed to the two stable pH treatments were statistically the same. However, horn samples from the reduced/high fluctuating treatment were significantly less hard than those from either of the stable pH treatments by 15 - 22% (p ≤ 0.045). The reduced/low fluctuating treatment horns were lower, but not significantly different, from those in the ambient/stable pH treatment (p=0.050), although they were significantly lower than those in reduced/stable pH treatment (p=0.004). There were no significant differences in horn stiffness among the treatments (ANOVA: F_3,36 = 0.077, p=0.972; Figure 5D; Table 2).

Antenna stiffness

Mean antennal flexural stiffness (EI) at the more flexible distal location did not differ among lobsters in any treatments (Kruskal-Wallis rank sum test, χ^{2 =} 3.52, df=3, p=0.318) (ambient/stable: 0.004 ± 0.002 Nm²; reduced/stable: 0.003 ± 0.002 Nm²; reduced/low fluctuating: 0.003 ± 0.001 Nm²; reduced/high fluctuating: 0.004 ± 0.002 Nm²) (Figure 6). At the stiffer proximal location, EI was significantly lower in lobsters exposed to reduced/stable pH treatments than lobsters exposed to reduced/high fluctuating pH (Welch’s ANOVA: F_3,22.55= 3.15, p=0.045; Games-Howell posthoc p=0.035). There were no significant differences among the other treatments (p≥0.20), although EI was slightly lower in antennae from the reduced/stable and reduced/low fluctuating treatment than those in ambient conditions (ambient/stable: 0.023 ± 0.012 Nm²; reduced/stable: 0.018 ± 0.004 Nm²; reduced/low fluctuating: 0.020 ± 0.004 Nm²; reduced/high fluctuating: 0.026 ± 0.010 Nm²).

Figure 6 Flexural stiffness of the antennae at proximal (striped shades) and distal (solid shades) locations. The low stiffness at the distal location was similar across all treatments. While antennae from the reduced/high fluctuating treatment were stiffer at the proximal base than lobsters in the reduced/stable pH treatment, there were no significant differences in flexural stiffness from the ambient/stable treatment. Box boundaries are the first and third quartiles, whiskers are 1.5 times the interquartile range, the center line is the median, and black and white diamonds are the mean. Values on the bottom of panels denote sample size.

Discussion

Predation is thought to be the largest source of mortality for juvenile California spiny lobster, P. interruptus, until they reach a size refuge against most predators at four to five years of age (Engle, 1979). Mortality rates during this critical life stage could be exacerbated if anti-predation strategies and defenses are compromised by changing ocean carbonate chemistry. OA-like conditions are known to affect the composition of calcified materials (Chan and Connolly, 2013; Kroeker et al., 2013; Davis et al., 2021), potentially disrupting the mechanical properties, though this link is not often studied in crustaceans (Taylor et al., 2015). Furthermore, mineralization effects are generally measured on one or two crustacean body regions, typically the carapace and chelae (Dutil et al., 2000; Boßelmann et al., 2007; Amato et al., 2008), but the strong localization of mineralization and specialization of different structures can lead to differential responses to OA. Lastly, exposure to stable reduced pH experimental conditions does not adequately assess the vulnerability of animals that live in dynamic pH environments (Wahl et al., 2015). Here, we tested the hypotheses that reduced stable pH conditions would affect the mineralization and mechanical properties of exoskeletal defense structures, that the effects would differ in magnitude across structures, and that fluctuating reduced pH conditions would mitigate these effects. After exposing juvenile P. interruptus to diurnal conditions that mimic their natural environment under near-future seasonal upwelling and/or OA, we found that most exoskeletal properties were unaffected by the treatment conditions, but that the integrity of some important exoskeletal defenses were compromised in variable ways and not as hypothesized. In particular, changes in mineralization or material properties were observed across key defensive structures (i.e., carapace, rostral horns, and antennae) in ways that were site specific, but absent of typical mineral-mechanical properties relationships. These effects were induced primarily by fluctuating reduced pH conditions and, to a limited extent, stable reduced pH conditions.

Exoskeletal defense integrity

Rostral horns

Spiny lobster rostral horns must be sufficiently hard and stiff to anchor animals in their dens when predators attempt to expel them. The spine-like shape of the horns is effective at locking into the interstices of hard substrate, while high hardness and stiffness limit deformation and yield stress. Structurally, rostral horns differ from the typical cuticle template by having only two layers, an outer and inner core (Lowder, 2019). The outer layer of the horn is poorly mineralized in comparison to typical crustacean cuticle, with just 1.4% Ca and 0.3% Mg. Neither Ca or Mg percentages were affected by reduced pH conditions in either layer, but hardness of the horn tips was lower in the reduced pH fluctuating conditions (both low and high fluctuating). The magnitude of this change in hardness is slight (<0.1 GPa, or 15%) compared to the decrease in microhardness observed in the chelae of juvenile blue and red king crabs (~44 VHN, or 46%) (Coffey et al., 2017). While red king crabs exhibited this decrease after exposure to pH 7.8, blue king crabs did so under a more extreme pH of 7.5; however, both species were exposed for a longer period of time (1 year) during which these effects manifested. It is unknown if the magnitude of change observed in spiny lobsters after three months of exposure would result in notable effects on lobster anchoring ability, but it is possible that longer exposure encompassing exoskeleton regeneration through one or more molt cycles could result in more dramatic effects.

We could not separate the horn layers for analysis of Ca and Mg content using the ICP-MS method, though EDX is sufficient for documenting differences in elemental composition (Fabritius et al., 2012; Taylor et al., 2015). The rostral horns have high amounts of halogens, such as chlorine, in the outer ring (Lowder, 2019), similar to the dactyl tips of crabs that scrape algae from rocks (Cribb et al., 2009). The replacement of calcium carbonate with halogenated compounds provides structures with greater toughness rather than hardness, which is beneficial for avoiding wear and breakage on rocky surfaces (Cribb et al., 2009; Schofield et al., 2009). This would be functionally important for the rostral horns that must lodge into the rocky underhangs of their shelter. It is interesting though that horn hardness was affected by reduced, fluctuating pH conditions despite its uncalcified state, hinting at a potentially complex interplay of elements under different ocean carbon chemistry parameters.

Antenna

Lobsters use their antennae as cantilevers to push away animals and maintain their distance (Kanciruk, 1980; Spanier and Zimmer-Faust, 1988). Sufficient flexural stiffness of the antennae is necessary to produce and maintain forces that hold predators at bay. In general, antennal cuticle has a thinner exocuticle layer compared to carapace cuticle, but both have similar Ca and Mg content. Flexural stiffness of the antenna is derived by the greater second moment of area associated with its hollow, tubular design, but it is also dependent on cuticle thickness and stiffness. Neither cuticle thickness nor mineral content of the antennae were affected by reduced pH conditions, unlike the carapace. Though there were no detectable changes in structure and composition, stiffness at the proximal base of the antennae in reduced/stable pH conditions was significantly reduced, by 31%, relative to the reduced pH/high fluctuating conditions. The mechanism underlying this change in bending stiffness is unclear, but may stem from weakened articulations between antennal annuli. Antennae breakage tends to occur in these regions (Lowder, pers. obs.), which warrants detailed morphological analysis.

Carapace

Should lobsters fail to keep predators at bay, their last line of exoskeletal defense is their armored carapace. The spiny lobster carapace is covered with thick, flattened spines and is more heavily mineralized than the abdominal segments (Lowder, 2019).

These attributes enhance hardness and stiffness, consequently providing greater resistance to the piercing and crushing forces of predator attacks, such as by the California sheephead. The abdomen, in contrast, facilitates the tail flip escape response and is not exposed to predators when the animal is sheltering. While there was no change to material properties after medium-term exposure to reduced pH/fluctuating conditions, carapace mineralization (both Ca and Mg) decreased and there was no change within the abdominal segment. The higher concentration of mineral as calcium carbonate and magnesian calcite in the carapace may have made it more favorable to draw from for potential hemolymph buffering under reduced pH conditions (Henry et al., 1981; Spicer et al., 2007). Thus, structures with greater mineral availability, and which benefit the most from the hardness that mineralization provides, may be more vulnerable to losing it in reduced pH/high pCO₂ conditions.

Discrepancies between elemental composition measured with ICP-MS and EDX are likely due to differences in the sensitivities of the techniques and/or in the cuticle regions sampled. Here, ICP-MS analysis revealed decreases to mineralization in reduced pH/fluctuating conditions, whereas no decrease was identified by EDX. Measuring elements over a 2-dimensional surface with EDX, compared to a volume of material with ICP-MS, captures more variation within and between cuticle samples, introducing greater variability, as observed in our data. EDX is also sensitive to surface morphology and sample preparation, which affects electron scattering and element detection. Furthermore, using EDX we only measured elements across regions of the endo- and exo-cuticle layers, whereas with ICP-MS we measured all layers, including the epicuticle. The distribution of mineral across cuticle layers in spiny lobsters is not well-documented, but the American lobster cuticle exhibits a defined calcitic layer near the outermost region of the exocuticle, bordering the epicuticle (Al-Sawalmih et al., 2008; Kunkel, 2013). If spiny lobster cuticles are similar, then changes to the mineralization of this layer would more likely be detected through our ICP-MS measurements of the whole cuticle, rather than our EDX sampling that did not capture this particular region.

Tangled links between composition and properties

Among mineralized organisms, changes to Ca and Mg content after exposure to reduced pH are assumed to belie changes to the mechanical function of calcified structures, but the limited studies testing this link so far defy expectations. In one study, increases in Mg content in the raptorial appendage of mantis shrimp exposed to ocean acidification-like conditions did not result in significant changes to the hardness or stiffness of the cuticle (deVries et al., 2016). In another study, Ca content increased in the carapace and chela of blue king crabs and red king crabs, respectively, but only the chela mechanical properties in both species changed (Coffey et al., 2017). For tanner crabs, an increase in carapace Mg and decrease in Ca did not result in any changes to microhardness, but chelae without significant alterations in mineralization experienced reductions in microhardness (Dickinson et al., 2021). In the present study, both Ca and Mg content decreased in the carapace of California spiny lobsters, but there was no significant change in cuticle hardness or stiffness. In all of these studies of crustaceans, changes in mineral content were less than 20%, which may not be of sufficient magnitude to induce measurable changes in nanomechanical properties. Several studies have established coarse correlations between the degree of mineralization and mechanical properties in crustacean cuticle (Currey at al., 1982; Raabe et al., 2005; Taylor and Kier, 2006; Sachs et al., 2006; Fabritius et al., 2012; Weaver et al., 2012; Rosen et al., 2020), but its complex structure and variation among species make it difficult to parse out fine-scale relationships and potential thresholds.

In contrast, cuticle mechanical properties may change independently from adjustments to mineral content, as observed here in the rostral horns of the spiny lobster that decreased in hardness without any corresponding changes to mineral content. Both cuticle structure (i.e., total thickness, layer thickness and stacking density) and mineral phase contribute to cuticle mechanical properties (Raabe et al., 2005). Changes to the former have not been observed in most crustaceans exposed to reduced pH conditions (Taylor et al., 2015; deVries et al., 2016; Lowder et al., 2017; Coffey et al., 2017), including species whose mechanical integrity was compromised (Coffey et al., 2017). Structural level changes, thus, are unlikely to be the mechanism driving observable effects in mechanical properties. The present study, as well as other crustacean OA studies, implicitly assume that the protein-chitin matrix of the cuticle remains unchanged, so that changes to the mineral:matrix ratios governing stiffness are driven by mineral content alone. What has yet to be widely observed in OA research are potential changes to the mineral phase in crustaceans. Crustaceans can manipulate calcium carbonate forms and use amorphous calcium carbonate, magnesian calcite, and calcite (Huner et al., 1979; Vigh and Dendinger, 1982; Addadi et al., 2003; Dillaman et al., 2005; Boßelmann et al., 2007), which may allow species to compensate for mineralization changes due to OA while maintaining structural function. Indeed, Tanner crabs that were post terminal molt appeared to shift amorphous calcium carbonate to calcite after two years of exposure to reduced pH conditions (Dickinson et al., 2021).

Crustacean cuticle is built with multiple levels of hierarchical structure (Raabe et al., 2005) that adds complexity to the large-scale link between structure, composition, and mechanical function under OA-like conditions. This link becomes clouded in other marine calcifiers as well, including mussels that develop more brittle and less stiff shell regions without any direct change in structure or composition (Fitzer et al., 2015). Consequently, direct relationships between compositional changes under OA-like conditions and material properties have not been consistently identified, and the perceived implications of observed changes to mineralization should be tested when possible.

Reduced pH affects intermolt exoskeleton

The effects of OA-like conditions on crustacean mineralization should be maximally expressed during molting, when animals shed and rapidly build a new exoskeleton. Mineral deposition peaks during the early postmolt period, but is otherwise considered stable throughout the duration of intermolt, except when mineral is potentially drawn away to be stored in gastroliths prior to the next molt (Shechter et al., 2008). As only two individuals (one from the reduced/stable pH treatment and one from the reduced/high fluctuating pH treatment) showed indications of entering the premolt stage (softened exoskeleton and apolysis) at the conclusion of the experiment, variations in molt stage do not likely explain variation in measurements. Juvenile California spiny lobsters of the size used in this study molt annually (California Department of Fish and Wildlife, 2004), and only a few lobsters from each treatment molted over the three month duration of this experiment. Yet, significant changes in mineralization and material properties occurred in the absence of molting, which has also been observed in Tanner crabs past their terminal molt (Dickinson et al., 2021). We found significant decreases in both Ca and Mg in the carapace, hardness in the rostral horns, and flexural stiffness in the antennae of intermolt lobsters, demonstrating that the fully-formed exoskeleton of intermolt animals is sensitive to reduced pH conditions.

Crustaceans use bicarbonate to regulate internal acid-base chemistry under external high pCO₂ conditions and the carbonate minerals in the cuticle can be one source of the ion (Henry et al., 1981). A significant decrease in calcium concentration in the carapace was only found in the reduced/high fluctuating treatment, potentially indicating that the lowest pH values of the experiment reached in that treatment at night (7.57) may have triggered calcium carbonate dissolution from the carapace to use for buffering. Magnesium concentration showed significant decreases from ambient conditions in all three reduced pH treatments. Magnesium is likely in the form of magnesian calcite, a more soluble form of calcite, and may have been drawn out first (Davis et al., 2000).

In other juvenile crustaceans, mineralization is often not affected by exposure to pH conditions expected within the next century. For example, there was no change in either [Mg] or [Ca] at moderate (0.25 pH unit) decreases in pH for the early juvenile stage of the American lobster, H. americanus, after five weeks of exposure (Menu-Courey et al., 2018). Juvenile blue swimming crabs and hermit crabs exposed to reduced pH conditions incurred no changes in Ca content, but hermit crabs exhibited small increases in Mg (Glandon et al., 2018; Ragagnin et al., 2018). The above results were all found in animals that had molted during the experiment. In Homarus gammarus that didn’t molt, there was no change in Mg and Ca after five weeks of exposure to a 0.4 unit decrease in pH (Small et al., 2016), although longer-term exposure might induce detectable changes.

Fluctuating pH conditions may exacerbate effects of reduced pH

Our results indicate that future scenarios of fluctuating pH conditions are not beneficial to P. interruptus, despite the current natural fluctuations they experience in kelp forests on diurnal and seasonal scales (Frieder et al., 2012; Kapsenberg and Hofmann, 2016; Kekuewa, 2020). While each reduced pH treatment had the same mean pH value, fluctuations above this value were hypothesized to provide a short, daily refuge period when processes like calcification can mainly be carried out (Wahl et al., 2018). Indeed, we found distinct differences in exoskeletal defenses between fluctuating and stable reduced pH conditions, but they were contrary to expectations. Mineralization (e.g., carapace Mg concentration) and material properties (e.g., rostral horn hardness) were significantly reduced in fluctuating conditions, but were unaffected by stable pH exposure. At the same time, antennae flexural stiffness was significantly lower in reduced, stable pH conditions.

Like some organisms, juvenile spiny lobsters may find fluctuating low pH conditions stressful; fluctuating reduced pH conditions raised oxidative stress levels in the mussel Mytilus edulis (Mangan et al., 2017). In this study, negative effects tended to be greater under reduced pH with high fluctuations than with low fluctuations. Potentially, the negative effects may result from the pH dropping below a threshold value at night around 7.55-7.60 pH units, which lobsters do not appear to encounter frequently, if at all, in the kelp forest habitat where they were collected (Frieder et al., 2012). Thus, OA may breach the sensitivity thresholds of organisms in fluctuating environments more rapidly than in relatively stable habitats, necessitating consideration of dynamics when developing predictions for when impacts on species and ecosystems will become more apparent.

Conclusions

In this study on juvenile California spiny lobsters, we found that fluctuating reduced pH conditions affect the mineral composition and mechanical properties of some of their exoskeletal defense structures, thereby adding P. interruptus to the list of species sensitive to OA-like conditions. This is the first study to establish spiny lobsters’ exoskeleton sensitivity to reduced pH, adding to the evidence of sensory impacts to this group (Ross and Behringer, 2019; Gravinese et al., 2020; Boco et al., 2021) that also has implications for predator defenses. The effects were variable across exoskeletal structures, demonstrating localized responses within an individual that should be accounted for when assessing vulnerability across crustacean species. While the magnitudes of these effects were not appreciably large, it is noteworthy that significant changes to the exoskeleton occurred during the intermolt phase, when the cuticle is generally considered to be stable. This observation reveals that the crustacean exoskeleton can be continuously responsive to environmental pH/pCO₂ conditions, and not just during the peak post molt period. Furthermore, the correlation between mineralization and mechanical properties is not always straightforward under OA-like conditions such that changes in mineralization are not necessarily a good indicator of impacts on mechanical properties. Thus, assumptions about mechanical integrity and function cannot be made based on mineral content alone.

Worldwide, 19 Panulirus species (Ptacek et al., 2001) live in shallow, dynamic habitats, such as coral reefs, rocky reefs, and kelp forests. Our study on the temperate P. interruptus shows that the natural daily fluctuations in pH typical of these habitats does not necessarily afford spiny lobsters with a helpful respite from greater pH variations associated with climate change predictions. Rather, reduced/fluctuating pH conditions induced exoskeletal changes not observed under stable reduced pH. Dynamic pH conditions should therefore be considered when assessing the vulnerabilities of other species of spiny lobsters, as well as other crustaceans occupying these habitats. Under future ocean conditions, the California spiny lobster may engage with predators using altered exoskeletal defenses that could potentially affect the outcomes of these interactions.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Data and metadata are available through PANGAEA: Lowder, Kaitlyn; deVries, Maya S; Hattingh, Ruan; Day, James M D; Andersson, Andreas J; Zerofski, Phillip; Taylor, Jennifer (2022): Exoskeletal predator defenses of juvenile California spiny lobsters (Panulirus interruptus) are affected by fluctuating ocean acidification. PANGAEA, https://doi.org/10.1594/PANGAEA.945362.The R code is accessible at:https://github.com/kblowder/Spiny-lobster-exoskeleton-OA.

Author contributions

KL and JT conceptualized the research and wrote the manuscript. KL carried out the experimental study, water chemistry measurements, hardness and stiffness nanoindentation measurements, SEM and EDX measurements, statistical analyses, and aided in animal collection. JT carried out flexural stiffness and ICP-MS measurements. MdV helped carry out SEM and EDX measurements, and RH aided in ICP-MS measurement and analysis. JD contributed to ICP-MS analysis. AA aided in designing water chemistry treatments and analyzing water chemistry results. PZ led animal collection and aided in experimental setup. JD, AA, and MdV edited the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was partially supported through access and utilization of the NanoEngineering Materials Research Center (NE-MRC) and the Nano3 facility at UC San Diego. This work was supported by the National Science Foundation Graduate Research Fellowship Program, the UC San Diego Frontiers of Innovation Scholars Program, and the PADI Foundation to KL.

Acknowledgments

Many volunteers were critical for carrying out this work: Cierra Kelly, Jack Shurtz, Duc Tran, Steven Ly, Eric Young, Grace Chan, Kyle Perdue, Zoe Sebright, and Alex Hill. Travis Buck, Kevin Hovel, Brett Pickering, Lyndsay Sutterley, Matt Costa, and Georgie Zelanak either provided traps, helped set them up, or helped recover trapped lobsters. Sebastian Davis, Yaamini Venkataraman, and Summer Webb helped collect water samples from La Jolla seagrass beds. David Cervantes from the Dickson Laboratory provided the sampling protocol and analyzed water samples. Scripps Coastal & Open Ocean BiogeochemistrY lab members, including Alyssa Griffin, Sam Kekuewa, and Ariel Pezner, helped train and troubleshoot while analyzing experimental water samples. Sabine Faulhaber provided training at the NanoEngineering Materials Research Center (NE-MRC). We also thank the editor and reviewers for their thoughtful comments. Portions of this work are contained in the dissertation of KL (Lowder, 2019).

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.

Supplementary material

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

References