Intra-Specific Variation Reveals Potential for Adaptation to Ocean Acidification in a Cold-Water Coral from the Gulf of Mexico

Kurman, Melissa D.; Gómez, Carlos E.; Georgian, Samuel E.; Lunden, Jay J.; Cordes, Erik E.

doi:10.3389/fmars.2017.00111

ORIGINAL RESEARCH article

Front. Mar. Sci., 04 May 2017

Sec. Global Change and the Future Ocean

Volume 4 - 2017 | https://doi.org/10.3389/fmars.2017.00111

This article is part of the Research Topic The Pathway to Solutions: New Frontiers in Climate Change Adaptation & Mitigation View all 13 articles

Intra-Specific Variation Reveals Potential for Adaptation to Ocean Acidification in a Cold-Water Coral from the Gulf of Mexico

$\r\nMelissa D. Kurman*&#x;$ Melissa D. Kurman¹^*^†

¹Department of Biology, Temple University, Philadelphia, PA, USA
²Marine Conservation Institute, Seattle, WA, USA
³Department of Biological Sciences, Lock Haven University of Pennsylvania, Lock Haven, PA, USA

Ocean acidification, the decrease in seawater pH due to the absorption of atmospheric CO₂, profoundly threatens the survival of a large number of marine species. Cold-water corals are considered to be among the most vulnerable organisms to ocean acidification because they are already exposed to relatively low pH and corresponding low calcium carbonate saturation states (Ω). Lophelia pertusa is a globally distributed cold-water scleractinian coral that provides critical three-dimensional habitat for many ecologically and economically significant species. In this study, four different genotypes of L. pertusa were exposed to three pH treatments (pH = 7.60, 7.75, and 7.90) over a short (2-week) experimental period, and six genotypes were exposed to two pH treatments (pH = 7.60 and 7.90) over a long (6-month) experimental period. Their physiological response was measured as net calcification rate and the activity of carbonic anhydrase, a key enzyme in the calcification pathway. In the short-term experiment, net calcification rates did not significantly change with pH, although they were highly variable in the low pH treatment, including some genotypes that maintained positive net calcification in undersaturated conditions. In the 6-month experiment, average net calcification was significantly reduced at low pH, with corals exhibiting net dissolution of skeleton. However, one of the same genotypes that maintained positive net calcification (+0.04% day⁻¹) under the low pH treatment in the short-term experiment also maintained positive net calcification longer than the other genotypes in the long-term experiment, although none of the corals maintained positive calcification for the entire 6 months. Average carbonic anhydrase activity was not affected by pH, although some genotypes exhibited small, insignificant, increases in activity after the sixth month. Our results suggest that while net calcification in L. pertusa is adversely affected by ocean acidification in the long term, it is possible that some genotypes may prove to be more resilient than others, particularly to short perturbations of the carbonate system. These results provide evidence that populations of L. pertusa in the Gulf of Mexico may contain the genetic variability necessary to support an adaptive response to future ocean acidification.

Introduction

Climate change is dramatically altering the Earth's terrestrial and aquatic environments, with ocean acidification expected to elicit some of the more severe alterations to marine environments (Keeling et al., 1976; Parmesan and Yohe, 2003; Feely et al., 2004; Hoegh-Guldberg and Bruno, 2010). Carbon dioxide is currently being emitted into the atmosphere at a rate of approximately 49 ± 4.5 Gt CO₂ year⁻¹, and approximately one-third of these emissions are ultimately absorbed by the world's oceans (Sabine et al., 2004; IPCC, 2013). The absorption of this excess CO₂ drives a series of chemical reactions that alter the carbonate chemistry of the seawater causing significant reductions in seawater pH and carbonate ion concentration ( ${CO}_{3}^{2 -}$ ) (Kleypas and Langdon, 2006). It is estimated that ocean pH will drop by 0.3–0.5 pH units within the century, with estimates as high as 0.7 pH units (Caldeira and Wickett, 2003). Consequently, ocean acidification lowers the saturation state (Ω) of calcium carbonate, such that some marine calcifiers are unable to maintain calcification rates that exceed dissolution, or do so at much greater energetic cost. This can compromise the biological and ecological functions of these species, including reduced extent and density of reef framework (Gattuso et al., 1998; Kleypas et al., 1999; Langdon and Atkinson, 2005), failed larval development (Kurihara, 2008), impairment of predator-prey interactions (Cripps et al., 2011; Gaylord et al., 2014), and changes in acid-base regulation (Pörtner et al., 2004).

Cold-water corals are one of the most important foundation species in the deep sea, providing habitat for a large number of associated species including many commercially important fish and invertebrates (Rogers, 1999; Roberts et al., 2006). In recent years, there has been great interest in predicting the physiological response of cold-water corals to future levels of ocean acidification, mainly because these species inhabit deep-sea waters where aragonite saturation states (Ω_ar) are naturally low, and thus will be experiencing corrosive waters sooner than shallow-water organisms (Orr et al., 2005; Guinotte et al., 2006). There is evidence that deep-sea corals in the Gulf of Mexico live close to the aragonite saturation horizon (Ω_ar ~ 1) with [ ${CO}_{3}^{2 -}$ ] about 94 μmol Kg⁻¹ (Georgian et al., 2016a). Moreover, these species naturally exhibit slow growth and calcification rates, which range between ~2 to 9 mm year⁻¹ (Brooke and Young, 2009; Larcom et al., 2014; Lartaud et al., 2014) making them particularly vulnerable to the effects of future ocean acidification. Nevertheless, recent evidence from experimental studies suggests that cold-water corals can be resistant to low aragonite saturation states, and may even be able to calcify in undersaturated conditions (Form and Riebesell, 2012; Maier et al., 2013; Hennige et al., 2015), though at a reduced rate (Maier et al., 2009; Lunden et al., 2014a). The ability to maintain calcification under low-pH conditions has been explained by the capacity of cold-water corals to upregulate the pH of their internal calcifying fluid by up to 0.8 units (McCulloch et al., 2012; Wall et al., 2015).

The rapid pace of ongoing and projected ocean acidification is expected to outpace the natural capacity of marine species to survive through acclimatization or adaptation (Hoegh-Guldberg, 2014); this may prove especially true for cold-water corals which are generally long-lived, slow growing, and have a low fecundity. However, a growing body of literature supports the concept that acclimation or genetic adaptations may confer resilience to reduced pH conditions (e.g., Langer et al., 2009; Parker et al., 2011; Pistevos et al., 2011; Sunday et al., 2011; Pančić et al., 2015; Vargas et al., 2017). The evolution of adaptive capacities within a population may occur from novel mutations, but will likely depend heavily on the standing variation within the gene pool. Therefore, identifying existing variation in the response of single species to acidified seawater is a critical first step toward assessing their adaptive potential to future ocean acidification. Standing variation to acidification has been investigated in a number of marine organisms, including bryozoans (Pistevos et al., 2011), oysters (Parker et al., 2012), coccolithophores (Langer et al., 2009; Lohbeck et al., 2012), and sea urchins (Kelly et al., 2013). Notably for cold-water corals, Lunden et al. (2014a) found that Lophelia pertusa individuals from the Gulf of Mexico exhibited a highly variable calcification response to reduced pH conditions, suggesting that certain acidification-resistant genotypes may provide such an adaptive capacity.

In the present study, we performed both a short-term (2-week) and long-term (6-month) assessment of the physiological performance of L. pertusa collected in the Gulf of Mexico after exposure to different levels of reduced pH. In order to assess the differential response to altered carbonate chemistry at a genotype level, six different genotypes were included in the long-term experiment, and four of these six were also used in the short-term experiment. We hypothesized that overall, decreased pH would elicit a negative response in the physiological performance of L. pertusa, but further hypothesized that some genotypes may be more resilient to the effects of ocean acidification. If upheld, this would imply a capacity for an evolutionary response to ocean acidification and an increased probability of persistence for scleractinian corals as they face future global climate and ocean change.

Materials and Methods

Study Species and Sample Collection

Lophelia pertusa is a scleractinian cold-water coral that is found throughout the world's oceans, typically between 200 and 1,000 m (Rogers, 1999), although it occurs as shallow as 40 m in Norwegian fjords due to the upwelling of deep, cold water (Strømgren, 1971). In the Gulf of Mexico, L. pertusa is most commonly observed between depths of 300 and 600 m, where it creates habitat that greatly increases local diversity (Cordes et al., 2008). Globally, its distribution is best explained by water temperature (Davies and Guinotte, 2011), while in the Gulf of Mexico, its distribution can best be predicted by depth, local topography, and the availability of hard substrata (Georgian et al., 2014).

L. pertusa were collected in the Gulf of Mexico in April and May 2014 on the R/V Atlantis using the DSV Alvin (Woods Hole Oceanographic Institution). All corals were collected at similar depths, temperature and carbonate chemistry conditions (Table 1) from the Viosca Knoll 826 (VK826) site (Cordes et al., 2008), named for the oil lease block governed by the U.S. Bureau of Ocean Energy Management (BOEM). Each collection was carefully placed in a separate temperature-insulated compartment on the submersible, and kept separate throughout the relocation and experiments. All live corals were kept in natural seawater and transported in insulated containers from Gulfport, Mississippi to Temple University, Philadelphia, Pennsylvania.

TABLE 1

Table 1. Collection sites for each of the different genotypes used in the present study with its respective in situ carbonate chemistry conditions.

At the time of arrival, all corals were kept in a 550 L recirculating aquarium system with artificial seawater (ASW) prepared using Instant Ocean^® (aquarium sea-salt mixture) and maintained at a temperature of ~8°C, salinity 35 ppt, and total alkalinity (TA) ~2,300 μmol kg⁻¹ (for full description of the system see Lunden et al., 2014b). The corals were fed 20 mL of artificial MarineSnow^® (Two Little Fishies) 3 days a week, as per comparable experimental procedures (Lunden et al., 2014a,b; Hennige et al., 2015; Georgian et al., 2016b). Corals were allowed to acclimate to laboratory aquaria conditions for a period of approximately 3 months prior to the start of experiments.

Paired water samples were collected directly above (<2 m) each of the six coral collection sites using vehicle-mounted Niskin bottles. Samples were measured for pH on the total scale (pH_T) using an Orion 5 Star pH meter and glass electrode (ROSS Ultra pH/ATC Triode 8107BNUMD). Total alkalinity was measured by acid-titration on an autotitrator (Mettler-Toledo DL15, 0.1 mol L⁻¹ HCl), with accuracy checked by measuring certified reference material (Dickson Lab, batches 138 and 141; Dickson et al., 2007). CO2Calc software (Robbins et al., 2010) was used to calculate pCO₂, ${HCO}_{3}^{-}$ , ${CO}_{3}^{2 -}$ , and Ω aragonite (Ω_ar) using pH and total alkalinity as input variables with the dissociation constants for boric acid and K1 and K2 from Mehrbach et al. (1973) refit by Dickson and Millero (1987), KHSO₄ from Dickson (1990), total boron from Lee et al. (2000), and pH on the total scale (pH_T). In situ salinity, pressure, and temperature were measured using a vehicle-mounted CTD (Seabird SBE 49). For the full methodology of water sampling and analysis, see Georgian et al. (2016a).

Determination of Coral Genotypes

Corals were collected from six discontinuous reef patches within VK826 to increase the probability that different genotypes were sampled. Previous work indicates relatively low levels of clonality among distinct coral colonies at this site as compared to L. pertusa populations in the North Atlantic (Waller and Tyler, 2005; Lunden et al., 2014a). To ensure that distinct genotypes were used in our experiments, each coral fragment was genotyped according to methods used previously in a similar study of L. pertusa (Lunden et al., 2014a). Briefly, coral genotypes were determined using 10 L. pertusa microsatellite loci (Morrison et al., 2008; Molecular Ecology Resources database entry 51059) from total DNA extracted from coral polyps. Of the six collections used in the live coral experiments, all of the microsatellite loci amplified, except for one collection in which only 6 of the 10 amplified. In all pairwise comparisons, the highest number of shared microsatellite genotypes was 2, the average was 0.8, and in 7 of the 15 pairwise comparisons, there were no single-locus genotypes shared between the two individuals. These results indicate that each of the samples represents a distinct genotype.

Manipulation of Seawater Chemistry

Because Instant Ocean^® at salinity 35 yields seawater with a total alkalinity of ~3,600 μmol kg⁻¹, 12.1 N HCl was added to previously prepared artificial seawater (ASW) in order to reduce the total alkalinity to ~2,300 μmol kg⁻¹ before adjusting the pH (Lunden et al., 2014b). The manipulation of the pH was accomplished by bubbling pure CO₂ gas into the different treatment tanks using an automated CO₂ injection system (American Marine Inc, PINPOINT pH Monitor). To reduce pH, pure CO₂ was mixed into a reaction chamber that dissolved the CO₂ into seawater before being delivered into the respective tanks. If pH levels fell below the desired threshold, ambient air was automatically bubbled into the tank. The pH meters underwent a two-point Tris-HCl and AMP-HCl calibration weekly (Dickson et al., 2007). In both experiments the pH_T of the experimental tanks was gradually brought down to the desired treatment conditions at a rate of ~0.1 pH units day⁻¹. Total alkalinity (TA) was measured twice a week as described above. Total alkalinity and pH_T were used to calculate pCO₂, ${HCO}_{3}^{-}$ , ${CO}_{3}^{2 -}$ , and Ω_ar in CO2Calc software (Robbins et al., 2010) as described above. Salinity was measured daily using a handheld refractometer, and temperature was recorded every 5 min using a temperature logger (Onset HOBO Pendant^®). A 10% water change was performed weekly for the long-term experiment, and twice weekly for the short-term experiment.

Short-Term Experiment

The short-term experiment began in August 2014. Three pH_T treatments were selected: 7.90, which approximates the average in situ pH_T at the VK826 collection site (Lunden et al., 2013); 7.75, which results in seawater near the aragonite equilibrium point (Ω_ar = 1); and 7.6, which generates undersaturated seawater with respect to aragonite (Ω_ar < 1). Each treatment contained 12 fragments from four different genotypes (36 fragments in total). These treatments were established in three separate 50 L tanks housed in a cold room, which was maintained at 8°C. Each of these fragments was weighed at the beginning and end of the 2-week pH exposure to determine net calcification rates.

Long-Term Experiment

In July 2014, the six collections were randomly split between two identical 550 L recirculating units (fully described in Lunden et al., 2014b). Two pH_T treatments were selected: 7.90 as the control, and 7.60 as the treatment. The pH was gradually reduced in August 2014 to treatment conditions (control = pH 7.90, acidified = pH 7.60). Fragments (n = 4) from each genotype (n = 6) of five to eight live polyps were set-aside in each tank for buoyant weighing every 2 weeks (48 fragments in total).

Net Calcification

The buoyant weight (BW) technique (Jokiel et al., 1978; Davies, 1989) was used to measure net calcification using a Denver Instruments SI-64 closed analytical balance (Precision of 0.1 mg). All fragments from each treatment were weighed before the gradual pH decline, when the pH reached the desired treatment level, after 14 days, and in the long-term experiment, subsequently every 2 weeks for a 6-month time period. The balance was closed within a plexiglass chamber to prevent any outside disturbance from affecting the weight. Each fragment was removed from its treatment tank and transported to the balance in a 4 L beaker without exposure to air. Fragments were gently attached to a monofilament fishing line, and suspended within the water from the bottom of the balance. After 10 s of stabilization, the weight was recorded. Each of the fragments in both the short- and long-term experiments was weighed in triplicate during each buoyant weight trial to minimize error.

Salinity and temperature of the water at time of weighing was recorded, and used to calculate the density of the water and to ensure that the treatment conditions were maintained throughout the weighing process. The dry weight of each fragment (W_a) was calculated using the measured buoyant weight (W_w), the density of the seawater during measurement (D_W), and the coral skeletal density (SD, 2.82 g cm⁻³) (Lunden et al., 2013).

\begin{array}{l} W_{a} = \frac{W_{w}}{1 - \frac{D_{w}}{S D}} \end{array}

Net calcification of L. pertusa was calculated as the change in weight over each 2-week interval, expressed as % net calcification d⁻¹. The average weight of the fragments at the start of the experiment was 15.81 ± 7.39 g. For the long-term experiment, net calcification rates were averaged monthly to yield a total of six time points. Percent net calcification per day was calculated by the equation described in Lunden et al. (2014a):

\begin{array}{l} N_{t} = 100 \times \frac{W_{w 2} - W_{w 1}}{W_{w 1} (T_{2} - T_{1})} \end{array}

where W_w2 and W_w1 is the final and initial standardized buoyant weight respectively, and T₁ and T₂ equal time 1 and time 2, respectively.

Carbonic Anhydrase Activity

At the beginning and end of the long-term experiment, coral fragments from each genotype were selected to assess the potential effects of ocean acidification on enzyme activity. Collected corals were frozen in liquid nitrogen and stored in −80°C until analyses were performed. Frozen coral polyps (2-3 polyps) from each genotype were extracted on liquid nitrogen, ground with a mortar to disrupt membrane cells, and homogenized with 3.5 mL of sterilized artificial seawater (ASW). The homogenate was centrifuged at 4°C for 3 min at 1,500 g. 200 μL was reserved for total protein quantification using bicinchoninic acid assay (Pierce™ BCA assay kit), and 1.5 mL of the remaining homogenate was transferred into a 7 mL scintillation vial and diluted 1:1 with 25 mM veronal buffer. The in-vitro assay for carbonic anhydrase was performed according the protocol described by Weis et al. (1989), and was measured by the decrease in pH resulting from the hydration of CO₂ to ${HCO}_{3}^{-}$ and H⁺ after the addition of substrate (CO₂−enriched distilled H₂O, pH = 3.5, 4°C). For each assay, 1 mL of buffered homogenate was transferred to a scintillation vial equipped with a magnetic stir bar and further diluted 1:1 with the same veronal buffer. The homogenate was maintained at ~3.5°C before running the assay. The mixture was constantly stirred, and 1 mL of substrate was added to the sample when the pH stabilized and temperature reached 13°C. From this time, the drop in pH was recorded every 10 s for a total of 120 s. As a negative control for CA activity, a subsample of the homogenate (~1.7 mL) was boiled for 5 min to denature the enzyme and treated as above. Enzyme activity was calculated as the difference of the decrease in pH between the sample with the active enzyme and the control with the denatured enzyme, and is given as enzyme units (EU) standardized to total protein content.

Data Analysis

For the short-term experiment, a two-way ANOVA was used with net calcification (change in mass as % d⁻¹) as the dependent variable, pH as a fixed factor (three levels: pH = 7.9, 7.75, 7.6) and genotype as a random effect (6 levels). Multiple pair-wise comparisons were then applied to investigate the differential response to pH at the genotype level, using fragments within genotype (n = 12) as replicates. When considered together, the net calcification rate data were not normally distributed, and there were no simple transformations that resulted in a normal distribution. However, when considered separately, the data within each pH level were normally distributed. Furthermore, ANOVA is known to be very robust with respect to violations of this assumption (Schmider et al., 2010), and therefore we proceeded with the analysis as planned.

In the long-term experiments, a three-way factorial ANOVA was performed with growth (% d⁻¹) as the dependent variable, pH and time as fixed factors and genotype as a random effect, using each individual fragment as a replicate (Quinn and Keough, 2002). Tukey's HSD test was used to test for pairwise differences in the response of each genotype to the different pH treatments, where fragments within genotypes (n = 24) were used as replicates. Additionally, to investigate potential differences in the amount of time that a given genotype could maintain positive calcification in undersaturated conditions, survival curves were generated using Kaplan-Meier “time to event” analysis using a Mantel-Cox log rank test followed by a post-hoc Wilcoxon exact test to determine if significant differences were found (Kaplan and Maier, 1958). In this case, the “event” in question was the change from positive to negative calcification rate, and the time when negative calcification was first observed in each of the fragments of a given genotype within a treatment was used as the data in the analysis.

Differences in enzyme activity of L. pertusa at the end of the long-term experiment were assessed with a two-way ANOVA with pH (two levels) and time (initial and final time) as fixed factors. If significant differences were found, genotype was further analyzed by multiple pair-wise comparison. All analyses were applied to the mean standardized values of enzyme activity among replicate fragments of each genotype. For all the analyses, results were considered significantly different with p < 0.05. All data are presented as average ± standard deviation unless otherwise indicated.

Results

Aquaria Conditions

Mean values of temperature, salinity and all carbonate chemistry parameters for the different treatments in both short and long-term experimental aquaria were relatively stable during the course of the two experiments. All parameters remained within a very narrow range around the target values of 8°C, 35 ppt, and pH values of 7.9, 7.75, and 7.6 (Table 2). These parameters, along with a total alkalinity near 2,300 resulted in Ω_ar values of approximately 1.5, 1, and 0.8 in the three treatment levels.

TABLE 2

Table 2. Summary of seawater chemistry parameters from each of the treatments for both long-and short-term experiments.

Short-Term Experiment

Net Calcification Rates

Corals within the control (pH = 7.9) and the slightly acidified treatments (pH = 7.75) displayed the highest net calcification rate averaging 0.004 ± 0.011 and 0.006 ± 0.003% d⁻¹ respectively (Table 3). In the lowest pH treatment, skeletal dissolution was observed at a rate of −0.033 ± 0.032% day⁻¹ (Table 3, Figure 1), although there was no significant effect due to pH [ANOVA F_{(2, 24)} = 1.03, P = 0.412] due to the high variance observed in net calcification rate measurements. There was also no significant difference in the overall net calcification rate of the different genotypes when averaged among treatments [ANOVA F_{(3, 24)} = 1.119, P = 0.413]. However, there was a significant interaction term between pH and genotype [ANOVA F_{(6, 24)} = 2.575, P = 0.045] (Table 4), indicating that different genotypes had different responses to the three pH treatments. If the two extreme values from G1 are dropped from this analysis, then the significance level of the interaction is P = 0.052, indicating that the hypothesis that different genotypes responded differently to changes in pH was still moderately supported. On average, genotype 1 (G1) and G3 experienced net dissolution (negative calcification) in the highly acidified treatment (pH = 7.6) with −0.144 ± 0.0156 and −0.0008 ± 0.001 % net calcification d⁻¹ respectively, whereas G4 experienced significant positive net calcification at a rate of 0.013 ± 0.001 % d⁻¹ in the same treatment (Figure 1). In particular, within the highly acidified treatment all three of the fragments from G4 had the highest positive calcification of any of the individual nubbins. However, pairwise post-hoc comparisons of net calcification rates indicated that the only genotype with a significantly different net calcification rate over the treatments was G1 in the low pH treatment (Tukey HSD, P < 0.034 for all 11 pairwise comparisons). This difference was mainly driven by the highly negative net calcification rate measured in two of the three fragments in undersaturated conditions (Figure 1). The third individual had a net calcification rate of +0.003% d⁻¹, which is within the range of the other fragments at higher pH.

TABLE 3

Table 3. Short- and long-term net calcification rates measured with the buoyant weight technique.

FIGURE 1

Figure 1. Short-term calcification response of individual genotypes of L. pertusa corals in the different ocean acidification conditions. G1 through G4 refers to each genotype. Values are given as mean ± S.D

TABLE 4

Table 4. Short-term experiment ANOVA.

Long-Term Experiment

Net Calcification Rates

The average net calcification rate in the control treatment (pH = 7.90) was 0.0057 ± 0.006% net calcification d⁻¹, while in the low pH (pH = 7.60) treatment it was −0.0078 ± 0.0043% net calcification d⁻¹. Although variability in the response of different individuals was evident, the overall trend of positive calcification at Ω_ar > 1 and dissolution at Ω_ar < 1 was still observed (Table 3, Figures 2A,B). All genotypes in the pH 7.6 treatment started to exhibit negative calcification rates as the experiment progressed, with net dissolution in 100% of the fragments by the end of the experiment (Figure 3). The three-way ANOVA indicated significant differences in the main effects of overall net calcification rate in the two treatments [F_{(1, 199)} = 49.09, P = 0.001], and average net calcification rate over time [F_{(5, 199)} = 8.44, P < 0.001], but no significant difference in overall net calcification rate among genotypes [F_{(5, 199)} = 1.01, P = 0.495]. In addition, the interaction terms between pH treatment and time [F_{(5, 199)} = 4.21, P = 0.006] as well as treatment and genotype were significant [F_{(5, 199)} = 2.73, P = 0.04] (Table 5).

FIGURE 2

Figure 2. (A) Long-term time-series of the calcification response of L. pertusa standardized by the initial weight for each time point, (B) calcification response standardized by initial weight t₀ at the beginning of the experiment. Values are given as mean ± S.D.

FIGURE 3

Figure 3. Long-term time-series of the calcification response by individual genotype standardized by initial weight for each time-point. Values are given as mean ± S.D.

TABLE 5

Table 5. Long-term experiment three-way factorial ANOVA with calcification as the dependent variable showing the effect of time, pH treatment, and genotype.

The treatment x genotype interaction is more complex and was explored using post-hoc pairwise comparisons. In this analysis, every pairwise comparison between the overall average net calcification rate of a given genotype in the pH 7.9 treatment and the same genotype's average net calcification rate at pH 7.6 was significant (Tukey HSD, P < 0.021 for 5 of the 6 comparisons), except for G5, which did not show a significant difference between treatments (Tukey HSD, P = 0.828). However, the net calcification rate of G5 in the low pH treatment was also not significantly different from any of the other genotypes in the low pH treatment (Tukey HSD, P > 0.645 for all 5 comparisons). The lack of significance in these tests can primarily be attributed to the high intra-genotype variance in net calcification rate, as one of the G5 fragments exhibited highly negative calcification at two of the time points in the low pH treatment (time 2 and time 4, Figure 3).

To further examine the potential differences in the response of genotypes to low pH and carbonate ion availability, a time-to-event survival analysis was used to determine which genotypes were capable of maintaining positive calcification the longest (Figure 4). Notably, out of the six genotypes being used in the experiment, two (G3 and G4) maintained positive net calcification (0.0115 and 0.00964 % g d⁻¹ respectively for the first 90 days) in the acidified (pH 7.6) treatment, significantly longer than the other genotypes (Kaplan-Meier survival analysis Mantel-Cox X² = 19.148, p < 0.01; post-hoc Wilcoxon exact test p < 0.05; Figure 4, Table 6).

FIGURE 4

Figure 4. Kaplan-Meier cumulative survival curves for the long-term experiment showing the month when the different genotypes reached negative net calcification (dissolution). Each color represents a different genotype, with no significant difference among pairwise comparisons of the times that they could maintain positive calcification indicated by the lines over the numbers in the key.

TABLE 6

Table 6. Post-hoc pair-wise comparisons for Kaplan-Meier analysis between different genotypes for low pH treatment in the long-term experiment.

Carbonic Anhydrase Activity

No significant differences were observed in the average carbonic anhydrase activity for corals grown in the pH 7.9 vs. pH 7.6 [F_{(1, 19)} = 0.32, P = 0.59]. The average activity at the initial time point was 0.010 ± 0.003 EU for corals in the pH 7.9 and 0.006 ± 0.003 EU for corals in pH 7.6 (Figure 5). After 6 months the average activity was similar for control conditions (pH 7.9) with an average 0.010 ± 0.006 EU, however, it increased in the low pH treatment (7.6) with an average 0.013 ± 0.012 EU (Figure 5), although this was not significant due to the high variability in the enzyme activity [F_{(1, 19)} = 0.32, P = 0.32]. The maximum average activity was measured at the low pH final time point for G2 0.031 EU, followed by G5 at the control conditions with a value of 0.020 EU (Figure 5).

FIGURE 5

Figure 5. Carbonic anhydrase activity for the long-term experiment shown for the average of all genotypes (general) and for individual genotypes (G1 through G6). T_i is the initial time point and T_f is the final time point after the six-month experiment. Values are given as mean ± S.D.

Discussion

The present study examined the potential effects of short- (2 weeks) and long-term (6 months) low-pH exposure on the physiological response of cold-water coral L. pertusa from the Gulf of Mexico. This vital ecosystem is particularly vulnerable to ongoing climate and ocean change as many of these habitats are rapidly approaching undersaturated carbonate conditions (Lunden et al., 2013; Georgian et al., 2016a), and L. pertusa appears to be living near the edge of its viable niche in the Gulf of Mexico (Lunden et al., 2014a; Georgian et al., 2016b). As cold-water coral habitats are becoming increasingly acidified, the calcification of new skeleton is expected to slow and existing skeletal structures may start dissolving (Andersson et al., 2008). This would result in a loss of ecologically significant primary habitat and regional habitat heterogeneity along continental slopes worldwide (Rogers, 1999; Roberts et al., 2006). In this study, we found an overall trend of decreased calcification rate under exposure to undersaturated conditions, but this response took between 2 weeks and 3 months to manifest in the genotypes examined here. In addition, there was a relatively high degree of variance in the response of individual genotypes across experiments, particularly in the short-term experiment. These findings provide evidence that the existing intra-specific variability in Gulf of Mexico L. pertusa populations may confer a higher degree of resistance to ocean acidification than was previously expected.

Short-Term Calcification Response

Our results did not show a statistically significant response in the short-term (2-week) exposure to low pH in L. pertusa; however, this was primarily due to the extremely high variability in the response of different genotypes within the pH 7.6 treatment. Short-term experiments of ocean acidification with Gulf of Mexico L. pertusa populations have found reductions in net calcification rate with the same pH treatments (Lunden et al., 2014a; Georgian et al., 2016b). Lunden et al. (2014a) found that short-term low pH exposure had an overall negative effect on calcification, although some coral fragments were able to maintain positive calcification rates even at undersaturated levels of aragonite. Moreover, Georgian et al. (2016b) found a 48% decline in net calcification of L. pertusa grown at pH 7.75, and the net dissolution of L. pertusa grown at pH of 7.6. In the low pH treatment (pH = 7.6) in this study, average net dissolution was also observed. However, we observed extremely high variability among genotypes, but consistent responses within genotype. Genotype G4 maintained significantly positive calcification rates in all individuals during the short-term exposure to undersaturated conditions.

Other short-term experiments with cold-water corals including different populations of Lophelia pertusa from the North Atlantic and Mediterranean Sea have shown contrasting results. Despite these inconsistencies, there is evidence for declines in calcification as pH decreases (Maier et al., 2009; Form and Riebesell, 2012). For example, a study using the reef-building cold-water coral Madrepora occulata, found that calcification rates are 50% slower than in pre-industrial times, and suggesting that ocean acidification has already impaired cold-water coral calcification (Maier et al., 2012). Moreover, Maier et al. (2009) found that lowering pH by 0.15 to 0.3 units (from an initial value of pH 8.1) according to future scenarios of ocean acidification resulted in a decrease in L. pertusa's gross calcification by 30 and 56% respectively.

In our study, coral fragments in the control treatment (pH = 7.9) and the slightly acidified treatment (pH = 7.75) maintained similar calcification rates. The similar calcification rates between the control and the 7.75 treatment can likely be attributed to the fact that both treatments were still supersaturated with respect to aragonite (Ω_ar), therefore positive calcification was still energetically favored and dissolution was not likely to occur. These results are similar to the in situ deep-water (McCulloch et al., 2012; Wall et al., 2015) and shallow-water (Cohen and McConnaughey, 2003; Venn et al., 2013; Holcomb et al., 2014) studies that have found corals growing in low pH conditions, presumably due to specialized cellular mechanisms for maintaining elevated pH and carbonate ion concentration within the calcifying fluid.

In the short-term experiment, two of the three replicate fragments from G1 displayed a very high net dissolution rate, while all of the replicates of G4 had positive net calcification within the low pH treatment. It should be noted that the overall net calcification rates reported here for the high pH treatments, designed to simulate ambient in situ conditions, were lower than those found in the previous studies of Gulf of Mexico L. pertusa populations, although within the range of these and other studies (Table 3). Since these studies were conducted in the same laboratory using the same methodology by many of the same scientists, it is possible that some of the coral fragments examined here were in a poorer condition at the time of collection than those collected in previous years, or that they included some genotypes that simply do not perform as well under ambient conditions.

The fact that genotype G4 was able to calcify under low pH and undersaturated conditions suggests that these corals may have started with larger energy reserve, such as carbohydrates, lipids, or proteins, allowing them to allocate additional energy to the calcification process, and as has been found in shallow-water calcifying corals that maintain calcification rates under environmental stress (Rodrigues and Grottoli, 2007; Schoepf et al., 2013). When comparing Gulf of Mexico and Norwegian L. pertusa populations, Georgian et al. (2016b) attributed the higher calcification rates of the Norwegian corals to higher food supply at the shallower Norwegian sites and an elevated metabolic response of the corals. However, increased food supply did not confer an increased resilience to lower pH in Madrepora oculata from the Mediterranean Sea (Maier et al., 2016). It is also possible that certain genotypes are more resilient to low-pH stress via mechanisms not yet understood, but plausibly related to the up-regulation of the calcification pathway or certain stress genes and pathways (Carreiro-Silva et al., 2014) that warrant further investigation.

The finding of a variable response of L. pertusa genotypes to the undersaturated treatment reflects the disparate results from the literature. In the short-term experiment here, two of the genotypes exhibited net dissolution (G1 and G3), while one of the genotypes (G4) was capable of maintaining significant positive net calcification in undersaturated conditions (Ω_ar = 0.81, range 0.76–0.85). This is similar to previous experiments with L. pertusa from the Gulf of Mexico (Lunden et al., 2014a), where a few of the genotypes tested were capable of positive calcification in undersaturated conditions, although the former experiment was not designed to include genetic replicates. Maier et al. (2009) found that L. pertusa from two different populations in the North Atlantic were able to calcify in slightly undersaturated levels of aragonite (Ω_ar = 0.97) albeit at reduced calcification rates. However, the authors used ⁴⁵Ca to measure gross calcification over a 24-h period, which does not take into account possible dissolution effects, while the buoyant-weight technique used here integrates both biologically mediated calcification and chemical dissolution. Therefore, is plausible that all of our coral fragments in the highly acidified treatment (pH = 7.6) were still able to calcify, but at slower rates that are insufficient to compensate for skeletal dissolution.

Long-Term Calcification Response

This is the first long-term experiment performed on L. pertusa from the Gulf of Mexico, although long-term experiments have been conducted using L. pertusa from the North Atlantic (Form and Riebesell, 2012; Hennige et al., 2015), and the Mediterranean sea (Maier et al., 2013; Movilla et al., 2014). Our results showed a significant effect of time and pH, where corals in the control treatment (pH 7.9) maintained positive net calcification, while, on average, corals in the acidified treatment (pH = 7.6) experienced declining net calcification rate over time and net dissolution by the end of the 6-month period. On average, corals were not able to maintain calcification even after only 1 month in the low pH treatments, however there was a variable response among genotypes over time. Genotype 3 maintained positive net calcification for 2 months, and G4 displayed positive net calcification for approximately 3 months before net dissolution started to occur, longer than any of the other genotypes (Figure 4). This provides evidence for an adaptive capacity in this population of L. pertusa that would improve its probability of surviving predicted levels of ocean acidification.

These results contrast with what has been found so far in other long-term experiments where there has not been a significant negative response in calcification for L. pertusa (Form and Riebesell, 2012; Maier et al., 2013; Movilla et al., 2014; Hennige et al., 2015) and for other cold-water coral species (Carreiro-Silva et al., 2014; Rodolfo-Metalpa et al., 2015; Gori et al., 2016). However, the majority of the studies performed in cold-water corals to date have used pH levels that still yielded Ω_ar saturation states at or near 1. Form and Riebesell (2012) found that in the long-term (six months), corals exposed to low pH conditions (pH: 7.75; Ω_ar: 0.93) experienced no significant change in calcification compared to control conditions suggesting an acclimation potential. Maier et al. (2013) and Movilla et al. (2014) in different experiments with L. pertusa and M. oculata, did not find any significant difference in net calcification after 9 and 6 month experiments respectively, however, all of the treatments in those studies were saturated with respect to aragonite (Ω_ar > 1). Likewise, Carreiro-Silva et al. (2014) performed a long-term experiment (8 months) on the cold-water coral Desmophylum diantus subjected to two levels of pH that simulated ambient (pH 8) and acidified conditions (pH 7.7) and found no significant differences in average net calcification, however their experiment were never subjected to undersaturation for aragonite. In this study, none of the L. pertusa genotypes examined were capable of maintaining positive net calcification beyond 90 days in undersaturated conditions. Although we did not measure gross calcification in the present study, it is clear that dissolution is an important factor controlling overall calcification and net calcification rates in the natural environment, and measurements of gross calcification should be incorporated into future studies (Rodolfo-Metalpa et al., 2015).

The most recent direct comparison to the present study, because of the inclusion of a treatment with Ω_ar < 1, measured net calcification in L. pertusa populations from the North Atlantic subjected to undersaturated conditions (Hennige et al., 2015). Although they did not find evidence for a significant response in net calcification rate, their data showed a negative trend with calcification rates approaching zero at undersaturated levels of aragonite. However, it is important to bear in mind that this maintenance of calcification under ocean acidification scenarios comes with a significant energetic cost. Norwegian populations of L. pertusa were found to be able to maintain calcification only by increasing their feeding and respiration rates, while the Gulf of Mexico populations did not exhibit this physiological response (Georgian et al., 2016b).

Population and Genotypic Variability in Response to Ocean Acidification

Environmental differences between the regions may account for some of the variability among experiments conducted in different regions, as relatively isolated L. pertusa populations (Morrison et al., 2011; Lunden et al., 2014a) may have acclimated or adapted to local pH and carbonate conditions. The corals used here were collected from 400 to 500 m at an ambient pH of approximately 7.9 (Table 1), as opposed to the North Atlantic and Norwegian corals that were collected from a depth of 70–150 m and a pH of approximately 8.0 to 8.1 (Maier et al., 2009; Form and Riebesell, 2012). Positive calcification rates were maintained in an experiment on Norwegian corals due to an increase in feeding rate and higher metabolic rate (Georgian et al., 2016b), and these mechanisms will remain effective as long as food supplies are sufficient. Model simulations indicate that L. pertusa reefs induce downwelling of organic-carbon rich waters, which would serve to ensure a consistent food supply (Soetaert et al., 2016), at least while surface productivity rates remain consistent. Even though overall net calcification is positive in lab experiments at low pH and undersaturation of aragonite, dissolution of coral fragments has been noted, particularly in areas that were not covered by tissue (Hennige et al., 2015). There was also a decrease in the breaking strength of new skeletal structure, changes in the crystallographic organization, as well as thinner corallites grown in high CO₂ (Hennige et al., 2015). Furthermore, Gulf of Mexico L. pertusa populations (Georgian et al., 2016b) and Mediterranean M. oculata populations (Maier et al., 2016) do not appear to be able to increase feeding rate to compensate for the increased energetic cost of calcification at low saturation states. These changes have the capacity to alter the stability of reef frameworks in the future.

Standing genetic variation and local adaptation may play a role in species persistence in a changing ocean. Kelly et al. (2013) tested the responses of larval Strongylocentrotus purpuratus cultured under varying pCO₂ conditions, and observed differences in larval length among individuals collected from two different populations characterized by different CO₂ regimes. The authors found that urchins from naturally high-CO₂ environments were less sensitive to OA, implying a capacity for adaptation to ocean acidification. In a related study using several populations of S. purpuratus, Pespeni et al. (2013) observed rapid selection of genes encoding proteins involved in calcification under acidification, demonstrating a potential for rapid evolution to occur in this species. Generally, in slow-growing organisms with long-generation times such as cold-water corals, persistence in a changing ocean will likely hinge on the maintenance of genetic diversity and the success of resilient genotypes, the potential for individual acclimatization, or local adaptation among populations in the presence of both selection and gene flow (Palumbi and Sotka, 2006). While it is unknown if the observed resilience within certain genotypes to climate change stressors in the present study of L. pertusa is heritable, the general consistency in the response of different nubbins of the same genotype within and among different experiments provides evidence for the adaptive capacity of this species.

Identification of the variability of responses due to genetics may lead to future studies of the mechanisms underlying such variation. One mechanism responsible for this could be differential expression of calcification genes. A recent study using next-generation sequencing techniques observed significant effects of OA on gene expression in larvae and primary polyps of the scleractinian coral Acropora millepora (Moya et al., 2012). Notably, genes involved in CaCO₃ deposition showed patterns of differential expression under OA, including down-regulation of genes encoding proteins that support the skeletal organic matrix. It is unclear if this study employed true biological replicates, and different genotypes will likely have different thresholds for their transcriptomic response. An interesting question for future work would be to investigate the potential for different individuals or populations to balance these processes under acidification, as has been observed in natural coral populations responding to thermal stress (Barshis et al., 2013).

Carbonic Anhydrase Activity and Calcification

We did not detect any significant change in the overall activity for carbonic anhydrase in the long-term experiment. Lophelia pertusa maintained stable activity after a 6-month exposure, but exhibited higher variability in the low pH treatment. While our measured values (between 0.002 and 0.031 units of enzyme activity) are low compared to studies of zooxanthellate corals, it has been found that CA activity in zooxanthellate corals is 29 times higher than in azooxanthellate species (Weis et al., 1989). Reported values for CA activity in the azooxanthellate coral Tubastrea aurea were 0.08 units of enzyme activity Weis et al. (1989) reported similar values of 0.06 units for different species of azooxanthellates, which are close to the values obtained here.

In the low pH treatment, there was a slight increase in overall CA activity due to two genotypes that showed increased activity despite a decrease in net calcification. Calcification in scleractinian corals occurs in a biologically mediated environment (Allemand et al., 2011) and carbonic anhydrase actively participate in the calcification process (Tambutté et al., 2007; Moya et al., 2008; Bertucci et al., 2013). In symbiotic scleractinians, reductions in calcification are correlated with decreased enzyme activities (Goreau, 1959; Weis et al., 1989), where calcification is tightly linked with photosynthesis and CA acts as a CO₂ supplier for the endosymbiont (Weis et al., 1989). However, metabolic CO₂ is not the main source of DIC for calcification in deep-water corals (Adkins et al., 2003; Mueller et al., 2013). An increase in CO₂ will cause an increase in DIC, specifically by producing ${HCO}_{3}^{-}$ and H⁺, which ultimately will alter the internal chemistry and lower the pH inside the cells (Allemand et al., 1998). It is plausible that the role of the isoform of αCA found in corals (Bertucci et al., 2013), is to regulate the build-up of H⁺ ions that result from the hydration of CO₂ to HCO³⁻ in order to avoid internal acidosis. This has been suggested for the scleractinians Stylophora pistillata (Moya et al., 2008), and the cold water coral Desmophyllum dianthus (Carreiro-Silva et al., 2014), where an over-expression of the genes that encode for CA was found when high levels of H⁺ were present in the internal compartments of the cells.

Conclusions

These data suggest that further ocean acidification will negatively impact L. pertusa reefs in the Gulf of Mexico, although transient exposures may not have a significant effect, and some genotypes might be more resilient than others. It has been previously stated that ecosystems in the Mediterranean Sea are likely to show the first signs of ocean acidification (Calvo et al., 2011), but our results coupled with previous work (Lunden et al., 2014a; Georgian et al., 2016b) suggest that the deep-water coral ecosystems of the Gulf of Mexico are also extremely vulnerable. Reductions in calcification rate and increased dissolution of standing dead skeleton would cause instability and degradation of the deep reef structures, which would lead to biodiversity loss and cascading effects throughout the larger Gulf of Mexico ecosystem. However, the finding that some genotypes are able to maintain calcification in undersaturated waters for between 2 weeks and 3 months may offer hope for the population as a whole. It is possible that these resistant genotypes will survive the onset of predicted levels of acidification and provide an avenue for adaptive resilience to future climate change. Future work would benefit from examining the genotypic response to acidification and other stressors, with a particular focus on the genetic and physiological mechanisms underlying the resilience of certain genotypes of cold-water corals.

Author Contributions

MK conducted the experiments, performed data analysis, and wrote the manuscript. CG designed and conducted the experiments, performed data analysis, and edited the manuscript. SG collected the corals, designed the experiments, and edited the manuscript. JL conceived of the study, carried out the field work, and edited the manuscript. EC conceived of the study, designed the experiments, led the research cruise, and edited the manuscript.

Conflict of Interest Statement

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.

Acknowledgments

The authors wish to acknowledge Dana Russell and Adam Butler who helped in the set-up and tank maintenance, and Alex Barkman who helped in running the carbonic anhydrase assays. Thanks to Cheryl Morrison for conducting the microsattelite analysis of the coral collections. The authors are grateful for the support and efforts of the pilots of the DSV Alvin and the crew and officers of the R/V Atlantis who made this work possible. We also wish to acknowledge all of the support of the scientists on cruise AT26-14. Special thanks to Ivan Hürzeler and the rest of the Acid Horizon team for helping to spread the word. This research was supported by NSF BIO-OCE grant #1220478. CEG was supported by the Fulbright-Colciencias scholarship program.

References

Adkins, J. F., Boyle, E. A., Curry, W. B., and Lutringer, A. (2003). Stable isotopes in deep-sea corals and a new mechanism for “vital effects.” Geochim. Cosmochim. Acta 67, 1129–1143. doi: 10.1016/S0016-7037(02)01203-6