rspb.royalsocietypublishing.org

Research Cite this article: Levy O, Karako-Lampert S, Ben-Asher HW, Zoccola D, Page`s G, Ferrier-Page`s C. 2016 Molecular assessment of the effect of light and heterotrophy in the scleractinian coral Stylophora pistillata. Proc. R. Soc. B 283: 20153025. http://dx.doi.org/10.1098/rspb.2015.3025

Received: 17 December 2015 Accepted: 31 March 2016

Subject Areas: ecology, environmental science Keywords: gene expression, coral feeding, bleaching, oxidative stress, microarray, heterotrophy

Authors for correspondence: Oren Levy e-mail: [email protected] Christine Ferrier-Page`s e-mail: [email protected]

Electronic supplementary material is available at http://dx.doi.org/10.1098/rspb.2015.3025 or via http://rspb.royalsocietypublishing.org.

Molecular assessment of the effect of light and heterotrophy in the scleractinian coral Stylophora pistillata Oren Levy1, Sarit Karako-Lampert1, Hiba Waldman Ben-Asher1, Didier Zoccola2, Gilles Page`s3 and Christine Ferrier-Page`s2 1

The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan 5290002, Israel Centre Scientifique de Monaco, 8 quai Antoine 1er, 98000 Monaco, Monaco 3 Universite´ de Nice-Sophia Antipolis, Institut de Recherche sur le Cancer et le Veillissement, IRCAN, CNRS-UMR 7284, INSERM U1081, Nice, France 2

OL, 0000-0002-5478-6307; CF-P, 0000-0002-0357-4486 Corals acquire nutrients via the transfer of photosynthates by their endosymbionts (autotrophy), or via zooplankton predation by the animal (heterotrophy). During stress events, corals lose their endosymbionts, and undergo starvation, unless they increase their heterotrophic capacities. Molecular mechanisms by which heterotrophy sustains metabolism in stressed corals remain elusive. Here for the first time, to the best of our knowledge, we identified specific genes expressed in heterotrophically fed and unfed colonies of the scleractinian coral Stylophora pistillata, maintained under normal and lightstress conditions. Physiological parameters and gene expression profiling demonstrated that fed corals better resisted stress than unfed ones by exhibiting less oxidative damage and protein degradation. Processes affected in lightstressed unfed corals (HLU), were related to energy and metabolite supply, carbohydrate biosynthesis, ion and nutrient transport, oxidative stress, Ca2þ homeostasis, metabolism and calcification (carbonic anhydrases, calciumtransporting ATPase, bone morphogenetic proteins). Two genes (cp2u1 and cp1a2), which belong to the cytochrome P450 superfamily, were also upregulated 249 and 10 times, respectively, in HLU corals. In contrast, few of these processes were affected in light-stressed fed corals (HLF) because feeding supplied antioxidants and energetic molecules, which help repair oxidative damage. Altogether, these results show that heterotrophy helps prevent the cascade of metabolic problems downstream of oxidative stress.

1. Introduction Reefs are degrading rapidly worldwide owing to increased frequency and intensity of anthropogenic disturbances [1]. Under stress events, many coral species bleach, i.e. lose their symbiotic dinoflagellates of the genus Symbiodinium [2]. Because endosymbionts sustain most of the host metabolism through the transfer of photosynthetic products, their loss induces host starvation and can lead to the death of the coral colony [3]. Mass bleaching events are expected to increase in frequency in the future [4], so it is urgent to understand the physiological changes induced by the loss of endosymbionts and how corals can survive without this nutritional source. Whether bleaching leads to mortality depends on several factors: the duration and intensity of the stress event, the tissue thickness, the endosymbiont contribution to the coral metabolism and the heterotrophic capacities of the coral host [5–8]. Heterotrophy—the capture of planktonic prey and the uptake of dissolved organic matter—is one of the main mechanisms by which some species can survive bleaching. Species able to increase their heterotrophic input were shown to be more resilient to bleaching and to recover faster [3,9], because a particulate food supply compensates for the loss of autotrophic products. It also allows for the maintenance of higher symbiont and pigment concentrations within the tissue,

& 2016 The Author(s) Published by the Royal Society. All rights reserved.

(a) Experimental set-up Four colonies of the scleractinian coral Stylophora pistillata, collected in the Red Sea under CITES permit no. DCI/89/32, were used. They were divided into 64 nubbins (16 nubbins per colony), and equally dispatched (four nubbins per colony) into four aquaria (30 l). All colonies were in symbiosis with Symbiodinium clade A1, identified using large, partial chloroplast subunit (23S)-rDNA sequences [15]. Aquaria were kept at 268C + 0.58C with 500 W heaters (Rena) connected to temperature controllers (ElliWell PC 902/T), and were provided with light (photosynthetically active radiation, PAR) from overhead metal halide lamps (OSRAM HQI 400 W) at an intensity of 200 mmoles photons m22 s21. Seawater was renewed at a rate of 2 l h21, mixed with a Rena pump, and no food was supplied, in order to keep the corals under autotrophy. The physico-chemical conditions in the aquaria were closely regulated as shown in electronic supplementary material, figure S1. Nubbins were maintained for seven weeks to monitor their rates of calcification and photosynthesis under these control conditions. For each measurement, one nubbin per colony and one per aquarium (n ¼ 16) were selected. Nubbins were then frozen for determination of protein, chlorophyll (Chla) and symbiont density, as described below. After these first measurements, two tanks remained at 200 mmoles photons m22 s21, whereas PAR was gradually increased over a period of 5 days in the two other aquaria to 500 mmoles photons m22 s21 to induce a light stress. Under each of the two light conditions, nubbins in one aquarium were kept unfed, whereas the others were fed every day for a period of 3 h with Artemia salina nauplii (ca 2000 nauplii l21). Four conditions were achieved: low light unfed (LLU), low light fed (LLF), high light unfed (HLU) and high light fed (HLF). The experiment was performed for 15 days, until a significant decrease in maximal photosynthetic efficiency (Fv/Fm) could be detected by pulse amplitude modulated (PAM) fluorometry [16]. Four nubbins per

(b) Measurements Calcification rates were determined using the buoyant weight method [17]. The same nubbins were then used for photosynthesis and respiration measurements. Rates of respiration (R) in the dark, and net or gross photosynthesis (Pn and Pg) at 200 and 500 mmol photons m22 s21 were assessed for each nubbin using the respirometry [18]. The maximum quantum yield of photosystem II (PSII) fluorescence (Fv/Fm) was also assessed using a diving-PAM fluorometer (Walz GmbH, Germany), as described in electronic supplementary material. Samples were then frozen for determination of total chlorophyll (Chla þ c2), protein and symbiont concentrations, according to Godinot et al. [19]. Data were normalized to the skeletal surface area of each fragment (mmol O2 cm22 h – 1) measured using the wax-dipping technique [20]. Data were checked for normality using the Kolmogorov–Smirnov test with the Lilliefors correction, and for variance homoscedasticity using the Levene test. The effect of treatment was tested using a factorial analysis of variance (ANOVA) with two factors (light and feeding). When there were significant interactions, the analyses were followed by the a posteriori Tukey’s test.

(c) Gene expression analysis Total RNA was extracted with TRIzolw reagent (Invitrogen). DNAse treatment of samples was applied to avoid genomic DNA contamination (DNAse I, Invitrogen). RNA concentrations were measured using a NanoDrop spectrophotometer (ND-1000), and sample quality was checked using a bioanalyser (Agilent). Microarray experiments were performed using 200 ng of labelled RNA samples (Agilent low-input quick Amp labelling kit) hybridized against the microarray in a custom Agilent two-colour gene expression platform with 8  15 K probes per slide. Oligonucleotide probes (60 mers) were designed based on approximately 12 000 genes predicted to encode proteins retrieved from a de novo assembly of 454-sequenced EST libraries of S. pistillata [21]. The intensity of the emitted fluorescence from each target spot on the array was detected using an Agilent G2565BA microarray scanner. Raw and processed data were deposited under the Gene Expression Omnibus, with accession number GSE53661 (http:// www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=wjudmswwtxm vxqr&acc=GSE53661). Stylophora pistillata expressed sequencing tag (EST) data are also stored in the Cnidarian Database of Centre Scientifique de Monaco: http://data.centrescientifique.mc/CSM data-home.html.

(d) Microarray data analysis Measurements from all arrays were subjected to correction and LOESS within-array normalization using Agilent FEATURE EXTRACTION software (Agilent Technologies, Santa Clara, CA). The remaining analyses were performed using Partekw GENOMICS SUITE software (v. 6.6, Copyright 2012, Partek Inc., St Louis, MO). Data from the four biological replicates and two technical replicates were used to perform two-way ANOVA. Genes with significantly up- or downregulated expression (false discovery rate (FDR) p-value , 0.05) were identified, with a cut-off of at least a 1.5-fold change. We then focused on expressed genes or pathways enriched as indicated by DAVID enrichment analysis program [22]. The gene ontology (GO) classification analysis was retrieved from the DAVID database and compared according to the different treatments. The DAVID GO fold change is defined

2

Proc. R. Soc. B 283: 20153025

2. Material and methods

treatment (one per colony, total of n ¼ 16) were snap frozen for gene expression analysis. The remaining two nubbins per colony and treatment (n ¼ 32) were used to assess rates of calcification, photosynthesis, respiration and tissue biomass (Chla, endosymbiont density and protein).

rspb.royalsocietypublishing.org

and thus facilitates the conservation of energy reserves [10,11]. Nutritional status will thus not only determine the resistance and resilience of corals to stress, but also their capacity to grow and reproduce. Although studies have started to assess some of the mechanisms by which heterotrophy sustains coral metabolism, others remain to be understood. In particular, the underlying mechanisms and molecular pathways involved in the physiological changes are unknown. To date, the changes in gene expression in response to feeding under normal or stressful conditions have not been investigated. The aim of this study was to identify genes differentially expressed both in fed and unfed coral colonies maintained under normal and light-stress conditions. Solar bleaching is widespread throughout the tropics, especially on shallow reef flats [12], causing colourless lesions on coral surfaces directly facing the sun. Here we seek to find changes in gene expression which will target specific biological functions that are enhanced or repressed by irradiance levels and feeding behaviour. To complement this, we have combined an array of physiological tools, such as calcification, photosynthesis and tissue analyses. No study to date has explored the association between transcriptome expression response and different feeding and irradiance levels, although heterotrophy is important for coral health, and gene expression has been used to understand the coral metabolic response to environmental changes [13,14]. Such analysis is timely owing to the impact of environmental changes on the future of coral reefs.

(b) 6

10

a

a

3 2

c

1

b 6 4 c

2

HLF

a a

2

b

1

c

0 LLU

LLF

HLU

LLU

0

LLU

b 0.4

0.2 LLF

HLU

HLF

protein concentration (mg cm–2)

a

0.6

LLU

HLF

LLF

HLU

HLF

b –0.5

a

a a

–1.0 –1.5 –2.0 a

(f) a

HLU

HLF

(e) a

LLF

(d) respiration (mmole O2 h–1 cm–2)

HLU

500 mmoles photon m–2 s–1 200 mmoles photon m–2 s–1

3

0.8

a

0 LLF

(c)

skeletal growth (%)

a 8

1.2

a a

0.8

a

0.4

0

LLU

LLF

HLU

HLF

Figure 1. Changes in physiological parameters after 15 days under the different culture conditions (LL, low light; HL, high light; F, fed; U, starved/unfed). Significant differences between treatments are shown with different letters. For Pg, differences were tested at 500 mmol photons m22 s21. All data are represented. as the ratio between the proportion of the submitted list to the proportion of the background one.

(e) Validation of microarray data Reverse transcription was performed using RTIII superscript (Invitrogen) on 2 mg of DNAse-treated RNA using random hexamer primers. Validations were performed using real-time quantitative PCR (qPCR) as described in Moya et al. [23] (electronic supplementary material).

3. Results (a) Coral holobiont physiology Under low light, food provision did not significantly change the physiology of S. pistillata, except that symbiont density was higher in fed corals (figure 1). After exposing corals to high light for 15 days, several parameters were affected, especially in the unfed (HLU) treatment (electronic supplementary

material, table S1). Symbiont density and total chlorophyll content of HLU corals were 80% lower compared with LLU nubbins (figure 1a,b), although the chlorophyll content per symbiont cell remained unchanged (ca 2.3  1026 mg total chl symbiont21; data not shown). The maximal photosynthetic efficiency changed from 0.650 + 0.010 in LLU corals to 0.438 + 0.019 in HLU corals (Tukey’s test, p ¼ 0.001). In addition, rates of net photosynthesis, and respiration per surface area (figure 1c,d), as well as skeletal growth (% increase per day, figure 1e) were 80%, 42% and 45% lower in HLU compared with LLU corals. Only protein content remained unchanged (figure 1f). Conversely to HLU corals, HLF corals decreased their symbiont density and chlorophyll content by only 26% (figure 1a,b), and their rates of net photosynthesis by 45% (figure 1c) compared with LLF corals. The maximal photosynthetic efficiency of the symbionts was slightly reduced from 0.680 + 0.013 in LLF to 0.571 + 0.031 (Tukey’s test, p ¼ 0.02) in HLF corals. Thus, chlorophyll per symbiont cell (ca 1.6  1026 mg total chl symbiont21; results not shown), respiration rates (figure 1d), skeletal growth (figure 1e)

Proc. R. Soc. B 283: 20153025

0 LLU

net photosynthesis (mmole O2 h–1 cm–2)

Chl concentration (mg cm–2)

b

5 4

3

rspb.royalsocietypublishing.org

symbiont density (106 × cell cm–2)

(a)

Venn diagram HLF versus LLF (956)

(a)

4

rspb.royalsocietypublishing.org

113

73

578

Proc. R. Soc. B 283: 20153025

192

167

1784 928

HLU versus HLF (1360)

HLU versus LLF (3482)

(b)

9.50

LLF

LLU

LLF

LLU

LLF

LLU

LLF

LLU

HLF

HLU 8.00

(c)

hierarchical clustering

HLU

LLF

HLF

11.00

Figure 2. (a) Venn diagram shows the gene expression profiles under the two conditions (HLU and HLF versus LLF), where LLF is used as a control for comparative expression. (b) Hierarchical clustering shows heat map of the expression prolife analysed based on the Venn diagram results with biological repeats. (c) The heat map of expression after averaging the biological repeats. and protein content (figure 1f ) remained unchanged between LLF and HLF corals. Therefore, HLF corals maintained an equivalent metabolism compared with LLU corals, except for their total chlorophyll content (figure 1b) and their rates of photosynthesis (figure 1c), which were slightly lower.

(b) Gene expression The principal component analysis captured 37.3% of the total variability for the different treatments (electronic supplementary material, figure S2). Replicate samples were grouped in three clusters. The first cluster included all samples maintained

2712 unique

5

0

–1

–2

–3 –4 log10 (p-value)

–5

–6

rspb.royalsocietypublishing.org

metabolic process cellular ketone metabolic process protein folding ion transmembrane transport ATP metabolic process generation of precursor metabolites and energy cellular process

770 common

0

–0.5 –1.0 –1.5 –2.0 –2.5 –3.0 –3.5 –4.0 –4.5 log10 (p-value)

Figure 3. Variations in global gene expression between treatments explored through GO analysis terms. A total of 2712 genes that significantly changed only in HLU corals compared with control conditions were grouped in seven GO categories, and 770 genes were commonly changed under HLU and HLF in comparison to the control conditions.

under low light (LLF and LLU). HLF samples were grouped in a second cluster, whereas HLU samples constituted a third cluster. There was no significant difference except for two genes with less than twofold in gene expression between LLF and LLU (ANOVA FDR, p-value . 0.05), in agreement with the physiological observations. Therefore, we have decided that each treatment sample will be normalized to the control LLF condition, which represents the best growth condition for the corals, and presented in our Venn diagram (figure 2a), and in the heat map (figure 2b,c). The Venn diagram (figure 2a) of all differentially expressed genes (having an FDR p-value , 0.05) highlighted 770 genes that commonly changed in HLU and HLF versus the control condition. The HLU versus LLF comparison showed a much higher number of differentially expressed genes (total of 3482) than the HLF versus LLF comparison (total of 956 genes; ANOVA FDR, p-value , 0.05, figure 2c). While changing the control condition to HLF, 1360 genes were differentially expressed in the HLU (ANOVA FDR, p-value . 0.05) from which about 87% were already expressed in the LLF control condition expressed lists (figure 2a). These findings suggest that high light induced more cellular changes in starved than in fed corals. Finally, 1616 genes were changed between HLU versus LLU. Variations in global gene expression were further explored through GO analysis (figure 3 and electronic supplementary material, table S2). Six enriched GO categories were changed both in HLU and HLF treatments compared with the control condition, with more than a sixfold change in processes such as electron transport chain, cellular respiration (with genes related to succinate and NADH dehydrogenases (cytochromes b and c)), and carbohydrate biosynthesis (with genes related to glycosyl transferase, glycogen synthase and glucose phosphate isomerase). Other processes affected were oxidoreduction (threefold change, with genes of thioredoxine, NADH, phosphogluconate and succinate dehydrogenases/reductases (cytochromes c and P450)), ketone metabolism and other metabolic processes (1.5-fold change). The above changes in GO

terms were linked to the decrease in rates of photosynthesis observed in both coral groups (HLU and HLF; figure 1). The 2712 genes that significantly changed only in HLU corals compared with the control condition were grouped in seven other GO categories. The highest fold change was observed in the processes of ion transmembrane transport (8.36-fold change), in ATP-metabolic processes (5.20-fold change), protein-folding (4.37-fold change), and in metabolites and energy (3.05-fold change in aco1, aconitase 1 (uptake and sequestration of iron), aifm3 (apoptosis-inducing factor) and ptgs1 (regulation of cell proliferation). Finally, cellular and metabolic processes included 260 and 217 genes, which represented a 1.5-fold change. Most of these processes were linked to the supply of energy by the symbionts, whose density had decreased by more than 80% in the HLU corals, inducing an equivalent decrease in the photosynthesis rate. To gain further insights into the cellular processes occurring during light stress, we investigated the genes that were the most up- or downregulated (cut-off more than twofold change). The analysis pointed to a group of upregulated genes involved in stress response and DNA- or protein-repair processes (figure 4 and table 1): redox regulation (thioredoxin (txnl), peroxidoxin ( prdx5), superoxide dismutase (sodc), methionine sulfoxide (mrsa), quinone oxidoreductase (qor) and hydroxyacylglutathione hydrolase (glo2)), molecular chaperones (heat shock protein 40, 60, 70 (dnaj, hsp60, hsp70)), protein degradation (serine protease: prs23), DNA damage (topoisomerase 2-beta: top2b) and metabolism (enolase (eno2), phosphoglucomutase ( pgm1) and Ca2þ transporting ATPase (atp2b1)). HLU corals exhibited the highest number and fold change of these regulated genes (total of 17) compared with the HLF corals (total of 7). The largest modification in expression level was observed in HLU corals, with an upregulation of two genes (cp2u1 and cp1a2) that belong to the cytochrome P450 superfamily. The fold change increase was 249 for cp2u1 and 10 for cp1a2 in the HLU versus LLF condition, respectively. Conversely, lower values were observed

Proc. R. Soc. B 283: 20153025

oxidation reduction metabolic process elcetron transport chain cellular respiration cellular carbohydrate biosynthetic process cellular ketone metabolic process

HLF versus LLF HLU versus LLF expression 10.0 5.75

–2.75 –7.0

Figure 4. Fold-change values of selected upregulated and downregulated genes. The colour-scale legend indicating the relative fold change is shown on the right. This figure was prepared using EXPANDER software.

under the HLF versus LLF treatment (1.9 for cp2u1 and a nonsignificant 21.3 for cp1a2). In addition, 20 significantly differently expressed genes linked to calcification (calm, fkb14, ca2d2, atp2b1, lrp1, hpcl1, pf2r), calcium pumps (bmpr2, ca2d2, cah1, cah2, cah6, cah12, cah13, cah14, caha) and solute carrier family (s23a2, sc6a1, ctr4, s22af ) were downregulated under HLU versus LLF conditions, whereas a slight decrease (cah14, cah13) or increase (sc6a1) of these genes was observed under LLF and HLF conditions, respectively ( post hoc tests, p-value , 0.05). Three candidate genes that were more downregulated under HLU than under HLF conditions were evaluated using qPCR assays for the validation of the microarray results (electronic supplementary material, figure S3). The expression of calm, s22af and sc6a1 (compared with the housekeeping gene L40A) was significantly higher (ANOVA, p-value , 0.05) in the HLF compared with the HLU.

This study was designed to determine and evaluate the physiological and molecular responses of a scleractinian coral, maintained under autotrophy or both auto- and heterotrophy, to an increase in light intensity affecting its nutritional status. Light-stress-driven changes in gene expression were well correlated with physiological changes, and documented a much higher response in unfed than in fed colonies of S. pistillata, as summarized in figure 5. Comparison of the two feeding conditions clearly showed that unfed colonies were affected by oxidative stress owing to a decrease in metabolic and energy processes. This condition prevented DNA repair and degradation of unfolded protein. Unfed colonies thus experienced a disruption in intracellular Ca2þ homeostasis, leading to a significant decrease in calcification. Conversely, well-fed heterotrophic corals better resisted the stress, because feeding protected them from oxidative damage by supplying antioxidants and energy-rich molecules that contributed to protein/ DNA repair. Bleaching was the physiological response of unfed coral colonies to light stress. These colonies rapidly lost ca 80% of their symbiont density and chlorophyll concentrations, as previously observed in the field [12], reducing photosynthesis (and autotrophic carbon acquisition) to very low levels. Photosystem II was the principal target for damage in these corals, with chronic photoinhibition [24]. Respiration rates were also lower under the HLU condition, and resulted either from a decline in lipid or protein levels, or from an acclimation response to compensate for reduced primary production [11,25]. On the contrary, fed corals lost only 15% of their photosynthetic pigments and slightly decreased their photosynthetic efficiency under high light. Therefore, they responded to irradiance increase with slight bleaching or adaptation to a new light environment [16]. Analysis of the gene expression changes identified molecular pathways involved in the coral responses to light stress. Overall, the fold change in gene expression, and the number of genes differentially expressed, were much lower in HLF than in HLU corals, confirming that feeding allowed coral colonies to avoid the cascade of molecular, and then phenotypic, changes in response to light stress. Genes linked to lipids or carbohydrates were consistently changed in HLU corals: diglyceride acetyltransferase 1 (Dgat1), glyceraldehyde-3-phosphate dehydrogenases (gapdhs) and enolase (enol2). These changes suggest that glycolysis might have been impacted by light stress in HLU corals, which relied on their storage lipids during stress to produce energy and compensate for a decrease in carbon acquisition by the symbionts. Conversely, the provision of external food to HLF corals may have enhanced the nocturnal recovery of the symbiont photosynthetic machinery, allowing the maintenance of high photosynthetic rates and energy supply to the coral host [10,26]. Other results from the gene expression of HLU corals also point to a cellular stress response, which was similar to that monitored in heat-stressed S. pistillata [27] or other coral species and coral larvae [14,27–30]. The processes enriched in lightstressed, unfed corals were involved in oxidative stress, DNA repair, redox regulation and molecular chaperones. The observed changes suggest that high light led to reactive oxygen species formation in the electron transport chain of the symbionts of HLU corals. In addition, in HLU corals, we observed the upregulation of two genes (cp2u1 and cp1a2)

Proc. R. Soc. B 283: 20153025

1.5

6

rspb.royalsocietypublishing.org

FKB14 CAH14 CAH6 CAH13 CAHA CAH2 CTR4 LRP1 CAH1 S23A2 SC6A1 HPCL1 CAH12 BMPR2 CA2D2 PTPRC PF2R CALM S22AF TXNL1 PRDX5 SODC MsrA Se1B QOR GLO2 TOP2B HSP70 DnaJ Hsp60 PRS29 ENOG PGM1 ATP2B1 CP2U1 CP1A2

4. Discussion

Table 1. List of genes significantly up- or downregulated in HLF and HLU corals (fold change .2) compared with those in the control condition. putative function

FKBP14

peptidyl-prolyl cis – trans isomerase

accelerates protein folding

CAH6, 13, 14 CAHA

carbonic anhydrases carbaryl hydrolase

reversible hydration of carbon dioxide hydrolase, carbon-nitrogen ligase

CAH1, 2 CTR4

carbonic anhydrases 1, 2 copper transporter complex

reversible hydration of carbon dioxide copper ion transmembrane transporter

LRP1 SLC23A2

lipoprotein receptor solute carrier family

lipid homeostasis, clearance of apoptotic cells ascorbic acid transporter

SLC6A1

solute carrier family 6 member 1

HPCL1 CAH12

hyppocalcin-like 1 carbonic anhydrase 12

calcium-binding protein reversible hydration of carbon dioxide

BMPR2 CA2D2

bone morphogenic protein receptor 2 a2d calcium channel subunit

bone formation mediates entry of calcium ion into cells

PTPRC

protein tyrosine phosphatase receptor

regulation of cellular processes

PF2R CALM

prostaglandin F receptor calmodulin 1

receptor for prostaglandin calcium-binding protein

S22aF TXNL1

solute carrier family 22 thioredoxin-like 1

organic-ion transporter oxidoreductase activity

PRDX5 SOD

peroxidoxin superoxide dismutase

antioxidant activity antioxidant activity

MsrA

methionine sulfoxide reductase

antioxidant activity

SelB QOR

selenocysteinyl-tRNA-specific quinone oxydoreductase

insertion of selenocystein into proteins quinone reductase activity

GLO2 TOP2B

hydroxyacyl glutathione hydrolase popoisomerase 2

Tiolesterase DNA repair

HSP60, 70

heat shock protein 60, 70

chaperones

DnaJ PRS29

heat shock protein 40 putative serine protease 29

chaperones enzyme-cutting peptide bonds

Enog PGM1

enolase 2 phosphoglucomutase 1

magnesium ion binding breakdown and synthesis of glucose

CP2U1, CP1A2

cytochrome P450

steroid metabolism, lipid biosynthesis

belonging to the cytochrome P450 superfamily, involved in cell detoxification and protection from oxidative stress [31]. cp2u1 (or family 2, subfamily U, polypeptide 1) encodes for an enzyme that catalyses the hydroxylation of arachidonic acid (AA) and other long-chain fatty acids, which protect from oxidative stress by activating peroxisome receptor gamma [32]. The other upregulated cyp gene (family 1, subfamily A, polypeptide 2) or cp1a2 also oxidizes steroids, fatty acids and xenobiotics, and might play the same role as cp2u1. The fact that cp2u1 was upregulated more than 200 times in HLU corals versus two times in HLF ones is consistent with the oxidative stress observed under high light. Because HLU corals were highly stressed by light, genes responsible for sensing and repairing DNA damage were finally activated. Therefore, while oxidative damage began to accumulate in the cells of HLU corals, thioredoxin and peroxiredoxin expression patterns increased to repair proteins and to regulate cellular redox balance. Superoxide dismutase (SOD) and molecular chaperone expressions were also upregulated probably to limit oxidative stress and to assist in protein folding.

Overall, unfed corals experienced a stress response to high light that was similar to the coral response to thermal or acidification stress [27,33,34], suggesting that the same cellprotection mechanisms are conserved among coral species. The decrease in the expression of genes related to energy metabolism and the large decrease in the photosynthetic rates of HLU corals suggest that they were not able to afford the energetic demand of cellular processes such as protein degradation, refolding and DNA repair. Fed corals presented fewer differentially expressed genes (figure 4), and no important bleaching or PSII photoinhibition was observed, suggesting that they were more resistant to high light. The food supply might have brought antioxidant products to corals. A. salina is, indeed, well known for being rich in vitamin E, which is present in unsaturated fatty acids and other antioxidant molecules [35]. This is also true of many other natural plankton species [36]. The involvement of vitamins in the protection against antioxidant stress in fed versus unfed corals is further supported by the downregulation of the solute carrier s23a2 by twofold under the HLU versus HLF corals. s23a2

Proc. R. Soc. B 283: 20153025

annotation

rspb.royalsocietypublishing.org

clone ID

7

8

cell survival

disruption Ca homeostasis

Ca homeostasis

loss of ATP production

ATP production

loss of zoox

conservation of zoox low ROS production

high ROS production light

light

food intake (vitamins, antioxidant) Ca2+

Cah ER

Ca2+

SLC

symbiont

Ca2+

ER

Ca2+

ROS M

M

protein folding ion transport ATP process metabolites energy ketone process electron transport chain carbohydrate synthesis oxido-reduction

Cah

SLC

Ca2+

ROS

symbiont

ATP nutrients

ATP nutrients

HLU

electron transport chain carbohydrate synthesis oxido-reduction

HLF

Figure 5. A schematic conceptual model summarizing the molecular and biological changes induced by high light in unfed and fed corals. Green represents mild stress, whereas red represents severe stress. ER, endoplasmic reticulum; M, mitochondrion; ROS, reactive oxygen species; zoox, zooxanthellae (symbionts). encodes for a sodium/ascorbate co-transporter, which accounts for the specific uptake of vitamin C [37]. In addition to oxidative stress, feeding might have supplied energetic molecules to maintain coral metabolism and to enhance protein repair, which is a process incurring high metabolic costs and energy expenditure [38,39]. Altogether, these genotypic and phenotypic patterns suggest that HLF corals avoided severe bleaching by decreasing oxidative stress and by having more energy to repair photo-oxidative damage. In HLU colonies, oxidative stress also induced a disruption in intracellular Ca2þ homeostasis. This disruption is marked by a sustained increase in intracellular Ca2þ, which led to cytoskeletal rearrangement, cell adhesion changes, decreased calcification and the initiation of cell death. This was evidenced by the downregulation in HLU corals of several genes related to calcium pathways: calmodulin (calm) and fkb14 ( peptidyl-prolyl cis–trans isomerase). Because calm interacts with members of all families of Ca2þ channels at the plasma membrane, endoplasmic reticulum and mitochondria, many other calcium-related genes were also downregulated in HLU corals: at2b1 (a plasma membrane calcium pump that regulates the calcium level in the cells [40]); pf2r (a receptor whose activity is mediated by a phosphatidylinositol –calcium messenger system, and that is involved in the regulation of apoptotic processes [41]); and lrp1 (which encodes a receptor involved in intracellular signalling and clearance of apoptotic cells [42]). A similar disruption in Ca2þ homeostasis was observed in corals experiencing heat-stress [14,43,44], again suggesting that this process is consistently affected among coral species during stress events. This disruption in Ca2þ homeostasis and the almost complete loss of endosymbionts in HLU corals were linked to the severely depressed calcification. We believe that the loss of endosymbionts, indeed, induced a loss of nutrients used to build the

organic matrix [45], whereas the disruption in Ca2þ homeostasis affected the calcium pumps used to pump calcium from the external medium to the site of calcification. Two other significant groups of genes involved in calcification were downregulated in HLU colonies. These genes are related to carbonic anhydrases (CAs), which catalyse the reversible hydration of carbon dioxide into bicarbonate [46], and those related to bone morphogenetic proteins (BMPs). BMPs are members of the transforming growth factor beta superfamily, and play a significant role in coral biomineralization [47]. As far as we know, this is the first report of the downregulation of CAs owing to light stress, although a similar pattern was observed in corals subjected to thermal-induced bleaching [38,48] or acidification stress [27]. Conversely, Vidal-Dupiol et al. [49] observed the upregulation of several CAs under acidification stresses at pH 7.8 and 7.4, but downregulation at a lower pH (7.2). Finally, in HLU corals, several genes related to the solute carrier proteins (SLC) were affected by high light stress and downregulated: SLC6A1 is an amino acid transporter, [50]; SLC26-A6, which functions as a bicarbonate/chloride exchanger in vertebrates [51], was also downregulated in corals experiencing heat stress [43]; SLC7A4 or CTR4 is a cationic amino acid transporter. Its downregulation can impact protein or membrane formation. Finally, PTPRC encodes for a protein that belongs to the tyrosine phosphatase (PTP) family. PTPs are known to be signalling molecules that regulate a variety of cellular processes, including cell growth, differentiation, mitosis and oncogenic transformation [52]. Overall, our results represent the first transcriptomic study comparing the stress response of heterotrophically fed and unfed colonies of the scleractinian coral Stylophora pistillata. Before gene expression can serve as a universal biomarker, we must understand its limitations, particularly with respect to physiological differences in stress tolerance. This case

Proc. R. Soc. B 283: 20153025

Ca2+

rspb.royalsocietypublishing.org

cell death

Data accessibility. The EST sequences reported in this paper were depos-

participated in data analysis and drafted the manuscript. Competing interests. We declare we have no competing interests.

Funding. This study was supported by the Scientific Center of Monaco, and was partially funded by the Israeli Science Foundation (grant no. 243/10) to O. Levy, as well as by the Groupement de Recherche International Coral Reefs.

Acknowledgements. We thank Professor Denis Allemand for the fruitful discussions on the results, Severine Sikorski and Renaud Grepin for performing the qPCR, and Dr Chaim Wachtel (head of the Microarray Unit) and Sivan Goren (Bioinformatics Unit), both from the Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, for performing the experimental part of the microarray and helping with the analysis of the data, respectively. We also thank Professor John Pandolfi as well as two anonymous reviewers for their fruitful comments on the manuscript.

References 1.

Carpenter KE et al. 2008 One-third of reef-building corals face elevated extinction risk from climate change and local impacts. Science 321, 560– 563. (doi:10.1126/science.1159196) 2. Hoegh-Guldberg O. 1999 Climate change, coral bleaching and the future of the world’s coral reefs. Mar. Freshw. Res. 50, 839–866. (doi:10.1071/ MF99078) 3. Grottoli AG, Rodrigues LJ, Palardy JE. 2006 Heterotrophic plasticity and resilience in bleached corals. Nature 440, 1186 –1189. (doi:10.1038/ nature04565) 4. Hoegh-Guldberg O et al. 2007 Coral reefs under rapid climate change and ocean acidification. Science 318, 1737 –1742. (doi:10.1126/science. 1152509) 5. Abrego D, Ulstrup KE, Willis BL, van Oppen MJ. 2008 Species-specific interactions between algal endosymbionts and coral hosts define their bleaching response to heat and light stress. Proc. R. Soc. B 275, 2273 –2282. (doi:10.1098/rspb. 2008.0180) 6. Day T, Nagel L, van Oppen MJ, Caley MJ. 2008 Factors affecting the evolution of bleaching resistance in corals. Am. Nat. 171, E72 –E88. (doi:10.1086/524956) 7. Fitt WK et al. 2009 Response of two species of Indo-Pacific corals, Porites cylindrica and Stylophora pistillata, to short-term thermal stress: the host does matter in determining the tolerance of corals to bleaching. J. Exp. Mar. Biol. Ecol. 373, 102–110. (doi:10.1016/j.jembe. 2009.03.011) 8. Palardy JE, Rodrigues LJ, Grottoli AG. 2008 The importance of zooplankton to the daily metabolic carbon requirements of healthy and bleached corals at two depths. J. Exp. Mar. Biol. Ecol. 367, 180–188. (doi:10.1016/j.jembe.2008.09.015) 9. Rodrigues LJ, Grottoli AG. 2007 Energy reserves and metabolism as indicators of coral recovery from bleaching. Limnol. Oceanogr. 52, 1874 –1882. (doi:10.4319/lo.2007.52.5.1874) 10. Ferrier-Page`s C, Rottier C, Beraud E, Levy O. 2010 Experimental assessment of the feeding effort of

11.

12.

13.

14.

15.

16.

17.

18.

three scleractinian coral species during a thermal stress: effect on the rates of photosynthesis. J. Exp. Mar. Biol. Ecol. 390, 118 –124. (doi:10.1016/j. jembe.2010.05.007) Fitt WK, Cook CB. 2001 The effects of feeding or addition of dissolved inorganic nutrients in maintaining the symbiosis between dinoflagellates and a tropical marine cnidarian. Mar. Biol. 139, 507 –517. (doi:10.1007/s002270100598) Brown BE, Dunne RP, Warner ME, Ambarsari I, Fitt WK, Gibb SW, Cummings DG. 2000 Damage and recovery of photosystem II during a manipulative field experiment on solar bleaching in the coral Goniastrea aspera. Mar. Ecol. Prog. Ser. 195, 117 –124. (doi:10.3354/meps195117) Bay LK, Guerecheau A, Andreakis N, Ulstrup KE, Matz MV. 2013 Gene expression signatures of energetic acclimatisation in the reef building coral Acropora millepora. PLoS ONE 8, e61736. (doi:10. 1371/journal.pone.0061736) Desalvo MK, Voolstra CR, Sunagawa S, Schwarz JA, Stillman TJH, Coffroth MA, Szmant M, Medina M. 2008 Differential gene expression during thermal stress and bleaching in the Caribbean coral Montastraea faveolata. Mol. Ecol. 17, 3952–3971. (doi:10.1111/j.1365-294X.2008.03879.x) Santos SR, Taylor DJ, Kinzie RA, Hidaka M, Sakai K, Coffroth MA. 2002 Molecular phylogeny of symbiotic dinoflagellates inferred from partial chloroplast large subunit (23S)-rDNA sequences. Mol. Phylogenet. Evol. 23, 97 –111. (doi:10.1016/ S1055-7903(02)00010-6) Warner ME, Chilcoat GC, McFarland FK, Fitt WK. 2002 Seasonal fluctuations in the photosynthetic capacity of photosystem II in symbiotic dinoglatellates in the Caribbean reef-building coral Montastraea. Mar. Biol. 141, 31 –38. (doi:10.1007/ s00227-002-0807-8) Jokiel PL, Maragos JE, Franzisket L. 1978 Coral growth: buoyant weight technique. In Coral reef: research methods (eds DR Stoddart, RE Johannes), pp. 379–396. Paris, France: Unesco. Hoogenboom MO, Rodolfo-Metalpa R, Ferrier-Pages C. 2010 Co-variation between autotrophy and

19.

20.

21.

22.

23.

24.

25.

26.

27.

heterotrophy in the Mediterranean coral Cladocora caespitosa. J. Exp. Biol. 213, 2399–2409. (doi:10. 1242/jeb.040147) Godinot C, Ferrier-Page`s C, Sikorski S, Grover R. 2013 Alkaline phosphatase activity of reef-building corals. Limnol. Oceanogr. 58, 227– 234. (doi:10. 4319/lo.2013.58.1.0227) Stimson J, Kinzie III RA. 1991 The temporal pattern and rate release of zooxanthellae from the reef coral Pocillopora damicornis (Linneaus) under nitrogenenrichment and control conditions. J. Exp. Mar. Biol. Ecol. 153, 63–74. (doi:10.1016/S0022-0981(05) 80006-1) Karako-Lampert S et al. 2014 Transcriptome analysis of the scleractinian coral Stylophora pistillata. PLoS ONE 9, e88615. (doi:10.1371/journal.pone. 0088615) Huang da W, Sherman BT, Lempicki RA. 2009 Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44 –57. (doi:10.1038/nprot.2008.211) Moya A, Tambutte´ S, Beranger G, Gaume B, Scimeca JC, Allemand D, Zoccola D. 2008 Cloning and use of a coral 36B4 gene to study the differential expression of coral genes between light and dark conditions. Mar. Biotechnol. (NY) 10, 653–663. (doi:10.1007/s10126-008-9101-1) Jones RJ, Hoegh-Guldberg O, Larkum AWD, Schreiber U. 1998 Temperature-induced bleaching of corals begins with impairment of the CO2 fixation mechanism in zooxanthellae. Plant Cell Environ. 21, 1219–1230. (doi:10.1046/j.1365-3040. 1998.00345.x) Nystro¨m M, Nordemar I, Tedengren M. 2001 Simultaneous and sequential stress from increased temperature and copper on the metabolism of the hermatypic coral Porites cylindrica. Mar. Biol. 138, 1225– 1231. (doi:10.1007/s002270100549) Beraud E, Gevaert F, Rottier C, Ferrier-Page`s C. 2013 The response of the scleractinian coral Turbinaria reniformis to thermal stress depends on the nitrogen status of the coral holobiont. J. Exp. Biol. 216, 2665– 2674. (doi:10.1242/jeb.085183) Kaniewska P, Campbell PR, Kline DI, RodriguezLanetty M, Miller DJ, Dove S, Hoegh-Guldberg O.

9

Proc. R. Soc. B 283: 20153025

ited with GenBank under accession numbers GSE53661 and are available through a website that we established for these data: http://data.centrescientifique.mc/.

Authors’ contributions. All authors carried out the laboratory work,

rspb.royalsocietypublishing.org

study highlights a different stress response for corals relying on heterotrophy, because food provision protected them against antioxidant stress and bleaching, while maintaining symbiont density and photosynthetic performance at a high level during stress (figure 5). Plankton concentration and particulate matter can be high in coral reefs, allowing corals to increase their heterotrophic input when needed [48,53]. Therefore, we hypothesize that nutrition is an important player in the resilience of coral reefs, and heterotrophic capacity will serve as a selective positive force in future environmental changes.

29.

31.

32.

33.

34.

35.

36.

38.

39.

40.

41.

42.

43.

44.

45. Cuif J-P, Dauphin Y, Freiwald A, Gautret P, Zibrowius H. 1999 Biochemical markers of zooxanthellae symbiosis in soluble matrices of skeleton of 24 Scleractinia species. Comp. Biochem. Physiol. 123A,269–278. (doi:10.1016/S1095-6433(99) 00059-8) 46. Moya A, Tambutte´ S, Bertucci A, Tambutte´ E, Lotto S, Vullo D, Supuran CT, Allemand D, Zoccola D. 2008 Carbonic anhydrase in the scleractinian coral Stylophora pistillata characterization, localization, and role in biomineralization. J. Biol. Chem. 283, 25 475 –25 484. (doi:10.1074/jbc.M804726200) 47. Zoccola D, Moya A, Beranger GE, Tambutte´ E, Allemand D, Carle GF, Tambutte´ S. 2009 Specific expression of BMP2/4 ortholog in biomineralizing tissues of corals and action on mouse BMP receptor. Mar. Biotechnol. (NY) 11, 260–269. (doi:10.1007/ s10126-008-9141-6) 48. Heidelberg KB, Sebens KP, Purcell JE. 2004 Composition and sources of near reef zooplankton on a Jamaican forereef along with implications for coral feeding. Coral Reefs 23, 263–276. (doi:10. 1007/s00338-004-0375-0) 49. Vidal-Dupiol J et al. 2013 Genes related to iontransport and energy production are upregulated in response to CO2-driven pH decrease in corals: new insights from transcriptome analysis. PLoS ONE 8, e58652. (doi:10.1371/journal.pone.0058652) 50. Ma B-F, Xie M-J, Zhou M. 2012 Bicarbonate efflux via GABAA receptors depolarizes membrane potential and inhibits two-pore domain potassium channels of astrocytes in rat hippocampal slices. Glia 60, 1761–1772. (doi:10.1002/glia.22395) 51. Gray MA. 2004 Bicarbonate secretion: it takes two to tango. Nat. Cell Biol. 6, 292 –294. (doi:10.1038/ ncb0404-292) 52. Tonks NK. 2013 Protein tyrosine phosphatases from housekeeping enzymes to master regulators of signal transduction. FEBS J. 280, 346– 378. (doi:10. 1111/febs.12077) 53. Yahel R, Yahel G, Genin A. 2005 Near-bottom depletion of zooplankton over coral reefs: I: diurnal dynamics and size distribution. In Coral reefs 24, 75– 85. (doi:10.1007/s00338-004-0449-z)

10

Proc. R. Soc. B 283: 20153025

30.

37.

compared to vitamins in Artemia franciscana nauplii and metanauplii. Mar. Biol. 136, 765– 772. (doi:10. 1007/s002270000280) Takanaga H, Mackenzie B, Hediger MA. 2004 Sodium-dependent ascorbic acid transporter family SLC23. Pflugers Arch. Eur. J. Physiol. 447, 677– 682. (doi:10.1007/s00424-003-1104-1) Downs CA, Fauth JE, Halas JC, Dustan P, Bemiss J, Woodley CM. 2002 Oxidative stress and seasonal coral bleaching. Free Radic. Biol. Med. 33, 533–543. (doi:10.1016/S0891-5849 (02)00907-3) Takahashi S, Nakamura T, Sakamizu M, van Woesik R, Yamasaki H. 2004 Repair machinery of symbiotic photosynthesis as the primary target of heat stress for reef-building corals. Plant Cell Physiol. 45, 251 –255. (doi:10.1093/pcp/pch028) Zoccola D, Tambutte´ E, Kulhanek E, Puverel S, Scimeca J-C, Allemand D, Tambutte´ S. 2004 Molecular cloning and localization of a PMCA P-type calcium ATPase from the coral Stylophora pistillata. Biochim. Biophys. Acta 1663, 117–126. (doi:10. 1016/j.bbamem.2004.02.010) Mu J, Kanzaki T, Si X, Tomimatsu T, Fukuda H, Shioji M, Murata Y, Sugimoto Y, Ichikawa A. 2003 Apoptosis and related proteins in placenta of intrauterine fetal death in prostaglandin f receptordeficient mice. Biol. Reprod. 68, 1968– 1974. (doi:10.1095/biolreprod.102.008029) Yahiro K et al. 2012 Low-density lipoprotein receptor-related protein-1 (LRP1) mediates autophagy and apoptosis caused by Helicobacter pylori VacA. J. Biol. Chem. 287, 31 104–31 115. (doi:10.1074/jbc.M112.387498) Kenkel CD, Meyer E, Matz MV. 2013 Gene expression under chronic heat stress in populations of the mustard hill coral (Porites astreoides) from different thermal environments. Mol. Ecol. 22, 4322–4334. (doi:10.1111/mec.12390) Edge SE, Morgan MB, Gleason DF, Snell TW. 2005 Development of a coral cDNA array to examine gene expression profiles in Montastraea faveolata exposed to environmental stress. Mar. Pollut. Bull. 51, 507 –523. (doi:10.1016/j.marpolbul.2005.07.007)

rspb.royalsocietypublishing.org

28.

2012 Major cellular and physiological impacts of ocean acidification on a reef building coral. PLoS ONE 7, e34659. (doi:10.1371/journal.pone.0034659) Barshis DJ, Ladner JT, Oliver TA, Seneca FO, TraylorKnowles N, Palumbi SR. 2013 Genomic basis for coral resilience to climate change. Proc. Natl Acad. Sci. USA 110, 1387 –1392. (doi:10.1073/pnas. 1210224110) DeSalvo MK, Sunagawa S, Fisher PL, Voolstra CR, Iglesias-Prieto R, Medina M. 2010 Coral host transcriptomic states are correlated with Symbiodinium genotypes. Mol. Ecol. 19, 1174 – 1186. (doi:10.1111/j.1365-294X.2010.04534.x) Kenkel CD et al. 2011 Development of gene expression markers of acute heat-light stress in reefbuilding corals of the genus Porites. PLoS ONE 6, e26914. (doi:10.1371/journal.pone.0026914) Goldstone JV. 2008 Environmental sensing and response genes in cnidaria: the chemical defensome in the sea anemone Nematostella vectensis. Cell Biol. Toxicol. 24, 483–502. (doi:10.1007/s10565-008-9107-5) Wang ZJ, Liang CL, Li GM, Yu CY, Yin M. 2006 Neuroprotective effects of arachidonic acid against oxidative stress on rat hippocampal slices. Chem. Biol. Interact. 163, 207–217. (doi:10.1016/j.cbi. 2006.08.005) Maor-Landaw K, Karako-Lampert S, Ben-Asher HW, Goffredo S, Falini G, Dubinsky Z, Levy O. 2014 Gene expression profiles during short-term heat stress in the red sea coral Stylophora pistillata. Glob. Change Biol. 20, 3026 –3035. (doi:10.1111/gcb.12592) DeSalvo MK, Estrada A, Sunagawa S, Medina M. 2012 Transcriptomic responses to darkness stress point to common coral bleaching mechanisms. Coral Reefs 31, 215–228. (doi:10.1007/s00338011-0833-4) Vismara R, Vestri S, Kusmic C, Barsanti L, Gualtieri P. 2003 Natural vitamin E enrichment of Artemia salina fed freshwater and marine microalgae. J. Appl. Phycol. 15, 75 –80. (doi:10.1023/ a:1022942705496) Maeland A, Ronnestad I, Fyhn HJ, Berg L, Waagbo R. 2000 Water-soluble vitamins in natural plankton (copepods) during two consecutive spring blooms