Salicornia europaea L. Functional Traits Indicate Its Optimum Growth

Salicornia europaea L. grows in areas periodically flooded by salty or brackish water. It has potential economic value, because it can be used as food, forage, or biofuel, and has potential in pharmaceuticals and cosmetics. Increasing interest in S. europaea is due to its extreme salt tolerance and well growth in marginal saline soils. However, the variation in its functional traits in response to environmental conditions is still poorly studied. There are still questions regarding the optimal level of salinity for different traits. Therefore, we worked to address the question if S. europaea traits from different scales are controlled by salinity level. Based on performed pot experiment, we found that almost all traits are salinity dependent but affected in different ways. We demonstrated that morphological, biomass, and anatomical properties indicate optimum growth between 200 and 400 mM NaCl and growth limitations at 0, 800, and 1000 mM NaCl. Moreover, we found the most affected traits which include photosynthetic pigments and protein content, plant surface area, peroxidase activity, and anatomic traits related to cell shape. Our results significantly expanded the knowledge about S. europaea functional traits variation in response to salinity, which can be important for discovering regulating processes and for possible future agricultural applications.


Introduction
Salicornia europaea L. belongs to the Amaranthaceae family (formerly Chenopodiaceae), subfam. Salicornioideae. The genus Salicornia is widespread in temperate and subtropical regions of the Northern Hemisphere but absent in South America and Australia [1]. Presence in habitats with changing seasonal and even daily dynamics has led to high physiological plasticity in plants [2]. This results in phenotypic variability and problems with establishing acceptable systematics [3]. Despite the high phenotypic variability, several common features can be distinguished [4,5]. S. europaea has an erect, highly branched stem. Secondary shoots are formed on the primary cylindrical shoots. The plant has strongly reduced leaves, and the assimilation area is located in the shoots. It is green most of its life, but at the end of life cycle, the stems turn red due to chlorophyll destruction which reduces photosynthesis and finally affects nutrient loss, biomass, and hydric balance [6]. Branches have spikes consisting of three flowers, one main and two laterals. Seeds are small, dark, ellipsoidal, and characterized by heteromorphism, i.e., color, shape, and size variability [7].
S. europaea grows in areas periodically flooded by salty or brackish water [8]. In Central Europe this species has been recognized as Salicornia ramosissima J. Woods (=S. herbacea L.) [9]. It has potential economic value because it is edible, either raw or cooked [10]. In addition, physiological traits responses [19,20,22,40]. The biochemical parameters, although useful, are black boxes in terms of anatomical and structural changes, which are not directly visualized. This is currently a gap in the literature for this species. To our knowledge, there are no comprehensive studies considering traits in the context of the adaptation and optimum growth of S. europaea in the salinity gradient. The level of salinity is still under question, which can be considered as optimal for S. europaea trait development. Moreover, it is still unknown which S. europaea functional traits are the most affected by salinity.
Therefore, to fill this gap in the knowledge the overarching question, we worked to address whether S. europaea traits are from different scales controlled by salinity level. To answer this question, we performed complex research on the morphological, anatomical, and physiological traits at different salinity levels. For morphological and anatomical assessments, we applied a novel image analysis method [5]. To present a complex example of the plant trait functional linkage, we applied similarity analysis between different salt treatments [41]. We also selected the most affected by salinity functional traits by the means of discriminant analysis [42]. We hypothesized that: (a) salinity affects plant morphological, anatomical, and physiological responses in different ways, and (b) plant trait responses can indicate optimum growth in the salinity gradient. Understanding complex mechanisms of salt stress adaptation of S. europaea, an extreme halophyte species, is possible only based on traits from different functional levels. The determination of key functional traits in salinity adaptations and S. europaea optimum growth is also important because of the possible future agricultural application perspectives.

Growth Responses to Different Salinity Levels
To investigate the morphological trait responses to salt stress, we measured plant height, number of branches, plant surface area, and shoot diameter. We found differences in the morphological features dependent on the salt concentrations (Figures 1 and 2). Plants which grew under extreme salinity (1000 mM NaCl) were smaller than those grown in other salt treatments. They had the smallest plant surface area (ca. 105 cm 2 ) and were significantly shorter (ca. 4.6 cm) and thinner (diameter ca. 0.248 cm) than the others ( Figure 2). Moreover, they had the lowest number of branches (ca. 9). There were not significant differences in plant height and shoot diameter between 0, 200, 400, and 800 mM NaCl treatments (Figure 2a,d). The plant surface area was higher in 400 (283 cm 2 ) and 800 mM NaCl (ca. 289 cm 2 ) ( Figure 2c). The highest number of branches (ca. 26) was noted at 800 mM, but there were no significant differences between 800, 400, and 200 mM treatments ( Figure 2b). Therefore, the relationship between this trait and salinity is not so clear.
For the effect of salinity on biomass accumulation, we measured shoot, root, and total fresh and dry weight. The highest values for fresh and dry weight (SFW 38.5 g, RFW 12.1 g, TFW 50.6 g, SDW 5.4 g, RDW 4.1 g, and TDW 9.5 g) were obtained in 400 mM NaCl, although for all these traits there were no statistically significant differences between 400 and 200 mM NaCl ( Figure 3). The lowest values for fresh and dry weights (SFW 2.03 g, RFW 0.44 g, TFW 2.47 g, SDW 0.26 g, RDW 0.08 g, and TDW 0.34 g) were obtained in 1000 mM NaCl. In general, for measurements of fresh and dry weights, there were not significant differences between 0, 800, and 1000 mM NaCl ( Figure 3). These results confirm morphological traits assessments, i.e., under 1000 mM NaCl, S. europaea cannot grow well. Moreover, growth limitations were detected at 0 and 800 mM NaCl.

Anatomical Responses to Different Salinity Levels
We based the effect of salinity on anatomical traits on measurements of area, perimeter, diameter, roundness, and aspect ratio of stem-cortex cells ( Figure 4). Our results showed that in the control and extremely saline condition (0 and 1000 mM NaCl), the plant cells area, cell perimeter, and cell diameter were lower than in moderate and high saline treatments (200, 400, and 800 mM NaCl) (Figure 5a-c). The plant cells area in the non-saline condition and 1000 mM NaCl was ca. 8394 and 9027 µm 2 respectively, the plant cell perimeter was ca. 345 and 346 µm, and the cell diameter was ca. 130 and 126 µm. There were not statistically significant differences between 200-800 mM NaCl treatments and cells area ranged between 13.360 and 15.967 µm 2 , cell perimeter between 431 and 479 µm, and cell diameter between 157 and 177 µm (Figure 5a  These results suggested that under non-saline and extreme saline conditions, stemcortex cells are smaller probably because of growth stress. We did not find statistically significant differences in anatomical traits between plants grown in 200-800 mM NaCl.

Photosynthetic Pigments and Soluble Protein Content
Photosynthetic pigments are traits that can affect photosynthetic performance, plant growth, and development. We measured concentrations of chlorophyll a and b and carotenoids. Results demonstrated decreasing photosynthetic pigments content with increasing NaCl concentrations (Figure 6a-c). The highest content was found in control condition (chlorophyll a: 0.715 mg/g FW; chlorophyll b: 0.427 mg/g FW and carotenoids 0.254 mg/g FW). The lowest content was found in 1000 mM NaCl (chlorophyll a: 0.241 mg/g FW; chlorophyll b: 0.151 mg/g FW and carotenoids 0.092 mg/g FW). Soluble protein content was more stable than pigments because, up to 400 mM NaCl, there were not significant changes in their content (ca. 6.0 mg/g FW) ( Figure 6d). However, a significant decrease was observed in 800 (ca. 2.39 mg/g FW) and 1000 mM NaCl (ca. 2.35 mg/g FW) (Figure 6d).

Hydrogen Peroxide (H 2 O 2 ) Content
Hydrogen peroxide is a strong ROS (reactive oxygen species) which is overproduced in plant cells under stressful conditions. To investigate how different salinity levels can affect ROS production we measured H 2 O 2 content. We observed that under 0 (ca. 4.4 nmol/g FW) to 800 mM NaCl (ca. 14.2 nmol/g FW) the level of H 2 O 2 did not differ significantly ( Figure 7a). However, we observed a statistically significant increase of the H 2 O 2 content in extreme salinity-1000 mM NaCl (ca. 26.2 nmol/g FW). It seems that S. europaea has a high tolerance to increasing salinity levels and can prevent H 2 O 2 overproduction.

Antioxidant Enzymes Activities and Proline Content
Antioxidant enzymes activities and osmolyte accumulation are two of the most important defense mechanisms against salinity. Therefore, the evaluation of these traits can help in a better understanding of S. europaea adaptive response. Thus, we measured activities of two antioxidant enzymes (CAT and POD) and the concentration of proline as an osmolyte accumulation marker. We found that with increasing NaCl concentration, CAT activity decreases. The highest activity of catalase was recorded in plants growing without salinity (36.7 µM min −1 mg −1 protein), and the lowest in 1000 mM NaCl (5.7 µM min −1 mg −1 protein) (Figure 7b). These results suggested that catalase in S. europaea is sensitive to salt concentration and in low salinity has high activity and efficiency. In contrast, POD activity ( Figure 7c) increased with increasing salinity. The lowest POD activity was in plants treated with 0 mM NaCl (3.7 µM min −1 mg −1 protein), and the highest activity in 1000 mM NaCl (39.2 µM min −1 mg −1 protein). In the case of proline, we found any statistically significant differences among all salinity levels ( Figure 7d). The proline content in plant tissues ranged from ca. 0.4 µmol/g FW in non-saline conditions to 0.7 µmol/g FW in 1000 mM NaCl.
The results suggested that S. europaea without NaCl in the medium had a stable proline concentration in tissue which helped it to cope with osmotic stress induced by soil salinity.

Discussion
Halophytes are plants that have evolved to survive under high salinity conditions, and many of them are thought to require salt exposure for better growth [43]. These plants can regulate their energy metabolism under saline conditions [43,44]. Thus, a better understanding of halophyte functional traits adaptations is important not only for understanding the natural salt marsh ecosystem function [45,46], but also for reconsidering salt tolerance in glycophyte plants to better manage salinity problems in the future agriculture [47]. As it was already mentioned, up to our knowledge, there is a lack of research involving salt stress responses of complex traits from different functioning levels. Addressing our research to answer the question if S. europaea traits from different scales are controlled by salinity level, we found that almost all of them were salinity dependent. However, we confirmed our first assumption, i.e., that traits are affected in different ways.
Analysis of growth responses expressed by morphological properties and biomass accumulation demonstrated that the extreme saline condition of 1000 mM negatively affected these traits. Moreover, biomass accumulation was also reduced at 0 and 800 mM compared with 200 and 400 mM NaCl. This is in line with research that already reported stimulation of S. europaea and other halophytes biomass accumulation at moderate salinity [48,49]. Research by Lv et al. [28] demonstrated limitation of plant height, stem diameter, and plant biomass in lower salinity than recorded by us, i.e., of ca. 700 mM NaCl. On the other hand, in line with our findings, they also proved S. europaea growth traits limitations in the substrate without salinity. Similar results were obtained by Piernik [27], who reported lower S. europaea height and biomass at a very low salinity of ca. 2 dSm −1 (ca. 20 mM NaCl) in field experiment. However, in this treatment, opposite to our findings, higher number of branches were noted. Additionally, the results of field research by Silva et al. [50] demonstrated a higher number of S. europaea branches at low salinity. It is worth emphasizing that morphological parameters under field conditions can be affected not only by salinity but also by plant density and competition between individuals [51]. Obtained results may suggest that plant biomass is more directly salt stress-dependent than morphological traits.
In the present study, an analysis of the anatomical traits demonstrated that the S. europaea stem-cortex cells under non-saline conditions and in 1000 mM NaCl were significantly smaller than in other salt treatments. However, we observed that under control conditions, plant cells were more elongated, while in higher salinities, they were rounder. We did not find statistically significant differences in most anatomical traits between plants grown in 200-800 mM NaCl. Our results are in line with other studies on halophytes that show the increasing cells size in water-absorbing tissues under moderately saline conditions [16]. The water storage in the parenchyma is the first S. europaea mechanism to overcome high Na + concentrations. This dilutes the accumulated salts and contributes to maintaining cellular turgor, allowing the plant to cope efficiently with high salinity [12,16]. In addition, it was reported that salinity induces vacuolization in many halophytes and may be responsible for the plant adaptation to salt stress [52]. For example, Akcin et al. [52] found that under stressful conditions the juice of halophytes increases. This response is associated with an increase in cell volume due to extra water being stored in the vacuoles for the plant to survive. Moreover, it was reported that S. europaea cells can remain turgid and continue proper cell function by ion compartmentalization in cell vacuoles. In this sense, the cytoplasm and organelles of the cell are protected from salt [12]. That is why our results also proved that S. europaea, by enlarging its cells under moderately saline conditions, can adapt to high salt stress. On the other hand, the reduction of cell size at extreme salinities can be explained by reduced ability to osmoregulate because of saturation of the system with Na + and Cl − or deficiency of Ca 2+ and K + , which are involved in almost all reactions related to plant development [2].
Our results showed that salinity also affects biochemical traits of S. europaea. First of all it negatively affects photosynthetic pigments. We observed that the highest concentration of photosynthetic pigments was under non-saline conditions, and the lowest under the highest salinity level (1000 mM NaCl). These results suggested that salinity significantly affects traits which are responsible for photosynthesis performance in plants. Reduction in photosynthetic pigments due to salt stress has been well documented in numerous papers [5,53,54]. This happens by the inactivation of enzymes involved in the synthesis of photosynthetic pigments [55,56]. Furthermore, this reduction can be due to ROS generation and increasing of chlorophyllase enzyme activity [56,57]. However, we found that photosynthetic pigment concentrations were relatively stable at levels above 400 mM NaCl, indicating that the plant can protect, to a certain extent, these traits under extreme saline conditions. High content of photosynthetic pigments in treatments without salinity was not related to high biomass accumulation, which suggests mechanisms of growth limitations other than those directly related to photosynthetic ability. We also calculated Chlorophyll a/b ratio to test if the rearrangement of chlorophyll contents can explain observed by the pattern in S. europaea growth performance, as it was reported, for example, in Arthrocnemum macrostachyum and Sarcocornia fruticosa by Ghanem et al. [58] (Supplementary Figure S1). Based on the obtained results, we can conclude that such a rearrangement was not a strategy of the investigated species.
We also observed that salinity negatively affects soluble protein content of S. europaea. This response may be due to the toxic effect of NaCl on protein synthesis, or the proteolysis of proteins caused by proteases induced by salt stress [59]. It is well documented that high concentrations of NaCl destroy the hydration layer of protein, causing its aggregation and denaturation [60]. However, soluble protein compounds were more stable than pigments and maintained a similar level up to 400 mM NaCl.
Hydrogen peroxide, as already stated, is one of the reactive oxygen species responsible for oxidative stress. On the other hand, H 2 O 2 is widely generated in many biological systems and mediates various physiological and biochemical processes in plants [61]. Our results showed that the H 2 O 2 content was the lowest in the control condition, and by increasing salt levels up to 800 mM NaCl, no significant increase was observed. Only at 1000 mM NaCl was its level significantly higher. The opposite pattern, i.e., increasing level of H 2 O 2 in plant tissues together with increasing salinity, has been observed for glycophytic species [54,62,63]. It seems that halophytic S. europaea can prevent H 2 O 2 overproduction. Nevertheless, in extreme saline conditions, impairment of H 2 O 2 production and scavenging can affect membrane structural integrity and peroxidation of lipids and limits plants growth and development [19].
In present study, we also monitored traits related to antioxidant enzymes activities and proline content changes. We found that the highest activity of catalase was under non-saline conditions and decreased together with increasing salt levels. Lower activity of CAT in high salt concentrations indicated CAT as a less effective scavenger of H 2 O 2 . Furthermore, the CAT has a poor affinity for H 2 O 2 and exhibits photo-inactivation and subsequent degradation [64,65]. A completely opposite activity had POD with the lowest values at 0 mM, and the highest at 1000 mM NaCl. In addition, up to 400 mM NaCl POD activity was relatively stable. Increase of the activity of POD at higher salinity levels indicated that this enzyme plays a key role in defense mechanisms of S. europaea by scavenging ROS from cells.
Although it is reported that accumulation of proline in S. europaea cells is an important adaptive response to salinity [18,43,66], in this study we did not prove a statistically significant increase in proline content by increasing salt levels. Some halophytes produce proline analogues, e.g., glycine betaine under salt stress to survive due to their ability to protect the protein turnover machinery, stabilize proteins, and prevent enzymes from denaturation [66,67]. It seems that S. europaea normally has stable proline levels in tissues, which helps it to cope with osmotic stress. However, as proved by Pellegrini et al. [68], proline could be involved in stress tolerance in the Salicornia genus regardless of the intensity and duration of the stress. For glycophytic species, a frequently increasing level of proline is observed together with increasing salinity [22,54].
We also referred biochemical traits to dry weight instead of fresh weight to test if differences observed between samples in terms of FW could be partly due to differential loss of water under stress ( Supplementary Figures S1 and S2). The tendency in response to salinity level was maintained, and therefore we can conclude that the differences between samples are due to the different levels of stress to which they were subjected.
To conclude on salinity impact, based on all investigated traits and to test our second assumption about indication of plant optimum growth, we performed non-metric multidimensional scaling (NMDS) with Bray-Curtis dissimilarity measure between treatments. Within the analysis we asked for original variables related to the resulting NMDS space as supplementary ones to make the interpretation of the results more clear (Figure 8). Supplementary response variables do not affect the definition of the ordination NMDS axes determined by similarity between samples and can be added to an existing ordination by projection, i.e., by regressing its data on the existing ordination axes [69]. NMDS results demonstrated that some traits can indicate the best S. europaea growth expressed by morphological, biomass, and anatomical properties between 200 and 400 mM NaCl ( Figure 8). The role of morphological traits adaptations were not stressed up to now. Plants grown in 200 and 400 mM NaCl were the most similar to each other. The first ordination axis explained 98% of the variability between these treatments located at the left side of the diagram and 0 and 1000 mM NaCl located at the right side ( Figure 8). The second ordination axis related to the differences between plant traits in 0 and 1000 mM NaCl explained only 1.7% of traits variability, related mostly to oxidative stress. This findings demonstrated by morphological and biomass traits are in line with the field studies of Piernik [27] and the laboratory observations of Lv et al. [28], Cárdenas-Pérez et al. [5], Muscolo et al. [2], and Rozema and Schatb [70] who reported S. europaea's optimum growth at moderate salinity and growth limitation under non-saline conditions. Based on field research and soil sampling on inland salt marshes in Central Europe, Piernik [71] reported the optimum growth of S. europaea at 38.1 dSm −1 , which corresponds to ca. 380 mM NaCl and is in the range of our findings. The growth stimulation at low and moderate salinity in halophytes such as S. europaea may be attributed to the improvement in shoot osmotic status because of increased ion uptake [72]. However, the very high salinity imposed a reduction of the growth, which is probably associated with the reduced turgor and the high energy cost of massive salt secretion and osmoregulation [20,48,73].
Based on discriminant analysis with a forward selection procedure and Monte Carlo permutation test, we selected the most affected traits by salinity. There were chlorophyll a, carotenoids, and protein content explaining ca. 23-25% of the variability between all treatments (Table 1), with the highest values related to low saline conditions. A similar amount of variability, i.e., ca. 25% was explained by plant surface area related to higher salinities. Photosynthetic area, represented in S. europaea by plant surface area, determines light interception and is an important parameter in determining plant productivity [74]. Of lower importance but statistically significant in the separation between treatments were POD activity related to the highest salinities, explaining ca. 0.8% of variability, and aspect ratio related to non-saline treatment, explaining 0.7%, and the number of branches and cell perimeter related to the moderate salinities, explaining respectively 0.6% and 0.2% of variability between treatments (Table 1). However, as was already mentioned and discussed, the number of branches can be strongly modified under field conditions not only by salinity but also by competition between individuals for resources [51]. We focused here on conditional effects, which exclude the effect of the most correlated traits. Conditional effects summarize the partial effect of each predictor, representing the variation (and its significance) explained by a predictor after accounting for the effect of predictors already selected [69]. The predictors were chosen in the order of their decreasing explained variation. We skipped simple effects, which summarize the independent effects of all explanatory variables, because they reflected significance assessed by already presented ANOVA. Our results demonstrated that though photosynthetic traits are the most affected by salinity, it does not guarantee high productivity of S. europaea in non-saline environments and does not affect productivity to the certain salinity level. This can explain why many halophytes are known as light-required species [75]-because they cope with reduced photosynthetic pigments under saline condition. Fortunately, a relatively higher S. europaea area under moderate salinity is advantageous for increasing photosynthetic capacity through capturing light at the expense of increased construction costs and produces large boundary layers responsible for reducing transpiration rates and thus heat exchange and carbon dioxide diffusion to the surrounding air [76]. Important findings in our research are also addressed to the anatomic traits most affected by salinity and related to cell shape. As proved by aspect ratio with increasing salinity, cells become rounder, optimizing their surface area and, as was evidenced by cell perimeter, they become enlarged. However, all regulatory mechanisms behind cell adaptations are still not recognized [21]. Our findings can be the key starting points for their future identification.  Finally, we would like to highlight that S. europaea is a species covering different ecotypes/subspecies, which can be specifically adapted to the local environments [5,12,21]. That is why we reported, in the introduction, Salicornia ramosissima J. Woods (= S. herbacea L.) from Central Europe [9] as a synonym/reference for the S. europaea population investigated by us.

Laboratory Plant Material and Plantlets Preparation
S. europaea seeds were collected from the nature reserve of halophytes "Ciechocinek" located in north-central Poland (52 • 53 N, 18 • 47 E; Central Europe) under permission WOP.6400.12.2020.JC. Matured plants spikes were harvested at the beginning of November, and the seeds were shaken from the spikes. After threshing, 250 seeds that were healthy and uniform in appearance were selected and planted in Petri dishes with Whatman no. 2 filter paper. The Petri dishes were watered with 5 mL distilled water and were placed in a growth chamber with 22 • C and 16 h light period for 9 days. After full germination, uniform plantlets were selected and planted into the pots that contained a mixture of vermiculite and sand (1:1). In total, 60 pots were prepared and divided into five groups (12 pots per each salt concentrations: 0 mM-distilled water, 200 mM, 400 mM, 800 mM, and 1000 mM NaCl). Before planting, each group of 12 pots was located on individual tray lacking drainage and were saturated to their full capacity with relevant NaCl solution (ca. 0.5 L per treatment, pot size: height 5.3 cm, diameter 5.5 cm). All pots were placed in a growth chamber with 22 • C and a 16 h light period and each group was irrigated through pouring distillate water once per week. To prevent nutrient deficiencies, plants were watered with Hoagland's solution once per two weeks [77]. In our previous research we focused on germination and early growth assessments of S. europaea [2]. Therefore, this time we focused on plantlets, which were grown for 60 days to the end of the experiment period.

Growth Analysis
Growth analysis covered morphological and biomass assessments. Two months after starting salt treatments, the plants were subjected to morphological analysis. For this purpose, plant images were taken by a Canon camera with resolution of 5 to 25 MP for morphological analysis. Images were taken from four sides of each sample (see Figure 1) in four replications. We used ImageJ version 1.47 program (National Institutes of Health, Bethesda, MD, USA) for image analysis to measure plant morphological parameters: height, number of branches, surface area, and shoot diameter. The results were calculated according to the number of pixels covering the plant converted to the appropriate metric units.
At the end of the experiment, shoot fresh weight (SFW), root fresh weight (RFW), and total fresh weight (TFW) were measured. After the samples were oven dried for 72 h at 72 • C, shoot dry weight (SDW), root dry weight (RDW), and total dry weight (TDW) were determined.

Anatomical Analysis
The anatomical analysis was performed based on the cross-sections of plant shoots (middle primary branch-fleshy segment shoot) taken from plants growing in five NaCl concentrations. Slices of fresh tissue were obtained by cutting them with a sharp bi-shave blade and slices of approximately 0.5 mm were used for analysis. The samples were stained with alkaline fuchsin and malachite green for microscopic observations. S. europaea crosssections photos for each NaCl concentration were taken using light microscopy. Photos were used to observe the cortex parenchyma cells, because water storage in the parenchyma is the first S. europaea mechanism to overcome high Na + concentrations [28]. We used ImageJ version 1.47 (National Institutes of Health, Bethesda, MD, USA) for image analysis. For each sample, more than 200 cells were analyzed to determine the following parameters: (1) cell area (A), which was calculated through the number of pixels inside the borderline; (2) cell perimeter, which was calculated due to the prescribed limits; (3) cell diameter, which was defined by the distance between two points; (4) aspect ratio (AR), which was defined as the quotient between the minimum and maximum diameter, determining the uniformity of the cell; and (5) cell roundness (R), which determines the circularness of a cell using Equation (1). In this equation, for perfectly round cells R = 1 [78].
where: A is cell area; MD is cell diameter.

Biochemical Analyzes
To evaluate the effect of different salinity levels on the biochemical traits of S. europaea plants, we measured traits related to photosynthetic activity, osmotic adjustment, oxidative stress, and enzymatic activities.

Photosynthetic Pigments Content
Chlorophylls (Ch a and Ch b) and carotenoids were extracted from fresh plant stems (100 mg) using 80% acetone for 6 h in darkness, and then centrifuged at 10,000 rpm, 10 min. Supernatants were quantified spectrophotometrically. The absorbances at 646, 663, and 470 nm wavelength were measured. The total content of chlorophyll a and b [79] and carotenoids [80], when 80% of acetone was used as dissolvent, were calculated according to Equations (2)-(4) and reported as mg-per-gram fresh weight.

Soluble Protein Content
Plant material (0.5 g) was placed in a mortar chilled with liquid nitrogen, then 7 mL of 50 mM phosphate buffer was added and homogenized. After fine grinding, the sample was centrifuged at 20,000× g for 25 min at 4 • C. The supernatant was collected for determination of protein content at 595 nm using the Bradford method [81]. The amount of protein was determined using a bovine serum albumin as a standard.

Hydrogen Peroxide (H 2 O 2 ) Content
Hydrogen peroxide was measured according to the methods described by Velikova et al. [82]. Three samples were prepared for each NaCl concentration. Stem tissues (0.5 g) were homogenized in an ice bath with 5 mL of 0.1% TCA (trichloroacetic acid). Then, the homogenate was centrifuged at 12,000× g for 15 min at 4 • C. The supernatant (0.5 mL) was transferred to a new tube and 0.5 mL of 10 mM potassium phosphate (pH 7) and 2 mL of 1 M KI were added. The solution was incubated in the dark for one hour and then the absorbance was read at 390 nm wavelength and the hydrogen peroxide concentration was given based on the standard curve from 0 to 40 nM and equation y = 0.0188x + 0.046, R 2 = 0.987 in nM-per-gram fresh weight.

Peroxidase (POD) Activity
The peroxidase activity was measured according to the method of Chance and Maehly [83]. For the measurement, a 3-mL reaction mixture was prepared with the following composition: 50 mM potassium phosphate buffer (pH 7), 20 mM guaiacol, 40 mM H 2 O 2 , and 0.1 mL of enzyme extract obtained from plant stems. The reaction was initiated by adding the enzyme extract. The increase in absorbance of the reaction solution was measured at 470 nm wavelength. The reading was recorded every 20 s. One unit of peroxidase activity was defined as an absorbance change of 0.01 units per minute. The enzyme activity was expressed on a protein basis [84].

Catalase (CAT) Activity
Catalase was determined by measuring residual H 2 O 2 with a titanium reagent [85]. The reaction mixture was 3 mL in total and consisted of: 1 mL of 6 mM H 2 O 2 and 1.9 mL of 0.1 M phosphate buffer (pH 7). To start the reaction, 0.1 mL of diluted enzyme extract was added to the test tube. The reaction was stopped after 5 min by adding 4 mL of titanium reagent, which formed a colored complex with residual H 2 O 2 . Mixtures without enzyme were used as control. The colored complex was centrifuged at 10,000× g for 10 min. The absorbance was measured at 415 nm. The residual H 2 O 2 content was calculated from the standard curve [86].

Proline Content
The proline content of the leaves was determined according to the method of Bates et al. [87]. Three samples were prepared for each NaCl concentration; 0.5 g of each sample were pulverized on ice and homogenized in a mortar with 5 mL of 3% aqueous sulfosalicylic acid. The homogenate was centrifuged at 18,000× g, 10 min at 4 • C, and the supernatant was collected. About 1 mL of the supernatant was added into a test tube and 2 mL of glacial acetic acid, 2 mL of ninhydrin reagent, and 1 mL of sulfosalicylic acid 3% were added to the tube. The reaction mixture was boiled in a 100 • C water bath for one hour. Then, the test tubes were put on ice and 4 mL of toluene were added to them. The samples were transferred to a separating funnel and after thorough mixing, the toluene containing the chromophore was separated. The absorbances were measured at 520 nm wavelength. The amount of proline was determined using the standard curve in the concentration range of 0 to 40 µg mL −1 and equation y = 0.0467x − 0.0734, R 2 = 0.963.

Statistical Analyzes
To compare all treatments, one-way analysis of variance (ANOVA) was performed followed by post hoc analysis using Tukey's test. The PAST 4.0 program [88] was used for statistical analyses. To discuss S. europaea's optimum growth based on the whole set of the traits, we applied non-metric multidimensional scaling (NMDS) with Bray-Curtis dissimilarity measure to demonstrate similarities between plants coming from different salinity treatments. To select the most affected traits by salinity treatments, we applied discriminant analysis with a forward-selection procedure and Monte Carlo permutation test. Only conditional effects were taken into account to exclude the effect of the most correlated traits [42]. For these analyses, the Canoco 5.0 program was applied [69].

Conclusions
Addressing our research to answer the question if S. europaea traits from different scales are controlled by salinity level, we found that almost all of them were salinity dependent. However, we proved that functional traits were affected by salinity in the different ways and demonstrated significant differences at different salinity levels. Moreover, we did not find a statistically significant relationship between proline levels and increasing salinity, which was not expected based on reported findings regarding glycophyte species. Based on analysis of all investigated traits, we demonstrated that morphological, biomass, and anatomical properties indicated optimum growth between 200 and 400 mM NaCl and growth limitations at 0, 800, 1000 mM NaCl. Moreover, we can conclude that the most affected traits by salinity include photosynthetic pigments, protein content, plant surface area, peroxidase activity, and anatomic traits related to cell shape. Our results significantly expand the knowledge about S. europaea functional traits variation in response to salinity, which can be important for discovering regulating processes and for possible future agricultural applications.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors. The data are not publicly available due to ongoing research.
Acknowledgments: This study was supported by funds for the Nicolaus Copernicus University in Toruń, Poland.

Conflicts of Interest:
The authors declare no conflict of interest.