# Isoprostanoid Profiling of Marine Microalgae

Claire Vigor, Camille Oger, Guillaume Reversat, Amandine Rocher, Bingqing Zhou, Amandyne Linares-Maurizi, Alexandre Guy, Valérie Bultel-Ponceé, Jean-Marie Galano, Joseph Vercauteren, et al.

► To cite this version:

Claire Vigor, Camille Oger, Guillaume Reversat, Amandine Rocher, Bingqing Zhou, et al. Isoprostanoid Profiling of Marine Microalgae. Biomolecules, 2020, 10 (7), pp.1073. 10.3390/biom10071073. hal-03085457

HAL Id: hal-03085457

https://hal.sorbonne-universite.fr/hal-03085457

Submitted on 2 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Isoprostanoid Profiling of Marine Microalgae

Claire Vigor<fnref n="1" />,*, Camille Oger<fnref n="1" />, Guillaume Reversat<fnref n="1" />, Amandine Rocher<fnref n="1" />, Bingqing Zhou<fnref n="1" />, Amandyne Linares-Maurizi<fnref n="1" />, Alexandre Guy<fnref n="1" />, Valérie Bultel-Poncé<fnref n="1" />, Jean-Marie Galano<fnref n="1" />, Joseph Vercauteren<fnref n="1" />, Thierry Durand<fnref n="1" />, Philippe Potin<fnref n="2" />, Thierry Tonon<fnref n="2" />,<fnref n="3" /> and Catherine Leblanc<fnref n="2" />

Received: 5 June 2020; Accepted: 14 July 2020; Published: 18 July 2020

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Algae result from a complex evolutionary history that shapes their metabolic network. For example, these organisms can synthesize different polyunsaturated fatty acids, such as those found in land plants and oily fish. Due to the presence of numerous double-bonds, such molecules can be oxidized nonenzymatically, and this results in the biosynthesis of high-value bioactive metabolites named isoprostanoids. So far, there have been only a few studies reporting isoprostanoid productions in algae. To fill this gap, the current investigation aimed at profiling isoprostanoids by liquid chromatography -mass spectrometry/mass spectrometry (LC-MS/MS) in four marine microalgae. A good correlation was observed between the most abundant polyunsaturated fatty acids (PUFAs) produced by the investigated microalgal species and their isoprostanoid profiles. No significant variations in the content of oxidized derivatives were observed for *Rhodomonas salina* and *Chaetoceros gracilis* under copper stress, whereas increases in the production of C18-, C20- and C22-derived isoprostanoids were monitored in *Tisochrysis lutea* and *Phaeodactylum tricornutum*. In the presence of hydrogen peroxide, no significant changes were observed for *C. gracilis* and for *T. lutea*, while variations were monitored for the other two algae. This study paves the way to further studying the physiological roles of isoprostanoids in marine microalgae and exploring these organisms as bioresources for isoprostanoid production.

**Keywords:** microalgae; PUFAs; isoprostanoids; oxidative stress; micro-LC-MS/MS

## 1. Introduction

Marine ecosystems account for approximately half of the global primary production, and unicellular eukaryotes, e.g., photosynthetic microalgae, as part of the phytoplankton, are major contributors to this ocean productivity [1]. These organisms also play critical roles in the biogeochemical cycle of many chemical elements, including carbon, nitrogen, sulfur, phosphorus and silica. Currently, more than 35,000 species of microalgae have been described, which likely represent only a small part of

<footnote n="1">¹ Institut des Biomolécules Max Mousseron, IBMM, Université de Montpellier, CNRS, ENSCM, Faculté de Pharmacie, 34093 CEDEX 5 Montpellier, France; camille.oger@umontpellier.fr (C.O.); guillaume.reversat@umontpellier.fr (G.R.); Amandine.rocher@umontpellier.fr (A.R.); bingqing.zhou@umontpellier.fr (B.Z.); amandyne.linares-maurizi@umontpellier.fr (A.L.-M.); alexandre.guy@umontpellier.fr (A.G.); valerie.bultel@umontpellier.fr (V.B.-P.); jean-marie.galano@umontpellier.fr (J.-M.G.); joseph.vercauteren@umontpellier.fr (J.V.); thierry.durand@umontpellier.fr (T.D.)</footnote>

<footnote n="2">² Integrative Biology of Marine Models, LBI2M (Sorbonne Université/CNRS), Station Biologique de Roscoff (SBR), 29680 Roscoff, France; philippe.potin@sb-roscoff.fr (P.P.); thierry.tonon@york.ac.uk (T.T.); catherine.leblanc@sb-roscoff.fr (C.L.)</footnote>

<footnote n="3">³ Centre for Novel Agricultural Products, Department of Biology, University of York, Heslington, YO10 5DD York, UK</footnote>

<footnote>* Correspondence: claire.vigor@umontpellier.fr; Tel.: +33-411-957-513</footnote>
the supposed biodiversity, since their number of species has been estimated to range between 200,000 and 800,000 [2].

Microalgae can grow mostly autotrophically but, also, heterotrophically or mixotrophically according to culture conditions and metabolic capacities. These are related to the different environments inhabited by these organisms, as well as their evolutionary history that shaped their network of biochemical pathways [3]. Some microalgae exhibit high contents in proteins, lipids, sugars and pigments, making them attractive for a number of biotechnological applications. Such potential has been investigated for the bio-based production of a wide range of compounds for the food, feed, energy, agriculture and health sectors [4].

Among the interesting compounds produced by microalgae are the omega-3 (ω-3 or n-3) long-chain polyunsaturated fatty acids (n-3 PUFAs), eicosapentaenoic acid (EPA; C20: 5n-3) and docosahexaenoic acid (DHA; C22:6n-3). While land plants and microalgae can produce the medium-chain α-linolenic acid (ALA; C18:3n-3), only the latter organisms can convert this precursor into EPA and DHA. These three fatty acids are considered as essential in human nutrition, because ALA cannot be synthesized de novo by humans, and the metabolic conversion efficiency of dietary ALA into EPA and DHA is low and insufficient to meet physiological demands [5]. These n-3 PUFAs have been shown to provide significant benefits on human health [6,7], notably in mitigating a number of pathological conditions, including cardiac diseases [8]. They are also important for the healthy development of the neural system [9,10], and as such, they are necessarily included in infant formula. Very recently, it has been proposed that dietary n-3 PUFAs selectively drive the expansion of adipocyte numbers to produce new fat cells and store saturated fatty acids, enabling the homeostasis of healthy fat tissue [11]. At present, marine fishes and fish oils are the main commercial sources of n-3 PUFAs. However, the suitability of these sources of PUFAs for human consumption has been questioned, notably because of biosafety (e.g., contents in heavy metals) and of overfishing. In addition, the current supply of n-3 PUFAs from these traditional sources is insufficient to satisfy human nutritional requirements [12]. Therefore, new sources of n-3 PUFAs have been investigated, such as wild-type and engineered microbes, including microalgae [13], and the extraction of fish oil from genetically modified crops [14].

It is well-established that PUFAs are highly reactive species sensitive to oxidation because of the presence of bis-allylic structures, in which α-hydrogen atoms are easily removed by the action of free radicals. Some of these free radicals, named reactive oxygen species, are produced under oxidative stress (OS) conditions and react with PUFAs to form, spontaneously through enzymatic reactions, oxidized derivatives of PUFAs. All of these oxidized metabolites are grouped under the term oxylipins. Most of the oxylipins studied so far are derived from the enzymatic transformation catalyzed by enzymes such as lipoxygenases or dioxygenases. During the last two decades, it has been shown that the nonenzymatic oxidation of PUFAs (NEO-PUFAs) leads to other valuable compounds. ALA are precursors of phytoprostanes (PhytoPs), arachidonic acid (AA; C20: 4n-6) of isoprostanes from the serie 2 (IsoPs, serie 2), EPA of isoprostanes from the serie 3 (IsoPs, serie 3), AdA (docosatetraenoic acid; C22: 4n-6) of dihomo-isoprostanes and dihomo-isofurans, and DPA<sub>n-6</sub> (docosapentaenoic acid; C22: 5n-6) and DHA of neuroprostanes (NeuroPs) (Figures 1 and 2). These NEO-PUFAs are considered to be very good markers of OS in plants and animals. They have also been shown to act as lipid mediators, with key functions in various cell-signaling pathways [15], and have been suggested to be potentially beneficial for human health [16].
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Structure of some isoprostanoids isomers derived from n-3 polyunsaturated fatty acids (PUFAs): ALA (α-linolenic acid), EPA (eicosapentaenoic acid) and DHA (docosahexaenoic acid).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Structure of some isoprostanoids isomers derived from n-6 PUFAs: AA (arachidonic acid), DPA<sub>n-6</sub> (docosapentaenoic acid) and AdA (adrenic acid).</figcaption></figure>

Recently [17], we have investigated variations in the isoprostanoid contents of red and brown macroalgae after exposure to oxidative (heavy metal) stress conditions [18]. In addition, changes in the production of isoprostanoids derived from C18, C20 and C22 fatty acids were observed in the microalga *Phaeodatylum tricornutum* subjected to oxidative stress by cultivation under increasing doses
of hydrogen peroxide (H₂O₂) [19]. This work suggested that nonenzymatic oxylipins in *P. tricornutum* may be involved in the control of important processes under various physiological and environmental conditions. In view of these findings, and to go further in the study of the potential production of NEO-PUFAs by algae, we were first interested in increasing our knowledge on the distribution of NEO-PUFAs in different lineages of marine microalgae by establishing qualitative and quantitative profiles under laboratory culture growth conditions. Based on previous analyses of fatty acids and lipid compositions in marine microalgae and, notably, the high production of EPA and DHA by some of them [20–22], we decided to select the four following species: the diatoms *Phaeodactylum tricornutum* and *Chaetoceros gracilis* known to exhibit high content of EPA, the haptophyte *Tisochrysis lutea* that has been shown to produce elevated amounts of ALA and DHA and the cryptophyte *Rhodomonas salina* that harbors similar and high levels of EPA and DHA. Our second objective was to assess changes in the isoprostanoid profiles of the selected microalgae under altered physiological conditions in relationship with the exposure to oxidative stress (copper and hydrogen peroxide treatments).

## 2. Materials and Methods

### 2.1. *Chemicals and Reagents*

All the NEO-PUFA analytical standards, as well as the internal standard (IS) mixture (C19-16-F₁t-PhytoP and C21-15-F₂t-IsoP) used to determine the calibration curve ratio, were synthesized according to previously described procedures [23–26]. NEO-PUFA analytical standards were as follows: phytoprostanes (9-L₁-PhytoP, *ent*-9-L₁-PhytoP, *ent*-16-*epi*-16-F₁t-PhytoP, 9-F₁t-PhytoP, 16-F₁t-PhytoP + 9-*epi*-9-F₁t-PhytoP, 16(*RS*)-16-A₁-PhytoP, 16-B₁-PhytoP and *ent*-16-B₁-PhytoP); phytofurans (*ent*-16(*RS*)-9-*epi*-ST-Δ¹⁴-10-PhytoF, *ent*-9(*RS*)-12-*epi*-ST-Δ¹⁰-13-PhytoF and *ent*-16(*RS*)-13-*epi*-ST-Δ¹⁴-9-PhytoF) coming from oxidation of the C18 n-3 ALA; isoprostanes derived from the C20 n-6 AA (15-F₂t-IsoP, 15-*epi*-15-F₂t-IsoP, 5-F₂t-IsoP, 5-*epi*-5-F₂t-IsoP and 5-F₂c-IsoP); isoprostanes coming from the oxidation of C20 n-3 EPA (8-F₃t-IsoP, 8-*epi*-8-F₃t-IsoP, 18-F₃t-IsoP and 18-*epi*-18-F₃t-IsoP); dihomo-isoprostanes and dihomo-isofurans derived from the C22 n-6 AdA (*ent*-7(*RS*)-7-F₂t-dihomo-IsoP and 7(*RS*)-ST-Δ⁸-11-dihomo-IsoF); neuroprostanes coming from the oxidation of C22 n-3 DHA (4(*RS*)-4-F₄t-NeuroP, 10-F₄t-NeuroP, 10-*epi*-10-F₄t-NeuroP, 20-F₄t-NeuroP and 20-*epi*-20-F₄t-NeuroP) and those derived from the oxidation of C22 DPAₙ-6 (4(*RS*)-4-F₃t-NeuroP₃Dₚₐₙ-6). The only exception is 16 (*RS*)-16-A₁-PhytoP that were purchased from Cayman Chemicals (Ann Arbor, MI, USA). Liquid chromatography – mass spectrometry (LC-MS) grade water, methanol, acetonitrile and chloroform were obtained from Fisher Scientific (Loughborough, UK). Hexane (CHROMASOLV for high performance liquid chromatography - HPLC), formic and acetic acid, ammonia and potassium hydroxide (Fluka for mass spectrometry) were purchased from Sigma-Aldrich (Saint Quentin Fallavier, France). Ethyl acetate (HPLC grade) was acquired from VWR (Fontenay-sous-bois, France). Solid-phase extraction (SPE) cartridges Oasis MAX with mixed polymer phase (3 mL, 60 mg) were obtained from Waters (Milford, MA, USA).

### 2.2. *Microalgal Species*

The four microalgae (*Tisochrysis lutea* RCC 1349, *Phaeodactylum tricornutum* RCC 69, *Chaetoceros gracilis* and *Rhodomonas salina* RCC 20) used in this study were obtained from the Roscoff Culture Collection (RCC) and from the EMBRC Roscoff culture facilities for *C. gracilis*. This latter strain is a kind gift from the Experimental Mollusc Hatchery of Ifremer at Argenton (France) and is cultivated for larvae feeding (Robert et al., 2004, https://archimer.ifremer.fr/doc/2004/rapport-1546.pdf).

### 2.3. *Cultivation of Microalgae and Oxidative Stress Treatments*

Microalgae were grown in Conway medium [27], commonly used in aquaculture, at a temperature of 18 °C and under a continuous light intensity of 300 µmoles m⁻² s⁻¹ for biomass production [28,29]. The volume of the culture was gradually brought, by successive subculture in increasing volumes of the medium, to a final volume of 10 L in Nalgene flasks placed under constant aeration. Cells were harvested
after reaching the stationary phase. After centrifugation (5000 rpm for 25 min), the supernatant was removed, and the pellet resuspended in 2 L of 0.45 µm filtered and autoclaved natural seawater (FSW) collected offshore at Roscoff (at a site with no direct chemical influence from the shore) and free of organic matter. After agitation to ease resuspension, cells were spun down again for 25 min at 5000 rpm, then resuspended as described above in a final volume of 900 mL. This suspension was split into nine glass flasks previously washed overnight with 1% HCl to limit the Cu adsorption and rinsed with mqH₂O and FSW. The 100 mL cell volume was brought to 1 L with FSW. Three flasks were considered as the control, three were used for copper stress (Cu(II) as CuCl₂) and three for incubation in the presence of hydrogen peroxide (H₂O₂).

Oxidative stress was triggered by adding Cu(II) as CuCl₂ (Merck, Germany) at a final concentration of 0.3 µM or H₂O₂ at 1 mM. After 24 h of incubation under the conditions described above, cells were harvested as explained in the previous section. Supernatants were discarded, and cells were washed one time using FWS and centrifugated before freezing in liquid nitrogen and freeze-drying. Algal material was stored at -20 °C until analysis.

## 2.4. Preparation of Algal Samples for Lipidomic Analysis

During the preparation of samples for such analysis, we made two important observations. First, we noticed that one sample of *P. tricornutum* obtained under H₂O₂ stress condition contained some water after lyophilization. This sample was not considered further for extraction. In addition, one sample of *T. lutea* obtained after copper stress showed significant differences in color and texture during the extractive process when compared to the other samples. The data acquired in LC-MS/MS for this latter sample showed numerous outliers (Grubbs' statistical test; data not shown), which were discarded for subsequent analysis.

A protocol similar to what was described for our previous work on macroalgae was applied for lipidomic analysis [18]. Freeze-dried microalgal samples were coarsely reduced to powder using, first, a Mixer Mill MM400 (Retsch®) bench top unit. Then 100 mg of powder was added in grinding matrix tubes (lysing matrix D, MP Biochemicals, Illkirch, France) with 25 µL of BHT (butylated hydroxytoluene 1% in water) and 1 mL of MeOH. Tubes were placed in a FastPrep-24 (MP Biochemicals), and samples were ground for 30 s at a speed of 6.5 m/s. Suspensions were transferred into a 15 mL centrifuge tube, and 1 mL of MeOH, 4 µL of IS (1 ng/µL) and 1.5 mL of phosphate buffer (50 mM, pH 2.1, prepared with NaH₂PO₄ and H₃PO₄) saturated in NaCl were added. Tubes were then stirred for 1 h at 20 °C. Subsequently, the mixture was vortexed and centrifuged at 5000 rpm for 5 min at room temperature. The organic phase was recovered in Pyrex tubes, and the solvent was dried under a stream of nitrogen at 40 °C. Afterward, lipids were hydrolyzed with 950 µL of KOH for 30 min at 40 °C. After incubation, 1 mL of formic acid (FA; 40 mM, pH 4.6) was added before running the SPE separation. First, SPE Oasis MAX cartridges were conditioned with 2 mL of MeOH and equilibrated with 2 mL of formic acid (20 mM, pH 4.5). After loading the sample, the cartridges were successively washed with 2 mL of NH₃ (2% (v/v)), 2 mL of a mixture of formic acid (20 mM):MeOH (70:30, v/v), 2 mL of hexane and, finally, 2 mL of a hexane:ethanol:acetic acid (70:29.4:0.6, v/v/v) mixture. Lastly, all samples were evaporated to dryness under a nitrogen flow at 40 °C for 30 min and reconstituted in 100 µL of a mobile phase (solvent A: water with 0.1% (v/v) of formic acid, solvent B: ACN:MeOH, 8:2 v/v with 0.1% (v/v) of formic acid and A:B ratio 83:17) for injection.

## 2.5. Preparation of Samples for Analysis of Extraction Yield and Matrix Effect

Parameters related to extraction yield (EY) and matrix effect (ME) were determined for a better description of microalgal isoprostanoid profiles. To this aim, three sets of samples were prepared. The first one was obtained by addition of 6.4 µL of two different concentrations of a standard mixture (36 PUFA oxidized metabolites at 0.5 and 8 ng/mL) into 100 mg of freeze-dried microalgae at the beginning of the extraction process described above to reach concentrations of 32 and 512 ng/g, respectively. This corresponds to the "pre-spiked samples". For the second set of samples, extraction was
done as explained in the previous section on 100 mg of microalgae up to the elution step. Then, eluates
were spiked with 6.4 µL of the two different concentrations of standard mixture used for the “pre-spiked
samples” and processed to complete the algal sample preparation protocol. These samples were named
“post-spiked samples”. The third set of samples consisted of standard solutions (final concentrations of
32 and 256 ng/mL) prepared in 100 µL of the mobile phase H₂O:ACN:FA (83:17:0.1, v/v/v). All sets
of samples were analyzed using the LC-MS/MS system described below. The EY was calculated as
the percentage difference between peaks areas of standards in pre-spiked and post-spiked samples.
The ME was determined as the percentage difference between peak areas of standard added to the
extracted samples (post-spiked sample) and pure standards diluted into the mobile phase. The ME
and EY were calculated for each isoprostanoid and for each species.

## 2.6. *Micro-LC-MS/MS Analysis*

All LC-MS analyses were carried out using an Eksigent® MicroLC 200 Plus (Eksigent Technologies, CA, USA) on a HALO C18 analytical column (100 * 0.5 mm, 2.7 µm; Eksigent Technologies, CA, USA) kept at 40 °C. The mobile phase consisted of a binary gradient of solvent A (water with 0.1% (v/v) of formic acid) and solvent B (ACN:MeOH, 8: 2, v/v with 0.1% (v/v) formic acid). The elution was performed at a flow rate of 0.03 mL/min using the following gradient profile (min/%B): 0/17, 1.6/17, 2.85/21, 7.3/25, 8.8/28.5, 11/33.3, 15/40, 16.5/95 and 18.9/95 and then returned to the initial conditions. Under these conditions, no sample contamination or sample-to-sample carryover was observed.

Mass spectrometry analyses were performed on an AB SCIEX QTRAP 5500 (Sciex Applied Biosystems, ON, Canada). The ionization source was electrospray (ESI), and it was operated in the negative mode. The source voltage was kept at -4.5 kV, and N₂ was used as the curtain gas. The multiple ion monitoring (MRM) of each compound was predetermined by MS/MS analysis to define the two transitions for quantification (T₁) and specification (T₂) (Table S1). The analysis was conducted by monitoring the precursor ion to the product ion (T₁). Peak detection, integration and quantitative analysis were performed by MultiQuant 3.0 software (Sciex Applied Biosystems). The quantification of the isoprostanoids was based on calibration curves obtained from the analyte to the IS area under the curve ratio. Linear regression of six concentrations of standards mixture (16, 32, 64, 128, 256 and 512 pg/µL) of each standard were calculated. The sensitivity of the method was evaluated through limit of detection (LOD) and limit of quantification (LOQ) parameters, which were defined as the lowest concentration with a signal to noise ratio above 3 and 10, respectively.

## 2.7. *Statistical Analysis*

All statistical analyses were performed with R [30], all graphics were created with different functions of the tidyverse package [31] and all the tables with kableExtra package [32]. Analyte concentrations were compared by one-way analysis of variance (ANOVA) and post-hoc (Tukey) test for multiple comparison using rstatix package [33]. For all analyses, the significance threshold was 0.05 for the *p*-value resulting from the statistical test used.

# 3. Results

## 3.1. *Analysis of Extraction Yield and Matrix Effect*

The analysis of NEO-PUFAs in natural matrices is extremely challenging, requiring highly sensitive and specific methods for their profiling and characterization. Therefore, a protocol relying on the specific extraction of lipophilic compounds (Folch's extraction), combined with a step of SPE to eliminate potentially interfering substances, was implemented to obtain an extract enriched in NEO-PUFAs. Such a protocol has proven to be efficient for similar analyses in the past [34,35]. Isoprostanoids were subsequently separated, identified and quantified using a micro-LC-MS/MS method validated by previous studies [36–38]. Identification relied on retention times observed during spiked experiments, determination of molecular masses and the analysis of specific MS/MS transitions. Calibration curves
for the calculation of the concentrations were established for 32 compounds (Table S2), as well as LODs and LOQs. Values were found to be dependent of the type of isoprostanoids and ranged between 0.16 and 0.63 pg injected for the LODs and between 0.16 and 1.25 pg injected for the LOQs. In addition, based on previous experiments done on macroalgae (Vigor et al., 2018), we decided to assess the influence of the matrix effect (ME) on the extraction protocol, since this can affect the extraction yields (EY) and/or mass ionization. Therefore, algal samples spiked with two different concentrations of a standard mixture (SM32 or SM256) were analyzed to calculate the EY and the ME, which subsequently enabled the determination of the efficiency of the sample processing (Table S3). The extraction yield, a parameter specific to each compound (standards and IS), allowed the evaluation of product losses that could happen by retention on the SPE cartridge and/or by partial elution during the washing steps. For the majority of analytes of *C. gracilis* and *R. salina*, the apparent loss of compounds during SPE was between 10% and 20%. The results were most often similar for the two spiked concentrations (32 and 512 ng/g). Regarding the type of compounds (PhytoPs, PhytoFs, IsoPs or NeuroPs), no specific trend was noticed. As far as *P. tricornutum* and *T. lutea*, the calculated extraction yield was more than 100% for some analytes, corresponding probably to the coelution of a compound that presents the same MRM transition. Note in the table the values of two or even three units considered to be outliers. To complete this validation, the matrix effect, corresponding to an ion-suppression/enhancement of coeluted matrix compounds, was evaluated. As for EY, ME is specific to each isoprostanoid, and there was no similar behavior across the same class of compounds or across selected species.

For the sake of clarity, results are presented species-by-species in the next sections. In addition, Table 1 provides a summary of the relative percentage distribution of each type of isoprostanoid (ALA, AA and EPA; AdA, EPA, DPA and DHA) in the four species studied.

<table><caption>**Table 1.** Relative percentage distribution of each type of nonenzymatic oxidation-polyunsaturated fatty acids (NEO-PUFAs) in *C. gracilis*, *P. tricornutum*, *T. lutea* and *R. salina*: ALA (α-linolenic acid), AA (arachidonic acid), AdA (adrenic acid), EPA (eicosapentaenoic acid), DPA<sub>n-6</sub> (docosapentaenoic acid) and DHA (docosahexaenoic acid).</caption><thead><tr><th>Microalgal<br>Species</th><th>Metabolites<br>of ALA</th><th>Metabolites<br>of AA</th><th>Metabolites<br>of AdA</th><th>Metabolites<br>of EPA</th><th>Metabolites<br>of DPA</th><th>Metabolites<br>of DHA</th></tr></thead><tbody><tr><td><strong>C. gracilis</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CTL</td><td>0.8%</td><td>6.6%</td><td>0.4%</td><td>89.0%</td><td>0.0%</td><td>3.1%</td></tr><tr><td>Cu<sup>2+</sup></td><td>1.1%</td><td>7.4%</td><td>0.5%</td><td>87.9%</td><td>0.0%</td><td>3.0%</td></tr><tr><td>H<sub>2</sub>O<sub>2</sub></td><td>0.8%</td><td>11.2%</td><td>0.4%</td><td>83.4%</td><td>0.0%</td><td>4.2%</td></tr><tr><td><strong>P. tricornutum</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CTL</td><td>65.5%</td><td>4.0%</td><td>1.7%</td><td>28.2%</td><td>0.0%</td><td>0.6%</td></tr><tr><td>Cu<sup>2+</sup></td><td>58.1%</td><td>4.5%</td><td>2.0%</td><td>34.8%</td><td>0.0%</td><td>0.5%</td></tr><tr><td>H<sub>2</sub>O<sub>2</sub></td><td>44.8%</td><td>4.3%</td><td>1.5%</td><td>48.7%</td><td>0.0%</td><td>0.7%</td></tr><tr><td><strong>T. lutea</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CTL</td><td>69.5%</td><td>0.2%</td><td>0.6%</td><td>0.1%</td><td>1.9%</td><td>27.7%</td></tr><tr><td>Cu<sup>2+</sup></td><td>73.6%</td><td>0.3%</td><td>0.8%</td><td>0.1%</td><td>2.2%</td><td>23.0%</td></tr><tr><td>H<sub>2</sub>O<sub>2</sub></td><td>67.9%</td><td>0.2%</td><td>0.7%</td><td>0.1%</td><td>1.7%</td><td>29.4%</td></tr><tr><td><strong>R. salina</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CTL</td><td>71.2%</td><td>2.4%</td><td>0.5%</td><td>18.2%</td><td>0.1%</td><td>7.4%</td></tr><tr><td>Cu<sup>2+</sup></td><td>66.7%</td><td>2.7%</td><td>0.4%</td><td>21.4%</td><td>0.2%</td><td>8.6%</td></tr><tr><td>H<sub>2</sub>O<sub>2</sub></td><td>79.1%</td><td>3.5%</td><td>1.5%</td><td>12.4%</td><td>0.1%</td><td>3.4%</td></tr></tbody></table>

### 3.2. *Rhodomonas Salina*

Analysis of the isoprostanoid profile of this species revealed the presence of 35 isoprostanoids (Table 2 and Figure S1). The concentrations of metabolites were comprised between 13.4 ng/g for the epimers 4(RS)-4-F<sub>3t</sub>-NeuroP and 2 µg/g for 16-B<sub>1</sub>-PhytoP. The total amount of isoprostanoids in *R. salina* was 10.6 µg/g.
<table><caption>**Table 2.** Quantification of NEO-PUFAs in *R. salina* incubated under control, copper or H₂O₂ stress conditions. Data are mean ± sd (n = 3) expressed as ng/g dry weight. NaN stands for Not a Number, because it is an impossible value.</caption><thead><tr><th rowspan="2">Component Name</th><th colspan="2">CTL</th><th colspan="2">Cu<sup>2+</sup></th><th colspan="2">H<sub>2</sub>O<sub>2</sub></th></tr><tr><th>Conc.</th><th>sd</th><th>Conc.</th><th>sd</th><th>Conc.</th><th>sd</th></tr></thead><tbody><tr><td>10-epi-10-F<sub>4t</sub>-NeuroP</td><td>49.6</td><td>18.1</td><td>81.2</td><td>51.5</td><td>13.2</td><td>5.03</td></tr><tr><td>10-F<sub>4t</sub>-NeuroP</td><td>40.0</td><td>11.9</td><td>62.5</td><td>34.3</td><td>15.1</td><td>5.99</td></tr><tr><td>13-epi-13-F<sub>4t</sub>-NeuroP</td><td>113</td><td>35.7</td><td>160</td><td>76.8</td><td>39.2</td><td>13.4</td></tr><tr><td>13-F<sub>4t</sub>-NeuroP</td><td>183</td><td>65.2</td><td>279</td><td>153</td><td>41.0</td><td>NaN</td></tr><tr><td>14(RS)-14-F<sub>4t</sub>-NeuroP</td><td>51.0</td><td>14.8</td><td>89.0</td><td>55.0</td><td>12.3</td><td>4.78</td></tr><tr><td>15-epi-15-F<sub>2t</sub>-IsoP</td><td>26.0</td><td>7.17</td><td>37.1</td><td>15.4</td><td>31.7</td><td>9.39</td></tr><tr><td>15-F<sub>2t</sub>-IsoP</td><td>14.0</td><td>3.78</td><td>21.9</td><td>10.1</td><td>17.3</td><td>6.42</td></tr><tr><td>16-B<sub>1</sub>-PhytoP</td><td>1960</td><td>96.9</td><td>2190</td><td>208</td><td>1410</td><td>844</td></tr><tr><td>18-F<sub>3t</sub>-IsoP</td><td>711</td><td>239</td><td>1100</td><td>575</td><td>343</td><td>142</td></tr><tr><td>18-epi-18-F<sub>3t</sub>-IsoP</td><td>240</td><td>63.6</td><td>393</td><td>203</td><td>174</td><td>54.9</td></tr><tr><td>20-epi-20-F<sub>4t</sub>-NeuroP</td><td>67.0</td><td>19.7</td><td>97.2</td><td>48.0</td><td>37.3</td><td>16.2</td></tr><tr><td>20-F<sub>4t</sub>-NeuroP</td><td>88.8</td><td>28.7</td><td>143.0</td><td>81.6</td><td>33.1</td><td>11.4</td></tr><tr><td>4(RS)-4-F<sub>3t</sub>-NeuroP</td><td>13.4</td><td>4.57</td><td>22.7</td><td>16.8</td><td>9.24</td><td>1.90</td></tr><tr><td>4(RS)-4-F<sub>4t</sub>-NeuroP</td><td>194</td><td>45.0</td><td>326</td><td>187</td><td>90.1</td><td>9.40</td></tr><tr><td>5-epi-5-F<sub>3t</sub>-IsoP</td><td>457</td><td>137</td><td>717</td><td>378</td><td>278</td><td>77.1</td></tr><tr><td>5(RS)-5-F<sub>2t</sub>-IsoP</td><td>70.1</td><td>17.7</td><td>107</td><td>49.0</td><td>95.7</td><td>22.3</td></tr><tr><td>5-F<sub>3t</sub>-IsoP</td><td>424</td><td>107</td><td>713</td><td>377</td><td>193</td><td>43.3</td></tr><tr><td>5-F<sub>2c</sub>-IsoP</td><td>149</td><td>31.0</td><td>222</td><td>96.2</td><td>143</td><td>34.9</td></tr><tr><td>7(RS)-ST-Δ<sup>18</sup>-11-dihomo-IsoF</td><td>57.1</td><td>3.55</td><td>62.9</td><td>3.56</td><td>125</td><td>22.4</td></tr><tr><td>8-epi-8-F<sub>3t</sub>-IsoP</td><td>38.9</td><td>11.5</td><td>61.4</td><td>33.0</td><td>18.9</td><td>7.81</td></tr><tr><td>8-F<sub>3t</sub>-IsoP</td><td>57.3</td><td>17.0</td><td>82.0</td><td>40.2</td><td>17.0</td><td>4.21</td></tr><tr><td>9-epi-9-F<sub>1t</sub>-PhytoP</td><td>514</td><td>115</td><td>851</td><td>455</td><td>668</td><td>119</td></tr><tr><td>9-F<sub>1t</sub>-PhytoP</td><td>584</td><td>113</td><td>888</td><td>404</td><td>687</td><td>122</td></tr><tr><td>9-L<sub>1</sub>-PhytoP</td><td>1510</td><td>76.6</td><td>1660</td><td>172</td><td>1110</td><td>619</td></tr><tr><td>ent-16-epi-16-F<sub>1t</sub>-PhytoP</td><td>440</td><td>96.7</td><td>715</td><td>362</td><td>550</td><td>118</td></tr><tr><td>ent-16-F<sub>1t</sub>-PhytoP</td><td>311</td><td>66.4</td><td>520</td><td>284</td><td>417</td><td>81.3</td></tr><tr><td>ent-16(RS)-9-epi-ST-Δ<sup>14</sup>-10-PhytoF</td><td>1790</td><td>173</td><td>2180</td><td>411</td><td>1290</td><td>NaN</td></tr><tr><td>ent-9(RS)-12-epi-ST-Δ<sup>10</sup>-13-PhytoF</td><td>434</td><td>NaN</td><td>570</td><td>124</td><td>404</td><td>149</td></tr></tbody></table>

Considering the four PhytoPs corresponding to epimers 9-F₁t-PhytoP, 9-epi-9-F₁t-PhytoP, 16-F₁t-PhytoP and 16-epi-16-F₁t-PhytoP, plus two other derivatives (16-B₁-PhytoP and 9-L₁-PhytoP) and two pairs of phytofuranoid forms (ent-16(RS)-9-epi-ST-Δ¹⁴-10-PhytoF and ent-9(RS)-12-epi-ST-Δ¹⁰-13-PhytoF), ALA is considered as the main source of isoprostanoids in *R. salina*. This is confirmed when assessing the amounts of oxidized derivatives produced for each potential precursor. Those from ALA represent an average value of 7.6 µg/g of algal dry weight, with the 16-B₁-PhytoP and the ent-16(RS)-9-epi-ST-Δ¹⁴-10-PhytoF being the most abundant (2 µg/g and 1.8 µg/g, respectively). DHA was also inferred to produce a wide range of compounds, with up to ten stereoisomeric NeuroPs that could be arranged by pairs. The sum of the DHA derivatives was 0.8 µg/g, i.e., approximately ten times less than the sum of the ALA derivatives. *R. salina* also synthesized six EPA derivatives, again as epimers that could be classified by pairs, which corresponded to an amount of 1.9 µg/g. Therefore, compared to the DHA-derived isoprostanoids, the EPA products had slightly less structural diversity but accumulated at a higher content. Among the other molecular species of interest, it is worth mentioning those derived from AA: five representatives (15-epi-15-F₂t-IsoP, 15-F₂t-IsoP, 5(RS)-5-F₂t-IsoP and 5-F₂c-IsoP) were identified, for a total content of 0.26 µg/g. To complete this description, other isoprostanoids were observed, including AdA derivatives (7(RS)-ST-Δ⁸-11-dihomo-IsoF at the level of 0.06 µg/g) and DPA<sub>n-6</sub> derivatives (4(RS)-4-F<sub>3t</sub>-NeuroP<sub>DPAn-6</sub> at the level of 0.01 µg/g). Based on this analysis, it is interesting to note that, while the cryptophyte *R. salina* is known to produce high amounts of EPA and DHA [39–41], the most abundant isoprostanoids were derived from the C18 ALA (71.5 % of the total amount of PUFA-oxidized derivatives).
After oxidative stress, no modification in the diversity of the molecules identified could be noticed. All the compounds observed under the control condition were still present under the OS condition, regardless of the type of stress applied. Few significant changes were observed in the content of the 35 NEO-PUFAs measured initially. In fact, based on the statistical analysis, the amount of only one compound showed a significant increase (Table S4) between the control condition and H₂O₂ stress. Indeed, it was observed that the content of the two 7(RS)-ST-Δ⁸-11-dihomo-IsoF epimers doubled, from 0.06 µg/g to 0.13 µg/g, after OS (Table 1). In addition, when comparing profiles obtained after copper and peroxide hydrogen additions, differences were noted for three compounds: 8-F₃t-IsoP, ent-16(RS)-9-epi-ST-Δ¹⁴-10-PhytoF and 7(RS)-ST-Δ⁸-11-dihomo-IsoF (Figure 3).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Changes in contents of selected isoprostanoids for the cryptophyte *Rhodomonas salina* between the control condition (CTL) and oxidative stress (Cu²⁺ and H₂O₂) conditions. Statistically relevant responses between the control and stress conditions (one-way ANOVA) are indicated by asterisks: \* *p* < 5 × 10⁻² and \*\* *p* < 5 × 10⁻³; ns, not significant.</figcaption></figure>

### 3.3. *Tisochrysis Lutea*

*T. lutea* has a greater diversity of compounds identified compared to the other species investigated, with 38 isoprostanoids observed. This consisted of 16 derivatives coming from ALA oxidation, ten derivatives from DHA, three from EPA, three from AA, four from AdA and, finally, two from DPA<sub>n-6</sub>. Despite this increased diversity, the total amount of isoprostanoids measured in this haptophyte was lower than what was measured in the cryptophyte *R. salina*, i.e., 7 µg/g. The details of the oxidized PUFA derivatives grouped by family are as follow: 4.8 µg/g from ALA, 1.9 µg/g from DHA, 0.008 µg/g from EPA, 0.012 µg/g from AA, 0.045 µg/g from AdA and 0.135 µg/g from DPA<sub>n-6</sub>. The levels of individual metabolites were comprised between 1.24 ng/g for 8-epi-8-F₃t-IsoP and 0.988 µg/g for 16-B₁-PhytoP, which represents a 1000-fold difference and indicated that isoprostanoids can be produced at very different ranges in this microalga species (Table 3).

The production of these molecules is in accordance with the fatty acid profile of *T. lutea*, since ALA and DHA are the most abundant PUFAs measured in this alga [39,42], and the contents of their oxidized derivatives represented 97% of the total amount of identified isoprostanoids. Taking a closer look at the four main families of metabolites, the amounts of ALA derivatives ranged from 0.07 µg/g (ent-9-D<sub>1t</sub>-PhytoP) to 1 µg/g (16-B<sub>1</sub>-PhytoP) and those of AA derivatives from 5 ng/g (5(RS)-5-F<sub>2t</sub>-IsoP) to 6 ng/g (5-F<sub>2c</sub>-IsoP). The contents of the EPA derivatives went from 1.2 ng/g for the lowest (8-epi-8-F₃t-IsoP) to 4.5 ng/g for the highest (18-epi-18-F₃t-IsoP) and, for DHA derivatives, from 0.1 µg/g for the lowest (10-epi-10-F₄t-NeuroP) to 0.5 µg/g for the highest (4(RS)-4-F₄t-NeuroP).
<table><caption>**Table 3.** Quantification of NEO-PUFAs in *T. lutea* incubated under control, copper or H₂O₂ stress conditions. Data are mean ± sd (n = 3, except for copper stress: n = 2) expressed as ng/g dry weight. NaN stands for Not a Number, because it is an impossible value.</caption><thead><tr><th rowspan="2">Component Name</th><th colspan="2">CTL</th><th colspan="2">Cu<sup>2+</sup></th><th colspan="2">H<sub>2</sub>O<sub>2</sub></th></tr><tr><th>Conc.</th><th>sd</th><th>Conc.</th><th>sd</th><th>Conc.</th><th>sd</th></tr></thead><tbody><tr><td>10-epi-10-F<sub>4t</sub>-NeuroP</td><td>137</td><td>12.4</td><td>190</td><td>NaN</td><td>175</td><td>51.7</td></tr><tr><td>10-F<sub>4t</sub>-NeuroP</td><td>105</td><td>10.3</td><td>153</td><td>NaN</td><td>132</td><td>37.5</td></tr><tr><td>13-epi-13-F<sub>4t</sub>-NeuroP</td><td>282</td><td>15.4</td><td>312</td><td>NaN</td><td>334</td><td>71.2</td></tr><tr><td>13-F<sub>4t</sub>-NeuroP</td><td>433</td><td>33.3</td><td>505</td><td>NaN</td><td>543</td><td>125</td></tr><tr><td>14(RS)-14-F<sub>4t</sub>-NeuroP</td><td>132</td><td>15.8</td><td>159</td><td>NaN</td><td>167</td><td>72.9</td></tr><tr><td>16-B<sub>1</sub>-PhytoP</td><td>988</td><td>99.0</td><td>1010</td><td>NaN</td><td>1050</td><td>169</td></tr><tr><td>16(RS)-16-A<sub>1</sub>-PhytoP</td><td>324</td><td>30.2</td><td>587</td><td>NaN</td><td>348</td><td>39.3</td></tr><tr><td>18-F<sub>3t</sub>-IsoP</td><td>4.54</td><td>1.13</td><td>4.60</td><td>NaN</td><td>2.46</td><td>NaN</td></tr><tr><td>20-epi-20-F<sub>4t</sub>-NeuroP</td><td>116</td><td>15.3</td><td>213</td><td>NaN</td><td>144</td><td>33.2</td></tr><tr><td>20-F<sub>4t</sub>-NeuroP</td><td>203</td><td>20.8</td><td>326</td><td>NaN</td><td>256</td><td>60.0</td></tr><tr><td>4(RS)-4-F<sub>3t</sub>-NeuroP</td><td>135</td><td>17.5</td><td>261</td><td>NaN</td><td>138</td><td>15.3</td></tr><tr><td>4(RS)-4-F<sub>4t</sub>-NeuroP</td><td>515</td><td>38.7</td><td>826</td><td>NaN</td><td>601</td><td>102</td></tr><tr><td>5(RS)-5-F<sub>2t</sub>-IsoP</td><td>5.05</td><td>0.635</td><td>13.1</td><td>NaN</td><td>6.10</td><td>1.09</td></tr><tr><td>5-F<sub>2c</sub>-IsoP</td><td>6.33</td><td>NaN</td><td>19.6</td><td>NaN</td><td>9.03</td><td>0.925</td></tr><tr><td>7(RS)-ST-Δ<sup>18</sup>-11-dihomo-IsoF</td><td>41.8</td><td>3.50</td><td>78.5</td><td>NaN</td><td>49.5</td><td>12.2</td></tr><tr><td>8-epi-8-F<sub>3t</sub>-IsoP</td><td>1.24</td><td>0.344</td><td>2.43</td><td>NaN</td><td>1.90</td><td>0.939</td></tr><tr><td>8-F<sub>3t</sub>-IsoP</td><td>2.27</td><td>0.448</td><td>4.21</td><td>NaN</td><td>3.64</td><td>0.97</td></tr><tr><td>9-epi-9-F<sub>1t</sub>-PhytoP</td><td>237</td><td>43.4</td><td>575</td><td>NaN</td><td>324</td><td>71.2</td></tr><tr><td>9-F<sub>1t</sub>-PhytoP</td><td>407</td><td>33.6</td><td>730</td><td>NaN</td><td>466</td><td>65.9</td></tr><tr><td>9-L<sub>1</sub>-PhytoP</td><td>727</td><td>82.2</td><td>1300</td><td>NaN</td><td>759</td><td>128</td></tr><tr><td>ent-16-epi-16-F<sub>1t</sub>-PhytoP</td><td>315</td><td>28.8</td><td>610</td><td>NaN</td><td>361</td><td>50.7</td></tr><tr><td>ent-16-F<sub>1t</sub>-PhytoP</td><td>381</td><td>20.4</td><td>682</td><td>NaN</td><td>409</td><td>66.6</td></tr><tr><td>ent-16(RS)-13-epi-Δ<sup>14</sup>-9-PhytoF</td><td>188</td><td>18.6</td><td>402</td><td>NaN</td><td>210</td><td>21.6</td></tr><tr><td>ent-16(RS)-9-epi-ST-Δ<sup>14</sup>-10-PhytoF</td><td>859</td><td>84.5</td><td>1810</td><td>NaN</td><td>1020</td><td>107</td></tr><tr><td>ent-7(RS)-7-F<sub>2t</sub>-dihomo-IsoP</td><td>3.26</td><td>0.602</td><td>13.7</td><td>NaN</td><td>4.47</td><td>1.48</td></tr><tr><td>ent-9-D<sub>1t</sub>-PhytoP</td><td>65.5</td><td>10.9</td><td>172</td><td>NaN</td><td>94.6</td><td>35.6</td></tr><tr><td>ent-9-epi-9-D<sub>1t</sub>-PhytoP</td><td>115</td><td>24.9</td><td>263</td><td>NaN</td><td>142</td><td>38.1</td></tr><tr><td>ent-9(RS)-12-epi-ST-Δ<sup>10</sup>-13-PhytoF</td><td>221</td><td>18.9</td><td>451</td><td>NaN</td><td>251</td><td>27.4</td></tr></tbody></table>

When assessing the impact of oxidative conditions, the oxidized metabolite diversity remained unchanged. Furthermore, no variation in the isoprostanoid content of *T. lutea* was monitored after incubation in the presence of H₂O₂. In contrast, the copper treatment had a strong effect; the contents of 17 among the 38 oxidized derivatives increased under this stress condition. Two compounds were very significantly impacted, as shown by the *p*-values adjusted for multiple comparisons: 5-F₂c-IsoP (*p* < 0.0005) and ent-16-epi-16-F₁t-PhytoP (*p* < 0.001). To a lesser extent, changes in contents of ent-16(RS)-9-epi-ST-Δ¹⁴-10-PhytoF, ent-9(RS)-12-epi-ST-Δ¹⁰-13-PhytoF and 4(RS)-4-F₃t-NeuroP<sub>DPAn-6</sub> were also statistically supported (*p* < 0.005) (Table S4). This concerned four of the six families of NEO-PUFAs identified, i.e., those derived from ALA, EPA, DPA<sub>n-6</sub> and AdA, while the contents of the derivatives from EPA and DHA did not change significantly. Comparing the contents of each of the four families between the control and copper stress condition indicated an increase by 160%, 184%, 104% and 93% for the derivatives of ALA, AA, AdA and for DPA<sub>n-6</sub>, respectively, i.e., a two to three-fold increase in isoprostanoid contents (Figure 4). The concentration of ent-9-D₁t-PhytoP, the metabolite from ALA with the highest content in the control condition (65 ng/g), reached a value of 172 ng/g after cupric stress. Similarly, the content of 5-F₂c-IsoP increased from 6 ng/g to 20 ng/g, of ent-7(RS)-F₂t-dihomo-IsoP from 3 ng/g to 14 ng/g and, of 4(RS)-4-F₃t-NeuroP<sub>DPAn-6</sub>, from 135 ng/g to 260 ng/g. None of the inventoried compounds shown a decrease in contents after applying any of the two oxidative stresses.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Changes in contents of selected isoprostanoids for the haptophyte *Tisochrysis lutea* between the control condition (CTL) and oxidative stress (Cu²⁺ and H₂O₂) conditions. Statistically relevant responses between the control and stress conditions (one-way ANOVA) are indicated by asterisks: \* *p* < 5 × 10⁻², \*\* *p* < 5 × 10⁻³ and \*\*\* *p* < 5 × 10⁻⁴; ns, not significant.</figcaption></figure>
### 3.4. Chaetoceros Gracilis

In this diatom, 28 different isoprostanoids were identified, derived from all the PUFAs mentioned above, except DPA<sub>n-6</sub>, and accounted for 2.45 µg/g (Table 4).

<table><caption>**Table 4.** Quantification of NEO-PUFAs in *C. gracilis* incubated under control, copper or H<sub>2</sub>O<sub>2</sub> stress. Data are mean ± sd (n = 3) expressed as ng/g dry weight. NaN stands for Not a Number, because it is an impossible value.</caption><thead><tr><th rowspan="2">Component Name</th><th colspan="2">CTL</th><th colspan="2">Cu<sup>2+</sup></th><th colspan="2">H<sub>2</sub>O<sub>2</sub></th></tr><tr><th>Conc.</th><th>sd</th><th>Conc.</th><th>sd</th><th>Conc.</th><th>sd</th></tr></thead><tbody><tr><td>10-epi-10-F<sub>4t</sub>-NeuroP</td><td>5.24</td><td>0.437</td><td>6.23</td><td>1.32</td><td>5.68</td><td>NaN</td></tr><tr><td>10-F<sub>4t</sub>-NeuroP</td><td>3.28</td><td>0.217</td><td>3.97</td><td>0.503</td><td>8.23</td><td>8.04</td></tr><tr><td>13-epi-13-F<sub>4t</sub>-NeuroP</td><td>11.0</td><td>0.779</td><td>14.3</td><td>1.51</td><td>24.7</td><td>22.8</td></tr><tr><td>13-F<sub>4t</sub>-NeuroP</td><td>14.6</td><td>1.88</td><td>16.6</td><td>3.60</td><td>31.8</td><td>25.8</td></tr><tr><td>15-epi-15-F<sub>2t</sub>-IsoP</td><td>13.2</td><td>0.443</td><td>16</td><td>0.24</td><td>19.3</td><td>14.1</td></tr><tr><td>15-F<sub>2t</sub>-IsoP</td><td>9.32</td><td>0.449</td><td>10.8</td><td>1.24</td><td>14.7</td><td>11.0</td></tr><tr><td>16-B<sub>1</sub>-PhytoP</td><td>4.38</td><td>0.913</td><td>7.04</td><td>2.10</td><td>4.22</td><td>1.29</td></tr><tr><td>18-F<sub>3t</sub>-IsoP</td><td>635</td><td>33.5</td><td>767</td><td>NaN</td><td>1040</td><td>792</td></tr><tr><td>18-epi-18-F<sub>3t</sub>-IsoP</td><td>362</td><td>16.1</td><td>428</td><td>NaN</td><td>335</td><td>NaN</td></tr><tr><td>20-epi-20-F<sub>4t</sub>-NeuroP</td><td>9.38</td><td>0.976</td><td>10.3</td><td>NaN</td><td>15.8</td><td>11.5</td></tr><tr><td>20-F<sub>4t</sub>-NeuroP</td><td>10.7</td><td>2.15</td><td>10.3</td><td>2.79</td><td>17.7</td><td>13.4</td></tr><tr><td>4(RS)-4-F<sub>4t</sub>-NeuroP</td><td>22.7</td><td>1.62</td><td>23.2</td><td>4.11</td><td>23.9</td><td>NaN</td></tr><tr><td>5-epi-5-F<sub>3t</sub>-IsoP</td><td>603</td><td>26.6</td><td>661</td><td>92.3</td><td>569</td><td>NaN</td></tr><tr><td>5(RS)-5-F<sub>2t</sub>-IsoP</td><td>35.1</td><td>3.19</td><td>41.4</td><td>7.19</td><td>61.1</td><td>52.8</td></tr><tr><td>5-F<sub>3t</sub>-IsoP</td><td>471</td><td>22.2</td><td>477</td><td>80.4</td><td>493</td><td>NaN</td></tr><tr><td>5-F<sub>2c</sub>-IsoP</td><td>105</td><td>3.71</td><td>140</td><td>20.1</td><td>248</td><td>180</td></tr><tr><td>7(RS)-ST-Δ<sup>18</sup>-11-dihomo-IsoF</td><td>10.8</td><td>0.754</td><td>14.2</td><td>0.802</td><td>12.9</td><td>2.66</td></tr><tr><td>8-epi-8-F<sub>3t</sub>-IsoP</td><td>57.5</td><td>3.01</td><td>63.7</td><td>7.63</td><td>56.0</td><td>NaN</td></tr><tr><td>8-F<sub>3t</sub>-IsoP</td><td>52.4</td><td>2.05</td><td>62.9</td><td>7.24</td><td>56.6</td><td>NaN</td></tr><tr><td>9-F<sub>1t</sub>-PhytoP</td><td>2.06</td><td>0.136</td><td>2.43</td><td>0.0889</td><td>2.83</td><td>1.29</td></tr><tr><td>9-L<sub>1</sub>-PhytoP</td><td>3.21</td><td>0.587</td><td>5.40</td><td>1.53</td><td>3.27</td><td>0.92</td></tr><tr><td>ent-16-epi-16-F<sub>1t</sub>-PhytoP</td><td>1.44</td><td>0.164</td><td>1.84</td><td>0.229</td><td>2.04</td><td>0.952</td></tr><tr><td>ent-16(RS)-9-epi-ST-Δ<sup>14</sup>-10-PhytoF</td><td>4.44</td><td>0.354</td><td>7.08</td><td>1.33</td><td>5.51</td><td>1.59</td></tr><tr><td>ent-9-epi-9-D<sub>1t</sub>-PhytoP</td><td>3.96</td><td>0.347</td><td>6.51</td><td>1.02</td><td>7.54</td><td>5.75</td></tr></tbody></table>

Seven phytoprostanoids and phytofuranoides derived from ALA were observed and represented 0.019 µg/g. The concentrations ranged from 1.4 ng/g (ent-16-epi-16-F<sub>1t</sub>-PhytoP) to 4.4 ng/g (ent-16(RS)-9-epi-ST-Δ<sup>14</sup>-10-PhytoF). Eight neuroprostanoids from DHA (0.077 µg/g), six EPA derivatives (2.2 µg/g), five AA derivatives (0.16 µg/g) and two AdA derivatives (0.011 µg/g) were also identified. DHA derivatives had contents ranging from 3.3 ng/g (10-F<sub>4t</sub>-NeuroP) to 14.6 ng/g (13B(RS)-13-F<sub>4t</sub>-NeuroP), while EPA derivatives accumulated from 52.4 ng/g (8-epi-8-F<sub>3t</sub>-IsoP) to 535 ng/g (18-epi-18-F<sub>3t</sub>-IsoP) and AA derivatives from 9.3 ng/g (15-F<sub>2t</sub>-IsoP) to 105 ng/g (5-F<sub>2c</sub>-IsoP). As it could be expected based on the high content of EPA found in *C. gracilis* [39,43,44], the most abundant isoprostanoids identified in this species were derived from this PUFA, notably the diasteroisomer pair 5(RS)-5-F<sub>3t</sub>-IsoP that accounted for approximately 42% (1.1 µg/g) of the total amount of oxidized metabolites measured.

For this alga, the qualitative profile remained mostly unchanged, similarly to the amounts of the individual molecules, between the control and the two oxidative stress conditions tested. The only exception was for the compound 18-F<sub>3t</sub>-IsoP; it showed a slight significant difference in its concentration under the Cu<sup>2+</sup> stress condition, increasing from 362 ng/g to 428 ng/g (p = 0.015) (Figure 5 and Table S4).
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Changes in content of selected isoprostanoids for the diatom *Chaetoceros gracilis* between the control condition (CTL) and oxidative stress (Cu²⁺ and H₂O₂) conditions. Statistically relevant responses between the control and stress conditions (one-way ANOVA) are indicated by asterisks: \* *p* < 5 × 10⁻² and \*\* *p* < 5 × 10⁻³; ns, not significant.</figcaption></figure>

### 3.5. *Phaeodactylum Tricornutum*

In this diatom, 21 different oxidized metabolites were identified and quantified, for a total of 0.32 µg/g. No derivatives of DPA<sub>n-6</sub> were observed under any conditions. ALA derivatives represented the main isoprostanoids in term of diversity with 12 metabolites (six PhytoPs and six PhytoFs) and, also, in terms of the content (0.21 µg/g, i.e., 66% of the total amount). Concentrations of ALA-oxidized metabolites were comprised within a range of 2.6 ng/g (with ent-16(RS)-13-epi-ST-Δ<sup>14</sup>-9-PhytoF) to 44.2 ng/g (9-F<sub>1t</sub>-PhytoP). The second-most abundant derivatives were produced from EPA, with four metabolites that accounted for 0.09 µg/g. The dynamic range was from 1.2 ng/g for 8-epi-8-F<sub>3t</sub>-IsoP to 47 ng/g for (5(R)-5-F<sub>3t</sub>-IsoP. One single isoprostanoid from AA was identified (5-F<sub>2c</sub>-IsoP; 13 ng/g), two from AdA (7(RS)-ST-Δ<sup>8</sup>-11-dihomo-IsoF epimers; 5 ng/g) and, finally, two from DHA (4(RS)-4-F<sub>4t</sub>-NeuroP epimers; 2 ng/g). No oxidized derivatives of DPA<sub>n-6</sub> were found, as previously stated for the other diatom, *C. gracilis* (Table 5). As already mentioned for *T. lutea* and *C. gracilis*, *P. tricornutum* produced isoprostanoids in accordance with its PUFA profile that contained mainly ALA and EPA [39,45].

<table><caption>**Table 5.** Quantification of NEO-PUFAs in *P. tricornutum* incubated under control, copper or H₂O₂ stress conditions. Data are mean ± sd (n = 3) expressed as ng/g dry weight. NaN stands for Not a Number, because they are impossible values.</caption><thead><tr><th rowspan="2">Component Name</th><th colspan="2">CTL</th><th colspan="2">Cu<sup>2+</sup></th><th colspan="2">H<sub>2</sub>O<sub>2</sub></th></tr><tr><th>Conc.</th><th>sd</th><th>Conc.</th><th>sd</th><th>Conc.</th><th>sd</th></tr></thead><tbody><tr><td>16-B<sub>1</sub>-PhytoP</td><td>15.1</td><td>1.38</td><td>42.7</td><td>5.94</td><td>3.80</td><td>1.23</td></tr><tr><td>4(RS)-4-F<sub>4t</sub>-NeuroP</td><td>1.84</td><td>0.0457</td><td>2.91</td><td>0.169</td><td>2.31</td><td>0.518</td></tr><tr><td>5-epi-5-F<sub>3t</sub>-IsoP</td><td>47.1</td><td>2.74</td><td>99.9</td><td>14.0</td><td>86.6</td><td>49.4</td></tr><tr><td>5-F<sub>3t</sub>-IsoP</td><td>33.7</td><td>2.07</td><td>72.9</td><td>8.36</td><td>57.9</td><td>32.4</td></tr><tr><td>5-F<sub>2c</sub>-IsoP</td><td>13.0</td><td>1.02</td><td>25.2</td><td>3.05</td><td>14.2</td><td>2.78</td></tr><tr><td>7(RS)-ST-Δ<sup>18</sup>-11-dihomo-IsoF</td><td>5.41</td><td>0.554</td><td>11.1</td><td>1.66</td><td>5.14</td><td>1.69</td></tr><tr><td>8-epi-8-F<sub>3t</sub>-IsoP</td><td>6.72</td><td>0.759</td><td>12.0</td><td>1.83</td><td>10.8</td><td>4.19</td></tr><tr><td>8-F<sub>3t</sub>-IsoP</td><td>3.92</td><td>0.304</td><td>7.99</td><td>0.911</td><td>6.92</td><td>3.48</td></tr><tr><td>9-F<sub>1t</sub>-PhytoP</td><td>44.2</td><td>4.36</td><td>52.5</td><td>3.05</td><td>37.2</td><td>3.93</td></tr><tr><td>9-L<sub>1</sub>-PhytoP</td><td>12.3</td><td>1.13</td><td>34.0</td><td>4.63</td><td>3.25</td><td>0.846</td></tr><tr><td>ent-16-epi-16-F<sub>1t</sub>-PhytoP</td><td>31.5</td><td>3.19</td><td>37.1</td><td>2.10</td><td>26.3</td><td>2.54</td></tr><tr><td>ent-16-F<sub>1t</sub>-PhytoP</td><td>81.0</td><td>8.05</td><td>90.9</td><td>4.69</td><td>67.3</td><td>7.15</td></tr><tr><td>ent-16(RS)-13-epi-Δ<sup>14</sup>-9-PhytoF</td><td>2.64</td><td>0.312</td><td>7.01</td><td>NaN</td><td>1.06</td><td>0.268</td></tr><tr><td>ent-16(RS)-9-epi-ST-Δ<sup>14</sup>-10-PhytoF</td><td>20.2</td><td>2.46</td><td>44.0</td><td>4.61</td><td>5.78</td><td>1.94</td></tr><tr><td>ent-9-epi-9-D<sub>1t</sub>-PhytoP</td><td>0.918</td><td>0.182</td><td>2.40</td><td>0.392</td><td>0.783</td><td>0.426</td></tr><tr><td>ent-9(RS)-12-epi-ST-Δ<sup>10</sup>-13-PhytoF</td><td>4.53</td><td>0.744</td><td>11.7</td><td>1.65</td><td>3.58</td><td>0.882</td></tr></tbody></table>

The isoprostanoid profile of the diatom *P. tricornutum* was strongly influenced by the copper treatment, in contrast to what was observed for the other diatom, *C. gracilis* (Figure 6).
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Changes in contents of selected isoprostanoids for the diatom *Phaeodactylum tricornutum* between the control condition (CTL) and oxidative stress (Cu²⁺ and H₂O₂) conditions. Statistically relevant responses between the control and stress conditions (one-way ANOVA) are indicated by asterisks: \* *p* < 5 × 10⁻², \*\* *p* < 5 × 10⁻³, \*\*\* *p* < 5 × 10⁻⁴ and \*\*\*\* *p* < 1 × 10⁻⁴; ns, not significant.</figcaption></figure>

A significant increase in the content of 14 metabolites among the 21 identified was observed (Table S4). This was particularly significant for *ent*-16(RS)-13-*epi*-ST-Δ¹⁴-9-PhytoF (*p* < 0.00005), with a concentration of 2.6 ng/g in the control condition and 7 ng/g under the copper stress condition. Compounds 16-B₁-PhytoP, 9-L₁-PhytoP, *ent*-16(RS)-9-*epi*-ST-Δ¹⁴-10-PhytoF and *ent*-9(RS)-12-*epi*-ST-Δ¹⁰-13-PhytoF were also strongly impacted by copper, with significant modifications in the contents (*p* < 0.0005). The concentrations of these different metabolites increased from 15 ng/g to 43 ng/g, 12 ng/g to 34 ng/g, 20 ng/g to 44 ng/g and 4.5 ng/g to 12 ng/g, respectively. The amounts
of the derivatives of ALA were enhanced by 94%, of AA by 94%, of AdA by 105% and of DHA by 58%, representing, on average, a two-fold increase in the isoprostanoid contents. No changes in the isoprostanoids for which EPA is the precursor was observed. No metabolite showed a decrease in contents.

Surprisingly, we monitored the lower levels of the isoprostanoids under H₂O₂ stress compared to the control condition for *P. tricornutum*. The contents of the two epimer series derived from ALA oxidation, *ent*-16(RS)-13-*epi*-ST-Δ¹⁴-9-PhytoF and *ent*-16(RS)-9-*epi*-ST-Δ¹⁴-10-PhytoF, decreased by factors of two and four, respectively (p < 0.005). In the same vein, the amounts of 16-B₁-PhytoP and 9-L₁-PhytoP were four times lower under the H₂O₂ stress condition compared to the control (p < 0.05).

**4. Discussion**

There is an increasing interest in studying oxylipin metabolism in marine microalgae. This is supported by recent publications describing profiles of enzymatically produced oxidized derivatives of PUFAs and their potential physiological roles [17,46–49]. So far, little emphasis has been put on biosynthesis, by eukaryotic phytoplankton, of isoprostanoids, i.e., oxylipins produced nonenzymatically by a reaction of ROS with the double-bonds of PUFAs. In this context, and to our knowledge, the current study is the first to report the production of phytoprostanes, phytofurans, isoprostanes (serie 2 and 3) and neuroprostanes, all derived from PUFA precursors that include ALA, AdA, EPA, DPA<sub>n-6</sub> and DHA, in the cryptophyte *R. salina*, the haptophyte *T. lutea* and the diatom *C. gracilis*. In addition, it extends the repertoire of isoprostanoids recently published for another diatom, *P. tricornutun* [19]. Under laboratory culture growth conditions, a good correlation between the presence of PUFAs and the biosynthesis of NEO-PUFAs was observed. We have already noticed this in our previous study on macroalgae, notably with the *Rhodophyta* species *Grateloupia turuturu* Yamada, known to be rich in AA and which produced significant amounts of oxidized AA metabolites [18]. The four microalgae investigated in the present study exhibited different levels of diversity, as well as marked differences in the amounts of isoprostanoids produced. A high content of ALA derivatives was quantified in *R. salina*. The diatom *P. tricornutum*, which appeared to contain lower amounts of oxidized derivatives compared to the other microalgae used, was mostly rich in AA derivatives. The diatoms *C. gracilis* showed profiles rich in NEO-PUFAs produced from EPA. High contents of ALA and DHA derivatives were identified in *T. lutea*. After exposure to oxidative stress conditions, changes in the diversity and amounts of isoprostanoids produced were species and stress-dependent. Under copper stress, no strong variations were observed in *R. salina* and *C. gracilis*, whereas a significant increase in the production of C18-, C20- and C22-derived isoprostanoids was monitored in *T. lutea* and *P. tricornutum*. H₂O₂ stress had different impacts. NEO-PUFA concentrations remained unchanged for *C. gracilis* and *T. lutea*, whereas profiles and contents were altered in *R. salina* and *P. tricornutum*, notably for the ALA-oxidized derivatives. Changes in phytoprostanes derived from ALA have been recently observed in this latter alga under a H₂O₂ treatment slightly different from the condition considered in our analysis (1 mM for 24 h), i.e., 0.25 and 0.75 mM of H₂O₂ applied during 48 h [19]. Interestingly, this study identified a number of isoprostanoids derived from ALA, ARA, EPA and DHA, which levels were differentially affected after oxidative stress. The authors have studied the influence of nine synthetic isoprostanoids, applied in the micromolar range, on the physiology and lipid metabolism of *P. tricornutum*. They observed an induction of the accumulation of triacylglycerols (storage lipids) and a reduction of growth without the alteration of photosynthesis. Such a study, describing the characterization of nonenzymatic oxylipins in *P. tricornutum* and suggesting physiological roles for these molecules, paved the way to better understand their importance in the biology of marine microalgae.

In the context of research on microalgal biorefinery, numerous studies were conducted on culture parameters and have been shown to impact the production of PUFAs, such as light [50], macronutrient depletion [51,52], temperature [53] or salinity [54]. Nitrogen depletion or salinity stress, for instance, were shown to induce oxidative stress and significant changes in PUFA productions [51], but little is known about the impacts of these abiotic parameters on oxidized PUFA derivatives. Under laboratory
culture growth conditions, our study showed that isoprostanoid profiles present good correlations with PUFA contents, and that their production could be increased in *T. lutea*, *P. tricornutum* and *R. salina* by applying direct oxidative stress, either through a copper or H₂O₂ addition. According to previous studies on lipid metabolism regulation, these results suggest that culture condition manipulation could also be an interesting field to be explored for improving the biotechnological production of microalgal isoprostanoids.

From a more methodological point of view, it is worth underlining the sensitivity of the measurements, since we have managed to measure metabolites present in very small quantities. Indeed, considering the case of *P. tricornutum* and of its lowest represented isoprostanoid (8-epi-8-F₃t-IsoP) presents at approximately 1 ng/g, it is satisfying to detect and reliably quantify molecules at such very low levels. With LODs ranging between 0.16 ng/g and 0.63 ng/g and LOQs comprised between 0.16 ng/g and 1.25 pg/g, we can consider the method as sensitive. Interestingly, these low detection limits enabled to detect and quantify a great diversity of metabolites. To our knowledge, the 38 oxidized metabolites detected in *T. lutea* represent the highest diversity of isoprostanoids identified from a given organism so far, including plants and animals. This number of 38 isoprostanoids is close to the number of molecules for which we currently have standards for our targeted lipidomics method (47 metabolites). However, we cannot ensure that these undetected compounds were absent. Some of the missing metabolites may be present in amounts lower than our LODs/LODs or may be produced but masked by other molecules from the matrix. To overcome these two issues, it is possible to foresee improving the sample preparation and, also, the chromatographic procedure by working on the choice of column, solvents and gradient. Furthermore, it is important to emphasize that more than 38 isoprostanoids may be present in the microalgae studied but could not be identified, because our analysis is based on a targeted lipidomics approach and, thus, only detects the metabolites present in the analytical method. Therefore, we think it is important to consider extending the library of isoprostanoid/oxylipin standards through the synthesis of new molecules by chemists, as well as to adopt an untargeted lipidomics method [55,56] to expand the investigation of algal isoprostanoids in the future. Another important observation is that there is a very large difference in the concentrations of NEO-PUFAs in the four species investigated, although this can be partly smoothed out by applying the correction factors of the extraction yield and matrix effect.

Finally, recent studies have shown promising biological activities for PhytoPs, IsoPs and NeuroPs [16,57]. For instance, Minghetti et al. showed the ability of the phytoprostane B₁-PhytoP, through novel mechanisms involving PPAR-γ, to specifically affect immature brain cells, such as neuroblasts and oligodendrocyte progenitors, thereby conferring neuroprotection against oxidant injury and promoting myelination [58]. Duda et al. showed the role, also as lipid mediator, of some phytoprostanes in the immediate effector phase of allergic inflammation [38]. More recently, the work of Lee et al. put forward the hypothesis of the neuroprotective effect of 4-F₄t-NeuroP in cellular and animal models [59]. Early studies on the cardiovascular system demonstrated that AA-oxidized derivatives induced platelet aggregation or showed a hypertensive effect [60,61]. More recent studies showed that IsoPs and NeuroPs have beneficial effects in cardiovascular disease. Indeed, Leguennec et al. revealed that the lipid mediator 4-F₄t-NeuroP derived from the nonenzymatic peroxidation of DHA has an antiarrhythmic effect in ventricular cardiomyocytes and in post-myocardial infarcted mice [62]. They also demonstrated the capability of such derivatives to prevent and protect rat myocardium from reperfusion damages following occlusion (ischemia) [63]. Due to high amounts quantified in some of the tested microalgae, especially after copper exposure, it may be worth exploring these organisms as a potential natural resource for the production of isoprostanoids. The extraction of these NEO-PUFAs from marine microalgae could be an interesting alternative to current productions by complex chemical syntheses, as are macroalgae. In this context, further works should focus on assessing how culture conditions alter the isoprostanoid contents and diversity in selected algae towards enhancing productions for future extractions from natural resources.
## 5. Conclusions

The current investigation aimed at profiling isoprostanoids by micro-LC-MS/MS in selected marine microalgae belonging to different lineages: the cryptophyte *Rhodomonas salina*, the haptophyte *Tisochrysis lutea* and the diatoms *Phaeodactylum tricornutum* and *Chaetoceros gracilis*. To our knowledge, this is the first report of such a wide variety of NEO-PUFAs produced in microalgae. For instance, our analysis allowed the detection of PUFA-oxidized derivatives never reported so far, and we detected no less than 38 different metabolites in *T. lutea*. Our study is also the first to establish a link between significant changes in the isoprostanoid profiles of some selected microalgae and heavy metal stress. It also highlights the impact of hydrogen peroxide stress on NEO-PUFAs in some cases. Based on recent studies showing promising biological activities for NEO-PUFAs and due to high amounts quantified in some of the tested microalgae, further works should focus on assessing how manipulating culture conditions could enhance the production of isoprostanoids in selected species, notably by targeting the PUFA biosynthetic pathways.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2218-273X/10/7/1073/s1: Figure S1: Chromatogram of selected reaction monitoring (SRM) of metabolites detected in *Chaetoceros gracilis*. Figure S2: Distribution of oxidized metabolites (sum of concentrations) classified according to the original PUFAs in the four species studied. Table S1: Selected reaction monitoring (SRM) of the isoprostanes derived from polyunsaturated fatty acids. Table S2: Standards calibration curves. Table S3: Determination of the matrix effect and extraction efficiency for isoprostanoid extractions from the four species studied: (a) *P. tricornutum* and *T. lutea* and (b) *C. gracilis* and *R. salina*. Each data point is the mean of six replicates. Table S4: Statistical results expressed with one-way ANOVA and a post-hoc test analysis for the four species studied. Statistically relevant responses between the control and stress conditions (one-way ANOVA) are indicated by asterisks: \* *p* < 5 × 10⁻², \*\* *p* < 5 × 10⁻³ and \*\*\* *p* < 5 × 10⁻⁴; ns, not significant.

**Author Contributions:** C.V., C.O., T.D., P.P., T.T. and C.L. designed and managed the study. C.V. and T.T. wrote the manuscript. C.O., J.V., T.D., J.-M.G., A.G., V.B.-P., P.P. and C.L. contributed to the manuscript corrections. A.R., B.Z., G.R. and A.L.-M. analyzed the NEO-PUFA contents. C.V., G.R., P.P., T.T. and C.L. conducted the study at the experimental facilities. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work received financial support from the EMBRC France (European Marine Biological Resource Centre) through a call of projects 2015. This project was also partly funded by IDEALG (France: ANR-10-BTBR-04) and by the European Joint Programming Initiative “A Healthy Diet for a Healthy Life” (JPI HDHL—ANR-15-HDHL-0003).

**Acknowledgments:** We thank our colleagues of the Marine Biological Resource Centre at the Station Biologique de Roscoff and, especially, Ronan Garnier for providing technical help to implement the microalgal cultivations and those from Merabomer-Corsaire for access to the sample extraction facilities. We also thank our colleagues of the MAMMA platform.

**Conflicts of Interest:** The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

## References

1. Worden, A.Z.; Follows, M.J.; Giovannoni, S.J.; Wilken, S.; Zimmerman, A.E.; Keeling, P.J. Environmental science. Rethinking the marine carbon cycle: Factoring in the multifarious lifestyles of microbes. *Science* **2015**, *347*, 1257-1259. [CrossRef]

2. Ebenezer, V.; Medlin, L.K.; Ki, J.S. Molecular detection, quantification, and diversity evaluation of microalgae. *Mar. Biotechnol.* **2012**, *14*, 129-142. [CrossRef]

3. Brodie, J.; Chan, C.X.; de Clerck, O.; Cock, J.M.; Coelho, S.M.; Gachon, C.; Grossman, A.R.; Mock, T.; Raven, J.A.; Smith, A.G.; et al. The Algal Revolution. *Trends Plant Sci.* **2017**, *22*, 726-738. [CrossRef]

4. Venkata, M.S.; Hemalatha, M.; Chakraborty, D.; Chatterjee, S.; Ranadheer, P.; Kona, R. Algal biorefinery models with self-sustainable closed loop approach: Trends and prospective for blue-bioeconomy. *Bioresour. Technol.* **2020**, *295*, 122128. [CrossRef] [PubMed]

5. Castro, L.F.; Tocher, D.R.; Monroig, O. Long-chain polyunsaturated fatty acid biosynthesis in chordates: Insights into the evolution of Fads and Elovl gene repertoire. *Prog. Lipid Res.* **2016**, *62*, 25-40. [CrossRef] [PubMed]
6. Calder, P.C. Very long chain omega-3 (n-3) fatty acids and human health. *Eur. J. Lipid Sci. Technol.* **2014**, *116*, 1280-1300. [CrossRef]

7. Calder, P.C. Very long-chain n-3 fatty acids and human health: Fact, fiction and the future. *Proc. Nutr. Soc.* **2018**, *77*, 52-72. [CrossRef] [PubMed]

8. Delgado-Lista, J.; Perez-Martinez, P.; Lopez-Miranda, J.; Perez-Jimenez, F. Long chain omega-3 fatty acids and cardiovascular disease: A systematic review. *Br. J. Nutr.* **2012**, *107*, S201-S213. [CrossRef]

9. Innis, S.M. Dietary (n-3) fatty acids and brain development. *J. Nutr.* **2007**, *137*, 855-859. [CrossRef]

10. Campoy, C.; Escolano-Margarit, M.V.; Anjos, T.; Szajewska, H.; Uauy, R. Omega 3 fatty acids on child growth, visual acuity and neurodevelopment. *Br. J. Nutr.* **2012**, *10*, S85-S106. [CrossRef]

11. Hilgendorf, K.I.; Johnson, C.T.; Mezger, A.; Rice, S.L.; Norris, A.M.; Demeter, J.; Greenleaf, W.J.; Reiter, J.F.; Kopinke, D.; Jackson, P.K. Omega-3 Fatty Acids Activate Ciliary FFAR4 to Control Adipogenesis. *Cell* **2019**, *179*, 1289-1305. [CrossRef] [PubMed]

12. Tocher, D.R.; Betancor, M.B.; Sprague, M.; Olsen, R.E.; Napier, J.A. Omega-3 Long-Chain Polyunsaturated Fatty Acids, EPA and DHA: Bridging the Gap between Supply and Demand. *Nutrients* **2019**, *11*, 89. [CrossRef]

13. Diao, J.; Song, X.; Guo, T.; Wang, F.; Chen, L.; Zhang, W. Cellular engineering strategies toward sustainable omega-3 long chain polyunsaturated fatty acids production: State of the art and perspectives. *Biotechnol. Adv.* **2020**, *40*, 107497. [CrossRef] [PubMed]

14. Napier, J.A.; Olsen, R.-E.; Tocher, D.R. Update on GM canola crops as novel sources of omega-3 fish oils. *Plant. Biotechnol. J.* **2019**, *17*, 703-705. [CrossRef] [PubMed]

15. Jahn, U.; Galano, J.-M.; Durand, T. Beyond Prostaglandins—Chemistry and Biology of Cyclic Oxygenated Metabolites Formed by Free-Radical Pathways from Polyunsaturated Fatty Acids. *Angew. Chem. Int. Ed.* **2008**, *47*, 5894-5955. [CrossRef]

16. Galano, J.M.; Lee, Y.Y.; Oger, C.; Vigor, C.; Vercauteren, J.; Durand, T.; Giera, M.; Lee, J.C. Isoprostanes, neuroprostanes and phytoprostanes: An overview of 25years of research in chemistry and biology. *Prog. Lipid Res.* **2017**, *68*, 83-108. [CrossRef]

17. Barbosa, M.; Collado-Gonzalez, J.; Andrade, P.B.; Ferreres, F.; Valentao, P.; Galano, J.M.; Durand, T.; Gil-Izquierdo, A. Nonenzymatic alpha-Linolenic Acid Derivatives from the Sea: Macroalgae as Novel Sources of Phytoprostanes. *J. Agric. Food Chem.* **2015**, *63*, 6466-6474. [CrossRef]

18. Vigor, C.; Reversat, G.; Rocher, A.; Oger, C.; Galano, J.M.; Vercauteren, J.; Durand, T.; Tonon, T.; Leblanc, C.; Potin, P. Isoprostanoids quantitative profiling of marine red and brown macroalgae. *Food Chem.* **2018**, *268*, 452-462. [CrossRef]

19. Lupette, J.; Jaussaud, A.; Vigor, C.; Oger, C.; Galano, J.M.; Reversat, G.; Vercauteren, J.; Jouhet, J.; Durand, T.; Marechal, E. Non-enzymatic synthesis of bioactive isoprostanoids in the diatom Phaeodactylum following oxidative stress. *Plant. Physiol.* **2018**, *178*, 1344-1357. [CrossRef]

20. Kovac, D.J.; Simeunovic, J.B.; Babic, O.B.; Misan, A.C. Algae in food and feed. *Food Feed Res.* **2013**, *40*, 21-32.

21. Chauton, M.S.; Reitan, K.I.; Norsker, N.H.; Tveterås, R.; Kleivdal, H.T. A techno-economic analysis of industrial production of marine microalgae as a source of EPA and DHA-rich raw material for aquafeed: Research challenges and possibilities. *Aquaculture* **2015**, *436*, 95-103. [CrossRef]

22. Kumari, P.; Kumar, M.; Reddy, C.R.K.; Jha, B. Algal lipids, fatty acids and sterols. In *Functional Ingredients from Algae for Foods and Nutraceuticals*; Dominguez, H., Ed.; Woodhead: Cambridge, UK, 2013; pp. 87-134.

23. Oger, C.; Brinkmann, Y.; Bouazzaoui, S.; Durand, T.; Galano, J.M. Stereocontrolled access to isoprostanes via a bicyclo[3.3.0]octene framework. *Org. Lett.* **2008**, *10*, 5087-5090. [CrossRef] [PubMed]

24. Oger, C.; Bultel-Ponce, V.; Guy, A.; Durand, T.; Galano, J.M. Total synthesis of isoprostanes derived from AdA and EPA. *Eur. J. Org. Chem.* **2012**, *2012*, 2621-2634. [CrossRef]

25. Guy, A.; Oger, C.; Heppekausen, J.; Signorini, C.; de Felice, C.; Fürstner, A.; Durand, T.; Galano, J. Oxygenated metabolites of n-3 polyunsaturated fatty acids as potential oxidative stress biomarkers: Total synthesis of 8-F3t-IsoP, 10-F4t-NeuroP and [D4]-10-F4t-NeuroP. *Chem. Eur. J.* **2014**, *20*, 6374-6380. [CrossRef]

26. Cuyamendous, C.; Leung, K.S.; Durand, T.; Lee, J.C.; Oger, C.; Galano, J.M. Synthesis and discovery of phytopurans: Metabolites of alpha-linolenic acid peroxidation. *Chem. Commun.* **2015**, *51*, 15696-15699. [CrossRef]

27. Lananan, F.; Jusoh, A.; Ali, N.; Lam, S.S.; Endut, A. Effect of Conway Medium and f/2 Medium on the growth of six genera of South China Sea marine microalgae. *Bioresour. Technol.* **2013**, *141*, 75-82. [CrossRef]
28. Heydarizadeh, P.; Boureba, W.; Zahedi, M.; Huang, B.; Moreau, B.; Lukomska, E.; Couzinet-Mossion, A.; Wielgosz-Collin, G.; Martin-Jézéquel, V.; Bougaran, G.; et al. Response of CO2-starved diatom Phaeodactylum tricornutum to light intensity transition. *Philos. Trans. R. Soc. B Biol. Sci.* **2017**, *372*, 20160396. [CrossRef]

29. del Pilar Sánchez-Saavedra, M.; Maeda-Martínez, A.N.; Acosta-Galindo, S. Effect of different light spectra on the growth and biochemical composition of Tisochrysis lutea. *J. Appl. Phycol.* **2016**, *28*, 839–847. [CrossRef]

30. R Core Team. *R: A Language and Environment for Statistical Computing*, R version 3.5.4 Patched; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: https://www.R-project.org/ (accessed on 11 March 2019).

31. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; D'Agostino McGowan, L.; François, R.; Grollemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the tidyverse. *J. Open Source Softw.* **2019**, *4*. [CrossRef]

32. Zhu, H.; Travison, T.; Tsai, T.; Beasley, W.; Xie, Y.; Yu, G.; Laurent, S.; Shepherd, R.; Sidi, Y.; Salzer, B.G.; et al. KableExtra: Construct Complex Table with 'kable' and Pipe Syntax. R package version 1.1.0. Available online: https://CRAN.R-project.org/package=kableExtra (accessed on 1 January 2019).

33. Kassambara, K. Rstatix: Pipe-Friendly Framework for Basic Statistical Tests. R package version 0.4.0. Available online: https://CRAN.R-project.org/package=rstatix (accessed on 1 January 2020).

34. Leung, K.S.; Chen, X.; Zhong, W.; Yu, A.C.; Lee, C.Y. Microbubble-mediated sonoporation amplified lipid peroxidation of Jurkat cells. *Chem. Phys. Lipids* **2014**, *180*, 53–60. [CrossRef]

35. Dupuy, A.; Le Faouder, P.; Vigor, C.; Oger, C.; Galano, J.M.; Dray, C.; Lee, J.C.; Valet, P.; Gladine, C.; Durand, T.; et al. Simultaneous quantitative profiling of 20 isoprostanoids from omega-3 and omega-6 polyunsaturated fatty acids by LC-MS/MS in various biological samples. *Anal. Chim. Acta* **2016**, *921*, 46–58. [CrossRef]

36. Yonny, M.E.; Rodriguez Torresi, A.; Cuyamendous, C.; Reversat, G.; Oger, C.; Galano, J.M.; Durand, T.; Vigor, C.; Nazareno, M.A. Thermal Stress in Melon Plants: Phytoprostanes and Phytofurans as Oxidative Stress Biomarkers and the Effect of Antioxidant Supplementation. *J. Agric. Food Chem.* **2016**, *64*, 8296–8304. [CrossRef]

37. Ruesgas-Ramón, M.; Figueroa-Espinoza, M.; Durand, E.; Suárez-Quiroz, M.; González-Ríos, O.; Rocher, A.; Reversat, G.; Vercauteren, J.; Oger, C.; Galano, J.; et al. Identification and quantification of phytoprostanes and phytofurans of coffee and cocoa by- and co-products. *Food Funct.* **2019**, *10*, 6882–6891. [CrossRef] [PubMed]

38. Gonzalez Roldan, N.; Engel, R.; Dupow, S.; Jakob, K.; Koops, F.; Orinska, Z.; Vigor, C.; Oger, C.; Galano, J.M.; Durand, T.; et al. Lipid Mediators From Timothy Grass Pollen Contribute to the Effector Phase of Allergy and Prime Dendritic Cells for Glycolipid Presentation. *Front. Immunol.* **2019**, *10*, 974. [CrossRef] [PubMed]

39. Jónasdóttir, S.H. Fatty Acid Profiles and Production in Marine Phytoplankton. *Mar. Drugs* **2019**, *17*, 151. [CrossRef] [PubMed]

40. Fernández-Reiriz, M.J.; Perez-Camacho, A.; Ferreiro, M.J.; Blanco, J.; Planas, M.; Campos, M.J.; Labarta, U. Biomass production and variation in the biochemical profile (total protein, carbohydrates, RNA, lipids and fatty acids) of seven species of marine microalgae. *Aquaculture* **1989**, *83*, 17–37. [CrossRef]

41. Guevara, M.; Bastardo L Fau-Cortez, R.; Cortez R Fau-Aredondo-Vega, B.; Arredondo-Vega B Fau-Romero, L.; Romero L Fau-Gómez, P.; Gómez, P. Rhodomonas salina (Cryptophyta) pastes as feed for Brachionus plicatilis (Rotifera). *Rev. Biol. Trop.* **2011**, *59*, 1503–1515.

42. Huang, B.; Marchand, J.; Thiriet-Rupert, S.; Carrier, G.; Saint-Jean, B.; Lukomska, E.; Moreau, B.; Morant-Manceau, A.; Bougaran, G.; Mimouni, V. Betaine lipid and neutral lipid production under nitrogen or phosphorus limitation in the marine microalga Tisochrysis lutea (Haptophyta). *Algal Res.* **2019**, *40*, 101506. [CrossRef]

43. Volkman, J.K.; Jeffrey, S.W.; Nichols, P.D.; Rogers, G.I.; Garland, C.D. Fatty acid and lipid composition of 10 species of microalgae used in mariculture. *J. Exp. Mar. Biol. Ecol.* **1989**, *128*, 219–240. [CrossRef]

44. Pratiwi, A.R.; Syah, D.; Hardjito, L.; Panggabean, L.M.G.; Suhartono, M.T. Fatty Acid Synthesis by Indonesian Marine Diatom, Chaetoceros gracilis. *HAYATI J. Biosci.* **2009**, *16*, 151–156. [CrossRef]

45. Cui, Y.; Thomas-Hall, S.R.; Schenk, P.M. Phaeodactylum tricornutum microalgae as a rich source of omega-3 oil: Progress in lipid induction techniques towards industry adoption. *Food Chem.* **2019**, *297*, 124937. [CrossRef] [PubMed]

46. Di Dato, V.; Orefice, I.; Amato, A.; Fontanarosa, C.; Amoresano, A.; Cutignano, A.; Ianora, A.; Romano, G. Animal-like prostaglandins in marine microalgae. *ISME J.* **2017**, *11*, 1722–1726. [CrossRef] [PubMed]
47. Di Dato, V.; Barbarinaldi, R.; Amato, A.; Di Costanzo, F.; Fontanarosa, C.; Perna, A.; Amoresano, A.; Esposito, F.; Cutignano, A.; Lanora, A.; et al. Variation in prostaglandin metabolism during growth of the diatom Thalassiosira rotula. *Sci. Rep.* **2020**, *10*, 5374. [CrossRef] [PubMed]

48. Meyer, N.; Rettner, J.; Werner, M.; Werz, O.; Pohnert, G. Algal Oxylipins Mediate the Resistance of Diatoms against Algicidal Bacteria. *Mar. Drugs* **2018**, *16*, 486. [CrossRef] [PubMed]

49. Russo, E.; d'Ippolito, G.; Fontana, A.; Sarno, D.; D'Alelio, D.; Busseni, G.; Ianora, A.; von Elert, E.; Carotenuto, Y. Density-dependent oxylipin production in natural diatom communities: Possible implications for plankton dynamics. *ISME J.* **2020**, *14*, 164–177. [CrossRef]

50. He, Q.; Yang, H.; Wu, L.; Hu, C. Effect of light intensity on physiological changes, carbon allocation and neutral lipid accumulation in oleaginous microalgae. *Bioresour. Technol.* **2015**, *191*, 219–228. [CrossRef]

51. Zhang, Y.-M.; Chen, H.; He, C.-L.; Wang, Q. Nitrogen Starvation Induced Oxidative Stress in an Oil-Producing Green Alga Chlorella sorokiniana C3. *PLoS ONE* **2013**, *8*, e69225. [CrossRef]

52. Yilancioglu, K.; Cokol, M.; Pastirmaci, I.; Erman, B.; Cetiner, S. Oxidative Stress Is a Mediator for Increased Lipid Accumulation in a Newly Isolated Dunaliella salina Strain. *PLoS ONE* **2014**, *9*, e91957. [CrossRef]

53. Renaud, S.M.; Thinh, L.-V.; Lambrinidis, G.; Parry, D.L. Effect of temperature on growth, chemical composition and fatty acid composition of tropical Australian microalgae grown in batch cultures. *Aquaculture* **2002**, *211*, 195–214. [CrossRef]

54. Church, J.; Hwang, J.H.; Kim, K.T.; McLean, R.; Oh, Y.K.; Nam, B.; Joo, J.C.; Lee, W.H. Effect of salt type and concentration on the growth and lipid content of Chlorella vulgaris in synthetic saline wastewater for biofuel production. *Bioresour. Technol.* **2017**, *243*, 147–153. [CrossRef]

55. Watrous, J.D.; Niiranen, T.J.; Lagerborg, K.A.; Henglin, M.; Xu, Y.-J.; Rong, J.; Sharma, S.; Vasan, R.S.; Larson, M.G.; Armando, A.; et al. Directed Non-targeted Mass Spectrometry and Chemical Networking for Discovery of Eicosanoids and Related Oxylipins. *Cell Chem. Biol.* **2019**, *26*, 433–442. [CrossRef] [PubMed]

56. Martínez, E.; Cosnahan, R.K.; Wu, M.; Gadila, S.K.; Quick, E.B.; Mobley, J.A.; Campos-Gómez, J. Oxylipins mediate cell-to-cell communication in Pseudomonas aeruginosa. *Commun. Biol.* **2019**, *2*, 66. [CrossRef] [PubMed]

57. Wang, W.; Yang, H.; Johnson, D.; Gensler, C.; Decker, E.; Zhang, G. Chemistry and biology of ω-3 PUFA peroxidation-derived compounds. *Prostag. Other Lipid Mediat.* **2017**, *132*, 84–91. [CrossRef] [PubMed]

58. Minghetti, L.; Salvi, R.; Lavinia Salvatori, M.; Antonietta Ajmone-Cat, M.; De Nuccio, C.; Visentin, S.; Bultel-Poncé, V.; Oger, C.; Guy, A.; Galano, J.-M.; et al. Nonenzymatic oxygenated metabolites of α-linolenic acid B1- and L1-phytoprostanes protect immature neurons from oxidant injury and promote differentiation of oligodendrocyte progenitors through PPAR-γ activation. *Free Radic. Biol. Med.* **2014**, *73*, 41–50. [CrossRef] [PubMed]

59. Lee, Y.Y.; Galano, J.-M.; Leung, H.H.; Balas, L.; Oger, C.; Durand, T.; Lee, J.C.-Y. Nonenzymatic oxygenated metabolite of docosahexaenoic acid, 4(RS)-4-F4t-neuroprostane, acts as a bioactive lipid molecule in neuronal cells. *FEBS Lett.* **2020**, *594*, 1797–1808. [CrossRef]

60. Song, W.L.; Paschos, G.; Fries, S.; Reilly, M.P.; Yu, Y.; Rokach, J.; Chang, C.T.; Patel, P.; Lawson, J.A.; Fitzgerald, G.A. Novel eicosapentaenoic acid-derived F3-isoprostanes as biomarkers of lipid peroxidation. *J. Biol. Chem.* **2009**, *284*, 23636–23643. [CrossRef]

61. Praticò, D.; Smyth Em Fau-Violi, F.; Violi F Fau-FitzGerald, G.A.; FitzGerald, G.A. Local amplification of platelet function by 8-Epi prostaglandin F2alpha is not mediated by thromboxane receptor isoforms. *J. Biol. Chem.* **1996**, *271*, 14916–14924. [CrossRef]

62. Roy, J.; Oger, C.; Thireau, J.; Roussel, J.; Mercier-Touzet, O.; Faure, D.; Pinot, E.; Farah, C.; Taber, D.F.; Cristol, J.-P.; et al. Guennec, Nonenzymatic lipid mediators, neuroprostanes, exert the antiarrhythmic properties of docosahexaenoic acid. *Free Radic. Biol. Med.* **2015**, *86*, 269–278. [CrossRef]

63. Roy, J.; Fauconnier, J.; Oger, C.; Farah, C.; Angebault-Prouteau, C.; Thireau, J.; Bideaux, P.; Scheuermann, V.; Bultel-Poncé, V.; Demion, M.; et al. Guennec, Non-enzymatic oxidized metabolite of DHA, 4(RS)-4-F(4t)-neuroprostane protects the heart against reperfusion injury. *Free Radic. Biol. Med.* **2017**, *102*, 229–239. [CrossRef]

<figure><img src="image_10.png" /></figure>