Phytoplankton Cell Lysis Associated with Polyunsaturated Aldehyde Release in the Northern Adriatic Sea

Diatoms are able to react to biotic and abiotic stress, such as competition, predation and unfavorable growth conditions, by producing bioactive compounds including polyunsaturated aldehydes (PUAs). PUAs have been shown to act against grazers and either enhance or inhibit the growth of different phytoplankton and bacteria both in culture and in the field. Presence of nanomolar concentrations of dissolved PUAs in seawater has been reported in the North Adriatic Sea (Mediterranean), suggesting that these compounds are released in seawater following diatom cell lysis. However, the origin of the PUAs and their effects on natural phytoplankton assemblages remain unclear. Here we present data from four oceanographic cruises that took place during diatom blooms in the northern Adriatic Sea where concentrations of particulate and dissolved PUAs were monitored along with phytoplankton cell lysis. Cell lysis was positively correlated with both concentrations of particulate and dissolved PUAs (R = 0.69 and R = 0.77, respectively), supporting the hypothesis that these compounds are released by cell lysis. However, the highest concentration of dissolved PUAs (2.53 nM) was measured when cell lysis was high (0.24 d−1) but no known PUA-producing diatoms were detected, suggesting either that other organisms can produce PUAs or that PUA-producing enzymes retain activity extracellularly after diatom cells have lysed. Although in situ concentrations of dissolved PUAs were one to three orders of magnitude lower than those typically used in laboratory culture experiments, we argue that concentrations produced in the field could induce similar effects to those observed in culture and therefore may help shape plankton community composition and function in the oceans.


Introduction
Diatoms are ubiquitous organisms, responsible for about one fifth of the photosynthesis on Earth [1]. They have developed sophisticated regulatory mechanisms to perceive changes in environmental conditions and respond accordingly [2]. Some of these mechanisms rely upon the production of bioactive molecules, the most well-known being the polyunsaturated fatty acid derivatives, called oxylipins [3][4][5].
In higher plants, oxylipins play an important role in acclimation to environmental stresses due to wounding, grazing by predators or competition for limiting nutrients [6] and can trigger several physiological responses, ranging from inducing cell division to cell lysis depending on the levels of oxylipins produced [7]. Similar to plants [8], diatom oxylipins such as polyunsaturated aldehydes (PUAs) are associated with a chemical defense system that impairs the reproductive success of copepods and other invertebrates [3,5], inhibits the growth of cultured phytoplankton [9,10] and affects growth in bacteria communities both in culture [11,12] and in the field [13].
Diatoms produce PUAs in at least three steps. First, chloroplastmembrane-localized glycolipids and plasma-membrane-localized phospholipids are hydrolyzed to generate free polyunsaturated fatty acids (PUFAs) [14,15]. Then, the free PUFAs are acted on by lipoxygenases that generate hydroxyperoxy fatty acids. Finally, hydropeoxy fatty acids are transformed into PUAs [16,17]. In contrast to plants, diatom-derived PUAs are almost never found in intact cells but are synthesized mainly after membrane disruption [18]. PUA production takes place as long as enzymes are in contact with precursor free fatty acids and is not inhibited by the amount of PUAs produced [19].
Diatom PUAs appear to be part of a nitric-oxide-based stress surveillance system that mediates intracellular communication involved in regulating stress response to unfavorable environmental conditions [4]. PUA production and release in the water remains low under optimal growth conditions, but increases as conditions become poorer, e.g., nutrients become limiting [20,21]. Such a surveillance mechanism depends upon the amount of PUA produced. High PUA exposure triggers cell lysis in exponentiallygrowing diatoms while lower exposures induce resistance, leading to decreasing effects in successive exposures [4]. In contrast, PUAs can accelerate cell lysis in diatoms in the late stationary phase of growth, even if they have been exposed to PUAs previously [22]. In nature, this surveillance mechanism would confer an advantage to PUA-producing diatoms at the peak of the bloom, when nutrient conditions progressively become limiting, by inhibiting the growth of competitors. Later on, when resources become depleted and cells begin to lyse, large PUAs releases could induce a synchronized cell lysis, leading to the sharp declines of diatoms commonly observed in nature [23].
The relevance of culture experiments has been questioned because PUA concentrations are typically orders of magnitude higher than those predicted in the field [24]. However, Dittami et al. [25] have shown that repeated application of low doses (20 nM) of octadienal caused cell lysis of a non PUA-producing Thalassiosira rotula strain (CCMP 1018), a concentration that is in the same range as that reported by Vidoudez et al. [26] in surface water during blooms of the diatom Skeletonema marinoi in the northern Adriatic Sea (Mediterranean). At present our understanding of the dynamics of PUA releases and effects on natural phytoplankton assemblages is very poor.
Here we present data from an extended field survey of PUAs in the northern Adriatic Sea (Mediterranean) that was carried out during March 2002March , 2004March , 2005March and 2006. We measured phytoplankton cell lysis along with the concentrations of PUAs produced after chemical or mechanical disruption of the cells (hereafter termed ''particulate PUAs'') as well as dissolved PUAs present in seawater. Our results show that cell lysis is positively correlated with particulate and dissolved concentrations of PUAs, supporting the hypothesis that these compounds are released in seawater following cell lysis.

Methods
Four oceanographic cruises were conducted in the North Adriatic Sea in March 2002March , 2004March , 2005 and 2006 aboard the RV Urania and RV Dallaporta. Samples were taken at the surface using 5 or 10-L Niskin bottles from stations located along three transects (Fig. 1

Oceanographic data
Dissolved inorganic nutrients (phosphate, silicic acid, nitrite and nitrate) were analyzed by colorimetry with an autoanalyzer (Systea-Alliance) according to Grasshoff et al. (1999) [27]. Because nitrite concentrations were always very low (,0.3 mM), dissolved inorganic nitrogen (DIN) concentrations were expressed as the sum of the nitrite and nitrate. Samples for particulate C and N (POC and PON, respectively) were filtered on acidified precombusted (450uC, 24 h) GF/F filters and stored at 280uC until analysis with a CHN elemental analyzer (Perkin Elmer 2400). Phytoplankton cell numbers and species composition were determined on fixed samples (buffered formaldehyde 2%) after concentration by sedimentation [28], using an inverted microscope (Zeiss Axiovert 35). Samples for chlorophyll a analysis were obtained by gentle filtration (,100 mm Hg) onto GF/F filters and storage at 280uC until further analysis. The filters were ground in 100% methanol and the extract injected into a Beckman System Gold HPLC. Monthly averages of satellite-determined chlorophyll a distribution were obtained from SeaWIFS data (http://gdata1. sci.gsfc.nasa.gov/).

Measurements of cell lysis
Phytoplankton lysis rates were estimated following an improved version [29] of the esterase method [30]. The method is based on the assumption that esterases are strictly intracellular enzymes that are released only by cell lysis; esterase activity is measured as increased fluorescence due to cleavage of fluorescein diacetate (FDA). Dissolved esterase activity (DEA) was measured in 0.22 mm  immediately to the sample after the measurement at time t = 60 min. Esterase activities were measured in three replicates (coefficient of variation averaged 5%). The cell lysis rate was calculated as the decrease in PEA with time due to the production of DEA during cell lysis (produced on an hourly basis) using the following equations [29]: where DEAprod~D EA T 1=2 |0:5 and T 1/2 represents the DEA turnover time, calculated from the esterase decay rate measured in each sample.

In situ quantification of PUAs
Since PUAs are particularly reactive and volatile, they are extracted by procedures based on the transformation of PUAs into more stable products, such as ethyl ester by the Wittig reaction with carbetoxyethylidene-triphenylphosphorane (CET-TPP) [31,32], or pentafluorobenzyl-oxime derivative by O-(2,3,4,5,6pentafluorobenzyl) hydroxylamine hydrochloride (PFBHA?HCl) [24]. The quantification of particulate PUA concentration was performed following the method of d'Ippolito et al. [31] in   1  2  3  Total  1  2  3  Total  1  2  3  Total  1  2  3  Total   Bacillariophyceae  -6118  -6118  27  139  498  115  305  673  4377  1221  208  514  13194 1737 Chaetoceros affinis Alexandrium minutum --  [24]. Particulate concentration of PUAs was calculated by dividing the amount of PUAs produced after mechanical cell disruption by the number of known-PUA producing diatom cells (mainly Skeletonema marinoi) [33]. The quantification of dissolved PUAs was performed in March 2006 following the method of Vidoudez and Pohnert [22]. One liter of surface seawater was gently filtered through 1-2 cm layer of dry commercial washed sea sand using a Büchner funnel with two layers of 25 mm Whatman filter paper. The filtered seawater, together with 5 mL of internal standard (1 mM Benzaldehyde in methanol), was run through a C18 cartridge column (Chromabond C18 ec, Machery-Nagel, Düren, DE) and the PUAs eluted by adding 6 ml of the derivatizing reagent (25 mM PFBHA in methanol). The samples were incubated 1 h at room temperature and then frozen in liquid nitrogen and stored at 280uC until further analysis. Concentrations of PUAs were then analyzed by gas chromatography-mass spectrometry following Wichard et al. method [24].

Statistical Analyses
Lysis rates were compared with PUA concentrations using reduced major axis model II regression [34]. A permutation test was performed to determine the significance of the slopes and to calculate the Pearson coefficients of correlation. Statistical analyses were performed using the R computing environment (R Core Team, 2013), and judging significance at the 95% confidence level.

Results
In all cruises, satellite-based chlorophyll a concentrations indicate that the highest phytoplankton biomass in the North Adriatic Sea occurred in the Po river plume along the western coastline of transects 2 and 3 (Fig. 1). The same was true for inorganic nutrients, with the highest concentrations observed in the Po river plume (  of inorganic nitrogen and phosphate were recorded (,1.6 mM and ,0.01 mM, respectively). The concentrations of total particulate PUAs were positively correlated with lysis rates (R = 0.69, P, 0.001) (Fig. 3A), and not with chlorophyll a concentrations, nutrient concentrations or cell abundances of S. marinoi (P.0.1 in all cases) (data not shown).
Dissolved PUAs were detected at each station in March 2006, with concentrations ranging from 0.06 to 2.53 nM (Fig. 4d and Table 3). Consistent with particulate PUAs, dissolved PUAs were mostly composed of octadienal ( Table 3). The highest concentrations of dissolved PUAs were observed along transect 2 (1.09 nM), located at the Po river mouth, and coincided with high lysis rates (0.21 d 21 ) and high concentration of particulate PUAs (1.43 fmol cell 21 ). A good correlation was observed between the concentrations of total dissolved PUAs and lysis rates (R = 0.77, P = 0.005) (Fig. 3B). In general, higher concentrations of dissolved PUAs were associated with periods when S. marinoi dominated the phytoplankton community, though the highest concentration of dissolved octadienal (2.40 nM) was measured when cell lysis was high (0.24 d 21 ) but neither S. marinoi cells nor particulate PUAs were detected. No correlation was observed between concentrations of total dissolved and particulate PUAs (P.0.1) (data not shown).  Fig. 5A and 5B) over the diel cycle, driven by changes in dissolved esterase activity ( Fig. S1A and S1B). No significant differences in lysis rates were observed between dark and light periods in either experiment (Student's t-test, P.0.1 and P.0.9 for 2005 and 2006, respectively), though the statistical tests were of low power because there we relatively few data points. Particulate and dissolved PUAs were composed of 70% and 66% octadienal, respectively, with concentrations ranging over almost an order of magnitude (from 0.54 to 4.56 fmol cell 21 and 0.04 to 0.44 nM for concentrations of particulate and dissolved PUAs, respectively) ( Fig. 5A and 5B). Concentrations of total particulate and dissolved PUAs were positively correlated with lysis rates (R = 0.69, P = 0.02 and R = 0.60, P = 0.04, respectively) (Fig 4C  and 4D). No significant differences related to the time of day were found for particulate and dissolved concentrations of PUAs (Student's t-test, P.0.05 and P.0.5, respectively).

Discussion
Our results indicate that high phytoplankton lysis rates (up to 0.38 d 21 ) can occur during diatom blooms in the North Adriatic, which is in good agreement with previous phytoplankton loss rate estimates based on 14 C uptake [35]. These rates are very high, since phytoplankton net growth rates in the region range from 0.1 to 0.5 d 21 in spring [35,36]. Although esterase assays are prone to overestimation of phytoplankton lysis rates (because a small but significant fraction of both DEA and PEA can originate from heterotrophs, rather than phytoplankton [37]), our results support the idea that phytoplankton cell lysis is an important loss factor driving the dynamics of phytoplankton dynamics. This has also been noted for individual species, for example, winter blooms of the prymnesiophyte Phaeocystis spp. in the North Sea [29,30,38].
Average cell lysis rates were three times higher during the diatom blooms in the Adriatic Sea (0.16 d 21 ) than those measured during non-diatom blooms in the North Sea (0.05 d 21 ) [30,38]. Numerous factors could be responsible for this difference. Several studies have shown that diatom cell lysis can be triggered by nutrient depletion [39][40][41] or exposure to high light [42], conditions which existed during our studies. However, no significant correlations were found between lysis rates and the levels of inorganic nutrients, POC, molar particulate C:N ratio or light exposure during the survey. Field studies have shown that diatom cell lysis can also be exogenously triggered by biotic factors, such as viruses [43,44] and algicidal bacteria [45,46]. In addition, there is growing evidence that the production of algaederived toxic compounds, such as PUAs, could play a significant role in determining algal bloom dynamics and fate [47,48]. However, direct links between environmental factors and phytoplankton cell lysis rates are very difficult to demonstrate in nature; these factors are dynamic and they can interact with each other [49]. Moreover, Brussaard & Riegman [49] showed that the magnitude of cell lysis induced by nutrient limitation can be modified by the presence of bacteria.
Concentrations of total particulate PUAs varied widely across locations (from 0 to 5.37 fmol cell 21 ) but were remarkably similar among the four cruises despite significant changes in diatom compositions and environmental conditions (Table 1 and 2). At  stations where S. marinoi dominated the diatom community, particulate PUAs were mainly composed of heptadienal, and, to a minor extent, octadienal and octatrienal (Table 3), in good agreement with field observations in the region [50]. In March 2006, particulate PUAs along transect 2 were mainly composed of octadienal, with trace of decatrienal, which may be due to the high numbers of Chaetoceros compressus (Table 2), a known octadienal and decatrienal-producing species [33].
Total dissolved PUAs ranged from 0 to 2.53 nM (Table 3), which is about 2 orders of magnitude lower than concentrations measured in S. marinoi cultures (.200 nM ) [22] but in good agreement with previous in situ measurements in the region [50]. The highest concentrations of dissolved PUAs were measured when cell lysis was high (0.24 d 21 ) but at these times, no particulate PUAs were detected. This result is consistent with previous work during S. marinoi blooms showing that traces of dissolved heptadienal and octadienal were found at several stations where no particulate PUAs were detected [26]. In the latter study, one possible explanation is that not all PUA-producing organisms were collected; samples for particulate PUAs were collected using GF/C filters (nominal pore size 1.2 mm). Thus smaller cells such as picoplankton and bacteria were under-sampled, and both groups have been reported to contain small amounts of particulate heptadienal [26]. An alternative explanation is that, at the point of sampling, all PUA producing species had already lysed, leaving dissolved PUAs as the only pool. Finally, preliminary results show that lipoxygenases, key enzymes for PUA production, are still active 48 h after being released in seawater [51], suggesting that free enzymes could still produce PUAs after cells have disappeared.
Positive correlations between both concentrations of particulate and dissolved PUAs and lysis rates along transects and over the 27h time course experiments (Fig. 3) suggest a tight coupling between the three variables. We hypothesize that the positive correlation between concentrations of particulate PUAs and lysis rates results because particulate PUAs in diatoms increase under physiological stress, as observed in cultures [20,25], and diatom cell lysis similarly occurs in response to stress (e.g., nutrient limitation). Together with the fact that the composition of particulate and dissolved PUAs is very similar, the correlation between dissolved PUAs and lysis rates support the hypothesis that PUAs are released in seawater following cell lysis. An alternative hypothesis is that grazing results in dissolved PUA in the seawater. However, micro-scale examination has failed to detect production of dissolved PUAs around the mouths of grazing copepods [52].
Dissolved PUA concentrations and cell lysis are also correlated, raising the question of whether concentrations of dissolved PUAs are high enough to induce phytoplankton cell lysis. Although diatom cell growth can be inhibited by nanomolar PUA concentrations in culture [25], cell lysis has been shown to occur only when cells are exposed to micromolar PUA concentrations [4,9,10]. The concentrations of dissolved PUAs during the survey would therefore be at least one order of magnitude lower than the toxic levels established in cultures, suggesting that PUA concentrations were probably not high enough to induce cell lysis. However, bulk measurements of PUAs in seawater may not reflect the local concentrations in the vicinity of the PUA-releasing cells. Based on a simple diffusion model, micromolar concentrations of PUAs would be found at 1 mm-distance from a single S. marinoi cells (see [10] for more details). Methods for the quantification of microscale chemical gradients, such as the surface plasmon resonance detection method coupled with antibodies [53], will be needed to determine the ecological consequences of PUAs on the structure and dynamics of diatom populations.