Seasonal Photoacclimation and Vulnerability Patterns in the Brown Macroalga Lessonia spicata (Ochrophyta)

: Fluctuations in solar radiation are one of the key factors affecting productivity and survival in habitat forming coastal macroalgae, in this regard, photoacclimation has a direct impact on the vulnerability and the capacity of seaweed to withstand, for instance, radiation excess. Here, we study ecophysiological responses through photosynthetic activity measurements under time-dependent (one year) ﬂuctuations in solar radiation in the brown macroalga L. spicata . The responses presented seasonal patterns, with an increase in photosynthetic capacity during summer, expressed in greater maximal electron transport rate (ETR max ) and diminished thermal dissipation (NPQ max ). Moreover, we studied photoprotective compounds (phenolic compounds) and total antioxidant capacity, which demonstrated an increase during periods of high solar radiation. In addition, content of photosynthetic pigment (Chl a , Chl c and Carotenoids) increased under greater solar irradiance. spicata can accumulate as reservoir photoprotective and antioxidant substances to withstand periods of high solar irradiance. All ecophysiological and biochemical responses in L. spicata indicate high photoacclimation and low vulnerability in the species, especially during with greater levels of solar irradiance.


Introduction
In the last century, according to the Intergovernmental Panel on Climate Change [1] each decade has been successively warmer than the previous one, illustrating one of the main consequences of global climate change. In addition to greater atmospheric CO 2 , there have been increasing patterns in the surface solar radiation resulting from the destruction of the ozone layer by chlorofluorocarbons (CFCs) and other substances that deplete atmospheric ozone [2]. In this context, the study of the photosynthetically active radiation (PAR, λ = 400-700 nm) and ultraviolet radiation (UVB, λ = 280-315 nm and UVA, λ = 315-400 nm) are important factors influencing primary production (photosynthetic organisms) and the ecological balance of marine organisms [3].
In central Chile, solar radiation varies throughout the year, with levels increasing from winter to summer [4]. Although there is some data on radiation levels in the area, there is no information on how photoautotrophic organisms respond upon increases and the effects on ecological balance, and how they will be affected in the medium to long term [5]. Radiation in central Chile can reach levels of up to 680.9 MJm −2 in summer [5]. Moreover, in the coastal waters off Central Chile, the coastal upwelling occurs in seasonal cycles and the strong southerly upwelling favors winds prevailing during spring-summer [6]. Therefore, Central Chile can be a helpful natural laboratory to carry out cutting-edge investigations on photobiology, photo-accumulation, and radiation tolerance in photosynthetic organisms. Oscillations in solar radiation can affect responses in photoacclimation or in the vulnerability of the sessile species as macroalgae [7,8].
It is known that the excess of radiation can disrupt the electron transport, causing an overproduction of reactive oxygen species (ROS), and in some cases, the presence of the solar radiation in low tidal zones, can induce oxidative stress and, in extreme cases, damage and cellular death [8][9][10]. Thus, it has been observed that one of the first metabolic processes after exposure to excess solar radiation in brown macroalgae corresponds to an impact on photosynthetic processes and on the induction in the biosynthesis of photoprotective and antioxidant compounds [11]. In this regard, macroalgae are considered vulnerable with respect to high ultraviolet radiation as it induces an imbalance between photodamage, photoprotection, and photo-repair mechanisms [9,11]. The concept of the vulnerability in macroalgae is associated with their photoacclimation and photo-accumulation abilities, defined as the concentration of photoprotective compounds synthesized under higher irradiance, which usually occurs during periods of low tide [7,12].
Stress caused by high solar radiation in brown macroalgae has been observed to enhance the production of different photoprotective compounds as polyphenols (also called phenolic compounds or phlorotannins) and carotenoids [5]. Polyphenols and carotenoids play an important role in dissipating energy and as antioxidants. The antioxidant properties of polyphenols are related to phenolic rings, which serve as electron donors for ROS as superoxide anions and hydroxyl radical [13]. Interestingly, the polyphenols in brown macroalgae have been found to have even stronger antioxidant capacities than in terrestrial plants [11,14]. In this context, it is well known that the use of photoprotective compounds and antioxidants extracted from brown seaweeds can be used for different biotechnological purposes, considering the pharmaceutical, cosmeceutical and food industries [15]; thus, it is certainly relevant to increase our knowledge on different macroalgae species to address basic responses and biotechnological potential.
Macroalgae have proven to be good biomonitoring organisms and represent ecosystem status, especially due to their sessile nature, widespread distribution, and responses to environmental disturbances [16]. However, despite the latter, the focus of most studies in Chile have been to develop evaluations associated with contamination [5,10,17,18], but there is scarce information on their physiological and biochemical responses with respect to stress, eventually caused by naturally fluctuating environmental conditions. Hence, in this investigation we studied the physiological and biochemical responses of Lessonia spicata (Surh) Santelices (Laminariales, Phaeophyceae, Ochrophyta), highly abundant brown macroalgae in the western Pacific Coast, and greatly relevant as key habitat-forming species in intertidal coastal ecosystems [5,19]. Therefore, we studied L. spicata vulnerability and photoacclimation, as well as the photosynthetic production, accumulation, and induction of photoprotection compounds and antioxidant capacity throughout seasonal fluctuations in solar irradiance.

Species and Sampling Location
Nine L. spicata thalli were collected at least 2 m distance among the thalli at 0.1-0.4 m above Chart Datum at 10:00 local Chilean time, on the 20th of each month from March 2019 to February 2020, throughout summer, autumn, winter, and spring. The samples were collected from the rocky shore at Cochoa Beach, Viña del Mar, Chile (32 • 57 19.0" S, 71 • 32 52.4" W). For photosynthetic evaluations, live material was transported in cooler boxes. For biochemical assessments, L. spicata blades were immediately frozen in liquid nitrogen, transported to the laboratory, and stored at −80 • C for further analyses.

Abiotic Parameters
The photosynthetically active radiation (PAR, 2019 to February 2020, throughout summer, autumn, winter, and spring. The samples were collected from the rocky shore at Cochoa Beach, Viña del Mar, Chile (32°57′19.0″ S, 71°32′52.4″ W). For photosynthetic evaluations, live material was transported in cooler boxes. For biochemical assessments, L. spicata blades were immediately frozen in liquid nitrogen, transported to the laboratory, and stored at −80 °C for further analyses.

Abiotic Parameters
The photosynthetically active radiation (PAR, ʎ = 400-700 nm), and ultraviolet A radiation (UVA, ʎ = 320-400 nm) were measured in the air and in the water, column using QSO-SUN 2.5V and USB-SU 100 cosine-corrected sensors (Apogee Instruments, Logan, UT, USA), respectively. Each PAR and UVA sensor was inside a polycarbonate box (Otter Box 3000) [20,21]. Water temperature was measured monthly with a waterproof temperature data logger (HOBO Pendant ® Temperature, Onset Data Loggers, Bourne, MA, USA). Salinity, pH, and conductivity were measured using a multiparameter water quality meter (HI 98194, Hanna Instruments, Woonsocket, RI, USA). In addition, PAR diffuse attenuation coefficient (KdPAR) for downward irradiance in the water column, was calculated measuring PAR radiation just below the water surface and every 0.5 m depth in the water column. The Beer-Lambert law was used to calculate the attenuation coefficients; the k values were expressed in m −1 .

Photosynthesis and Energy Dissipation Measurements
In vivo chlorophyll a fluorescence of Photosystem II was measured in blades using portable pulse amplitude modulated fluorometer (MINI-PAM II-Walz GmbH, Germany).
L. spicata basal fluorescence (Ft) and maximal fluorescence (Fm') were measured under light conditions to obtain the effective quantum yield (ΔF/Fm') being ΔF =Fm' − Ft/Fm'. To determine maximal quantum yield (Fv/Fm), algal thalli was incubated in a dark chamber with fresh seawater for 15 min before measuring rapid light curves to fully photo-reduce all reaction centers. Fv is the variable fluorescence as the difference between Fm and Fo, being Fm the maximal fluorescence and Fo the basal fluorescence of dark-adapted thalli.
The electron transport rate (ETR) was determined by exposing the tissue for 20 s to twelve irradiances of actinic light in the rapid light curve (RLC), E1 = 25, E2 = 45, E3 = 66, E4 = 90, E5 = 125, E6 = 190, E7 = 285, E8 = 420, E9 = 625, E10 = 845, E11= 1150 and E12 = 1500 µmol photons m −2 s −1 , by using the following formula by [22]: where ΔF/Fm' is the effective quantum yield, E is the PAR on each pulse of light expressed in µmol m −2 s −1 , A is the blade absorptance equivalent to the fraction of the incident irradiance that is absorbed by the algae [23], and FII is the portion of chlorophyll a related to PSII (400-700 nm), being 0.8 in brown macroalgae [9]. The ETR parameters as maximum electron transport rate (ETRmax), and the initial slope of ETR versus irradiance function (αETR) as index of photosynthetic efficiency, were calculated from the tangential model by [24]. Saturation irradiance (EkETR) was calculated from the intercept of ETRmax and αETR.
Non photochemical quenching (NPQ) was calculated as proposed by [25]: Finally, maximal non photochemical quenching (NPQmax), the initial slope of NPQ versus irradiance function (αNPQ) and the saturation of irradiance of the EKNPQ, were obtained from the hyperbolic tangent fitting models, according to [26] of NPQ versus irradiance. The yield of losses, Y(NO) and Y(NPQ) are determined according to [27] and [28], = 400-700 nm), and ultraviolet A radiation (UVA, 2019 to February 2020, throughout summer, autumn, winter, and spring. The samples were collected from the rocky shore at Cochoa Beach, Viña del Mar, Chile (32°57′19.0″ S, 71°32′52.4″ W). For photosynthetic evaluations, live material was transported in cooler boxes. For biochemical assessments, L. spicata blades were immediately frozen in liquid nitrogen, transported to the laboratory, and stored at −80 °C for further analyses.

Abiotic Parameters
The photosynthetically active radiation (PAR, ʎ = 400-700 nm), and ultraviolet A radiation (UVA, ʎ = 320-400 nm) were measured in the air and in the water, column using QSO-SUN 2.5V and USB-SU 100 cosine-corrected sensors (Apogee Instruments, Logan, UT, USA), respectively. Each PAR and UVA sensor was inside a polycarbonate box (Otter Box 3000) [20,21]. Water temperature was measured monthly with a waterproof temperature data logger (HOBO Pendant ® Temperature, Onset Data Loggers, Bourne, MA, USA). Salinity, pH, and conductivity were measured using a multiparameter water quality meter (HI 98194, Hanna Instruments, Woonsocket, RI, USA). In addition, PAR diffuse attenuation coefficient (KdPAR) for downward irradiance in the water column, was calculated measuring PAR radiation just below the water surface and every 0.5 m depth in the water column. The Beer-Lambert law was used to calculate the attenuation coefficients; the k values were expressed in m −1 .

Photosynthesis and Energy Dissipation Measurements
In vivo chlorophyll a fluorescence of Photosystem II was measured in blades using portable pulse amplitude modulated fluorometer (MINI-PAM II-Walz GmbH, Germany).
L. spicata basal fluorescence (Ft) and maximal fluorescence (Fm') were measured under light conditions to obtain the effective quantum yield (ΔF/Fm') being ΔF =Fm' − Ft/Fm'. To determine maximal quantum yield (Fv/Fm), algal thalli was incubated in a dark chamber with fresh seawater for 15 min before measuring rapid light curves to fully photo-reduce all reaction centers. Fv is the variable fluorescence as the difference between Fm and Fo, being Fm the maximal fluorescence and Fo the basal fluorescence of dark-adapted thalli.
The ETR parameters as maximum electron transport rate (ETRmax), and the initial slope of ETR versus irradiance function (αETR) as index of photosynthetic efficiency, were calculated from the tangential model by [24]. Saturation irradiance (EkETR) was calculated from the intercept of ETRmax and αETR. Non photochemical quenching (NPQ) was calculated as proposed by [25]: Finally, maximal non photochemical quenching (NPQmax), the initial slope of NPQ versus irradiance function (αNPQ) and the saturation of irradiance of the EKNPQ, were obtained from the hyperbolic tangent fitting models, according to [26] of NPQ versus irradiance. The yield of losses, Y(NO) and Y(NPQ) are determined according to [27] and [28], = 320-400 nm) were measured in the air and in the water, column using QSO-SUN 2.5V and USB-SU 100 cosine-corrected sensors (Apogee Instruments, Logan, UT, USA), respectively. Each PAR and UVA sensor was inside a polycarbonate box (Otter Box 3000) [20,21]. Water temperature was measured monthly with a waterproof temperature data logger (HOBO Pendant ® Temperature, Onset Data Loggers, Bourne, MA, USA). Salinity, pH, and conductivity were measured using a multiparameter water quality meter (HI 98194, Hanna Instruments, Woonsocket, RI, USA). In addition, PAR diffuse attenuation coefficient (K dPAR ) for downward irradiance in the water column, was calculated measuring PAR radiation just below the water surface and every 0.5 m depth in the water column. The Beer-Lambert law was used to calculate the attenuation coefficients; the k values were expressed in m −1 .

Photosynthesis and Energy Dissipation Measurements
In vivo chlorophyll a fluorescence of Photosystem II was measured in blades using portable pulse amplitude modulated fluorometer (MINI-PAM II-Walz GmbH, Germany).
L. spicata basal fluorescence (F t ) and maximal fluorescence (F m' ) were measured under light conditions to obtain the effective quantum yield To determine maximal quantum yield (F v /F m ), algal thalli was incubated in a dark chamber with fresh seawater for 15 min before measuring rapid light curves to fully photo-reduce all reaction centers. F v is the variable fluorescence as the difference between F m and F o , being F m the maximal fluorescence and F o the basal fluorescence of dark-adapted thalli.
The electron transport rate (ETR) was determined by exposing the tissue for 20 s to twelve irradiances of actinic light in the rapid light curve (RLC), E1 = 25, E2 = 45, E3 = 66, E4 = 90, E5 = 125, E6 = 190, E7 = 285, E8 = 420, E9 = 625, E10 = 845, E11= 1150 and E12 = 1500 µmol photons m −2 s −1 , by using the following formula by [22]: where ∆F/F m is the effective quantum yield, E is the PAR on each pulse of light expressed in µmol m −2 s −1 , A is the blade absorptance equivalent to the fraction of the incident irradiance that is absorbed by the algae [23], and F II is the portion of chlorophyll a related to PSII (400-700 nm), being 0.8 in brown macroalgae [9]. The ETR parameters as maximum electron transport rate (ETR max ), and the initial slope of ETR versus irradiance function (α ETR ) as index of photosynthetic efficiency, were calculated from the tangential model by [24]. Saturation irradiance (Ek ETR ) was calculated from the intercept of ETR max and α ETR .
Non photochemical quenching (NPQ) was calculated as proposed by [25]: Finally, maximal non photochemical quenching (NPQ max ), the initial slope of NPQ versus irradiance function (α NPQ ) and the saturation of irradiance of the EK NPQ , were obtained from the hyperbolic tangent fitting models, according to [26] of NPQ versus irradiance. The yield of losses, Y(NO) and Y(NPQ) are determined according to [27] and [28], where Y(NO) is the fraction of energy passively dissipated in form of heat and fluorescence, mainly due to closed PSII reaction centers. High values indicate inability to protect itself against photodamage by excess radiation.
Y(NPQ) is the fraction of energy dissipated in form of heat via regulated photoprotective NPQ mechanisms.

Biochemical Variables
The pigments were extracted from 20 mg of nitrogen-grounded fresh tissue with 1.5 mL of 90% acetone. Samples were incubated overnight at 4 • C in the darkness and then centrifuged at 16,200 g (NU-C200R, NuWind, Plymouth, MN, USA) for 10 min at 4 • C. Finally, for chlorophyll a (Chla) and c (Chlc), the absorbance was determined at 510, 630, 664, and 750 nm, according to [29], and for the total carotenoids the absorbance at 480 nm as proposed by [30]. These analyses were conducted using a microplate spectrophotometer (SPECTROstar Nano, BMG Labtech, Offenburg, Germany). Pigment concentrations were expressed as mg g −1 dry weight after determining the fresh to dry weight ratio (2.84 for L. spicata).
The polyphenols or phenolic compounds (PC) were measured using 0.25 grs of frozen fresh tissue and mixed with 2.5 milliliters (mL) of 80% methanol in 15 mL conical tubes. Mixtures were incubated overnight at 4 • C in the darkness with vigorous agitation and then centrifuged at 2253 g for 30 min at 4 • C. The supernatant was collected to measure PC content colorimetrically using Folin-Ciocalteu reagent [31] and phloroglucinol (1,3,5trihydroxybenzene, Sigma P-3502) as standard. Finally, the absorbance was determined at 760 nm with the microplate reader [7]. Phenolic concentration was expressed as mg g −1 dry weight.
Antioxidant activity DPPH (2,2-diphenyl-1-picrylhydrazyil) assay (i.e., EC 50 ) was estimated by reducing the stable free radical DPPH [32]. In this regard, the supernatant for PC can also be used for this analysis [7,12]. DPPH solution was prepared in 90% methanol in distilled water (90MeOH: 10H 2 O) to a concentration of 1.27 mM. The reaction was complete after 30 min in a dark room at 20 • C, and the absorbance was read at 517 nm using the microplate reader. A calibration curve made with DPPH was used to calculate the remaining concentration of DPPH in the reaction mixture after incubation. Finally, Trolox (6-Hydroxy-2,5,7,8-tetramethylchromane-2-carboxylic acid) (0 to 50 µM) was applied as the reference antioxidant. The results obtained were expressed as µmol TEAC (Trolox equivalent antioxidant capacity) g −1 DW.

Statistical Analysis
The significant effects between physiological and biochemical variables were analyzed using ANOVA according to [33]. This test was performed for L. spicata including month (one-way) as a fixed factor with 12 levels (mean ± SE, n = 9); the probability applied in the statistical analysis was p < 0.05. Homogeneity of variance was tested using Cochran tests and by visual inspection of the residuals. Student Newman-Keuls tests (SNK) were performed after significant ANOVA interactions. Homogeneity of variance and homoscedasticity were tested using Cochran's C-tests tests and by visual inspection of the residuals. All data conformed to homogeneity of variance. Analyses were performed by using SPSS v.21 (IBM, Endicott, New York, NY, USA). The general variation patterns between physiological and biochemical variables measured in L. spicata were explored using a multivariate approach. A principal coordinates analysis (PCO) was performed for this purpose based on Euclidean distance using PERMANOVA+ within PRIMER v.6 software package [34]. Each one of the variables was represented by an arrow in the ordination plot pointing to the samples displaying the highest amount of that compound. Pearson correlation coefficients were calculated and tested between all measured dependent variables using Sigma plot 12.0 (Systat software, Slough, Berkshire, UK).

Environmental Conditions
The seawater temperature ranged from 12.25-15.21 • C, pH from 7.90-8.01, salinity between 31.74 and 33.98 PSU (Table 1)  pointing to the samples displaying the highest amount of that compound. Pearson correlation coefficients were calculated and tested between all measured dependent variables using Sigma plot 12.0 (Systat software, Slough, Berkshire, UK).

Physiological Responses and Biochemical Variables
Principal coordinates analysis (PCO diagram; Figure 3) along the year, shows a positive correlation of the first axis (34.9% of total variation) with NPQmax, Chla and Chlc, Car and AlfaNPQ, being highest in the samples from springtime. On the other hand, ETRmax, PC, AA and EkETR were highest in the samples collected in summer and autumn, and they were positively correlated with this axis (23.5% of total variation, Figure 3). The small angles between the arrows are indicative of high positive correlation between the variables (Figure 3). ETR versus irradiance (rapid light curves) in the months reaching maximal (April) and minimal (July) values are presented in Figure 4A. In April, the maximal value of 88 µmol m −2 s −1 was reached, whereas in July only 30 µmol m −2 s −1 was reached ( Figure 4A). The NPQ was lower in April and increased in the wintertime or July ( Figure 4B). Yield II decreased with the irradiance, reaching the minimal values in April around 0.2, whereas the minimal values in July were around 0.1 ≈ 0.05 ( Figure 4C,D). In the same context, yield loss as Y(NO) was the same trend in both months, with high values of 0.58 in July ( Figure  4C,D). Finally, the maximal values of Y(NPQ) were 0.55 in July ( Figure 4C) and 0.7 in April ( Figure 4D). At irradiances lower than 100 µmol photons m −2 s −1 , values of Y(NO) were higher than Y(NPQ), whereas at irradiances higher than 400 µmol photons m −2 s −1 were the reverse.
The Fv/Fm varied significantly depending on time (p < 0.05) ( Figure 5A and Table S1). Fv/Fm increased in samples collected in wintertime. ETRmax increased in summer and autumn ( Figure 5B and Table S1). Also, αETR was significantly higher in winter, summer and

Physiological Responses and Biochemical Variables
Principal coordinates analysis (PCO diagram; Figure 3) along the year, shows a positive correlation of the first axis (34.9% of total variation) with NPQ max , Chla and Chlc, Car and Alfa NPQ , being highest in the samples from springtime. On the other hand, ETR max , PC, AA and Ek ETR were highest in the samples collected in summer and autumn, and they were positively correlated with this axis (23.5% of total variation, Figure 3). The small angles between the arrows are indicative of high positive correlation between the variables (Figure 3).
Water 2021, 13, x FOR PEER REVIEW 7 of 15 autumn ( Figure 6A and Table S1). The EkETR increased significantly in summertime (Figure 6B and Table S1). The NPQmax was significantly different throughout the year (p < 0.05) ( Figure 7 and Table S1). The NPQmax increased under low temperature and low irradiance conditions in wintertime (Figure 7).  ETR versus irradiance (rapid light curves) in the months reaching maximal (April) and minimal (July) values are presented in Figure 4A. In April, the maximal value of  Figure 4A). The NPQ was lower in April and increased in the wintertime or July ( Figure 4B). Yield II decreased with the irradiance, reaching the minimal values in April around 0.2, whereas the minimal values in July were around 0.1 ≈ 0.05 ( Figure 4C,D). In the same context, yield loss as Y(NO) was the same trend in both months, with high values of 0.58 in July ( Figure 4C,D). Finally, the maximal values of Y(NPQ) were 0.55 in July ( Figure 4C) and 0.7 in April ( Figure 4D). At irradiances lower than 100 µmol photons m −2 s −1 , values of Y(NO) were higher than Y(NPQ), whereas at irradiances higher than 400 µmol photons m −2 s −1 were the reverse.   The F v /F m varied significantly depending on time (p < 0.05) ( Figure 5A and Table S1). F v /F m increased in samples collected in wintertime. ETR max increased in summer and autumn ( Figure 5B and Table S1). Also, α ETR was significantly higher in winter, summer and autumn ( Figure 6A and Table S1). The Ek ETR increased significantly in summertime ( Figure 6B and Table S1). The NPQ max was significantly different throughout the year (p < 0.05) (Figure 7 and Table S1). The NPQ max increased under low temperature and low irradiance conditions in wintertime (Figure 7).
Chla and Chlc contents were significantly different throughout the sampling period (p < 0.05) (Figure 8 and Table S2). The significant differences revealed clear differences along time, which were that Chla, Chlc and carotenoids were higher in springtime compared to the rest to the year (Figure 8).     Chla and Chlc contents were significantly different throughout the sampling period (p < 0.05) (Figure 8 and Table S2). The significant differences revealed clear differences along time, which were that Chla, Chlc and carotenoids were higher in springtime compared to the rest to the year (Figure 8).
The PC and AA presented significant differences throughout the year (Figure 9, Table  S2). PC decreased in wintertime and increased in spring, presenting the highest values in summer ( Figure 9A). Similar values presented the antioxidant capacity with the greatest levels in summertime ( Figure 9B).
The Pearson correlation showed a positive correlation between PC and AA, throughout all seasons and months (Table S3). Furthermore, a positive correlation was found between αETR, EkETR, Fv/Fm, and ETRmax, and between ETRmax and antioxidant activity. Finally, NPQmax, pigment content and photoprotective and antioxidant activity were positively correlated (Table S3).  The PC and AA presented significant differences throughout the year (Figure 9, Table S2). PC decreased in wintertime and increased in spring, presenting the highest values in summer ( Figure 9A). Similar values presented the antioxidant capacity with the greatest levels in summertime ( Figure 9B).

Discussion
Photophysiological, biochemical and biological responses in habitat-forming species L. spicata presented seasonally dependent variation caused most likely due to fluctuations in solar radiation. We found that the maximum electron transport rate (ETRmax), phenolic compounds, antioxidant activity, and saturation of irradiance in ETR (EkETR) increased in summer, whereas NPQmax, Chla, Chlc, and carotenoids increased in springtime (Figure 3). Our results suggest that the light stress in summertime can provoke photobiological and biochemical changes in L. spicata specimens, related to effective photoprotection of photosynthesis. The PAR diffuse attenuation coefficient (KdPAR) for downward irradiance in the water column, was 0.360 m −1 in autumn, 0.288 m −1 in winter, 0.304 m −1 in spring and 0.336 m −1 in summer, indicating less light penetration in the water column in the months with greater solar irradiance, such as a spring and summer; in contrast the highest light penetration occurred in wintertime. These variations can be explained by seasonal cycles with strong southerly upwelling due to increased winds prevailing during winter, and because nutrient content is limited in summer and can be also consumed by the phytoplankton [e.g., 6].
In this study, the maximum quantum yield (Fv/Fm) was higher at the end the winter, with values around 0.73, and at the beginning of springtime; in this regard, Fv/Fm is a useful parameter for a quick assessment of the photosynthetic conditions [9,35] and is used as an indicator of photoinhibition on primary producers [36]. [5] show values of Fv/Fm in L. spicata of 0.76. Here, L. spicata specimens displayed a reduction around 15% in the Fv/Fm, in spring and summer, indicating that the alga suffered photoinhibition as it has been The Pearson correlation showed a positive correlation between PC and AA, throughout all seasons and months (Table S3). Furthermore, a positive correlation was found between α ETR , Ek ETR , F v /F m , and ETR max, and between ETR max and antioxidant activity. Finally, NPQ max , pigment content and photoprotective and antioxidant activity were positively correlated (Table S3).

Discussion
Photophysiological, biochemical and biological responses in habitat-forming species L. spicata presented seasonally dependent variation caused most likely due to fluctuations in solar radiation. We found that the maximum electron transport rate (ETR max ), phenolic compounds, antioxidant activity, and saturation of irradiance in ETR (E kETR ) increased in summer, whereas NPQ max , Chla, Chlc, and carotenoids increased in springtime (Figure 3). Our results suggest that the light stress in summertime can provoke photobiological and biochemical changes in L. spicata specimens, related to effective photoprotection of photosynthesis. The PAR diffuse attenuation coefficient (K dPAR ) for downward irradiance in the water column, was 0.360 m −1 in autumn, 0.288 m −1 in winter, 0.304 m −1 in spring and 0.336 m −1 in summer, indicating less light penetration in the water column in the months with greater solar irradiance, such as a spring and summer; in contrast the highest light penetration occurred in wintertime. These variations can be explained by seasonal cycles with strong southerly upwelling due to increased winds prevailing during winter, and because nutrient content is limited in summer and can be also consumed by the phytoplankton e.g., [6].
In this study, the maximum quantum yield (F v /F m ) was higher at the end the winter, with values around 0.73, and at the beginning of springtime; in this regard, F v /F m is a useful parameter for a quick assessment of the photosynthetic conditions [9,35] and is used as an indicator of photoinhibition on primary producers [36]. [5] show values of F v /F m in L. spicata of 0.76. Here, L. spicata specimens displayed a reduction around 15% in the F v /F m , in spring and summer, indicating that the alga suffered photoinhibition as it has been reported in other brown macroalgae species [7,9,37]. Photoinhibition of photosynthesis is usually linked to a reduced photosynthetic production or maximal electron transport rate (ETR max ) [37,38], as well as related to limited availability of nutrients in seawater [39]. Here, the photosynthetic activity (ETR max ), as index of productivity, which was higher in summertime, with almost 60-100 µmol m −2 s −1 . These results indicate that under higher irradiance of summer there may be decreased excitation transfer to the reaction centers due to photoinhibition. Similar results were found in Cystoseira tamariscifolia [7,12] in spring and summer, with higher values in respect to the wintertime. The latter suggests that macroalgae may withstand radiation derived stress in springtime and decrease the vulnerability in wintertime. Indeed [21], showed similar results in Gelidium corneum plants, which were transplanted from areas with low to high irradiances (i.e., 3 to 5 m of depth in the water column).
In the same trend, the α ETR and Ek ETR increased in months with high solar irradiance, indicating a photoacclimation pattern typical of sun affinity algae [12]. In addition, the decrease in ETR max observed in L. spicata throughout the year could be partially explained by increased energy demand from other photoprotection cellular processes [12,40,41]. In this context, it is known that photosynthetic activity can enhance the accumulation of PC in other brown macroalgae [42] as well as in the green macroalga Ulva rigida [43]. Indeed, the NPQ max dissipation in the form of heat, used as a photoprotection indicator, also associated to photoinhibition and the concentrations of the pigments, presented the highest trend of increase in springtime, this suggests, an increase of the nutrient fertilization, product by seasonal upwelling [6].
The energy dissipation pattern or yield loss in L. spicata at higher irradiances than 400 µmol photons m −2 s −1 was produced under regulated photoprotective mechanisms, i.e., Y(NPQ) values were higher than Y(NO). In the summer period, the ETR presented the high values of the productivity index and the NPQ was lower for an increase in wintertime (July-August). The NPQ represents a loss of potential energy by the system probably involving several factors such as thermodynamic losses as waste heat and losses in the xanthophyll cycle of photosynthesis in some phototrophs [26]; thus, the contribution of non-photoregulated dissipation Y(NO) in the yield losses increased compared to April, but it remains lower than that of Y(NPQ). Y(NPQ) was higher under nutrient-enriched periods compared to low conditions. The high values of Y(NPQ) are indicative of a high photoprotective capacity, whereas high values of Y(NO) reflect the inability of a plant to protect itself against damage by excess illumination [44]. Indeed, similar trends were shown by [44] in the green alga Ulva rigida. The authors evidenced that the differences in energy losses, and the photosynthetic conditions are more favorable under high nutrients conditions, than low carbon and nitrate. A decrease in chlorophylls and carotenoids was measured in L. spicata individuals in the summertime, which may also be related to low levels of environmental nitrogen as it has been suggested for other macroalgae [45].
In our study, we found that higher internal content of the photoprotectors such as phenolics compounds in summer-season has a relation with photoacclimation patterns. The phenol and antioxidant capacity of this seaweed increased in spring, a period concomitant with high solar irradiance and nutrient in surface waters due to upwelling. L. spicata could accumulate nitrogen during winter and spring as a reservoir for periods of high irradiance and need of photoprotective mechanisms; for instance, as observed in summer in C. tamariscifolia, [7,9,12]. The PC, as photoprotective substances, presented a positive correlation with antioxidant capacity; this response suggests an increase in the process of energy dissipation and antioxidant control under great solar irradiance. High levels of PC were also detected in C. tamariscifolia in the Mediterranean Sea during high solar irradiance in springtime, which was attributed to increase of PAR and UV radiation [7]. Also considering that the process of acclimation, due to high levels of radiation is dependent on greater non-photochemical quenching, photoinhibition, and the accumulation of UV screening photoprotector [9,46].
In addition, it has been observed that L. spicata displays an increase in PC upon excess solar radiation, evidenced after 48 h of UV and PAR exposure in the laboratory [5]; similar observations have been made in Macrocystis pyrifera, a cosmopolitan macroalga and a species with a key ecological role [47][48][49]. Records by [7] in the brown macroalga C. tamariscifolia revealed higher total antioxidant capacity during summer (average 10,165 kJ m −2 PAR, 1051 kJ m −2 UVA, and 57.5 kJ m −2 UVB) compared to daily doses of winter (4000 kJm −2 PAR, 400 kJ m −2 UVA, and 15 kJ m −2 UVB). In this respect, we found that the highest antioxidant capacity in L. spicata during summer is accordance with PAR and UVA levels of 8000-9000 kJ m −2 and 500-600 kJ m −2 , respectively, also accompanied by higher content of PC. This information suggests a proficient photoacclimation strategy to reduce physiological vulnerability. Overall, the photoprotective substances as phenolic compound content was correlated with antioxidant capacity, suggesting an important role of PC in the photoprotection, not only as UV-screening substances, but also as antioxidant compounds [50][51][52].
General responses in L. spicata demonstrate the development of high photoacclimation adaptations, as evidenced during months with increased solar irradiance. This must be considered when reducing the potential negative effects of solar irradiance on photosynthetic activity. The data showed that L. spicata can produce UV-screening substances and photosynthetic pigments and present high phenotypical plasticity subject to yearly fluctuations in solar irradiance. Thus, the real consequence in these oscillations directly provoke responses inside and they can modulate their photoacclimation or the vulnerability of these sessile species such as other macroalgae [7,12].

Conclusions
In this study, L. spicata displayed a high photoacclimation capacity, evidenced in increased productivity and dynamic photoinhibition, together with greater UV photoprotections systems, especially during the high irradiance periods of spring and summer. Thus, the phenolics compounds and antioxidant capacity were generally greater throughout these seasons, confirming the photoacclimation strategy to avoid vulnerability in L. spicata.