# A 13CO2 Enrichment Experiment to Study the Synthesis Pathways of Polyunsaturated Fatty Acids of the Haptophyte Tisochrysis lutea

Marine Remize, Frédéric Planchon, Matthieu Garnier, Ai Ning Loh, Fabienne Le Grand, Antoine Bideau, Christophe Lambert, Rudolph Corvaisier, Aswani Volety, Philippe Soudant

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

## To cite this version:

Marine Remize, Frédéric Planchon, Matthieu Garnier, Ai Ning Loh, Fabienne Le Grand, et al. A 13CO2 Enrichment Experiment to Study the Synthesis Pathways of Polyunsaturated Fatty Acids of the Haptophyte Tisochrysis lutea. Marine drugs, 2022, 20, pp.22. 10.3390/md20010022. hal-04056518

HAL Id: hal-04056518

https://hal.science/hal-04056518

Submitted on 6 Apr 2023

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_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Article

# A ¹³CO₂ Enrichment Experiment to Study the Synthesis Pathways of Polyunsaturated Fatty Acids of the Haptophyte *Tisochrysis lutea*

Marine Remize <fnref n="1" />,<fnref n="2" />,* ID, Frédéric Planchon ¹ID, Matthieu Garnier <fnref n="3" />, Ai Ning Loh <fnref n="4" />, Fabienne Le Grand ¹ID, Antoine Bideau <fnref n="1" />, Christophe Lambert ¹ID, Rudolph Corvaisier <fnref n="1" />, Aswani Volety <fnref n="5" /> and Philippe Soudant <fnref n="1" />,* ID

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

check for
updates

**Citation:** Remize, M.; Planchon, F.; Garnier, M.; Loh, A.N.; Le Grand, F.; Bideau, A.; Lambert, C.; Corvaisier, R.; Volety, A.; Soudant, P. A ¹³CO₂ Enrichment Experiment to Study the Synthesis Pathways of Polyunsaturated Fatty Acids of the Haptophyte *Tisochrysis lutea*. *Mar. Drugs* **2022**, *20*, 22. https://doi.org/10.3390/md20010022

Academic Editor: Ipek Kurtboke

Received: 2 November 2021

Accepted: 14 December 2021

Published: 24 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The production of polyunsaturated fatty acids (PUFA) in *Tisochrysis lutea* was studied using the gradual incorporation of a ¹³C-enriched isotopic marker, ¹³CO₂, for 24 h during the exponential growth of the algae. The ¹³C enrichment of eleven fatty acids was followed to understand the synthetic pathways the most likely to form the essential polyunsaturated fatty acids 20:5n-3 (EPA) and 22:6n-3 (DHA) in *T. lutea*. The fatty acids 16:0, 18:1n-9 + 18:3n-3, 18:2n-6, and 22:5n-6 were the most enriched in ¹³C. On the contrary, 18:4n-3 and 18:5n-3 were the least enriched in ¹³C after long chain polyunsaturated fatty acids such as 20:5n-3 or 22:5n-3. The algae appeared to use different routes in parallel to form its polyunsaturated fatty acids. The use of the PKS pathway was hypothesized for polyunsaturated fatty acids with n-6 configuration (such as 22:5n-6) but might also exist for n-3 PUFA (especially 20:5n-3). With regard to the conventional n-3 PUFA pathway, Δ6 desaturation of 18:3n-3 appeared to be the most limiting step for *T. lutea*, "stopping" at the synthesis of 18:4n-3 and 18:5n-3. These two fatty acids were hypothesized to not undergo any further reaction of elongation and desaturation after being formed and were therefore considered "end-products". To circumvent this limiting synthetic route, *Tisochrysis lutea* seemed to have developed an alternative route via Δ8 desaturation to produce longer chain fatty acids such as 20:5n-3 and 22:5n-3. 22:6n-3 presented a lower enrichment and appeared to be produced by a combination of different pathways: the conventional n-3 PUFA pathway by desaturation of 22:5n-3, the alternative route of ω-3 desaturase using 22:5n-6 as precursor, and possibly the PKS pathway. In this study, PKS synthesis looked particularly effective for producing long chain polyunsaturated fatty acids. The rate of enrichment of these compounds hypothetically synthesized by PKS is remarkably fast, making undetectable the ¹³C incorporation into their precursors. Finally, we identified a protein cluster gathering PKS sequences of proteins that are hypothesized allowing n-3 PUFA synthesis.

**Keywords:** long-chain PUFA synthesis; desaturases; elongases; PKS pathway; 20:5n-3 (EPA); 22:6n-3 (DHA); *Tisochrysis lutea*; ¹³C artificial enrichment

## 1. Introduction

Long chain polyunsaturated fatty acids (LC-PUFA) such as 20:5n-3 (EPA) and 22:6n-3 (DHA) are important compounds for most marine metazoans for their growth, reproduction, and development. They are not able to synthetize them in sufficient quantities and thus have to acquire them from their diet. On the basis of the food web, protists are the

<footnote n="1">¹ UMR 6539 LEMAR, CNRS, IRD, Ifremer, University of Brest, 29280 Plouzane, France; Frederic.Planchon@univ-brest.fr (F.P.); fabienne.legrand@univ-brest.fr (F.L.G.); Antoine.Bideau@univ-brest.fr (A.B.); christophe.lambert@univ-brest.fr (C.L.); rudolph.corvaisier@univ-brest.fr (R.C.)</footnote>

<footnote n="2">² GREENSEA, Promenade du Sergeant Navarro, 34140 Meze, France</footnote>

<footnote n="3">³ PBA, Ifremer, Rue de l'Ile d'Yeu, BP 21105, CEDEX 03, 44311 Nantes, France; Matthieu.Garnier@ifremer.fr</footnote>

<footnote n="4">⁴ Center for Marine Science, Department of Earth and Ocean Sciences, University of North Carolina Wilmington, 5600 Marvin K. Moss Ln, Wilmington, NC 28403, USA; lohan@uncw.edu</footnote>

<footnote n="5">5 50 Campus Drive, Elon University, Elon, NC 27244, USA; avolety@elon.edu</footnote>

<footnote>* Correspondence: marineremize@greensea.fr (M.R.); philippe.soudant@univ-brest.fr (P.S.)</footnote>
main producers of these fatty acids and present a key role in marine ecosystem functioning. 20:5n-3 and 22:6n-3 are also particularly important in human nutrition. They are known to have beneficial effects on cardiovascular diseases or diabetes. However, due to high demand for human nutrition and aquaculture of carnivore species, a shortage of these two compounds found in fish oil is predicted to occur by 2050 [1]. Despite their economic and ecologic interests, biological and ecological processes responsible for their synthesis are still under investigation. It is, then, of first concern to understand how 20:5n-3 and 22:6n-3 are produced at the basis of the food webs, and how global changes could affect their availability at higher trophic levels.

In phytoplankton and microzooplankton, fatty acids are synthetized via different metabolic pathways [2–4]. The most “conventional” pathway is the fatty acid synthase (FAS) pathway, followed by the elongation and front-end desaturation steps of the n-3 and n-6 pathways. Starting with the initial formation of acetyl-CoA and then malonyl-CoA in aerobic conditions, these pathways produce more complex fatty acids by progressive addition of two atoms of carbon (elongation steps) or desaturations of precursors such as 16:0 or 18:0 [5–7]. These two pathways can be connected by the so-called ω-3 desaturase (or methyl end desaturase) pathway. Within the n-3 and n-6 pathways, an alternative route of Δ8 desaturation can also bypass the Δ6 desaturation step and has already been identified in Haptophyte [8]. These routes allowed the synthesis of 20:5n-3 as well as 22:6n-3 (Figure 1).

<figure><img src="image_6.png" /><figcaption>**Figure 1.** Microalgae fatty acid synthesis pathways. Desaturases are noted with “ΔX” (yellow arrows) and “ωY-des (ΔX)” (blue arrows), where X refers to the location of carbon holding the newly formed double bond from the front end (or carboxyl end) and Y its position from the methyl end. Elo: elongase, FAS: fatty acid synthase.</figcaption></figure>

An alternative O₂ independent pathway, called the polyketide synthase (PKS) pathway, is responsible of long chain PUFA synthesis such as 20:5n-3 and 22:6n-3 [3,9,10]. It has been found in bacteria and protists such as thraustochytrids, dinophytes, and haptophytes [11–15]. The PKS pathway relies on the same four basic enzymatic reactions (condensation, reduction, dehydration, and reduction) as the FAS pathway. Opposed to the conventional pathway, the PKS pathway is less energy consuming, because it requires fewer reduction and dehydration steps than “conventional” pathways [3]. The metabolites used to form the carbon chain are simultaneously desaturated and elongated, creating long-chain PUFA [3,16,17].

Even if some microalgae species share all or part of O₂-dependent n-3 and n-6 pathways and O₂-independent PKS pathways, PUFA composition of primary producers varies greatly according to species. Diatoms synthetize more 20:5n-3 as well as C16 PUFA, while
dinophytes or haptophytes contain more 22:6n-3 or C₁₈ PUFA. Other groups such as cyanobacteria or some chlorophytes classes are unable to build 20:5n-3 or 22:6n-3 or only in very low proportions (<1%) [18–21].

To improve knowledge of the synthesis routes and production of 20:5n-3 and 22:6n-3, studies have focused on identifying genes coding for the different elongase and desaturase [10,22–25]. Moreover, in recent years, the use of ¹³C substrate allowed monitoring the incorporation of labelled substrates into targeted organic macromolecules. Different metabolic intermediates or end products such as fatty acids [26–29] are monitored and quantified. This has already been applied to E. coli [30], yeast [29], and microalgae [31–35]. The development of technics such as gas chromatography coupled to mass spectrometry (GC-c-IRMS) consists of a noticeable improvement in the direct resolution of isotopic composition of organic macromolecules the so-called compounds specific isotope analysis (CSIA) including fatty acids [36–39].

The present study aims at investigating the synthesis pathways of essential PUFA 20:5n-3 and 22:6n-3 using stable isotope (¹³C) labelling experiment of the haptophyte *Tisochrysis lutea*. *T. lutea* is intensively used in aquaculture (hatchery) and industry [40]. The incorporation of ¹³C was monitored in 11 FA during 24 h at a high temporal resolution (each 0.5 to 2 h). Progressive accretion of the ¹³C-labelled CO₂ into FA (from precursors to PUFA of interest) allowed us to constrain FAS, elongase/desaturase, and PKS involvement in 20:5n-3 and 22:6n-3 production by *Tisochrysis lutea*. In parallel to the monitoring of ¹³C incorporation into FA, growth, physiological status, and other cellular parameters (morphology, viability, esterase activity, and lipid content) were monitored by flow cytometry analysis.

## 2. Results

### 2.1. *Algae Physiology and Biochemistry* during the 24 h Experiment

Cell abundance for *Tisochrysis lutea* during the 24 h of experiment increased sharply from $t_0$ to $t_{24}$. The experiment allowed cell concentration to double for the three balloons (Figure 2A). Despite the attention given to homogenization at the time of subculture from inoculum, the second enriched balloon had a cell concentration twice higher than Tl1 and TIT. This difference remained constant during the entire experiment. Cell abundance varied from $4.3 \times 10^6$ cells·mL⁻¹ to $9.5 \times 10^6$ cells·mL⁻¹ for the most concentrated balloon Tl2 and from, on average, $2.6 \times 10^6$ cells·mL⁻¹ to $6.4 \times 10^6$ cells·mL⁻¹ for the two others (Figure 2A). Despite the concentration differences, the general slopes for the three balloons were very similar ($0.22 \times 10^6$ cells·mL⁻¹·h⁻¹ for Tl2 and $0.15 \times 10^6$ cells·mL⁻¹·h⁻¹ for Tl1 + TIT) (Figure 2A). Bacteria were also found in higher abundance Tl2 (Figure 2B), almost five times higher than in Tl1 and TIT. However, bacteria increased only by a factor of 1.2 for Tl2 between $t_0$ and $t_{24}$ versus a factor 3.5 in average for Tl1 and TIT. Bacteria concentration was around 6.6 times higher than algae concentration for Tl1 and TIT on average over the 24 h of the experiment. For Tl2, bacteria concentration was 16 times higher than algae concentration at the beginning of the experiment, but this ratio decreased progressively until $t_{24}$ (8.7 times higher) (Figure 2B).

FSC and SSC were, respectively, considered a proxy of cell size and cell complexity, and FL3 was considered a proxy of chlorophyll content. SSC and red fluorescence (FL3) did not significantly vary during the entire experiment (Bartlett tests, *p* > 0.05 and ANOVA *p* > 0.05) (Table S1). FSC increased slightly with time for the three balloons (Bartlett test, *p* > 0.05 and ANOVA *p* = 0.03) (Table S1). The percentage of dead microalgae (as measured by SYTOX staining assay) remained below 7% for the 24h of the experiment (Table S1).

Particulate organic carbon concentration increased similarly to the cell abundance for the two labelled balloons from 3.6 to 9.0 mmolC·L⁻¹ for Tl1, from 4.5 to 11.6 mmolC·L⁻¹ for Tl2, and from 3.4 to 8.2 mmolC·L⁻¹ for TIT (Figure 3A). Increase in POC, as a function of experiment duration ($R^2$ = 0.81, *p* < 0.0001), occurred at a relatively constant rate of 0.2 mmolC·L⁻¹·h⁻¹ for Tl1, Tl2, and TIT considered together. Total fatty acid (TFA) concentration also increased for the three balloons (between 0.26 mmolC·L⁻¹ to 0.60 mmolC·L⁻¹ for Tl1, from 0.25 mmolC·L⁻¹ to 0.80 mmolC·L⁻¹ for Tl2, and from 0.17 mmolC·L⁻¹
to 0.55 mmolC·L⁻¹ for TIT) (Figure 3B). The slope was 15 µmolC·L⁻¹·h⁻¹ (R² = 0.72, *p* < 0.0001). POC was significantly correlated to cell concentration (R² = 0.80 *p* < 0.0001). The slope of the relation between POC and cell abundance is a proxy of carbon content per cell for *T. lutea*, which was, on average for the three balloons, equal to 1.07 fmolC·cell⁻¹ (Figure 3C). TFA concentration was linearly and positively correlated with POC concentration (R² = 0.73, *p* < 0.0001). The slope of the regression between TFA and POC concentration indicates that TFA represent in average 7.7% of bulk POC (Figure 3D).

<figure><img src="image_7.png" /><figcaption>**Figure 2.** Temporal dynamics of cell concentrations of the two enriched balloons of *Tisochrysis lutea* (TI1 and TI2, filled black circles and filled gray triangles, respectively) and of the control balloon (TIT, empty squares) (A) and corresponding bacteria concentrations (B) during the 24 h experiment.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 3.** Particulate organic carbon (POC) concentration (A) of the two enriched balloons of *Tisochrysis lutea* (TI1 and TI2, filled black circles and filled gray triangles, respectively) and of the control balloon (TIT, empty squares) and total fatty acids (TFA) concentration according to culture age in hours (B). POC concentration according to algae concentration (C). Relation between total fatty acids and POC concentration (D). All regressions (dotted lines) have been calculated with data from the three balloons combined together.</figcaption></figure>
## 2.2. <sup>13</sup>C Atomic Enrichment (AE) of Particulate Organic Carbon and Dissolved Inorganic Carbon

Dissolved inorganic carbon (DIC) progressively enriched in the two balloons with <sup>13</sup>CO<sub>2</sub> (Figure 4A). The enrichment trends were similar for the two balloons after t<sub>4</sub>, with an important increase in the DIC atomic enrichment until t<sub>20</sub> (up to 58.1 and 61.1% for Tl1 and Tl2, respectively). The increase in AE tended to stabilize after t<sub>20</sub>. Final levels of enrichment were 61.5 and 64.6% for Tl1 and Tl2, respectively (Figure 4A). Atomic enrichment (AE<sub>POC</sub>) increased sharply after t<sub>4</sub> for the two balloons until the end of the experiment. Enrichment levels at t<sub>24</sub> were 25.2 and 34.7% for Tl1 and Tl2, respectively (Figure 4B).

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Atomic enrichment of the dissolved inorganic carbon (DIC) (A) and particulate organic carbon (POC) (B) of the two enriched balloons of *Tisochrysis lutea* (Tl1 and Tl2, filled black circles and filled gray triangles, respectively).</figcaption></figure>

## 2.3. Fatty Acid Composition in Neutral and Polar Lipids in *T. lutea*

Neutral lipids and polar lipids represented, respectively, 37% and 63% of TFA on average for the three balloons (Figure 5A). The proportions of individual fatty acid in NL and PL did not vary throughout the experiment. Total bacteria fatty acids (iso15:0, ante15:0, iso16:0, iso17:0, 15:0, 17:0, 21:0, 15:1n-5—Tables S2–S4) remained below 1% for both NL and PL fractions during the 24 h. Branched fatty acids were only present in trace amounts (Tables S2–S4). Concentrations in µg·L<sup>-1</sup> and µmol·C·L<sup>-1</sup> as well as proportions in% of all identified and quantified FA in neutral and polar lipid fractions according to sampling time are available in the Supplementary Files (Tables S2–S4).

We focus the presentation of the results on the polar lipid fraction, as it is the predominant fraction containing FA (Figure 5B). During the experiment, thirty two fatty acids (FA), as listed in the Material and Methods section, were identified and quantified for *T. lutea*, with 12 being over 1% of the TFA in PL (14:0, 16:0, 18:0, 16:1n-7, 18:1n-9, 18:1n-7, 18:2n-6, 18:3n-3, 18:4n-3, 18:5n-3, 22:5n-6, and 22:6n-3) (Figure 5A). Although under 1% for PL, the 16:3n-6, 20:5n-3, and 22:5n-3 were also presented due to their potential synthesis significance (Figure 5A). PUFA (in average 30%) and SFA (21%) were the main FA categories for polar lipids (PL) during the 24 h. PUFA n-3 represented 25% of the TFA, PUFA n-6 5%.

In PL, 14:0 and 22:6n-3 (respectively, 21% and 18% on average over the 24 h) were the most abundant, followed by 18:1n-9, 18:4n-3, and 16:0 (11–13% of TFA). Finally, 18:3n-3, 18:5n-3, and 22:5n-6 ranged from 3 to 4% of the TFA (Figure 5B). Patterns observed for NL are available in Supplementary Files (Figure S1).
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Proportions (%) of NL vs. PL (A) and proportions (%) of fifteen fatty acids in the PL fraction in average over the 24 h (B).</figcaption></figure>

## 2.4. Fatty Acid ¹³C Atomic Enrichment

Figure 6 shows the atomic enrichment (AE) of the eleven fatty acids over time. Despite the different timing and level of enrichment between the two balloons, the temporal dynamic of fatty acids enrichment remained similar. The 18:2n-6 and 18:1n-9 + 18:3n-3 had the highest AE during the entire experiment. 22:5n-6, 16:0, and finally, 20:5n-3 were next. The less enriched fatty acids in the polar lipid fraction were in decreasing order 22:5n-3, 22:6n-3, 18:4n-3, 18:5n-3, and finally, 18:0 (Figure 6). For the NL (Supplementary Files, Figure S2), 20:5n-3, 22:5n-6, 18:1n-9, 16:0, and 18:2n-6 were always the most enriched. The sequence for the other fatty acids remained close to that of polar lipids. It has to be noted that enrichments of 20:5n-3 and, to a lesser extent, of 22:5n-3 were higher in NL than in PL.

Table 1 explored FA synthesis pathways with regard to their most expected direct precursor. Most ratios were below 1, except for the 20:5n-3/18:5n-3 ratio, which was above 1. Similar patterns were observed in NL (Supplementary Files, Table S5).

<table><caption>**Table 1.** Mean ratio of atomic enrichment (AE) for pairs of FA (FAA vs. FAB) in the polar lipids (PL) (mean ± SD, n = 9 sampling dates t₈ to t₂₄) for the two enriched balloons (Tl1, Tl2, Tl = *Tisochrysis lutea*).</caption><thead><tr><th rowspan="3">Fatty Acid B/Fatty Acid A</th><th colspan="4">Polar Lipids</th></tr><tr><th colspan="2">Tl1</th><th colspan="2">Tl2</th></tr><tr><th>Mean *</th><th>SD</th><th>Mean *</th><th>SD</th></tr></thead><tbody><tr><td>18:5n-3/18:4n-3</td><td>0.78</td><td>0.10</td><td>0.78</td><td>0.14</td></tr><tr><td>20:5n-3/18:5n-3</td><td>2.00</td><td>0.92</td><td>1.77</td><td>0.51</td></tr><tr><td>22:5n-3/20:5n-3</td><td>0.98</td><td>0.20</td><td>0.97</td><td>0.06</td></tr><tr><td>22:6n-3/22:5n-3</td><td>0.88</td><td>0.03</td><td>0.87</td><td>0.02</td></tr><tr><td>22:6n-3/20:5n-3</td><td>0.86</td><td>0.17</td><td>0.84</td><td>0.03</td></tr><tr><td>22:6n-3/22:5n-6</td><td>0.61</td><td>0.05</td><td>0.59</td><td>0.06</td></tr></tbody></table>

<footnote>* If the AE of the product (B) exceeds the AE of the reactant (A), ratio > 1, then it is necessary to consider another formation process for B. If the ratio is <1, transformation of A into B is considered possible. If the ratio is close to 1, the fatty acids A and B are at the equilibrium in terms of label incorporated, implying B is then synthesized simultaneously or very rapidly from A.</footnote>
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Atomic enrichment of 11 main fatty acids in the polar lipid (PL) fraction during a 24 h ¹³C labelling experiment of the two enriched balloons of *Tisochrysis lutea* (TI1 and TI2, filled circles and filled triangles, respectively).</figcaption></figure>

## 2.5. Identification of Candidate Proteins for PKS Synthesis in *T. lutea*

Thirty sequences of potential candidate proteins involved in *T. lutea* PUFA synthesis have been identified and are presented in Supplementary File (Table S6). Only fourteen presented the four main domains potentially coding for the enzymes used in PKS PUFA synthesis pathways: ketoacyl reductase (KR), polyketide synthase (KS), dehydrase/dehydrogenase (DH), and enoyl reductase (ER). Among these sequences, four sequences (TISO_14962, TISO_14968, TISO_14975, and TISO_14977) were part of the same cluster (group of homologous proteins) and presented multiple KS, KR, ER, and DH domains as well as phosphopantetheine (PP)-binding domains (Figure 7). TISO_14962 also possessed methyltransferases and thioesterase domains (Figure 7). TISO_14977 presented a domain acknowledged to be involved in acetyl-CoA synthesis. Within this cluster, TISO_14973 was also selected, as it contains an atypical domain, specifically recognized as being involved in n-3 PUFA synthesis. Nine other sequences (TISO_04539, TISO_06404, TISO_06537, TISO_08047, TISO_11097, TISO_16495, TISO_27353, TISO_37260, and TISO_37631) were also found, containing the four main domains (up to 18 for KR in TISO_08047). Except TISO_37631, these sequences also have thioesterase, sulfotransferase, or peptide-synthesis-related domains, and thus they might be in charge of the synthesis of more complex lipids.
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Cluster of candidate proteins suspected involved in PKS PUFA synthesis pathway in *T. lutea*. The name of each protein is annotated with TISO_ (for *Tisochrysis lutea*) and associated number. In the legend, the text written in bold italic correspond to domain names as shown in NCBI conserved domain database, followed by its suspected role.</figcaption></figure>

**3. Discussion**

This study investigated long-chain PUFA synthesis pathways in the haptophyte *Tisochrysis lutea* using the incorporation of ¹³CO₂. Addition of ¹³CO₂ did not affect *T. lutea* physiology. Cell viability remained above 93% during the experiment, while cell complexity and chlorophyll content did not vary significantly according to sampling time. Cell size (as attested by FSC) increased slightly during the 24 h experiment. *T. lutea* produced FA to a level of 7% of POC; predominantly in the form of PL (66%).

Major FA of *T. lutea* were similar in proportions to those found in other prymnesio-
phycea (Haptophytes), i.e., 14:0, 16:0, 18:1n-9, 18:4n-3, and 22:6n-3 [41–45]. As reported
before in Huang et al. (2019) [46], *T. lutea* had a low content of neutral lipids during exponen-
tial phase, and PUFA were mainly found in the polar fraction. *Tisochrysis lutea* accumulates
neutral lipids mainly during stationary phase or under nutritive limitations [45].

The final level of atomic enrichment (AE) into the different FA witnessed active
synthesis, as most fatty acids had a higher AE than that of POC (30% on average for the two
balloons). 22:5n-6 was the most enriched long chain PUFA (LC-PUFA) in the PL fraction.
22:5n-6 and 18:2n-6 were the only ¹³C labelled n-6 fatty acids detectable by GC-c-IRMS.
None of the known synthesis intermediates (18:3n-6, 20:3n-6, 20:4n-6, and 22:4n-6) between
18:2n-6 and 22:5n-6 [4] had measurable ¹³C-labelling and were below 1% in the FA profile
during our experiment. It is then difficult to hypothesize the pathway used to create 22:5n-6
with this missing information. However, even though the different intermediates were
undetectable, 18:2n-6 and 22:5n-6 atomic enrichments being very close cannot exclude them
to be related to each other. While studying the existence of an alternative Δ8 desaturase in Haptophyte, Qi et al. (2002) [8] noticed the absence of intermediates of the n-6 Δ8 desaturase pathway (20:2n-6, 20:3n-6 and 20:4n-6)) in *Isochrysis galbana*. It was attributed to relatively high active enzymes that could form the end-product 22:5n-6 with a rapid flow through these n-6 intermediates. Our results agree with this, as ¹³C enrichment of n-6 intermediates could not be detected by compound specific isotope analysis. To demonstrate the existence of these pathways, it would be interesting to combine functional analysis of desaturases by expression in yeast and GC-c-IRMS monitoring of the intermediates after ¹³C labelling of their precursors.

However, it is also possible that another pathway not involving “classical” n-6 FA intermediates exist in *T. lutea*. Previous studies showed the existence of PKS genes in various species of the prymnesiophytes including *Isochrysis galbana* [47], closely phylogenetically related to *Tisochrysis lutea*. We identified five candidates; proteins potentially involved in PKS synthesis pathway in *T. lutea*. Even if their function has not been verified, it is possible that at least one of the proteins presented in Figure 7 was responsible for the formation of n-6 PUFA in the haptophyte. Thus, our hypothesis is that an n-6 PKS pathway might also exist in *T. lutea* (Figure 8). Finally, PKS and “classical” n-6 routes might not be completely independent and could interact in the synthesis of 22:5n-6 in *T. lutea*.

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

**Figure 8.** Hypothesized pathways for 22:5n-6 synthesis in *T. lutea* in the PL. Numbers in the boxes correspond to final AE value. The triangles symbolize the desaturases (front-end in yellow and methyl-end in purple), the circles the enzymes involved PKS pathway (KR: 3-ketoacyl synthase, KS: 3-ketoacyl-ACP-reductase, DH: dehydrase, 2.2I: 2-trans, 2-cis isomerase, 2.3I: 2-trans, 2-cis isomerase, ER: enoyl reductase), and the squares the elongases.

Despite being one of the most abundant FA, 18:4n-3 showed a low ¹³C-enrichment (23%). The synthesis of 18:4n-3 from 18:3n-3 by Δ6 desaturase had already been described by *Isochrysis* sp. [48]. We assume that such activity also exists in *Tisochrysis*, phylogenetically close to *Isochrysis*. However, as 18:3n-3 co-elute with 18:1n-9, it was not possible to measure its AE and to assess whether this could be a limiting step in n-3 pathway (Figure 9). The 18:5n-3 had the lowest enrichment, and the ratio 18:5n-3/18:4n-3 was below the threshold value (R = 0.78), indicating a feasible transformation of 18:4n-3 into 18:5n-3. The existence of Δ3 desaturase that could support the production of 18:5n-3 (18:5Δ3,6,9,12,15) from 18:4n-3
(18:4Δ6,9,12,15) had been suggested by Joseph (1975) [49] to explain the presence of this unusual FA in dinophytes. A more recent study by Ahman et al. (2011) [23] showed in *Ostreococcus lucimarinus* that a Δ4 desaturase was surprisingly able to add a double bond in 18:4n-3 at the Δ3 position leading to the formation of 18:5n-3 when the gene was expressed in yeast cell and supplemented by 18:4n-3 as substrate. With our results and the discovery of Ahman et al. (2011) [23], we proposed that a Δ4 desaturase of *T. lutea* might be able to act as a Δ3 desaturase on 18:4n-3 to produce 18:5n-3 (Figure 9). Desaturation of 18:4n-3 into 18:5n-3 had been previously hypothesized by Kotajima et al. (2014) [50] in the prymnesiophyte *Emiliania huxleyi*.

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

**Figure 9.** Hypothesized pathways to produce 18:5n-3 in *T. lutea*. Numbers in boxes correspond to final mean AE value, and number in the yellow box the mean value of ratio of the two surrounding fatty acids. The triangles symbolize the desaturases (front-end in yellow and methyl-end in purple), the circles the enzymes involved PKS pathway (KR: 3-ketoacyl synthase, KS: 3-ketoacyl-ACP-reductase, DH: dehydrase, 2.2I: 2-trans, 2-cis isomerase, 2.3I: 2-trans, 2-cis isomerase, ER: enoyl reductase). The directions with dashed arrows cannot be proven with the enrichment dynamics.

The 18:5n-3 was also described as an intermediate of 22:6n-3 synthesis by PKS pathway [4]. However, its low enrichment, as compared to 22:6n-3, appeared not compatible with a hypothetical production through this pathway. Nevertheless, one may speculate that there are two separated PKS pathways, one for the 22:6n-3 and one for the 18:5n-3, as these two PUFA are localized in different cell compartments. The 18:5n-3 is generally associated with chloroplastic glycolipids, while the 22:6n-3 is predominant in the other cellular compartments [51–53].

Surprisingly, 20:5n-3 in PL was more enriched than 18:4n-3, its precursor in the n-3 pathway [4]. As AE of 20:5n-3 is higher than AE of 18:4n-3, it seems very unlikely that 20:5n-3 was produced via the pathway involving 18:4n-3 elongation and 20:4n-3 Δ5 desaturation. The existence of the alternative Δ8 desaturase pathway have been studied before in *Isochrysis galbana* and *Pavlova lutheri* [8,54,55]. However, as for the n-6 PUFA, intermediates (20:3n-3 and 20:4n-3) of the alternative Δ8 pathway were not detected by fatty acid analysis of *Isochrysis galbana* [8]. Similarly, in our study, intermediates (20:3n-3 and 20:4n-3) of this pathway to synthesize 20:5n-3 have not been found in sufficient amount to be measured by CSIA. As proposed by Qi et al. (2002) [8] for *Isochrysis galbana*, the synthesis of 20:5n-3 via 20:3n-3 and 20:4n-3 by *Tisochrysis lutea* might be very rapid, explaining why these two intermediates were only found in trace amounts (0.12% and 0.02% in PL and 0.05% and 0.33% in NL, respectively).
Due to their lower enrichments, 18:4n-3 and 18:5n-3 seemed unlikely involved in long chain PUFA synthesis such as 20:5n-3 and 22:6n-3. Based on the enrichment dynamics, elongation of 20:5n-3 into 22:5n-3 and further desaturation into 22:6n-3, respectively, by Δ5 elongase and Δ4 desaturase could be possible in *Tisochrysis lutea*. Ratio 22:5n-3/20:5n-3 in PL was within the threshold, indicating a simultaneous enrichment of both 20:5n-3 and 22:5n-3 in *T. lutea*. Such enzymes have been evidenced in haptophytes [54,56,57].

Considering the diversity of PKS gene in haptophytes [47], the possibility of produc-
tion of 22:6n-3 directly by PKS PUFA synthesis pathway might be possible, as previously
shown with thraustochytrids [10,58]. Synthesis of 22:6n-3 by PKS pathway might be at
play in parallel with the n-3 pathway. Indeed, we identified a protein cluster gathering
the four main domains potentially coding for the enzymes used in PKS PUFA synthesis
pathways: ketoacyl reductase (KR), polyketide synthase (KS), dehydrase/dehydrogenase
(DH), and enoyl reductase (ER). Protein clusters are groups of similar proteins that most
likely shared the same or similar functions [59]. By considering this cluster (candidate
proteins TISO_14962, TISO_14968, TISO_14968, TISO_14973, TISO_14975, and TISO_14977,
Figure 8), it could be possible that these proteins act together and allow n-3 PUFA synthesis
via PKS pathway. Interestingly, protein TISO_14973, while possessing only two of the four
domains of interest (KS and DH), presented a specific n-3 domain. This protein might
act concomitantly with the other proteins of the same cluster and allow the access to the
missing reductase activities (KR and ER). Finally, the ratio 22:5n-6/22:6n-3 was below the
threshold value making possible the conversion of 22:5n-6 into 22:6n-3 if we assumed that
ω3-desaturase might exist in haptophyte. Synthesis of 22:6n-3 by both n-3 and n-6 pathway
might be feasible in *Tisochrysis lutea* (Figure 10). These different ways to produce 22:6n-3
might contribute to betaine lipids synthesis. Indeed, betaine lipids are generally highly
unsaturated in C20 and C22 PUFA, especially in 22:6n-3 in haptophytes [60–63].

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

**Figure 10.** Hypothesized pathways to produce DHA in *T. lutea*. Numbers in boxes correspond to final mean AE value, and number in the yellow box is the mean value of ratio of the two surrounding fatty acids. The triangles symbolize the desaturases (front-end in yellow and methyl-end in purple), the circles the enzymes involved in PKS pathway (KR: 3-ketoacyl synthase, KS: 3-ketoacyl-ACP-reductase, DH: dehydrase, 2.2I/2.3I: 2-trans, 2-cis or 2-trans, 3-cis isomerases, ER: enoyl-reductase), and the squares the elongases.
# 4. Material and Methods

## 4.1. Algal Culture and ¹³C Labelling

This study was conducted following the experimental design described by Remize et al. (2020) [34]. The marine prymnesiophyte *Tisochrysis lutea* (T-iso, CCAP 927/14) was cultured in 2 L batch condition in balloons under continuous light (24 h light cycle, 100 µmoles photons m⁻²·s⁻¹) at 20 °C and with pH regulation at 7.50 ± 0.05 by CO₂ injection. Filtered seawater was previously enriched with Conway medium [64] and inoculated with *T. lutea* preculture in a growing stage (exponential phase, four days old). The experimental setup was composed of two cultures (Tl1 and Tl2) receiving the labelling ¹³C-CO₂ gas (Sigma-Aldrich, <3%atom ¹⁸O, 99.0%atom ¹³C) and one culture (TlT) receiving petrochemical CO₂ gas. ¹³C-incorporation in Tl1 and Tl2 began after inoculation (t₀) and was maintained for 24 h (t₂₄).

During the first hours of the experiment, the ¹³CO₂ injection tube of balloon Tl1 had been temporarily disconnected from the system. Consequently, balloon Tl2 received earlier the ¹³CO₂ and thus started to incorporate ¹³C before balloon Tl1.

## 4.2. Samples Collection

Sampling was performed as described in Remize et al. (2020) [34], i.e., at 30 min, 1 h, 2 h, 3 h, 4 h, and then every 2 h. A total of 16 samples was collected during the 24 h monitoring. At each sampling time, a total volume of 30 to 70 mL was collected for (i) flow cytometry analysis of cellular parameters, (ii) bulk isotopic analysis of particulate organic carbon (¹³C-POC) and dissolved inorganic carbon (¹³C-DIC) by EA-IRMS, (iii) fatty acid (FA) analysis in neutral lipids (NL) and polar lipids (PL) by GC-FID, and (iv) compound specific isotope analysis (CSIA) of FA (¹³C-FA) by GC-IRMS, as described in the following paragraphs.

## 4.3. Flow Cytometry Analysis

Algae growth cellular variables were measured using an Easy-Cyte Plus 6HT flow cytometer (Guava Merck Millipore®, Darmstadt, Germany) equipped with a 488 nm blue laser, detectors of forward (FSC) and side (SSC) light scatters, and three fluorescence detectors: green (525/30 nm), yellow (583/26 nm), and red (680/30 nm). The protocol, the variables studied, and the probes used for this cytometry following are described in Remize et al. (2020) [34]. Briefly, forward scatter (FSC), side scatter (SSC), and red fluorescence (FL3, red emission filter long pass, 670 nm) are used to study, respectively, cell size, complexity, and chlorophyll content. The fluorescent probe (SYTOX, Molecular Probes, Invitrogen, Eugene OR, USA, final concentration of 0.05 µM) was used to assess cell viability on FL1 detector (green fluorescence). The concentration of bacteria was also monitored by using SYBR® Green (Molecular Probes, Invitrogen, Eugene, OR, USA, #S7563) on FL1 detector. Concentrations of algae and bacteria were given cells per mL, and cellular variables were expressed in arbitrary units (a.u).

## 4.4. POC Concentration and Bulk Carbon Isotopic Composition

For particulate organic carbon (POC) and stable isotopic composition measurements, 30–70 mL of samples were filtered through pre-combusted 0.7 µm nominal pore-size glass fiber filters (Whatman GF/F, Maidstone, UK). The filter was processed, subsampled, and encapsulated as described in Remize et al. (2020) [34]. POC concentrations of all samples were measured using a CE Elantech NC2100 (ThermoScientific, Lakewood, NJ, USA) according to protocol by Remize et al. (2020) [34]. Bulk carbon isotopic composition (¹³C-POC) was analyzed by continuous flow on an Elemental Analyzer (EA, Flash 2000; Thermo Scientific, Bremen, Germany) coupled to a Delta V+ isotope ratio mass spectrometer (Thermo Scientific). Calibration was performed with international standards and in-house standard described in Table 2.
<table><caption>**Table 2.** List of international and in-house standards used for EA-IRMS and GB-IRMS analysis.</caption><thead><tr><th>Description</th><th>Nature</th><th>Analysis</th><th>δ<sup>13</sup>C (‰)</th><th>SD</th></tr></thead><tbody><tr><td>IAEA-CH<sub>6</sub></td><td>Sucrose (C<sub>12</sub>H<sub>22</sub>O<sub>11</sub>)</td><td>13C-POC</td><td>-10.45</td><td>0.03</td></tr><tr><td>IAEA-600</td><td>Caffeine (C<sub>8</sub>H<sub>10</sub>N<sub>4</sub>O<sub>2</sub>)</td><td>13C-POC</td><td>-27.77</td><td>0.04</td></tr><tr><td>Acetanilide</td><td>Acetanilide (C<sub>8</sub>H<sub>9</sub>NO)</td><td>13C-POC</td><td>+29.53</td><td>0.01</td></tr><tr><td>CA21 (in-house std)</td><td>Calcium carbonate (CaCO<sub>3</sub>)</td><td>13C-DIC</td><td>+1.476</td><td></td></tr><tr><td>Na<sub>2</sub>CO<sub>3</sub> (in-house std)</td><td>Sodium carbonate</td><td>13C-DIC</td><td>-6.8805</td><td></td></tr><tr><td>NaHCO<sub>3</sub> (in-house std)</td><td>Sodium bicarbonate</td><td>13C-DIC</td><td>-5.9325</td><td></td></tr></tbody></table>

## 4.5. DIC Concentration and Bulk Carbon Isotopic Composition

Samples for dissolved inorganic carbon (DIC) concentration and stable isotopic composition were collected from the filtrate of POC samples and processed as described in Remize et al. (2020) [34]. Analyses were conducted in a gas bench coupled to a Delta Plus mass spectrometer from Thermo Fisher Scientific, Bremen, Germany (GB-IRMS).

## 4.6. Isotopic Data Processing

We used the atomic proportion of ¹³C in percent (%atom of ¹³C) to express the results instead of the δ notation due to ¹³C-labelling. Conversion between δ notation and %atom¹³C notation can be done as follow [65]:

$$ \%\text{atom}^{13}\text{C} = 100 \times \frac{\left(\frac{\delta^{13}\text{C}}{1000} + 1\right) \times \left(\frac{^{13}\text{C}}{^{12}\text{C}}\right)_{\text{VPDB}}}{1 + \left(\frac{\delta^{13}\text{C}}{1000} + 1\right) \times \left(\frac{^{13}\text{C}}{^{12}\text{C}}\right)_{\text{VPDB}}} \quad (1) $$

where (¹³C/¹²C)PDB = 0.0112372, the ratio of ¹³C to ¹²C in the international reference VPDB standard.

Atomic enrichment (AE) of POC and DIC is then calculated from atom%¹³C-POC correction by POCcontrol values (i.e., corrected by 1.08%) and from atom%¹³C DIC corrected by control values (DICcontrol = 1.12%), respectively, according to the following equations:

$$ AE_{\text{POC}} = \%\text{atom}^{13}\text{C} - \text{POC}_{\text{control}} \quad (2) $$

$$ AE_{\text{DIC}} = \%\text{atom}^{13}\text{C} - \text{DIC}_{\text{control}} \quad (3) $$

## 4.7. Fatty Acids Analysis

### 4.7.1. Fatty Acids Analysis by Gas Chromatography Flame Ionisation Detector (GC-FID)

Lipid extraction, separation of neutral and polar lipid fractions, and transesterification processes are described elsewhere [34]. Fatty acids methyl esters (FAME) samples were analyzed by gas chromatography on a Varian CP8400 gas chromatograph (Agilent, Santa Clara, CA, USA) and separated concomitantly on two columns: one polar (ZB-WAX: 30 mm × 0.25 mm ID × 0.2 µm, Phenomenex, Torrance CA, USA) and the other apolar (ZB-5HT: 30 m × 0.25 mm ID × 0.2 µm, Phenomenex, Torrance CA, USA). The FAME of *T. lutea* were quantified using C23:0 as an internal standard (2.3 µg in each lipid fraction prior transmethylation) and were identified by comparison of their retention times with commercial standards (Supelco 37 component FAME mix, the PUFA No. 1 and No. 3 and the Bacterial Acid Methyl Esther Mix from Sigma-Aldrich, Darmstadt, Germany) and in-house standards mixtures. FA concentrations were reported as µg C·L⁻¹ and as % of total fatty acids from each lipid fraction. Thirty two fatty acids (FA) were thus identified and quantified: iso15:0, anteiso15:0, 14:0, 15:0, 16:0, 18:0, 22:0, 24:0, 14:1n-5, 16:1n-9, 16:1n-7, 17:1n-1, 18:1n-9, 18:1n-7, 16:2n-7, 16:2n-4, 16:4n-3, 18:2n-6, 18:3n-6, 18:3n-3, 18:4n-3, 18:5n-3, 20:2n-6, 20:3n-6, 20:4n-6, 20:4n-3, 20:5n-3, 22:2n-6, 22:4n-6, 22:5n-6, 22:5n-3, and 22:6n-3. Individual fatty acid and total fatty acid concentrations (as the sum of both fractions, named thereafter TFA) obtained in µg·L⁻¹ by GC-FID were also expressed in µmolC·L⁻¹
($\mu\text{g}\cdot\text{L}^{-1}$/molecular weight of individual fatty acid $\times$ carbon number of individual fatty acid) to ease the comparison with POC concentrations expressed in $\mu\text{molC}\cdot\text{L}^{-1}$ as well.

### 4.7.2. Fatty Acid Compound-Specific Isotope Analysis and Processing

Samples for compound-specific isotope analyses (CSIA) of FAME were performed on a Thermo Fisher Scientific GC ISOLINK TRACE ULTRA (Bremen, Germany) using the same apolar column as mentioned above for FAME analysis. Only the fatty acids with the highest concentrations, as measured by GC-FID analyses, were considered for CSIA (namely 14:0, 16:0, 18:0, 18:1n-9, 18:2n-6, 18:3n-3, 18:4n-3, 18:5n-3, 20:5n-3, 22:5n-6, 22:5n-3, and 22:6n-3). The other FA presenting a too low signal amplitude (<800 mV) on the GC-c-IRMS did not allow precise isotope ratio analysis. Additionally, 18:1n-9 and 18:3n-3 co-eluted for GC-c-IRMS on the apolar column, but most of the isotopic signature for neutral lipids (NL) is attributed to 18:1n-9. However, in the polar fraction (PL), 18:1n-9 and 18:3n-3 are in relatively similar proportion and so were considered together. Additionally, $^{13}$C enrichment of 18:5n-3 could only be measured in the polar lipid fraction, but its concentration was too low in neutral lipid fraction to measure its isotope composition.

To evidence FA conversion of fatty acid A into fatty acid B in *T. lutea*, we calculated the $AE_{FA}$ ratio (R) of product B over expected precursor A. R was defined with a confidence interval calculated at $\alpha = 0.1$ (defined arbitrarily) as follows:

$$R = \frac{AE_{FA(B)}}{AE_{FA(A)}} \qquad (4)$$

where A is the fatty acid hypothesized to be a precursor to fatty acid B, and $AE_{FA(A)}$ and $AE_{FA(B)}$ are their respective atomic enrichments at each sampling time.

If the AE of the product (B) exceeds the AE of the reactant (A), ratio $> 1$, then it is necessary to consider another formation process for B, since any molecule formed from A would have the same AE as A or below. If the ratio is $< 1$, transformation of A into B is considered possible. If the ratio is close to 1, the fatty acids A and B are at equilibrium in terms of label incorporated, implying B is then synthesized simultaneously or very rapidly from A.

### 4.8. Identification within the Silico Proteome of *T. lutea* of PKS Enzymes Involved PUFA Synthesis

The in silico proteome generated from last annotated version of the genome of *T. lutea* was used to identify putative proteins involved in n-3 PUFA PKS pathways [66]. We used the PKS previously identified in the haptophyte *Chrysochromulina tobin* as query for a BLASTp analysis, using e-value $< 10^{-3}$ as threshold [67]. Analysis of conserved domain was performed using the NCBI CD database V3,18 with e-value $< 10^{-2}$ as threshold. The genome location of genes encoding selected proteins was identified to evaluate genes' proximity and occurrence of gene clusters. Proteins and cluster of proteins containing the four domains ketoacyl reductase (KR), dehydrase/dehydrogenase (DH), enoyl reductase (ER), and polyketide synthase (KS) were selected as candidates.

### 4.9. Statistical Analysis

To assess the potential effect of time and difference between balloons during algae development and of $^{13}$CO$_2$ incorporation, Bartlett tests and ANOVA were performed on physiological and biochemical parameters, as well as PERMANOVA analysis on FA percentage separately in NL and PL. All statistical analyses were performed using R software.

## 5. Conclusions

The synthesis of long-chain PUFA in *Tisochrysis lutea* appeared to involve multiple pathways (Figures 8–10). First, the assumption of the use of PKS pathway for 22:5n-6 (DPA-6) was attested regarding the fast enrichment observed for this FA as well as the absence of detectable intermediates more or equally enriched. PKS pathway appeared to
be particularly efficient in *T. lutea* and induced a strong incorporation of the ¹³C-marker. However, the possibility of use of the conventional n-6 PUFA pathway should not be excluded, as 18:2n-6 presented a similar level of enrichment as 22:5n-6. It would only endorse that following desaturation and elongation steps to form 22:5n-6 were particularly dynamic and thus did not allow the accumulation of the ¹³C-label into n-6 intermediates. Within n-3 PUFA pathway, the Δ6-desaturase route seemed slower than the n-6 pathway in *T. lutea* in producing the two C₁₈ polyunsaturated fatty acids 18:4n-3 and 18:5n-3. We assumed 18:4n-3 and 18:5n-3 were unlikely synthesis intermediates of 20:5n-3 and 22:6n-3, as their enrichments were lower than the latter. Although 22:6n-3 was present in higher proportion than 22:5n-6, it was not enriched as fast, possibly because its synthesis may be more complex. Indeed, 22:6n-3 could be synthesized by *Tisochrysis lutea* via a combination multiples pathway: from 22:5n-6 via ω-3 desaturase pathway, from desaturation and elongation of 20:5n-3 and 22:5n-3, and via PKS pathway. Further studies are needed to better constrain the plausible routes taken by this prymnesiophyte to produce long chain PUFA.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/md20010022/s1, Table S1: Cellular parameters of *Tisochrysis lutea* (morphology (FSC and SSC), viability (FL1-SYTOX), and chlorophyll content (FL3) using flow cytometry analysis (mean ± SD of the 3 balloons). Values for SYTOX are in%, values for FL3/SSC/FSC in arbitrary unit (a.u), Table S2: Concentrations in µg·L⁻¹ of all identified and quantified FA in neutral and polar lipid fractions according to sampling time, Table S3: Concentrations in µmol·C·L⁻¹ of all identified and quantified FA in neutral and polar lipid fractions according to sampling time, Table S4: Concentrations in% of all identified and quantified FA in neutral and polar lipid fractions according to sampling time. Table S5: Mean ratio of atomic enrichment (AE) for pairs of FA (FAA vs. FAB) in the neutral lipids (NL) (mean ± SD, *n* = 9 sampling dates t8 to t24) for the two enriched balloons (Tl1, Tl2, Tl = *Tisochrysis lutea*). If the ratio is equal to or close to 1, A and B are assumed at equilibrium, and B is synthesized quickly from A; if the ratio is below 1, the transformation of B from A is possible but slow. Finally, if the ratio is above 1, A is not a main precursor of B, which has to be synthesized by a different pathway. Table S6: List of the potential candidate protein sequences involved in *Tisochrysis lutea* PKS synthesis pathway. The suspected function of each protein has been assumed using the NCBI conserved domain database (CDD) (Marchler-Bauer et al., 2017) by identifying the role of each domain recognized in the sequence. In columns KS/KR/DH/ER are written the number of domain corresponding to these functions in the studied sequences. ACS: acetyl-CoA synthetase, A_NRPS: adenylation domain of the non-ribosomal peptide synthetase (NRPS), Croto: crotonase/enoyl-CoA hydratase, EntF: enterobactin non-ribosomal peptide synthetase or thioesterase domain of Type I PKS, FAAL: fatty acyl-AMP ligase, GrsT: alpha/beta hydrolase, HM: hydroxymethylglutaryl-CoA synthase, MT: methyltransferase, PP: phosphopantetheine-binding (="swinging arm"), Sulf: sulfo-transferase, Thio: thioesterase. Figure S1: Proportions (%) of NL vs. PL (A) and proportions (%) of fifteen fatty acids in the NL fraction in average over the 24 h (B). Figure S2: Atomic enrichment of 11 main fatty acids in the polar lipid (NL) fraction during a 24h ¹³C labelling experiment. Tl: *Tisochrysis lutea*. Reference [68] is cited in the supplementary materials.

**Author Contributions:** Conceptualization, M.R., F.P. and P.S.; data curation, M.R.; Formal analysis, F.P.; funding acquisition, M.R., A.N.L., A.V. and P.S.; investigation, M.R., F.P., F.L.G., C.L., M.G. and P.S.; methodology, M.R., F.P., A.N.L., F.L.G., A.B., M.G. and P.S.; project administration, M.R. and P.S.; resources, R.C.; software, M.R., F.L.G. and A.B.; supervision, M.R., F.P., A.N.L. and P.S.; validation, M.R. and C.L.; visualization, M.R.; writing—original draft, M.R.; writing—review and editing, M.R., F.P., A.N.L., F.L.G., A.B., C.L., M.G, R.C., A.V. and P.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Université de Bretagne Occidentale (UBO) and the Center for Marine Sciences (CMS) at the University of North Carolina Wilmington (UNCW), the Interdisciplinary School for the Blue Planet (ISblue), and the Walter-Zellidja grant of the Académie Française.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.
**Acknowledgments:** We would like to thank Philippe Miner, Nelly Le Goic, Margaux Mathieu-Resuge, Elodie Fleury, Adeline Bidault, Korydwen Terrasson, and Corentin Baudet for their help during sampling and Oanez Lebeau for her support during the isotopic analyses.

**Conflicts of Interest:** The authors declare no conflict 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. FAO. *The State of World Fisheries and Aquaculture 2020*; FAO: Quebec, QC, Canada, 2020; ISBN 978-92-5-132692-3.

2. Arao, T.; Yamada, M. Biosynthesis of Polyunsaturated Fatty Acids in the Marine Diatom, Phaeodactylum Tricornutum. *Phytochemistry* **1994**, *35*, 1177–1181. [CrossRef]

3. Metz, J.G.; Roessler, P.; Facciotti, D.; Levering, C.; Dittrich, F.; Lassner, M.; Valentine, R.; Lardizabal, K.; Domergue, F.; Yamada, A.; et al. Production of Polyunsaturated Fatty Acids by Polyketide Synthases in Both Prokaryotes and Eukaryotes. *Science* **2001**, *293*, 290–293. [CrossRef] [PubMed]

4. Bell, M.V.; Tocher, D.R. Biosynthesis of polyunsaturated fatty acids in aquatic ecosystems: General pathways and new directions. In *Lipids in Aquatic Ecosystems*; Kainz, M., Brett, M.T., Arts, M.T., Eds.; Springer: New York, NY, USA, 2009; pp. 211–236. ISBN 978-0-387-88607-7.

5. Gurr, M.I.; Harwood, J.L.; Frayn, K.N. Fatty acids structure and metabolism: Fatty acids biosynthesis. In *Lipid Biochemistry: An Introduction*; Wiley-Blackwell: New York, NY, USA, 2002; pp. 21–59. ISBN 0-632-05409-3.

6. Guschina, I.A.; Harwood, J.L. Lipids and Lipid Metabolism in Eukaryotic Algae. *Prog. Lipid Res.* **2006**, *45*, 160–186. [CrossRef]

7. Harwood, J.L.; Guschina, I.A. The Versatility of Algae and Their Lipid Metabolism. *Biochimie* **2009**, *91*, 679–684. [CrossRef]

8. Qi, B.; Beaudoin, F.; Fraser, T.; Stobart, A.K.; Napier, J.A.; Lazarus, C.M. Identification of a cDNA encoding a novel C18-Δ9 polyunsaturated fatty acid-specific elongating activity from the docosahexaenoic acid (DHA)-producing microalga, Isochrysis galbana. *FEBS Lett.* **2002**, *510*, 159–165. [CrossRef]

9. Hauvermale, A.; Kuner, J.; Rosenzweig, B.; Guerra, D.; Diltz, S.; Metz, J.G. Fatty Acid Production in Schizochytrium Sp.: Involvement of a Polyunsaturated Fatty Acid Synthase and a Type I Fatty Acid Synthase. *Lipids* **2006**, *41*, 739–747. [CrossRef] [PubMed]

10. Ye, C.; Qiao, W.; Yu, X.; Ji, X.; Huang, H.; Collier, J.L.; Liu, L. Reconstruction and Analysis of the Genome-Scale Metabolic Model of Schizochytrium limacinum SR21 for Docosahexaenoic Acid Production. *BMC Genom.* **2015**, *16*, 799. [CrossRef]

11. Ratledge, C. Fatty Acid Biosynthesis in Microorganisms Being Used for Single Cell Oil Production. *Biochimie* **2004**, *86*, 807–815. [CrossRef] [PubMed]

12. John, U.; Beszteri, B.; Derelle, E.; Van de Peer, Y.; Read, B.; Moreau, H.; Cembella, A. Novel Insights into Evolution of Protistan Polyketide Synthases through Phylogenomic Analysis. *Protist* **2008**, *159*, 21–30. [CrossRef]

13. Freitag, M.; Beszteri, S.; Vogel, H.; John, U. Effects of Physiological Shock Treatments on Toxicity and Polyketide Synthase Gene Expression in Prymnesium parvum (Prymnesiophyceae). *Eur. J. Phycol.* **2011**, *46*, 193–201. [CrossRef]

14. Eichholz, K.; Beszteri, B.; John, U. Putative Monofunctional Type I Polyketide Synthase Units: A Dinoflagellate-Specific Feature? *PLoS ONE* **2012**, *7*, e48624. [CrossRef] [PubMed]

15. Armenta, R.E.; Valentine, M.C. Single-Cell Oils as a Source of Omega-3 Fatty Acids: An Overview of Recent Advances. *J. Am. Oil Chem. Soc.* **2013**, *90*, 167–182. [CrossRef]

16. Wallis, J.G.; Watts, J.L.; Browse, J. Polyunsaturated Fatty Acid Synthesis: What Will They Think of Next? *Trends Biochem. Sci.* **2002**, *27*, 467–473. [CrossRef]

17. Uttaro, A. Biosynthesis of Polyunsaturated Fatty Acids in Lower Eukaryotes. *IUBMB Life* **2006**, *58*, 563–571. [CrossRef] [PubMed]

18. Caramujo, M.-J.; Boschker, H.T.S.; Admiraal, W. Fatty Acid Profiles of Algae Mark the Development and Composition of Harpacticoid Copepods. *Freshw. Biol.* **2007**, *53*, 77–90. [CrossRef]

19. Pereira, H.; Barreira, L.; Figueiredo, F.; Custódio, L.; Vizetto-Duarte, C.; Polo, C.; Rešek, E.; Engelen, A.; Varela, J. Polyunsaturated Fatty Acids of Marine Macroalgae: Potential for Nutritional and Pharmaceutical Applications. *Mar. Drugs* **2012**, *10*, 1920–1935. [CrossRef]

20. de Carvalho, C.; Caramujo, M.-J. Fatty Acids as a Tool to Understand Microbial Diversity and Their Role in Food Webs of Mediterranean Temporary Ponds. *Molecules* **2014**, *19*, 5570–5598. [CrossRef] [PubMed]

21. Hixson, S.M.; Arts, M.T. Climate Warming Is Predicted to Reduce Omega-3, Long-Chain, Polyunsaturated Fatty Acid Production in Phytoplankton. *Glob. Chang. Biol.* **2016**, *22*, 2744–2755. [CrossRef] [PubMed]

22. Courchesne, N.M.D.; Parisien, A.; Wang, B.; Lan, C.Q. Enhancement of Lipid Production Using Biochemical, Genetic and Transcription Factor Engineering Approaches. *J. Biotechnol.* **2009**, *141*, 31–41. [CrossRef] [PubMed]

23. Ahmann, K.; Heilmann, M.; Feussner, I. Identification of a Δ4-Desaturase from the Microalga Ostreococcus lucimarinus. *Eur. J. Lipid Sci. Technol.* **2011**, *113*, 832–840. [CrossRef]

24. Vaezi, R.; Napier, J.; Sayanova, O. Identification and Functional Characterization of Genes Encoding Omega-3 Polyunsaturated Fatty Acid Biosynthetic Activities from Unicellular Microalgae. *Mar. Drugs* **2013**, *11*, 5116–5129. [CrossRef]
25. Kimura, K.; Okuda, S.; Nakayama, K.; Shikata, T.; Takahashi, F.; Yamaguchi, H.; Skamoto, S.; Yamaguchi, M.; Tomaru, Y. RNA Sequencing Revealed Numerous Polyketide Synthase Genes in the Harmful Dinoflagellate Karenia mikimotoi. *PLoS ONE* **2015**, *10*, e0142731. [CrossRef] [PubMed]

26. de Swaaf, M.E.; de Rijk, T.C.; van der Meer, P.; Eggink, G.; Sijtsma, L. Analysis of Docosahexaenoic Acid Biosynthesis in Cryptthecodinium cohnii by 13C Labelling and Desaturase Inhibitor Experiments. *J. Biotechnol.* **2003**, *103*, 21–29. [CrossRef]

27. Van den Meersche, K.; Middelburg, J.J.; Soetaert, K.; van Rijswijk, P.; Boschker, H.T.S.; Heip, C.H.R. Carbon-Nitrogen Coupling and Algal-Bacterial Interactions during an Experimental Bloom: Modeling a ¹³C Tracer Experiment. *Limnol. Oceanogr.* **2004**, *49*, 862–878. [CrossRef]

28. Ecker, J.; Scherer, M.; Schmitz, G.; Liebisch, G. A Rapid GC-MS Method for Quantification of Positional and Geometric Isomers of Fatty Acid Methyl Esters. *J. Chromatogr. B* **2012**, *897*, 98–104. [CrossRef] [PubMed]

29. Zhang, H.; Wu, C.; Wu, Q.; Dai, J.; Song, Y. Metabolic Flux Analysis of Lipid Biosynthesis in the Yeast Yarrowia lipolytica Using 13C-Labeled Glucose and Gas Chromatography-Mass Spectrometry. *PLoS ONE* **2016**, *11*, e0159187. [CrossRef] [PubMed]

30. Zhao, J.; Shimizu, K. Metabolic Flux Analysis of Escherichia coli K12 Grown on 13C-Labeled Acetate and Glucose Using GC-MS and Powerful Flux Calculation Method. *J. Biotechnol.* **2003**, *101*, 101–117. [CrossRef]

31. Xiong, W.; Liu, L.; Wu, C.; Yang, C.; Wu, Q. 13C-Tracer and Gas Chromatography-Mass Spectrometry Analyses Reveal Metabolic Flux Distribution in the Oleaginous Microalga Chlorella protothecoides. *Plant Physiol.* **2010**, *154*, 1001–1011. [CrossRef] [PubMed]

32. Martzolff, A.; Cahoreau, E.; Cogne, G.; Peyriga, L.; Portais, J.-C.; Dechandol, E.; Le Grand, F.; Massou, S.; Gonçalves, O.; Pruvost, J.; et al. Photobioreactor Design for Isotopic Non-Stationary 13C-Metabolic Flux Analysis (INST 13C-MFA) under Photoautotrophic Conditions. *Biotechnol. Bioeng.* **2012**, *109*, 3030–3040. [CrossRef] [PubMed]

33. Cui, J.; Diao, J.; Sun, T.; Shi, M.; Liu, L.; Wang, F.; Chen, L.; Zhang, W. 13C Metabolic Flux Analysis of Enhanced Lipid Accumulation Modulated by Ethanolamine in Cryptthecodinium cohnii. *Front. Microbiol.* **2018**, *9*, 956. [CrossRef] [PubMed]

34. Remize, M.; Planchon, F.; Loh, A.N.; Le Grand, F.; Bideau, A.; Le Goic, N.; Fleury, E.; Miner, P.; Corvaisier, R.; Volety, A.; et al. Study of Synthesis Pathways of the Essential Polyunsaturated Fatty Acid 20:5n-3 in the Diatom Chaetoceros muelleri Using 13C-Isotope Labeling. *Biomolecules* **2020**, *10*, 797. [CrossRef] [PubMed]

35. Remize, M.; Planchon, F.; Loh, A.N.; Le Grand, F.; Lambert, C.; Bideau, A.; Bidault, A.; Corvaisier, R.; Volety, A.; Soudant, P. Identification of Polyunsaturated Fatty Acids Synthesis Pathways in the Toxic Dinophyte Alexandrium minutum Using 13C-Labelling. *Biomolecules* **2020**, *10*, 1428. [CrossRef]

36. Grosse, J.; van Breugel, P.; Boschker, H.T.S. Tracing Carbon Fixation in Phytoplankton-Compound Specific and Total ¹³C Incorporation Rates: ¹³C Uptake into Macromolecules. *Limnol. Oceanogr. Methods* **2015**, *13*, 288–302. [CrossRef]

37. Menzel, R.; Ngosong, C.; Ruess, L. Isotopologue Profiling Enables Insights into Dietary Routing and Metabolism of Trophic Biomarker Fatty Acids. *Chemoecology* **2017**, *27*, 101–114. [CrossRef]

38. Wei, X.; Shi, B.; Koo, I.; Yin, X.; Lorkiewicz, P.; Suhail, H.; Rattan, R.; Giri, S.; McClain, C.J.; Zhang, X. Analysis of Stable Isotope Assisted Metabolomics Data Acquired by GC-MS. *Anal. Chim. Acta* **2017**, *980*, 25–32. [CrossRef]

39. Lima, V.F.; de Souza, L.P.; Williams, T.C.R.; Fernie, A.R.; Daloso, D.M. Gas Chromatography–Mass Spectrometry-Based 13C-Labeling Studies in Plant Metabolomics. In *Plant Metabolomics*, António, C., Ed.; Methods in Molecular Biology; Springer New York: New York, NY, USA, 2018; Volume 1778, pp. 47–58. ISBN 978-1-4939-7818-2.

40. Eikrem, W.; Medlin, L.K.; Henderiks, J.; Rokitta, S.; Rost, B.; Probert, I.; Throndsen, J.; Edvardsen, B. *Haptophyta. Handbook of the Protists*; Archibald, J.M., Simpson, A.G.B., Slamovits, C.H., Margulis, L., Melkonian, M., Chapman, D.J., Corliss, J.O., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 1–61. ISBN 978-3-319-32669-6.

41. 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]

42. Pernet, F.; Tremblay, R.; Demers, E.; Roussy, M. Variation of Lipid Class and Fatty Acid Composition of *Chaetoceros muelleri* and *Isochrysis* sp. Grown in a Semicontinuous System. *Aquaculture* **2003**, *221*, 393–406. [CrossRef]

43. Patil, V.; Källqvist, T.; Olsen, E.; Vogt, G.; Gislerød, H.R. Fatty Acid Composition of 12 Microalgae for Possible Use in Aquaculture Feed. *Aquacult Int.* **2007**, *15*, 1–9. [CrossRef]

44. Gladyshev, M.I.; Sushchik, N.N.; Makhutova, O.N. Production of EPA and DHA in Aquatic Ecosystems and Their Transfer to the Land. *Prostaglandins Other Lipid Mediat.* **2013**, *107*, 117–126. [CrossRef]

45. da Costa, F.; Le Grand, F.; Quéré, C.; Bougaran, G.; Cadoret, J.P.; Robert, R.; Soudant, P. Effects of Growth Phase and Nitrogen Limitation on Biochemical Composition of Two Strains of Tisochrysis lutea. *Algal. Res.* **2017**, *27*, 177–189. [CrossRef]

46. 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]

47. Kohli, G.S.; John, U.; Van Dolah, F.M.; Murray, S.A. Evolutionary Distinctiveness of Fatty Acid and Polyketide Synthesis in Eukaryotes. *ISME J.* **2016**, *10*, 1877–1890. [CrossRef] [PubMed]

48. Wang, S.; Zheng, L.; Cui, Z.; Chen, J.; Yang, B.; Han, X.; Liu, C. Cloning and Molecular Characterization of a Delta-6 Fatty Acid Desaturase Gene from *Isochrysis* sp. CCMM5001. *J. Appl. Phycol.* **2016**, *28*, 921–929. [CrossRef]

49. Joseph, J.D. Identification of 3, 6, 9, 12, 15-Octadecapentaenoic Acid in Laboratory-Cultured Photosynthetic Dinoflagellates. *Lipids* **1975**, *10*, 395–403. [CrossRef] [PubMed]
50. Kotajima, T.; Shiraiwa, Y.; Suzuki, I. Functional Screening of a Novel Δ15 Fatty Acid Desaturase from the Coccolithophorid Emiliania huxleyi. *Biochim. Biophys. Acta (BBA)-Mol. Cell Biol. Lipids* **2014**, *1841*, 1451–1458. [CrossRef] [PubMed]

51. Heydarizadeh, P.; Poirier, I.; Loizeau, D.; Ulmann, L.; Mimouni, V.; Schoefs, B.; Bertrand, M. Plastids of Marine Phytoplankton Produce Bioactive Pigments and Lipids. *Mar. Drugs* **2013**, *11*, 3425–3471. [CrossRef]

52. Nalder, T.D.; Miller, M.R.; Packer, M.A. Changes in Lipid Class Content and Composition of Isochrysis Sp. (T-Iso) Grown in Batch Culture. *Aquacult Int.* **2015**, *23*, 1293–1312. [CrossRef]

53. Zhong, Y.; Li, Y.; Xu, J.; Cao, J.; Zhou, C.; Yan, X. Isolation of Chloroplasts from Marine Microalgae Isochrysis Galbana Parke Suitable for Organelle Lipid Composition Analysis; In Review; Creative Commons: Mountain View, CA, USA, 2020.

54. Zhou, X.-R.; Robert, S.S.; Petrie, J.R.; Frampton, D.M.F.; Mansour, M.P.; Blackburn, S.I.; Nichols, P.D.; Green, A.G.; Singh, S.P. Isolation and Characterization of Genes from the Marine Microalga Pavlova salina Encoding Three Front-End Desaturases Involved in Docosahexaenoic Acid Biosynthesis. *Phytochemistry* **2007**, *68*, 785–796. [CrossRef]

55. Fraser, T.C.M.; Qi, B.; Elhussein, S.; Chatrattanakunchai, S.; Stobart, A.K.; Lazarus, C.M. Expression of the Isochrysis C18-Δ9 Polyunsaturated Fatty Acid Specific Elongase Component Alters Arabidopsis Glycerolipid Profiles. *Plant Physiol.* **2004**, *135*, 859–866. [CrossRef]

56. Tonon, T.; Harvey, D.; Larson, T.R.; Graham, I.A. Identification of a Very Long Chain Polyunsaturated Fatty Acid Δ4-Desaturase from the Microalga Pavlova lutheri. *FEBS Lett.* **2003**, *553*, 440–444. [CrossRef]

57. Robert, S.S.; Petrie, J.R.; Zhou, X.-R.; Mansour, M.P.; Blackburn, S.I.; Green, A.G.; Singh, S.P.; Nichols, P.D. Isolation and Characterisation of a Δ5-Fatty Acid Elongase from the Marine Microalga Pavlova salina. *Mar. Biotechnol.* **2009**, *11*, 410–418. [CrossRef]

58. Zhang, K.; Li, H.; Chen, W.; Zhao, M.; Cui, H.; Min, Q.; Wang, H.; Chen, S.; Li, D. Regulation of the Docosapentaenoic Acid/Docosahexaenoic Acid Ratio (DPA/DHA Ratio) in Schizochytrium limacinum B4D1. *Appl. Biochem. Biotechnol.* **2017**, *182*, 67–81. [CrossRef]

59. Zaslavsky, L.; Ciufo, S.; Fedorov, B.; Tatusova, T. Clustering Analysis of Proteins from Microbial Genomes at Multiple Levels of Resolution. *BMC Bioinform.* **2016**, *17*, 276. [CrossRef]

60. Kato, M.; Sakai, M.; Adachi, K.; Ikemoto, H.; Sano, H. Distribution of Betaine Lipids in Marine Algae. *Phytochemistry* **1996**, *42*, 1341–1345. [CrossRef]

61. Eichenberger, W.; Gribi, C. Lipids of Pavlova lutheri: Cellular Site and Metabolic Role of DGCC. *Phytochemistry* **1997**, *45*, 1561–1567. [CrossRef]

62. Meireles, L.A.; Guedes, A.C.; Malcata, F.X. Lipid Class Composition of the Microalga Pavlova lutheri: Eicosapentaenoic and Docosahexaenoic Acids. *J. Agric. Food Chem.* **2003**, *51*, 2237–2241. [CrossRef] [PubMed]

63. Armada, I.; Hachero-Cruzado, I.; Mazuelos, N.; Ríos, J.L.; Manchado, M.; Cañavate, J.P. Differences in Betaine Lipids and Fatty Acids between Pseudoisochrysis paradoxa VLP and Diacronema Vlkianum VLP Isolates (Haptophyta). *Phytochemistry* **2013**, *95*, 224–233. [CrossRef] [PubMed]

64. Walne, P.R. Experiments on the Large-Scale Culture of the Larvae Ostrea edulis. *J. Fish. Invest. Min. Agric. Fish. Lond.* **1966**, *2*, 25–53.

65. Hunter, W.R.; Jamieson, A.; Huvenne, V.A.I.; Witte, U. Sediment Community Responses to Marine vs. Terrigenous Organic Matter in a Submarine Canyon. *Biogeosciences* **2013**, *10*, 67–80. [CrossRef]

66. Berthelier, J.; Casse, N.; Daccord, N.; Jamilloux, V.; Saint-Jean, B.; Carrier, G. A Transposable Element Annotation Pipeline and Expression Analysis Reveal Potentially Active Elements in the Microalga Tisochrysis lutea. *BMC Genom.* **2018**, *19*, 378. [CrossRef]

67. John, U.; Beszteri, S.; Glöckner, G.; Singh, R.; Medlin, L.; Cembella, A.D. Genomic Characterisation of the Ichthyotoxic Prymnesiophyte Chrysochromulina polylepis, and the Expression of Polyketide Synthase Genes in Synchronized Cultures. *Eur. J. Phycol.* **2010**, *45*, 215–229. [CrossRef]

68. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C.J.; Lu, S.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; et al. CDD/SPARCLE: Functional classification of proteins via subfamily domain architectures. *Nucleic Acids Res.* **2017**, *45*, D200–D203. [CrossRef]