Layout 1 INTRODUCTION Marine diatoms are some of the most successful groups of unicellular eukaryotic photosynthetic organ- isms, forming massive blooms in coastal and open waters and having a key role in global marine primary production (Armbrust, 2009). Elucidating the mechanisms that regu- late diatom bloom dynamics and mainly their decay de- terminants is therefore fundamental to understand nutrient and carbon cycling so as energy flow at the global as well as the local scale. Programmed cell death (PCD) is an active, self-con- trolled death mechanism well known in higher organisms which has been also detected in unicellulars (Bidle et al., 2007; Segovia, 2011). However, its occurrence in unicel- lular organisms is still considered a matter of ongoing de- bate (Deponte, 2008; Jimenez et al., 2009), and the specific mode of action in which PCD is initiated remains to be elu- cidated. PCD has been suggested to act in phytoplankton populations to eliminate aged cells and to reduce population size in order to cope with nutrient limitation, for instance, in the final stages of blooms (Bidle and Falkowski, 2004; Vardi et al., 1999, 2006, 2007). After a certain threshold of limitation of environmental factors (mainly nutrients), PCD induces the massive decay of the bloom, sometimes through massive production of resting or resistance stages (Ellegaard et al., 2013; McQuoid et al., 2002). These fac- tors include nitrogen and light limitation (Berges and Falkowski, 1998), iron starvation (Bidle and Bender, 2008; Thamatrakoln et al., 2012), virus infection (Bidle et al., 2007; Vardi et al., 2009), CO2 limitation (Vardi et al., 1999) and several secondary metabolites (Cervia et al., 2009; Costas et al., 1993). Reactive oxygen species (ROS) are often involved in inducing PCD in reaction to a stress factor in a wide range of organisms and cell lines, from bacteria to mammalians and plants (Andrianasolo et al., 2007; Jones, 2008; Kuwabara et al., 2008; Lam, 2008; Rocken- feller and Madeo, 2008; Scherz-Shouval and Elazar, 2007), mainly at intermediate concentrations. In fact, several cell reactions to ROS strongly depend upon the ROS concen- tration and may range from production of antioxidant com- pounds to necrotic cell death (Chandra et al., 2000), or PCD (Costas et al., 1993; Vardi et al., 2009). Apart from ROS concentrations, several variables come into play by modu- lating and regulating the cell response to ROS. These in- clude the cell type, the stress type, its magnitude and its duration in time (Martindale and Holbrook, 2002). Different case studies report the co-occurrence of ROS and PCD in marine phytoplankton and suggest a correlation between them (Butow et al., 1997; Segovia and Berges, 2009; Vardi et al., 1999). In the dinoflagel- late Peridinium gatunense, PCD leading to the collapse Advances in Oceanography and Limnology, 2015; 6(1/2): 13-20 ORIGINAL ARTICLE DOI: 10.4081/aiol.2015.5466 Expression of death-related genes and reactive oxygen species production in Skeletonema tropicum upon exposure to the polyunsaturated aldehyde octadienal Alessandra A. Gallina,1# Chih-Ching Chung,2,3 Raffaella Casotti1* 1Stazione Zoologica Anton Dohrn di Napoli, Villa Comunale, 80121 Naples, Italy; 2Institute of Marine Environmental Chemistry and Ecology, National Taiwan Ocean University, Keelung 20224, Taiwan, China; 3Center of Excellence for the Oceans, National Taiwan Ocean University, Keelung 20224,Taiwan, Republic of China #Present address: Department of Biology, Colorado State University, Fort Collins, CO 80523, USA *Corresponding author: raffaella.casotti@szn.it ABSTRACT The effects of 4E/Z-octadienal (OCTA) on ScDSP-1 and ScDSP-2 gene expression and reactive oxygen species (ROS) production were investigated in the marine diatom Skeletonema tropicum (formerly costatum) using qRTPCR and flow cytometry. ScDSP-1 and ScDSP-2 genes have been previously shown to be involved in cell death in ageing cells and in response to photosynthetic stress. OCTA in- duced a differential, concentration-dependent DSP gene expression associated to ROS production, 821.6 and 97.7 folds higher for ScDSP- 1 and ScDSP-2, respectively. Among the concentrations tested, only 8 μM OCTA, which caused a reduction of 50% in cell concentrations at 24 h, was able to elicit an expression pattern consistent with a signalling role. Interestingly, only intermediate levels of reactive oxygen species (ROS) (i.e., 1.5±0.1 increase) were observed to be elicited by such concentration. These results suggest that ROS are key components of the molecular cascade triggered by polyunsaturated aldehydes (PUA) and leading to cell death. This could have implications for bloom final stages at sea, where PUA may act as effectors of diatom population dynamics through ROS acting as modulators. Key words: Diatom; chemical ecology; bloom evolution; Programmed cell death; ROS. Received: July 2015. Accepted: October 2015. No n- co mm er cia l u se on ly 14DSP expression in response to PUA of the algal population is mediated by oxidative stress in- duced by carbon dioxide (CO2) limitation (Vardi et al., 1999). Similarly, ROS-mediated cell death was observed at the end of a bloom of the same species, despite an in- crease in antioxidant activity (Butow et al., 1997). PCD- like cell death concomitant with accumulation of ROS is also reported to be induced by the toxic cyanobacterium Microcystis sp. (Vardi et al., 2002), or H2O2 (Vardi et al., 2007). In the chlorophyte Dunaliella tertiolecta, dark-in- duced PCD, involving the activity of caspase-like en- zymes has been linked to an increase in ROS production (Segovia and Berges, 2009). In the prymnesiophyte Emil- iana huxley viral infection has been associated to an en- hanced production of ROS, possibly leading to PCD (Evans et al., 2006). Interestingly, the secondary metabolite Euplotin C de- rived from the marine ciliated protist Euplotes crassus, is found to induce PCD in the congeneric Euplotes vannus, which does not produce this metabolite (Cervia et al., 2009). This suggests that specific secondary metabolites can play an ecological role in broadening phytoplankton niche size through different mechanisms, including PCD (Cervia et al., 2009). In the marine diatom Phaeodactylum tricornutum, the exposure to diatom-derived secondary metabolites, polyunsaturated aldehydes (PUA), results in an altered ex- pression of metacaspases through the overexpression of a gene associated to NO generation (Vardi et al., 2008). In the marine diatom Thalassiosira weissflogii, PUA were shown to induce PCD-like cell death, suggesting for a role in activating the cell death cascade, therefore regulating population dynamics and species succession (Casotti et al., 2005). In Phaeodactylum tricornutum it is speculated that a sophisticated stress surveillance system exists, in- volving NO and Calcium, in response to PUA, which could lead to bloom termination and population-level cell death at specific PUA concentrations (Vardi et al., 2006). Even PUA precursors, the polyunsaturated fatty acids (PUFAs), have been reported to inhibit growth of Candida albicans and Candida dubliniensi by inducing ROS ac- cumulation and consequent PCD (Thibane et al., 2012). In the marine diatom Skeletonema tropicum (formerly costatum) two genes have been identified, coding for two identical Ca-regulated proteins, named S. costatum Death- Specific Proteins (ScDSP-1 and ScDSP-2) (Chung et al., 2005), which are both composed of a transmembrane do- main and a pair of EF-hand motifs (Luan et al., 2002; Means and Dedman, 1980). Both ScDSP gene expressions were shown to be strongly upregulated in aging cells and also during light stress. Their expression is correlated to DNA fragmentation, suggesting for a possible role of the encoded protein in the signal transduction of stress to the cell death machinery (Chung et al., 2005, 2008). In partic- ular, ScDSP expression has been related to the blockage of the electron flow between PSII and cytochrome b6f and it was found to be NO dependent (Chung et al., 2008). More- over, two DSP-like proteins (TpDSP1 and TpDSP2), have been recently identified in the marine diatom Thalassiosira pseudonana and have been shown to be upregulated both under iron limitation and starvation, suggesting for a role in Fe-induced stress and consequent death (Thamatrakoln et al., 2012). However, they have also recently been shown to enhance growth under Fe limitation and to be controlled by light intensity (Thamatrakoln et al., 2013). In this study we report the gene expression of the two death-related genes ScDSP-1 and ScDSP-2 in ST exposed to the polyunsaturated aldehyde octadienal, and relate it to ROS production. The overall aim is to relate cell death to chemical stress induced by PUA and to advance hy- potheses on bloom decay mechanisms at sea. METHODS Experimental setup and cultures A unialgal culture of Skeletonema tropicum Cleve (for- merly costatum) was grown at 23°C on a 12h-12h light- dark cycle under a photon flux density of 110 µmol quanta m–2 s–1 (Chung et al., 2005). Natural seawater, amended with f/2 nutrients (Guillard, 1975) was used as medium. For gene expression measurements, S. tropicum was grown in 10 L polycarbonate carboys (Nalgene). When the culture reached the late exponential phase of growth (LE, Fig. 1), cells were transferred into 2 L polycarbonate bot- tles (Nalgene) and treated with different concentrations of 4E/Z-octadienal (OCTA) or H2O2. Samples were collected at 1,2,3,4,5 and 1,2,3,4,12 days since the beginning of the Fig. 1. Effect of OCTA and H2O2 on cell concentrations of S. tropicum. LE is the time point indicated as Late Exponential growth phase. Data are means of triplicates ±SD. When the bar is not visible, it is smaller than the symbol. No n- co mm er cia l u se on ly 15 A.A. Gallina et al. experiment for cell counts and total RNA extraction, re- spectively. Cell concentrations were determined by using a Sedgwick-Rafter counting chamber (Hausser Scientific Partnership). At least 200 cells from each sample were counted using a light microscope (BX60; Olympus, Center Valley, USA) at 200× magnification. Growth rates were calculated using the following equation: µ=ln [(N1/N0)/t] (eq. 1) where N0 and N1 are cell density at the start and the end of the growth period, and t is the time between measure- ments (in days). Since no dead vs alive cells assay has been used, it is possible that the cells counts reported in the study might include damaged or dying cells. Gene expression experiments were performed using three OCTA concentrations, namely the concentration which caused a reduction of 50% in cell concentration at 24 h, half this concentration and twice this concentration (4, 8 and 16 µM). For ROS detection, three independent cultures were grown in 1 L glass flasks until the late ex- ponential growth phase. Cells were harvested and directly exposed to the same concentrations of OCTA used for gene expression measurements and analysed through flow cytometry as described below. The controls were dye- loaded samples processed exactly as the other samples ex- cept for PUA inoculation. Each experiment was performed at least twice, but generally three times, i.e. the data presented are the results of biological replicates. Chemicals The PUA used in this study was 2E,4E/Z-octadienal (OCTA, from Sigma Aldrich Inc. Milan, Italy). This was chosen because it is the main and most commonly pro- duced PUA by diatoms, and, in particular, by Skeletonema spp. (Wichard et al., 2005). The concentrations used were chosen accounting for multiples of the concentration caus- ing a reduction of 50% in cell concentrations at 24h, as empirically determined for S. tropicum (data not shown) resulting to be 8 μM. The reference value of 8 μM has been chosen also based on Ribalet et al., (2007), where 8 μM represented the EC50 for growth in the congeneric species S. marinoi. OCTA working solutions were dis- solved in HPLC-grade pure Methanol (JT Baker, Phillips- burg, NJ, USA) at room temperature. Since Methanol can be toxic, it was tested on growing cultures and cultures showed no reaction up to 7 microliters per mL. Subse- quently, and conservatively, no more than 4 microliters of Methanol plus OCTA were ever inoculated into the cul- tures. Before incubation, the exact OCTA concentration in the working solution was assessed spectrophotometri- cally (Hewlett-Packard 8453- Hewlett-Packard Co., Palo Alto, CA, USA) by measuring absorption at 274 nm and a specific molar absorption coefficient of 31,000 (Pippen and Nonaka, 1958). Hydrogen peroxide (H2O2, 30 wt. % from Sigma Aldrich, Milan, Italy) was used as positive control, after dilution 1/200 in filtered seawater and added to samples as a 50 μL:1 mL ratio. The ROS scavenger, 4-Hydroxy-2,2,6,6-tetram- ethylpiperidine 1-oxyl (TEMPOL, Sigma-Aldrich, Milan, Italy) was used as negative control, so as sodium di- ethyldithiocarbamate trihydrate (DETC, Sigma-Aldrich), an inhibitor of superoxide dismutase (SOD) at a final con- centration of 5 mM and 1 mM, respectively, in filtered sea water (FSW). Samples with DETC were incubated for 120 min prior to PUA addition. The ROS-sensitive dye dihydrorodamine 123 (DHR, Molecular Probes, Leiden, NL) (5 mM stock in DMSO) was used to monitor ROS production at the single-cell level. DHR can be oxidized by different ROS, including hydrogen peroxide (H2O2) and peroxynitrite (ONOO–), to form the fluorescent de- rivative rhodamine 123 (excitation WL 492-495 nm, emission WL=517-527 nm), and is commonly used to in- vestigate oxidative stress in vivo (Ischiropoulos et al., 1999). DHR has been previously used in different mi- croalgal species to measure oxidative burst (Jamers et al., 2009; Vardi et al., 1999). Samples were incubated with both 10 μM DHR and PUA in the dark at room tempera- ture and analysed using flow cytometry (see below). Flow cytometry To quantify ROS production in vivo, a flow cytometer (FACScalibur, Becton Dickinson, Palo Alto, CA, USA) equipped with a 488-nm laser as excitation source and fil- tered sea water as sheath was used. A 530/30-nm beam- pass emission filter was used for detection of DHR-derived green fluorescence. Both the sheath fluid and the sample velocity (65 μLmin–1) were kept constant during all the ex- periments. As a trigger signal, red fluorescence was used with a threshold at channel 52. Red fluorescence was col- lected through a 650-nm long-pass filter and was also used as a proxy for chlorophyll cell fluorescence. Fluorescent beads (Coulter FlowSet Fluorospheres, Beckman Coulter, Fullerton, CA, USA) were used as internal standard. Green fluorescence values were normalized to the green fluores- cence of the beads, and all data were expressed as fold changes of green fluorescence relative to the control val- ues. Data acquisition (104 cells on average for each sample) was performed using CellQuest software (Becton-Dikin- son). Data analysis was performed using FCS4 Express (De Novo Softwares, Los Angeles, CA, USA). Total RNA extraction Two-hundred mL of culture were collected through a fil- tration system using a 2-µm pore size polycarbonate filter (Nucleopore; Whatman, Maidstone, UK), placed over a No n- co mm er cia l u se on ly 16DSP expression in response to PUA GF/F glass microfiber filter (Whatman). Filtration never lasted longer than 30 min. Cells were then resuspended in 0.7 mL guanidine isothiocyanate buffer (RLT buffer, Qia- gen, Valencia, CA, USA) containing 1% β-mercaptoethanol (Sigma Aldrich, St. Louis, MO, USA), immediately frozen at -80°C, and stored for a few days. After disrupting the cells by sonication on ice (Sonicator ultrasonic processor XL, Heat System Ultrasonics, Farmingdale, NY, USA), total RNA was extracted using the silica-membrane spin column included in the RNeasy Plant Mini Kit (Qiagen). The resid- ual genomic DNA was removed by an additional treatment with the RNase-Free DNase Set (Qiagen) according to the manufacturer’s instruction. Total RNA concentration and purity were determined with a spectrophotometer (ND-100; NanoDrop Technologies, Wilmington, DE, USA) by meas- uring absorption at 260 and 280 nm. Reverse transcription and Real time quantitative polymerase chain reaction One μg of total RNA was treated with DNase I and re- verse transcribed into first-strand complementary DNA (cDNA) fragments by using random hexamer oligonu- cleotides and the High Capacity cDNA Reverse Transcrip- tion Kit (Applied Biosystems, Foster City, CA, USA). The conditions for the retro-transcription were set as in Chung et al. (2005): 25°C for 10 min, 37°C for 120 min, and 85°C for 5 min. Real time quantitative polymerase chain reac- tions (Q-PCR) were performed using 2× SYBR Green PCR master mix (Applied Biosystems) with 250 nM of each of the forward and reverse primers. The nucleotide sequences of all primer pairs used in the Q-PCRs are indi- cated in Tab. 1. The GeneAmp 7000 sequence detection system (Applied Biosystems) was used to carry out all the reactions. The threshold cycle at which the fluorescence intensity exceeded a preset threshold was used to calculate the target gene mRNA and 18S rRNA expression levels. The RNA molar ratio of ScDSP mRNA and 18S rRNA was calculated using the formula described in Chung et al., (2005), and the data normalized to the control values. The levels of ScDSP mRNA expression were expressed as fold changes of mRNA abundance relative to the control. All data reported are means of three technical replicates but the experiments were repeated at least twice. Statistical analyses Pairwise comparisons among treatments and control were assessed by Student’s t-test using the Excel spread- sheet (Microsoft, Redmond, Washington, MA, USA). RESULTS AND DISCUSSION Cultures were affected by OCTA exposure. 16 μM OCTA reduced cell concentrations to 51.17% of initial val- ues after 24 h (Fig. 1; Tab.2). H2O2 also strongly affected cell concentrations which decreased to 50% as early as 2 h after inoculation and to zero 24 h later (Fig. 1). ScDSP-1 and ScDSP-2 gene expression and ROS lev- els increased in S. tropicum cultures exposed to different concentrations of OCTA. A peak in ScDSP-1 and ScDSP- 2 expression was visible as early as 1 h after exposure and its intensity for cells in LE exposed to the concentration of 8 μM OCTA was ca. 822 and 98 folds that in the control for ScDSP-1 and ScDSP-2, respectively (Fig. 2 a,b). Gene expression decreased afterward and showed no further peak until the end of the experiment. In the same strain, an increase in ScDSP expression has been previously re- ported when cultures entered senescence (Chung et al., 2005) and also in response to a light stress involving the blockage of electron flow at the level of the photosystem II (Chung et al., 2008). The latter also involved NO in- crease in the gene regulation and involvement of a PCD machinery (Chung et al., 2008). Overall, these results pro- vide evidence that the DSP genes have a role in regulating the molecular mechanism of mortality in S. tropicum under stress. In our experiments, gene expression was highest at the concentrations where ROS production was intermediate, suggesting that there exists a threshold for the activation of the two death-related genes investigated. A peak in ROS was observed 20 min after inoculation with OCTA (Fig. 3), followed by a decrease. The intensity of the peak Tab. 1. Primer pairs used for the qRT PCR reactions, as in Chung et al. (2008). Target gene Primer name Nucleotide sequence (5’-end to 3’-end) ScDSP-1 ScDSP-SG-F GAACA AGCAA ACTGC ACTCG TC ScDSP-SG-R GTCAA GAATG TTGGT CGTCG CG ScDSP-2 ScDSP-SG-F GAACA AGCAA ACTGC ACTCG TC ScDSP2-SG-R GTAGG CATCT GCTAT TCTTT CTG 18S rRNA Ske-18S-F GAATT CCTAG ATATC GCAGT TCATC Ske-18S-R GCTAA TCCAC AATCT CGACT CCTC Tab. 2. Effect of OCTA on growth rate of S. tropicum cultures after 24 h. Initial concentration (t0, i.e., LE) was 2.53 105 cell mL–1. Growth rate (day–1) Concentration (µM) LE OCTA° 0 0.30±0.10 4 0.03±0.08 8* -0.63±0.39 16 -0.80±0.22- 16H2O2 2.70±0.35 LE, Late exponential growth phase (as in Fig. 1); °2E,4E/Z-octadienal; *concentration causing a reduction of 50% in cell concentrations at 24 h (see text). Data are means ±SD, n=3. No n- co mm er cia l u se on ly 17 A.A. Gallina et al. was PUA concentration-dependent, and was highest for 16 μM OCTA (1.9±0.2 increase; P<0.001, n=3). Both the ROS scavenger Tempol and the SOD inhibitor DETC re- duced the intensity of the peak, therefore supporting the conclusion that the peak was really due to an increase of ROS. No increase in DHR-derived green fluorescence was observed when methanol only was added to the cul- tures, excluding an artifact due to the solvent used (data not shown). Hence, it is suggested that cell death is trig- gered in S. tropicum only by specific PUA concentrations, which are those causing a reduction of 50% in cell con- centrations at 24 h (8 μM OCTA). This is supported by the observation that 8 μM OCTA induced intermediate levels of ROS production (1.5±0.1 increase in LE cul- tures), while a higher concentration (16 μM OCTA) elicited a considerably lower ScDSP expression, even though it induced a higher ROS production (1.9±0.2 in- crease). This could be related to different molecular sig- nals triggered in the two cases. In our experiments, while the addition of H2O2 to the cultures induced an increased DSP gene expression, the levels were not as high as those elicited by OCTA, even though the DHR-ROS derived green fluorescence was the highest observed (7.5±0.9 times the control at 40 min after exposure, not shown). This suggests that ROS could be important modulators in the pathway leading to cell death. Even though evidence for the occurrence of apoptotic-like features in S. tropicum in response to PUA was not pro- vided in this study, it has been previously reported that PUA are able to induce apoptotic-like cell death in the ma- rine diatom Thalassiosira weissflogii (Casotti et al., 2005). In addition, in the marine diatom Phaeodactylum tricornutum PUA exposure induced overexpression of a gene associated to NO production (designated as PtNOA), that resulted in an altered expression of superoxide dis- mutase and an increased metacaspase activity, which are known factors in the stress responses and PCD pathway (Vardi et al., 2008). Indeed, in our OCTA-exposed S. trop- icum cultures, PUA had a toxic effect drastically affecting growth and inducing cell death. ROS are generally indicated in the literature as PCD- inducing factors (Aldsworth et al., 1999; Simon et al., 2000). However, we suggest that a difference exists be- tween the intracellular effectors of a PCD-like process (e.g., NO and ROS), and other factors and stimuli that, causing an alteration of NO and ROS, are able to modulate the response leading to cell death. Based on the results re- ported here, it is proposed that PUA act as such factors, which induce modifications of the intracellular ROS lev- els, which then result in a differential activation of ScDSP genes. The successive reaction depends then on the amount of ROS produced and possibly on the cell physiological and growth status. This scenario is supported by the ob- servation that the increase in DSP gene expression in our experiments was not linearly correlated to PUA concen- trations. While it is possible that the experimental design was not appropriate to detect sudden peaks occurring at intermediate time points, it is unlikely that these lasted for shorter than 1 h. Indeed, when a further time point at 30 min after inoculation was included, the expression levels were not even comparable to the highest expression meas- ured during the experiments (data not shown). Therefore, Fig. 2. ScDSP1 (a) and ScDSP2 (b) mRNA expression in S. tropicum cultures exposed to different OCTA concentrations and H2O2 in LE cultures. Data are expressed as mean of fold changes of mRNA abundance (relative to control) from technical triplicates. The ex- periments were repeated at least twice. No n- co mm er cia l u se on ly 18DSP expression in response to PUA considering that a burst in ROS production always oc- curred 20 min after exposure to all the different PUA tested, it is possible to conclude that the concentration of 8 μM OCTA represents a threshold value below which the cell activates defences mechanisms to cope with the re- lated oxidative stress. (Chung et al., 2005, 2008). Based on the results reported here, it is suggested that in the case of PUA-induced chemical stress, DSP genes ex- pression is regulated by ROS. Although we cannot exclude a role for NO in the process, data from the congeneric S. marinoi show that ROS and not NO have a main role in the response to PUA in this diatom and probably in the con- generic S. tropicum (Gallina et al., 2014). In addition, con- sidered that DSP is reported to be a Ca2+ regulated protein (Chung et al., 2005), it is then possible that in S. tropicum ROS are acting as Ca2+ mobilizing agents, as also happens in other biological systems (Gonzalez et al., 2002; Mallilankaraman et al., 2011; Wendehenne et al., 2004). CONCLUSIONS The reported findings are only preliminary and need to be substantiated by further experiments aimed at prov- ing a clear connection between ROS and DSP expression. In addition, direct data on PCD activation have to be ac- quired, possibly using the genomic approach used by Lau- ritano et al., (2015). However, based on these preliminary results and combining them on the vast literature on PUA in diatoms, we propose that in Skeletonema spp. the ROS downstream response to PUA activates either a protection pathway (i.e., antioxidant defence) or a cell death cascade either through activation of specific genes or through non- genetically controlled cell death, depending on the PUA used, its concentration, the time of exposure, the concen- tration of ROS produced and the physiological state of the cell. This has important implications for bloom regulation and population dynamics and induction of cell lysis, as instance, immediately after the peak stage of a diatom bloom, when cell concentrations are highest but environ- mental factors become limiting and do not support further growth (Vanboekel et al., 1992). It has been proposed that the evolution of PCD in uni- cellular organisms as a byproduct of selection might be explained by the possibility that PCD could be a maladap- tive stress-induced process that can be selected under cer- tain conditions to benefit only some organisms or part of a population (Nedelcu et al., 2010). As related to this con- cept, if an intrapopulation signal (e.g., PUA) triggers the production of PCD-inducing signal (e.g., ROS), and if this intrapopulation signal-dependent death is beneficial for a part of the population (for instance, by fitting cell density to the environmental conditions), then PCD could be adaptively selected for and consequently evolve into al- truistic adaptation (Nedelcu et al., 2010). This hypothesis needs to be tested in situ, where PCD has been shown to exist, for instance in reaction to a viral attack (Bidle et al. 2007), supporting its role as population regulatory mech- anism, or, triggering the production of survival forms such as cysts or resting stages, using high ROS levels as signals (Vardi et al., 1999), so to assure perpetuation of the species later on, when the environment conditions become favourable again (McQuoid et al., 2002). ACKNOWLEDGMENTS AAG deeply acknowledges SH Hung and SF. Wang for their assistance during experiments and technical sup- port. Travel fund for AAG has been supported by the Na- tional Science Council of Taiwan (NSC) through the Summer Program for Italian Graduate Students fellow- ship. AAG has been supported by a SZN-OU PhD fellow- ship. Materials for this study were partly supported by the Taiwan National Science Council (NSC 99-2611-M-019 -014 -MY2). REFERENCES Aldsworth TG, Sharman RL, Dodd CER, 1999. Bacterial suicide through stress. Cell. Mol. Life Sci. 56: 78-383. Andrianasolo EH, Haramaty L, Degenhardt K, Mathew R, White E, Lutz R, Falkowski P, 2007. Induction of apoptosis by diterpenes from the soft coral Xenia elongata. J. Nat. Prod. 70:1551-1557. Fig. 3. ROS production in S. tropicum exposed to different con- centrations of OCTA. Data are values of DHR123-derived green fluorescence normalized to the control values. Data are means ±SD from three biological replicates. No n- co mm er cia l u se on ly 19 A.A. Gallina et al. Armbrust EV, 2009. The life of diatoms in the world’s oceans. Nature 459:185-192. Berges JA, Falkowski PG, 1998. Physiological stress and cell death in marine phytoplankton: Induction of proteases in re- sponse to nitrogen or light limitation. Limnol. Oceanogr. 43:129-135. Bidle KD, Bender SJ, 2008. Iron starvation and culture age ac- tivate metacaspases and programmed cell death in the ma- rine diatom Thalassiosira pseudonana. Eukaryot. Cell 7:223-236. Bidle KD, Falkowski PG, 2004. Cell death in planktonic, photo- synthetic microorganisms. Nat. Rev. Microbiol. 2:643-655. Bidle KD, Haramaty L, Ramos JBE, Falkowski P, 2007. Viral activation and recruitment of metacaspases in the unicellular coccolithophore, Emiliania huxleyi. P. Natl. Acad. Sci. USA 104:6049-6054. Butow BJ, Wynne D, TelOr E, 1997. Antioxidative protection of Peridinium gatunense in Lake Kinneret: seasonal and daily variation. J. Phycol. 33:780-786. Casotti R, Mazza S, Brunet C, Vantrepotte V, Ianora A, Miralto A, 2005. Growth inhibition and toxicity of the diatom alde- hyde 2-trans, 4-trans-decadienal on Thalassiosira weissflogii (Bacillariophyceae). J. Phycol. 41:7-20. Cervia D, Di Giuseppe G, Ristori C, Martini D, Gambellini G, Bagnoli P, Dini F, 2009. The secondary metabolite euplotin C Induces apoptosis-like death in the marine ciliated protist Euplotes vannus. J. Eukaryot. Microbiol. 56:263-269. Chandra J, Samali A, Orrenius S, 2000. Triggering and modula- tion of apoptosis by oxidative stress. Free Radical Biology and Medicine 29, 323-333. Chung CC, Hwang SPL, Chang J, 2005. Cooccurrence of ScDSP gene expression, cell death, and DNA fragmentation in a marine diatom, Skeletonema costatum. Appl. Environ. Microbiol. 71:8744-8751. Chung CC, Hwang SPL, Chang J, 2008. Nitric oxide as a sig- naling factor to upregulate the death-specific protein in a marine diatom, Skeletonema costatum, during blockage of electron flow in photosynthesis. Appl. Environ. Microbiol. 74:6521-6527. Costas E, Aguilera A, Gonzalezgil S, Lopezrodas V, 1993. Con- tact inhibition: also a control for cell proliferation in unicel- lular algae? Biol. Bull. 184:1-5. Deponte M, 2008. Programmed cell death in protists. BBA Mol. Cell Res. 1783:1396-1405. Ellegaard M, Ribeiro S, Lundholm N, Andersen TJ, Berge T, Ekelund F, Harnstrom K, Godhe A, 2013. Using the sedi- ment archive of living dinoflagellate cysts and other protists resting stages to study dynamics through time, pp. 149-153. In: F. Marrett, J.M. Lewis and L.R. Bradely (eds.), Biolog- ical and Geological perspectives of dinoflagellates. Geolog- ical Society, London. Evans C, Malin G, Mills GP, Wilson WH, 2006. Viral infection of Emiliania huxleyi (Prymnesiophyceae) leads to elevated pro- duction of reactive oxygen species. J. Phycol. 42:1040-1047. Gallina AA, Brunet C, Palumbo A, Casotti R, 2014. The Effect of polyunsaturated aldehydes on Skeletonema marinoi (Bacillariophyceae): the involvement of reactive oxygen species and nitric oxide. Mar. Drugs 12:4165-4187. Gonzalez A, Schmid A, Salido GM, Camello PJ, Pariente JA, 2002. XOD-catalyzed ROS generation mobilizes calcium from intracellular stores in mouse pancreatic acinar cells. Cell. Signal. 14:153-159. Guillard RRL, 1975. Culture of phytoplankton for feeding ma- rine invertebrates, pp. 29-60. In:, W. L. Smith and M.H. Chanley (eds.), Cultures of marine invertebrate animals. Plenum Publishing Co., New York. Ischiropoulos H, Gow A, Thom SR, Kooy NW, Royall JA, Crow JP, 1999. Detection of reactive nitrogen species using 2,7- dichlorodihydrofluorescein and dihydrorhodamine 123. In: J.N. Abelson, M.I. Simon and H. Sies (eds.). Methods in en- zymology - Nitric oxide Part C - Biological and antioxidant activities, vol. 301. Academic Press, San Diego. Jamers A, Lenjou M, Deraedt P, Van Bockstaele D, Blust R, de Coen W, 2009. Flow cytometric analysis of the cadmium- exposed green alga Chlamydomonas reinhardtii (Chloro- phyceae). Eur. J. Phycol. 44, 541-550. Jimenez C, Capasso JM, Edelstein CL, Rivard CJ, Lucia S, Breusegem S, Berl T, Segovia M, 2009. Different ways to die: cell death modes of the unicellular chlorophyte Dunaliella viridis exposed to various environmental stresses are mediated by the caspase-like activity DEVDase. J. Exp. Bot. 60:815-828. Jones DP, 2008. Radical-free biology of oxidative stress. Am. J. Physiol. Cell Physiol. 295:C849-C868. Kuwabara M, Asanuma T, Niwa K, Inanami O, 2008. Regulation of cell survival and death signals induced by oxidative stress. J. Clin. Biochem. Nutr. 43:51-57. Lam E, 2008. Programmed cell death in plants: orchestrating an intrinsic suicide program within walls. Crit. Rev. Plant Sci. 27:413-423. Lauritano C, Orefice I, Procaccini G, Romano G, Ianora A, 2015. Key genes as stress indicators in the ubiquitous diatom Skeletonema marinoi. BMC Genomics 16:411. Luan S, Kudla J, Rodriguez-Concepcion M, Yalovsky S, Gruis- sem W, 2002. Calmodulins and calcineurin B-like proteins: calcium sensors for specific signal response coupling in plants. Plant Cell 14:S389-S400. Mallilankaraman K, Gandhirajan RK, Hawkins BJ, Madesh M, 2011. Visualization of vascular Ca2+ signaling triggered by paracrine derived ROS. J. Vis. Exp. DOI: 10.3791/3511 Martindale JL, Holbrook NJ, 2002. Cellular response to oxida- tive stress: signaling for suicide and survival. J. Cell. Phys- iol. 192:1-15. McQuoid MR, Godhe A, Nordberg K, 2002. Viability of phyto- plankton resting stages in the sediments of a coastal Swedish fjord. Eur. J. Phycol. 37:191-201. Means AR, Dedman JR, 1980. Calmodulin - an intracellular cal- cium receptor. Nature 285:73-77. Nedelcu AM, Driscoll WW, Durand PM, Herron MD, Rashidi A, 2010. On the paradigm of altruistic suicide in the unicel- lular world. Evolution 65:3-20. Pippen EL, Nonaka M, 1958. A convenient method for synthe- sizing normal aliphatic 2,4-dienals. J. Org. Chem. 23:1580- 1582. Ribalet F, Berges JA, Ianora A, Casotti R, 2007. Growth inhibition of cultured marine. phytoplankton by toxic algal-derived polyunsaturated aldehydes. Aquat. Toxicol. 85:219-227. Rockenfeller P, Madeo F, 2008. Apoptotic death of ageing yeast. Exp. Gerontol. 43: 876-881. Scherz-Shouval R, Elazar Z, 2007. ROS, mitochondria and the No n- co mm er cia l u se on ly 20DSP expression in response to PUA regulation of autophagy. Trends Cell Biol. 17:422-427. Segovia M, 2011. Cell death or survival in the unicellular chlorophyte Dunaliella: the role of metacaspases and cas- pase-like proteases. Eur. J. Phycol. 46:32-33. Segovia M, Berges JA, 2009. Inhibition of caspase-like activities prevents the appearance of reactive oxygen species and dark-induced apoptosis in the unicellular chlorophyte Dunaliella tertiolecta. J. Phycol. 45:1116-1126. Simon HU, Haj-Yehia A, Levi-Schaffer F, 2000. Role of reactive oxygen species (ROS) in apoptosis induction. Apoptosis 5:415-418. Thamatrakoln K, Bailleul B, Brown CM, Gorbunov MY, Kustka AB, Frada M, Joliot PA, Falkowski PG, Bidle KD, 2013. Death-specific protein in a marine diatom regulates photo- synthetic responses to iron and light availability. P. Natl. Acad. Sci. USA 110:20123-20128. Thamatrakoln K, Korenovska O, Niheu AK, Bidle KD, 2012. Whole-genome expression analysis reveals a role for death- related genes in stress acclimation of the diatom Thalas- siosira pseudonana. Environ. Microbiol. 14:67-81. Thibane VS, Ells R, Hugo A, Albertyn J, van Rensburg WJJ, Van Wyk PWJ, Kock JLF, Pohl CH, 2012. Polyunsaturated fatty acids cause apoptosis in C. albicans and C. dubliniensis biofilms. Biochim. Biophys. Acta 1820.1463-1468. Vanboekel WHM, Hansen FC, Riegman R, Bak RPM, 1992. Lysis-induced decline of a Phaeocystis spring bloom and coupling with the microbial food web. Mar. Ecol.-Progr. Ser. 81:269-276. Vardi A, Berman-Frank I, Rozenberg T, Hadas O, Kaplan A, Levine A, 1999. Programmed cell death of the dinoflagellate Peridinium gatunense is mediated by CO2 limitation and ox- idative stress. Curr. Biol. 9:1061-1064. Vardi A, Bidie KD, Kwityn C, Hirsh DJ, Thompson SM, Callow JA, Falkowski P, Bowler C, 2008. A diatom gene regulating nitric-oxide signaling and susceptibility to diatom-derived aldehydes. Curr. Biol. 18:895-899. Vardi A, Eisenstadt D, Murik O, Berman-Frank I, Zohary T, Levine A, Kaplan A, 2007. Synchronization of cell death in a dinoflagellate population is mediated by an excreted thiol protease. Environ. Microbiol. 9:360-369. Vardi A, Formiggini F, Casotti R, De Martino A, Ribalet F, Mi- ralto A, Bowler C, 2006. A stress surveillance system based on calcium and nitric oxide in marine diatoms. PLoS Biol. 4:411-419. Vardi A, Schatz D, Beeri K, Motro U, Sukenik A, Levine A, Ka- plan A, 2002. Dinoflagellate-cyanobacterium communica- tion may determine the composition of phytoplankton assemblage in a mesotrophic lake. Curr. Biol. 12:1767-1772. Vardi A, Van Mooy BAS, Fredricks HF, Popendorf KJ, Ossolin- ski JE, Haramaty L, Bidle KD, 2009. Viral glycosphin- golipids induce lytic infection and cell death in marine phytoplankton. Science 326:861-865. Wendehenne D, Durner J, Klessig DF, 2004. Nitric oxide: a new player in plant signalling and defence responses. Curr. Opin. Plant Biol. 7:449-455. Wichard T, Poulet SA, Halsband-Lenk C, Albaina A, Harris R, Liu DY, Pohnert G, 2005. Survey of the chemical defence potential of diatoms: screening of fifty one species for alpha,beta,gamma,delta-unsaturated aldehydes. J. Chem. Ecol. 31:949-958. No n- co mm er cia l u se on ly