# H3K27me3 and EZH Are Involved in the Control of the Heat-Stress-Elicited Morphological Changes in Diatoms

Mhammad Zarif, Ellyn Rousselot, Bruno Jesus, Leïla Tirichine, Céline Duc

► To cite this version:

Mhammad Zarif, Ellyn Rousselot, Bruno Jesus, Leïla Tirichine, Céline Duc. H3K27me3 and EZH Are Involved in the Control of the Heat-Stress-Elicited Morphological Changes in Diatoms. International Journal of Molecular Sciences, 2024, 25 (15), pp.8373. 10.3390/ijms25158373. hal-04678502

HAL Id: hal-04678502

https://hal.science/hal-04678502v1

Submitted on 27 Aug 2024

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

H3K27me3 and EZH Are Involved
in the Control of the Heat-Stress-
Elicited Morphological Changes in
Diatoms

Mhammad Zarif, Ellyn Rousselot, Bruno Jesus, Leïla Tirichine and Céline Duc

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

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

Article

# H3K27me3 and EZH Are Involved in the Control of the Heat-Stress-Elicited Morphological Changes in Diatoms

Mhammad Zarif <fnref n="1" />, Ellyn Rousselot <fnref n="1" />, Bruno Jesus <fnref n="2" /> ID, Leïla Tirichine <fnref n="1" />,<fnref n="3" /> ID and Céline Duc <fnref n="1" />,* ID

**Abstract:** Marine water temperatures are increasing due to anthropogenic climate change, constituting a major threat to marine ecosystems. Diatoms are major marine primary producers, and as such, they are subjected to marine heat waves and rising ocean temperatures. Additionally, under low tide, diatoms are regularly exposed to high temperatures. However, physiological and epigenetic responses to long-term exposure to heat stress remain largely unknown in the diatom *Phaeodactylum tricornutum*. In this study, we investigated changes in cell morphology, photosynthesis, and H3K27me3 abundance (an epigenetic mark consisting of the tri-methylation of lysine 27 on histone H3) after moderate and elevated heat stresses. Mutants impaired in PtEZH—the enzyme depositing H3K27me3—presented reduced growth and moderate changes in their PSII quantum capacities. We observed shape changes for the three morphotypes of *P. tricornutum* (fusiform, oval, and triradiate) in response to heat stress. These changes were found to be under the control of PtEZH. Additionally, both moderate and elevated heat stresses modulated the expression of genes encoding proteins involved in photosynthesis. Finally, heat stress elicited a reduction of genome-wide H3K27me3 levels in the various morphotypes. Hence, we provided direct evidence of epigenetic control of the H3K27me3 mark in the responses of *Phaeodactylum tricornutum* to heat stress.

**Keywords:** *Phaeodactylum tricornutum*; H3K27me3; EZH; heat stress; photosynthesis; morphotype

## 1. Introduction

Diatoms are an important group of unicellular phytoplankton, primarily photosynthetic. They contribute to ~20% of the total primary production [1]. Among phytoplankton, the coastal species *Phaeodactylum tricornutum* is distinguished from other diatoms by its frustule, which is either devoid of silica or poorly silicified. The silica deficiency provides this species with cell morphological plasticity. In addition to its morphological flexibility, *P. tricornutum* possesses an extensive array of molecular tools, making it an ideal model to investigate cellular morphological plasticity. This diatom has three major morphotypes: fusiform (FM), oval (OM), and triradiate (TM) [2]. FM is the most abundant and stable morphotype. OM cells are more resistant to stresses such as low salinity or low light, while TM is favored in unstressed conditions [3]. Moreover, several conditions, such as temperature, are known to impact morphotype conversion. However, the impact of high temperatures on marine diatoms as well as the molecular players controlling morphotype conversion, remain mostly unknown. *P. tricornutum* is able to adapt to various environments, and its morphogenesis appears to be independent of its genotype [4,5].

*P. tricornutum* is a marine photosynthetic organism using solar photons to produce organic matter, the light energy being converted into chemical energy through electron transfer reactions. The photosystem II, hereafter referred to as PSII, is a multi-protein

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

check for
updates

**Citation:** Zarif, M.; Rousselot, E.; Jesus, B.; Tirichine, L.; Duc, C. H3K27me3 and EZH Are Involved in the Control of the Heat-Stress-Elicited Morphological Changes in Diatoms. *Int. J. Mol. Sci.* **2024**, *25*, 8373. https://doi.org/10.3390/ijms25158373

Academic Editor: Fucheng Lin

Received: 28 June 2024

Revised: 23 July 2024

Accepted: 27 July 2024

Published: 31 July 2024

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

**Copyright:** © 2024 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/).

<footnote n="1">¹ Nantes Université, CNRS, US2B, UMR 6286, F-44000 Nantes, France; mhammad.zarif@univ-nantes.fr (M.Z.); leila.tirichine@univ-nantes.fr (L.T.)</footnote>

<footnote n="2">² Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France; bruno.jesus@univ-nantes.fr</footnote>

<footnote n="3">³ Institute for Marine and Antarctic Studies (IMAS), Ecology and Biodiversity Centre, University of Tasmania, Hobart, TAS 7004, Australia</footnote>

<footnote>* Correspondence: celine.duc@univ-nantes.fr</footnote>
complex that absorbs light energy. In PSII, there are two antenna proteins, PsbB and PsbC, that transfer energy to the PsbA and PsbD proteins of the reaction center. These four proteins are large transmembrane subunits. At its lumenal side, PSII also contains an oxygene-evolving complex (OEC) made of five proteins (PsbO, PsbQ, PsbU, PsbV, and Psb31) [6]. The OEC performs the water oxidation reaction. Finally, PSII contains small transmembrane proteins such as PsbM [6]. Electron transfer from PSII to PSI is mediated by the plastoquinone and then by the cytochrome *b6-f* complex. This complex is a plastoquinone-plastocyanin reductase. It is composed of several subunits, such as the cytochromes *f* (petA) and *b6* (petB), as well as the petC protein, which is an iron-sulfur cluster protein known as Rieske protein [7]. The plastocyanin electron carrier is absent in *P. tricornutum* and replaced by the cytochrome *c6* [8], which transfers electrons to photosystem I (PSI). In diatoms, PSI is composed of seven trans-membrane subunits (PsaA, PsaB, PsaF, PsaI, PsaJ, PsaL, and PsaM) and three stromal subunits (PsaC, PsaD, and PsaE) [9]. Finally, electrons are transferred from PSI to ferredoxin, which is attached to the PSI through the PsaA-PsaC-PsaE pocket [9]. Diatoms are able to dissipate excess light energy as heat through the binding of xanthophyll pigments to the light-harvesting complex (LHC). Under low light, zeaxanthin, via its antheraxanthin intermediate, is converted to violaxanthin by the zeaxanthin epoxidase (ZEP) enzymes (putatively ZEP1 to 3 in *P. tricornutum*). On the contrary, under excessive light, the reverse reaction is catalyzed by the violaxanthin de-epoxidase (VDE) enzymes (putatively VDE, VDL1, and VDL2 in *P. tricornutum*). A second xanthophyll cycle is constituted by the reversible conversion of diatoxanthin to diadinoxanthin, using the same enzymes [10]. These xanthophyll cycles protect diatoms from high light-induced damage caused to the photosynthetic apparatus. Diatoms are, thus, able to grow under a wide range of light intensities.

*P. tricornutum* is found in coastal waters, including tidal areas where temperature can be highly variable due to changing water levels. Besides, temperature is known to impact *P. tricornutum* morphology [3]. This diatom is atypical in that it is the only species known to have three major morphotypes. Anthropogenic climate change is notably triggering the warming and acidification of marine waters. These changing conditions impact the growth and physiology of diatoms [11]. The polycomb repressive complex 2 (PRC2) is involved in cell differentiation in various eukaryotes [12]. The PRC2 core proteins are the enhancer of zeste homologs (EZH), suppressor of zeste 12 (SUZ12), embryonic ectoderm development (EED), and retinoblastoma binding proteins 4/7 protein (RBBP4/7, also referred to as p48/p55) [13]. This complex deposits tri-methylation on the lysine 27 of the histone H3 (H3K27me3), an epigenetic mark considered as transcriptionally repressive. In *P. tricornutum*, the homolog of the catalytic subunit of PRC2 is the enhancer of zeste homolog protein (PtEZH). It controls cell differentiation in this species [14]. Previous studies have investigated the role of the H3K27me3 mark in thermo-morphogenesis in *Arabidopsis thaliana*. Indeed, intensive studies reported an epigenetic repression of *FLOWERING LOCUS C* (FLC), notably via PRC2 and H3K27me3 [15], and a genome-wide accumulation of H3K27me3 was reported in high temperatures [16]. However, the link between the H3K27me3, its methylase, and the heat-stress morphological responses has never been investigated in *P. tricornutum*.

To date, we know that cell differentiation in *P. tricornutum* is (i) independent of the genotype [2,4], (ii) controlled by PtEZH and H3K27me3 in standard conditions [14], and (iii) modulated by temperature [3]. Hence, we hypothesized that increasing water temperature due to global warming, especially in coastal regions and under low tide, will impact *P. tricornutum* morphology through epigenetic control, such as a modulation of H3K27me3 levels. To test our hypothesis, we used two different heat stress conditions (30 °C and 37 °C) and evaluated diatom responses at the physiological, morphological, molecular, and epigenetic levels. Besides, we used mutants impaired in PtEZH to investigate if and how H3K27me3 controls heat stress response in *P. tricornutum*. Overall, most *Ptezh* mutants have a reduced maximum specific growth rate in optimal growth conditions and moderate changes in their photo-physiological responses. We showed that the three morphotypes of
*P. tricornutum* (fusiform, oval, and triradiate) were affected in their cell shape during heat stress (after being grown under standard temperature conditions), a feature lost in *Ptezh* mutants. In response to heat stress conditions (30 °C and 37 °C), several genes encoding proteins involved in photosynthesis were differentially expressed in wild type and *Ptezh* mutants, with different expression patterns. Moreover, there was a strong reduction of genome-wide levels of the H3K27me3 epigenetic mark in response to heat stress.

## 2. Results

### 2.1. Analysis of the Physiological Responses in the Ptezh Mutants under Standard Conditions

*Ptezh* mutant growth was reported to be slightly slower compared to WT [14]. Therefore, we first characterized how *Ptezh* mutants are affected by the absence of H3K27me3 in standard temperature conditions before analyzing the heat stress responses. To do so, we monitored cell growth and photosynthetic features of the *Ptezh* mutants in the FM and TM genetic backgrounds. We used two mutants in each morphotype. FM-*del6*, TM-*M2.10*, and TM-*M2.11* are knock-out mutants [14]. These three mutants produce a truncated version of the PtEZH protein (Figure S1a,b). The fourth mutant, FM-*M2.3*, is a knock-down mutant [14] presenting an insertion of a phenylalanine residue in the EZH1/2-like domain (Figure S1c). We monitored cell growth at 19 °C in *Ptezh* mutants for 25 days (Figure S1e,f) and calculated the maximum specific growth rate (µmax). FM-*del6*, FM-*M2.3*, and TM-*M2.11* mutants presented a reduced µmax (Figure 1a,b), while TM-*M2.10* displayed a higher one (Figure 1b). However, TM-*M2.10* and TM-*M2.11* reached a stationary phase at a lower cell concentration compared to TM (Figure S1f). Hence, growth is perturbed under normal conditions by the absence of a functional PtEZH protein.

The chloroplast has a central role in the growth of photosynthetic organisms. Since most *Ptezh* mutants exhibited a decreased growth rate, we examined whether the loss of PtEZH could affect photosynthesis in *P. tricornutum*. To do so, we measured several parameters: the α light-limited initial slope for rapid light curves (RLC), the maximal relative electron transport rate (rETRmax), Ek (light saturation coefficient) that corresponds to the photon irradiance at the onset of light saturation coefficient, and maximal non-photochemical quenching (NPQmax) in *Ptezh* mutants. The α values were lower in FM-*M2-3* (Figure 1c) and higher in TM-*M2.11* and TM-*M2.10* (Figure 1d) compared to WT. Only the FM-*del6* mutant displayed a lower photosynthetic performance based on the rETRmax parameter compared to FM (Figure 1e,f). Moreover, the Ek coefficient, which is the light level at which ETR reaches a maximum and the relation between PAR irradiances and ETR stops to be linear, remained unchanged in *Ptezh* mutants compared to WT (Figure S1g,h). Furthermore, NPQmax (Figure 1g,h) and E50 (Figure 1i,j), which is the light level at which 50% of the NPQ occurs, were both reduced in FM-*M2.3* (Figure 1g-i) and TM-*M2.10* (Figure 1h-j) while the other mutants remained unaffected. This suggests that FM-*M2.3* and TM-*M2.10* mutants are less tolerant to higher light intensity.

NPQ enables the dissipation of excessive energy as heat and xanthophylls are important players in this photoprotective mechanism [17]. Since some *Ptezh* mutants had a reduced NPQmax, we investigated transcript levels of the genes encoding enzymes involved in the xanthophyll cycles: (i) the zeaxanthin epoxidase enzymes (ZEP1 to 3) and the violaxanthin de-epoxidase enzymes (VDE, VDL1 and VDL2). Moderate changes in transcript abundance were observed in the *Ptezh* mutants (Figure 2a,b). However, ZEP1 and ZEP2 transcript abundance was increased in FM-*del6* and FM-*M2.3*, respectively (Figure 2a). Therefore, the production of enzymes involved in the xanthophyll synthesis pathway did not seem to be altered by the loss of PtEZH.
<figure><img src="image_8.png" /></figure>

**Figure 1.** Growth and photosynthesis analysis of *Ptezh* mutants in FM and TM lines. (a,b) Growth rates are represented with error bars indicating standard deviations based on triplicate cultures in FM (a) and TM (b) lines. Maximum growth rates were determined from the growth curve obtained by cell counting performed every 2 days for 25 days on cultures grown at 19 °C (12/12 light–dark period; light intensity of 55 µmol photons.m⁻².s⁻¹). (c–j) Photosynthesis analysis performed when Fv/Fm was at its maximum for the *Ptezh* mutants in FM and TM lines grown at 19 °C (12/12 light–dark period; light intensity of 55 µmol photons.m⁻².s⁻¹). (c,d) α slope for rapid light curves in FM (c) and TM (d) lines. (e,f) rETRmax is the PSII maximum relative electron transfer rate in FM (e) and TM (f) lines. (g,h) NPQmax is the maximum non-photochemical quenching in FM (g) and TM (h) lines. (i,j) E50 is the light level at which 50% of the NPQ occurs in FM (i) and TM (j) lines. Student's *t*-test: *p* < 0.05; ns, non-significant.
<figure><img src="image_9.png" /><figcaption>**Figure 2.** Expression of genes encoding zeaxanthin epoxidase (ZEP), violaxanthin de-epoxidase (VDE), and violaxanthin de-epoxidase-like (VDL). (a,b) Relative transcript levels of genes encoding ZEP1, ZEP2, ZEP3, VDE, VDL1, and VDL2 enzymes in the FM (a) and TM (b) lines. The ZEP1, ZEP2, and ZEP3 enzymes putatively convert (i) zeaxanthin via the antheraxanthin intermediate to violaxanthin and (ii) diadinoxanthin in diatoxanthin. The VDE, VDL1, and VDL2 enzymes putatively catalyze the reverse reactions [10]. Levels are measured by qRT-PCR on three biological replicates consisting of 5-day-old cultures at 19 °C (12/12 light–dark period; light intensity of 55µmol photons.m⁻².s⁻¹) in (a,b). Student's t-test: * *p* < 0.05; ns, non-significant.</figcaption></figure>

## 2.2. Gene Expression Analysis of Photosynthesis Proteins in Ptezh Mutants

The measurements described above relied on the fluorescence mainly emitted by PSII. To complete our analysis of *Ptezh* physiology and, more specifically, of photosynthesis in standard conditions, we measured the transcript abundance of various genes encoding proteins involved in photosynthesis in wild type and *Ptezh* mutants. Regarding the PSII complex, transcript abundance of *PsbC*, *PsbA*, and *PsbD* increased in FM-*del6* compared to FM (Figure 3a). Transcript abundance of all PSII tested genes (with the exception of *PsbD*) increased in both TM-*M2.10* and TM-*M2.11* mutants (Figure 3b). Moreover, *petC1* and *Ferredoxin* transcript abundance increased in the four *Ptezh* mutants (Figure 3c). However, a *petC2* increased transcript abundance was monitored in both *Ptezh* TM mutants (Figure 3c). One should note that two genes were identified to encode Rieske proteins, *petC1* (Phatr3_J46657) and *petC2* (Phatr3_J13358) in *P. tricornutum*. Finally, for PSI, transcript abundance of *PsaA* and *PsaB* was increased in both *Ptezh* FM mutants, while *TM-M2.10* showed an increased *PsaC* transcript level (Figure 3d). Therefore, genes encoding proteins involved in photosynthesis were differentially expressed in *Ptezh* mutants, with different levels of expression in the FM and TM lines (Figure S3a).
<figure><img src="image_10.png" /><figcaption>**Figure 3.** Expression of genes encoding proteins of the photosynthesis apparatus. (a,b) Relative transcript levels of genes encoding proteins of PSII (antenna: PsbB and PsbC; reaction center: PsbA and PsbD; OEC: PsbO and PsbU; small transmembrane protein: PsbM) in the FM (a) and TM (b) lines. (c) Relative transcript levels of genes encoding proteins of the cytochrome b6-f complex (Rieske proteins: petC1 and petC2) and ferredoxin in the FM (upper panel) and TM (lower panel) lines. (d) Relative transcript levels of genes encoding proteins of PSI (trans-membrane subunits: PsaA, PsaB; stromal subunit: PsaC) in the FM (upper panel) and TM (lower panel) lines. Levels are displayed with a logarithmic scale (a–c) and measured by qRT-PCR on three biological replicates consisting of 5-day-old cultures grown at 19 °C (12/12 light–dark period; light intensity of 55 µmol photons.m⁻².s⁻¹) in (a–d). Student's t-test: * *p* < 0.05; ns, non-significant.</figcaption></figure>

## 2.3. Morphotype Conversion in Response to Prolonged Heat Stress

To decipher the relationship between temperature and the deposition of H3K27me3 by PtEZH, we applied two types of prolonged heat stress: (i) “moderate heat stress” (MHS) applied for 5 days at 30 °C on 7-day-old cultures pre-acclimated at 19 °C; (ii) “elevated heat stress” (EHS) for 2 days at 37 °C on 5-day-old cultures pre-acclimated at 19 °C.
### 2.3.1. Morphotype Conversion in Ecotypes of Pure Morphotypes

Temperature impact on morphotypes has previously been assessed in different ecotypes of *P. tricornutum* [3]. However, cultures were analyzed in response to a transition from cold (15 °C) to warm temperature (28 °C) in Pt1 8.6, Pt3, and Pt8 ecotypes. These lines presented a mixture of morphotypes, making it difficult to distinguish between morphotype conversion elicited by heat stress and a differential lethality between morphotypes within the same ecotype. To address this, we purified Pt1 8.6, Pt3, and Pt16 ecotypes to exclusively contain FM (Pt1 8.6-FM), OM (Pt3-OM), or TM (Pt16-TM) cells, respectively, at 19 °C to have representative lines for each morphotype. We used Pt9, a tropical ecotype presenting a mixture of FM and OM cells at 19 °C, as a control.

In both stress conditions, we observed the appearance of OM cells in Pt1 8.6-FM and FM cells in Pt3-OM (MHS, Figure 4a; EHS, Figure 4c). However, FM cells appeared in Pt16-TM only in response to EHS (Figure 4c). The tropical line Pt9 did not show any changes in morphotype abundance during the first two days in MHS (Figure 4a) and in EHS (Figure 4c). As a control, we assessed the relative abundance of morphotypes for each line maintained at 19 °C and did not observe any changes over time (Figure S2a,c), indicating that morphotype conversion was due to heat stress. We then compared cell growth for each morphotype in order to understand if changes in morphotype abundance could be due to changes in cell viability in response to heat stress. In Pt1 8.6-FM, FM cells display a reduced growth after 2 days at 30 °C (Figure 4b), while OM cells have an increased growth at 30 °C (Figure 4b) and 37 °C (Figure 4d). In Pt3-OM, OM cells displayed a reduced growth after 3 days at 30 °C while FM cells had a constant growth (Figure 4b). In Pt16-TM, TM cells displayed a reduced growth after 2 days at 37 °C (Figure 4d). In Pt9, a similar growth pattern was observed for FM and OM cells in MHS (Figure 4b) and EHS (Figure 4c). Therefore, the changes observed in the relative abundance of morphotypes in response to heat stress were the consequence of decreased growth of the pre-existing morphotype simultaneously with increased growth of the appearing morphotype.

<figure><img src="image_11.png" /><figcaption>**Figure 4.** *Cont.*</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 4.** The impact of prolonged-heat stress exposures in representative ecotypes of *P. tricornutum*. (a,c) Relative abundance of each morphotype in Pt1 8.6-FM, Pt3-OM, and Pt16-TM that are representative of the FM, OM, and TM morphotypes, respectively. The Pt9 tropical strain was used as a control. (b,d) Cell concentration measured for each morphotype (FM and OM) for Pt1 8.6-FM, Pt3-OM and Pt9 subjected to MHS (b) and for Pt1 8.6-FM, Pt3-OM, Pt16-TM, and Pt9 subjected to EHS (d). Relative abundance presented in panels a and c were inferred from these measures. A moderate heat stress (MHS) was applied for 5 days at 30 °C to 7-day-old cultures pre-acclimated at 19 °C (a,b). An elevated heat stress (EHS) was applied for 2 days at 37 °C to 5-day-old cultures pre-acclimated at 19 °C (c,d). Measures were carried out before transferring the cultures to heat-stress conditions (T0). They were then performed every day for 5 days in MHS and 2 days in EHS after transfer under heat stress.</figcaption></figure>

Furthermore, morphotype conversion can occur with or without cell division. To determine whether morphotype conversion observed in response to prolonged heat stress aligns with cell division, we synchronized cells through a prolonged darkness period of 36 h [18]. We then monitored morphotype changes in response to EHS since we noticed changes for the three morphotypes only at 37 °C. We observed that OM cells appeared in Pt1 8.6-FM after 4 h at 37 °C and the abundance is stable after 6 h (Figure 5). FM cells appeared in Pt3-OM after 4 h at 37 °C and then the relative abundance of FM cells were increasing until +8 h (Figure 5). Regarding Pt16-TM, TM cells appeared only after 10 h at 37 °C (Figure 5). *P. tricornutum* divides approximately once per day and cell synchronization decreased after 10 h following illumination [18]. Therefore, morphotype conversion in Pt1 8.6-FM seemed to occur without cell division while, in Pt3, it resulted from OM cells converting to FM cells without cell division and then with cell division. Finally, morphotype conversion in Pt16-TM seemed to occur mainly through cell division.
<figure><img src="image_13.png" /><figcaption>**Figure 5.** Kinetics of the morphological changes observed during exposure of 5-day-old cultures to 37 °C in representative ecotypes of *P. tricornutum*. Synchronized cells were transferred to 37 °C at T0, simultaneously with illumination to release cell cycle progression from G1. The relative abundance of each morphotype in Pt1 8.6-FM (**left**), Pt3-OM (**middle**), and Pt16-TM (**right**) was then assessed every 2 h for 14 h.</figcaption></figure>

### 2.3.2. Regulation of Morphotype Conversion by PtEZH during a Prolonged Heat Stress

We then analyzed how H3K27me3 depletion would impact the heat stress response in *P. tricornutum*. To do so, we monitored the morphotype abundance in various *Ptezh* mutants under both MHS and EHS. *Ptezh* mutants were impaired in PtEZH function and depleted in H3K27me3 [14]. We observed the appearance of OM cells in Pt1 8.6-FM but not in the *FM-del6* and *FM-M2.3* mutants in both MHS (Figure 6a) and EHS (Figure 6c). The *TM-M2.10* and *TM-M2.11* were generated in the Pt8 ecotype [14] presenting a mixture of TM and FM cells [2]. Therefore, we purified Pt8 to have only TM cells (Pt8-TM), as previously carried out for Pt16-TM, the other triradiate ecotype used above. FM cells appeared in the Pt8-TM line in response to EHS (Figure 6d) but not to MHS (Figure 6b). The *TM-M2.10* and *TM-M2.11* mutants have fusiform cells, and no morphotype changes were observed in either MHS (Figure 6b) or EHS (Figure 6d). Therefore, the loss of PtEZH and the associated H3K27me3 depletion led to the loss of temperature-dependent morphotype conversion in all the investigated *Ptezh* mutants.

<figure><img src="image_14.png" /><figcaption>**Figure 6.** *Cont.*</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 6.** Impact of prolonged-heat stress exposures on the morphology of *Ptezh* mutants. (a,b) Relative abundance of FM (a) and TM (b) lines subjected to a moderate heat stress (MHS). (c,d) Relative abundance of FM (c) and TM (d) lines subjected to an elevated heat stress (EHS). A moderate heat stress (MHS) was applied for 5 days at 30 °C on 7-day-old cultures pre-acclimated at 19 °C (a,b). An elevated heat stress (EHS) was applied for 2 days at 37 °C on 5-day-old cultures pre-acclimated at 19 °C (c,d). Measures were carried out before transferring cultures to heat stress conditions (T0). They were then performed every day for 5 days in MHS and 2 days in EHS after transfer under heat stress.</figcaption></figure>

## 2.4. Analysis of the Effects of Prolonged Heat Stresses on Photosynthesis Gene Expression

Besides the heat-triggered morphological conversion, we also investigated the physiological response of *Ptezh* mutants in response to heat stress. To do so, we analyzed the transcript abundance of genes encoding (i) enzymes of the xanthophyll cycle, (ii) proteins involved in photosynthesis, and (iii) heat shock proteins.

Firstly, we investigated the transcript abundance of the genes encoding enzymes involved in the xanthophyll cycles. In response to MHS, the *ZEP1*, *ZEP2*, *ZEP3*, *VDE*, and *VDL1* transcript abundance were reduced in FM and *FM-del6* (Figure 7a). In response to EHS, *FM-M2.3* exhibited increased transcript abundances for all genes (with the exception of *ZEP3*, Figure 7b). On the contrary, FM presented reduced transcript abundances for all genes in EHS (with the exception of *VDL1* and *VDL2*, Figure 7b). The *ZEP2* and *ZEP3* transcript abundance was reduced in *FM-del6* while *VDE* transcript abundance increased in EHS (Figure 7b). Therefore, a prolonged exposure to 30 °C down-regulates the expression of genes encoding enzymes involved in the xanthophyll cycle in the wild type and the *FM-del6* mutant. Surprisingly, a prolonged exposure to 37 °C up-regulates the expression of most genes encoding enzymes involved in the xanthophyll cycle in the knock-down mutant *FM-M2.3*.
<figure><img src="image_16.png" /><figcaption>**Figure 7.** Expression of genes encoding zeaxanthin epoxidase (ZEP), violaxanthin de-epoxidase (VDE), and violaxanthin de-epoxidase-like (VDL) in response to prolonged heat stresses. (a,b) Relative transcript level of genes encoding ZEP1, ZEP2, ZEP3, VDE, VDL1, and VDL2 enzymes in MHS (a) and EHS (b) lines. The ZEP1, ZEP2, and ZEP3 enzymes putatively convert (i) zeaxanthin *via* the antheraxanthin intermediate to violaxanthin and (ii) diadinoxanthin in diatoxanthin. The VDE, VDL1, and VDL2 enzymes putatively catalyze the reverse reactions [10]. Levels are displayed with a logarithmic scale and measured by qRT-PCR on three biological replicates consisting of 7-day-old cultures pre-acclimated at 19 °C transferred to 30 °C for 4 days (a) and 5-day-old cultures pre-acclimated at 19 °C transferred to 37 °C for 2 days (b). Student's *t*-test; \* *p* < 0.05.</figcaption></figure>

Secondly, we measured transcript abundance for a set of photosynthetic genes in FM *Ptezh* mutants subjected to either MHS or EHS. In response to MHS, the *PsbB*, *PsbC*, and *PsbD* transcript abundance increased in FM-*del6* and FM-*M2.3* while the *PsbA*, *PsbO*, *PsbU*, and *PsbM* transcript abundance decreased in wild type and *Ptezh* mutants (Figure 8a). In response to EHS, transcript abundance of most genes related to PSII also decreased (Figure 9a). The *petC2* transcript level decreased in wild type and *Ptezh* mutants under both MHS (Figure 8b) and EHS conditions (Figure 9b). The *Ferredoxin* transcript abundance increased in MHS in wild type and *Ptezh* mutants (Figure 8b) and only in wild type in EHS (Figure 9b). For the PSI complex, *PsaA* and *PsaB* transcript abundances increased in wild type in MHS (Figure 8c) and EHS (Figure 9c) while *PsaC* abundance decreased in both conditions (Figures 8c and 9c). The *Ptezh* mutants exhibited changes in *PsaA*, *PsaB*, and *PsaC* transcript abundance only in response to EHS in (Figure 9c). This suggests that heat stress impacts gene expression of proteins involved in photosynthesis to different extents in wild type and *Ptezh* mutants, depending on the applied temperature.
<figure><img src="image_17.png" /></figure>

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

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

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

**Figure 8.** Expression of genes encoding proteins of the photosynthesis apparatus in response to MHS. (a) Relative transcript levels of genes encoding proteins of PSII (antenna: PsbB and PsbC; reaction center: PsbA and PsbD; OEC: PsbO and PsbU; small transmembrane protein: PsbM). (b) Relative transcript levels of genes encoding proteins of the cytochrome b6-f complex (Rieske proteins: petC1 and petC2) and ferredoxin. (c) Relative transcript levels of genes encoding proteins of PSI (trans-membrane subunits: PsaA, PsaB; stromal subunit: PsaC). (d) Relative transcript levels of genes encoding heat shock proteins (HSP, heat shock protein; HSF3, heat shock factor protein 3). Levels are displayed with a logarithmic scale and measured by qRT-PCR on three biological replicates consisting of 7-day-old cultures pre-acclimated at 19 °C transferred to 30 °C for 4 days. Student's t-test; * *p* < 0.05.
<figure><img src="image_21.png" /><figcaption>**Figure 9.** Expression of genes encoding proteins of the photosynthesis apparatus in response to EHS. (a) Relative transcript level of genes encoding proteins of PSII (antenna: PsbB and PsbC; reaction center: PsbA and PsbD; OEC: PsbO and PsbU; small transmembrane protein: PsbM). (b) Relative transcript level of genes encoding proteins of the cytochrome b6-f complex (Rieske proteins: petC1 and petC2) and ferredoxin. (c) Relative transcript level of genes encoding proteins of PSI (trans-membrane subunits: PsaA, PsaB; stromal subunit: PsaC). (d) Relative transcript level of genes encoding heat shock proteins (HSP, heat shock protein; HSF3, heat shock factor protein 3). Levels are displayed with a logarithmic scale and measured by qRT-PCR on three biological replicates consisting of 5-day-old cultures pre-acclimated at 19 °C transferred to 37 °C for 2 days. Student's *t*-test; \* *p* < 0.05.</figcaption></figure>

Thirdly, since heat-shock protein production is modulated by various stresses such as heat, we measured the transcript abundance of a set of genes encoding heat-shock proteins (HSP70G, HSF3, HPSP20, and HPSP20A). Transcript abundances decreased in both wild type and *Ptezh* mutants for the four tested genes in response to MHS (Figure 8d) and EHS (Figure 9d). Therefore, prolonged exposure to increased temperatures greatly decreased the expression of heat-shock proteins.

## 2.5. *Analysis of the Effects of Prolonged Heat Stress on H3K27me3 Levels*

We showed above that the heat-triggered morphological changes are lost in the absence of a functional PtEZH protein that deposits the H3K27me3 epigenetic mark. Therefore, we
investigated whether prolonged heat stress modulates the abundance of H3K27me3 on a genome-wide scale. To do so, we extracted proteins as described by [19] and determined the amount of the H3K27me3 mark by Western blotting under MHS. This experiment was carried out in the Pt1 8.6-FM, Pt3-OM, and Pt16-TM ecotypes, which harbor only FM, OM, or TM cells, respectively, as well as in the Pt9 tropical ecotype. In response to MHS, we observed a complete loss or a strong reduction of H3K27me3 levels in Pt1 8.6-FM, Pt3-OM, and Pt16-TM ecotypes, which represent the FM, OM, and TM morphotypes, respectively (Figure 10). The tropical ecotype Pt9 for which, morphotype abundance remained unchanged in response to MHS (Figure 4a) also exhibited a loss of H3K27me3 (Figure 10). Therefore, exposure to prolonged heat stress leads to a strong reduction of the H3K27me3 epigenetic mark at a genome-wide scale.

<figure><img src="image_22.png" /><figcaption>**Figure 10.** Genome-wide levels of the epigenetic mark H3K27me3 in response to MHS. The H3K27me3 levels were quantified by Western blotting in 7-day-old cultures pre-acclimated at 19 °C transferred to 30 °C for 4 days. H4 was used as a loading control. Presented blots are representative of several blots and biological replicates. Two biological replicates (referred to as 1 and 2) consisting of independent cultures collected at the same time are presented for each ecotype.</figcaption></figure>

**3. Discussion**

In this study, we investigated how some epigenetic players control morphological changes and differentiation in *P. tricornutum*. We hypothesized that PtEZH, the catalytic subunit of the PRC2 complex that deposits the H3K27me3 epigenetic mark, is involved in the heat-elicited morphotype changes. We first carried out a characterization of *Ptezh* mutants to investigate how depletion of H3K27me3 could impact the diatom physiology besides the previously described impact on morphology and differentiation [14]. Second, we applied two prolonged heat stresses (30 °C and 37 °C) to reveal how temperature increase could impact morphological changes and the expression of photosynthesis-related genes in *P. tricornutum* wild type and *Ptezh* mutants. Finally, we investigated how PtEZH could control heat-triggered morphotype changes in *P. tricornutum* and how H3K27me3 levels are modulated by a prolonged heat stress.

In a previous study [14], *Ptezh* mutants showed slower growth monitored over a period of 17 days. Moreover, only *FM-del6* and *TM-M2.11* mutants were analyzed. Here, we monitored the growth for 25 days of the two mutants available in the FM (*FM-del6* and *FM-M2.3*) and TM (*TM-M2.10* and *TM-M2.11*) lines. The maximum specific growth rate was reduced in all mutants, with the exception of *TM-M2.10*. Since *TM-M2.10* showed higher α and lower NPQ values at lower light levels, this mutant might be more efficient in converting light energy under the low light intensities used in our study, which might explain why it grew quicker (Figure 1b) than the other lines. Furthermore, with a longer monitoring time, we observed that only *FM-M2.3*, which is a knock-down mutant, displayed a higher cell concentration at 25 days compared to its wild type (Figure S1e). Besides, only *FM-del6* showed a lower photosynthetic performance based on the rETRmax parameter compared to FM. Both *FM-M2.3* and *TM-M2.10* mutants seemed to have an increased sensitivity to excess light energy based on the reduced NPQmax and E50 parameters. Therefore, *Ptezh* mutants show moderate changes in their photosynthetic capacities in standard conditions.

*P. tricornutum* has an optimal growth at 19 °C and a suboptimal growth above 25 °C. This diatom is considered to be eurythermal. To investigate heat stress responses, we chose three purified ecotypes, Pt1 8.6-FM, Pt3-OM, and Pt16-TM, representing the fusiform, oval, and triradiate morphotypes, respectively. The growth rate of OM cells in Pt3-OM
was reported to be decreased at 28 °C [3]. We confirmed this finding and found that the conversion of OM cells to FM cells occurs without cell division in the first 4h at 37 °C then in the course of cell division (Figure 5). Besides, the conversion of FM to OM cells in Pt1 8.6-FM at 37 °C occurs without cell division while in Pt16-TM, TM cells convert to FM cells with cell division. In this study, we chose to apply two different temperatures (30 °C–MHS and 37 °C–EHS). The Pt8 ecotype (presenting a mixture of TM and FM cells as described in [3]) did not show any changes in morphotype abundance during a transition from 15 °C to 28 °C [3]. Therefore, we applied a higher temperature to see if we could trigger a heat-elicited morphotype change in a triradiate ecotype. Indeed, both pure triradiate lines used in our study (Pt16-TM and Pt8-TM) showed conversion of TM to FM at 37 °C but not at 30 °C (summarized in Figure S3c). It is interesting to notice that the Pt1 8.6-FM, Pt3-OM, and Pt16-TM ecotypes originally contain a mixture of FM/OM, OM/FM, and TM/FM. Thus, heat stress carried out on pure ecotypes (Pt1 8.6-FM, Pt3-OM, and Pt16-TM) triggered the re-appearance of the secondary morphotype. Moreover, this secondary morphotype shows a higher growth rate under heat stress, suggesting its potential role in providing a better adaptation to changing environmental conditions. Given the predominant shift to the fusiform morphotype in response to higher temperatures and its general prevalence, it is plausible that the fusiform morphotype has evolved in response to the rising temperatures experienced over the last century. This morphotype might show an optimized surface-to-volume ratio that confers advantages in challenging environmental conditions and is likely to increase nutrient uptake [20]. Moreover, elongated cells dominate during the summer period [21] when water temperatures are higher, suggesting that higher temperatures favor the fusiform morphotype.

In land plants, photosynthesis is reported to be a highly sensitive physiological process. Indeed, photosynthesis inhibition occurs under moderate heat stress, while the photosynthesis apparatus is damaged under elevated heat stress. More specifically, the photosystem II and its OEC are the most sensitive components [22]. Expression of some genes encoding proteins of petC, PSII, and the light-harvesting complex of the PSI were shown to be down-regulated in warming-adapted populations of *P. tricornutum* [23]. Here, we reported a reduction of *PsbO* and *PsbU* transcript abundance in response to moderate and elevated heat stresses in wild type and *Ptezh* mutants (Figure S3b). This likely leads to a reduction in PsbO and PsbU protein levels. These two proteins are cofactors of PSII, and their loss triggers the inactivation of the plant PSII during heat stress [22]. Thus, this strongly suggests a PSII impairment in response to heat stress in *P. tricornutum*. Besides, transcript abundance for genes encoding PSI proteins increased in response to moderate and elevated heat stresses in wild type and *Ptezh* mutants (Figure S3b). Hence, *P. tricornutum* might present a higher PSI activity under heat stress. We also observed that transcript abundance for genes encoding cytochrome *b6-f* subunits is affected by prolonged heat stress (Figure S3b). Therefore, electron and energy flow from PSII to PSI might follow a different distribution under heat stress compared to optimal growth conditions in the diatom *P. tricornutum*. This is likely associated with a different balance between xanthophyll pigments, based on the changes we observed in the expression of genes encoding enzymes involved in the xanthophyll cycle. Furthermore, the heat-shock proteins participate in cell protection during heat stress and contribute to heat tolerance, notably in plants. It has been recently reported in *P. tricornutum* that HSP70A might be involved in the folding of PsbA and PsbD, the reaction center proteins of PSII [24]. Additionally, *HSP70A* expression was strongly induced after 1 h at 26 °C, although it returned to levels similar to optimal conditions after 5 h [24]. In our study, we observed that prolonged heat stress triggered a decrease in the expression of genes encoding the heat-shock proteins HSP70G, HSF3, HPSP20, and HPSP20A (Figure S3b). Therefore, the PSII impairment may result from reduced protection by these molecular chaperones.

Epigenetic control of heat stress responses has been reported in various eukaryotes (reviewed in [25]). In *P. tricornutum*, specific sets of long non-coding RNAs were identified in each morphotype (FM, OM, or TM), suggesting a role of these epigenetic players in
morphotype identity [26]. Furthermore, global DNA methylation levels in CHG and CHH contexts were reported to decrease within gene bodies under warming conditions [23]. Importantly, the triradiate cell shape is abolished in *Ptezh* mutants TM-M2.10 and TM-M2.11, which lack the epigenetic mark H3K27me3 [14]. Our study also reported the role of epigenetics in the control of morphogenesis. Using a reverse genetic approach, we showed that the heat-elicited morphotype changes are lost in the *FM-del6*, *FM-M2.3*, *TM-M2.10*, and *TM-M2.11* mutants (summarized in Figure S3c). This led us to conclude that PtEZH controls morphotype changes observed during prolonged heat stresses. It will be interesting to complete this analysis with *Ptezh* mutants in the OM morphotype when available. We also observed a strong decrease or a loss of the H3K27me3 epigenetic mark when 7-day-old cultures of *P. tricornutum*, pre-acclimated at 19 °C, were subjected to prolonged heat stress at 30 °C for 4 days. Thus, *Ptezh* mutants that are already depleted in the H3K27me3 epigenetic mark cannot modulate the levels of this epigenetic mark, leading to an absence of heat-triggered morphotype changes. Therefore, we could postulate that in response to heat stress, PtEZH controls morphological changes via the H3K27me3 epigenetic mark. Further studies would be pertinent in analyzing morphological changes in response to heat stress in mutants impaired in the H3K27me3 demethylase, an enzyme that has not yet been identified [27]. Furthermore, the H3K27me3 levels are more elevated in TM compared to FM [14]. Therefore, a higher temperature might be necessary to elicit the morphotype conversion and a complete H3K27me3 depletion, explaining why morphological changes in triradiate lines Pt8-TM and Pt16-TM were observed at 37 °C, not at 30 °C. Finally, we observed a genome-wide depletion of the H3K27me3 epigenetic mark in response to a moderate heat stress in *P. tricornutum*. However, in *Arabidopsis thaliana*, an opposite observation was reported. Growth under warm temperature (27 °C) led to a genome-wide accumulation of H3K27me3 in this plant [16]. The epigenetic control of heat stress responses might thus be different in marine and terrestrial photosynthetic organisms. Hence, this study opened new avenues regarding the investigation of heat stress responses that are under the control of PtEZH and H3K27me3 in *P. tricornutum*. As further perspectives, it will be relevant to investigate the genome-wide distribution of H3K27me3 and PtEZH in each morphotype (FM, OM, or TM) upon heat stress. In parallel, a transcriptomic analysis will be pertinent to identify key regulators involved in the control of heat-elicited morphotype changes, based on changes in their expression concomitantly with their H3K27me3 levels. Besides, analyzing the cross talk between various epigenetic marks such as the transcription permissive and repressive ones (namely H3K4me3 and H3K36me3 vs. H3K9me3 and H3K27me3) as well as DNA methylation will be valuable to understand how epigenetics triggers heat-elicited morphotype changes in the model diatom *P. tricornutum*.

To conclude, this study reports the first evidence of heat-triggered morphological changes in *P. tricornutum* controlled by PtEZH—the diatom enhancer of zeste homolog that is the catalytic subunit of PRC2—and the epigenetic mark H3K27me3. Each morphotype of *P. tricornutum* adapts differently to a given environment. OM cells, for instance, have better acclimation in benthic environments such as sediment surfaces, showing strong adherence and survival abilities to limiting growth conditions. Conversely, FM and TM cells are better adapted to unstressed growth conditions and pelagic environments [28]. Marine temperature warming will, therefore, modulate morphotype abundance in the diatom *P. tricornutum* and likely in other primary marine producers. This might result in an imbalanced ecosystem impacting community structure and function.

# 4. Materials and Methods

## 4.1. Used Ecotypes and Genetically-Modified Lines

Several ecotypes of *P. tricornutum* were used in the current study. For the following ecotypes, pure morphotypes were isolated as follows: FM for clone Pt1 8.6 (CCAP 1055/1) originated from an estuary near Blackpool (UK), OM for Pt3 (CCAP 1052/1B, CCMP2558), an ecotype originally isolated from the ecotype Pt2 (Plymouth, UK), TM for Pt16 (RCC641) originated from Helgoland (Helgoland, Atlantic Ocean, North Sea, Germany). To isolate
pure morphotypes, cultures presenting a mixture of morphotypes were streaked on EASW-Agar (1/2 EASW [29], 1% agar) plates. After 21 days of culture at 19 °C, under 12/12 light-dark period with a light intensity of 55 µmol photons.m⁻².s⁻¹, isolated colonies were resuspended in EASW. Serial dilutions were performed in 24-well plates to isolate dilutions containing only cells of one single morphotype. All lines were grown axenically using enhanced artificial sea water (EASW) [29] in batch cultures at 19 °C, under 12/12 light-dark period with a light intensity of 55 µmol photons.m⁻².s⁻¹. For the heat stress experiments, we used the tropical ecotype Pt9 (CCMP633) collected in the Territory of Guam (Northern Mariana, Islands, Micronesia), which presents a mixture of fusiform and oval cells at 19 °C [2]. In order to analyze the effects of PtEZH depletion, the following published lines [14] were used: for FM, the Pt1 8.6 ecotype without (referred to as FM-WT) and with the Cas9 control vector introduced by bombardment (referred to as FM-Cas9) and the *Ptezh* knock-out lines FM-Del6 and FM-M2.3; for TM, the Pt8 (CCMP2560) ecotype (described in [2]) collected in the Jericho Beach (Canada) without (referred to as TM-WT) and with the Cas9 control vector introduced by bombardment (referred to as TM-Cas9) and the *Ptezh* knock-out lines TM-M2.10 and TM-M2.11. To assess growth of wild type and *Ptezh* mutants in FM and TM lines, cells were grown in 30 mL of EASW with an initial concentration of 10⁵ cells/mL. Cell counts were measured on triplicate cultures using flow cytometry (CytoFLEX, Beckman Coulter Life Sciences, Brea, CA, USA). Counting was performed every 2 days for 25 days with 1 mL of sample taken from each line. Maximum growth rates (µmax, d⁻¹) were determined by fitting growth kinetic data with a Gompertz model [30].

## 4.2. *Heat Stress Experiments*

The effect of moderated (MHS) and elevated (EHS) heat stress on morphotype conversion and viability was assessed in Pt1 8.6-FM (pure FM), Pt3-OM (pure OM), Pt16-TM (pure TM) and Pt9 (FM + OM) lines. For MHS, cultures were transferred to 30 °C (under 12/12 light-dark period with a light intensity of 55µmol photons.m⁻².s⁻¹) after a 7-day pre-acclimation at 19 °C (under 12/12 light-dark period with a light intensity of 55 µmol photons.m⁻².s⁻¹) of cultures initially inoculated at a concentration of 10⁵ cells.mL⁻¹. Then, the cell concentration and abundance of each morphotype were analyzed daily for 5 days with a Malassez chamber for cell counting. Briefly, cells were diluted, and counting was performed on ten chambers, the Malassez cell containing 100 chambers for a total volume of 1µL. For EHS, cultures were transferred to 37 °C (under 12/12 light-dark period with a light intensity of 55 µmol photons.m⁻².s⁻¹) after a 5-day pre-acclimation at 19 °C (under 12/12 light-dark period with a light intensity of 55 µmol photons.m⁻².s⁻¹) of cultures initially inoculated at a concentration of 10⁵ cells.mL⁻¹. Then, the cell concentration and abundance of each morphotype were analyzed daily for 2 days with a Malassez chamber for cell counting. For each ecotype, three replicates were monitored and MHS and EHS experiments were repeated three and two times, respectively.

## 4.3. *Heat Stress Experiments on Synchronised Cells*

*P. tricornutum* cells were arrested in the G1 phase by prolonged darkness (36 h). Cells were released from this G1 checkpoint by illumination, and they were simultaneously transferred to 37 °C with a light intensity of 55 µmol photons.m⁻².s⁻¹ for heat stress. Cell counting and morphotype abundance were assessed every 2 h during 14 h, starting at illumination (T0). For each ecotype, six replicates were monitored.

## 4.4. *PAM Fluorescence Measurements*

For PAM fluorescence measurements, cells were grown at 19 °C (under 12/12 light-dark period with a light intensity of 55 µmol photons.m⁻².s⁻¹) in 20 mL of EASW with an initial concentration of 10⁵ cells/mL. Fluorescence analyses were performed every day and data analyses (presented in Figure 1c-j and Figure S1g,h) were performed when Fv/Fm was at its maximum. The Fv/Fm parameter represents the maximum PSII quantum efficiency and was calculated as follows: $Fv/Fm = (Fm - Fo)/Fm$, with Fm being the maximum
imposed with the Chimera software v1.16 [36]. Protein domains were determined from InterProScan (5.69-101.0) [37].

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25158373/s1.

**Author Contributions:** Conceptualization, C.D., B.J. and L.T.; experiments and data acquisition, C.D., E.R. and M.Z.; formal analysis, C.D., B.J. and L.T.; data interpretation, C.D., B.J. and L.T.; writing—original draft preparation, C.D., B.J. and L.T.; E.R. contributed to the text editing; writing—review and editing, C.D. and L.T.; visualization, C.D. and B.J.; supervision, C.D.; project administration, C.D.; funding acquisition, C.D., B.J. and L.T. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by RÉGION PAYS DE LA LOIRE (France), grant "Etoiles Montantes" (to C.D.), by the region of Pays de la Loire (ConnecTalent EPIALG project) and the Epicycle ANR project (ANR-19-CE20-0028-02) to LT.

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

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** All data are included in this article.

**Acknowledgments:** We are grateful to Samuel Menicot and Carine Pruvost for their technical assistance.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Falkowski, P.G.; Barber, R.T.; Smetacek, V. Biogeochemical Controls and Feedbacks on Ocean Primary Production. *Science* **1998**, *281*, 200–206. [CrossRef] [PubMed]

2. De Martino, A.; Meichenin, A.; Shi, J.; Pan, K.; Bowler, C. Genetic and Phenotypic Characterization of *Phaeodactylum tricornutum* (Bacillariophyceae) Accessions. *J. Phycol.* **2007**, *43*, 992–1009. [CrossRef]

3. De Martino, A.; Bartual, A.; Willis, A.; Meichenin, A.; Villazán, B.; Maheswari, U.; Bowler, C. Physiological and Molecular Evidence That Environmental Changes Elicit Morphological Interconversion in the Model Diatom *Phaeodactylum tricornutum*. *Protist* **2011**, *162*, 462–481. [CrossRef]

4. Chaumier, T.; Yang, F.; Manirakiza, E.; Ait-Mohamed, O.; Wu, Y.; Chandola, U.; Jesus, B.; Piganeau, G.; Groisillier, A.; Tirichine, L. Genome-Wide Assessment of Genetic Diversity and Transcript Variations in 17 Accessions of the Model Diatom *Phaeodactylum tricornutum*. *ISME Commun.* **2024**, *4*, ycad008. [CrossRef] [PubMed]

5. Rastogi, A.; Vieira, F.R.J.; Deton-Cabanillas, A.F.; Veluchamy, A.; Cantrel, C.; Wang, G.; Vanormelingen, P.; Bowler, C.; Piganeau, G.; Hu, H.; et al. A Genomics Approach Reveals the Global Genetic Polymorphism, Structure, and Functional Diversity of Ten Accessions of the Marine Model Diatom *Phaeodactylum tricornutum*. *ISME J.* **2020**, *14*, 347–363. [CrossRef]

6. Feng, Y.; Li, Z.; Li, X.; Shen, L.; Liu, X.; Zhou, C.; Zhang, J.; Sang, M.; Han, G.; Yang, W.; et al. Structure of a Diatom Photosystem II Supercomplex Containing a Member of Lhcx Family and Dimeric FCPII. *Sci. Adv.* **2023**, *9*, eadi8446. [CrossRef]

7. Grouneva, I.; Rokka, A.; Aro, E.M. The Thylakoid Membrane Proteome of Two Marine Diatoms Outlines Both Diatom-Specific and Species-Specific Features of the Photosynthetic Machinery. *J. Proteome Res.* **2011**, *10*, 5338–5353. [CrossRef]

8. Sandmann, G.; Reck, H.; Kessler, E.; Böger, P. Distribution of Plastocyanin and Soluble Plastidic Cytochrome c in Various Classes of Algae. *Arch. Microbiol.* **1983**, *134*, 23–27. [CrossRef]

9. Nagao, R.; Kato, K.; Ifuku, K.; Suzuki, T.; Kumazawa, M.; Uchiyama, I.; Kashino, Y.; Dohmae, N.; Akimoto, S.; Shen, J.-R.; et al. Structural Basis for Assembly and Function of a Diatom Photosystem I-Light-Harvesting Supercomplex. *Nat. Commun.* **2020**, *11*, 2481. [CrossRef]

10. Coesel, S.; Oborník, M.; Varela, J.; Falciatore, A.; Bowler, C. Evolutionary Origins and Functions of the Carotenoid Biosynthetic Pathway in Marine Diatoms. *PLoS ONE* **2008**, *3*, e2896. [CrossRef]

11. Doney, S.C.; Ruckelshaus, M.; Emmett Duffy, J.; Barry, J.P.; Chan, F.; English, C.A.; Galindo, H.M.; Grebmeier, J.M.; Hollowed, A.B.; Knowlton, N.; et al. Climate Change Impacts on Marine Ecosystems. *Ann. Rev. Mar. Sci.* **2012**, *4*, 11–37. [CrossRef]

12. Schuettengruber, B.; Bourbon, H.-M.; Di Croce, L.; Cavalli, G. Genome Regulation by Polycomb and Trithorax: 70 Years and Counting. *Cell* **2017**, *171*, 34–57. [CrossRef]

13. Glancy, E.; Ciferri, C.; Bracken, A.P. Structural Basis for PRC2 Engagement with Chromatin. *Curr. Opin. Struct. Biol.* **2021**, *67*, 135–144. [CrossRef]

14. Zhao, X.; Rastogi, A.; Deton Cabanillas, A.F.; Ait Mohamed, O.; Cantrel, C.; Lombard, B.; Murik, O.; Genovesio, A.; Bowler, C.; Bouyer, D.; et al. Genome Wide Natural Variation of H3K27me3 Selectively Marks Genes Predicted to Be Important for Cell Differentiation in *Phaeodactylum tricornutum*. *New Phytol.* **2021**, *229*, 3208–3220. [CrossRef]
15. De Lucia, F.; Crevillen, P.; Jones, A.M.E.; Greb, T.; Dean, C. A PHD-Polycomb Repressive Complex 2 Triggers the Epigenetic Silencing of FLC during Vernalization. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 16831-16836. [CrossRef]

16. Kim, J.; Bordiya, Y.; Xi, Y.; Zhao, B.; Kim, D.H.; Pyo, Y.; Zong, W.; Ricci, W.A.; Sung, S. Warm Temperature-Triggered Developmental Reprogramming Requires VIL1-Mediated, Genome-Wide H3K27me3 Accumulation in *Arabidopsis*. *Development* **2023**, *150*, dev201343. [CrossRef] [PubMed]

17. Goss, R.; Lepetit, B. Biodiversity of NPQ. *J. Plant Physiol.* **2015**, *172*, 13-32. [CrossRef] [PubMed]

18. Huysman, M.J.; Martens, C.; Vandepoele, K.; Gillard, J.; Rayko, E.; Heijde, M.; Bowler, C.; Inzé, D.; Van De Peer, Y.; De Veylder, L.; et al. Genome-Wide Analysis of the Diatom Cell Cycle Unveils a Novel Type of Cyclins Involved in Environmental Signaling. *Genome Biol.* **2010**, *11*, R17. [CrossRef] [PubMed]

19. Lin, X.; Tirichine, L.; Bowler, C. Protocol: Chromatin Immunoprecipitation (ChIP) Methodology to Investigate Histone Modifications in Two Model Diatom Species. *Plant Methods* **2012**, *8*, 48. [CrossRef]

20. Hattich, G.S.I.; Jokinen, S.; Sildever, S.; Gareis, M.; Heikkinen, J.; Junghardt, N.; Segovia, M.; Machado, M.; Sjöqvist, C. Temperature Optima of a Natural Diatom Population Increases as Global Warming Proceeds. *Nat. Clim. Chang.* **2024**, *14*, 518-525. [CrossRef]

21. Kléparski, L.; Beaugrand, G.; Edwards, M.; Schmitt, F.G.; Kirby, R.R.; Breton, E.; Gevaert, F.; Maniez, E. Morphological Traits, Niche-Environment Interaction and Temporal Changes in Diatoms. *Prog. Oceanogr.* **2022**, *201*, 102747. [CrossRef]

22. Allakhverdiev, S.I.; Kreslavski, V.D.; Klimov, V.V.; Los, D.A.; Carpentier, R.; Mohanty, P. Heat Stress: An Overview of Molecular Responses in Photosynthesis. *Photosynth. Res.* **2008**, *98*, 541-550. [CrossRef]

23. Wan, J.; Zhou, Y.; Beardall, J.; Raven, J.A.; Lin, J.; Huang, J.; Lu, Y.; Liang, S.; Ye, M.; Xiao, M.; et al. DNA Methylation and Gene Transcription Act Cooperatively in Driving the Adaptation of a Marine Diatom to Global Change. *J. Exp. Bot.* **2023**, *74*, 4259-4276. [CrossRef] [PubMed]

24. Yang, W.; Gao, S.; Bao, M.; Li, X.; Liu, Z.; Wang, G. HSP70A Promotes the Photosynthetic Activity of Marine Diatom *Phaeodactylum tricornutum* under High Temperature. *Plant J.* **2024**, *118*, 2085-2093. [CrossRef]

25. Murray, K.O.; Clanton, T.L.; Horowitz, M. Epigenetic Responses to Heat: From Adaptation to Maladaptation. *Exp. Physiol.* **2022**, *107*, 1144-1158. [CrossRef] [PubMed]

26. Debit, A.; Charton, F.; Pierre-Elies, P.; Bowler, C.; Cruz de Carvalho, H. Differential Expression Patterns of Long Noncoding RNAs in a Pleiomorphic Diatom and Relation to Hyposalinity. *Sci. Rep.* **2023**, *13*, 2440. [CrossRef] [PubMed]

27. Rastogi, A.; Lin, X.; Lombard, B.; Loew, D.; Tirichine, L. Probing the Evolutionary History of Epigenetic Mechanisms: What Can We Learn from Marine Diatoms. *AIMS Genet.* **2015**, *02*, 173-191. [CrossRef]

28. Zhao, P.; Gu, W.; Wu, S.; Huang, A.; He, L.; Xie, X.; Gao, S.; Zhang, B.; Niu, J.; Peng Lin, A.; et al. Silicon Enhances the Growth of *Phaeodactylum tricornutum* Bohlin under Green Light and Low Temperature. *Sci. Rep.* **2014**, *4*, 3958. [CrossRef] [PubMed]

29. Vartanian, M.; Desclés, J.; Quinet, M.; Douady, S.; Lopez, P.J. Plasticity and Robustness of Pattern Formation in the Model Diatom *Phaeodactylum tricornutum*. *New Phytol.* **2009**, *182*, 429-442. [CrossRef]

30. Jauffrais, T.; Jesus, B.; Méléder, V.; Turpin, V.; Russo, A.D.P.G.; Raimbault, P.; Jézéquel, V.M. Physiological and Photophysiological Responses of the Benthic Diatom *Entomoneis Paludosa* (Bacillariophyceae) to Dissolved Inorganic and Organic Nitrogen in Culture. *Mar. Biol.* **2016**, *163*, 115. [CrossRef]

31. Platt, T.; Gallegos, C.; Harrison, W. Photoinhibition of Photosynthesis in Natural Assemblages of Marine Phytoplankton. *J. Mar. Res.* **1980**, *38*, 687-701.

32. Serôdio, J.; Lavaud, J. A Model for Describing the Light Response of the Nonphotochemical Quenching of Chlorophyll Fluorescence. *Photosynth. Res.* **2011**, *108*, 61-76. [CrossRef] [PubMed]

33. Siaut, M.; Heijde, M.; Mangogna, M.; Montsant, A.; Coesel, S.; Allen, A.; Manfredonia, A.; Falciatore, A.; Bowler, C. Molecular Toolbox for Studying Diatom Biology in *Phaeodactylum tricornutum*. *Gene* **2007**, *406*, 23-35. [CrossRef] [PubMed]

34. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, Scalable Generation of High-quality Protein Multiple Sequence Alignments Using Clustal Omega. *Mol. Syst. Biol.* **2011**, *7*, 539. [CrossRef] [PubMed]

35. Kelley, L.A.; Mezulis, S.; Yates, C.M.; Wass, M.N.; Sternberg, M.J.E. The Phyre2 Web Portal for Protein Modeling, Prediction and Analysis. *Nat. Protoc.* **2015**, *10*, 845-858. [CrossRef] [PubMed]

36. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera? A Visualization System for Exploratory Research and Analysis. *J. Comput. Chem.* **2004**, *25*, 1605-1612. [CrossRef]

37. Paysan-Lafosse, T.; Blum, M.; Chuguransky, S.; Grego, T.; Pinto, B.L.; Salazar, G.A.; Bileschi, M.L.; Bork, P.; Bridge, A.; Colwell, L.; et al. InterPro in 2022. *Nucleic Acids Res.* **2023**, *51*, D418-D427. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.