Seasonal Acclimation Modulates the Impacts of Simulated Warming and Light Reduction on Temperate Seagrass Productivity and Biochemical Composition

Increases in seawater temperature and reduction in light quality have emerged as some of the most important threats to marine coastal communities including seagrass ecosystems. Temperate seagrasses, including Zostera marina, typically have pronounced seasonal cycles which modulate seagrass growth, physiology and reproductive effort. These marked temporal patterns can affect experimental seagrass responses to climate change effects depending on the seasons of the year in which the experiments are conducted. This study aimed at evaluating how seasonal acclimatization modulates productivity and biochemical responses of Zostera marina to experimental warming and irradiance reduction. Seagrass shoots were exposed to different temperatures (6, 12, 16, 20, and 24°C), combined with high (180 μmol photons m–2 s–1) and low (60 μmol photons m–2 s–1) light conditions across four seasons (spring: April, summer: July, and autumn: November 2015, and winter: January 2016). Plants exhibited similar temperature growth rates between 16 and 20°C; at 24°C, a drastic reduction in growth was observed; this was more accentuated in colder months and under low irradiance conditions. Higher leaf growth rates occurred in winter while the largest rhizomes were reached in experiments conducted in spring and summer. Increases in temperature induced a significant reduction in polyunsaturated fatty acids (PUFA), particularly omega-3 (n-3 PUFA). Our results highlight that temperate seagrass populations currently living under temperature limitation will be favored by future increases in sea surface temperature in terms of leaf and rhizome productivity. Together with results from this study on Z. marina from a temperate region, a wider review of the reported impacts of experimental warming indicates the likely reduction in some compounds of nutritional importance for higher trophic levels in seagrass leaves. Our results further demonstrate that data derived from laboratory-based studies investigating environmental stress on seagrass growth and acclimation, and their subsequent interpretation, are strongly influenced by seasonality and in situ conditions that precede any experimental exposure.

Increases in seawater temperature and reduction in light quality have emerged as some of the most important threats to marine coastal communities including seagrass ecosystems. Temperate seagrasses, including Zostera marina, typically have pronounced seasonal cycles which modulate seagrass growth, physiology and reproductive effort. These marked temporal patterns can affect experimental seagrass responses to climate change effects depending on the seasons of the year in which the experiments are conducted. This study aimed at evaluating how seasonal acclimatization modulates productivity and biochemical responses of Zostera marina to experimental warming and irradiance reduction. Seagrass shoots were exposed to different temperatures (6,12,16,20, and 24 • C), combined with high (180 µmol photons m −2 s −1 ) and low (60 µmol photons m −2 s −1 ) light conditions across four seasons (spring: April, summer: July, and autumn: November 2015, and winter: January 2016). Plants exhibited similar temperature growth rates between 16 and 20 • C; at 24 • C, a drastic reduction in growth was observed; this was more accentuated in colder months and under low irradiance conditions. Higher leaf growth rates occurred in winter while the largest rhizomes were reached in experiments conducted in spring and summer. Increases in temperature induced a significant reduction in polyunsaturated fatty acids (PUFA), particularly omega-3 (n-3 PUFA). Our results highlight that temperate seagrass populations currently living under temperature limitation will be favored by future increases in sea surface temperature in terms of leaf and rhizome productivity. Together with results from this study on Z. marina from a temperate region, a wider review of the reported impacts of experimental warming indicates the likely reduction

INTRODUCTION
Seagrasses play a vital role in coastal zones in terms of productivity, support of biodiversity, constitution as a primary food source of nutrition for several marine organisms, protection of the coastline, carbon storage and nutrient retention (Hemminga and Duarte, 2000;Carruthers et al., 2007;Nordlund et al., 2018). However, over the last few decades, global seagrass distribution has been reduced by around 30%. Particularly Zostera marina is lost at a rate of 1.3% per year due to anthropogenic pressures and climate change effects (Orth et al., 2006;Waycott et al., 2009). Noteworthy, a recent study revealed signs of recovery of European seagrass meadows due to conservation efforts (de los Santos et al., 2019).
Zostera marina L. is the dominant habitat-forming seagrass species in the northern hemisphere, currently distributed from subtropical regions in the Pacific coast of Mexico (24.3 • N) and in the Mediterranean Sea (35.1 • N), to sub-Arctic regions in Greenland (64.2 • N) and in the northern coast of Russia (74.3 • N) covering a latitudinal distribution of ∼50 • (Short et al., 2006). Reflecting its broad distribution, this species is adapted to a wide range of environmental conditions including annual average temperatures ranging from ∼1 to 25 • C (Tyberghein et al., 2012). Temperate seagrasses typically have pronounced seasonal cycles which modulate seagrass growth, physiology and flowering (e.g., Orth and Moore, 1986;Alcoverro et al., 2001). Temperate Z. marina populations usually grow faster and develop larger photosynthetic structures in summer months under more suitable climate conditions, storing high reserves of energetic compounds in the rhizomes such as carbohydrates. In less favorable environmental conditions, plants reduce their aboveground structures to reduce respiratory demands, and use the energetic reserves (e.g., Olesen and Sand-Jensen, 1993;Lee et al., 2005).
Temperature is one the most important environmental factors controlling seagrass physiological, growth and reproductive processes (Lee et al., 2007;Qin et al., 2020b). Global warming has emerged as an important threat for marine coastal species including seagrass ecosystems (e.g., Duarte, 2002). Seagrass responses to predicted increases in temperature, including Z. marina, depends on population-specific thermal tolerance and their capacity to adapt to new local climate conditions (Short and Neckles, 1999). As a result, expected global warming may drive different effects within Z. marina populations across their latitudinal distribution range. For instance, some subarctic populations were reported to exhibit more favorable responses to temperature increases Olesen et al., 2015;Beca-Carretero et al., 2018a). On the contrary, warming may negatively affect southern Z. marina populations when summer temperatures can exceed their physiological optima (e.g., Reusch et al., 2005;Lopez-Calderon et al., 2016). Physiological and biochemical responses of centrally distributed populations to warming events remain more unclear (Beca-Carretero et al., 2018a;Aoki et al., 2020;Dubois et al., 2020).
Irradiance is the most relevant environmental driver affecting photosynthetic activity and the vertical distribution of seagrasses (Ralph et al., 2007). The predicted increase in the frequency and intensity of storms alongside continuous coastal modifications will likely cause a reduction in underwater light and thus, cause low light stress in seagrass populations (Duarte, 2002;Silva et al., 2013;Bertelli and Unsworth, 2018). While some species can tolerate low light for short periods, longer exposure can cause chronic physiological stress and eventual seagrass loss (e.g., Short and Wyllie-Echeverria, 1996;Nguyen et al., 2021). Global warming is expected to enhance light limitation stress in seagrasses as minimum irradiance requirements, to maintain adequate growth and physiological performance, are increased at higher temperatures (Lee et al., 2007). Experimental studies of Z. marina populations distributed at their central distribution range reported that plants exposed to light-limited conditions (50 µmol photons m −2 s −1 ) were highly vulnerable to temperatures above 20 • C, causing growth inhibition and negative carbon budgets (Beca-Carretero et al., 2018b). Recently, in situ studies provide evidence for largescale declines in Z. marina meadows due to combined effects of light reduction due to precipitation event sand sediment runoff and anomalous warming episodes of temperatures above 28 • C (Johnson et al., 2021).
Warming episodes are well documented for summer seasons, however, with the substantial change in climate in northern latitudes, they also occur during cold seasons in higher latitudes where Z. marina is present (Short et al., 2007;Bokhorst et al., 2009). Anomalous climate events throughout the year can alter seasonal rhythms, uncoupling plant internal clocks and, consequently, affecting their vegetative development, population dynamics and reproductive cycles (Alcoverro et al., 2001;Bokhorst et al., 2009;Bjerke et al., 2011). As a result, different growth and physiological seagrass responses to climate change effects can be expected depending on the seasons of the year in which these events occur. However, the majority of the controlled temperature exposure experiments using seagrasses were conducted within a single season (e.g., review in Beca-Carretero et al., 2018a). Some studies reported differential sensitivities to temperature and irradiance treatments across seasons in the temperate species Z. marina (Staehr and Borum, 2011), which makes it highly relevant to assess the potential effects of warming during different seasons. While several studies have reported effects of warming or change in light climate on seagrass ecosystems (review in Koch et al., 2013;review in Roca et al., 2016;review in Nguyen et al., 2021), only few authors assessed the potential synergic effects of multiple stressors (e.g., Koch et al., 2007;Ontoria et al., 2019b).
Irish coasts are characterized by the presence of Z. marina meadows in subtidal areas reaching maximum depths of 4-6 m, and in lesser extension in the lower intertidal (Madden et al., 1993;Beca-Carretero et al., 2019b). In Ireland, there is an overall lack of seagrass research including their distribution, ecology and healthy status (Dale et al., 2007;Wilkes et al., 2017). Recent studies reported that large subtidal Z. marina meadows remain undisturbed along the Irish coast (Beca-Carretero et al., 2019b;Cott et al., 2021), and it is expected that their distribution is wider than previously known (Beca-Carretero et al., 2020;Hastings et al., 2020). Ireland's climate is defined as a temperate oceanic climate, characterized by defined climatic seasons with a lack of extreme cold or warm temperatures. Western Irish seagrasses are exposed to minimum temperatures of 4-5 • C in winter and a maximum of 16-17 • C reached in summer. With optimal temperature for growth of worldwide Z. marina populations ranging from 15.3 to 24 • C (review in Lee et al., 2007;review in Beca-Carretero et al., 2018b), it can be expected that an increase of 2-3 • C in Irish nearshore systems will favor seagrasses growth and production (Beca-Carretero et al., 2019b).
Seagrasses constitute an important source of food for several marine organisms; important nutritional compounds include proteins, carbohydrates or lipids and fatty acids (review in Kim et al., 2021). Despite the ecological importance of seagrasses as a food source to higher tropic levels, only few studies have explicitly assessed the role of environmental changes on their nutritional value. Recently, studies of the seagrass species Z. marina, P. oceanica and Halophila stipulacea reported changes in essential nutritional compounds including reductions of polyunsaturated fatty acids (PUFAs) relative to saturated fatty acids (SFAs), lower leaf fiber and proteins, and an increase in leaf necrosis in response to increasing temperatures have been documented (e.g., Beca-Carretero et al., 2018bHernán et al., 2019;Nguyen et al., 2020a). Such biochemical changes render seagrasses exposed to stress conditions a less desirable food source for herbivores (Hernán et al., 2017) although it is currently not clear how warming can modify biochemical composition and associated nutritional value. Here, we assessed how combined temperature and irradiance exposure treatments affect centrally distributed Z. marina populations descriptors (morphological, production and biochemical composition), and how such responses are modulated by seasonal acclimatization. Specifically, we investigated, (i) the effects of a temperature increase by 3-4 • C above in situ temperatures on seagrass performance; (ii) the synergistic effects of thermal stress and light reduction on seagrass traits and survival; (iii) the extent to which any potential seasonal acclimatization affects the experimental optima for growth and productivity patterns and; finally (iv) the potential impact of experimental increases in temperature on the biochemical composition and nutritional value of seagrass leaves reported in the literature in comparison with results from this study. To address these questions, Z. marina shoots from Ireland were exposed to different temperature treatments (ranging from 6 to 24 • C) combined with high (180 µmol photons m −2 s −1 ) versus low irradiance (60 µmol photons m −2 s −1 ) levels across four seasons (spring: April, summer: July, and autumn: November 2015, and winter: January 2016). We hypothesize that, (i) predicted warming will favor growth of centrally distributed Z. marina populations, such as Irish meadows, likely growing under temperature limitation, (ii) plants exposed to low irradiance are more vulnerable to thermal stress than plants exposed to high light, and (iii) plants incubated in different seasons display diverse thermal responses and distinct optimum of temperature for growth. Finally, (iv) we expect that warming will reduce PUFA contents and particularly omega-3, as high temperatures dismiss the requirements of unsaturation levels in photosynthetic structures.

Plant Collection
Plants were collected from a monospecific Zostera marina L. meadow situated in western Ireland at Finavarra (FV) (53.149437N, −9.133993W), southern Galway Bay (Figure 1). This area is classified as a Special Area of Conservation (SAC), Natural Heritage Area (NHA), and described as an "unpolluted" area by the Irish Environmental Protective Area (EPA) (NPWS, 2014) (more details in Beca-Carretero et al., 2019a. Mature apical shoots were manually harvested by SCUBA diving or snorkeling at a shallow depth of 2-3 meters at intervals of 5-10 meters along a 50 m transect line to prevent potential resampling of the same genotypes. At this depth range, the meadow is denser and more homogenous than at shallower or deeper areas (Beca-Carretero et al., 2019a). We selected shoots with similar weights (average of 2.7 ± 0.6 g DW shoot −1 ) and size (average of 32.3 ± 8.3 cm shoot −1 ) and the number of mature rhizome segments (3-4). Shoots were transferred to cooling tanks filled with ambient seawater and directly transported to the laboratory at the National University of Ireland Galway campus (Supplementary Figure 1) within 1 h. The collected samples were carefully cleaned, e.g., any remnants of sediment, epiphytic organisms or dead plant materials were removed from the shoots.
The experimental temperature range represented the natural annual range of sea surface temperature (SST) in western Ireland (∼5-17 • C) (Supplementary Figure 1); the 20 • C treatment represented the predicted increases of 2.6-3.1 • C of SST for the study region by the end of 2100 (IPCC, 2014;Hobday et al., 2016;Tinker and Howes, 2020); 24 • C temperature represented an elevation by +3-4 • C above maximum in situ SST at the collection site (∼ 19−20 • C, Stengel et al., 1999) that may be reached during a summer heat-wave events in the Atlantic Ocean (Feudale and Shukla, 2011). The duration of the temperature treatment was selected to simulate the duration of previous warming episodes in temperate regions (Hobday et al., 2016).
For each temperature treatment, two irradiance levels (180 µmol photons m −2 s −1 and 60 µmol photons m −2 s −1 ) were tested. Light levels were chosen based on in situ values recorded in September 2014 (maximum annual temperatures of 16 -17 • C) during (typical) cloudy days at the location where plants were collected for the experiment. It is noteworthy that in western Ireland 63% of days have some precipitation and more than 96% days have more than 20% cloud cover 1 . Therefore, we considered that using light levels recorded during cloudy days were ecologically more relevant for Irish coastal ecosystems than those during sunny days. Irradiances measured using a DIVING-PAM-II 2 along a 50 m transect long were 174.7 ± 42.4 µmol photons m −2 s −1 [n = 8]) at intermediate depths of Irish seagrass meadows (2-2.5 m), which is ∼70% higher than at deeper regions of the meadow (4-5 m; 57.6 ± 25.4 µmol photons m −2 s −1 [n = 8]) (Beca-Carretero et al., 2019b). Similar light levels were also previously used in previous Z. marina experiments (e.g., Evans et al., 1986;Höffle et al., 2011;Staehr and Borum, 2011;Beca-Carretero et al., 2018a).  Table 1).
Before starting the experiment, shoots were progressively acclimated from in situ temperature to the target experimental temperature (6, 12, 16, 20, and 24 • C) by gradually increasing or reducing chamber temperatures by 1 • C every 24-48 h (Figure 2). Also, plants were progressively acclimated to the target irradiance applied (60 or 180 µmol photons m −2 s −1 ) by increasing or reducing ∼30 µmol photons m −2 s −1 every 5-7 days. The experiment was started when all treatments reached their target temperatures and irradiance to ensure that all plants were pre-acclimated to lab-conditions for the same number of days (10-14 days).
For each temperature and irradiance treatment, three individual plants were incubated in each of the three individual transparent cylinders (n = 3) consisting of cylindrical perspex bottles (ID 10 cm, height 35 cm) with a volume of ∼3.5 L. Shoots were loosely tied to a weighed-down plastic net at the bottom of each container to maintain them in a vertical orientation. Air was supplied through pumps ensuring constant mixing and CO 2 supply. 80% of the seawater in the cylinders was replaced every two to 3 days. Irradiance levels (180 and 60 µmol photons m −2 s −1 ) were measured inside the cylinder at the upper position of the seagrass leaves; the light was provided by Lumilux cool daylight fluorescent lamps (OSRAM L18W/865, Germany) TABLE 1 | Geographic coordinates at Finavarra (FV), western Ireland, where Z. marina shoots were collected, water depth relative to the mean water level, in situ sea surface water temperature (SST) ( • C) measured at times when Z. marina biomass was collected, and averaged SST and daylength (h) (https://www.timeanddate.com) across 14 days at collection time (collection date 7 ± days).
Latitude (Beca-Carretero et al., 2018a). We continuously recorded water temperature using HOBO loggers (UA-002-64, Onset) installed in one cylinder per temperature and irradiance treatment. After 15-16 days, we measured a set of response parameters (see following).

Morphological Descriptors
Once the Z. marina shoots were acclimated to the target temperatures and before starting the experiment, all individual shoots were morphologically characterized. We measured: shoot fresh weight (g FW), total length (cm shoot −1 ), number of leaves (shoot −1 ) and the length of the rhizomes (cm). We standardized leaves and rhizome length across the all experiments to avoid potential impacts of morphology in plant responses by adjusting all shoots to two-three mature rhizome segments (2-3 cm total rhizome length) and 3-4 healthy leaves. At the end of each experimental treatment, the number of new leaves and new rhizome segments were counted, and rhizome and leaf elongation rates were assessed. To measure leaf growth (elongation) rates over time, two holes were pierced into leaves, one above the other with a distance of 2 mm above the basal meristem of the plant with a hypodermic needle (Sand-Jensen, 1975;Short and Duarte, 2001). We quantified the leaf formation rate (leaves shoot −1 d −1 ) by identifying the new number of leaves without punched holes, divided by the number of incubation days. Leaf elongation rate (cm shoot −1 d −1 ) was calculated as the length of new leaf material produced during the incubation time (in days), measured (i) from the base of the meristem to the punched holes, and (ii) the length of the newly produced leaves (Short and Duarte, 2001). The relative growth rate (RGR) was calculated according to: where Bi is the initial weight and Bf, the final weight of the shoot and t is the incubation period in days. We haphazardly selected a group of 10 Z. marina shoots with similar size at the time of collection from the same seagrass meadow to assess the ratio of initial fresh weight (FW): dry weight (DW). Moreover, we evaluated the shoot mortality rate by the following equation: Nf is the final shoot population (n = 9; 3 shoots per cylinder); Ni is the initial shoot population and t is the duration of the experiment (in days).

Biochemical Responses
After 15-16 days of incubation, only healthy tissues (avoiding epiphytes or damaged parts) of the youngest and second youngest leaves were selected for biochemical analysis. Selected biomass was cleaned with distilled water prior to processing. Samples were frozen at −20 • C and 48 h later freeze-dried. The samples were then kept at −20 • C. 24 h before conducting the experiments samples were again freeze-dried to remove the potential humidity.

Fatty Acid Analysis
We determined the fatty acid content and composition of Z. marina leaf biomass by applying the protocol previously used for macroalgae and seagrasses (Schmid et al., 2014;Beca-Carretero et al., 2018b. Fatty acid methyl esters (FAME) were obtained by direct transmethylation of ∼ 20-30 mg of powdered leaf biomass with dry methanol containing 2% (v/v) H 2 SO 4 . To prevent oxidation, vials were closed with nitrogen gas before being heated at 80 • C for 2 h under continuous stirring conditions. After transmethylation, we added 1 mL of Milli-Q water and later extracted the FAME using 0.5 mL of n-hexane. Analysis of FAME was conducted using an Agilent 7890A/5975C Gas Chromatograph/mass selective detector (GC/MSD) Series (Agilent Technologies, United States) equipped with a flame ionization detector and a fused silica capillary column (DB-WAXETR, 0.25 mm × 30 m × 0.25 µm, Agilent Technologies, Catalog No.: 122-7332). Identification of FAME was achieved by co-chromatography with authentic commercially available FAME standard of fish oil (Menhaden Oil,catalog no. 47116,Supelco). Total and individual fatty acid contents were quantified by comparison with a known quantity of added pentadecanoic acid 15:0 (99%, catalog no. A14664-09, Alfa Aesar, United Kingdom) as an internal standard. We added the 15:0 standard (10 µl, 5 mg mL −1 ) before starting the direct transmethylation and expressed the results as the mean values of 3 replicates (n = 3) for each treatment.
Fatty acids (FAs) content and composition were chosen as seagrass indicator due to their proved sensitivity to environmental fluctuations including temperature and irradiance alongside its importance as nutritional compound (e.g., Falcone et al., 2004;Sanina et al., 2008;Beca-Carretero et al., 2020).

Pigment Extraction
Chlorophylls and total carotenoids were determined following two consecutive extractions, using 5 mL of 80% acetone each time to extract pigments from ∼20-30 mg of powdered leaf biomass (Beca-Carretero et al., 2019a. The first extraction was conducted over 20 h, and the second over 4 h. To ensure optimal extraction and to avoid pigment oxidation, both extractions were performed in darkness at 4 • C and with continuous stirring, and the vials were closed under nitrogen gas. After the first extraction, samples were centrifuged (10,000 rpm) for 60 s, the supernatant (5 mL) kept in darkness at 4 • C, and the remaining biomass was used for the second extraction. The supernatants of both extractions were combined to a final volume of 10 mL, which was then used for pigment analysis (n = 3). We quantified chlorophyll a, b and carotenoids by spectrophotometric absorbance (CARY 50 Scan UV-Visible Spectrophotometer), following the equations of Lichtenthaler and Wellburn (1983).

Literature Review
A thorough literature survey evaluated previous seagrass studies on biochemical components with a potential nutritional value in response to controlled temperature exposure. Any data from field experiments or in situ observations were excluded. For the literature review selected compounds include lipids, fatty acids, proteins, carbohydrates, pigments and carbon and nitrogen composition among others. We included reports on necrosis as a parameter that may affect nutritional value. Table 2 contains details of leaf compounds which varied more than ± 5% when seagrasses were grown under predicted warming conditions compared to in situ summer temperature. In addition, Supplementary Tables 13-14 (Excel file) include all results regardless of the degree of variation reported.

Statistics
Data of morphological descriptors, growth rates, FA content and composition and photosynthetic pigments were Ln transformed, checked for homogeneity of variance using Bartlett's test and for normality applying Kolmogorov-Smirnov Frontiers in Marine Science | www.frontiersin.org Carbohydrates; chemical elements; fatty acids; proteins; pigments; other compounds; necrosis. % of variation was calculated according to the relative percentage of change from the in situ summer temperature and the predicted warming temperature of the temperature experiments, C. = carbon. Chl. a/b = Chlorophyll a/b; Phaeo. a/b = Phaeophytin a/b.
test. As data did not meet the criteria, the non-parametric test PERMANOVA analyses based on a similarity matrix created from the Euclidean distances were implemented. We performed two PERMANOVA models. The first model aimed at investigating differences in seagrass descriptors responses to month (April, July, November 2015 and January 2016), temperature (6, 12, 16, 20, and 24 • C) and irradiance levels (180 µmol photons m −2 s −1 vs. 60 µmol photons m −2 s −1 ). We used "temperature, " "irradiance, " and "month" as fixed factors. A pairwise test was applied to identify the treatments that differed significantly (p < 0.05). The second model was implemented to assess responses of seagrass traits to temperature and irradiance treatments each separate month. In this model "temperature" and "irradiance" were fixed factors. A pairwise test was applied to identify the treatments that differed significantly (p < 0.05). All data treatments and statistical analysis were performed using the PRIMERandPERMANOVA 6 statistical package. Growth rates, FAs and pigment concentrations are reported as means and standard deviation (SD). Pearson correlation analysis was used to assess potential relationships between FA composition (e.g., PUFA, SFA) and experimental temperatures (6,12,16,20, and 24 • C) including both irradiance levels (180 µmol photons m −2 s −1 vs. 60 µmol photons m −2 s −1 ).

Effect of Temperature on Plant Traits
Temperature treatment significantly affected growth of Z. marina plants in terms of leaf elongation, leaf formation, rhizome elongation, internode formation and RGR (PERMANOVA, p < 0.01). Overall, most growth descriptors displayed bellshaped patterns, with 16-20 • C representing more favorable temperatures, and lowest rates reached at coldest (6 • C) and warmest temperatures (24 • C) at both irradiance levels (Figure 3 and Supplementary Tables 1-5).
Most biochemical descriptors in Z. marina leaves including fatty acid content and composition and photosynthetic pigments were significantly affected by temperature treatments (Figure 4 and Supplementary Table 1). The average total fatty acid (TFA) content of Z. marina leaves was 1.41 ± 0.32% of DW, however, no clear pattern in response to temperature increase was observed. The most abundant FAs were polyunsaturated fatty acids (PUFAs), which accounted for an average of 57.45 ± 6.6% of TFA across all seasons (April, July, November and January). The most abundant PUFA was α-linolenic acid (ALA, 18:3 n-3), followed by linoleic acid (LA, 18:2 n−6) and hexadecatrienoic acid (HTA, 16:3 n−3). The second most abundant group of FAs were saturated fatty acids (SFAs) which accounted for an average value of 35.1 ± 6.6% of TFA, where palmitic acid (16:0) was the most common SFA, followed by stearic acid (18:0). Finally, monounsaturated fatty acids (MUFAs) accounted for an average of 4.3 ± 1.5% of TFA (Supplementary Tables 6-9).
Common to all seasonal experiments, the observed reductions in PUFAs following high temperature exposure were mostly attributed to changes in n−3 PUFAs (18:3 n−3 and 16:3 n−3), and occurred at both high and low light; lowest percentages were detected at 24 • C (Figure 5). For example, in plants incubated at higher irradiances, a significant reduction (Pearson correlation, p < 0.05, n = 30) of n−3 PUFAs from 6 to 24 • C was observed with a decrease of 0.9% by 1 • C increase. On the other hand, a significant increase in SFA and n−6 PUFA was observed with temperature ( Figure 5).

Effect of Low Light on Plant Traits
Leaf and rhizome growth rates of plants in low light was significantly reduced in comparison to plants exposed to high light (PERMANOVA, p < 0.01) (Figure 3 and Supplementary Table 2). For instance, plants acclimated at light-limited conditions produced new leaves and rhizome segments 21.1 and 18.7% slower than plant acclimated to high irradiances (PERMANOVA, p < 0.05). At temperatures above the optimum range (16-20 • C), light limitation induced a significant, abrupt decline in both leaf and rhizome growth compared to plants grown at high light. Additionally, at high temperatures, mortality rates were higher in plants exposed to low light than in plants incubated at high light; at low light, 30% of the incubated shoots died when grown at 24 • C (Supplementary Table 12).  Reduction in irradiance at optimal temperatures for growth (16-20 • C) generated an increase in FA unsaturation levels in leaves. The n−3/n−6 PUFA ratio in Z. marina exposed to low light significantly increased by 9.4% in comparison to plants exposed to high light (PERMANOVA, p < 0.05) (Figure 4 and Supplementary Table 2). Concentrations of photosynthetic pigments were also higher in low-light than in high lightacclimated plants (Figure 4 and Supplementary Table 10).

Effect of Seasonal Acclimation on Plant Responses to Temperature and Irradiance Treatments
At all seasons, temperature increases of 4 • C or more above in situ temperatures significantly enhanced growth of both rhizomes and leaf structures. For instance, in the experiment conducted in April, leaf growth at 16 • C was 44.9% higher than in plants exposed to the in situ temperature of 12 • C (PERMANOVA, p < 0.05, Figure 3). In summer, rhizome segments of plants incubated at 20 • C grew 58% faster than those of plants incubated at in situ temperatures (16 • C). Similar trends were observed in experiments conducted in January and November (PERMANOVA, p < 0.05, Supplementary Table 2).
Leaf and rhizome descriptors displayed different growth patterns depending on the seasonal acclimatization of the plants. Significantly higher leaf elongation and leaf formation rates were observed when the experiment was conducted in the coldest month, January (in situ temperatures of 5-6 • C), compared to April or July (PERMANOVA, p < 0.05, Figure 3 and Supplementary Table 2). By contrast, significantly higher rhizome elongation and internode formation rates were observed in plants incubated in warmer months (April and July) than in plants collected in January and November. For instance, at optimal temperatures for growth (16-20 • C), rhizome elongation of plants incubated in July (0.9 ± 0.01 cm d −1 ) was 71 and 99% higher than in plants incubated in January and November, respectively. Noticeably, both leaf and rhizome growth descriptors were significantly lower in November in comparison with plants incubated during other months (PERMANOVA, p < 0.05). At 24 • C, plants incubated in November and January displayed negative RGRs, but RGRs of plants incubated in April and July were positive. Mortality occurred in experiments conducted in April, November and January at 24 • C but all plants survived in all experimental treatments in July (Supplementary Table 11), and was more accentuated in plants incubated at low light in January (HL: 0.039 ± 0.02 dead shoot d −1 ; LL: 0.054 ± 0.02 dead shoot d −1 ). Interestingly, 22% of the plants incubated in January at 12 • C displayed the first stages of reproductive structures with the development of leaf sheaths corresponding to reproductive shoots (Figure 6). There was no evidence of reproductive structures in any other seasonal experiment.
TFA and photosynthetic pigment contents of Z. marina leaves differed significantly across seasons, however, there was no apparent clear pattern in response to temperature treatment (Figure 4) (PERMANOVA, p < 0.05). For instance, in plants incubated in April, the highest TFA contents were observed at lowest temperatures (6 and 12 • C) with average values of 2.05 ± 0.1% of DW, and lowest contents at highest experimental temperatures (24 • C, 1.3 ± 0.1% of DW) (PERMANOVA, p < 0.05). While in November, a contrary pattern with the largest TFA contents observed at 24 • C (1.7 ± 0.2% of DW) and lowest TFA were detected at lowest temperatures (12 • C, 1.3 ± 0.05% of DW) (Supplementary Table 10). Nearly all seasonal results revealed significant interactive effects between temperature, irradiance and time on growth descriptors (Supplementary Table 2) (PERMANOVA, p < 0.05). Particularly, growth responses to temperature treatments were differently modulated by light, depending on season. For example, in January, as temperature increased, low light significantly reduced RGR in comparison to high light; on the contrary, plants incubated in April and July where not affected by light treatment (PERMANOVA, p < 0.05). Additionally, rhizome elongation of plants incubated in April did not display negative responses to warming (20-24 • C) at either light treatment (Figure 3). However, a different pattern was observed in July, with low light significantly reducing rhizome elongation.

DISCUSSION
To the best of our knowledge, this is the first study attempting to assess the potential impacts of seasonality on outcomes of controlled exposure experiments using seagrasses. Our results indicate that seasonal acclimation significantly conditioned growth and biochemical responses of seagrass to thermal and low light stress. Besides, the reported outcomes pointed out that warming will likely favor growth and production of centrally distributed Z. marina populations which represent the broader habitat range of this species. These outcomes are of high relevance in predicting the effects of climate change on seagrasses and highlight the importance of seasonal conditioning on seagrass responses.

Effects of Temperature Increases on Growth Performance
Results confirm our hypothesis that predicted warming episodes of ∼3-4 • C above in situ current conditions (IPCC, 2014), including colder periods, in temperate regions may stimulate annual increases in growth of seagrass populations currently living under temperature-limited conditions (e.g., Krause-Jensen and Duarte, 2014).
In this study, regardless of light treatment, experimental optimum growth temperatures (between 16 and 20 • C) for Z. marina from western Ireland were lower than previously reported for other centrally located populations of the species (Lee et al., 2007;review in Beca-Carretero et al., 2018a). Such different latitudinal responses may be explained by the local adaptations of seagrass thermal tolerance to in situ temperature conditions (Short and Neckles, 1999;Bennett et al., 2015). In the context of climate mitigation potential of seagrasses in temperate regions, the expected increases in productivity can be interpreted as a likely higher capacity to uptake and store carbon in the sediments, therefore reinforcing their role as a carbon sink (Macreadie et al., 2019). Negative effects of thermal stress were observed in most of the seagrass descriptors at 24 • C. These results are in accordance with previous studies which correlated increases in respiration, growth inhibition and shoot death with anomalous high temperatures ranging from 20 to 30 • C (e.g., review in Lee et al., 2007;Collier and Waycott, 2014;Gao et al., 2017).
The observed warming stimulation in seagrass productivity under lower irradiances can favor their capacity to colonize deeper areas (more than 4-6 m in Irish coasts), as warming enhances photosynthetic efficiency of temperaturelimited seagrass populations (Krause-Jensen et al., 2012;Krause-Jensen and Duarte, 2014). In agreement with this, warm-adapted Atlantic, Pacific and Mediterranean Z. marina populations can colonize depths of 10-20 m (Morita et al., 2007;Ondiviela et al., 2014;Boutahar et al., 2020), while the vertical distribution of central and northern Z. marina populations is restricted to 4-10 m (e.g., Phillips, 1974;Krause-Jensen et al., 2011;Beca-Carretero et al., 2020). It is noteworthy that other factors such as greater light penetration can also explain the colonization of deeper areas within its southern distribution range such as in Mediterranean regions. In addition, both marine and freshwater macrophytes were reported to extend their vertical distribution to greater depths in response to increases in temperature (Duarte and Kalff, 1987;Rooney and Kalff, 2000;Jorda et al., 2020). As expected, plants growing at limiting irradiances had lower biomass production in both above and belowground tissues, probably explained by lower photosynthetic activity and carbon acquisition (e.g., Lee and Dunton, 1997;Ruiz and Romero, 2003).

Synergistic Effects of Thermal Stress and Light Reduction on Seagrass Traits and Survival
Results confirm that plants living at reduced irradiances or are experiencing episodes of high turbidity in the water column may be more vulnerable to warming episodes (e.g., Longstaff and Dennison, 1999;Krause-Jensen et al., 2005). In our study, at higher temperatures, growth was reduced at low light compared to high light, and there were more visible signs of leaf necrosis alongside higher mortality rates in most seasonal experiments. A recent study reported largescale declines of Z. marina populations after a combined effect of light reduction due to precipitation events alongside anomalous warming episodes (Johnson et al., 2021). The authors also found a rapid demographic recovery, mostly explained by the high abundance of the seedling. Studies in the Mediterranean Sea show that extreme summer temperatures provoked high mortality in deep-adapted (17-25 m) populations of Posidonia oceanica which were not able to recover in subsequent years, whereas shallow populations were less affected and recovered fully (Marbà and Duarte, 2010;Aoki et al., 2020). Additionally, projected sea-level rise will further reduce the irradiance available to Z. noltei plants distributed at their deeper range (∼2-2.8 m), resulting in more extreme conditions and physiological stress (Ondiviela et al., 2020).
Interestingly, proportions of n−3 PUFA relative to n−6 PUFA in plants receiving less light increased significantly, thus increasing the level of unsaturation of FA in leaves. Similar results were previously observed in the tropical seagrass H. stipulacea (Beca-Carretero et al., 2019a) and also in other marine primary producers (Guihéneuf and Stengel, 2013). This regulation of FA unsaturation levels in photosynthetic structures under light limitation may represent a photoadaptative mechanism enabling optimal photosynthetic and physiological performance (Gombos et al., 1994;Sanina et al., 2008;Solovchenko et al., 2008). Regarding pigmentation, plants exposed to low light accumulated higher contents of chlorophyll than those exposed to higher irradiances. This is a well-described mechanism of seagrasses to favor photosynthetic efficiency in low-light environments (e.g., Olesen et al., 2002;York et al., 2013).

Effects of Seasonal Modulation on Seagrass Responses to Temperature Increases and Irradiance Reduction
Experimental growth rates measured across seasons at respective in situ temperatures (6 • C in January, 12 • C in April and November and 16 • C in August) were highly similar to those measured in the field in shallow seagrass meadows in western Ireland (53.3272 • N,9.6168 • W) (Supplementary Table 12). These results confirm the robustness of our experimental conditions and resultant response measurements but also highlight the importance of capturing seasonal responses of temperate seagrasses to in situ conditions.
Unexpectedly, optimal conditions for growth (16-20 • C) were consistent across all seasons but critically, seasonality clearly modulated growth patterns of above and belowground compartments. On the one hand, higher rates of rhizome elongation and segment formation were observed in spring and summer when, in situ, plants were exposed to optimal environmental conditions for photosynthesis. This finding may be related to the seasonal capacity of temperate seagrasses to accumulate and store excess of energy produced from photosynthesis as carbohydrates in rhizomes in more suitable environmental conditions (e.g., Hemminga, 1998;Alcoverro et al., 2001). Such high-energy compounds can later support flower and seed production or act as an internal carbon source in less favorable environmental conditions in winter (e.g., Govers et al., 2015). On the other hand, the significantly higher leaf production observed in plants incubated in January-February might represent a mechanism to adjust photosynthetic structures to more suitable environmental conditions in spring. In contrast to this, plants incubated in November displayed the lowest leaf and rhizome growth rates which may represent a physiological acclimation to less suitable growth conditions in winter (e.g., Orth and Moore, 1986;Olesen and Sand-Jensen, 1993). There was evidence that seasonality modulated thermal and irradiance responses. For example, the growth performance of plants collected in April (in situ temperature at 12 • C) was less negatively impacted by extreme warming (24 • C) than plants collected in summer (in situ temperature 16 • C), particularly under low light conditions. Two further aspects suggest that seasonal acclimation modulates the observed experimental thermal responses: (i) highest mortality was observed in January (winter) while in July (summer) all plants survived; (ii) there were contrasting growth responses observed in plants incubated in April and November: in situ environmental conditions at these times of the year are very similar, however, growth rates in autumn were nearly half of those observed during experiments conducted in spring.
The flowering of Z. marina shoots observed in January at 12 • C suggests that warming in cold seasons may advance flowering events in temperate seagrasses. Recent field studies reported that flowering in Irish Z. marina populations in situ currently starts in early March-April (11-12 • C) and ends in October-November (Beca-Carretero et al., 2019b). Previous studies also reported that anomalous warming can alter reproductive processes by increasing seed production in seagrasses Qin et al., 2020a). Overall, year-round effects of climate change can lead to an internal decoupling of physiological and biochemical adjustments in marine primary producers (Bokhorst et al., 2009). For example, the accumulation of energy reserves such as non-structural carbohydrates may be affected, limiting their accumulation and impairing plant survival in cold winter periods (Alcoverro et al., 2001). Similarly, changes in temperature can affect flower production, flower fertilization and consequently seed production and germination (Qin et al., 2020a,b).
Detailed experimental studies that have investigated environmental stress on plants across different seasons are scarce and are mostly limited to terrestrial systems. For example, De Boeck et al. (2011) reported different sensitivities and responses to extreme events in terrestrial plants, dependent on the seasonal timing of the experiment. Our findings clearly demonstrate that, for temperate seagrasses, the time of year when temperature exposure occurs contributes to the observed growth, biochemical and survival responses. These results have relevant implications for future experimental design and interpretation of climate change effects on seagrasses. Nevertheless, caution is needed when extrapolating laboratory-derived results to natural habitats due to the complexity and the synergic interaction of several environmental parameters co-occurring in natural environments (De Boeck et al., 2010). Clearly, additional investigations are necessary to capture seasonal responses more comprehensively.

Effects of Experimental Warming on the Nutritional Value of Seagrasses Leaves
Particular nutritionally valuable compounds in seagrasses comprise essential fatty acids, carbohydrates, proteins or pigments, alongside secondary metabolites like phenolic compounds (e.g., Wang and Zheng, 2001;Pradheeba et al., 2011;Ahmed and Stepp, 2016).
The detailed review of the literature, together with experimental evidence from the present study, reveal sometimes contradicting trends pointing toward differential responses of specific biochemical compounds to increasing temperature. For instance, in some studies (e.g., Marín-Guirao et al., 2018;Artika et al., 2020), carbohydrates and photosynthetic pigments significantly increased in seagrass leaves at elevated temperatures; however, other studies report a significant decrease (e.g., Egea et al., 2018). Similar contrasting observations are documented for proteins, carbon, nitrogen and phosphorous contents (e.g., Hernán et al., 2017;Nguyen et al., 2020a;Viana et al., 2020; Table 2 and Supplementary Tables 14, 15). These differential patterns are likely to be linked to individual species-specific thermal tolerances, different degrees and duration of warming experiments or experimental time of year; they highlight an important gap in research into the effect of warming on the nutritional value of seagrasses. Common to all studies was the increase in necrotic leaf area under increased temperatures.
Our experimental results clearly demonstrate that warming significantly reduces omega-3 (n−3) PUFA levels, and favors the accumulation of SFA (see also Beca-Carretero et al., 2018b, as previously reported for other marine primary producers such as micro-and macroalgae (e.g., Gladyshev et al., 2013;Aussant et al., 2018). PUFAs play a key role in regulating the membrane fluidity, electron transport and photosynthetic activity in the thylakoids of the chloroplast; hence, the reported lower levels of FA saturation in photosynthetic structures with increases in temperature could be expected (e.g., Sanina et al., 2008;Gosch et al., 2015;Sijil et al., 2019). Changes in FA levels and composition have previously been implemented as a sensitive indicator of responses of seagrasses, saltmarshes and algae to different environmental settings including temperature, irradiance or salinity (Duarte et al., 2018;Kostetsky et al., 2018;Beca-Carretero et al., 2020).
Climate change is expected to increase metabolic demands of seagrass herbivores (Burnell et al., 2013) and thus affect their nutritional preferences (de los Santos et al., 2012;Hernán et al., 2016;Pagès et al., 2018). Recent studies have demonstrated that changes in the biochemical and nutritional composition of terrestrial and marine primary producers induced by climate change effects are impacting trophic interactions between plants and consumers and their feeding behavior (O'Connor, 2009;Prado et al., 2010;Hoover et al., 2012;Hernán et al., 2016). Herbivore feeding experiments should be undertaken to evaluate the potential effects of diets with reduced n−3 PUFA contents, simulating future warming conditions.

Ecological Significance
Under global warming scenarios marine heatwaves are predicted to increase in duration, intensity and frequency but also can occur across different seasons including in colder periods (Belkin, 2009;Frölicher et al., 2018). Our experimental results are thus ecologically highly relevant as they indicate that warming will not only affect seagrass productivity and physiology but also their reproductive effort. Conventionally, climate experiments and particularly warming assessment (review in Beca-Carretero et al., 2018b), have been conducted in warmer months, however, our study highlighted the importance of conducting climate change experiments across all seasons. In addition, our outcomes showed the higher vulnerability of plants exposed to low levels of irradiance under warming episodes. This is particularly important in the context of expected increases and variability of extreme events such as storms and extreme warming events associated with global change, resulting in synergetic stresses in shallow marine habitats (Lefcheck et al., 2017). Nevertheless, our experiments only represent short-time simulations such as marine heatwaves, and complementary experiments simulating long-term effects of warming and light stress are necessary to complement such an approach. Additionally, the magnitude of the selected range of temperatures, irradiance treatments and the duration of the exposure would significantly affect seagrass responses to experimental warming (Nguyen et al., 2021). Overall, our results have important ecological implications for the adequate management of seagrass ecosystems to ensure a balance of optimal temperature and irradiance conditions; they also point toward the critical importance to develop all-year monitoring programs in order to capture more comprehensively the potential effects of climate change events across all seasons.
Lastly, our results demonstrated that warming will dismiss the capacity of temperate seagrasses to produce and accumulate the essential omega-3 PUFA. These observations can have critical biological consequences as omega-3 are essential FAs that cannot be synthetized by herbivores and are transferred throughout the trophic chain (Sargent et al., 1999;Tocher, 2003;Parrish, 2009). Besides the direct ecological impact on herbivores, a reduction in omega-3 PUFA in seagrasses in response to climate change may also reduce the nutritional quality of some marine food sources (fishes or crustaceans) for human consumption (e.g., Kang, 2011;Hixson and Arts, 2016).

CONCLUSION
Our results demonstrate that responses of temperate seagrasses to experimental temperature and light stress are partially modulated by the seasonal acclimatization during an annual cycle: the time of year at which experiments are conducted strongly influences the observed growth, physiological and survival responses. These finding have highly relevant implications for experimental design, and in the understanding and interpretation of the effects of experimental warming on seagrasses. Increases in sea surface temperature (SST) predicted for the end of this century are expected to favor the growth and production of temperate seagrasses growing under temperature-limited conditions; by contrast, light-limited plants will be more vulnerable to thermal stress. Additionally, warming will affect the biochemical composition and associated nutritional quality of seagrass leaves with a potential effect on herbivores and its associated trophic change.

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/s.

AUTHOR CONTRIBUTIONS
PB-C and DS developed the experimental design. PB-C, MJ-M, CS, and FG collected the samples and conducted the laboratory experiments at NUIG. TA-G, FG, CS, and PB-C performed the laboratory analysis. PB-C performed the statistical analysis, prepared a first draft of the manuscript. DS provided resources and funding for the project. All co-authors commented on and provided edits to finalize the original manuscript. All coauthors are in agreement with the submission of the final manuscript. All authors contributed to the article and approved the submitted version.