Impact of El Niño-Southern Oscillation 2015-2016 on the soluble proteomic profile and cytolytic activity of Millepora alcicornis (“fire coral”) from the Mexican Caribbean

View article
Aquatic Biology

Introduction

Coral reefs play a critical role in marine ecology and human sustainability (Jain, Sonawane & Mandrekar, 2008; Mumby & Van Woesik, 2014; Anthony, 2016; Williams et al., 2016). Moreover, the organisms that constitute these ecosystems are considered a rich source of novel bioactive agents with considerable pharmaceutical and biotechnological potential, since they produce a great variety of molecules with unique structural features (Jain, Sonawane & Mandrekar, 2008; Cragg & Newman, 2013; Mayer et al., 2016; Li, Ding & Li, 2017). Organisms of the genus Millepora, which belong to the class Hydrozoa, are considered the second most important reef forming species, after scleractinian hard corals of the class Anthozoa (Lewis, 2006; Ruiz-Ramos, Weil & Schizas, 2014). In addition to their ecological importance, these organisms, commonly known as “fire corals,” induce local and systemic toxic effects in humans (Lewis, 2006; Rojas-Molina et al., 2012).

Unfortunately, coral reefs are extremely vulnerable to the stress related to greenhouse gas emissions, mainly ocean warming (Lesser, 2007; Baird et al., 2009; Mumby & Van Woesik, 2014; Lough, 2016). Elevated seawater temperature provokes disturbances that can seriously affect and break down the homeostatic capacity of coral reefs to overcome stressors (Mora, Graham & Nyström, 2016). One of the most devastating consequences of global warming is coral bleaching, in which corals and hydrocorals lose their photosynthetic symbiotic algae of the genus Symbiodinium or their pigments, which exposes the white exoskeleton composed of calcium carbonate (Hoegh-Guldberg, 1999; Lesser, 2006; Bonesso, Leggat & Ainsworth, 2017; Neal et al., 2017).

The frequency and seriousness with which this phenomenon occurs have increased in recent years. Massive bleaching events have been recorded in all the tropical regions of the world in 1987, 1998, 2003, 2005, and 2010 (Eakin et al., 2010; Hoegh-Guldberg & Bruno, 2010; Heron et al., 2016). In fact, during 2014–2017, the worst documented bleaching event was observed (Eakin et al., 2017; Hughes et al., 2017b, 2018). Numerous investigations have demonstrated that bleaching is deleterious to coral reefs. Bleaching events have caused massive damage to coral reefs around the world, with very severe effects on the balance of biodiversity in marine tropics. This phenomenon causes a general deterioration in reef health, as it provokes an increase in coral diseases, a decay in reef calcification, breakdown of reef framework by bioeroders, and loss of essential habitat for related reef organisms (Baker, Glynn & Riegl, 2008; Lesser, 2011; Grottoli et al., 2014; Okazaki et al., 2016; Hughes et al., 2017a). Some climate models predict that if CO2 emissions continue at the current rate, bleaching events will increase its frequency and severity, seriously threatening the survival of coral reefs (Houlbreque & Ferrier-Pagès, 2009; Lesser, 2011; Mumby & Van Woesik, 2014; Lough, 2016; Neal et al., 2017; Hughes et al., 2018).

Studies that have addressed the etiology and effects of bleaching have focused on Anthozoan species. To date there are no reports regarding the impact of thermal stress on the expression of genes and proteins of Hydrozoan species. Genomic, transcriptomic, and proteomic approaches have revealed that heat stress causes a differential pattern in the expression of genes and proteins from various cellular processes, including oxidative stress response, Ca2+ homeostasis, cytoskeletal organization, protein synthesis, apoptosis, endo-exophagocytosis, and immune response. Also, the differential expression of heat shock proteins and transcription factors has been observed (Bay & Palumbi, 2015; Pinzón et al., 2015; Raina et al., 2015; Seneca & Palumbi, 2015; Weston et al., 2015; Maor-Landaw & Levy, 2016; Oakley et al., 2016, 2017; Ricaurte et al., 2016; Sloan & Sawyer, 2016; Huang et al., 2017; Ruiz-Jones & Palumbi, 2017; Traylor-Knowles et al., 2017; Mayfield et al., 2018a).

It is generally accepted that cnidarian-Symbiodinium symbiosis is fundamental for the formation of coral structures, since algae carry out photosynthesis and transfer more than 50% of their photosynthetic products to the cnidarian host (Houlbreque & Ferrier-Pagès, 2009; Davy, Allemand & Weis, 2012; Fransolet, Roberty & Plumier, 2012; Douglas, 2003; Venn, Loram & Douglas, 2008; Yellowlees, Rees & Leggat, 2008). Up to date, the role of heterotrophic carbon (C) input in the resilience and recovery of bleached reef forming cnidarians is still controversial. It has been observed that during bleaching certain species can increase their heterotrophic feeding to maintain and restore energy reserves (Grottoli, Rodrigues & Palardy, 2006; Aichelman et al., 2016). However, some studies suggest that heterotrophically derived fixed carbon does not completely attenuate C budget imbalance (Levas et al., 2016; Tremblay et al., 2016).

Undoubtedly, after a bleaching event the central metabolism of both symbiotic partners is significantly altered (Ricaurte et al., 2016; Hillyer et al., 2017, 2018; Oakley et al., 2017; Mayfield et al., 2018a), and the contribution of heterotrophy in the response of cnidarian-Symbiodinium symbiosis to elevated temperature is still unclear (Grottoli, Rodrigues & Palardy, 2006; Levas et al., 2016; Aichelman et al., 2016; Tremblay et al., 2016). Considering that toxins are essential in prey capture and digestion, it is very likely that thermal stress, which provokes bleaching, alters their expression in cnidarians. Therefore, assessment of the effect of thermal stress on the expression of cnidarian toxins used to capture preys is fundamental to understand if heterotrophy constitutes a significant survival mechanism of reef-forming cnidarians after a bleaching phenomenon. In this context, the aim of the present study was to determine changes in the soluble proteomic profile and cytolytic activity of Millepora alcicornis (“fire coral”) that suffered bleaching during El Niño-Southern Oscillation (ENSO) 2015–2016 in the Mexican Caribbean Sea.

Materials and Methods

Specimen collection

Hydrocoral fragments were collected by removing small portions from the edges of three normal and three visibly bleached M. alcicornis colonies in November 2016, within the third global-scale mass bleaching event ever documented (Hughes et al., 2017b, 2018). Hydrocorals were obtained at depths of 4–10 m. The possibility of sampling genetic clones was reduced by collecting fragments that were at least 10 m apart in the area known as “La Bocana Chica” that belongs to the Parque Nacional Arrecife de Puerto Morelos (Quintana Roo, México). Collection permit was issued by Secretaría de Agricultura, Desarrollo Rural, Pesca y Alimentación (permit number PPF/DGOPA-139/15). Specimens were frozen in liquid nitrogen and transported to our laboratory at the Universidad Autónoma de Querétaro.

Determination of cell density in host tissues

Bleached and normal fragments of M. alcicornis were cut into one cm2 squares, subsequently, tissues were fixed with 10% buffered formalin during 3 days. Thereafter, samples were decalcified in a 5% solution of HCl for 5 days. The decalcifying solution was replaced daily. After this treatment, tissues were homogenized in a Glass-Col tissue homogenizer system during 2 min at 70 rpm, and centrifuged at 2,400 rpm. The resulting pellet was resuspended in 70% ethanol. A total of 10 μl-aliquots of this homogenate were used to count Symbiodinium cells using a Neubauer hemocytometer. Cell counts were performed in triplicate and were standardized in relation to the total surface area of one cm2.

Determination of chlorophyll content in Symbiodinium

Chlorophyll a and c2 contents were quantified employing the method of Shirur et al. (2014). Briefly, normal and bleached hydrocoral fragments were incubated in the dark for 24 h at 4 °C with a mixture of acetone: dimethyl sulfoxide (95:5 v/v). Chlorophylls a (Chl a) and c (Chl c2) contents were calculated using the equations of Jeffrey & Humphrey (1975) and were then standardized to M. alcicornis surface area.

Soluble proteome extraction

Nematocyst discharge and soluble proteome extraction from normal and bleached M. alcicornis fragments were accomplished by stirring the hydrocoral fragments in deionized water (pH 7) at 4 °C for 24 h. The aqueous extracts obtained, containing the soluble proteomes, were centrifuged at 12,000 rpm for 15 min at 4 °C. Subsequently, the supernatants were filtered through a 0.45 μm pore filter (Millipore Ltd, Hertfordshire, England). The filtrates were lyophilized and stored at −70 °C. For proteomic analysis, the lyophilized filtrates were dissolved in 100 μl of deionized water employing microtubes for centrifuge. Thereafter, proteins were cleaned using a Ready Prep 2-D Cleanup Kit (Bio-Rad, Hercules, CA, USA). The protein concentration in both soluble proteomes was determined using a 2D-Quant Kit (GE Healthcare Life Sciences, Uppsala, Sweden).

SDS–PAGE

Sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS–PAGE) was carried out as previously described (Schägger & Von Jagow, 1987). Samples of normal and bleached M. alcicornis-soluble proteomes containing 80 μg of protein were run in 18% polyacrylamide gels at 150 V for 2 h at 4 °C, using tris-glycine as buffer. Protein bands were visualized with BioSafe Coomassie blue G-250 dye (Bio-Rad, Hercules, CA, USA). Molecular weights were determined using Precision Plus Protein™ standards (Bio-Rad, Hercules, CA, USA).

First-dimension step, isoelectric focusing

One milligram of protein from the soluble proteomes of normal and bleached hydrocorals was mixed with rehydration buffer (8M Urea, 2% CHAPS, 50 mM dithiothreitol (DTT), 0.2% BioLyte® 3/10 Ampholyte) overnight on immobilized pH gradient using 11 cm ReadyStrip™ IPG strips pH 3–10 (Bio-Rad, Hercules, CA, USA), which were then subjected to first-dimension isoelectric focusing (IEF) in a PROTEAN® i12™ IEF System (Bio-Rad, Hercules, CA, USA) for a total of 20,000 Vh, as previously described (Ricaurte et al., 2016).

High-resolution two-dimensional gel electrophoresis

After completion of the IEF, and prior to running the second dimension, strips were firstly submerged in equilibration buffer I (urea 6M, 2% SDS, Tris–HCl 0.05M, pH 8.8, 50% glycerol, and 2% (w/v) DTT) to reduce the S–S bonds, and subsequently alkylated in equilibration buffer II (urea 6M, 2% SDS, Tris–HCl 0.05M, pH 8.8, 50% glycerol and 2.5% (w/v) Iodoacetamide). The second-dimensional separation was carried out in 18% SDS-polyacrylamide gels at 150 V for 2 h at 4 °C, using tris-glycine buffer. Gels containing the soluble proteomes of normal and bleached M. alcicornis were stained with BioSafe™ Coomassie blue G-250 dye (Bio-Rad, Hercules, CA, USA).

Image analysis

Images of the stained gels were captured in a ChemiDoc™ MP (Bio-Rad, Hercules, CA, USA). imaging system at 600 dpi resolution. Spot detection, matching, and fold changes were determined with the PD-Quest™ software (Bio-Rad, Hercules, CA, USA), version 8.0.1. Protein spots with fold change ≥2 between proteomes of normal and bleached specimens were regarded as differentially expressed. These spots were selected and marked for excision. All experiments were performed in triplicates.

Protein in-gel digestion and data analysis

Selected spots were manually excised and faded with ACN: 50 mM NH4HCO3 (50:50 v/v). Protein spots were processed for protein identification at Instituto Nacional de Medicina Genómica, México. Protein samples were subjected to enzymatic digestion for 18 h at 37 °C with mass spectrometry-grade-trypsin (V528A; Promega, Madison, WI, USA). Peptides obtained from the digestion (ACN:H2O:formic acid 50:45:5 v/v) were extracted and the volume of samples was reduced in a concentrator Eppendorf 5301 (Eppendorf, Hamburg, Germany). Samples were desalted by a C18 column (ZipTipC18 15 μm; Millipore, Bedford, MA, USA). Subsequently, a MALDI-TOF/TOF tandem mass spectrometry analysis was carried out to identify the selected proteins. Afterward, each protein sample was co-crystallized in plates using an α-Cyano-4-hydroxycinnamic acid matrix. MALDI-TOF/TOF experiments were carried out on a 4800 MALDI-TOF/TOF Analyzer (Applied Biosystems, Foster City, CA, USA) equipped with a 200 Hz, 355 nm Nd:YAG laser, linked to 4,000 series Explorer software version 3.5.3. The mass spectrometer was operated in positive ion mode and externally calibrated using a peptide mass standard. The laser source was set to 2,500–2,800 for MS and 3,500–3,800 for MS/MS acquisition. MS spectra were recorded in positive ion reflector mode with 25 laser shots. Precursor ions were automatically selected for fragmentation, which was carried out employing a collision energy of two kV, using air as the collision gas (pressure, 2 × 10–6 torr), with the accumulation of 400 shots. Mass spectra were acquired within a mass range of 800–4,000 Da (m/Z). The parental ion of Glu1-Fibrino-Peptide B, diluted in the matrix (1.3 pmol/μl/spot), was employed for internal calibration at m/z = 1,570.690 Da. The most intense ion signals per spot position having an S/N >20 were chosen for MS/MS acquisition. Following MALDI-TOF/TOF analysis, a search was performed, employing the Protein Pilot™ software version 2.0 (Applied Biosystems, Foster City, CA, USA). Carbamidomethylation of cysteine was set as fixed modification. The Basic Local Alignment Search Tool (BLASTp) (https://blast.ncbi.nlm.nih.gov/Blast.cgi) was used for identification of homolog sequences. Database was set to non-redundant protein sequences, using Cnidaria as “organism parameter search.” Similar sequences were aligned with the multiple sequence alignment program ClustalW2 (https://www.ebi.ac.uk/Tools/msa/clustalw2/). The criterion for the selection of protein spots considered the comparison of experimental isoelectric point (pI), molecular weight, and sequence coverage of proteins with that of the theoretical pI and molecular weight of the proteins detected on the 2DE-PAGE gels.

Evaluation of cytolytic activity of the soluble proteomes from normal and bleached M. alcicornis specimens

Comparative hemolytic activity of the soluble proteomes from normal and bleached M. alcicornis specimens

Hemolytic activity of the soluble proteomes from normal and bleached M. alcicornis specimens was monitored according to a method previously described with some modifications (Hernández-Matehuala et al., 2015). Briefly, samples for the assay contained a mixture (one mL) of Alsever’s solution (120 mM D-glucose, 30 mM sodium citrate, seven mM NaCl, and two mM citric acid; pH 7.4) with 50 μl of a 1% suspension of erythrocytes from rat and the required volume of each soluble proteome. These samples were incubated at 37 °C for 30 min. After centrifugation at 1,430 g for 4 min at 4 °C, the absorbance at 415 nm of the supernatant fluid containing the hemoglobin released from lysed erythrocytes was measured in a Benchmark Plus microplate spectrophotometer (Bio-Rad Laboratories, Hercules, CA, USA). Each experiment was normalized with respect to complete hemolysis, which was measured by diluting the erythrocyte sample in deionized water instead of Alsever’s buffer. One hemolytic unit (HU50) was defined as the amount of protein sample required to cause 50% hemolysis. The hemolytic activity was plotted in GraphPad Prism 6.0 (GraphPad Software, San Diego, CA, USA).

Comparative phospholipase A2 activity of the soluble proteomes from normal and bleached M. alcicornis specimens

The phospholipase A2 (PLA2) activity of the soluble proteomes from normal and bleached M. alcicornis specimens was determined using a secretory PLA2 colorimetric assay kit (Cayman Chemical, Ann Arbor, MI, USA). The assay uses the 1,2-dithio analog of di-heptanoyl phosphatidylcholine as substrate. Free thiols generated upon hydrolysis of the thioester bond at the sn-2 position by PLA2 were detected using DTNB (5,5′-dithio-bis-(2- nitrobenzoic acid)). Color changes were monitored by a Benchmark Plus microplate spectrophotometer (Bio-Rad Laboratories, Hercules, CA, USA) at 414 nm, sampling every minute for 10 min. The PLA2 activity was expressed as micromoles of hydrolyzed substrate per minute per milligram of protein. The PLA2 activity was calculated using GraphPad Prism 6.0 (GraphPad Software, San Diego, CA, USA).

Comparative proteolytic activity of the soluble proteomes from normal and bleached M. alcicornis specimens

The proteolytic activity of the soluble proteomes from normal and bleached M. alcicornis specimens was determined according to a method previously described with some modifications (Nagaraju et al., 2006). Briefly, soluble proteomes (50 μg protein) were incubated separately with one mL of substrate (2% casein in 50 mM Tris, 100 mM NaCl, five mM CaCl2, pH 8.8) for 2 h and 30 min at 37 °C. Undigested casein was precipitated by adding 1.5 mL of 0.44M trichloroacetic acid. The digested casein in the supernatant (one mL) was determined using the Folin–Ciocalteu’s reagent. One unit of activity was defined as the amount of enzyme required to cause an increase in the optical density by 0.01 at 660 nm/min. Data obtained from the proteolytic assay were plotted using GraphPad Prism 6.0 (GraphPad Software, San Diego, CA, USA).

Data analysis and statistics

Results are expressed as mean ± S.E.M. from n = 3 experiments. In the case of the hemolytic activity, concentration-response curves (CRC) were repeated three times employing different rats. CRC were plotted and fitted to the Boltzmann equation, using the data analysis and graphics software GraphPad Prism 6.0 (GraphPad Software, San Diego, CA, USA). Means ± S.E.M. obtained from cell density, chlorophyll content, hemolytic, PLA2, and caseinolytic assays were compared using an unpaired Student’s t-test. In all cases, statistical significance is indicated by p < 0.05.

Results

Normal (Fig. 1A) and visibly discolored (Fig. 1B) fragments of the “fire coral” M. alcicornis were collected in the Mexican Caribbean in November 2016. Bleaching caused significant decreases in chlorophyll a and chlorophyll c2 content (Fig. 1C), and cell density of symbionts (Fig. 1D).

Normal (A) and bleached (B) fragments from M. alcicornis.

Figure 1: Normal (A) and bleached (B) fragments from M. alcicornis.

(C) Chlorophyll a (Chl a) and chlorophyll c (Chl c2) content per cm2. (D) Symbiont´s density quantified for normal and bleached specimens. Data are mean ± SEM. *Indicate significant difference. Photographs from normal and bleached specimens were taken by Norma Olguín-Lopez and Víctor Hugo Hernández-Elizárraga.

Protein content in lyophilized soluble proteomes from normal and bleached specimens were 21.12 and 13.96 μg of protein per mg of lyophilizate, respectively. The SDS–PAGE analysis revealed that proteomes from normal and bleached hydrocorals have different protein profiles and contained proteins with a broad range of molecular weights, ranging in size from 10 to 200 kDa (Fig. 2).

Electrophoresis profiles of the soluble proteomes from normal (A) protein molecular weight marker (B) and bleached (C) specimens of M. alcicornis.

Figure 2: Electrophoresis profiles of the soluble proteomes from normal (A) protein molecular weight marker (B) and bleached (C) specimens of M. alcicornis.

Samples (80 μg of protein) were separated by one dimensional SDS–PAGE using 18% (w/v) polyacrylamide under non-reducing conditions. Protein bands were stained using Coomassie Blue R-250. Arrows indicate protein bands in which differences were observed.

2DE-PAGE analysis allowed the detection of 52 and 75 protein spots in the proteomes of normal (Fig. 3A) and bleached hydrocorals (Fig. 3B), respectively. 2DE-PAGE gels showed well-resolved protein spots with molecular masses between 10 and 50 kDa and pI values ranging from 4 to 8. A total of 17 proteins were differentially expressed in bleached hydrocorals with respect to normal specimens, six of these proteins were up-regulated after bleaching (Fig. 3A), while eleven were down-regulated (Fig. 3B). Analysis by MALDI-TOF/TOF, the ProteinPilot software, and BLASTp was carried out with eight of the differentially expressed proteins (Table 1). The other protein spots were not examined due to their low concentration.

Two-dimensional gels of the soluble proteomes from normal (A) and bleached (B) M. alcicornis specimens.

Figure 3: Two-dimensional gels of the soluble proteomes from normal (A) and bleached (B) M. alcicornis specimens.

Up-regulated proteins are marked with green circles and down-regulated proteins are marked with red circles. Protein spots with differential expression were selected based on a fold change ≥2.
Table 1:
M. alcicornis proteins identified by MALDI-TOF/TOF, the ProteinPilot software, and BLASTp analysis.
Spot Protein name Accession no. Peptides No of peptides matched m/z Mw/pI Sequence coverage
201 Calmodulin XP_020603574.1 EAFSLFDK
DGNGFISAAELR
VFDKDGNGFISAAELR
EAFSLFDKDGDGTITTK
ADQLTEEQIAEFK
EADIDGDGQVNYEEFVK
6 956.5549
1,249.6982
1,738.9722
1,845.0011
1,563.8358
1,927.9585
16.89/3 24.5%
1002 Disintegrin and metalloproteinase domain-containing protein 7 KII70418.1 TLGENGKEFSK
EFLLIFR
2 1,223.7189
937.6454
10.47/4.02 7.5%
3402 DELTA-actitoxin-Ate1a-like XP_020897118.1 FLQTEMDK
ELGVDYFLERR
EAHLILFLVNAR
RIKLKYAHLGSSNPP
4 1,028.5859
1,378.8389
1,378.8082
1,713.0261
18.55/4.39 22.7%
4102 Phospholipase A2 XP_004208832.1 TLGLPTAPYVEALSFAR
QDDERANYAK
2 1,806.0854
1,211.6534
11.61/4.99 19.8%
7102 ATP-dependent Zinc metalloprotease KII60716.1 NRARALINIAHPK
QFISSLEDAETIKSFNNLR
KNAGKAINAHLIPK
QKAENILSFLQK
ALLLLKDFDSVPPVVR
LSVNTPNNRQPLTITVSTKK
6 1,474.918
2,212.2581
1,474.922
1,418.8232
1,799.0536
2,211.2761
12.56/5.98 68.9%
7805 Component 3 exocytosis complex KII70395.1 NFAKLSQSLETLDSQLNLR
TSGDKGKTSPR
2 2,211.3179
1,182.6781
70.09/6.15 3.9%
5301 Actin XP_020600429.1 SYELPDGKVITIGNER 1 1,792.0223 32.02/5.41 5.6%
3203 Collagen alpha XP_020603217.1 GSAGPPGATGFPGAAGR 1 1,459.8524 14.41/4.83 4.3%
DOI: 10.7717/peerj.6593/table-1

Note:

The protein spots marked in red showed down-regulation in bleached Millepora alcicornis, while the protein spots marked in green showed up-regulation expression.

Our results showed that calmodulin (CaM), actin, and two potential toxins (DELTA-actitoxin-Ate1a-like and a PLA2) were up-regulated after bleaching, whereas component 3 exocytosis complex, collagen alpha, a disintegrin metalloproteinase domain-containing protein 7, and an ATP-dependent zinc metalloprotease were down-regulated. All proteins were identified on the basis of manual annotations using BLASTP to identify proteins with homology to Cnidarian proteins described in the NCBI database.

On the other hand, the soluble proteomes from normal and bleached specimens elicited a concentration-dependent hemolysis. In addition, both soluble proteomes displayed PLA2 and caseinolytic activities. No statistically significant difference was observed between the hemolytic and enzymatic activities produced by the proteome from normal specimens and those induced by the proteome from bleached hydrocorals (Table 2).

Table 2:
Protein content, hemolytic, PLA2, and caseinolytic activities of the soluble proteomes from normal and bleached M. alcicornis specimens.
Extracts Protein content (μg protein/mg lyophilized) Hemolytic activity (Hu50, μg protein/ml) PLA2 activity (μmol/min/mg) Caseinolytic activity (U/min/μg protein)
Normal 21.12 ± 1.1 0.1343 ± 0.02 41.37 ± 2.1 27.57 ± 1.43
Bleached 13.96 ± 1.4* 0.09061 ± 0.005 49.11 ± 0.9 30.82 ± 2.05
DOI: 10.7717/peerj.6593/table-2

Notes:

Each value is expressed as the mean ± S.E.M. (n = 3).

p < 0.05 compared to control.

PLA2, phospholipase A2.

Discussion

In recent decades, increases in the frequency and intensity of coral bleaching events have resulted in declines in coral cover worldwide (Roff, Zhao & Mumby, 2015; Lough, 2016). Accordingly, these damaged ecosystems have suffered a dramatic “phase shift,” in which habitats dominated by reef-forming organisms are now dominated by macroalgae and low-relief corals (Jackson et al., 2014; Cramer et al., 2017). In fact, coral reefs are considered the most degraded ecosystems on the planet due to climate change (Eakin et al., 2010; Roff, Zhao & Mumby, 2015; Cramer et al., 2017). During 2014–2017, record high temperatures triggered the third global scale event of coral bleaching ever registered (Eakin et al., 2017; Hughes et al., 2017b, 2018). It is estimated that 30% of the reef areas that were monitored after the ENSO 2015–2016 showed severe bleaching (Hughes et al., 2018), and the Caribbean coral reefs were seriously affected. Hydrocorals of the genus Millepora spp. are among the organisms most affected by bleaching (Dias & Gondim, 2016). Particularly, M. alcicornis has been considered very vulnerable to temperature rise, since it has suffered severe bleaching episodes in different times and locations, such as in the Great Barrier Reef (Marshall & Baird, 2000), where an 85% mortality of this species was observed during the summer of 2005; in the Florida Keys reef in the summers of 2006 and 2007 (Wagner, Kramer & Van Woesik, 2010), and in Puerto Rico and the Caribbean during 1987, 1993, 1995, 1998, 2003, 2005, 2009, and 2010 (Dias & Gondim, 2016).

In the present study, normal and visibly bleached fragments of the “fire coral” M. alcicornis, (Figs. 1A and 1B) were collected in the Parque Nacional Arrecife de Puerto Morelos (Quintana Roo, México) in November 2016, the warmest year in NOAA’s 137-year series (NOAA, 2017). Since bleaching is provoked by a decrease in symbionts population and a reduction in the concentration of their photosynthetic pigments (Hoegh-Guldberg, 1999; Lesser, 2006, 2011; Obura, 2009), chlorophyll contents, and symbiont’s density were determined in M. alcicornis. As was expected, chlorophyll a and chlorophyll c2 concentrations per cm2 of hydrocoral were significantly lower in the bleached specimens (Fig. 1C), which indicates that warming provoked breakdown of the hydrocoral-algae symbiosis due to photoinhibition, in a similar way as that observed in Anthozoan species (Warner, Fitt & Schmidt, 1999; Takahashi & Murata, 2008; Lesser, 2011). Additionally, bleached M. alcicornis specimens showed a decrease of 40% in the density of symbionts per square centimeter (Fig. 1D), which corresponds to a “moderate bleaching” according to the categories used in the ReefBase (Oliver, Berkelmans & Eakin, 2009). Unlike what has been previously observed, the M. alcicornis colonies examined in the present study after the bleaching event of 2015–2016, did not show severe bleaching. This could indicate that this Caribbean hydrocoral is developing a certain degree of thermotolerance, which might be associated with the switching of symbionts, as reported in other species (Silverstein, Cunning & Baker, 2015; Bay et al., 2016; Swain et al., 2018). However, this hypothesis should be proven in future studies.

Comparative proteomic studies have been used to study the effect of heat stress on the total proteome from some Anthozoan species that include: Acropora palmata, Acropora microphthalma, Pocillopora acuta, Seriatopora hystrix, and Aiptasia pallida. The results derived from such investigations have shown that thermal stress, which is responsible for coral bleaching, provokes the differential expression of several proteins involved in different processes essential for the survival of the cnidarians, including cytoskeleton organization, thermal and UV stress response, redox state, immunity, calcium homeostasis, transcription factors, exocytosis, and metabolism (Weston et al., 2015; Ricaurte et al., 2016; Oakley et al., 2017; Mayfield et al., 2018a, 2018b). However, to date there are no reports about the effect of breakdown of the cnidarian-Symbiodinium symbiosis due to thermal stress on the expression of proteins, including toxins, in an organism of the class Hydrozoa. Therefore, in this study the differential expression of proteins in the soluble proteomes from normal and bleached M. alcicornis specimens was analyzed.

The extraction of high-quality proteins from corals with calcareous exoskeleton is not an easy task (Cheng et al., 2018). The most common methods reported for protein extraction include the use of trizol (Garcia et al., 2016; Mayfield et al., 2018a), glass bead-assisted extraction (Weston et al., 2015), and sonication-assisted extraction with rehydration buffer (Ricaurte et al., 2016). However, the method chosen for the extraction of M. alcicornis proteins involved osmotic shock in bidistilled water. This method causes the discharge of the nematocysts content (Ibarra-Alvarado et al., 2007; García-Arredondo et al., 2015, 2015, 2016; Morabito et al., 2012; Hernández-Matehuala et al., 2015).

Sodium dodecyl sulfate–polyacrylamide gel electrophoresis evidenced that soluble proteomes from normal and bleached M. alcicornis contain proteins with a broad range of molecular weights (Fig. 2). Most of the bands detected in both electrophoretic profiles are within the range of 10–200 kDa. A similar electrophoretic pattern was observed for the proteome of isolated nematocysts from Hydra magnipapillata, (a well-known member of the Hydrozoa class), in which most of the protein bands are in the range of 10–70 kDa (Balasubramanian et al., 2012). Some differences were observed in the SDS–PAGE profiles obtained from normal and bleached hydrocorals, indicating a differential protein expression between both types of specimens.

Moreover, 2DE-PAGE showed that the soluble proteomes from normal and bleached specimens of M. alcicornis resolved 58 and 75 proteins, respectively. It is important to mention that in the case of the present investigation, we examined the effect of elevated seawater temperature on the soluble proteome of M. alcicornis, unlike what has been done in previous studies carried out on Anthozoan species, in which the impact of thermal stress on the total proteomes from Acropora palmata, A. microphthalma, Pocillopora acuta, Seriatopora hystrix, and Aiptasia pallida was analyzed (Weston et al., 2015; Ricaurte et al., 2016; Oakley et al., 2017; Mayfield et al., 2018a, 2018b). After bleaching, the differential expression of 17 proteins was observed in M. alcicornis, six proteins were up-regulated, while 11 were down-regulated.

Calmodulin, actin, and collagen are some of the proteins differentially expressed in bleached specimens of M. alcicornis. Previous studies have reported the differential expression of these same proteins, which are involved in calcium homeostasis, cytoskeleton, and extracellular matrix (ECM) in Anthozoan species subjected to bleaching (Bay & Palumbi, 2015; Pinzón et al., 2015; Raina et al., 2015; Seneca & Palumbi, 2015; Weston et al., 2015; Maor-Landaw & Levy, 2016; Oakley et al., 2016, 2017; Ricaurte et al., 2016; Sloan & Sawyer, 2016; Huang et al., 2017; Ruiz-Jones & Palumbi, 2017; Traylor-Knowles et al., 2017; Mayfield et al., 2018a).

Calmodulin is a key Ca2+ sensor, whose signaling is important in numerous cellular processes, such as cell cycle and calcium homeostasis (Haeseleer et al., 2002; Desalvo et al., 2008; Reyes-Bermudez, Miller & Sprungala, 2012). Several studies have demonstrated that expression of genes and proteins involved in Ca2+ homeostasis is modified after bleaching in Acropora palmata, Orbicella faveolata, Acropora millepora, and Galaxea astreata (Desalvo et al., 2008, 2010; Rodriguez-Lanetty, Harii & Hoegh-Guldberg, 2009; Moya et al., 2012; Ricaurte et al., 2016; Huang et al., 2018). In the case of M. alcicornis, CaM was up-regulated after bleaching, which implies that changes are occurring in biomineralization and other calcium-dependent processes related to the development of the hydrocorals. It has been shown that CaM overexpression and a rupture in Ca2+ homeostasis are linked to an abnormal reorganization of the actin cytoskeleton (Desalvo et al., 2008; Reyes-Bermudez, Miller & Sprungala, 2012).

Actin was another protein that was up-regulated in bleached M. alcicornis specimens, in a similar way to what was found in scleractinian corals such as Acropora palmata, Stylophora pistillata, and Montastraea faveolata (Desalvo et al., 2008; Kenkel, Meyer & Matz, 2013; Maor-Landaw et al., 2014; Ricaurte et al., 2016; Louis et al., 2017). Actin is implicated in the construction of filaments and supports the majority of motile events in eukaryotic cells (Fletcher & Mullins, 2010). The increase in the expression of actin observed in bleached hydrocorals might be related to the relocation of symbionts as a mechanism of photoprotection (Petrou, Ralph & Nielsen, 2017) and to the increase of heterotrophic feeding to mitigate the energetic imbalance by the departure of symbionts (Wooldridge, 2014; Tremblay et al., 2016), which requires a significant contribution by the cytoskeleton to collect and transport nutrients.

Among the proteins that showed a down-regulation in bleached hydrocorals was collagen. This result is in agreement with what was found in specimens of Acropora palmata, Stylophora pistillata, and Anemonia viridis that were subjected to bleaching conditions (Moya et al., 2012; Maor-Landaw & Levy, 2016; Ricaurte et al., 2016). Collagen, abundantly found in the calicoblastic space, is central during the biomineralization process since it provides the structural support to which proteins of the subfamily CARP and analog proteins bind in nucleation and mineral growth sites (Mohamed et al., 2016). Moreover, a decline in symbiont´s density has been related to changes in the active volume of calcification due to a reduction in dissolved inorganic carbon (Drake et al., 2013; D’Olivo & McCulloch, 2017). Thus, down-regulation of collagen and reduced symbiont´s density suggest that thermal stress induces a decrease in the rate of calcification in the hydrocoral M. alcicornis.

The exocytosis multiprotein complex provides spatial targeting of exocytotic vesicles to the plasma membrane (Picco et al., 2017), and therefore is related to the process of symbionts expulsion (Bieri et al., 2016). In fact, the exocytosis complex component 4 was up-regulated in Porites astreoides, and has been proposed as a biomarker of coral thermal stress (Kenkel, Meyer & Matz, 2013; Kenkel et al., 2014). In contrast, exocytosis complex component 3 was down-regulated in bleached M. alcicornis, as previously reported for other Endo-exo phagocytosis-related proteins in Acropora palmata (Ricaurte et al., 2016). Thus, the proposal that up-regulation of proteins that are part of the exocytosis complex represents a biomarker of thermal stress is still controversial.

Three of the soluble proteins from M. alcicornis, which were differentially expressed in bleached hydrocorals, showed amino acid sequence similarity to potential toxins. Despite the ecological relevance of reef forming hydrocorals, the information regarding the chemical structure and mechanism of action of their toxins is scarce. Previous studies carried out by our research group showed that the aqueous extract of M. alcicornis caused hemolysis of rat erythrocytes and showed PLA2 activity (Hernández-Matehuala et al., 2015). Zymographic analysis of this extract revealed that it contained ∼28–30 kDa cytolysins with PLA2 activity; ∼200 kDa cytoloysins, which do not have PLA2 activity (Hernández-Matehuala et al., 2015) and proteins with proteolytic activity at ∼25, ∼40, and 80–200 kDa (Olguín-López Norma, Hernández-Elizarraga Víctor, Hernández-Matehuala Rosalina, Cruz-Hernández Andrés, Guevara-González Ramón, Caballero-Perez Juan, Ibarra-Alvarado César , and Rojas-Molina Alejandra, 2018, unpublished data). Interestingly, we found that bleaching induced up-regulation of proteins that showed sequence homology with DELTA-actitoxin-Ate1a-like protein from the sea anemone Exaiptasia pallida (Sunagawa et al., 2009) and a putative PLA2 from Hydra magnipapillata (Sher et al., 2005). In contrast, bleaching elicited down-regulation of a protein that has sequence similarity with a ATP-dependent Zinc metalloprotease from Aiptasia pallida (Sunagawa et al., 2009).

Cnidarian toxins are mainly proteins and peptides (Anderluh & Maček, 2002; Nevalainen et al., 2004; Mariottini et al., 2015; Podobnik & Anderluh, 2017) and can be classified into three categories based on their mechanism of action: enzymes (PLA2 and metalloproteinases); pore-forming proteins (actinoporins, jellyfish toxins, hydralysins-related toxins, and membrane attack complex-perforins), and neurotoxins (Jouiaei et al., 2015b). The most studied pore forming toxins of cnidarians are actinoporins, that have been found in many species of marine anemones (Anderluh & Maček, 2002; Kristan et al., 2009; Šuput, 2009). One of the proteins whose expression was up-regulated after bleaching displayed sequence similarity with actinoporin DELTA-actitoxin-Ate1a-like protein. This finding suggests that M. alcicornis produce PFTs, whose structure might be related to that of the anemone actinoporins, which is interesting considering that this type of toxins have been found in only one species of the class Hydrozoa, Hydra magnipapillata (Jouiaei et al., 2015b).

On the other hand, the expression of a protein that exhibited sequence similarity to a putative PLA2 of 15.7 kDa from Hydra magnipapillata, which was found by a large-scale search at the Hydra magnipapillata EST Database, was up-regulated after bleaching. This Hydra PLA2 showed homology with PLA2s from Apis mellifera, the lizard Heloderma uspectum (“Gila Monster”), and the scorpions Mesobuthus tamulus and Pandinus imperator, mainly in the calcium binding domain and in the catalytic site. All of these PLA2s display neurotoxic activity and belong to the secreted PLA2 (sPLA2) group III (sPLA2-III) (Sher et al., 2005). In contrast to snake venom sPLA2s (Costa, Camargo & Antunes, 2015; Kordiš & Križaj, 2017; Tonello & Rigoni, 2017), cnidarian sPLA2s have been poorly characterized and their enzymatic action greatly differs between species (Nevalainen et al., 2004).

The present study also evidenced that proteins that displayed sequence similarity to disintegrin and metalloproteinase domain-containing protein 7 and an ATP-dependent Zinc metalloprotease showed down-regulation after bleaching. Metalloproteinases have been detected in some cnidarian nematocyst venoms, such as those of Stomolophus meleagris (Li et al., 2014), Olindias sambaquiensis (Knittel et al., 2016), Pelagia noctiluca (Frazão et al., 2017), and Chironex fleckeri (Jouiaei et al., 2015a). This class of toxins shows a great diversity of effects. Snake venom metalloproteinases (SVMPs) affect the ECM in various ways, including the release of ECM-derived biologically-active peptides, which exerts either reparing or damaging effects on tissues (Gutiérrez et al., 2016). These enzymes are also involved in the proteolytic processing of other venom proteins (Da Silveira et al., 2007) and are capable of causing severe inflammation by disrupting capillary vessels and tissues (Weston et al., 2013). However, the different kinds of cnidarian metalloproteinases and their effects are not yet clear. It has been demonstrated that SVMPs undergo various post-translational modifications, which contribute to their great functional diversity (Moura-da-Silva et al., 2016). This could be the case of cnidarian metalloproteinases. If this is so, then it could be hypothesized that the imbalance of energy produced by symbionts expulsion or death induce cnidarians to optimize their energy expenditure, which impairs the metalloproteases synthesis that represents a great energy investment.

Cnidarians employ toxins as a defense against predators and in prey capture to obtain heterotrophic nutrition (Anderluh & Maček, 2002; Anderluh et al., 2011; Houlbrèque, Rodolfo-Metalpa & Ferrier-Pagès, 2015; Ben-Ari, Paz & Sher, 2018). Particularly, both, PFTs and PLA2s produce cytolysis and play a very important role in defensive and offensive actions, since they provoke the destruction of cell membranes (Rojko et al., 2016; Macrander & Daly, 2016; Podobnik & Anderluh, 2017). It has been documented that certain corals, such as Oculina arbusculata, can increase their heterotrophic nutrition to maintain and restore energy reserves after a bleaching event (Grottoli, Rodrigues & Palardy, 2006; Aichelman et al., 2016). In the case of the present study, the differential expression of potential toxins suggests that thermal stress alters heterotrophic competences to mitigate the energy imbalance in hydrocorals. However, heat stress did not affect the hemolytic, PLA2, and proteolytic activity of M. alcicornis soluble proteome, which suggests that the cytolysis induced by the toxins of this hydrocoral is mainly produced by enzymes and PFTs (Nevalainen et al., 2004; Lee et al., 2011; Anderluh et al., 2011; Jouiaei et al., 2015a; Rojko et al., 2016; Knittel et al., 2016), whose expression was not affected after suffering moderate bleaching. It is worth mentioning that in a previous study carried out on M. alcicornis specimens collected in the Mexican Caribbean in 2008, we found that thermal stress did not significantly modify the PLA2 activity of the aqueous extract prepared from bleached M. alcicornis (García-Arredondo et al., 2015), which agrees with what we found in the present study with specimens collected in 2016. These findings seem to indicate that M. alcicornis subjected to moderate bleaching does not lose its capability to synthesize cytolysins. Evidently, it is important to continue investigating the impact of elevated ocean temperatures that provoke bleaching on the synthesis of toxins produced by reef forming cnidarians and the significance of heterotrophic feeding as a mechanism to counteract the deleterious effects of this phenomenon.

Conclusions

The present study represents the first report of the effect of thermal stress on the proteomic profile of a reef forming cnidarian of the class Hydrozoa. The decrease in symbiont’s density and chlorophyll a and c2 levels suggests that ENSO 2015–2016 induced a moderate bleaching in colonies of M. alcicornis that inhabit the Mexican Caribbean. Proteins involved in various important cellular processes, such as calcium homeostasis, exocytosis, and cytoskeleton organization were differentially expressed in bleached hydrocorals. Three of the differentially expressed proteins showed amino acid sequence similarity to potential toxins, however, the hemolytic, PLA2, and proteolytic activity of bleached M. alcicornis soluble proteome was not modified. These results suggests that the cytolytic effect induced by this hydrocoral is produced by toxins whose synthesis is not altered after bleaching, which allow us to hypothesize that M. alcicornis is able to prey heterotrophically when suffering from moderate bleaching, giving it a better chance to face the effects of high temperature.

Supplemental Information

13 Citations   Views   Downloads