Rising Temperature Is a More Important Driver Than Increasing Carbon Dioxide Concentrations in the Trait Responses of Enhalus acoroides Seedlings

: Increasing temperature and CO 2 concentration are among the most important factors affecting marine ecosystems under climate change. We investigated the morphological, biochemical, and physiological trait responses of seedlings of the tropical seagrass Enhalus acoroides under experimental conditions. Trait responses were greater under temperature effects than increasing CO 2 concentration. Seedlings under rising temperatures showed enhanced leaf growth, lower leaf nutrient content, and stimulated down-regulating mechanisms in terms of photo-physiology. Increasing CO 2 concentrations did not show any signiﬁcant effects independently. There was a signiﬁcant interaction for some of the trait responses considered, such as leaf number and carbon content in the roots, and trends of higher starch concentrations in the leaves and lower rETR max under combined enriched CO 2 and high temperature, even though none of these interactions were synergistic. Understanding the single and interactive trait responses of seagrass seedlings to increasing temperature and CO 2 concentration is of importance to determine the relative responses of early life stages of seagrasses, which may differ from adult plants, in order to form a more holistic view of seagrass ecosystem health under climate change.


Introduction
Seagrasses are important ecosystem engineers that are globally distributed and provide important ecosystem services [1]. They are one of the most threatened marine ecosystems due to their presence in shallow coastal areas resulting in accelerated seagrass biomass loss in recent decades [2]. While there is some space for optimism in temperate and subtropical areas [3][4][5], regional losses can still be high in the tropics, such as in Indonesia, where a 75% loss has been reported in the last five years [6].
Before seagrass loss occurs, changes in environmental conditions can be evaluated by seagrass trait responses, which can be used as indicators of environmental change in coastal management [7,8]. They have been described as valuable bioindicators since they integrate environmental impacts over measurable and definable timescales [9]. The relative sensitivity of seagrass trait responses to specific stressors or the suitability of different species is unknown. However, its understanding is critical in the selection of seagrass indicators for the assessment of different stressors [8]. As a solution, laboratory experiments have helped to isolate the responses of seagrass traits to single and multiple stressors [10][11][12].
Single effects of climate change have been widely studied in adult seagrasses (see Reference [13] for a review about this topic) although less information is available for seedlings, especially for tropical species. In seedlings, temperature responses varied from increased mortality due to inhibition of the photosynthetic system [14,15] or an increase on germination speed [16], and sublethal effects on photo-physiological performance and leaf or root development [17] to positive effects in terms of enhanced above-ground biomass [18,19]. The previously mentioned studies have also shown that the effects of temperature on seedlings vary depending on the age of the seedling, the temperature tested, or the length of the exposure. In contrast, CO 2 enrichment increases the rate of the dark reaction in P/S and, thus, growth and storage capacity, increasing carbon reserves in the below-ground tissues and lowering nitrogen content in leaves, as shown in Posidonia oceanica seedlings [20].
In the context of global climate change, seagrasses are rarely affected by just one stressor [11][12][13]. Species-specific trait responses to interactive effects of climate change stressors have been examined on adult seagrass plants [10,[21][22][23][24][25], but little information is available on seagrass early life stages, and even less in tropical species. From these studies on adult seagrasses, we know that it is possible that CO 2 enrichment could ameliorate or enhance the effects of warming in seagrass seedlings, but more studies are needed to provide a better understanding.
Enhalus acoroides is a common tropical Indo-Pacific foundation seagrass species widely distributed within the region. It is found from shallow intertidal to subtidal areas, being especially exposed to local and global disturbances, including climate change. E. acoroides is the biggest seagrass species in terms of its architectural properties [26] with shoots up to 1 m long [27]. It shows both horizontal clonal growth and production of new plants from seeds that are genetically unique individuals (genets). Sexual reproduction is a major sink of resources in E. acoroides compared to other co-habiting species [28,29], and seeds are large and full of carbon reserves [18]. The photosynthetic process in E. acoroides potentially has a very unique carbon concentrating mechanism more dependent on the uptake of CO 2 rather than HCO 3 [30] and is potentially C i limited at current concentrations [31]. As large seagrasses show low growth rates and biomass turnover [32], environmental changes could be detected more quickly in the seedling phase. Moreover, seedling development and growth are the most critical steps in seagrass survival [33] and could be more sensitive to changes than adult individuals from the same species. However, there is a lack of research conducted on responses of E. acoroides to environmental changes, especially on seeds and seedlings [18,29,30,[34][35][36][37][38][39]. The understanding of seedling trait responses under different environmental conditions may help to understand the effect of those changes, but also enhance future restoration plans of these threatened ecosystems. The conservation and management of these ecosystems requires knowledge about how they perform at the individual plant level in both early and adult life stages.
This study aims to measure the morphological, biochemical, and physiological trait responses of E. acoroides in the seedling phase under an increased average temperature and increased CO 2 concentrations. We conducted a laboratory experiment using seedlings of E. acoroides to test the hypothesis that seagrass seedling development would be enhanced (synergistic effect) by independent and combined increases in temperature and CO 2 concentrations in terms of growth and photosynthetic performance. The results of this study will help to understand the success of seedlings of E. acoroides under extended exposure to high temperatures and increasing CO 2 concentrations, which are predicted to occur under global climate change scenarios.

Collection and Maintenance of Seagrass Seedlings
Fruits of E. acoroides were collected on 12 January, 2017 on the southwest side of Barrang Lompo Island, South Sulawesi, Indonesia (S 5 • 03 05, E 119 • 19 37) where E. acoroides is abundant in a depth range of 1-3 m. Annual seawater temperature range in the dry season varies in this area between 26 and 32 • C [40], and light intensity ranges between 100 and 250 µmol photons m −2 s −1 [41]. The ripe seagrass fruits were opened, packed in a Styrofoam box with wet cotton, and then transported to the Marine Experimental facilities (MAREE) at the Leibniz Centre for Tropical Marine Research (ZMT) in Bremen (Germany) in less than 24 h. Once at the MAREE, seeds were planted directly in polypropylene trays previously filled with silicate sediment of at least 10 cm in depth.
The trays were kept in 250-l aquaria filled with artificial seawater (Red Sea Salt, Red Sea Europe Company, 27130 Verneuil d'Avre et d'Iton, France) under controlled fluorescent light (200 ± 30 µmol photons m −2 s −1 , 12:12 h light:dark), temperature (26 ± 1 • C), and salinity (35) conditions. Seagrass seeds were maintained under these controlled conditions for approximately 8 weeks until the experiment began. During this time, seeds developed into seedlings with both short leaves and roots (see below for average size).

Experimental Design and Setup
We conducted a fully crossed two-factorial experiment ( Figure 1) to study the effects of temperature, CO 2 , and their interaction on the trait responses of E. acoroides seedlings. We ran the experiment using an indoor flow-through system with 24 aquaria (10 × 29 × 50 cm) in a temperature-controlled room set at 26 • C (MAREE, ZMT, Bremen, Germany).

Collection and Maintenance of Seagrass Seedlings
Fruits of E. acoroides were collected on 12 January, 2017 on the southwest side of Barrang Lompo Island, South Sulawesi, Indonesia (S 5°03′05, E 119°19′37) where E. acoroides is abundant in a depth range of 1-3 m. Annual seawater temperature range in the dry season varies in this area between 26 and 32 °C [40], and light intensity ranges between 100 and 250 µmol photons m −2 s −1 [41]. The ripe seagrass fruits were opened, packed in a Styrofoam box with wet cotton, and then transported to the Marine Experimental facilities (MAREE) at the Leibniz Centre for Tropical Marine Research (ZMT) in Bremen (Germany) in less than 24 h. Once at the MAREE, seeds were planted directly in polypropylene trays previously filled with silicate sediment of at least 10 cm in depth.
The trays were kept in 250-l aquaria filled with artificial seawater ( (35) conditions. Seagrass seeds were maintained under these controlled conditions for approximately 8 weeks until the experiment began. During this time, seeds developed into seedlings with both short leaves and roots (see below for average size).

Experimental Design and Setup
We conducted a fully crossed two-factorial experiment ( Figure 1) to study the effects of temperature, CO2, and their interaction on the trait responses of E. acoroides seedlings. We ran the experiment using an indoor flow-through system with 24 aquaria (10 × 29 × 50 cm) in a temperature-controlled room set at 26 °C (MAREE, ZMT, Bremen, Germany). Twenty-four small clear glass aquariums (striped box) were placed inside the water baths under the selected temperatures and were independently supplied with seawater either enriched with CO2 or at ambient CO2 using a flow-through system attached to a peristaltic pump (black oval). Water was pre-conditioned with CO2 from water reservoirs (striped or dotted circles) of the desired CO2 concentrations. In addition, CO2 bubbling was also added directly in the small aquariums to help maintain CO2 treatments within the ETs.
In each aquarium, three seagrass seedlings were carefully planted by hand. Each seedling had leaves and roots fully developed with an average (±standard error, SE) of 1.93 ± 0.06 cm in leaf length, 6.27 ± 0.22 cm of root length, 1.10 ± 0.02 cm of seed diameter, and 0.9-5.5 g fresh weight (FW). Each aquarium was previously filled with pre-cleaned silicate sediments (~7 cm depth) and 10 L of natural seawater (NSW, North Sea). The NSW was pre-filtered before storage to remove plankton and avoid organic matter reminerali- Figure 1. Experimental setup in this study. Experimental tanks (ETs) were used as a water bath (white box) set at either regimes of 27 or 31 • C. Twenty-four small clear glass aquariums (striped box) were placed inside the water baths under the selected temperatures and were independently supplied with seawater either enriched with CO 2 or at ambient CO 2 using a flow-through system attached to a peristaltic pump (black oval). Water was pre-conditioned with CO 2 from water reservoirs (striped or dotted circles) of the desired CO 2 concentrations. In addition, CO 2 bubbling was also added directly in the small aquariums to help maintain CO 2 treatments within the ETs.
In each aquarium, three seagrass seedlings were carefully planted by hand. Each seedling had leaves and roots fully developed with an average (±standard error, SE) of 1.93 ± 0.06 cm in leaf length, 6.27 ± 0.22 cm of root length, 1.10 ± 0.02 cm of seed diameter, and 0.9-5.5 g fresh weight (FW). Each aquarium was previously filled with precleaned silicate sediments (~7 cm depth) and 10 L of natural seawater (NSW, North Sea). The NSW was pre-filtered before storage to remove plankton and avoid organic matter remineralization. A pump was placed in each aquarium to ensure the water aeration and homogeneous mixing by moving water from the bottom to the top, reducing the diffusive boundary layer.
We applied two levels of temperature: low and high (27 and 31 • C) and two levels of CO 2 concentrations (ambient, 400 µatm, and enriched, 700 µatm) and all the possible combinations, yielding in four different treatments: low temperature and ambient CO 2 concentrations, low temperature and enriched CO 2 concentrations, high temperature and ambient CO 2 concentrations, and high temperature and enriched CO 2 concentrations. Each of the 24 aquaria were independent and served as replicates (n = 6) of the four treatments. Experimental temperature levels were selected based on lower and higher average temperatures in seagrass meadows in the region. The levels of CO 2 concentrations were selected according to scenarios forecasted by the Intergovernmental Panel on Climate Change (IPCC). Therefore, the selection of experimental temperature and CO 2 concentration levels was based on obtaining realistic future trait responses rather than to obtain stress-related responses.
The experimental temperatures were obtained by placing aquaria in eight experimental tanks (ET) (42 × 87 × 58 cm) of 250 L (three aquaria in each ET), which served as water baths maintaining a constant water temperature. The water temperature was continuously controlled in each ET using heaters (EHEIM GmbH, Germany) connected to a digital controller and probe (±0.2 • C). A pump was also placed in each ET to ensure that the water was mixed to keep the temperature constant.
Each of the two CO 2 concentrations, termed ambient and enriched, were applied to two 115-l water reservoirs. We achieved the desired concentration in each water reservoir by constantly bubbling NSW with CO 2 from the gas cylinder's main tank. Each water reservoir was used to provide an equal and constant water flow to aquaria at a rate of~4 mL min −1 by using a 24-channel peristaltic pump (ISMATEC, Cole-Parmer GmbH, Germany) ensuring total water renovation inside each aquarium every~1.5 days. Therefore, each water reservoir pumped NSW to half of the total aquaria. Water reservoirs were manually emptied from any remaining water and refilled with fresh NSW every other day. A stock nutrient solution was added and gently mixed with NSW in each water reservoir during refill to achieve a final concentration of 1 µM NH 4 NO 3 and 0.1 µM KH 2 PO 4 (Merck, Germany) to avoid nutrient limitation during the experiment. The desired CO 2 concentrations (ambient or enriched) were maintained by also directly bubbling each aquarium with CO 2 from the gas cylinder's main tank. The loss of CO 2 in the aquaria was avoided by covering each ET with a transparent PVC lid that also reduced water evaporation. Water constantly overflowed from the aquaria to the water bath of the ETs, ensuring water renewal. At the same time, ETs were drained of the surplus water flowing out of the aquaria. The light was provided with two LED (light emitting diode) lamps (Hydra Fifty-two HD, AquaIllumination ® , Iowa) placed at the same height at the top part of each ET, providing the aquaria with a light intensity of 200 ± 20 µmol photons m −2 s −1 set on a 12:12 h light:dark cycle with sunrise and sunset simulation.
Seedlings were acclimated in the experimental setup at 27 • C and no CO 2 additions for five days. During this phase, all ETs were controlled by the same circulation system. Therefore, they were all connected to ensure the same initial conditions for all seedlings. Physico-chemical parameters were regularly controlled to ensure that the desired conditions were maintained. After this acclimation phase, ETs were separated according to the experimental treatments described above. The temperature of four random ETs (i.e., 12 aquaria) was increased by 1 • C per day from 27 to 31 • C, while the other four ETs were maintained at the same temperature as during the acclimation phase. Once the temperature of each aquarium was stable, CO 2 enrichment began and the experimental treatments lasted 35 days (20 March 2017 to 24 April 2017).

Water Sampling
Water parameters, including pH, dissolved oxygen, temperature, conductivity, and salinity, were measured in each aquarium every three days using a multi-parameter probe Appl. Sci. 2021, 11, 2730 5 of 18 (WTW Multiprobe). Water samples for silicate and phosphate were taken every week from random aquaria of each treatment (n = 2) and from both water reservoirs with a syringe, immediately filtered (0.45 µm pore size) in pre-rinsed polyethylene bottles, and frozen (−20 • C). Analysis was performed using a continuous flow injection analyzing system (Skalar SAN ++ -System) following standard methods for nutrient analysis [42].
For alkalinity measurements, water was sampled directly from the aquaria and immediately filtered (0.45 µm in pore size) and stored in the fridge (4 • C) in the dark. Samples were analyzed within two days, and mercury chloride was added to samples that were not analyzed within this period. Alkalinity was photometrically measured with Metrohm Dosimat ® using an Open Vessel Titration [43]. pCO 2 values were estimated based on salinity, conductivity, alkalinity, pH, and dissolved phosphate and silicate concentrations [44].

Morphological Traits
At the end of the experiment, seedlings were removed from the aquaria and separated with a glass spatula into three parts (leaves, seeds, and roots) and were weighed separately to get a biomass of each plant tissue. The number, maximum length, width, and surface area (SA) of the leaves were measured. Leaf length measurements were taken from the leaf base to the tip. The leaf SA was calculated by multiplying the leaf length by the width. In the seeds, the diameter and height were estimated while, in the roots, the number and maximum length were measured.

Biochemical Traits
The nutrient content was analyzed in the pooled material of the three seedlings from each aquarium. The nitrogen and carbon content (%N and %C, respectively) were analyzed in leaves, seeds, and roots that were previously dried (60 • C for 48 h) and ground to a fine powder with mortar and pestle. Aliquots of the ground sample were weighed into tin capsules using an analytical scale prior to analysis (Euro EA3000 Elemental Analyzer).
The concentrations of soluble sugars (sucrose) and starch were measured on leaf and seed material that was frozen (−80 • C), subsequently freeze-dried (48 h), and ground to a fine powder with mortar and pestle. Sucrose was extracted from plant tissue by heating (80 • C) the sample in EtOH 95%. The ethanol extracts were then evaporated by bubbling the samples with N 2 , and the remaining residues were dissolved in deionized water. Starch was extracted after having the ethanol-insoluble residue in 0.1 N NaOH for 24 h. The sucrose and starch concentrations were determined spectrophotometrically (486 and 640 nm, respectively) using a F200-Pro TECAN© plate reader. Resorcinol and anthrone assays were respectively used for sucrose and starch determination, while sucrose was used as the standard for calibration curves [45,46]. Results were reported in sucrose equivalents g dry weight (DW). Current testing of this method has shown that NaOH extracts both starch and cellulose, which can confound the results. Regarding the sucrose determination, this method only determines ketoses (as fructose), so we are ignoring glucose, which is the other component of sucrose, which may lead to an underestimation of the final concentrations (M. Birkicht, personal communication).

Physiological Traits
Seedling growth rate was measured using the leaf marking method [47]. At the beginning of the experiment, leaves were perforated close to the seed using a needle. At the end of the experiment, the elongation of the leaf, known as the length from the base of the leaf to the mark, was measured. Leaf SA growth rates were obtained by multiplying the leaf elongation by the leaf width and dividing the result by the number of days since the seagrass leaves were marked. Leaves without holes were considered new leaves and were measured from the leaf base to the tip of the leaves.
Responses of photosynthetic performance were measured at the end of the experiment using a PAM-2500 chlorophyll fluorometer (Walz, Germany) with rapid light curves (RLC).
The PAM optical cable was attached to the mid part of the second leaf of the seedling by using a leaf clip, and~3 mm in distance from the tissue. This position was chosen since it represents similar distances from the water surface (and, thus, from the light source) among plants with different leaf lengths [48]. Before measurements, the leaves were adapted in the dark for 5 min while maintained in NSW. The RLC consisted of 12 saturating light pulses (separated by 30-s intervals) increasing photosynthetic active radiation (PAR) between pulses until 3001 µmol photons m −2 s s −1 . From the RLC data, several parameters were calculated. The maximum quantum yield (Fv/Fm), the relative maximum electron transport rate (rETR max ), the light saturation coefficient (Ek), and the slope of the light limited part of the curve (α) calculated using Phytotools [49] under R software [50], following the model in Reference [51]. Efficient utilization of Fv/Fm was calculated by following Equation (1) [52].
Maximum quantum yield: where Fm is the maximum dark-adapted fluorescence and F o is the minimum dark-adapted fluorescence. The relative electron transport rate (rETR) was calculated for each step of the curve following Equation (2) [53].
where F m is the light-adapted maximum fluorescence, F is the fluorescence yield at a particular light level, and PAR is the photosynthetic active radiation. The maximum rETR value (rETR max ) was calculated as the inflection point of the fitted rETR curve.

Statistical Analysis
Average values of the three seedlings within each aquarium were calculated for physiological and morphological trait responses with no calculation of the deviation of the data. Descriptive statistics as mean values and standard error (SE) were estimated within treatments (n = 6). Prior to statistical analysis, the variance of the residuals was checked for normality (Shapiro-Wilk test) and homoscedasticity (Levene test). A two-way analysis of variance (ANOVA) was performed to detect the influence of temperature (two levels: 26 and 31 • C), CO 2 (two levels: ambient and enriched), and the interaction between these two factors on E. acoroides seedling performance. The statistical tests were carried out with SPSS (Statistical Package for the Social Sciences, IBM SPSS Statistics for Windows v.24, Armonk, NW, USA). All comparisons were considered significant at p-values < 0.05.

Water Parameters
Salinity and conductivity were constant across all experimental treatments (Table 1). Within each defined treatment, the temperature was nearly constant. The pCO 2 concentrations in the two water reservoirs were~485 and~972 µatm for the ambient and enriched treatments, respectively. In aquaria, carbonate system parameters of the ambient CO 2 treatments remained well within the target concentration of 400 µatm (calculated pCO 2 379 µatm). On the contrary, the concentrations in aquaria under enriched CO 2 treatments were lower than the target concentration (700 µatm) and the concentration in the water reservoir, but higher than the ambient treatments, with mean values of 464 and 624 µatm in the 27 • C-enriched CO 2 and 31 • C-enriched CO 2 treatments, respectively (Table 1). Among all treatments, pCO 2 and alkalinity values were the highest in the 31 • Cenriched CO 2 treatment. This is in line with the low pH observed in this treatment, showing that the CO 2 enrichment reduced the pH of the water. Dissolved oxygen was lower in the regimes of 31 • C when compared to the regimes of 27 • C. Phosphate and silicate concentrations remained very low throughout the experiment and were similar in all treatments during all the experimental periods.

Morphological Traits
The responses of leaves to the temperature and CO 2 treatments differed across morphological traits (Figure 2 and Table 2). The highest leaf number was found in the 31 • Cenriched CO 2 treatment with a mean of five leaves per seedling, showing a significant interaction effect (Figure 2a and Table 2).    2 and their interaction on the morphological leaf traits ( Figure 2) and growth rate (Figure 3).  The photosynthetic performance of the leaves was significantly higher under regimes of 27 °C (Figure 4), as indicated by higher rETRmax values calculated from the rapid light curves (Tables 3 and 4). In comparison, the maximum leaf length, leaf SA, and leaf weight values were significantly higher at a high temperature (31 • C) regardless of the CO 2 treatment (Figure 2b-d and Table 2). Even though CO 2 treatments did not show a significant effect in most morphological traits, some trends were observed in terms of leaf length, with longer leaves under ambient CO 2 compared to enriched CO 2 treatments (p-value = 0.065, Table 2). The leaf width did not show any significant difference among treatments (Figure 2e and Table 2). For the seeds and roots, there were no differences among any of the treatments ( Figure S1 and Table S1).

Physiological Traits
Leaf growth rates were significantly higher under the regimes of 31 • C while there was no influence of CO 2 treatments (Figure 3, Table 2). The photosynthetic performance of the leaves was significantly higher under regimes of 27 • C (Figure 4), as indicated by higher rETR max values calculated from the rapid light curves (Tables 3 and 4). The photosynthetic performance of the leaves was significantly higher under regimes of 27 °C (Figure 4), as indicated by higher rETRmax values calculated from the rapid light curves (Tables 3 and 4).   Additionally, there was an interaction (p-value = 0.078, Table 4) between the two factors, which are temperature and CO2 concentration, in rETRmax and Ek parameters, with the highest values observed under the 27 °C-enriched CO2 treatment (Table 3). In contrast,   Additionally, there was an interaction (p-value = 0.078, Table 4) between the two factors, which are temperature and CO 2 concentration, in rETR max and Ek parameters, with the highest values observed under the 27 • C-enriched CO 2 treatment (Table 3). In contrast, the 31 • C-enriched CO 2 treatment showed the lowest rETR max and Ek values (Table 3). Fv/Fm and α measures of the photosynthetic maximum quantum yield and efficiency did not show significant differences among the four treatments, nor any interaction between temperature regimes or CO 2 concentrations (Tables 3 and 4).

Biochemical Traits
Biochemical traits, namely carbon, nitrogen, sucrose, and starch content, showed different responses in leaves, seeds, and root tissues in the four experimental treatments considered in the study (Figures 5 and 6, Tables 5 and 6). With increasing temperatures, the percentage of carbon and nitrogen significantly decreased within the leaves (Figure 5a,b and Table 5).  nitrogen content, sucrose concentrations significantly responded to an increasing ature on seed tissues, which are less used in regimes of 31 °C (Figure 6c and Tab highest concentrations of sucrose in the seeds were found in the 31 °C-ambient C ment, whereas the highest starch concentration in the seeds was found in the 3 riched CO2 treatment, even though differences were not significant (Figure 6c,d a  6).  Table 6. Results of analysis of variance (two-way ANOVA test) of the effects of temperatu CO2 on the sucrose and starch content of the leaves and seeds ( Figure 6). p-values in bold tistically significant (<0.05) (TSS: Total sum of squares).

Temperature as the Main Driver of Seedling's Trait Responses
Temperature was the dominant driver of E. acoroides seedling's trait respons this experiment, and there were no synergistic effects in the interaction of both Figure 6. The effect of increasing temperature and CO 2 enrichment on sucrose (left column) and starch (right column) content (sucrose equivalents g −1 DW) on leaves (a,b) and seeds (c,d) in seedlings of Enhalus acoroides at the end of the experimental period. Values are mean ± SE (n = 3 for leaves and n = 6 for seeds). Table 5. Results of analysis of variance (two-way ANOVA test) on the effects of temperature and CO 2 on the carbon and nitrogen content of the leaf, seed, and root tissue ( Figure 5). p-Values in bold are statistically significant (<0.05) (TSS: Total sum of squares).  Table 6. Results of analysis of variance (two-way ANOVA test) of the effects of temperature and CO 2 on the sucrose and starch content of the leaves and seeds ( Figure 6). p-Values in bold are statistically significant (<0.05) (TSS: Total sum of squares). In contrast, there were no significant differences in single factor effects of temperature or CO 2 in the carbon and nitrogen content within the seeds (Figure 5c,d and Table 5) or in the roots (Figure 5e,f and Table 5). However, there was an interactive effect between temperature regimes and CO 2 concentration in the carbon content in the root tissue, which showed the highest concentrations in the 31 • C-ambient CO 2 and 27 • C-enriched CO 2 treatments (Figure 5e and Table 5).

Dependent
Leaf sucrose and starch concentrations were generally higher in enriched CO 2 treatments (Figure 6a,b and Table 6). However, these differences were not statistically significant. There was an interaction effect on starch concentrations in the leaves (p-value = 0.07), which was especially high in the 31 • C-enriched CO 2 treatment. In contrast to carbon or nitrogen content, sucrose concentrations significantly responded to an increasing temperature on seed tissues, which are less used in regimes of 31 • C (Figure 6c and Table 6). The highest concentrations of sucrose in the seeds were found in the 31 • C-ambient CO 2 treatment, whereas the highest starch concentration in the seeds was found in the 31 • C-enriched CO 2 treatment, even though differences were not significant (Figure 6c,d and Table 6).

Temperature as the Main Driver of Seedling's Trait Responses
Temperature was the dominant driver of E. acoroides seedling's trait responses within this experiment, and there were no synergistic effects in the interaction of both factors. Therefore, our initial hypothesis was rejected. In general, our results showed that seedlings of seagrass E. acoroides perform better at the higher temperature treatments, as bigger leaves with high maximum leaf lengths, leaf SA, and weight were always observed. This is in accordance with what was observed in the field in previous studies with adult individuals, as the temperature was also the major factor influencing the growth rate with increasing values above 30 • C [27]. This pattern of an enhanced size of leaves under an increasing temperature has also been observed in seedlings of Zostera marina [19] while detrimental or no changes were observed in above-ground traits of other persistent species with big-size seedlings, such as P. oceanica, at sub-lethal temperatures [15,17]. The optimum temperature for photosynthesis in adult E. acoroides from Philippines ranged between 24 and 33 • C, with an optimal temperature set at 27 • C [35] or 33 • C [30]. Thus, based on the photosynthetic potential of PSII (Fv/Fm) observed in this study, seedlings are still within the species' optimal range. However, rETR max (photosynthetic capacity) and Ek (saturation irradiance) values were higher in the seedlings' grown at regimes of 27 • C than in the ones grown at 31 • C, showing that some down-regulating mechanism might be potentially initiated at high target temperatures. These detrimental effects on photo-physiology were also observed in P. oceanica seedlings at increasing temperatures [14,17]. Even though, rETR values were higher than in younger seedlings of the same species [18], which indicates a higher independence in maternal reserves in comparison to the younger leaves.
Leaf nitrogen content under regimes of 27 • C was within the range observed in adult E. acoroides plants under natural conditions in South-East Asia meadows [34]. Contrary to this, carbon content was lower than in younger E. acoroides seedlings [18] showing that internal reserves were being used. This has also been observed in P. oceanica seedlings of a similar age [20]. Both leaf nitrogen and carbon contents, however, were also significantly affected by an increasing temperature. Lower leaf nutrient content under regimes of 31 • C was observed, which was likely related to nutrient reabsorption, often occurring under higher growth rates due to dilution of nutrients in leaf tissue [54,55]. On the contrary, higher carbon content in below-ground tissues were detected at the high target temperature, which has been shown in previous studies with P. oceanica and E. acoroides seedlings [17,18]. Higher sucrose concentrations observed in seeds might be a consequence of nutrient reallocation, which is a common mechanism in adult seagrasses under unfavorable environmental conditions [13].
Therefore, while higher temperatures are positive in terms of growth, performance in the regimes of 31 • C suggests that seedlings were close to suboptimal temperatures in terms of photo-physiology and nutrient content. To further test temperature tolerance of seedlings, experiments under higher temperatures (above 32 • C) would be necessary. Therefore, the small differences observed could suggest that higher average temperatures within the current climate change scenario are not detrimental for seedling performance as they can depend on their internal reserves from weeks until months [18,29,[56][57][58][59]. This physiological dependence on maternal seed reserves might help them to ameliorate the effects of environmental factors, as shown in this study.

Interactions of Temperature and CO 2 Concentrations on Seedling's Trait Responses
In our study, enriched CO 2 concentrations alone did not significantly affect performance of E. acoroides seedlings in terms of its photo-physiology or growth, nor did they affect the carbon balance of the seedlings, suggesting that they are still dependent on their internal carbon reserves. Previous studies have shown that, in many cases, seagrass are C ilimited at the current oceanic concentrations. Therefore, an increasing CO 2 concentration has often been related to an overall better performance, in terms of growth or photosynthetic rates, in both temperate and tropical species, including E. acoroides [31,[60][61][62][63]. For instance, Thalassia hemprichii, which is another persistent tropical seagrass species, can increase its growth rate under CO 2 enrichment (~1000 µmol L −1 ) [64] as well as Z. marina, which showed an increase in growth in CO 2 cultures enriched two-fold when compared to the control [65].
In terms of the biochemical traits, the absence of effects of CO 2 on leaf nutrient contents suggests, as previously observed with increasing nitrogen concentrations in seawater, that seedlings are still sustained by internal nutrient reserves [18]. Carbon balance in the body of a plant is a major factor determining the rate of growth, which also occurs in seagrass [66]. Changes in carbon balance increasing the storage in below-ground tissues have been observed in other seagrass species under enriched CO 2 concentrations for longer periods than in this study [20,67], suggesting it may be a long-term strategy of plant survival.
Therefore, contrary to the beneficial effects shown in other previous studies in marine primary producers [13], our results did not show any clear positive effect of CO 2 enrichment. We did, however, observe interactions with temperature that could benefit seagrass development in this initial seedling phase. This interaction of increasing CO 2 concentrations on thermal acclimation, even within sublethal temperatures, is still not well understood in seagrasses.
The study results showed enhanced photo-physiological activity in terms of higher rETR max and Ek values in the 27 • C-enriched CO 2 treatment, but not under the regimes of 31 • C, showing a weak interactive effect of both factors (p-value = 0.078). Even though CO 2 concentrations in the 31 • C-enriched treatment were greater (~600 µatm) than in the 27 • Cenriched CO 2 treatment (~465 µatm) due to increased CO 2 solubility with temperature, we did not observe an enhanced performance in the 31 • C-enriched CO 2 treatment. This higher rETR max (photosynthetic capacity) was not related to the biggest leaf morphological traits, which were observed under the 31 • C-enriched CO 2 treatment nor with higher internal nutrient contents. However, a trend of increasing concentrations was observed in enriched CO 2 treatment. This interaction could be more important in a tropical environmental as, even though some sites where E. acoroides grows are highly eutrophic [68,69], these systems are typically oligotrophic and dissolved nutrients are usually the limiting factor [34]. A better understanding of responses at a biochemical or morphological level has shown that the high photosynthetic performance observed in some studies is more a short-term response to C i limitation that is difficult to maintain over time due to the potential limitation of other factors [22,70,71]. This was also observed with P. oceanica seedlings, which showed higher photosynthetic performance after two months of study than at the end of the study period [20].
We also found a significant interactive effect in the carbon content in roots (p-value = 0.049), with the highest concentration in the 27 • C-enriched CO 2 treatment, but with no relation to root biomass. In this treatment, we also observed the highest nutrient content in leaves. Therefore, the relative contribution of root versus leaves in meeting carbon or nitrogen requirements might be influenced by the higher photosynthetic performance under enriched CO 2 treatments (i.e., the highest value of rETR max ).
These surprising results are becoming more common in recent studies [72,73], showing that seagrass responses to CO 2 enrichment are more complex than suggested in previous studies [74].

Different Trait Responses and Dependence on the Internal Reserves
From the discussion above, it is clear that there is no single trait that can clearly indicate the effect of different stressors, and results are even more confounding when multiple stressors are influencing seedling performance. Physiological and biochemical traits in seagrass have usually been highlighted to be the first to respond under different stressors [75], which is followed by morphological and community traits (i.e., biomass loss). After 35 days, most leaf traits have been influenced by temperature, enhancing their production in terms of size. However, physiological and biochemical traits were negatively influenced by this factor, including photo-physiological activity and nutrient content in leaves. This also suggests that it is important to look to a combination of traits to get a more holistic picture of plant responses.
Overall, leaf traits responded more quickly to the experimental treatments, as we observed no responses in seeds or roots, with the exception of carbon content in this latter tissue. Therefore, leaf tissue is a better early indicator of temperature or interactive CO 2 effects in seedlings. Although there is not much experimental evidence with adult E. acoroides individuals, based on their life-history traits, this species is classified as persistent with very large blades and slow growth. Therefore, using seedling traits as indicators of climate change effects in this species might be an advantage as faster responses could be obtained.
Seedlings have large maternally derived nutrient resources, so they do not appear to be limited by carbon, or other nutrients [18]. This lack of physiological independence likely buffers the seedlings from environmental stressors, leading to small differences in the responses. Additionally, seedlings are genetically unique individuals that have not been selected for by the natural environment via the "environmental sieve" (sensu Harper). Consequently, they may not respond the same as genotypes that have been selected by survival in the field (adult plants) and may not be representative samples of the population. In summary, we observed deleterious effects on seedlings, as lower nutrient content in leaves, photo-physiological performance, or the prioritization of the development of above-ground tissues. However, contrary to what was observed with adult seagrass, CO 2 enrichment did not help seedlings to cope with high temperatures. Further studies on the long-term ecological consequences of the changes in these traits are needed to better understand how seedling's development influences adult populations.
Supplementary Materials: The following are available online at https://www.mdpi.com/2076-341 7/11/6/2730/s1. Figure S1: The effect of increasing temperature and CO 2 enrichment on (a) seed diameter, (b) seed height, (c) seed weight, (d) root number, (e) maximum root length, and (f) root weight of Enhalus acoroides seedlings at the end of the experimental period. Values are mean ± SE (n = 6). Table S1: Results of analysis of variance (two-way ANOVA test) of the effects of temperature and CO 2 and their interaction on the morphological seed and root traits ( Figure S1). No statistical differences were detected (p-value < 0.05).  Institutional Review Board Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.