Light Induced Changes in Pigment and Lipid Profiles of Bryopsidales Algae

Bryopsidales (Chlorophyta) are cultured and consumed in several regions of the planet and are known for their high nutritional value and bioprospection potential, due to a high content of relevant polar lipids and polysaccharides. Among other characteristic features, these marine algae are known for possessing unique photosynthetic pigment-protein complexes and for the absence (in nearly all species investigated) of a functional xanthophyll cycle, a ubiquitous photoprotection mechanism present in most algae and plants. With the aim of characterizing the photophysiology of this atypical group of algae, we investigated the changes in pigment content and polar lipidome of two Bryopsidales species (Codium tomentosum and Bryopsis plumosa) exposed for 7 days to low or high irradiance (20 vs. 1,000 μmol photons m–2 s–1). Our results show that high light has a strong effect on the pigment composition, triggering the time-dependent accumulation of all-trans-neoxanthin (t-Neo) and violaxanthin (Viola). High light-acclimated macroalgae also displayed a shift in the characteristic polar lipidome, including a trend of accumulation of lyso-glycolipids, and highly unsaturated phospholipids and betaine lipids. We hypothesize that the observed shifts on the lipid composition could promote the interaction between t-Neo and Viola with the siphonaxanthin–chlorophyll–protein complexes (SCP) of photosystem II (PSII) within the thylakoid membranes of the chloroplasts. Light induced changes in pigment and lipid composition could contribute to the fitness of Bryopsidales algae by reducing damages to the photosynthetic apparatus under increased irradiance.

Bryopsidales (Chlorophyta) are cultured and consumed in several regions of the planet and are known for their high nutritional value and bioprospection potential, due to a high content of relevant polar lipids and polysaccharides. Among other characteristic features, these marine algae are known for possessing unique photosynthetic pigmentprotein complexes and for the absence (in nearly all species investigated) of a functional xanthophyll cycle, a ubiquitous photoprotection mechanism present in most algae and plants. With the aim of characterizing the photophysiology of this atypical group of algae, we investigated the changes in pigment content and polar lipidome of two Bryopsidales species (Codium tomentosum and Bryopsis plumosa) exposed for 7 days to low or high irradiance (20 vs. 1,000 µmol photons m −2 s −1 ). Our results show that high light has a strong effect on the pigment composition, triggering the time-dependent accumulation of all-trans-neoxanthin (t-Neo) and violaxanthin (Viola). High light-acclimated macroalgae also displayed a shift in the characteristic polar lipidome, including a trend of accumulation of lyso-glycolipids, and highly unsaturated phospholipids and betaine lipids. We hypothesize that the observed shifts on the lipid composition could promote the interaction between t-Neo and Viola with the siphonaxanthin-chlorophyll-protein complexes (SCP) of photosystem II (PSII) within the thylakoid membranes of the chloroplasts. Light induced changes in pigment and lipid composition could contribute to the fitness of Bryopsidales algae by reducing damages to the photosynthetic apparatus under increased irradiance.

INTRODUCTION
Bryopsidales is a monophyletic order of siphonous green algae with high ecological relevance in both tropical and temperate marine coastal ecosystems (Cremen et al., 2019). Recently, these organisms have received increasing attention due to their remarkable content of valuable polar lipids and polysaccharides with potential pharmaceutical and nutraceutical applications (e.g., Maeda et al., 2012;Arata et al., 2015;Santos et al., 2015;Fleurence and Levine, 2016). Furthermore, some species of the genera Codium and Caulerpa are amongst the most popular edible macroalgae and listed as delicacies in gourmet cuisine (Paul et al., 2014;Pérez-Lloréns et al., 2018).
Despite their relevance on both ecological and economical levels, the photophysiology of Bryopsidales is poorly understood. The photosynthetic antennae of these macroalgae are characterized by the unique association of the xanthophylls siphonaxanthin and its esterified form siphonaxanthindodecenoate (also referred to as siphonein) at the level of the photosynthetic antennae, together with chlorophylls a and b Yokohama et al., 1977;Nakayama and Okada, 1990). These siphonaxanthin-type light harvesting complexes are responsible for an increased absorption in the green and blue-green light regions, predominant in deep waters or shaded subtidal marine habitats (Anderson, 1983;Wang et al., 2013). In recent years, it has also been reported that several species of Bryopsidales lack a functional xanthophyll cycle (Cruz et al., 2015;Christa et al., 2017;Giovagnetti et al., 2018), a common photoprotection mechanism present in most algae and higher plants (Latowski et al., 2011;Goss and Lepetit, 2015). The most common xanthophyll cycle comprises the sequential de-epoxidation of the pigments violaxanthin (Viola) to antheraxanthin and zeaxanthin under high light, allowing excess energy dissipation as heat. On the other hand, the accumulation of the xanthophylls all-trans-neoxanthin (t-Neo) and Viola, putatively involved in photoprotection, have been reported in Bryopsidales macroalgae under high light (Uragami et al., 2014;Cartaxana et al., 2018). However, the mechanisms of photoprotection necessary for these algae to cope with increased irradiance, particularly relevant in their natural habitats during low tides or summer, remain largely unexplored.
In this work, we investigated photophysiological changes in two different Bryopsidales species, Codium tomentosum (Stackhouse, 1797) and Bryopsis plumosa (Agardh, 1823), with the aim of identifying key players involved in photoprotection and/or in enhancing the photosynthetic performance under increased irradiance, by (i) characterizing the time-dependent changes in the pigment composition in response to high light and (ii) analyzing the changes in the polar lipidome in high light acclimated algae, focusing on lipid classes with known relevance for the structure and functionality of chloroplast membranes and molecules possibly involved in photoprotection. Our results show that high irradiance triggers time-dependent accumulation of the putative photoprotective xanthophylls t-Neo and Viola and a shift in the polar lipidome that could correspond to a rearrangement of chloroplast membranes, possibly allowing a better photosynthetic performance under high light. Potential interactions between the theorized photoprotective model involving t-Neo (Uragami et al., 2014) and the observed alterations in the polar lipidome are discussed.

Algae Collection and Cultivation
Wild specimens of C. tomentosum were collected from the rocky shores of Aguda beach, Portugal (41 • 04 66.22 N,8 • 65 32.72 W) in May 2019 during low tide. B. plumosa was acquired from Kobe University Macroalgal Culture Collection (KU-0990, KU-MACC, Japan). Prior to treatments, macroalgae were maintained in a thermostatically controlled room (17 • C), with aeration and long day photoperiod (16/8 h, light/dark) and an irradiance of 130 µmol photons m −2 s −1 (FLUORA T8, L58 W/77, OSRAM). Medium f/2 supplemented with 1 mg L −1 GeO 2 was prepared according to Andersen (2005). For experiments, algae were precultured under Low Light (LL, 20 µmol photons m −2 s −1 ) for 1 week. Subsequently, they were exposed for seven days to High Light (HL, 1,000 µmol photons m −2 s −1 ) or LL conditions. Light was provided with REEF-SPEC bluewhite T5 lamps (Red Sea). Algal samples were flash-frozen with liquid nitrogen, freeze-dried, and stored at −80 • C until further analysis. Five independent samples were collected per treatment.

Pigment Analysis
Pigment extraction, separation, and identification were performed as described in detail by Mendes et al. (2007). Briefly, approximately 4 mg of freeze-dried macroalgal samples were extracted in 1 mL of methanol-based solvent (95% methanol, 2% ammonium acetate), sonicated and incubated at −20 • C for 20 min in the dark. Pigment extracts were filtered through 0.2 µm PTFE membrane filters and immediately injected into a Prominence-i LC2030C HPLC system (Shimadzu, Japan). Chromatographic separation was performed with a Supelcosil C18 column (250 mm length; 4.6 mm diameter; 5 µm particles; Sigma-Aldrich, United States) for reverse phase chromatography and a 35 min elution program. The solvent gradient followed Kraay et al. (1992), with an injection volume of 50 µL and a flow rate of 0.6 mL min −1 . Pigments were identified from retention times and absorption spectra. Concentrations were calculated from the signal in the photodiode array (PDA) detector, using calibration curves constructed with pure crystalline standards (DHI group, Hørsholm, Denmark).

Lipid Extraction
Lipids were extracted using modified Bligh and Dyer method (Bligh and Dyer, 1959;Rey et al., 2020a), where chloroform was replaced with the less hazardous dichloromethane (Cequier-Saìnchez et al., 2008). Freeze-dried macroalgal biomass (50 mg) was homogenized with a mortar and pestle to obtain small flakes and then mixed with methanol (2.5 mL) and dichloromethane (1.25 mL) in a glass centrifuge tube. After vigorous homogenization and short sonication (1 min), samples were incubated on ice for 150 min on an orbital shaker (Stuart equipment, Bibby Scientific, Stone, United Kingdom). After centrifugation at 2,000 rpm for 10 min, the organic phase was transferred to another glass tube, followed by addition of water (2.25 mL) and dichloromethane (1.25 mL). Organic phase was collected after centrifugation at 2,000 rpm for 10 min. The aqueous phase was re-extracted with dichloromethane (2 mL), and organic phase was collected after centrifugation and combined with the first. Extraction was repeated on the original biomass, and the resulting organic phase was combined with the previous. Lipid extracts were dried under a nitrogen stream and preserved at −20 • C until further analysis.
Acquisition of LC-MS and MS/MS data was performed in positive and negative modes (electrospray voltage 3.0 and −2.7 kV, respectively), using positive/negative switching toggles, with capillary temperature of 350 • C and sheath gas flow of 20 U. In MS experiments, a high resolution of 70,000, automatic gain control (AGC) target of 1 × 10 6 , and maximum injection time (IT) of 100 ms were applied. For MS/MS determinations, a resolution of 17,500, AGC target of 1 × 10 5 and IT of 50 ms were used. The cycles consisted in one full scan mass spectrum, and 10 data-dependent MS/MS scans that were repeated continuously throughout the experiments, with dynamic exclusion of 60 s and intensity threshold of 2 × 10 4 . Normalized collision energy TM (CE) ranged between 20, 25, and 30 eV. Data acquisition was performed using the Xcalibur data system (V3.3, Thermo Fisher Scientific, United States).
Polar lipid species were identified by the m/z ratio of the ions observed in LC-MS spectra, mass accuracy of ≤5 ppm and typical retention time. Confirmation of lipid species was done by detailed analysis of MS/MS spectra. MS/MS fragmentation patterns characteristic for the polar lipid classes, including the polar head groups and fatty acyl chains, were previously described (Rey et al., 2017(Rey et al., , 2020a. Ion peak integration and assignments were performed using MZmine version 2.53, based on an in-house lipid database. Previous studies on C. tomentosum (Rey et al., 2020a,b) were consulted as reference lipidome. Peaks with raw intensity lower than 1 × 10 4 were excluded. Integrated peaks were normalized by dividing the peak areas of the extracted ion chromatograms (XIC) of each lipid species of each lipid class by the peak area of the internal standard selected for the class (Octave, Matlab).

Statistical Analysis
Changes in the total pigment pool were assessed by principal component analysis (PCA). Variance of the content of single pigments was tested using repeated measures partly nested ANOVA. Box-plots, interaction-plots and Pearson's residuals analyses were used to verify normality and homogeneity of variance and absence of block-by/within-block interactions. Sphericity assumption was verified using Mauchly's test. Pairwise post hoc comparisons were performed using Welch's two-samples t-test. These analyses were carried out using R (R Core Team, 2018).
Lipidomic data assessment was performed using MetaboAnalyst (Xia et al., 2015). Data filtering was performed using relative standard deviation (RSD = SD/mean). Normalized data were glog transformed and auto scaled. PCA were performed to visualize the general 2D clustering. Heat maps were performed using Euclidean distance measure and Ward clustering algorithm, considering the top 25 variables that most contributed to variability (PLS-DA vip). Differences in lipid content were analyzed with Welch's two-samples t-test (R Core Team, 2018).

Light-Induced Changes in Pigment Content
A total of nine characteristic photosynthetic pigments were identified in C. tomentosum and B. plumosa: chlorophylls a and b (Chl a, Chl b), siphonaxanthin, and its ester siphonaxanthin dodecenoate (Siph, Siph-do) as main light harvesting pigments; t-Neo, 9 -cis-neoxanthin (c-Neo), and Viola as typical xanthophylls; and ε,εand β,ε-carotenes (referred hereafter as ε-Car and α-Car, respectively) ( Table 1). The general light response was consistent in both species, with HL acclimation leading to significant changes in the composition of chlorophylls and carotenoids (Figure 1). Samples of B. plumosa acclimated to HL presented a drastic decline of the main light harvesting pigments (significantly different for Chl a, Chl b, Siph, and Siph-do), while the effect was less pronounced in C. tomentosum, with no significant change for Chl a content. Regarding carotenoids, both species presented highly significant accumulation of t-Neo and Viola under high light. Only in B. plumosa a decrement of pigment concentration per dried weight for c-Neo, εand α-Car was observed.
Principal component analysis (Figure 2) showed a clear separation between HL and LL conditions in both species, with x-and y-axis explaining, respectively 69.1 and 18.3% of variance in C. tomentosum, 76.5, and 11.9% in B. plumosa. In particular, LL samples formed a unique cluster with no difference between time points, while HL macroalgae formed separate clusters that gradually separate during 7 days of acclimation, with the samples at day 0 being the closest and day 7 the farthest from the LL cluster. ANOVA analysis (Supplementary Table 1) showed that both factors (irradiance and time) and their interaction had a significant effect on the variance of most pigment (expressed as

Time-Dependent Accumulation of All-Trans-Neoxanthin and Violaxanthin
In both macroalgal species, HL had the strongest effect on the concentrations of xanthophylls t-Neo and Viola. As shown in Figure 3, these pigments were accumulated during 7 days of acclimation to HL, with a positive correlation between their concentrations and the time of exposure. In particular, Viola was accumulated in significant concentrations already after 1 day, while t-Neo started accumulating from day 2, and kept increasing throughout the HL acclimation period. This pattern was consistent between C. tomentosum and B. plumosa. The lightand time-dependent accumulation effect was not observed for c-Neo, with concentrations fluctuating around the initial level without showing a particular pattern (Supplementary Figure 1). Pigment to Chl a ratios of c-Neo, Siph, Siph-do, Chl b, εand α-Car are reported in Supplementary Figure 1.
Coherently with these results, in the PCA analysis (Figure 2), t-Neo and Viola resulted as the variables that most affected the trend of differentiation of HL samples within the week of acclimation (showed by the respective red arrows pointing in the direction of the HL acclimated clusters). Accordingly, as showed by ANOVA (Supplementary Table 1) the effect of all considered factors (irradiance, time of exposure, and their interaction) was greatest on Viola and t-Neo compared to other pigments (p value <0.001 for all factors and interaction) in both species.

Polar Lipids Identification
The present study focused on the identification of lipid species belonging to the three classes of polar lipids, namely glycolipids, phospholipids, and betaine lipids, as previously reported for C. tomentosum (Rey et al., 2017(Rey et al., , 2020a. We identified a total of 275 lipid species in C. tomentosum and 191 in B. plumosa, belonging to glycolipids (86 in C. tomentosum; 64 in B. plumosa), phospholipids (142 in C. tomentosum; 120 in B. plumosa), and betaine lipids (47 in C. tomentosum; 7 in B. plumosa) classes.
In particular, C. tomentosum lipidome included: (i) for glycolipids, monogalactosyl diacylglycerols (MGDG), monogalactosyl monoacylglycerols ( The same lipid classes (except for SQMG and LPI) were also present in B. plumosa, but with a lower number of lipids species for each lipid class. The complete list of all identified molecular species with relative quantification and fatty acyl chains is reported in Supplementary Table 2.

Light-Induced Changes in the Polar Lipidome
Considering the general lipid profile of C. tomentosum (Figure 4), few quantitative changes could be observed under HL. Although most lipid classes did not present significant content alterations (normalized XIC areas), we observed a significant increment of the lyso-glycolipids MGMG ( Figure 4B) and lyso-betaine lipids MGTS under HL (Figure 4F).
On the other hand, B. plumosa presented more pronounced quantitative changes to several lipid classes (Figure 5): different glycolipids, DGDG and SQDG (Figure 5A), and phospholipids classes, PC and PE (Figure 5C), increased significantly under HL, while other phospholipids decreased (i.e., PG and LPG) (Figures 5C,D). B. plumosa also presented a significant increase of some betaine lipids (i.e., DGTS) ( Figure 5E). However, in this macroalgal species betaine lipids were generally present with extremely low concentrations in both treatments (Figures 5E,F). Changes under HL within relevant molecular species, analyzed with multivariate analyses (PCA and heatmap/clustering) by lipid class (i.e., glycolipids, phospholipids, and betaine lipids), are described next.

Glycolipids
In HL acclimated C. tomentosum, despite a clear clustering was not observed in the PCA performed on the whole glycolipid dataset (Figure 6A), a trend of differentiation emerged from the heatmap/clustering analysis of the 25 glycolipid molecular species that most contributed to discriminate between treatments ( Figure 6B). In particular, several lyso-glycolipids containing  Table 2). Similarly, other relevant species for FIGURE 1 | Pigment concentrations of low light (LL) and high light (HL) acclimated macroalgae at the beginning (day 0) and at the end (day 7) of the light acclimation (mean ± SD; n = 5). (A) Codium tomentosum; (B) Bryopsis plumosa. Asterisks indicate significant differences between LL and HL (t-tests): * p < 0.05; * * p < 0.01; * * * p < 0.001. For pigment abbreviations see Table 1.  Table 2). In B. plumosa, the glycolipid pool presented a more drastic response to HL, with a clear differentiation of the HL samples cluster in the PCA ( Figure 6C) and a decrease of all 25 most relevant molecular species (mostly represented by MGDG) in the heatmap/clustering analysis ( Figure 6D). The trend of MGDG species observed was opposite to the quantitative increase of other glycolipid classes (DGDG and SQDG) previously reported under HL ( Figure 5A).

Phospholipids
Principal component analysis, performed with the complete molecular species dataset of phospholipids, did not show a separation of HL acclimated C. tomentosum ( Figure 7A). Nonetheless, heatmap/clustering analysis performed on the 25 molecular species that most contributed to differentiation, showed a clear separation of HL samples after 7 days of acclimation ( Figure 7B). The phospholipids involved belonged mostly to PG and PC classes. In particular, it was observed an increase of molecular species with polyunsaturated fatty acids, often presenting 18:  (Figures 7B,D).

Betaine Lipids
In C. tomentosum, while PCA of the complete betaine lipids dataset did not show evident separation of HL cluster (Figure 8A) (as for glycolipids and phospholipids), heatmap/clustering analysis showed a general trend of increase of all 25 molecular species that most contribute to differentiation after 7 days of HL acclimation (Figure 8B). This analysis showed not only the increase of several lyso-betaine lipids (MGTS), coherently with previous observations (Figure 4F) Given that only seven molecular species of betaine lipids were identified in B. plumosa, with considerably low concentrations (Figures 5E,F and Supplementary Table 2), the relative data was excluded from this analysis.

Pigment Contribution to Excess Energy Regulation
Codium tomentosum and B. plumosa presented the expected pigment pool for Bryopsidales algae Nakayama and Okada, 1990;Cruz et al., 2015). In both species, acclimation to higher irradiance determined: (i) a general decrease of light harvesting pigments (namely Chls a and b, Siph, Siph-do), a known high light acclimation response in plants and algae and (ii) a change in the xanthophyll composition, involving a time-dependent accumulation of t-Neo and Viola. These results are coherent with previously reported studies performed on different species of Bryopsidales (Uragami et al., 2014;Cartaxana et al., 2018) and demonstrate a trend of acclimation to higher irradiances. Furthermore, the detailed characterization of the time-dependent changes in the pigment composition of both C. tomentosum and B. plumosa indicate that the response to high light can be faster (within one to two days of high light exposure) than previously reported (1-2 weeks) under lower irradiances (Uragami et al., 2014;Cartaxana et al., 2018). The accumulation of t-Neo and Viola under high light has been previously correlated with the oligomerization of siphonaxanthin-chlorophyll-protein complexes (SCP) of photosystem II (PSII) in Codium intricatum, leading to the hypothesis that these pigments might contribute to the dissipation of excessive light energy (Uragami et al., 2014). Indeed, according to this model the oligomerization of SCP would determine a decrease of energy transfer to PSII, by adjusting the distance between energy donor (SCP) and acceptor (PSII core). A precedent study also reported a correlation between oligomerization of antenna complexes (isolated from Bryopsis corticulans) and the dissociation of Chl b, reducing the light absorption capabilities of SCP (Chen et al., 2005). Therefore, this process could also represent an additional sink of excitation energy in the t-Neo-mediated oligomerization of SCP proposed by Uragami et al. (2014). In most algae and plant taxawith the only exception of the ancient green alga lineages of Mesostigmatophyceae (Yoshii et al., 2003) and Bryopsidales (Qin et al., 2015)--the cis isomer of neoxanthin (c-Neo) represents the sole molecular form associated with antenna complexes, while t-Neo is not associated with light harvesting or photoprotection (Takaichi and Mimuro, 1998). Hence, the involvement of t-Neo in photoprotection could represent a unique mechanism of Bryopsidales . The presence of such photoprotective activity could represent a safety valve for the dissipation of excess excitation energy under higher irradiances, particularly relevant for these algae that lack a functional xanthophyll cycle (Franklin et al., 1996;Raniello et al., 2004Raniello et al., , 2006Cruz et al., 2015;Christa et al., 2017;Giovagnetti et al., 2018), an important and ubiquitous photoprotection mechanism in aquatic and terrestrial oxygenic phototrophs (Latowski et al., 2011;Goss and Lepetit, 2015).
Our results also show the co-occurrent accumulation of t-Neo and Viola, with a relevant pattern. During one week of acclimation to HL, Viola is already accumulated at day 1, while t-Neo begins to be significantly accumulated only at day 2 (Figure 3). The accumulation of both these pigments clearly contributed to the differentiation of samples during HL acclimation, as showed in PCA analysis (Figure 2). Viola is known as the biochemical precursor of t-Neo, within the ßbranch of carotenoid biosynthetic pathway (Lohr, 2011), and this biochemical relationship explains the contextual accumulation of these pigments. It has also been proposed that Viola could participate in the energy tuning at the level of SCP together with t-Neo, given that both xanthophylls present the same weak interactions with the SCP (Uragami et al., 2014), possibly allowing fast and reversible interactions with the antenna complexes. Therefore, despite the role of Viola in the theorized photoprotective model remains unclear, the present work shows that this pigment is clearly involved in this photoacclimation process and is essential for the accumulation of t-Neo. During HL acclimation we also observed a putative accumulation of other carotenoids belonging to the same ß-branch of carotenoid biosynthesis (namely antheraxanthin, zeaxanthin, and ß-carotene), but the relative HPLC signal was too close to the detection limit to ensure their correct identification and quantification (data not shown). We propose that the synthesis of these carotenoids might also increase under higher irradiances FIGURE 4 | Sum of the normalized extracted ion chromatograms (XIC) areas (mean ± SD; n = 5) of the peaks identified for the lipids species of glycolipid, phospholipid, and betaine lipid classes in Codium tomentosum exposed to low light (LL) and high light (HL), at the beginning (day 0) and at the end (day 7) of the light acclimation. or longer acclimation than that used in this study, due to their role as biochemical precursors of Viola and t-Neo in the pathway of carotenoid biosynthesis (Lohr, 2011).

Patterns of Lipid Differentiation Under High Light
While the lipidome characterization obtained for C. tomentosum was consistent with previous studies (Da Costa et al., 2015;Rey et al., 2017Rey et al., , 2020a, we found no available literature information on the lipidome of B. plumosa. The number of lipid species identified was higher in C. tomentosum (275 lipids species) than in B. plumosa (191 lipid species). High light acclimation caused more significant changes in the lipidome of B. plumosa than that of C. tomentosum (Figures 4-7), which could be linked to the greater complexity of the polar lipidome of C. tomentosum (i.e., higher number of identified lipid species) determining higher plasticity during acclimation. Most likely, in the presence of such rich and plastic polar lipidome, the imposed HL treatment (1,000 µmol photons m −2 s −1 for 1 week) was not sufficient for triggering a more drastic response in this species. Also, our target Bryopsidales are characterized by completely different macromorphologies: C. tomentosum displays a thick thallus, divided in an external photosynthetic region that hosts the chloroplasts and an internal colorless region (called medulla), while B. plumosa is formed by several thin filaments, all originating from a central segment and forming different layers (Graham and Wilcox, 2000;Lee, 2008). All these features possibly make the observed light acclimation response species-specific: the polar lipidome of wild specimens of C. tomentosum (like those collected and cultured for our study) could be more efficient in responding to environmental changes compared to cultured strains of B. plumosa. Nonetheless, heatmap/clustering analyses allowed to highlight several meaningful trends in both target species.
Under HL conditions, C. tomentosum accumulated some phospholipid species (i.e., PG and PC) and betaine lipids presenting fatty acyl chains with higher level of unsaturation (Figures 7, 8). The presence of polyunsaturated fatty acids could be linked with a photoprotective role within the thylakoid membranes, since they could contribute to the scavenging of reactive oxygen species (ROS; Mène-Saffrané et al., 2009;Schmid-Siegert et al., 2016). A similar trend of polyunsaturated fatty acid increase was observed in B. plumosa's glycolipids and phospholipids (Figures 6, 7). Additionally, in C. tomentosum we reported an increase of lyso-forms of glycolipids (MGMG) FIGURE 5 | Sum of the normalized extracted ion chromatograms (XIC) areas (mean ± SD; n = 5) of the peaks identified for the lipids species of glycolipid, phospholipid, and betaine lipid classes in Bryopsis plumosa exposed to low light (LL) and high light (HL), at the beginning (day 0) and at the end (day 7) of the light acclimation.  (Figures 4, 6). Cleavage of MGDG into MGMG could be related to their role in ROS scavenging (Schmid-Siegert et al., 2016). A similar trend was observed for lyso-betaine lipids: MGTS increased under HL (Figures 4, 8), indicating that acclimation to higher irradiance could lead to cleavage of a fatty acyl chain. In the case of MGMG, we observed a trend of increase of C16 fatty acids. It is known that these fatty acyl chains are esterified in the sn-2 position of glycolipids in the plastid membranes of green algae (Li-Beisson et al., 2015). For this reason, the increase of C16 lyso-glycolipids may represent a tendency of cleavage of fatty acids in sn-1 position under HL.
Monogalactosyl diacylglycerols, DGDG, SQDG, PG, and DGTS are known to be important components of chloroplast membranes in green algae (Hölzl and Dörmann, 2007;Kobayashi et al., 2016;Jüppner et al., 2017;Li-Beisson et al., 2019). In particular, MGDG, DGDG, SQDG, and PG are considered "the essential glycerolipid quartet" in all photosynthetic eukaryotes, given that their presence is thought to be critical for maintaining an optimal photosynthetic yield (Boudière et al., 2014). Alterations in the content of the lipid species from these classes allow us to conclude that high light seemed to stimulate a remodeling of chloroplast membranes, possibly allowing a better response to increased irradiance. Finally, despite some species-specific differences, common differentiation trends could be identified in both C. tomentosum and B. plumosa, leading to the conclusion that high-light acclimation response of the polar lipidome presents some shared patterns between different Bryopsidales species.

Pigment-Lipid Interactions Within the Thylakoid Membranes
As recently reviewed (Goss and Latowski, 2020), lipids are thought to be essential for xanthophyll-mediated photoprotection, in particular in relation with the xanthophyll cycle. In diatoms and higher plants, an increase of MGDG was correlated with high solubilization of hydrophobic xanthophyll cycle pigments and of their respective rate of de-epoxidation, therefore contributing to photoprotection (Latowski et al., 2004;Goss et al., 2005). MGDG are considered essential for the re-organization of the membrane environment in which antenna complexes are immersed, allowing the interaction between xanthophyll cycle pigments and enzymes, thus contributing to the FIGURE 6 | Principal component analysis (PCA) and heatmap/clustering analysis of glycolipids profile discriminating macroalgae at low light (LL) and high light (HL) at the beginning (day 0) and at the end (day 7) of the light acclimation. (A) PCA 2D score plot-Codium tomentosum; (B) heatmap/clustering-C. tomentosum; (C) PCA 2D score plot-Bryopsis plumosa (D) heatmap/clustering-B. plumosa. PCA scores plot were performed using the complete glycolipid dataset (normalized XIC areas) while for heatmaps/clustering analysis the 25 molecular species that most contributed to differentiation between treatments were selected (PLS-DA vip).
formation of the energy-dissipative quenching (Goss and Latowski, 2020). In this model, xanthophyll cycle pigments can be found either in association with the antenna complexes or immersed in the MGDG shield located around the complexes (Lepetit et al., 2010;Schaller et al., 2010).
Despite Bryopsidales were reported to lack a functional xanthophyll cycle (e.g., Christa et al., 2017), similarities can be found between the proposed model of t-Neo mediated oligomerization of SCP and the formation of quenching sites within functionally detached antennae oligomers of diatoms and higher plants (Miloslavina et al., 2008;FIGURE 7 | Principal component analysis (PCA) and heatmap/clustering analysis of phospholipids profile discriminating macroalgae at low light (LL) and high light (HL) at the beginning (day 0) and at the end (day 7) of the light acclimation. (A) PCA 2D score plot-Codium tomentosum; (B) heatmap/clustering-C. tomentosum; (C) PCA 2D score plot-Bryopsis plumosa (D) heatmap/clustering-B. plumosa. PCA scores plot were performed using the complete phospholipids dataset (normalized XIC areas) while for heatmaps/clustering analysis the 25 molecular species that most contributed to differentiation between treatments were selected (PLS-DA vip). Holzwarth et al., 2009;Giossi et al., 2020). For these reasons, the rearrangement of the MGDG pool observed in our study could reflect significant changes happening in the thylakoid membrane environment of Bryopsidales, possibly influencing the interaction between pigment and proteins within the SCP. We hypothesize that alterations of the biochemical properties of the lipid bilayer could promote the interaction between the excess of t-Neo FIGURE 8 | Principal component analysis (PCA) and heatmap/clustering analysis of betaine lipids profile discriminating Codium tomentosum at low light (LL) and high light (HL) at the beginning (day 0) and at the end (day 7) of the light acclimation. (A) PCA 2D score plot; (B) heatmap/clustering analysis. PCA scores plot were performed using the complete betaine lipids dataset (normalized XIC areas) while for heatmaps/clustering analysis the 25 molecular species that most contributed to differentiation between treatments were selected (PLS-DA vip). and Viola molecules, likely dissolved in the thylakoid membrane, with SCP.
Contextually, the reported trend of increase in relevant membrane lipids with higher number of unsaturations could also contribute to increase the fluidity of the membrane environment, possibly allowing faster interaction of t-Neo and Viola with the SCP and contributing to the instauration of a rapid and reversible oligomerization and separation of the complexes. It is possible that all of these processes would contribute to the proposed t-Neo-mediated photoprotective model (Uragami et al., 2014).

CONCLUSION
Due to scarce data on the general photophysiology of Bryopsidales, the function of high light-dependent accumulation of t-Neo and Viola in this group of algae remains to be elucidated. Nonetheless, we clearly show the co-occurrence of increased levels of these pigments with a shift in the polar lipidome composition that could reflect a rearrangement of the thylakoid membranes under high light, possibly contributing to the onset of the theorized neoxanthin-mediated photoprotection at the level of the light harvesting complexes. Observed changes in pigment and lipid composition could reduce damages to the photosynthetic apparatus under increased irradiance, thus contributing to the overall fitness of Bryopsidales algae.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
CG, SC, and PC: conceptualization. CG, SC, RM, and PC: methodology-algae collection, cultivation, and experimental setup. PC: methodology-pigment analysis. FR, TM, and MRD: methodology-lipidomics. CG, FR, and PC: data analysis. CG: investigation and first manuscript draft. SC, FR, MRD, and PC: supervision. PC: project administration and funding acquisition. All authors contributed to manuscript revision, read, and approved the submitted version.