# Colorimetric Detection of Caspase 3 Activity and Reactive Oxygen Derivatives: Potential Early Indicators of Thermal Stress in Corals

Authors: Mickael Ros; Mathieu Pernice; Sebastien Le Guillou; Martina A. Doblin; Verena Schrameyer; Olivier Laczka
Journal: Journal of Marine Biology
 (2016)

---

## Abstract

There is an urgent need to develop and implement rapid assessments of coral health to allow effective adaptive management in response to coastal development and global change. There is now increasing evidence that activation of caspase-dependent apoptosis plays a key role during coral bleaching and subsequent mortality. In this study, a “clinical” approach was used to assess coral health by measuring the activity of caspase 3 using a commercial kit. This method was first applied while inducing thermal bleaching in two coral species,Acropora millepora andPocillopora damicornis. The latter species was then chosen to undergo further studies combining the detection of oxidative stress-related compounds (catalase activity and glutathione concentrations) as well as caspase activity during both stress and recovery phases. Zooxanthellae photosystem II (PSII) efficiency and cell density were measured in parallel to assess symbiont health. Our results demonstrate that the increased caspase 3 activity in the coral host could be detected before observing any significant decrease in the photochemical efficiency of PSII in the algal symbionts and/or their expulsion from the host. This study highlights the potential of host caspase 3 and reactive oxygen species scavenging activities as early indicators of stress in individual coral colonies.


---

## Body

## 1. Introduction

Tropical marine ecosystems based on scleractinian corals are among the most productive and biologically diverse marine ecosystems on Earth [1]. At the heart of their success is the coral symbiosis with dinoflagellate algae,Symbiodinium sp. (also called zooxanthellae), which live within their tissues and provide corals with the bulk of their required energy [2]. The mutualistic association between corals and their symbiotic dinoflagellates largely explains the ability of corals to build the framework of coral reefs, which serve as a habitat for millions of other species [3]. This framework supports a biodiversity that ultimately underpins tourism, substance gathering, and fishing in at least 100 nations [4].Coral reef ecosystems are undergoing rapid and unprecedented decline worldwide, primarily driven by coastal development and global change [5]. An estimated 27% of the world’s reefs have already been lost and if the current pressures continue unabated, nearly 60% could be lost by the year 2030 [6]. Given the ecosystem services of coral reefs (e.g., $3.8 billion per year for the Great Barrier Reef only, [7]), losing coral from tropical marine ecosystems will have substantial social and economic impacts on coastal populations. There is, therefore, an urgent need to develop and implement rapid assessments of coral health to ensure that if stress is detected, steps can be taken to intervene in an adaptive approach.Current coral monitoring and assessment protocols generally involve major reef surveys of coral cover and species richness to reveal changes in community structure through the loss of individual organisms, populations, and species [8]. These surveys are time-consuming and severely limited by their inability to diagnose coral health in real-time and detect stress prior to the onset of coral mortality. Both factors hinder effective intervention and management to reverse the trajectory of decline. This lack of real-time diagnostic tools makes it difficult to develop reactive and effective management policies, practices, and implementation strategies [9].One common route of coral degradation around the world occurs via bleaching (i.e., the disruption of the coral-dinoflagellate symbiosis), a process known to be primarily caused by elevated sea surface temperature. The variation in susceptibility to coral bleaching is commonly explained by the occurrence of different clades or subclades of symbioticSymbiodinium and their different tolerance to thermal stress [10]. However, there is increasing evidence that the coral host can influence the sensitivity of the coral-dinoflagellate assemblage to stress by regulating a large arsenal of stress proteins [11]. Caspases are a family of cysteine proteases that are divided within three distinct subfamilies (the initiator, executioner, and cytokine activator caspases) and are normally present in the cells as inactive enzymes [12]. Caspase activation is known to trigger programmed cell death or apoptosis in a wide range of metazoans, from corals to humans [13, 14], and is therefore a potential broad indicator of stress in animals. Among the dozen existing caspases, caspases 8, 9, 12, and 3 are known to have pivotal roles in cell death and apoptosis processes, being cleaved and activated by one another under stressful conditions [15, 16]. This amplifying proteolytic cascade leads to the final destruction of cells as part of the cell death and apoptosis process [17]. These mechanisms, present in all multicellular organisms, play a key role in tissue health, development, growth, and responses to stress and pathogens [18]. In corals, the initiation of cell death has been shown to occur in response to ultraviolet radiation or diseases and there is now increasing evidence that activation of caspase-dependent apoptosis plays a key role during coral bleaching and subsequent mortality [19, 20].One of the causes of apoptotic pathway triggering is the presence of reactive oxygen species (ROS), acting as transduction signal in cells [20, 21]. During stressful conditions, ROS quantities overwhelm the capacity of antioxidants to handle them, causing oxidative stress in living cells [22]. Heat stress damages the photosystem II (PSII) in symbionts [23] directly related to ROS production [24–27]. In healthy corals and zooxanthellae, ROS are naturally transformed by antioxidants into less toxic components [22].In human medicine, clinical biomarkers are used to diagnose the condition of a patient and to develop the appropriate medical treatment. Cellular diagnosis is a systematic approach to defining cellular biomarkers based on their functionality within a cell and determining how alterations in the activity of cellular proteins/enzymes may reflect overall cellular performance [28]. In this study, we developed a diagnostic coral assay using a commercially available human caspase 3 detection kit. This assay was applied to detect the change in caspase 3 activity as a measure of coral host health during thermal stress in two species of reef-building corals,Acropora millepora andPocillopora damicornis. We then undertook further investigations onPocillopora damicornis, selected in view of its clear caspase response during the preliminary study, examining in more detail its response to thermal stress, followed by a recovery phase. We hypothesized that increasing levels of ROS can be detected before bleaching is observed based on the common assumption that oxidative stress can initiate coral bleaching. [27]. To test this, we determined the changes in the activity of ROS-scavenging enzymes including catalase [22] and glutathione [27, 29], in parallel with caspase 3 apoptotic enzyme activity [19, 20], photochemical efficiency [30], andSymbiodinium density [31] (Figure 1), during a thermal stress laboratory experiment. By screening the dynamics of this large suite of biomarkers, this study further aimed to provide diagnostic markers of sublethal stress in corals.Figure 1
Model of oxidative stress induced by heat stress in corals leading to bleaching events. Blue shapes stand for measures carried out in the study, using the following methods: total glutathione quantitation (1), caspase 3 (2) and catalase (3) activities measurements, PAM fluorometry (4), zooxanthellae density estimation (5), and visual bleaching monitoring (6). PSI: photosystem I; PSII: photosystem II; TM: thylakoid membrane;O
2
-: superoxide anion; H2O2: hydrogen peroxide; HO•: hydroxyl radical; SOD: superoxide dismutase; CAT: catalase; GSHp: glutathione peroxidase; GSH: glutathione; GSSG: glutathione disulphide; APX: ascorbate peroxidase.


## 2. Materials and Methods

### 2.1. Coral Collection and Maintenance

Two coral species representative of Indo-Pacific coral reef ecosystems,Acropora millepora andPocillopora damicornis, were collected on Heron Island reef (23°33′S, 151°54′E) in the southern Great Barrier Reef, Australia, in February 2013 and 2014. Three healthy-looking colonies ofAcropora millepora and three healthy colonies ofPocillopora damicornis were fragmented using wire cutters. These coral fragments were immediately transported to the aquarium facility at the University of Technology Sydney where they were divided into single upward-growing branch tips (6–8 cm long) and hung using conventional fishing line after a 6-week acclimation period. Coral fragments were then acclimatized in conditions simulating their natural environment (artificial seawater, 26°C; constant and optimal flow conditions within the different tanks; 110 μmol photons m−2 s−1 light intensity with 12 h : 12 h day : night cycle) for 4 weeks prior to the thermal stress experiment to make sure that the perturbation generated through the experimental manipulation of the organisms had no other source but thermal stress.


### 2.2. Experimental Treatments

#### 2.2.1. Thermal-Induced Stress onAcropora millepora andPocillopora damicornis

Two branches of each colony for each coral species were randomly placed into each of six aquaria (10 litres, 3 tank replicates per temperature treatment linked to a common sump controlled using a heater-chiller, total of 12 coral branches per aquarium) under the same conditions as for acclimation (artificial seawater, 26°C; 110μmol photons m−2 s−1 light intensity with 12 h : 12 h day : night cycle, light from 10 am to 10 pm). To explore the effect of temperature on caspase 3 activity, the systems were then subjected to the following experimental regimes: stress (0.5°C increase every 12 hours) and a control group (stable at 26°C, Supplementary Figure 1 in Supplementary Material available online at http://dx.doi.org/10.1155/2016/6825949) for 8 days. Water temperature was recorded every 2 minutes using temperature sensors (Hobo U22 Water Temp Pro v2, Onset Computer Corporation, MA, USA). Before coral branches were harvested for host protein analyses, the photosystem II (PSII) efficiency of symbionts was measured using a Pulse Amplitude Modulated (PAM) fluorometer (as described below) on a randomly selected subset (n
=
3) from the treatment and control tanks at the end of the night cycle (i.e., after 12 h of dark adaptation) on days 0, 2, 4, 6, and 8. Because branching corals often present a large spatial variation in tissue pigmentation and physiological condition along the branch, attention was paid to the selection of coral tissues with similar pigmentation (e.g., at similar distance to the branch tip) and to areas of the coral branch exposed to similar light levels for harvesting the coral branches. The same coral branches were then removed from the tanks, directly snap frozen using liquid nitrogen, and stored at −80°C for further analysis of caspase 3 activity and symbiont cell density.


#### 2.2.2. Thermal-Induced Stress and Recovery Phase onPocillopora damicornis

For each temperature treatment, 36 coral nubbins were placed into a system made of 6 aquaria (4 litres each, 6 nubbins per tank) linked to a 200-litre filtered sump, with artificial seawater flowing continuously under the same 7 conditions as for the acclimation period (Supplementary Figure 2). Experiments were conducted for 16 days under the following treatments: (i) control (stable at26
±
0.8°C (mean ± SD) for the whole experiment), (ii) heat stress (increasing 1
±
0.8°C every 24 hours from 26°C, then remaining at 32°C), and (iii) recovery (increasing 1
±
0.4°C every 24 hours from 26°C to 32°C, then decreasing 2°C every 24 hours to reach 26°C and remain constant for 8 days). Water temperature was recorded every hour with temperature sensors (Hobo U22 Water Temp Pro v2, Onset Computer Corporation, MA, USA) and profiles are shown in Supplementary Figure 3.


### 2.3. Coral Tissue Separation by Water-Picking

The air brush technique [32] was used to remove the living coral tissue from the skeleton in a 50 mL Falcon tube containing 5 mL of cold 0.45 μm filtered seawater. To dissociate the host tissue from the zooxanthellae, the samples were centrifuged twice at 4525 ×g for 8 minutes at 4°C. The supernatant containing the coral tissue was immediately aliquoted and used to determine caspase 3 activities and total protein content. Remaining aliquots of coral tissue extract were snap frozen in liquid nitrogen and kept at −80°C for analyses of ROS-scavenging compounds. Coral skeletons were frozen and kept at −20°C for further total surface area estimation. Zooxanthellae pellets were resuspended in 5 mL of Tris Buffer (10 mM Tris, 125 mM NaCl, 0.1 mM EDTA, pH 7.4). For the first thermal-induced stress experiment, the pellets were cryopreserved and kept at −80°C for cell density estimates whereas, for the recovery experiment, they were stored at 4°C for cell counting on the same day.


### 2.4. Zooxanthellae Counts and Coral Fragments Surface Area Estimation

Zooxanthellae were counted using a haemocytometer (Neubauer-improved cell counting chamber) and a compound light microscope where triplicate counts were made for each sample so that at least 400 cells could be counted, and algal densities could be estimated to within 10% counting error.Surface area of host living tissue was estimated using the wax weight method [33, 34] to consider corallite area and use it for further normalization of the data. To avoid any tissue and water residue, coral fragments were put in bleach for 10 minutes and then dried at 60°C overnight. Coral fragments were weighed before waxing and after being coated with one layer of paraffin wax held at 75°C, dipping in it for 5 seconds and leaving it for 5 minutes at room temperature. For standardisation, 8 cubic metallic objects of known surface area (10.34–60.72 cm2) were treated the same way as coral fragments to produce a standard curve (y
=
0.0601
x
-
0.1413, R
2
=
0.9974). Nubbin surface area could then be estimated using the differential paraffin wax weight.


### 2.5. Determination of Zooxanthellae Photochemical Efficiency

To determine the photochemical efficiency of PSII in zooxanthellae, the maximum quantum yield of PSII was measured with a PAM fluorometer. This noninvasive technique allows the determination of the photosynthetic parametersF
0 (minimum level of Chl a fluorescence) and F
M (maximum level of Chl a fluorescence) upon delivery of a saturating pulse of light. These parameters were then used to calculate F
V
=
F
M
-
F
0 and the maximum quantum yield (F
V
/
F
M) of photosystem II (PSII) [35, 36]. While other parameters such as NPQ, ETRmax, and alpha could have been used in this study, F
V
/
F
M has been widely employed in previous coral studies [30, 37]. Measurements were done after a dark adaptation period, at 6 locations on each coral branch (n
=
3) to capture both upward and downward facing surfaces. They were made with a 3 mm Diving-PAM fibre optic using the following settings: gain = 6, SP width = 0.7 s, and SP intensity ≥3000 μmol photons m−2 s−1. Care was taken to ensure a minimum fluorescence of 300. Measurements were undertaken at the same time-steps as for coral host measurements on days 0, 2, 4, 6, and 8.


### 2.6. Protein Content Quantification in Coral Tissue

In order to quantify the protein content of the coral host tissue (needed to normalize caspase 3 activity), a commercial total protein assay kit was used (BCA Protein Assay Kit, Pierce™, Rockford, USA). Following the manufacturer’s instructions, 25 μL of fresh coral tissue sample was introduced into each well of a Nunc MaxiSorp 96-well plate (Sigma Aldrich, Australia) and left to react for 30 minutes at 37°C with 200 μL of BCA reagent before reading the absorbance at 560 nm with a FLUOstar OPTIMA plate reader (BMG Labtech, Ortenberg, Germany). For each plate analysed, a protein standard curve was additionally obtained using bovine serum albumin (BSA) solutions at concentrations ranging from 0 to 2000 μg mL−1. Each coral sample was analysed in triplicate. The protein concentration in each sample was calculated using the standard curve and expressed in units of μg mL−1.


### 2.7. Detection of Caspase 3 Activity

Caspase 3 activity within the coral host tissue samples was measured using the Caspase 3 Assay Kit, Colorimetric (Sigma Aldrich, Australia). This assay is based on the monitoring of DEVD-p-nitroanilide cleavage in yellow-colored p-nitroaniline (pNA) through time [38]. For each assay, in addition to the coral host tissue samples tested, the following standards were undertaken: (i) a blank which consisted of assay buffer and substrate only, in order to verify that the coloration was not due to substrate degradation through time, (ii) a negative control for each sample containing a caspase 3 inhibitor to determine the signal baseline, and (iii) a positive control containing commercial caspase 3 provided with the kit, allowing the verification of the test efficiency. All samples and standards were loaded as quadruplicates. The temperature was maintained at 37°C during the assay and the kinetics of the reaction were followed for 3 h by measuring absorbance at 405 nm every minute on each well using the plate reader software, directly after addition of the caspase 3 substrate. To convert the absorbance values to the concentration of substrate degraded, a standard curve was produced relating the concentration of pNA with the obtained absorbance values. The reagents were added to the wells to a final volume of 100 μL, as described in Supplementary Table 1.


### 2.8. Catalase Activity Detection

Catalase activity was measured using the Catalase Fluorometric Detection Kit (Enzo Life Sciences, Australia). This assay is based on the measurement of hydrogen peroxide (H2O2) substrate leftovers after its degradation by the catalase enzyme through time [39]. In the presence of H2O2 and horseradish peroxidase (HRP), nonfluorescent ADHP (10-acetyl-3,7-dihydroxyphenoxazine) is oxidized to the highly fluorescent Resorufin which, after excitation at 570 nm, emits fluorescence at 590 nm and can be quantitatively measured and then correlated to catalase activity using a commercial catalase standard provided within the kit. Due to the high sensitivity of this test and the small linearity range of standard (0 to 1 U·mL−1), coral samples had to be diluted up to 250-fold. Each of them was analysed in triplicate.


### 2.9. Total Glutathione Detection

Reduced glutathione (GSH) and glutathione disulphide (GSSG) produced from its oxidation by glutathione peroxidase were quantified in host tissue using the Glutathione Assay Kit (Sigma Aldrich, Australia). The plate reader was set to measure absorbance of 2-nitro-5-thiobenzoate (TNB) continuously produced from 5,5′-dithiobis (2-nitrobenzoic acid) (DTNB) oxidation by GSH (every 100 seconds for 20 minutes at 405 nm). GSSG was also measured by reducing it into GSH by glutathione reductase [40]. Results were calculated from a GSH standard curve ranging from 0 to 50 μM, using technical quadruplicates.


### 2.10. Statistical Analysis

Statistical analyses were conducted using IBM SPSS Statistics 19 software. Due to collection permit limitations, the number of biological replicates used in this study was limited to 3. Kolmogorov-Smirnov and Levene’s tests were used first to test the normality and homoscedasticity of the data, respectively. When both of these conditions were met, an analysis of variance (ANOVA) comparing means across treatments was carried out. When Levene’s test did not confirm the homogeneity of variances, a nonparametric Kruskal-Wallis test was carried out and a nonparametric Mann-WhitneyU test was used to compare the control and thermal stress treatment. When the ANOVA test showed a significant difference, a post hoc Fisher Least Significant Difference (LSD) test was performed in order to determine which treatment was different. Student’s t-test was also performed for each sampling time allowing comparison between the values obtained for the control and for the thermally stressed samples. Results were considered significant at the 5% level.


## 2.1. Coral Collection and Maintenance

Two coral species representative of Indo-Pacific coral reef ecosystems,Acropora millepora andPocillopora damicornis, were collected on Heron Island reef (23°33′S, 151°54′E) in the southern Great Barrier Reef, Australia, in February 2013 and 2014. Three healthy-looking colonies ofAcropora millepora and three healthy colonies ofPocillopora damicornis were fragmented using wire cutters. These coral fragments were immediately transported to the aquarium facility at the University of Technology Sydney where they were divided into single upward-growing branch tips (6–8 cm long) and hung using conventional fishing line after a 6-week acclimation period. Coral fragments were then acclimatized in conditions simulating their natural environment (artificial seawater, 26°C; constant and optimal flow conditions within the different tanks; 110 μmol photons m−2 s−1 light intensity with 12 h : 12 h day : night cycle) for 4 weeks prior to the thermal stress experiment to make sure that the perturbation generated through the experimental manipulation of the organisms had no other source but thermal stress.


## 2.2. Experimental Treatments

### 2.2.1. Thermal-Induced Stress onAcropora millepora andPocillopora damicornis

Two branches of each colony for each coral species were randomly placed into each of six aquaria (10 litres, 3 tank replicates per temperature treatment linked to a common sump controlled using a heater-chiller, total of 12 coral branches per aquarium) under the same conditions as for acclimation (artificial seawater, 26°C; 110μmol photons m−2 s−1 light intensity with 12 h : 12 h day : night cycle, light from 10 am to 10 pm). To explore the effect of temperature on caspase 3 activity, the systems were then subjected to the following experimental regimes: stress (0.5°C increase every 12 hours) and a control group (stable at 26°C, Supplementary Figure 1 in Supplementary Material available online at http://dx.doi.org/10.1155/2016/6825949) for 8 days. Water temperature was recorded every 2 minutes using temperature sensors (Hobo U22 Water Temp Pro v2, Onset Computer Corporation, MA, USA). Before coral branches were harvested for host protein analyses, the photosystem II (PSII) efficiency of symbionts was measured using a Pulse Amplitude Modulated (PAM) fluorometer (as described below) on a randomly selected subset (n
=
3) from the treatment and control tanks at the end of the night cycle (i.e., after 12 h of dark adaptation) on days 0, 2, 4, 6, and 8. Because branching corals often present a large spatial variation in tissue pigmentation and physiological condition along the branch, attention was paid to the selection of coral tissues with similar pigmentation (e.g., at similar distance to the branch tip) and to areas of the coral branch exposed to similar light levels for harvesting the coral branches. The same coral branches were then removed from the tanks, directly snap frozen using liquid nitrogen, and stored at −80°C for further analysis of caspase 3 activity and symbiont cell density.


### 2.2.2. Thermal-Induced Stress and Recovery Phase onPocillopora damicornis

For each temperature treatment, 36 coral nubbins were placed into a system made of 6 aquaria (4 litres each, 6 nubbins per tank) linked to a 200-litre filtered sump, with artificial seawater flowing continuously under the same 7 conditions as for the acclimation period (Supplementary Figure 2). Experiments were conducted for 16 days under the following treatments: (i) control (stable at26
±
0.8°C (mean ± SD) for the whole experiment), (ii) heat stress (increasing 1
±
0.8°C every 24 hours from 26°C, then remaining at 32°C), and (iii) recovery (increasing 1
±
0.4°C every 24 hours from 26°C to 32°C, then decreasing 2°C every 24 hours to reach 26°C and remain constant for 8 days). Water temperature was recorded every hour with temperature sensors (Hobo U22 Water Temp Pro v2, Onset Computer Corporation, MA, USA) and profiles are shown in Supplementary Figure 3.


## 2.2.1. Thermal-Induced Stress onAcropora millepora andPocillopora damicornis

Two branches of each colony for each coral species were randomly placed into each of six aquaria (10 litres, 3 tank replicates per temperature treatment linked to a common sump controlled using a heater-chiller, total of 12 coral branches per aquarium) under the same conditions as for acclimation (artificial seawater, 26°C; 110μmol photons m−2 s−1 light intensity with 12 h : 12 h day : night cycle, light from 10 am to 10 pm). To explore the effect of temperature on caspase 3 activity, the systems were then subjected to the following experimental regimes: stress (0.5°C increase every 12 hours) and a control group (stable at 26°C, Supplementary Figure 1 in Supplementary Material available online at http://dx.doi.org/10.1155/2016/6825949) for 8 days. Water temperature was recorded every 2 minutes using temperature sensors (Hobo U22 Water Temp Pro v2, Onset Computer Corporation, MA, USA). Before coral branches were harvested for host protein analyses, the photosystem II (PSII) efficiency of symbionts was measured using a Pulse Amplitude Modulated (PAM) fluorometer (as described below) on a randomly selected subset (n
=
3) from the treatment and control tanks at the end of the night cycle (i.e., after 12 h of dark adaptation) on days 0, 2, 4, 6, and 8. Because branching corals often present a large spatial variation in tissue pigmentation and physiological condition along the branch, attention was paid to the selection of coral tissues with similar pigmentation (e.g., at similar distance to the branch tip) and to areas of the coral branch exposed to similar light levels for harvesting the coral branches. The same coral branches were then removed from the tanks, directly snap frozen using liquid nitrogen, and stored at −80°C for further analysis of caspase 3 activity and symbiont cell density.


## 2.2.2. Thermal-Induced Stress and Recovery Phase onPocillopora damicornis

For each temperature treatment, 36 coral nubbins were placed into a system made of 6 aquaria (4 litres each, 6 nubbins per tank) linked to a 200-litre filtered sump, with artificial seawater flowing continuously under the same 7 conditions as for the acclimation period (Supplementary Figure 2). Experiments were conducted for 16 days under the following treatments: (i) control (stable at26
±
0.8°C (mean ± SD) for the whole experiment), (ii) heat stress (increasing 1
±
0.8°C every 24 hours from 26°C, then remaining at 32°C), and (iii) recovery (increasing 1
±
0.4°C every 24 hours from 26°C to 32°C, then decreasing 2°C every 24 hours to reach 26°C and remain constant for 8 days). Water temperature was recorded every hour with temperature sensors (Hobo U22 Water Temp Pro v2, Onset Computer Corporation, MA, USA) and profiles are shown in Supplementary Figure 3.


## 2.3. Coral Tissue Separation by Water-Picking

The air brush technique [32] was used to remove the living coral tissue from the skeleton in a 50 mL Falcon tube containing 5 mL of cold 0.45 μm filtered seawater. To dissociate the host tissue from the zooxanthellae, the samples were centrifuged twice at 4525 ×g for 8 minutes at 4°C. The supernatant containing the coral tissue was immediately aliquoted and used to determine caspase 3 activities and total protein content. Remaining aliquots of coral tissue extract were snap frozen in liquid nitrogen and kept at −80°C for analyses of ROS-scavenging compounds. Coral skeletons were frozen and kept at −20°C for further total surface area estimation. Zooxanthellae pellets were resuspended in 5 mL of Tris Buffer (10 mM Tris, 125 mM NaCl, 0.1 mM EDTA, pH 7.4). For the first thermal-induced stress experiment, the pellets were cryopreserved and kept at −80°C for cell density estimates whereas, for the recovery experiment, they were stored at 4°C for cell counting on the same day.


## 2.4. Zooxanthellae Counts and Coral Fragments Surface Area Estimation

Zooxanthellae were counted using a haemocytometer (Neubauer-improved cell counting chamber) and a compound light microscope where triplicate counts were made for each sample so that at least 400 cells could be counted, and algal densities could be estimated to within 10% counting error.Surface area of host living tissue was estimated using the wax weight method [33, 34] to consider corallite area and use it for further normalization of the data. To avoid any tissue and water residue, coral fragments were put in bleach for 10 minutes and then dried at 60°C overnight. Coral fragments were weighed before waxing and after being coated with one layer of paraffin wax held at 75°C, dipping in it for 5 seconds and leaving it for 5 minutes at room temperature. For standardisation, 8 cubic metallic objects of known surface area (10.34–60.72 cm2) were treated the same way as coral fragments to produce a standard curve (y
=
0.0601
x
-
0.1413, R
2
=
0.9974). Nubbin surface area could then be estimated using the differential paraffin wax weight.


## 2.5. Determination of Zooxanthellae Photochemical Efficiency

To determine the photochemical efficiency of PSII in zooxanthellae, the maximum quantum yield of PSII was measured with a PAM fluorometer. This noninvasive technique allows the determination of the photosynthetic parametersF
0 (minimum level of Chl a fluorescence) and F
M (maximum level of Chl a fluorescence) upon delivery of a saturating pulse of light. These parameters were then used to calculate F
V
=
F
M
-
F
0 and the maximum quantum yield (F
V
/
F
M) of photosystem II (PSII) [35, 36]. While other parameters such as NPQ, ETRmax, and alpha could have been used in this study, F
V
/
F
M has been widely employed in previous coral studies [30, 37]. Measurements were done after a dark adaptation period, at 6 locations on each coral branch (n
=
3) to capture both upward and downward facing surfaces. They were made with a 3 mm Diving-PAM fibre optic using the following settings: gain = 6, SP width = 0.7 s, and SP intensity ≥3000 μmol photons m−2 s−1. Care was taken to ensure a minimum fluorescence of 300. Measurements were undertaken at the same time-steps as for coral host measurements on days 0, 2, 4, 6, and 8.


## 2.6. Protein Content Quantification in Coral Tissue

In order to quantify the protein content of the coral host tissue (needed to normalize caspase 3 activity), a commercial total protein assay kit was used (BCA Protein Assay Kit, Pierce™, Rockford, USA). Following the manufacturer’s instructions, 25 μL of fresh coral tissue sample was introduced into each well of a Nunc MaxiSorp 96-well plate (Sigma Aldrich, Australia) and left to react for 30 minutes at 37°C with 200 μL of BCA reagent before reading the absorbance at 560 nm with a FLUOstar OPTIMA plate reader (BMG Labtech, Ortenberg, Germany). For each plate analysed, a protein standard curve was additionally obtained using bovine serum albumin (BSA) solutions at concentrations ranging from 0 to 2000 μg mL−1. Each coral sample was analysed in triplicate. The protein concentration in each sample was calculated using the standard curve and expressed in units of μg mL−1.


## 2.7. Detection of Caspase 3 Activity

Caspase 3 activity within the coral host tissue samples was measured using the Caspase 3 Assay Kit, Colorimetric (Sigma Aldrich, Australia). This assay is based on the monitoring of DEVD-p-nitroanilide cleavage in yellow-colored p-nitroaniline (pNA) through time [38]. For each assay, in addition to the coral host tissue samples tested, the following standards were undertaken: (i) a blank which consisted of assay buffer and substrate only, in order to verify that the coloration was not due to substrate degradation through time, (ii) a negative control for each sample containing a caspase 3 inhibitor to determine the signal baseline, and (iii) a positive control containing commercial caspase 3 provided with the kit, allowing the verification of the test efficiency. All samples and standards were loaded as quadruplicates. The temperature was maintained at 37°C during the assay and the kinetics of the reaction were followed for 3 h by measuring absorbance at 405 nm every minute on each well using the plate reader software, directly after addition of the caspase 3 substrate. To convert the absorbance values to the concentration of substrate degraded, a standard curve was produced relating the concentration of pNA with the obtained absorbance values. The reagents were added to the wells to a final volume of 100 μL, as described in Supplementary Table 1.


## 2.8. Catalase Activity Detection

Catalase activity was measured using the Catalase Fluorometric Detection Kit (Enzo Life Sciences, Australia). This assay is based on the measurement of hydrogen peroxide (H2O2) substrate leftovers after its degradation by the catalase enzyme through time [39]. In the presence of H2O2 and horseradish peroxidase (HRP), nonfluorescent ADHP (10-acetyl-3,7-dihydroxyphenoxazine) is oxidized to the highly fluorescent Resorufin which, after excitation at 570 nm, emits fluorescence at 590 nm and can be quantitatively measured and then correlated to catalase activity using a commercial catalase standard provided within the kit. Due to the high sensitivity of this test and the small linearity range of standard (0 to 1 U·mL−1), coral samples had to be diluted up to 250-fold. Each of them was analysed in triplicate.


## 2.9. Total Glutathione Detection

Reduced glutathione (GSH) and glutathione disulphide (GSSG) produced from its oxidation by glutathione peroxidase were quantified in host tissue using the Glutathione Assay Kit (Sigma Aldrich, Australia). The plate reader was set to measure absorbance of 2-nitro-5-thiobenzoate (TNB) continuously produced from 5,5′-dithiobis (2-nitrobenzoic acid) (DTNB) oxidation by GSH (every 100 seconds for 20 minutes at 405 nm). GSSG was also measured by reducing it into GSH by glutathione reductase [40]. Results were calculated from a GSH standard curve ranging from 0 to 50 μM, using technical quadruplicates.


## 2.10. Statistical Analysis

Statistical analyses were conducted using IBM SPSS Statistics 19 software. Due to collection permit limitations, the number of biological replicates used in this study was limited to 3. Kolmogorov-Smirnov and Levene’s tests were used first to test the normality and homoscedasticity of the data, respectively. When both of these conditions were met, an analysis of variance (ANOVA) comparing means across treatments was carried out. When Levene’s test did not confirm the homogeneity of variances, a nonparametric Kruskal-Wallis test was carried out and a nonparametric Mann-WhitneyU test was used to compare the control and thermal stress treatment. When the ANOVA test showed a significant difference, a post hoc Fisher Least Significant Difference (LSD) test was performed in order to determine which treatment was different. Student’s t-test was also performed for each sampling time allowing comparison between the values obtained for the control and for the thermally stressed samples. Results were considered significant at the 5% level.


## 3. Results

### 3.1. Response ofAcropora millepora andPocillopora damicornis to Thermal-Induced Stress

#### 3.1.1. Detection of Caspase 3 Activity in Coral Hosts Using the Commercial Kit

Caspase 3 activity within the coral host tissue ofAcropora millepora was measured using the commercial kit. As shown in Figure 2, the kinetic response obtained from the coral samples is similar to the one obtained with the commercial caspase 3 (positive control). For the coral samples, the initial absorbance values are higher due to the natural absorbance of the coral tissue. To test for the potential contribution of symbionts to the signal, we investigated the presence of caspase 3 activity within the zooxanthellae extracted from the coral samples. The results showed that no significant activity was found in this fraction of the samples (Figure 2), confirming the specificity of the test towards the coral host.Figure 2
Detection of caspase 3 activity inA. millepora host tissue. Evolution of the absorbance recorded over time in coral host tissue (■) and its negative control (□); with commercial caspase 3 (♦) and its negative control (◊); in the zooxanthellae fraction (●) and its negative control (○); and in the blank (–).


#### 3.1.2. Detection of Caspase 3 and Photochemical Efficiency in Heat-Stressed Corals

Caspase activity showed an upward trend over time in heat-treatedA. millepora, but variability between replicates yielded no significant differences over time, just as for control samples (Figure 3(a)) (Kruskal-Wallis: p
=
0.331). A significant increase in caspase 3 activity was observed between control and stressed corals on day 6 (t-test: t
6
=
-
3.003, n
=
3, and p
=
0.041; average ± SD: caspase 3 activity was 3.8
±
0.3 and 4.3
±
0.2 U (nmol pNa) min−1 mg−1 for control and treatment corals, resp.). This increasing trend appeared to continue on day 8; however the increase was not significant, likely due to the variability among heat-stressed corals probably affecting the analysis (t-test: t
4
=
0.867, n
=
3, and p
=
0.473).Figure 3
Effect of heat stress on caspase 3 activity and photosynthetic efficiency inAcropora millepora (a) andPocillopora damicornis. (b). Caspase 3 activity is expressed as percentage of control; PSII efficiency (maximum quantum yield, F
V
/
F
M) as measuredin situ by Pulse Amplitude Modulated (PAM) fluorometry in corals incubated in control (filled circles) and heat stress treatments (open circles). Error bars represent SD; n
=
3 coral branches for caspase 3 activity and n
=
3 coral branches for PAM fluorometry. S
i
g
n
i
f
i
c
a
n
t
∗ difference (post hoc LSD analysis, p
=
0.05) between treatment and control for the same time point.
(a)
(b)There was a slight decrease in the maximum quantum yield (F
V
/
F
M) with time in stressed corals (average ± SD: 0.64
±
0.05 at day 0 and 0.59
±
0.04 at day 8) (Figure 3(a)). This decrease in the zooxanthellae PSII efficiency in stressed corals was significant over time (Kruskal-Wallis: p
=
0.05) while remaining stable within the control treatment. When comparing the PSII efficiency between the control and heated treatments at each time point, a significant difference was observed on day 6 (t-test: t
4
=
-
3.003, n
=
3, and p
=
0.041) showing a decrease of the F
V
/
F
M in heated corals (average ± SD F
V
/
F
M: control at day 6 = 0.651; treatment at day 6 = 0.576). This decreasing trend seemed to continue on day 8 but the differences were not statistically significant (t-test: t
4
=
-
3.389, n
=
3, and p
=
0.076).Caspase 3 activity trended upward over time in heatedPocillopora damicornis corals while it remained constant in the controls (Figure 3(b)). When comparing the two treatments at each time point, a significant difference appeared on day 4 (t-test at day 4: t
4
=
-
2.919, n
=
3, and p
=
0.043; t-test at day 6: t
4
=
-
6.629, n
=
3, and p
=
0.003) confirming an increase of the caspase 3 activity earlier than inA. millepora (average ± SD activity: control at day 4
=
3.2
±
0.3 U (nmol pNa) min−1 mg−1; treatment at day 4
=
3.8
±
0.3 U (nmol pNa) min−1 mg−1; control at day 6
=
3.1
±
0.4 U (nmol pNa) min−1 mg−1; treatment at day 6
=
5.2
±
0.3 U (nmol pNa) min−1 mg−1). No significant difference was observed on day 8 (t-test: t
4
=
-
1.480, n
=
3, and p
=
0.273).There was a slight decrease in the maximum quantum yield of PSII (F
V
/
F
M) over time in stressed corals (from 0.68
±
0.02 on day 0 to 0.62
±
0.03 on day 8) (Figure 3(b)). The ANOVA shows a statistically significant effect of time on the photosynthetic efficiency of the symbionts (ANOVA with 2 factors, interaction between treatment and time: F
1,20
=
1.003, p
=
0.329; time: F
4,20
=
2.922, p
=
0.047; treatment F
4,20
=
0.988, p
=
0.437), with a significant difference between day 8 and days 0 and 2 (post hoc Fisher LSD, p
<
0.05). When comparing the PSII efficiency between both treatments for the course of the experiment, no significant difference could be observed between stressed and control corals (t-test: t
>
0,05).


#### 3.1.3. Effect of Temperature Stress on Symbiont Cell Density

InAcropora millepora, zooxanthellae density did not change during the course of the experiment in either control or treatment specimens, with an average of 0.44
±
0.04 × 106 cells cm−2 (Figure 4(a)). However, comparison between treatments showed expectedly higher zooxanthellae density on day 2 (Mann-Whitney U test: p
=
0.05), followed by a decrease on day 4 showing significantly higher zooxanthellae densities in the control specimens on day 6 (Mann-Whitney U test: p
=
0.05). The zooxanthellae density ofPocillopora damicornis remained similar in control and treatment specimens for most of the experiment with an average of 0.22
±
0.02 × 106 cells cm−2 (Figure 4(b)). However, on day 8, there were higher zooxanthellae densities in control than in stressed corals (Mann-Whitney U test: p
=
0.05).Figure 4
Effect of heat stress on symbiont cell density inAcropora millepora (a) andPocillopora damicornis (b) incubated in control (filled triangles) and heat stress treatments (open triangles). Error bars represent SD; n
=
3 coral branches for caspase 3 activity and n
=
3 coral branches for PAM fluorometry. S
i
g
n
i
f
i
c
a
n
t
∗ difference (post hoc LSD analysis, p
=
0.05) between treatment and control for the same time point.
(a)
(b)


### 3.2. Response ofPocillopora damicornis to Thermal-Induced Stress and Recovery Potential

#### 3.2.1. Photochemical Efficiencies

Photochemical efficiency (Figure5), measured by maximum quantum yield (F
V
/
F
M), decreased through time (Kruskal-Wallis (K-W): p
=
0.012) for samples within the control (0.67
±
0.02 at day 0 and 0.54
±
0.07 at day 16), the stress (K-W: p
=
0.033) (0.70
±
0.01 at day 0 and 0.63
±
0.01 at day 6), and the recovery treatment (K-W: p
=
0.029) (0.69
±
0.01 at day 0 and 0.59
±
0.03 at day 10). An increasing trend of F
V
/
F
M was detected for samples within the recovery treatment after 12 days of the experiment (0.65
±
0.04), although not significant (p
=
0.127). Photochemical efficiency in the stress treatment was significantly higher than control at day 2 (Mann-Whitney (M-W): p
=
0.05). Photochemical efficiency in recovery treatment is lower than control in stress phase at day 4 (M-W: p
=
0.05) and higher than control during the recovery phase at day 16 (M-W: p
=
0.05).Figure 5
Photosynthetic efficiencies as quantum yield (F
V
/
F
M) measured by PAM fluorometry ofPocillopora damicornis corals (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


#### 3.2.2. Symbiont Densities

Symbiont density (Figure6) in control corals increased significantly at day 8 (M-W: p
=
0.05) (3.85
×
10
5
±
0.95
×
10
5 cells·cm−2 at day 0 and 5.25
×
10
5
±
0.51
×
10
5 cells·cm−2 at day 8 to 5.03
×
10
5
±
1.73
×
10
5 cells·cm−2 at day 16). No significant change in cell density occurred during the recovery treatment. The stress treatment did not show a significant change in symbiont density (K-W: p
=
0.282). Until day 6, there was no significant difference between control and both treatments (p
>
0.05) but at day 8 the difference of cell density became significantly higher for control compared to recovery treatment (M-W: p
=
0.05) until the end of the experiment except for day 12 (M-W: p
=
0.127).Figure 6
Symbiodinium densities in living tissue ofPocillopora damicornis coral fragments (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


#### 3.2.3. Caspase 3 Activities

In control corals, caspase 3 activities remained very low and stable without any significant difference during the experiment (K-W:p
=
0.114) (Figure 7). Stress treatment shows a significant increase (M-W: p
=
0.05) of caspase 3 activity in corals at day 4 (6.3
±
2.2 nmol pNa·min−1·mg protein−1 and 2.7
±
1.1 nmol pNa·min−1·mg protein−1 for the control). Both recovery and stress treatments show a substantial increase (however not statistically significant, M-W: p
=
0.083) in caspase 3 activity at day 6 (resp., 38.8
±
23.6 nmol pNa·min−1·mg protein−1 and 41.6
±
10.9 nmol pNa·min−1·mg protein−1). Recovery phase shows a decrease of caspase 3 activity in corals as soon as temperature goes down. Activity levels of corals in recovery phase stabilize to control activity levels at day 14 (M-W: p
=
0.248, 2.7
±
1.4 nmol pNa·min−1·mg protein−1 for the control and 3.9
±
0.11 nmol pNa·min−1·mg protein−1 for the recovery treatment).Figure 7
Caspase 3 activities measured as DEVD-p-nitroanilide cleavage inPocillopora damicornis coral host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


#### 3.2.4. Glutathione Quantities

Total glutathione levels (Figure8) in control changed significantly through time (K-W: p
=
0.013 increasing from day 0 (11.7
±
2.8 μmol GSH·μg protein−1) to day 6 (69.9
±
23.0 μmol GSH·μg protein−1; M-W: p
=
0.05)). Recovery treatment shows a significant (M-W: p
=
0.05) increase of GSH quantities detectable from day 4 (M-W: p
=
0.127) with very similar GSH levels to the stress treatment (67.00
±
24.0 μmol GSH·μg protein−1 and 65.8
±
28.5 μmol GSH·μg protein−1, resp.) to the next sampling point at day 6 (218.4
±
49.3 μmol GSH·μg protein−1 for the recovery treatment and 153.7
±
16.8 μmol GSH·μg protein−1 for the stress treatment). Recovery phase shows an immediate decrease of GSH levels scaling with temperature. GSH quantities for the recovery treatment in this phase are similar to control levels at day 10 (M-W: p
=
0.827).Figure 8
Total glutathione quantities in coralPocillopora damicornis host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6) during 16 days. Asterisks (∗) stand for a significant difference (α
=
0.05).


#### 3.2.5. Catalase Activities

Catalase activity (Figure9) in control corals was constant from day 0 to day 8 (K-W: p
=
0.402) but increased significantly (M-W: p
=
0.05) on day 10 (from 0.131
±
0.081 U·mL−1·μg protein−1 at day 8 to 0.399
±
0.083 U·mL−1·μg protein−1 at day 10). Compared to the control, no significant differences of catalase activities were observed in both recovery and stress treatments (M-W: p
>
0.05).Figure 9
Catalase activities measured as hydrogen peroxide leftovers inPocillopora damicornis coral host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6) during 16 days.


## 3.1. Response ofAcropora millepora andPocillopora damicornis to Thermal-Induced Stress

### 3.1.1. Detection of Caspase 3 Activity in Coral Hosts Using the Commercial Kit

Caspase 3 activity within the coral host tissue ofAcropora millepora was measured using the commercial kit. As shown in Figure 2, the kinetic response obtained from the coral samples is similar to the one obtained with the commercial caspase 3 (positive control). For the coral samples, the initial absorbance values are higher due to the natural absorbance of the coral tissue. To test for the potential contribution of symbionts to the signal, we investigated the presence of caspase 3 activity within the zooxanthellae extracted from the coral samples. The results showed that no significant activity was found in this fraction of the samples (Figure 2), confirming the specificity of the test towards the coral host.Figure 2
Detection of caspase 3 activity inA. millepora host tissue. Evolution of the absorbance recorded over time in coral host tissue (■) and its negative control (□); with commercial caspase 3 (♦) and its negative control (◊); in the zooxanthellae fraction (●) and its negative control (○); and in the blank (–).


### 3.1.2. Detection of Caspase 3 and Photochemical Efficiency in Heat-Stressed Corals

Caspase activity showed an upward trend over time in heat-treatedA. millepora, but variability between replicates yielded no significant differences over time, just as for control samples (Figure 3(a)) (Kruskal-Wallis: p
=
0.331). A significant increase in caspase 3 activity was observed between control and stressed corals on day 6 (t-test: t
6
=
-
3.003, n
=
3, and p
=
0.041; average ± SD: caspase 3 activity was 3.8
±
0.3 and 4.3
±
0.2 U (nmol pNa) min−1 mg−1 for control and treatment corals, resp.). This increasing trend appeared to continue on day 8; however the increase was not significant, likely due to the variability among heat-stressed corals probably affecting the analysis (t-test: t
4
=
0.867, n
=
3, and p
=
0.473).Figure 3
Effect of heat stress on caspase 3 activity and photosynthetic efficiency inAcropora millepora (a) andPocillopora damicornis. (b). Caspase 3 activity is expressed as percentage of control; PSII efficiency (maximum quantum yield, F
V
/
F
M) as measuredin situ by Pulse Amplitude Modulated (PAM) fluorometry in corals incubated in control (filled circles) and heat stress treatments (open circles). Error bars represent SD; n
=
3 coral branches for caspase 3 activity and n
=
3 coral branches for PAM fluorometry. S
i
g
n
i
f
i
c
a
n
t
∗ difference (post hoc LSD analysis, p
=
0.05) between treatment and control for the same time point.
(a)
(b)There was a slight decrease in the maximum quantum yield (F
V
/
F
M) with time in stressed corals (average ± SD: 0.64
±
0.05 at day 0 and 0.59
±
0.04 at day 8) (Figure 3(a)). This decrease in the zooxanthellae PSII efficiency in stressed corals was significant over time (Kruskal-Wallis: p
=
0.05) while remaining stable within the control treatment. When comparing the PSII efficiency between the control and heated treatments at each time point, a significant difference was observed on day 6 (t-test: t
4
=
-
3.003, n
=
3, and p
=
0.041) showing a decrease of the F
V
/
F
M in heated corals (average ± SD F
V
/
F
M: control at day 6 = 0.651; treatment at day 6 = 0.576). This decreasing trend seemed to continue on day 8 but the differences were not statistically significant (t-test: t
4
=
-
3.389, n
=
3, and p
=
0.076).Caspase 3 activity trended upward over time in heatedPocillopora damicornis corals while it remained constant in the controls (Figure 3(b)). When comparing the two treatments at each time point, a significant difference appeared on day 4 (t-test at day 4: t
4
=
-
2.919, n
=
3, and p
=
0.043; t-test at day 6: t
4
=
-
6.629, n
=
3, and p
=
0.003) confirming an increase of the caspase 3 activity earlier than inA. millepora (average ± SD activity: control at day 4
=
3.2
±
0.3 U (nmol pNa) min−1 mg−1; treatment at day 4
=
3.8
±
0.3 U (nmol pNa) min−1 mg−1; control at day 6
=
3.1
±
0.4 U (nmol pNa) min−1 mg−1; treatment at day 6
=
5.2
±
0.3 U (nmol pNa) min−1 mg−1). No significant difference was observed on day 8 (t-test: t
4
=
-
1.480, n
=
3, and p
=
0.273).There was a slight decrease in the maximum quantum yield of PSII (F
V
/
F
M) over time in stressed corals (from 0.68
±
0.02 on day 0 to 0.62
±
0.03 on day 8) (Figure 3(b)). The ANOVA shows a statistically significant effect of time on the photosynthetic efficiency of the symbionts (ANOVA with 2 factors, interaction between treatment and time: F
1,20
=
1.003, p
=
0.329; time: F
4,20
=
2.922, p
=
0.047; treatment F
4,20
=
0.988, p
=
0.437), with a significant difference between day 8 and days 0 and 2 (post hoc Fisher LSD, p
<
0.05). When comparing the PSII efficiency between both treatments for the course of the experiment, no significant difference could be observed between stressed and control corals (t-test: t
>
0,05).


### 3.1.3. Effect of Temperature Stress on Symbiont Cell Density

InAcropora millepora, zooxanthellae density did not change during the course of the experiment in either control or treatment specimens, with an average of 0.44
±
0.04 × 106 cells cm−2 (Figure 4(a)). However, comparison between treatments showed expectedly higher zooxanthellae density on day 2 (Mann-Whitney U test: p
=
0.05), followed by a decrease on day 4 showing significantly higher zooxanthellae densities in the control specimens on day 6 (Mann-Whitney U test: p
=
0.05). The zooxanthellae density ofPocillopora damicornis remained similar in control and treatment specimens for most of the experiment with an average of 0.22
±
0.02 × 106 cells cm−2 (Figure 4(b)). However, on day 8, there were higher zooxanthellae densities in control than in stressed corals (Mann-Whitney U test: p
=
0.05).Figure 4
Effect of heat stress on symbiont cell density inAcropora millepora (a) andPocillopora damicornis (b) incubated in control (filled triangles) and heat stress treatments (open triangles). Error bars represent SD; n
=
3 coral branches for caspase 3 activity and n
=
3 coral branches for PAM fluorometry. S
i
g
n
i
f
i
c
a
n
t
∗ difference (post hoc LSD analysis, p
=
0.05) between treatment and control for the same time point.
(a)
(b)


## 3.1.1. Detection of Caspase 3 Activity in Coral Hosts Using the Commercial Kit

Caspase 3 activity within the coral host tissue ofAcropora millepora was measured using the commercial kit. As shown in Figure 2, the kinetic response obtained from the coral samples is similar to the one obtained with the commercial caspase 3 (positive control). For the coral samples, the initial absorbance values are higher due to the natural absorbance of the coral tissue. To test for the potential contribution of symbionts to the signal, we investigated the presence of caspase 3 activity within the zooxanthellae extracted from the coral samples. The results showed that no significant activity was found in this fraction of the samples (Figure 2), confirming the specificity of the test towards the coral host.Figure 2
Detection of caspase 3 activity inA. millepora host tissue. Evolution of the absorbance recorded over time in coral host tissue (■) and its negative control (□); with commercial caspase 3 (♦) and its negative control (◊); in the zooxanthellae fraction (●) and its negative control (○); and in the blank (–).


## 3.1.2. Detection of Caspase 3 and Photochemical Efficiency in Heat-Stressed Corals

Caspase activity showed an upward trend over time in heat-treatedA. millepora, but variability between replicates yielded no significant differences over time, just as for control samples (Figure 3(a)) (Kruskal-Wallis: p
=
0.331). A significant increase in caspase 3 activity was observed between control and stressed corals on day 6 (t-test: t
6
=
-
3.003, n
=
3, and p
=
0.041; average ± SD: caspase 3 activity was 3.8
±
0.3 and 4.3
±
0.2 U (nmol pNa) min−1 mg−1 for control and treatment corals, resp.). This increasing trend appeared to continue on day 8; however the increase was not significant, likely due to the variability among heat-stressed corals probably affecting the analysis (t-test: t
4
=
0.867, n
=
3, and p
=
0.473).Figure 3
Effect of heat stress on caspase 3 activity and photosynthetic efficiency inAcropora millepora (a) andPocillopora damicornis. (b). Caspase 3 activity is expressed as percentage of control; PSII efficiency (maximum quantum yield, F
V
/
F
M) as measuredin situ by Pulse Amplitude Modulated (PAM) fluorometry in corals incubated in control (filled circles) and heat stress treatments (open circles). Error bars represent SD; n
=
3 coral branches for caspase 3 activity and n
=
3 coral branches for PAM fluorometry. S
i
g
n
i
f
i
c
a
n
t
∗ difference (post hoc LSD analysis, p
=
0.05) between treatment and control for the same time point.
(a)
(b)There was a slight decrease in the maximum quantum yield (F
V
/
F
M) with time in stressed corals (average ± SD: 0.64
±
0.05 at day 0 and 0.59
±
0.04 at day 8) (Figure 3(a)). This decrease in the zooxanthellae PSII efficiency in stressed corals was significant over time (Kruskal-Wallis: p
=
0.05) while remaining stable within the control treatment. When comparing the PSII efficiency between the control and heated treatments at each time point, a significant difference was observed on day 6 (t-test: t
4
=
-
3.003, n
=
3, and p
=
0.041) showing a decrease of the F
V
/
F
M in heated corals (average ± SD F
V
/
F
M: control at day 6 = 0.651; treatment at day 6 = 0.576). This decreasing trend seemed to continue on day 8 but the differences were not statistically significant (t-test: t
4
=
-
3.389, n
=
3, and p
=
0.076).Caspase 3 activity trended upward over time in heatedPocillopora damicornis corals while it remained constant in the controls (Figure 3(b)). When comparing the two treatments at each time point, a significant difference appeared on day 4 (t-test at day 4: t
4
=
-
2.919, n
=
3, and p
=
0.043; t-test at day 6: t
4
=
-
6.629, n
=
3, and p
=
0.003) confirming an increase of the caspase 3 activity earlier than inA. millepora (average ± SD activity: control at day 4
=
3.2
±
0.3 U (nmol pNa) min−1 mg−1; treatment at day 4
=
3.8
±
0.3 U (nmol pNa) min−1 mg−1; control at day 6
=
3.1
±
0.4 U (nmol pNa) min−1 mg−1; treatment at day 6
=
5.2
±
0.3 U (nmol pNa) min−1 mg−1). No significant difference was observed on day 8 (t-test: t
4
=
-
1.480, n
=
3, and p
=
0.273).There was a slight decrease in the maximum quantum yield of PSII (F
V
/
F
M) over time in stressed corals (from 0.68
±
0.02 on day 0 to 0.62
±
0.03 on day 8) (Figure 3(b)). The ANOVA shows a statistically significant effect of time on the photosynthetic efficiency of the symbionts (ANOVA with 2 factors, interaction between treatment and time: F
1,20
=
1.003, p
=
0.329; time: F
4,20
=
2.922, p
=
0.047; treatment F
4,20
=
0.988, p
=
0.437), with a significant difference between day 8 and days 0 and 2 (post hoc Fisher LSD, p
<
0.05). When comparing the PSII efficiency between both treatments for the course of the experiment, no significant difference could be observed between stressed and control corals (t-test: t
>
0,05).


## 3.1.3. Effect of Temperature Stress on Symbiont Cell Density

InAcropora millepora, zooxanthellae density did not change during the course of the experiment in either control or treatment specimens, with an average of 0.44
±
0.04 × 106 cells cm−2 (Figure 4(a)). However, comparison between treatments showed expectedly higher zooxanthellae density on day 2 (Mann-Whitney U test: p
=
0.05), followed by a decrease on day 4 showing significantly higher zooxanthellae densities in the control specimens on day 6 (Mann-Whitney U test: p
=
0.05). The zooxanthellae density ofPocillopora damicornis remained similar in control and treatment specimens for most of the experiment with an average of 0.22
±
0.02 × 106 cells cm−2 (Figure 4(b)). However, on day 8, there were higher zooxanthellae densities in control than in stressed corals (Mann-Whitney U test: p
=
0.05).Figure 4
Effect of heat stress on symbiont cell density inAcropora millepora (a) andPocillopora damicornis (b) incubated in control (filled triangles) and heat stress treatments (open triangles). Error bars represent SD; n
=
3 coral branches for caspase 3 activity and n
=
3 coral branches for PAM fluorometry. S
i
g
n
i
f
i
c
a
n
t
∗ difference (post hoc LSD analysis, p
=
0.05) between treatment and control for the same time point.
(a)
(b)


## 3.2. Response ofPocillopora damicornis to Thermal-Induced Stress and Recovery Potential

### 3.2.1. Photochemical Efficiencies

Photochemical efficiency (Figure5), measured by maximum quantum yield (F
V
/
F
M), decreased through time (Kruskal-Wallis (K-W): p
=
0.012) for samples within the control (0.67
±
0.02 at day 0 and 0.54
±
0.07 at day 16), the stress (K-W: p
=
0.033) (0.70
±
0.01 at day 0 and 0.63
±
0.01 at day 6), and the recovery treatment (K-W: p
=
0.029) (0.69
±
0.01 at day 0 and 0.59
±
0.03 at day 10). An increasing trend of F
V
/
F
M was detected for samples within the recovery treatment after 12 days of the experiment (0.65
±
0.04), although not significant (p
=
0.127). Photochemical efficiency in the stress treatment was significantly higher than control at day 2 (Mann-Whitney (M-W): p
=
0.05). Photochemical efficiency in recovery treatment is lower than control in stress phase at day 4 (M-W: p
=
0.05) and higher than control during the recovery phase at day 16 (M-W: p
=
0.05).Figure 5
Photosynthetic efficiencies as quantum yield (F
V
/
F
M) measured by PAM fluorometry ofPocillopora damicornis corals (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


### 3.2.2. Symbiont Densities

Symbiont density (Figure6) in control corals increased significantly at day 8 (M-W: p
=
0.05) (3.85
×
10
5
±
0.95
×
10
5 cells·cm−2 at day 0 and 5.25
×
10
5
±
0.51
×
10
5 cells·cm−2 at day 8 to 5.03
×
10
5
±
1.73
×
10
5 cells·cm−2 at day 16). No significant change in cell density occurred during the recovery treatment. The stress treatment did not show a significant change in symbiont density (K-W: p
=
0.282). Until day 6, there was no significant difference between control and both treatments (p
>
0.05) but at day 8 the difference of cell density became significantly higher for control compared to recovery treatment (M-W: p
=
0.05) until the end of the experiment except for day 12 (M-W: p
=
0.127).Figure 6
Symbiodinium densities in living tissue ofPocillopora damicornis coral fragments (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


### 3.2.3. Caspase 3 Activities

In control corals, caspase 3 activities remained very low and stable without any significant difference during the experiment (K-W:p
=
0.114) (Figure 7). Stress treatment shows a significant increase (M-W: p
=
0.05) of caspase 3 activity in corals at day 4 (6.3
±
2.2 nmol pNa·min−1·mg protein−1 and 2.7
±
1.1 nmol pNa·min−1·mg protein−1 for the control). Both recovery and stress treatments show a substantial increase (however not statistically significant, M-W: p
=
0.083) in caspase 3 activity at day 6 (resp., 38.8
±
23.6 nmol pNa·min−1·mg protein−1 and 41.6
±
10.9 nmol pNa·min−1·mg protein−1). Recovery phase shows a decrease of caspase 3 activity in corals as soon as temperature goes down. Activity levels of corals in recovery phase stabilize to control activity levels at day 14 (M-W: p
=
0.248, 2.7
±
1.4 nmol pNa·min−1·mg protein−1 for the control and 3.9
±
0.11 nmol pNa·min−1·mg protein−1 for the recovery treatment).Figure 7
Caspase 3 activities measured as DEVD-p-nitroanilide cleavage inPocillopora damicornis coral host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


### 3.2.4. Glutathione Quantities

Total glutathione levels (Figure8) in control changed significantly through time (K-W: p
=
0.013 increasing from day 0 (11.7
±
2.8 μmol GSH·μg protein−1) to day 6 (69.9
±
23.0 μmol GSH·μg protein−1; M-W: p
=
0.05)). Recovery treatment shows a significant (M-W: p
=
0.05) increase of GSH quantities detectable from day 4 (M-W: p
=
0.127) with very similar GSH levels to the stress treatment (67.00
±
24.0 μmol GSH·μg protein−1 and 65.8
±
28.5 μmol GSH·μg protein−1, resp.) to the next sampling point at day 6 (218.4
±
49.3 μmol GSH·μg protein−1 for the recovery treatment and 153.7
±
16.8 μmol GSH·μg protein−1 for the stress treatment). Recovery phase shows an immediate decrease of GSH levels scaling with temperature. GSH quantities for the recovery treatment in this phase are similar to control levels at day 10 (M-W: p
=
0.827).Figure 8
Total glutathione quantities in coralPocillopora damicornis host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6) during 16 days. Asterisks (∗) stand for a significant difference (α
=
0.05).


### 3.2.5. Catalase Activities

Catalase activity (Figure9) in control corals was constant from day 0 to day 8 (K-W: p
=
0.402) but increased significantly (M-W: p
=
0.05) on day 10 (from 0.131
±
0.081 U·mL−1·μg protein−1 at day 8 to 0.399
±
0.083 U·mL−1·μg protein−1 at day 10). Compared to the control, no significant differences of catalase activities were observed in both recovery and stress treatments (M-W: p
>
0.05).Figure 9
Catalase activities measured as hydrogen peroxide leftovers inPocillopora damicornis coral host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6) during 16 days.


## 3.2.1. Photochemical Efficiencies

Photochemical efficiency (Figure5), measured by maximum quantum yield (F
V
/
F
M), decreased through time (Kruskal-Wallis (K-W): p
=
0.012) for samples within the control (0.67
±
0.02 at day 0 and 0.54
±
0.07 at day 16), the stress (K-W: p
=
0.033) (0.70
±
0.01 at day 0 and 0.63
±
0.01 at day 6), and the recovery treatment (K-W: p
=
0.029) (0.69
±
0.01 at day 0 and 0.59
±
0.03 at day 10). An increasing trend of F
V
/
F
M was detected for samples within the recovery treatment after 12 days of the experiment (0.65
±
0.04), although not significant (p
=
0.127). Photochemical efficiency in the stress treatment was significantly higher than control at day 2 (Mann-Whitney (M-W): p
=
0.05). Photochemical efficiency in recovery treatment is lower than control in stress phase at day 4 (M-W: p
=
0.05) and higher than control during the recovery phase at day 16 (M-W: p
=
0.05).Figure 5
Photosynthetic efficiencies as quantum yield (F
V
/
F
M) measured by PAM fluorometry ofPocillopora damicornis corals (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


## 3.2.2. Symbiont Densities

Symbiont density (Figure6) in control corals increased significantly at day 8 (M-W: p
=
0.05) (3.85
×
10
5
±
0.95
×
10
5 cells·cm−2 at day 0 and 5.25
×
10
5
±
0.51
×
10
5 cells·cm−2 at day 8 to 5.03
×
10
5
±
1.73
×
10
5 cells·cm−2 at day 16). No significant change in cell density occurred during the recovery treatment. The stress treatment did not show a significant change in symbiont density (K-W: p
=
0.282). Until day 6, there was no significant difference between control and both treatments (p
>
0.05) but at day 8 the difference of cell density became significantly higher for control compared to recovery treatment (M-W: p
=
0.05) until the end of the experiment except for day 12 (M-W: p
=
0.127).Figure 6
Symbiodinium densities in living tissue ofPocillopora damicornis coral fragments (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


## 3.2.3. Caspase 3 Activities

In control corals, caspase 3 activities remained very low and stable without any significant difference during the experiment (K-W:p
=
0.114) (Figure 7). Stress treatment shows a significant increase (M-W: p
=
0.05) of caspase 3 activity in corals at day 4 (6.3
±
2.2 nmol pNa·min−1·mg protein−1 and 2.7
±
1.1 nmol pNa·min−1·mg protein−1 for the control). Both recovery and stress treatments show a substantial increase (however not statistically significant, M-W: p
=
0.083) in caspase 3 activity at day 6 (resp., 38.8
±
23.6 nmol pNa·min−1·mg protein−1 and 41.6
±
10.9 nmol pNa·min−1·mg protein−1). Recovery phase shows a decrease of caspase 3 activity in corals as soon as temperature goes down. Activity levels of corals in recovery phase stabilize to control activity levels at day 14 (M-W: p
=
0.248, 2.7
±
1.4 nmol pNa·min−1·mg protein−1 for the control and 3.9
±
0.11 nmol pNa·min−1·mg protein−1 for the recovery treatment).Figure 7
Caspase 3 activities measured as DEVD-p-nitroanilide cleavage inPocillopora damicornis coral host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6 and then stable at 32°C until day 16) during 16 days. Asterisks (∗) stand for a significant difference (p
≤
0.05).


## 3.2.4. Glutathione Quantities

Total glutathione levels (Figure8) in control changed significantly through time (K-W: p
=
0.013 increasing from day 0 (11.7
±
2.8 μmol GSH·μg protein−1) to day 6 (69.9
±
23.0 μmol GSH·μg protein−1; M-W: p
=
0.05)). Recovery treatment shows a significant (M-W: p
=
0.05) increase of GSH quantities detectable from day 4 (M-W: p
=
0.127) with very similar GSH levels to the stress treatment (67.00
±
24.0 μmol GSH·μg protein−1 and 65.8
±
28.5 μmol GSH·μg protein−1, resp.) to the next sampling point at day 6 (218.4
±
49.3 μmol GSH·μg protein−1 for the recovery treatment and 153.7
±
16.8 μmol GSH·μg protein−1 for the stress treatment). Recovery phase shows an immediate decrease of GSH levels scaling with temperature. GSH quantities for the recovery treatment in this phase are similar to control levels at day 10 (M-W: p
=
0.827).Figure 8
Total glutathione quantities in coralPocillopora damicornis host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6) during 16 days. Asterisks (∗) stand for a significant difference (α
=
0.05).


## 3.2.5. Catalase Activities

Catalase activity (Figure9) in control corals was constant from day 0 to day 8 (K-W: p
=
0.402) but increased significantly (M-W: p
=
0.05) on day 10 (from 0.131
±
0.081 U·mL−1·μg protein−1 at day 8 to 0.399
±
0.083 U·mL−1·μg protein−1 at day 10). Compared to the control, no significant differences of catalase activities were observed in both recovery and stress treatments (M-W: p
>
0.05).Figure 9
Catalase activities measured as hydrogen peroxide leftovers inPocillopora damicornis coral host fractions (n
=
3) in control (26°C); recovery treatment (ramping from 26°C at day 0 to 32°C at day 6 then cooling to 26° at day 9); and stress treatment (ramping from 26°C at day 0 to 32°C at day 6) during 16 days.


## 4. Discussion

In this study, we firstly show that caspase 3 activity could be monitored in two species of corals using a detection method designed for human medical applications such as breast cancer diagnosis [41] (Figure 3). In bothAcropora millepora andPocillopora damicornis, temperature stress resulted in a significant increase of caspase 3 activity. Previous studies clearly demonstrate the functional activation of caspase 3 in cnidarians in response to various stressors [19, 20, 42]. However, the time scale of this response has never been compared with other methods of stress detection such as PAM fluorometry, a proxy for symbiont health. In our initial experiment, the increase in caspase 3 activity could be detected in the host tissue after 4 days of thermal stress (i.e., temperature increase of 2°C from ambient to 32°C) inP. damicornis (39% increase in caspase 3 activity compared to the control at day 4 and 87% at day 6) and after 6 days of thermal stress inA. millepora (35% increase in caspase 3 activity compared to the control). In comparison, PAM fluorometry detected a decline in PSII efficiency of the symbionts (F
V
/
F
M) after 6 days inA. millepora with no significant decline inP. damicornis during the experiment (Figure 3). These results indicate that, by using a caspase 3 activity assay, heat stress can be detected before any loss of symbionts or change in PSII efficiency, especially inP. damicornis. Furthermore, even if not statistically significant yet, our data suggest that we may be able to detect the increase in caspase 3 activities after only 2 days of exposure to thermal stress.The lack of statistically significant differences observed in the caspase 3 activity at day 8 between treatments is likely to be related to elevated variability within the stressed samples (Figure3). In fact, this high variability could be explained by the large disparities observed between the coral fragments in terms of tissue cover, with some fragments showing higher sensitivity to thermal stress than others (some fragments had lost most of their tissue whereas others were still undergoing bleaching, showing high caspase 3 activities). According to our observations, these differences were not related to position in the tank or orientation of the coral branches.The recovery experiment showed that early response of the coralPocillopora damicornisto thermal stress could be monitored using biomarkers related to oxidative stress responses. Previous studies clearly demonstrate the role of glutathione [27] in stress response mechanisms allowing the corals to cope with stress through ROS scavenging, and with the same pattern of response [21]. However, the time scale of these responses has never been compared with other methods of stress detection such as PAM fluorometry, a proxy for symbiont health [43]. In our study, the increase in caspase 3 activity and GSH levels could be detected in the host tissue at 4 days of thermal stress (i.e., temperature increase of 4°C) (Figure 8). In comparison, PAM fluorometry showed a significant decline in photosynthetic efficiency of the symbionts (F
V
/
F
M) during the experiment, occurring in both the control and recovery treatment (Figure 5). Furthermore, a loss of symbionts does not induce a loss of photochemical efficiency (and vice versa), impeding the correlation of these values with each other (Figure 6). While these results indicate that F
V
/
F
M is not the optimal way to assess coral health in the case of this thermal stress experiment, other PAM parameters such as NPQ, ETR, alpha, and Ik need to be considered for optimum assessment using PAM fluorometry. Moreover, by using GSH quantitation and caspase 3 activity assays, heat stress response can be detected before any loss of symbionts or change in photochemical efficiency of PSII. This contrast between the two methods in the time needed to detect stress could reflect (i) physiological differences between the coral host and the symbiont and/or (ii) methodological disparity in the scale of measurement; that is, while cellular biomarker assays target the tissue of the whole coral branch (between 1 and 2 cm), PAM fluorometry only allows assessment of parts of the coral branch with intact photosynthetic tissue and has a significantly smaller footprint (0.5 cm diameter of illumination by the PAM fibre).Coral response to heat stress is heavily variable. Several cellular mechanisms lead to coral bleaching [44, 45]; in our experiments, we could visually observe host tissue detachment from its skeleton, which still had zooxanthellae associated with it. Considering the increase of caspase 3 activity with stress (Figure 7), this phenomenon seems to result from apoptosis (regulated cell death) and not necrosis (uncontrolled cell death). This latter type of bleaching does not allow a proper recovery phase for corals, since polyps are no longer present to integrateSymbiodinium cellsde novo. This could explain the constant loss of symbionts in living host tissue with no turn-back. Besides,Symbiodinium cells survival is highly compromised 5 days after expulsion [46] and thus cannot wait for the recolonization of the skeleton by new polyps, while in competition for substrate with other microalgae present in the system. Because the loss of zooxanthellae is delayed after the end of heat stress, density could not be used as a relevant bioindicator of coral bleaching. In contrast the control shows doubling of symbiont density, which has already been reported by Fitt et al. (2009) [47]. This increase could be the result of cell division inside the host compartment. At the same time, an increase of catalase activity occurs in the control, which suggests an involvement of the symbionts during the heat stress. Host and algal fractions were separated prior to analysis, but hydrogen peroxide can diffuse through the membrane and the cell walls of algae [22]. Moreover, even though catalase activity in symbiotic organisms such as the sea anemoneAnemonia viridis is mostly imputed to the host fraction [48], this trend may not be the same in corals such asP. damicornis. As the catalase activity in control is not consistent, and because the catalase response of stressed corals shifts through time, we cannot conclude that this biomarker could be used for coral health monitoring (Figure 9). To determine the origin of this response, catalase activity and hydrogen peroxide levels in symbionts should be measured in future studies.


## 5. Conclusion

In diagnostic application, nondestructive methods such as PAM fluorometry are very attractive. However, we showed that detectable changes in the photochemical efficiency of symbionts occur after changes in the concentration of host biomarkers. Clearly, if intervention to prevent coral loss is a desired management outcome, then development of rapid sublethal indicators of coral stress is imperative. We have demonstrated that GSH and caspase 3 assays can detect declines in coral health prior to the onset of thermal-related bleaching and mortality and therefore could be implemented as an early indicator of stress in corals, complementing traditional coral reef monitoring protocols. Nevertheless, the large disparities observed between the coral fragments in terms of tissue cover, with some fragments showing higher sensitivity to thermal stress than others (some fragments had lost most of their tissue whereas others were still healthy or in bleaching process, showing high caspase 3 activity and GSH levels), do not allow us to determine sublethal values of biomarkers that could be used universally. Furthermore, several challenges remain to be addressed before any possible integration of this sublethal detection tool in coral monitoring programs. Firstly, given that the methodology described here is invasive, the impact of sample collection on coral populations should be kept to a minimum. Then, caspase and GSH assays need to be performed on fresh living corals without handling stress, interfering with the specific thermal stress response (these proteins are not stable after freezing), which could be an issue in remote areas. In this context, follow-up studies to increase the sensitivity (thereby reducing the amount of coral tissue needed) using electrochemical sensing of caspase 3 activity [48] or near-infrared fluorescence to detect intracellular levels of GSH in vivo [49] show great promises. There is also a need to (i) target the variation among coral individuals using increased numbers of biological replicates; (ii) screen the involvement of different clades ofSymbiodinium in GSH levels and caspase 3 activity; (iii) improve resolution of results by studying response of other coral species; (iv) determine trigger values of caspase 3 activity in relation to a broad range of stressors (not only ROS) and assess their recovery following stress exposure; and (v) validate this technique in field studies. Indeed, in the field, the temperature and light levels are not constant and can fluctuate on a diel rhythm. Further, seasonal changes might help corals to prepare for elevated summer temperatures and light levels.While this new tool cannot substitute the need to monitor the effect of coral stress at the ecosystem-scale level, our work provides further evidence that it could be used to monitor health of individual coral colonies and help coral reef management authorities or public aquariums to detect and adaptively manage the effect of environmental stressors.



---

*Source: 6825949-2016-03-03.xml*