Physiological and Biogeochemical Responses of Super-Corals to Thermal Stress from the Northern Gulf of Aqaba, Red Sea

Grottoli, Andréa G.; Tchernov, Dan; Winters, Gidon

doi:10.3389/fmars.2017.00215

ORIGINAL RESEARCH article

Front. Mar. Sci., 04 July 2017

Sec. Coral Reef Research

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

This article is part of the Research Topic The Future of Coral Reefs Subject to Rapid Climate Change: Lessons from Natural Extreme Environments View all 16 articles

Physiological and Biogeochemical Responses of Super-Corals to Thermal Stress from the Northern Gulf of Aqaba, Red Sea

$\r\nAndra G. Grottoli*$ Andréa G. Grottoli¹^*

Dan Tchernov²

Gidon Winters³

¹School of Earth Sciences, Ohio State University, Columbus, OH, United States
²Leon H Charney School of Marine Science, University of Haifa, Haifa, Israel
³The Dead-Sea and Arava Science Center (ADSSC), Tamar Regional Council, Neve Zohar, Israel

Mass coral bleaching is increasing in frequency and severity, leading to the loss of coral abundance and diversity. However, some corals are less susceptible to bleaching than others and can provide a model for identifying the physiological and biogeochemical traits that underlie coral resilience to thermal stress. Corals from Eilat in the Gulf of Aqaba in the northern Red Sea do not bleach unless seawater temperatures are sustained at +6°C or higher above their average summer maximum. This extreme thermal tolerance qualifies these as super-corals, as most corals bleach when exposed to temperatures that are only +1–2°C above their thermal maximum. Here, we conducted a controlled bleaching experiment (+6°C) for 37 days (equivalent to 32° heating weeks) on three species of corals from Eilat: Stylophora pistillata, Pocillopora damicornis, and Favia favus. To assess the response of the holobiont to thermal stress, the following variables were measured on each coral: endosymbiotic algal cell density, Chlorophyll a, endosymbiotic mitotic cell division, total lipids, protein, carbohydrate, and the stable carbon (δ¹³C) and oxygen (δ¹⁸O) isotopic composition of the skeleton and the δ¹³C of the animal host tissue and endosymbiotic algae. While all three species appeared visibly bleached, their physiological and biogeochemical responses were species-specific. S. pistillata catabolized lipids but still maintained total energy reserves and biomass. Increases in both skeletal δ¹³C and δ¹⁸O indicates that calcification declined in this species. P. damicornis was the least affected by bleaching. It maintained its total energy reserves and biomass, and isotopic evidence suggests that it maintained calcification and was not dependent on heterotrophy for meeting metabolic demand when bleached. Finally, F. favus catabolized protein and carbohydrates, and suffered losses in total energy reserves and biomass. Nevertheless, isotopic evidence suggest that photosynthesis and calcification were maintained, and that this species has a high baseline heterotrophic capacity. Thus, just like their non-super-coral conspecific counterparts, maintaining energy reserves and biomass, and heterotrophic capacity appear to be traits that underlie the thermal tolerance of these super-corals from Eilat. Given the high thermal tolerance of these super-corals, these populations could provide viable seed stock for repopulating coral losses on other reefs.

Introduction

Increased atmospheric carbon dioxide (CO₂) is causing the oceans to warm and to become more acidic. Over the course of this century, tropical seawater temperatures are expected to rise by 1–3°C coupled with a drop in pH of 0.2–0.4 (IPCC, 2013). At this rate, tropical reefs are expected to experience significant declines in coral abundance, coral diversity, and reef growth during the second half of this century (e.g., Hoegh-Guldberg et al., 2007; Eakin et al., 2009; Veron et al., 2009; Hoegh-Guldberg, 2011). While coral calcification is projected to decline by 5% by 2034, severe bleaching will be annual by 2030 in some regions and globally by 2055 (van Hooidonk et al., 2014). Recent studies have shown that the warming of tropical oceans is a much more immediate threat to the survival of coral reefs than is ocean acidification (e.g., Cantin et al., 2010; Frieler et al., 2012; Chua et al., 2013). Corals can survive slight decreases in calcification, but it is not likely that they can physiologically survive severe bleaching on an annual basis.

Elevated temperatures can cause corals to bleach: a physiological response to sustained elevated temperatures that results in a breakdown of the coral symbiosis, the loss of the coral's algal endosymbiont partners, which results in a white appearance. The consequences of bleaching to the coral are reduced coral growth, disease, reduced fertility, and death (e.g., Brown, 1997; Hoegh-Guldberg, 1999; Omori et al., 1999; Buddemeier et al., 2004; Tchernov et al., 2011). The frequency of bleaching events appears to be increasing (Eakin et al., 2009; Heron et al., 2016) and back-to-back bleaching years have been observed on basin-wide scales (coralreefwatch.noaa.gov). Annual bleaching is predicted to occur by 2030 in some regions (van Hooidonk et al., 2013, 2014), resulting in large coral losses.

Interestingly, there are corals that are less susceptible to and/or recover more quickly from bleaching (e.g., Loya et al., 2001; Grottoli et al., 2006, 2014; Fine et al., 2013). Resilience to bleaching has been associated with high energy reserves, heterotrophic plasticity, and endosymbiont type (Rowan et al., 1997; Baker et al., 2004; Tchernov et al., 2004; Berkelmans and van Oppen, 2006; Grottoli et al., 2006, 2014; Rodrigues and Grottoli, 2007; Anthony et al., 2009; Guest et al., 2012). Corals with elevated levels of energy reserves (i.e., lipids, protein, carbohydrates) either take longer to bleach, bleach less severely, and/or recover more quickly from bleaching than corals with lower energy reserves because they have an energy source to buffer them against losses in photosynthetically derived fixed carbon (e.g., Rodrigues and Grottoli, 2007; Anthony et al., 2009; Grottoli et al., 2014; Schoepf et al., 2015). Corals that increase their intake of heterotrophic carbon (i.e., zooplankton, dissolved and particulate organic carbon) or increase the proportionate contribution of heterotrophic vs. photoautotrophic carbon in their tissues when bleached, are able to partially or fully supplement the deficit in their carbon budgets due to declines in photosynthesis when bleached, and are able to recover more quickly from bleaching (Grottoli et al., 2006, 2014; Rodrigues and Grottoli, 2006; Palardy et al., 2008; Levas et al., 2013, 2016). Corals that are able to shuffle or switch their dominant endosymbiont type for a thermally tolerant one are less sensitive to repeat exposures to thermal stress (e.g., Rowan et al., 1997; Baker et al., 2004; Berkelmans and van Oppen, 2006; Grottoli et al., 2014). In addition, there are regions such as the northern region of the Gulf of Aqaba, where entire communities of corals can already tolerate extreme warming conditions, thought the underlying traits that produce that extreme tolerance is poorly understood. These corals may provide seed stock for restoring dying reefs as well as be a model system for identifying the physiological and biogeochemical traits that underlie coral resilience to thermal stress.

Corals in the northern tip of the Gulf of Aqaba in the Red Sea do not bleach when exposed to temperatures that are ≥6°C above their average summer time maximum temperatures for 3–4 weeks (Fine et al., 2013) coupled with extreme salinities of 41ppt (Paldor and Anati, 1979). This is remarkable resilience to extreme increases in seawater temperature given that most corals bleach when seawater temperatures exceed summertime maximum temperature by as little as 1–2°C for as little as 2–4 weeks (e.g., Hoegh-Guldberg, 2011). Fine et al. (2013) proposed that this extreme thermal tolerance in these northern Red Sea corals is due to natural selection processes on the larvae transported to the northern section of the Gulf of Aqaba through the very hot straights of Bab el Mandeb in the southern Red Sea. However, the underlying physiological traits responsible for the resilience of these corals is not fully understood. We hypothesized that like their less thermally tolerant counterparts, maintenance of energy reserves and biomass, as well as heterotrophic plasticity, are all traits responsible for extreme temperature tolerance in super-corals of the northern Gulf of Aqaba. To test this hypothesis, bleaching was induced by thermal stress in an experiment conducted on three species of corals (Stylophora pistillata, Pocillopora damicornis, and Favia favus) from the Gulf of Aqaba near Eilat, Israel and the physiological and biogeochemical responses of the host and algal endosymbiont were assessed.

Materials and Methods

Two fragments from five healthy colonies of S. pistillata and P. damicornis (5–7 cm tall), and two fragments for seven healthy colonies of F. favus (6–8 cm in diameter) were collected between 2 m and 4 m depth in the Gulf of Eilat (northern Gulf of Aqaba), Israel (29° 30′ 5.55″ N, 34° 55′ 2.49″ E) in mid-February 2005 and glued to labeled gray PVC tiles using underwater epoxy putty. All three species are common in this region. S. pistillata is a branching coral, a reproductive brooder, and is dominated by the endosymbiont Symbiodinium clade A (Shlesinger and Loya, 1985; Karako-Lampert et al., 2004; Winters et al., 2009; Byler et al., 2013; Ziegler et al., 2017). P. damicornis is also a branching coral, is a member of a genus that reproduces by both broadcasting and brooding (though spawning has never been documented for this species in the Red Sea), and is dominated by the endosymbiont Symbiodinium clade C (Shlesinger and Loya, 1985; Karako-Lampert et al., 2004; Fine et al., 2013; Sawall et al., 2015; Ziegler et al., 2017). F. favus is a mounding coral, a broadcast spawner, and dominated by the endosymbiont Symbiodinium clade C (Shlesinger and Loya, 1985; Karako-Lampert et al., 2004; Fine et al., 2013). Symbiodinium clade D has never been reported for healthy or bleached corals in this region (Fine et al., 2013).

One fragment from each colony was placed in each of two 95-liter outdoor flow-through seawater tanks at the Interuniversity Institute for Marine Science of Eilat, Israel at the northern tip of the Gulf of Aqaba, Red Sea, and allowed to acclimate to the tanks and recover for 2 weeks. The tanks were shaded with black shading nets to simulate photosynthetically active radiation levels at collection depth, and the seawater was filtered through a small coarse gravel filter (70–200 μm pore size), allowing only small organic particulates and dissolved organics into the tanks. Thus, heterotrophic food sources available to the corals were limited to dissolved and particulate organic matter that came through the tank system. Coral were not supplemented with zooplankton. Flow rates were maintained in each tank at 2 l/min allowing each tank to fully flush every 47 min. On 28 February 2005, the seawater temperature was gradually increased over 3days using three 300W ATMAN aquarium heaters per tank yielding an average daily temperature of 32°C, while the second tank remained at ambient seawater temperature as controls with an average temperature of 26°C. The rate of temperature increase in the 32°C tank was chosen to minimize heat shock and also get the tank up to temperature relatively quickly. Since these corals are very tolerant of elevated temperatures, the rate of increase was not detrimental to these corals. The original temperature logger data has been lost, but the typical standard deviation of average daily temperatures in experimental tanks at the institute are ± 0.5°C. On 6 April 2005 (37 days later), the fragments were collected, frozen at −80°C, and shipped frozen to the Grottoli laboratory for analyses. In the lab, a quarter of each fragment was used for the following analyses: (1) chlorophyll a (Chl a) concentration, (2) total lipid, soluble animal protein, and soluble animal carbohydrate concentrations, (3) endosymbiotic algae concentrations and mitotic indexes, and (4) stable isotopic analyses of the host tissue, endosymbiotic algae, and underlying skeleton.

Algal Endosymbiont Analyses and Coral Energy Reserves

Chl a was extracted from ground samples in 100% acetone and the absorbance measured at 663 and 630 nm wavelengths using a Spectronic 20 Genesis visible spectrophotometer. The Chl a concentration was calculated using equations from Jeffrey and Humphrey (1975) corrected for turbidity with the absorbance value measured at 750 nm. Algal endosymbionts were separated from the animal tissue by centrifugation and three subsamples were counted with a hemocytometer and light microscope. The total number of algal cells and the number of mitotically dividing algal cells were recorded.

Total lipids were extracted in a 2:1 chloroform/methanol solution, washed in 0.88% KCl, and then dried to a constant weight (method modified from Grottoli et al., 2004). Soluble animal protein was extracted using the bicinchoninic acid method (Smith et al., 1985) with bovine serum albumin as a standard (Pierce BCA Protein Assay Kit). Soluble animal carbohydrate was extracted using the phenol-sulfuric acid method (Dubois et al., 1956) with glucose as a standard. All analyses were made from whole ground coral samples, converted into Joules (Gnaiger and Bitterlich, 1984), and standardized to both grams of ash free dry weight (gdw) and to surface area as measured using the aluminum foil technique (Marsh, 1970). However, normalization to ash free dry weight is more robust when comparing species with different tissue depths, morphology, and polyp structure than the same data normalized to surface area (Edmunds and Gates, 2002). Therefore, only results of data standardized to ash free dry weight are reported here (However, note that all statistics were also performed on the area normalized data and yielded almost identical patterns). Total energy reserves were calculated as the sum of total lipids, proteins, and carbohydrates. Total dry biomass was standardized to surface area. These analytical methods have been successfully used in corals before (e.g., Rodrigues and Grottoli, 2007).

Stable Isotopic Analyses

Coral tissue was completely removed from the skeleton using a Water-pik (Johannes and Wiebe, 1970). The tissue slurry was separated into animal host and algal endosymbiont fractions by centrifugation using established methods (Rodrigues and Grottoli, 2006). Briefly, coral slurry was homogenized to break open the animal cells, then centrifuged at 4,000 rpm for 5 min to separate the animal host and endosymbiotic algal fraction. The supernatant (i.e., animal host) was decanted and the algal pellet resuspended in deionized water and centrifuged again to remove any remaining animal host particles. The supernatants from both centrifugations were combined. Each organic fraction was isolated onto a pre-burned glass fiber filter under vacuum, and individually combusted in a Costech Elemental Analyzer. The resulting CO₂ gas was analyzed for δ¹³C [δ¹³C = permil deviation of ¹³C:¹²C relative to Vienna Peedee Belemnite Limestone standard (VPDB)] with a Thermo Delta Advantage IV via a Thermo ConFlow III open split in the Grottoli lab. At least 10% of all samples were run in duplicate. The standard deviation of repeated measurements of the USGS24 standard for δ¹³C was ±0.08 ‰. The host tissue isotope values were reported as δ¹³C_h while the algal endosymbiont isotope data were reported as δ¹³C_e. In general, as the incorporation of isotopically depleted zooplankton and other heterotrophic carbon sources (i.e., dissolved and particulate organic carbon) into coral tissues increases, the δ¹³C_h decreases (Rodrigues and Grottoli, 2006; Levas et al., 2013; Schoepf et al., 2015). As the rate of photosynthesis increases, isotopic fractionation decreases, and endosymbiotic algae incorporate relatively less ¹²C than ¹³C into their tissues, resulting in increased δ¹³C_e (Muscatine et al., 1989; Rodrigues and Grottoli, 2006). The difference between the δ¹³C_h and δ¹³C_e were also calculated to assess the proportionate contribution of heterotrophy and photoautotrophy in the carbon budget (sensu Muscatine et al., 1989). Specifically, decreases in δ¹³C_h–δ¹³C_e are interpreted to mean that the proportionate contribution of heterotrophic carbon to coral tissues has increased.

The skeleton was prepared for δ¹³C (δ¹³C_s) and δ¹⁸O (δ¹⁸O_s)[δ¹⁸O = permil deviation of ¹⁸O:¹⁶O relative to VPDB] analysis by shaving the top 100–200 μm of skeletal material off the tip of each coral fragment using a Dremel tool with a diamond-tipped drill bit. Approximately ~80 μg of skeletal material was acidified with 100% ortho-phosphoric acid in an automated Kiel III carbonate device, and the δ¹³C and δ¹⁸O values of the resulting CO₂ were measured using a Finnigan MAT 252 mass spectrophotometer at Stanford University. At least 10% of all samples were run in duplicate. The standard deviation of repeated measurements of an internal standard was ±0.03 ‰ for δ¹³C and ± 0.05 ‰ for δ¹⁸O. Skeletal values were reported as δ¹³C_s and δ¹⁸O_s. In general, as the rate of photosynthesis increases or as the incorporation of heterotrophic carbon into coral tissues decreases, the pool of inorganic C available for calcification through animal respiration becomes enriched resulting in increased δ¹³C_s (Grottoli and Wellington, 1999; Grottoli, 2002). In addition, variability in δ¹⁸O_s results from temperature-induced kinetic fractionation and seawater δ¹⁸O (Epstein et al., 1953). As temperature increases, δ¹⁸O_s decreases (Weber and Woodhead, 1972; Kim and O'Neil, 1997). Since seawater δ¹⁸O_sis constant in this study because the source seawater was the same in all tanks, it did not influence the δ¹⁸O_s values.

Statistical Analysis

Univariate two-way analysis of variance (ANOVA) was used to test the effects of species and temperature on Chl a, algal endosymbiont, mitotic index, lipid, protein, carbohydrate, total energy reserves, and tissue biomass concentrations as well as δ¹³C_s, δ¹⁸O_s, δ¹³C_h, δ¹³C_e, and δ¹³C_h–δ¹³C_e. Prior to statistical analysis, all data was tested for normality using a Shapiro-Wilk's test and homogeneity of variance was assessed with plots of expected vs. residual values. Any data failing to meet this assumption was log transformed and then met the assumptions of normal distribution. A posteriori slice tests (i.e., tests of simple effects, Winer, 1971) determined if bleached and control averages significantly differed within species. Bonferroni corrections were not used due to increased likelihood of false negatives (Quinn and Keough, 2002). With only one tank for treatment and one tank for controls, we were unable to assess tank effects independent of treatment effects. However, every aspect of both tanks that were measured were the same (i.e., seawater source, flow, light, position relative to the sun's arc) except for temperature. In addition, one bleached P. damicornis and two control F. favus fragments were lost during shipping or storage, yielding an unbalanced number of samples in the treatment and controls for these two species. Statistical analyses were generated using SAS software, Version 8.02 of the SAS System for Windows. Values of p < 0.05 were considered significant.

To get an overall assessment of the effect of treatment across all variables, non-parametric multi-variate statistics were employed. The non-parametric analyses were conducted on three sub-sets of the data: (1) all of the physiology data, (2) all of the isotopic data, and (3) all of the physiology and isotopic data combined. Since two samples did not have isotopic data, the multi-variate analyses that included the isotopes were comprised of 29 corals whereas the analyses of the physiological variables alone included 31 corals. Euclidean-distance based resemblance matrix was constructed using normalized data of all of the measured variables. Non-parametric multidimensional scaling (NMDS) plots were generated to graphically represent relationships between each coral sample across all treatments and species. Analysis of similarities (ANOSIM) was then used to test for the effect of species and treatment. ANOSIM simultaneously uses all variables measured, and tests whether all data (i.e., all physiological and/or isotopic) from one coral group differs from another coral group. The ANOSIM pairwise test statistic R ranges from 0 (no difference and complete overlap between groups) to 1 (maximum difference and no overlap between groups) and is a strong indicator of separation between groups, especially at low sample sizes such as in this study, as long as the overall ANOSIM test statistic R is significant (Clarke and Gorley, 2006). Multi-variate analyses were conducted using the software package Primer V6, and for the ANOSIMs a p < 0.05 was considered significant. Finally, since fragments pairs of each colony were included in both treatments, differences between treatment and control for any variable were due to treatment alone, independent of genotype.

Results

The temperature exposure of the treatment corals was equivalent to 32° heating weeks (DHW) (coralreefwatch.noaa.gov) relative to the controls. All three species in the elevated temperature treatment (i.e., 32°C) appeared pale or white and as such were bleached, whereas all of the corals in the control treatment (i.e., 26°C) were brown in color and appeared non-bleached.

Endosymbiotic Algae

Average endosymbiotic algal concentrations dropped to 5, 49, and 19% of control concentrations in S. pistillata, P. damicornis, and F. favus, respectively (Table 1, Figure 1A). At the same time, Chl a concentrations in bleached S. pistillata and P. damicornis decreased to 37 and 33% of control values, respectively (Table 1, Figure 1B), and mitotic divisions were 26-fold lower in bleached S. pistillata and 9-fold lower in bleached F. favus corals relative to their controls (Table 1, Figure 1C). While endolithic algae were not present in noticeable quantities in either P. damicornis or S. pistillata, F. favus showed evidence of harboring endolithic algae. Because all samples of F. favus showed approximately the same quantity of endolithic algae irrespective of treatment, it was assumed that their contribution to Chl a was relatively constant and that any differences in Chl a between the bleached and control fragments was due to the bleaching effect alone.

TABLE 1

Table 1. Results of three two-way ANOVAs testing for the effect of species (S), and temperature (T) on average chlorophyll a concentration, endosymbiotic algal density and mitotic index.

FIGURE 1

Figure 1. Average (A) endosymbiotic algal density, (B) chlorophyll a (Chl a), and (C) mitotic division density in control (black bars) and bleached (gray bars) Stylophora pistillata, Pocillopora damicornis, and Favia favus corals. Averages are plotted ± 1 standard error, sample sizes for all panels indicated in brackets in (A), and significant differences between average control and bleached values within each species are indicated with a *. Corresponding ANOVA results are given in Table 1. gdw, grams of ash free dry weight.

Energy Reserves

While lipids and protein were unaffected by bleaching overall, a posteriori slice tests revealed that lipid concentrations in bleached S. pistillata and protein concentrations in bleached F. favus were 49 and 48% lower, respectively, than in their controls (Table 2, Figures 2A,B). Carbohydrate concentrations dropped to 40 and 41% of control values in bleached P. damicornis and F. favus, respectively (Table 2, Figure 2C). Despite losses of at least one component of energy reserves in each species, only bleached F. favus had significantly lower total energy reserves (43% decrease) and experienced a 27% decrease in biomass compared to controls at the species level (Table 2, Figures 2D,E).

TABLE 2

Table 2. Results of five two-way ANOVAs testing for the effect of species (S), and temperature (T) on average lipid, protein, carbohydrate, total energy reserves, and total biomass.

FIGURE 2

Figure 2. Average (A) Lipid, (B) protein, (C) carbohydrate, (D) total energy reserves, and (E) total biomass in control (black bars) and bleached (gray bars) Stylophora pistillata, Pocillopora damicornis, and Favia favus corals. Averages are plotted ± 1 standard error, sample sizes for all panels indicated in brackets in (A), and significant differences between average control and bleached values within each species are indicated with a *. Corresponding ANOVA results are given in Table 2. J, Joules; Total Energy Res, total energy reserves which is the sum of lipids, protein, and carbohydrate values.

Stable Isotopes

Average δ¹³C_s values were significantly higher in bleached than in control corals overall (Table 3), but this difference was only significant at the species level for S. pistillata where bleached δ¹³C_s was 1.7 ‰ higher than control corals (Table 3, Figure 3A). Average δ¹⁸O_s was 0.9 ‰ higher in the bleached than in control S. pistillata (Table 3, Figure 3B). At the same time, average δ¹³C_h values were 2.2 ‰ and 1.8 ‰ higher in bleached than in control S. pistillata and F. favus (Table 3, Figure 3C), respectively, while average δ¹³C_e values were 1.0 ‰ higher in bleached than in control F. favus, and at least 2.4 ‰ higher in F. favus overall compared to the other two species (Table 3, Figure 3D). Finally, average δ¹³C_h–δ¹³C_e values were 0.8 ‰ higher in bleached than in control corals overall (Table 3), but this difference was only significant at the species level for S. pistillata where the average δ¹³C_h–δ¹³C_e was 1.5 ‰ higher in bleached than in control corals (Figure 4).

TABLE 3

Table 3. Results of five two-way ANOVAs testing for the effect of species (S), and temperature (T) on coral skeletal δ¹³C, skeletal δ¹⁸O, host tissue δ¹³C, endosymbiontic algal δ¹³C, and the difference between host and endosymbiontic algal δ¹³C (δ¹³C_h–δ¹³C_z).

FIGURE 3

Figure 3. Average stable isotopes in control (black circles) and bleached (white circles) coral skeleton (s), host tissue (h) and algal endosymbiont (e) fractions. (A) δ¹³C_s, (B) δ¹⁸O_s, (C) δ¹³C_h, and (D) δ¹³C_e in Stylophora pistillata, Pocillopora damicornis, and Favia favus corals. All data are reported in permil (‰) relative to V-PDB. Averages are plotted ± 1 standard error, sample sizes for each average are indicated in brackets, and significant differences between average control and bleached values within each species are indicated with a *. Corresponding ANOVA results are given in Table 3.

FIGURE 4

Figure 4. Average difference between host (δ¹³C_h) and endosymbiont (δ¹³C_e) stable carbon isotope values in bleached (open circles) and control (black circles) Stylophora pistillata, Pocillopora damicornis, and Favia favus corals from Figures 3C,D. Heterotrophy contributes more to the fixed carbon pool when the difference is <0, while photosynthesis contributes more when the difference is ≥0. Averages are plotted ±1 standard error calculated through error propagation and are reported in permil (‰) relative to V-PDB. Sample sizes for each average are in brackets. Significant differences between control and bleached within each species are indicated with *. Corresponding ANOVA results are given in Table 3.

Multi-Variate Assessments

Overall, both bleached S. pistillata and F. favus corals significantly differed from their controls, whereas bleached and control P. damicornis corals did not differ from each other (Table 4, Figure 5). This was true for all combinations of the data (i.e., physiological variables alone, isotopic variables alone, of both physiological and isotopic data combined). The degree of separation between the bleached and control corals of both S. pistillata and F. favus was similar with R values produced from the physiological and isotopic data combined of 0.66 and 0.64, respectively. In addition, all species significantly differed from each other, though F. favus differed the most from the other two species as indicated by R > 0.7 when compared to either of the other two species when both the physiological and isotopic data were considered (Table 5, Figure 5).

TABLE 4

Table 4. One-way ANOSIMs of physiological, isotopic, and both physiological and isotopic variables combined where each individual species and treatment was a factor allowing for pairwise tests of each species and treatment combination.

FIGURE 5

Figure 5. Non-parametric multi-dimensional scaling ordination using the physiological and isotopic variables from 29 of the 31 corals. Two corals were not included because they did not have complete isotope data sets. The NMDS plots and corresponding ANOSIMS patterns were similar when build from just the physiology data (all 31 corals) or just the isotopic data (29 corals). Corresponding ANOSIMs results are given in Tables 4, 5.

TABLE 5

Table 5. Two-way ANOSIMs of physiological, isotopic, and both physiological and isotopic variables combined where species (S. pistillata, P. damicornis, F. favus) and treatment (26°C, 32°C) were the factors.

Discussion

Overall, all treatment corals bleached (Figure 1), but the degree to which they were physiologically affected by the temperature stress varied among species (Figures 1–5). Here, we discuss the physiological and isotopic responses of each coral species, and explore the implications of those responses to the persistence of corals in the future.

Stylophora pistillata

S. pistillata bleached the most severely of all three species studied with the greatest percentage drop in endosymbiotic algal density (Figure 1A). Large drops in endosymbiotic algal density are associated with dramatic decreases in photosynthesis, resulting in corals that are unable to meet their total metabolic demand (Grottoli et al., 2006, 2014; Tremblay et al., 2012). Here, S. pistillata catabolized its lipid reserves (Figure 2A), presumably to compensate for losses in photosynthetically fixed carbon. The decline in lipids could have been further exacerbated by the absence of supplemental feeding, since heterotrophic carbon is critical to lipid synthesis in bleached corals (Baumann et al., 2014). Increases in δ¹³C_h (Figure 3C) suggest that the lipids catabolized during bleaching (Figure 2A) were isotopically depleted and/or that in the absence of isotopically depleted zooplankton in the diet, coral tissues became progressively enriched. Given that bleached Hawaiian corals were found to catabolize isotopically enriched lipids (Grottoli and Rodrigues, 2011), it is most likely that the enriched δ¹³C_h observed here in bleached corals is a result of a lack of zooplankton. In addition, increases in δ¹³C_h−e (Figure 4) suggests that S. pistillata incorporated less heterotrophic C into its tissues when bleached [i.e., increases (decreases) in δ¹³C_h−e are indicative of dramatic decreases (increases) in the incorporation of heterotrophic carbon in tissues (Rodrigues and Grottoli, 2006)]. This is consistent with the fact that corals in this study were not fed zooplankton during the experiment. However, even if they had been fed, this species is known to decrease feeding rates on zooplankton, pico-, and nano-plankton, as well as suffer increases in total organic carbon losses when bleached (Ferrier-Pages et al., 2010; Tremblay et al., 2012), and that it has limited heterotrophic plasticity capacity in response to a changing environment (Alamaru et al., 2009). Despite the loss of lipids and isotopic evidence of reduced heterotrophic carbon input, bleached S. pistillata only showed a trend of declining total energy reserves (p < 0.13) and maintained its biomass (Figures 2D,E). This may be because this species is able to rapidly curb apoptosis and acclimate to sustained thermal stress when exposed to elevated temperatures (Kvitt et al., 2014, 2016).

The increase in δ¹³C_s and δ¹⁸O_s when bleached (Figures 3A,B) is consistent with decreases in calcification during bleaching (Rodrigues and Grottoli, 2006) due to greater equilibration with seawater DIC when calcification rates drop (McConnaughey, 1989). Increases in δ¹³C_s could also indicate an increase in photosynthesis rates (Grottoli and Wellington, 1999; Grottoli, 2002). However, this is not likely since this species typically experiences significant declines in photochemical efficiency (Fv/Fm) and photosynthesis rates at elevated temperatures (Tremblay et al., 2012; Fine et al., 2013), and because it is not physiologically possible to increase photosynthesis rates when bleached. In addition, the increase in δ¹⁸O_s deceptively suggests that seawaters were cooler, not hotter, in the bleached coral tank (Figure 3B). Increases in δ¹⁸O_s in bleached corals are caused by isotopic equilibration with seawater as calcification dramatically slows (Rodrigues and Grottoli, 2006). Thus, the skeletal isotopic signature here indicates a decline in calcification in bleached S. pistillata. Decreasing calcification rates is one mechanism by which S. pistillata may reduce its energetic demand and maintain total energy reserves and biomass. In principle, maintenance of energy reserves and biomass should confer an advantage to S. pistillata and facilitate recovery from bleaching, as has been found for the branching Caribbean coral Porites divaricata (Grottoli et al., 2014).

Pocillopora damicornis

Of all three species studied here, P. damicornis was the least affected by bleaching with the smallest decline in endosymbiont cell density and the only species to maintain mitotic cell division (Figure 1). Fine et al. (2013) showed that photochemical efficiency of PSII in this species is among the least affected by 34°C temperatures compared to other corals from this site, despite hosting the typically thermally sensitive clade C Symbiodinium (Karako-Lampert et al., 2004). At the same time, P. damicornis preserved its energy reserves and biomass (Figure 2), showed no metabolic shifts in its isotopic signatures (Figures 3, 4), and overall was unchanged when bleached (Figure 5). The maintenance of its skeletal isotopic values when bleached (Figures 3A,B) suggests that calcification was maintained when bleached (Rodrigues and Grottoli, 2006). This is not an unreasonable interpretation as some populations of Pacific P. damicornis corals have optimal calcification rates at 31°C (Clausen and Roth, 1975)—only 1° cooler than the experimental temperatures. The lack of any shifts in δ¹³C_h, δ¹³C_e, and δ¹³C_h−e indicate that the proportionate contribution of photoautotrophic and heterotrophic C to tissues did not change with bleaching. Given that these corals were not fed zooplankton during the experiment, these results further suggest that P. damicornis may not be very dependent on heterotrophic C for meeting metabolic demand or may not be able to change its heterotrophic intake under bleaching stress. Work by Ziegler et al. (2014) shows that Pocillopora are not heterotrophically plastic in response to changes in light, which may indicate that they would not be heterotrophically plastic in response to temperature changes either.

Schoepf et al. (2013) also found no changes in energy reserves and biomass in P. damicornis under control and elevated temperature conditions (26.5 vs. 29°C). This is surprising since corals of this genus are often among the most sensitive to thermal stress in other regions (Hueerkamp et al., 2001; Obura, 2001; McClanahan et al., 2004; Guest et al., 2012; Foster et al., 2014; Pisapia et al., 2016) and even suffer lipid losses when bleached (Rodriguez-Troncoso et al., 2016). But there are also cases where Pocillopora corals do not readily bleach compared to other species (Sebastian et al., 2009; Guest et al., 2012). Some evidence suggests that Pocillopora corals are beginning to acclimatize or adapt to thermal stress in Australia (Ulstrup et al., 2006; Maynard et al., 2008) and may have high acclimatization potential in northern Red Sea populations (Sawall et al., 2015), though not through heterotrophic plasticity mechanisms (Ziegler et al., 2014). These findings suggest that P. damicornis in Eilat is a more resilient population to temperature-induced bleaching than those in other locations because (1) they can maintain their energy reserves, biomass, and calcification, and (2) appear to not be dependent on heterotrophy for meeting metabolic demand when bleached.

Favia favus

Bleached F. favus suffered declines in more physiological variables than either of the other two species with significant decreases in endosymbiotic algal density, mitotic divisions, protein, carbohydrates, total energy reserves, and biomass (Figures 1, 2), yet still maintained Chl a concentrations (Figure 1B). Even though photochemical efficiency (F_v/F_m) tends to decline at elevated temperatures for F. favus from Eilat (Winters et al., 2006; Kuguru et al., 2007; Fine et al., 2013), the observed increase in δ¹³C_e (Figure 3D) and maintained Chl a levels suggests that some level of photosynthesis persisted while bleached. At the same time, F. favus catabolized its protein and carbohydrate reserves, which underlie the significant declines in total energy reserves and biomass (Figures 2B–E), presumably to compensate for losses in photosynthetically fixed C. Increases in δ¹³C_h (Figure 3C) suggest that the protein and carbohydrates catabolized during bleaching (Figures 2B,C) were isotopically depleted, or that in the absence of isotopically depleted zooplankton in the diet, coral tissues became progressively enriched. Further research is needed to evaluate the isotopic composition of protein and carbohydrates. However, since these corals were not supplied zooplankton during the experiment, it is likely that the enrichment in δ¹³C_h observed here is at least in part a result of a lack of zooplankton availability. Although the δ¹³C_h−e did not change with bleaching in F. favus (Figure 4), it was lower overall than in S. pistillata and P. damicornis indicating a higher contribution of heterotrophic C to its tissues in general, irrespective of bleaching, than in the other two species (Figure 4). This is consistent with field-based findings showing that F. favus allocates a higher proportion of heterotrophic C to its lipids than does S. pistillata (Alamaru et al., 2009). In addition, some species are able to meet up to 36% of metabolic demand when bleached through the heterotrophic acquisition of dissolved and particulate organic carbon (Levas et al., 2013, 2016). The lower overall δ¹³C_h−e of F. favus suggests that when bleached it was able to take up dissolved and/or particulate organic carbon during the tank experiment to meet at least part of its metabolic demand heterotrophically. But any uptake in heterotrophic C was not enough to prevent large losses in protein and biomass. Had the corals been fed zooplankton, it is possible that F. favus might have maintained its energy reserves and biomass.

Finally, the skeletal isotopic signature in F. favus was unaffected by bleaching (Figures 3A,B), suggesting that calcification was maintained or at least did not decrease dramatically (Rodrigues and Grottoli, 2006). This is also consistent with the interpretation that some level of photosynthesis was maintained, since calcification is tightly linked to photosynthesis (see review by Allemand et al., 2011). Furthermore, in healthy corals, δ¹⁸O_s decreases by ~0.2 ‰ per 1°C seawater warming (Weber and Woodhead, 1972; Wellington et al., 1996). Here a non-significant trend of a 0.36 ‰ decline in δ¹⁸O_s (P < 0.085) suggests that F. favus skeletal isotopic signature recorded the seawater warming event, though not the ~1.2 ‰ decline that would have be expected if it had accurately recorded the full 6°C increase in temperature. Nevertheless, the fact that a warming trend was recorded is further support that F. favus calcified throughout the experiment.

Overall, F. favus suffered significant declines in energy reserves and biomass, but maintained Chl a, appears to have maintained calcification, and appears to have high baseline rates of heterotrophic C contributions to its tissues. The maintenance of Chl a and the apparent high baseline heterotrophic capacity of this coral could explain why this genus experiences little to no mortality following bleaching events (McClanahan, 2004; Sutthacheep et al., 2013). In Hawaii, the mounding coral Porites lobata can also bleach severely, but recovers quickly due to a high baseline contribution of heterotrophic carbon to the diet (Palardy et al., 2008; Levas et al., 2013). Thus, the high baseline heterotrophic capacity of the mounding coral F. favus may facilitate the rapid recovery and low mortality of this species following bleaching stress, and allow it to persist into the future, despite its initial dramatic energy reserve and biomass losses when bleached.

Summary

The corals found in the northern tip of the Gulf of Aqaba are able to tolerate temperatures well above their thermal summer maximum for extended periods of time without bleaching, making them temperature-tolerant super-corals. Unlike most corals that bleach after a 4-DHW exposure and experience significant mortality after 8-DHW, all three species of corals here bleached only after being exposed to the equivalent of 32-DHW. However, the physiological and biogeochemical responses to bleaching were species-specific. P. damicornis and S. pistillata maintained their total energy reserves and biomass while F. favus lost total energy reserves and biomass. Interestingly, total energy reserves in all three species of Eilat corals, irrespective of bleaching status, were higher than levels found in corals that were the most susceptible to repeat bleaching in the Caribbean (Grottoli et al., 2014) and consistent with findings that high energy reserves is a common trait found among the most thermally tolerant corals (Rodrigues and Grottoli, 2007; Anthony et al., 2009; Grottoli et al., 2014; Schoepf et al., 2015). While, we found no isotopic evidence of increases in heterotrophy as a mechanism for coping with the declines in photoautotrophically fixed carbon supply in bleached corals, the isotopic signature of F. favus suggests that it has a high baseline input of heterotrophic C, which could underlie its low mortality rate following bleaching events (McClanahan, 2004; Sutthacheep et al., 2013), despite its energy reserve and biomass losses. In addition, endosymbiont type switching or shuffling to acquire or increase thermally tolerant Symbiodinium clade D is probably not an underlying trait to Eilat super-coral thermal tolerance as Symbiodinium D has never been detected in this population of corals (Fine et al., 2013) even though it is present in more southern regions of the Red Sea (Baker et al., 2004; Ziegler et al., 2014). These findings support our hypothesis that energy reserves, biomass, and heterotrophic capacity are central traits to bleaching tolerance in all corals including northern Gulf of Aqaba super-corals, though corals need not have all of these traits to be thermally resilient. In addition, these resilience traits have been identified in a number of coral species throughout the world, suggesting that these phenotypic traits for thermal tolerance may be universal. But, Eilat coral can tolerate much more severe warming events than other populations of corals with similar phenotypic traits, supporting the finding that the Eilat coral thermal tolerance is also driven by the underlying genetics of these coral populations (Fine et al., 2013). Additional study is needed to confirm these findings where the corals are provided zooplankton throughout the experiment. Irrespective of the underlying mechanism for thermal tolerance, super-corals from the northern portion of the Gulf of Aqaba could be used for restocking reefs in other regions where bleaching events have dramatically reduced coral abundance. This may be particularly relevant for Stylophora and P. damicornis corals as they are among those most likely to die following bleaching events in some regions (e.g., Glynn, 1984; Glynn et al., 2001; McClanahan, 2004).

Author Contributions

AG conceived of the study and secured the funding. DT and GW executed the experiment. AG analyzed the samples, statistically analyzed the data, interpreted the results, and wrote the manuscript. DT contributed to the revising of the manuscript.

Funding

Partial funding was provided by startup funds from the University of Pennsylvania and the Ohio State University and by funding from National Science Foundation - Biological Oceanography grant #1459536.

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.

The reviewer MN and handling Editor declared their shared affiliation, and the handling Editor states that the process nevertheless met the standards of a fair and objective review.

Acknowledgments

We would like to thank Ryan Moyer, Dan Mucciarone, and Lindsey Swiaki for their assistance.

References

Alamaru, A., Loya, Y., Brokovich, E., Yam, R., and Shemesh, A. (2009). Carbon and nitrogen utilization in two species of Red Sea corals along a depth gradient: Insights from stable isotope analysis of total organic material and lipids. Geochim. Cosmochim. Acta 73, 5333–5342. doi: 10.1016/j.gca.2009.06.018