Polyunsaturated Aldehydes from Large Phytoplankton of the Atlantic Ocean Surface (42°N to 33°S)

Polyunsaturated aldehydes (PUAs) are organic compounds mainly produced by diatoms, after cell wounding. These compounds are increasingly reported as teratogenic for species of grazers and deleterious for phytoplanktonic species, but there is still scarce information regarding concentration ranges and the composition of PUAs in the open ocean. In this study, we analyzed the spatial distribution and the type of aldehydes produced by the large-sized (>10 μm) phytoplankton in the Atlantic Ocean surface. Analyses were conducted on PUAs released after mechanical disruption of the phytoplankton cells, referred to here as potential PUAs (pPUAs). Results show the ubiquitous presence of pPUA in the open ocean, including upwelling areas, as well as oligotrophic gyres. Total pPUA concentrations ranged from zero to 4.18 pmol from cells in 1 L. Identified PUAs were heptadienal, octadienal and decadienal, with heptadienal being the most common (79% of total stations). PUA amount and composition across the Atlantic Ocean was mainly related to the nitrogen:phosphorus ratio, suggesting nutrient-driven mechanisms of PUA production. Extending the range of trophic conditions considered by adding data reported for productive coastal waters, we found a pattern of PUA variation in relation to trophic status.


Introduction
Phytoplankton species produce a plethora of bioactive compounds [1]. A group of these chemicals are the polyunsaturated aldehydes (PUAs), an oxylipin-type compound commonly produced by marine diatoms [2]. PUAs from diatoms have been found to be highly reactive [3], but the nature and ecological significance of these compounds remain largely unknown.
PUAs are produced after lipoxidation of intracellular polyunsaturated fatty acids (PUFAs), mainly C16 and C20, once the cell membrane has been damaged [4][5][6]. In nature, processes producing the disruption of phytoplankton cells are viral infection, grazing or/and cell lysis during senescence [7,8]. As cells die, its content and derived compounds, such as these PUAs, are released into the surrounding medium. In the vicinity of the broken cell, the released PUAs create microzones, where they react with other dissolved chemicals and interact with the neighboring organisms.
To date, different types and grades of interactions of PUAs have been described, and different ecological roles of these compounds have been proposed. Released PUAs have been suggested as a chemical defense in diatoms against grazers [2,3]. Likewise, deleterious effects on different plankton species have also been reported [9][10][11][12]. On the other hand, several studies dispute the effects of PUAs on grazers [13,14]. Recently, experiments in mesocosm studies challenged the PUA effects on bacterial community structure or diversity [15], although a strain-dependent effect has been also reported for isolated oceanic bacteria [16]. Vardi et al. [17] provided evidence for a PUA-derived nitric oxide-based intercellular signaling system for the perception of stressed cells, potentially playing a role in the community succession of diatom blooms. These authors proposed that as conditions become increasingly detrimental during the final phase of the bloom, the PUA concentration exceeds a threshold value and could act as a bloom termination signal, triggering cell death.
Despite the evidence implicating the effects of PUAs at different ecological levels, basic questions regarding the natural distribution of PUA producers and the range of PUA concentrations in nature still remain open. Studies in nature have been restricted to coastal areas [18]. It has been reported that octadienal, heptadienal and octadienal seems to be the most common of the PUAs produced by diatoms in culture [19] and the most commonly observed in natural blooming events [18]. Dissolved PUAs quantified after a bloom event of Skeletonema marinoi in the northwestern Adriatic Sea [18] ranged from 0.05 nM to 128 nM, and potential PUA levels (pPUA), that is, PUA concentrations obtained by artificial cell disruption of phytoplankton, reached 18.1 nmol per cell (>1.2 μm) from 1 L [18]. These observed ranges of pPUA in nature are very far from the micromolar concentrations usually employed in experimental toxicity assays [10], where biomass concentrations exceed natural phytoplankton abundance by several orders of magnitude. As has been pointed out by others authors [20], further analysis of natural plankton samples under non-blooming conditions will help us to understand the fluctuations and rate of production of these chemicals in the ocean.
In this work, we analyze the pPUA derived from large phytoplankton in the nutrient-poor open ocean. The potential PUA concentrations were obtained by mechanical cell disruption of large-sized class phytoplankton (>10 μm) collected from the surface waters (0-3 m) of the Atlantic Ocean. In an ecological context, the levels of pPUA must be considered as the dissolved PUA concentration that would be released in 1 L of seawater at the sampled location if 100% of the collected phytoplanktonic cells were wounded naturally (e.g., by grazers, pathogens or lysis). A large phytoplankton fraction (>10 μm) is generally composed of eukaryotic phytoplankton, mostly diatoms, dinoflagellates and coccolithophorids. They usually represent a minority fraction of the phytoplankton biomass in the open ocean, but the large phytoplankton and, especially, diatoms appear to play a more important role as oceanic productivity increases [21,22]. Given the existing knowledge about PUA production in eutrophic coastal waters [18], here, we aim to examine the nature and amount of pPUA derived from large phytoplankton over a range of oceanic conditions.

Study Site
The cruises transected widely different areas in the Atlantic Ocean, embracing a range of environmental gradients in the open ocean ecosystem. Physicochemical characteristics of water masses (temperature, salinity, nutrients) ultimately determine the composition, size structure and cell physiological state of the phytoplankton community. Longhurst [23] proposed a subdivision of the ocean based on the above-mentioned water characteristics. This biogeographical regionalization of the Atlantic Ocean constitutes a useful tool for describing and comparing the pPUA distributions obtained in this wide sampling area. Following Longhurst's criteria, the various biogeographical provinces were identified for the four cruise tracks ( Figure 1, Table 1). Thus, latitudinal northern transects (T1 and T2 in Figure 1) visited mainly two biogeographical provinces: Caribbean Province (CARB) and North Atlantic Subtropical Gyre (NASE), passing through Sargasso Sea and the Azores region. The Sargasso Sea, with extremely low nutrient availability, denoted the boundary between the eastern and western section of this transect by defining a strong trophic gradient ( Table 1). The meridional transect (T3 in Figure 1) began in Cá diz (Spain) and went southward to Brazil, travelling through NASE, the North Atlantic and the Western Tropical provinces (NATR; WTRA) and the South Atlantic Gyral Province (SATL). This transect toured most of the ultraoligotrophic ocean and crossed the equatorial upwelling. Finally, the southern latitudinal transect (T4 in Figure 1) also visited the tropical South Atlantic oligotrophic gyre and delved into the Benguela upwelling region (BENG). In the South Atlantic Gyre, chlorophyll values are generally low throughout the year, with patches of higher phytoplankton concentrations inflowed from Benguela between October and December [23].
Our cruises sampled along a latitudinal gradient of surface temperature, with maximum temperature (>27 °C) at equatorial stations ( Figure 1). The thermal gradient was especially evident in transects T1 and T2 (Figure 1), ranging from 14 °C to 27 °C. Stations sampled in the summer (T2, T4 and the southern part of T3; Supplementary Table S1) showed relatively warm waters (Supplementary Figure S1 and Table S1).
Field measurements of total chlorophyll a (TChla) matched the pattern derived from MODIS ( Figure 1; Table 1 Table S1).

Surface pPUA Distribution and Biogeographical Patterns
The chemical analyses of the extracts obtained from the large phytoplankton fraction (>10 μm) showed the presence of three polyunsaturated aldehydes: heptadienal (HEPTA), octadienal (OCTA) and decadienal (DECA). Traces of octatrienal were also found in some samples, generally in very low concentrations, near the limit of detection, and ultimately have not been considered.
pPUA were present at 95% of the sampling sites ( Figure 2; Supplementary Table S1) denoting a widespread distribution of large-sized PUA producers in the surface waters of the open ocean during the sampling period. HEPTA had a wide distribution, being present at 79% of the sampling sites, while OCTA and/or DECA were present at 61% of the sites (Supplementary Table S1). It is known that PUA producers usually release a combination of different aldehydes [6], although less is known about the ratios of production of each species, i.e., the HEPTA:OCTA:DECA ratio, and its variation in relation to environmental conditions. On the basis of the available literature, HEPTA is the most common PUA produced by diatoms [19], in agreement with our observations (Figure 2). However, DECA also reached considerably high concentrations in the transect, T2 ( Figure 1; Table 1; Supplementary Table S1) Table 1. Biogeographical province main characteristics (following Longhurst, 1998). The transects are as in Figure 1. Data are expressed as the average ± SD. PUA, polyunsaturated aldehyde; Chla, chlorophyll a.   Supplementary Table S1). Increasing color darkness denotes increasing concentrations. See the legend below for detailed ranges. Biogeographical province acronyms are as used in Figure 1. A whisker-box plot of PUA concentrations grouped by bio-geographical ocean region is shown in Supplementary Figure S2.
Significant differences were also found when DECA was compared regionally (Kruskal-Wallis test; chi-squared = 18.48, p-value = 0.0024), being more abundant in northern tropical regions (NATR and WTRA) compared with southern ones, but also showed higher values on equatorial upwelling (included in the WTRA province) (Figure 2). OCTA and DECA followed a complementary pattern of distribution regionally (Supplementary Figure S2), with DECA being produced in the tropical Atlantic Ocean. OCTA and DECA distributions showed a significant and negative correlation (Spearman's rank correlation coefficient = −0.2808), while DECA and HEPTA were positively correlated (Spearman's rank correlation coefficient = 0.4569).
The latitudinal northern transects (T1 and T2 in Figure 1) included the provinces, CARB, NASE and NATR. Within these areas, plankton is subject to considerable seasonal variation in the environmental forcing [23,25]. Blooms of large diatoms usually develop in mid-April (>5 μm), and a transition to oligotrophic conditions occurs at the end of June. Transect T1 was sampled during March-April (boreal spring), and transect T2 was sampled during June-July (boreal summer); thereby, seasonal differences in the abundance and composition of phytoplanktonic community must be expected.
Phytoplankton biomass along T2 (averaged TChla: 0.07 ± 0.05 mg m −3 ) was much lower than that obtained for T1 (averaged TChla: 0.19 ± 0.14 mg m −3 ), and significant differences in nutrient availability were observed between western and eastern Sargassum (Figure 2, Table 1). The summer population of NASE (0.039 ± 0.021 mg Chla m −3 ) was associated with highly oligotrophic conditions (Table 1) and produced significantly higher amounts of PUA compared to spring or winter communities (Supplementary Figure S3) (Kruskal-Wallis test, chi-squared = 10.0766, p-value = 0.006485). When we consider each PUA chain length separately, DECA and HEPTA were higher in summer, compared to OCTA, which was not detected (Supplementary Table S1).
Our results show that NASE assemblages collected in summer showed higher pPUA production and a different PUA production ratio than that collected in spring or winter, with an apparent increase of DECA production. Moreover, as occurred for the comparison of different biogeographical regions (Supplementary Figure S2), DECA and OCTA production followed a complementary temporal pattern, with higher DECA in summer associated with lower OCTA (Supplementary Figure S3).
The production of PUA has been shown to be strain dependent [26]. Different strains of the same species, isolated in different seasons, produced different quantities of PUA in culture. Taylor and collaborators demonstrated that summer strains of Skeletonema marinoi had a higher production of PUA [26] and postulated that strains that were isolated during times of higher grazing pressure and less beneficial nutrient conditions had higher pPUA production. This is consistent with our results, thus extending to natural assemblages under physiological (nutrient) stress. This fact is especially interesting, because the observed patterns of variability of pPUA in nature shows the net effect of taxonomical, physiological and environmental effects, which are hard to reproduce in laboratory conditions. Seasonal environmental changes, such as day length or temperature, could also impact PUA persistence in water, toxicity or the signaling role. It has been experimentally demonstrated that temperature significantly affects the persistence of HEPTA, OCTA and DECA in natural Atlantic seawater, being more persistent at 10 °C than at either 15 °C or 20 °C. Moreover, DECA was consistently more persistent than OCTA and HEPTA at the higher temperature assayed [27]. The higher persistence of DECA at warmer waters compared with HEPTA or OCTA could imply less sensitivity to DECA of the low latitude strains (and species) compared with the native species of higher latitudes. Depending on the average levels of exposure of the cells to released PUAs, the type and quantity of these compounds could have a different effect on plankton community structure in temperate waters compared with, for instance, polar and sub-polar environments.
Our main interest in this work was not to identify PUA producers, but to describe the range and composition of PUAs potentially present in open oceanic areas. In fact, we should not rule out the possibility that other non-diatom large-sized class producers were responsible for the pPUA concentration observed. Most oligotrophic oceanic areas are dominated by picoplanktonic species [28], which would not be collected by a big-sized mesh, such as the one used in our sampling. However, large-sized colonies (>1 mm) of Trichodesmium (Cyanophyta) are especially abundant in oligotrophic tropical regions characterized by warm (>22° C) surface waters [25], the equatorial Atlantic region from 5°S to 15°N and also are present in the North Atlantic Subtropical Gyre and the western Sargasso Sea [28][29][30]. There are no reports regarding PUA production from Cyanophyta; however, Trichodesmium is able to synthesize polyunsaturated fatty acids (PUFA), such as eicosapentaenoic acid (EPA) and arachidonic acid (AA) [31]. These compounds may act as HEPTA and DECA precursors, respectively [5,6]. Other non-diatom DECA producers, such as the colony-forming Phaeocystis pouchetii (Prymnesiophyte), may also produce DECA [32]. Therefore, pPUA should not be strictly associated to diatoms, as other producers as Trichodesmium trichomes or Phaeocystis colonies are plausible.

pPUA Distribution and Trophic Gradients
The biogeographical regions described above embrace a wide range of ecosystems with diverse productivity, such as oligotrophic subtropical gyres, upwellings and THE mesotrophic shelf system. Such trophic gradation is evident in total chlorophyll and nutrient distributions (Figure 1, Supplementary Table S1). From ultraoligotrophic (TChla < 0.06 mg m −3 , sensu [33]) to mesotrophic conditions (TChla 0.1-0.3 mg m −3 , sensu [33]), large phytoplankton showed decreasing pPUA levels (Figure 3a), suggesting trophic control of PUA production. Examination of the relationship between pPUA and the availability of dissolved inorganic nutrients (nitrogen and phosphorus ratio, N:P) showed that the lowest (N:P < 1) and the highest (N:P > 20) available N:P were coincident with the lowest pPUA concentrations (Figure 3b). N:P ratios from two to 15 were related to the highest pPUA concentrations. . Nine data running averages of (a) pPUA over a gradient of total chlorophyll a (TChla), (b) pPUA and nutrients (nitrate and phosphate) over the gradient of the nitrate to phosphate ratio and (c) different chain lengths of aldehydes, heptadienal, octadienal and decadienal, over the gradient of the nitrate to phosphate ratio.
PUA production is dependent on taxonomic composition, cellular PUFA contents (PUA precursors) and enzymatic lipoxygenase activity. These variables are directly related to the cell physiological state and the trophic state of the surrounding environment and could explain the observed trend. Our results strongly suggest the control of nutrient availability on the regulation of pPUA production in large-sized phytoplankton assemblages of the Atlantic Ocean. This nutrient control has been previously observed by other authors in experimental conditions in mesocosms [34]. Several researchers showed that a single species can significantly change the type and the quantity of PUAs, depending on nutrient availability [35,36]. Ribalet and coworkers [34] proposed that the pPUA production of the diatom, Skeletonema marinoi (mainly heptadienal and octadienal), was limited by the amount of substrates (PUFA, polyunsaturated fatty acids) under P-limitation and by enzyme activity under N-limitation. This hypothesis could explain the trends represented in Figure 3, but the goal here is that the present results have been obtained for natural assemblages of mixed species of cells probably subjected to a chronic nutrient limitation.
Octadienal, but especially heptadienal and decadienal, showed an alternance in dominance along the N:P gradient (Figure 3c), and their distributions were negatively correlated (Spearman correlation coefficient: −0.280869), suggesting qualitative changes in the biochemical routes of PUA production in relation to nutrients variability. Recently, it has been described that in low-phosphorus marine systems, phytoplankton cells replace membrane phospholipids with non-phosphorus lipids [37]. Likewise, the ratios of PUAs that are produced could be affected by the differences in the biosynthetic pathways of the various PUA precursors. OCTA is mainly biosynthesized by glycolipids routes; DECA is mainly biosynthesized by phospholipids routes, and HEPTA can be synthesized by both glycolipids and phospholipids routes [5,38]. Hence, membrane remodeling under nutrient stress would affect the pool of PUA precursors, varying the pPUA ratios of the phytoplankton population.
The pPUA concentrations obtained from a phytoplankton fraction >1.2 μm in a bloom of Skeletonema marinoi in Adriatic Sea ranged from 18 to 25 nmol from cells in 1 L ( [18]; pPUA was referred to in this publication as particulate PUA). However, pPUA measured in the Atlantic Ocean did not reach concentrations of the nanomolar scale, differing three orders of magnitude from that found in the diatom bloom in the Adriatic Sea. Some possible reasons may be proposed for explaining such differences. First, the phytoplankton size fractions considered are significantly different in the two studies. pPUA in the coastal bloom of S. marinoi was related to the phytoplankton fraction larger than 1.2 μm, while our study was addressed on the size fraction larger than 10 μm. Second, cell densities during bloom events could reach 10 6 cell L −1 , which represents three to four orders of magnitude higher cell densities than that usually measured in open ocean areas (10 2 -10 3 cell L −1 ). Due to the low cell densities typically found in the open ocean, the pPUA concentrations measured for the large-sized phytoplankton fraction in open oceanic areas seem reasonable. Moreover, marine diatoms can span from a few micrometers to millimeters in length [39], and probably, lower-sized diatoms would also contribute significantly to the released PUA pool in the ocean.
We extended the trophic gradient of this study by including the published and unpublished data of phytoplankton samples collected in marine coastal areas (Figure 4). pPUA concentrations ranged from pM (equatorial upwelling) to micromole (Adriatic Sea and Shallow Bay) concentrations. The highest pPUA concentrations were measured in coastal areas with TChla over 1 mg m −3 and the lowest concentrations were found in subtropical and tropical Atlantic waters with TChla around 0.2 mg m −3 . This result can be explained from the increasing abundance of PUA producers towards eutrophic conditions. In contrast, pPUA was higher under ultra-oligotrophic conditions (TChla < 0.1 mg m −3 , sensu [33]). Since the fraction of large phytoplankton expressed as TChla is especially low at ultraoligotrophic conditions [40], the high pPUA values under oligotrophic conditions are even more evident when pPUAs are normalized by large-sized phytoplankton biomass when available (expressed as pmol PUA μg Chla −1 ) (Figure 4). We relate this result to the nutrient stress of pPUA producers, in agreement with Graneli et al. [41]. They hypothesized an increase in the release rate of chemical compounds associated with N:P unbalanced phytoplankton communities. pPUA normalized by FChla > 10 μm (pmol μg FChla −1 ) (red spots) concentrations measured at various oceanic areas (see the left axis). Data were grouped following the trophic gradient based on total chlorophyll concentration (see the sources below). (b) Box-plot for the ranges of total chlorophyll (TChla) in ecological areas. Ocean areas were arranged on the vertical axes in relation to mean Chla. Sources: 1 pPUA for phytoplankton > 1.2 μm (n = 8, [18]); total chlorophyll was derived from remotely sensedmeasurements ( [42]). 2 pPUA for phytoplankton > 10 μm; total Chla was measured in the field through standard fluorometric methods (n = 12, [43]). 3 pPUA for phytoplankton > 10 μm; total Chla was measured in the field through standard fluorometric methods (n = 68; [44]). 4 This work.
In this work, we report, for the first time, data on the pPUA distribution for open oceanic areas. The concentration ranges derived from large phytoplankton are related to trophic conditions. However, the sources of variability of pPUA are diverse, and data on natural ambient ones are still scarce. Future studies centered on the seasonal variation of pPUA and the depth profiles of dissolved PUA would be a useful next step toward understanding the ecological relevance of these compounds in nature. Physiological variables, such as cell lysis rates, should be also related to natural dissolved concentrations of PUA. Since PUA production prevails under certain trophic conditions and among certain taxa (e.g., large eukaryotes, such as diatoms), PUAs could also take part in determining the composition and size structure of the phytoplankton community.

Study Area
Samples for polyunsaturated aldehyde extractions were collected in several cruises carried out on board the B/O Hespé rides and B/O Sarmiento de Gamboa from December 2010 to July 2011 (Supplementary Table S1). A total of 71 sites were visited, and 137 viable surface samples of the phytoplankton fraction > 10 μm were collected at four different transects (T1 to T4 in Figure 1) Temperature, salinity and chlorophyll fluorescence were monitored using a CTD (SeaBird 9 plus, Bellevue, WA, USA).

Sampling
Water samples were collected by means of 12-L Niskin bottles mounted on a rosette frame. A volume of 250-500 mL was collected for the total and fractionated (>10 μm) Chla determination, filtered through Whatman ® glass microfiber filters (25 mm diameter, grade GF/F, Sigma-Aldrich, Buchs, Switzerland) and kept frozen for at least 6 hours until extraction. Pigments were extracted by placing them in 5-7 mL of 90% acetone at 4°C for 24 h, and the fluorescence of the extract was determined by means of a Turner Designs fluorometer [45] calibrated with a Chla standard (Sigma-Aldrich, Buchs, Switzerland). No phaeophytin correction was applied. Total chlorophyll (TChla) was measured in all sampling sites, but fractionated chlorophyll was only measured at T1 and T2.
Samples used for dissolved nutrients analyses were drawn from the Niskin bottles into polyethylene vials and immediately frozen (−20 °C) until analysis. Concentrations were measured spectrophotometrically with a Skalar autoanalyzer following standard procedures [46,47]. The detection limits of concentrations were 0.01 and 0.02 μM for NO 3 + NO 2 and PO 4 , respectively.
At each sampling site, two samples of a variable, but controlled, volume (from 50 to 450 L) of surface seawater (3 m) were collected for pPUA quantification. Natural seawater was pumped at 11-12 L min −1 . The sample volume was selected in relation to the CTD fluorescence signal. We generally collected a high water volume, considering that the open ocean is oligotrophic (TChla < 0.1 mg m −3 ). To quantify pPUA in large phytoplankton, water samples were firstly passed through a 200-μm mesh to eliminate the zooplankton fraction and then through a 10-μm mesh. Phytoplankton cells retained on the 10 μm mesh were concentrated into 100 mL vials by using 0.2 μm filtered seawater. Concentrated samples were filtered at low vacuum pressure onto polycarbonate filters (Nucleopore, 0.4 μm). The phytoplankton cells retained were rinsed from the filter with 1 mL of 25 mM of O- (2,3,4,5,6-pentafluorobenzyl) hydroxylamine hydrochloride (PFBHA, Fluka, Basel, Switzerland) solution in 100 mM Tris HCl, pH 7.2 [48], kept refrigerated or immersed in liquid nitrogen, depending on the cruise length, and transported to the lab. Once at the lab, samples were conserved at −80 °C until extraction.

Polyunsaturated Aldehydes Extraction and Quantification
Each sample was unthawed and immediately transferred to a conical glass vial, and 5 μL of internal standard were added (benzaldehyde, 1 mM in methanol, Sigma-Aldrich, Buchs, Switzerland). Samples were artificially disrupted by sonication on ice following the recommendations of Wichard et al. [48], and then, the samples were kept 1 h at ambient temperature to ensure a complete derivatization of aldehydes. For extraction, 0.5 mL of methanol and 1 mL of hexane were added, and the sample was vortexed for 1 min. The mixture was acidified by the addition of sulfuric acid and vortexed again [18]. The hexane phase was taken to dryness under vacuum and kept at −80 °C until analysis by gas chromatography-mass spectrometry (GC-MS). Before GC-MS analysis, the samples were dissolved in 50 μL of hexane and transferred into a 100-μL glass micro insert (Waters ® VWR International™, Radnor, PA, USA) closed with caps fitted with a PTFE septum.
GC-MS analysis was performed using an Agilent 7890A GC (Agilent Technologies Inc., Santa Clara, CA, USA) coupled to a Waters Micromass Quattro Micro (Milford, MA, USA) and equipped with a 30-m DB-5MS column (0.25 mm internal diameter, 0.25 μm stationary phase thickness, 95% dimethylpolysiloxane, 5% diphenylpolysiloxane). The GC oven temperature was programmed as follows: held at 60 °C for 1 min, then increased at 5 °C per minute to 300 °C and finally held at 300 °C for 1 min. The inlet temperature was 250 °C in splitless mode (1 μL). Helium was the carrier gas at a constant flow of 1 mL/min. The electron energy was 70 eV. The PUA-PFBHA derivatives were identified by comparison to external standards and the presence of their molecular ions (m/z 305 for heptadienal, 319 for octadienal and 347 for decadienal) and also the presence of major fragments ions, i.e., m/z 276 for the PUAs, 271 for the internal standard and 181 for all. Quantification was based on the ratio between the fragment at m/z 276 of the derivatized PUA and m/z 271 of the derivatized internal standard (that is, benzaldehyde). The calibration was performed using a synthetic standard of derivatized (25 mM PFBHA-solution in Tris-HCl at pH 7.2) 2E,4E-heptadienal (90%, Sigma-Aldrich Chemie GmbH, Steinheim, Germany), 2E,4E-octadienal (≥96%, Sigma-Aldrich Chemie GmbH) and 2E,4E-decadienal (85%, Sigma-Aldrich Chemie GmbH) in ranges from 0.5-15 nM, 15-200 nM and 200-8000 nM, to cover the wide range of molarities found in the concentrated extracts of the samples (a maximum value of 4200 nM). Calibration curves were repeated with every group of analyzed samples for three different ranges: 0.5-15 nM, 15-200 nM and 200-8000 nM of each aldehyde. The correlation coefficients obtained for each aldehyde at every range were: R 2 = 0.9884 ± 0.009, 0.9919 ± 0.011 and 0.9963 ± 0.004 for heptadienal; R 2 = 0.9875 ± 0.013, 0.9716 ± 0.037 and 0.9896 ± 0.011 for octadienal; R 2 = 0.9509 ± 0.046, 0.9823 ± 0.005 and 0.9915 ± 0.006 for decadienal, respectively.
The peak area ratios of each aldehyde (m/z 276), relative to the internal standard (m/z 271), was plotted against the concentration of the aldehyde. The obtained concentrations were finally normalized to the filtered volume of seawater, and the results were expressed as the molarity of PUA obtained from cells contained in 1 L of seawater.

Statistics
Statistical analysis of the data was performed by using statistical package R. All data were checked for normality and homoscedasticity before analysis by using the Shapiro-Wilk test (p-value > 0.05) and Levene's test (p-value > 0.05), respectively. The Kruskal-Wallis rank sum test and post hoc pairwise comparison by means of the Mann-Whitney test were applied for comparing pPUA levels.

Conclusions
In the present study, we provided the first assessment of pPUA concentrations in the open ocean and demonstrate that large-sized PUA producers are widely spread in the Atlantic Ocean. Heptadienal was the most widely distributed one. Other pPUAs, such as octadienal or decadienal, also dominated, depending on the site or season. Our results showed a relationship between pPUA in large-sized phytoplankton and nutrient availability, with higher pPUA concentrations under ultraoligotrophic conditions and at intermediate nitrate:phosphate ratios (from two to 10). Given that Chl:biomass ratio also varies in relation to the physiological status of phytoplankton, the PUA concentration in future studies should be related to conservative biomass variables, such as carbon instead of chlorophyll. Finally, we have found that the quantity of these compounds can range from the pM to μM level in the ocean, from ultraoligotrophic to eutrophic conditions.