Ca 2+ /Na + Ratio as a Critical Marker for Field Evaluation of Saline-Alkaline Tolerance in Alfalfa ( Medicago sativa L.)

: Current indices of saline-alkaline (SA) tolerance are mainly based on the traditional growth and physiological indices for salinity tolerance and likely affect the accuracy of alfalfa tolerance predictions. We determined whether the inclusion of soil alkalinity-affected indices, particularly Ca 2+ , Mg 2+ , and their ratios to Na + in plants, based on the traditional method could improve the prediction accuracy of SA tolerance in alfalfa, determine important indices for SA tolerance, and identify suitable alfalfa cultivars in alkaline salt-affected soils. Fifty alfalfa cultivars were evaluated for their SA tolerance under SA and non-SA field conditions. The SA-tolerance coefficient (SATC) for each investigated index of the alfalfa shoot was calculated as the ratio of SA to non-SA field conditions, and the contribution of SATC under different growth and physiological indices to SA tolerance was quantified based on the inclusion/exclusion of special alkalinity-affected indices. The traditional method, excluding the special alkalinity-affected indices, explained nearly all of the variation in alfalfa SA tolerance, and the most important predictor was the SATC of stem length. The new method, which included these special alkalinity-affected indices, had similar explanatory power but instead identified the SATC of shoot Ca 2+ /Na + ratio, followed by that of stem length, as key markers for the field evaluation of SA tolerance. Ca 2+ , Mg 2+ , and their ratios to Na + hold promise for enhancing the robustness of SA-tolerance predictions in alfalfa. These results encourage further investigation into the involvement of Ca 2+ in such predictions in other plant species and soil types under more alkaline salt-affected conditions. The [SL], [SDM], [SS], [PRO], [K + ], [Na + ], [Ca 2+ ], [Mg 2+ ], [K + /Na + ], [Ca 2+ /Na + ], and [Mg 2+ /Na + ] represent the SATC of stem length, the SATC of shoot dry mass, the SATC of shoot soluble sugar content, the SATC of shoot proline content, the SATCs of shoot Na + , K + , Ca 2+ , Mg 2+ , K + /Na + ratio, Ca 2+ /Na + ratio, and Mg 2+ /Na + ratio, respectively. Abbreviations: SL, stem length; SDM, shoot dry mass; SS, shoot soluble sugar content; PRO, shoot proline content.

Given the growing demand for food, animal feed, and fiber to accommodate an expanding human population with limited new productive land, it will become necessary to utilize salt-affected lands [14,15]. As one of the most vital perennial legume forages cultivated worldwide, alfalfa (Medicago sativa L.) is a prominent feature of salt-affected lands, such as saline soil [16], mainly because of its good nutritional forage quality for livestock husbandry and its ability to fix nitrogen; however, it is considered moderately sensitive to salinity stress and alkalinity stress [17] Different from saline soils, where the focus is on ionic osmotic stress due to high salinity [18][19][20][21][22], alfalfa grown in such SA conditions often experiences alkaline stress, which inhibits plant growth and reduces yield [23,24]. Alfalfa is commonly at risk of ion toxicity, osmotic stress, and even high pH stress [25,26]. Therefore, it is essential to identify suitable alfalfa cultivars to effectively use salt-affected lands.
To accomplish this objective, a crucial prerequisite is distinguishing the most informative indices for the salt tolerance of alfalfa [27,28]. Traditionally, the methods for identifying these important indices for plant SA tolerance are often derived from salinity tolerance [29,30] (Table S1) [17,19,20,[23][24][25]28,. This is especially true for NaCl tolerance, in which the underlying mechanism for the maintenance of adequate K + in plant tissues under salt stress seems to depend upon selective K+ uptake and selective cellular K + and Na + compartmentation and distribution in the shoots [21,27,55] Consequently, the conventional indices for SA tolerance in alfalfa, including its main growth indices (such as dry biomass [17], root parameters [17,31,32], and relative growth rate [23] at the temporal scale of both days and weeks under controlled solutions and soil conditions) as well as its major physiological indices, such as photosynthesis components [23], antioxidant enzyme activity and malondialdehyde (MDA) [31], proline [17,31], soluble sugar [31], soluble protein [31], and cation content [23,24,[31][32][33], are often relied upon to express SA tolerance under controlled laboratory conditions. Recent reports have illustrated that important indices for conveying saline tolerance or SA tolerance can be identified by quantifying the relationships between plant physiological indices and biomass [28,30,56], or between some physiological indices and comprehensive D-values based on statistical tools, such as linear regression analysis [57]. Notably, the relative values of the investigated parameters under stressed conditions and nonstressed conditions have been emphasized due to their ability to better reflect the salt tolerance of plants by eliminating the inherent differences among cultivars [29,57].
Interestingly, increased Na + , Ca 2+ , and Mg 2+ accumulation in plant roots and an altered Ca 2+ and Mg 2+ balance in entire alfalfa plants were observed under controlled alkaline solution conditions, indicating the critical roles of plant Na + , Ca 2+ , and Mg 2+ in conveying SA tolerance in alfalfa under alkaline salt-affected conditions [24]. Additionally, a high soil pH environment can cause soil Ca 2+ and Mg 2+ to precipitate (CaCO3 and MgCO3), leaving fewer ions available for plant uptake [7,58]. Similarly, a high pH environment outside the roots can lead to greater absorption of Na + , thus competing with the uptake of other nutrient cations, including Ca 2+ and Mg 2+ [58]. Wang et al. (24) suggested that the contributions of Ca 2+ and Mg 2+ to the maintenance of plant ion balance and the response to SA stress should not be ignored under SA conditions. Thus, it is reasonable to infer that the roles of plant Ca 2+ and Mg 2+ in explaining the SA tolerance of alfalfa may be considerable under soil conditions vis-à-vis their traditional roles in terms of salinity tolerance and SA solution conditions (Table S1) [17,19,20,[23][24][25]28,. Considering that most research to date on SA tolerance in alfalfa has been conducted under SA solution conditions rather than SA soil conditions (Table S1) [17,19,20,[23][24][25]28,, the respective contributions of Ca 2+ and Mg 2+ to the SA tolerance of alfalfa should be recognized and investigated thoroughly. Crucially, soil-mediated SA stress is a continuous stress that acts upon plants throughout their entire vegetative growth period, rather than a suddenly imposed stress that occurs under controlled solution conditions [17,30,59]. Moreover, simulating such stresses in soil conditions is especially complex and difficult under field conditions when compared with the use of solutions [30,59]. To the best of our knowledge, few studies have quantified the contributions of plant nutrient cations, especially Ca 2+ and Mg 2+ , to SA tolerance in alfalfa under field SA soil conditions. Therefore, this study was conducted with two objectives: firstly, to test the hypothesis that the inclusion of soil alkalinity-affected indices, especially Ca 2+ , Mg 2+ , and their ratios to Na + in the plant, based on the traditional method could improve the prediction accuracy of SA tolerance in alfalfa. To obtain predictions using this new method and the traditional method, we quantified the relationships between SA tolerance in alfalfa and investigated the growth and physiological indices of the alfalfa shoot from field SA and field non-SA conditions. The second objective was to identify important indices for SA tolerance in alfalfa and suitable alfalfa cultivars in alkaline salt-affected soils. Taken together, the results of this study will help provide theoretical and practical guidance for the precise evaluation of SA tolerance in alfalfa crops and serve as a valuable genetic resource for developing SA-tolerant alfalfa cultivars in alkaline salt-affected soils.

Plant Growth and Experimental Setup
The field experiment was conducted at the Da'an Sodic Land Experimental Station (45°36' N, 123°53' E, and 132.1 m) in the Songnen Plain, northeast China ( Figure S1). The climate is semiarid and semihumid, averaging 413.7 mm of precipitation annually [13]. The soil characteristics can be found in Table 1. The soil samples were collected from the 0-10 cm soil layer before planting. Electrical conductivity (EC) and pH were determined under a 1:5 soil-to-water-extract ratio. Organic matter was determined with the potassium dichromate method. Available N, P, and K were determined by the industry standards LY/T1228-2015, LY/T1232-2015, and LY/T1236-1999, respectively.  The field experiment consisted of two soil types: a crop field, which served as the control ("NSA" = non-SA), and the SA-field, which served as the treatment. A split-plot design with three replicates was used with SA stress treatments as the main plots and cultivars as the subplots. The main plots consisted of the SA stress treatment and control (NSA). The subplots of 3 m 2 per plot consisted of 50 cultivars and were randomly placed in each main plot.
The seeds of each cultivar were sown at a density of 10 g per 3 m 2 on 1 June 2018. The cultivars used and their sources are given in Table S2. The alfalfa seeds were sown about 2 cm deep and covered with mixtures of aeolian soil and the local SA soil according to the local practice to prevent the formation of a soil crust and promote plant emergence even after irrigation [3,7]. After sowing, the plots were then irrigated using irrigation water with an EC of 1.05 mS cm −1 and pH of 7.52 obtained from an 80 m deep well, followed by weeding practices carried out as needed by hand.

Experimental Sampling and Measurements
Plant samples were measured and analyzed during the seedling stage of the cultivar at approximately five weeks after the seeds were sown in 2018. In 2019, plant samples were measured and analyzed during the seedling stage of the cultivar after regreening. This seedling stage was selected because it is the critical period for the early screening of alfalfa for SA tolerance [17,25,26,60]. Each cultivar contained three biological replicates. In each plot, two or more plants per cultivar were randomly harvested by digging up the entire plant, removing the plant at the shoot base, and then averaging the aboveground weight of one alfalfa plant. The stem length (SL) of alfalfa was measured as the length of the longest stem on the plant [61]. Shoot dry mass (SDM) was obtained by oven-drying the fresh samples at 105 °C for 30 min (for the first recording) and then again at 75 °C to a constant weight. The soluble sugar (SS) content of each sample was measured following the modified method described by Liu et al. [62]. A ca. 50 mg dry-powder sample of the shoot was extracted in 15 mL of double-distilled H2O and boiled for 20 min. Next, 5 mL of anthrone reagent was added to the 1 mL extract and incubated at 95 °C for 20 min and then cooled to room temperature. The absorbance of the solution was measured at 620 nm with a spectrophotometer (T6, Puxi General Instrument Co., ltd. Beijing China). The proline (PRO) content was determined according to the modified method by Bates et al. [63]. Specifically, ca. 50 mg of a dry-powder sample from the shoot was homogenized in 5 mL of 3% (w/v) sulfosalicylic acid, following which the mixture was placed in a boiling water bath for 10 min. The 18 mL plastic tubes were cooled and then centrifuged at 3000 × g for 20 min. Approximately 2 mL of the supernatant was mixed with 2 mL of glacial acetic acid and 2 mL of acidic ninhydrin reagent, and then boiled for 30 min. The plastic tubes were cooled, 4 mL luene was added to each tube, and the absorbance at 520 nm was determined with a spectrophotometer (T6, Puxi General Instrument Co., ltd. Beijing China).
Finally, ca. 0.1 g dry powder from each sample was digested with a mixture of HNO3 and HClO4 (v/v = 2:1) and diluted to 100 mL. Using this mixture, the total cation contents of Na + , K + , Ca 2+ , and Mg 2+ were determined by inductively coupled plasma mass spectrometry (ICPS-7500, Shimadzu Corporation, Japan). The K + /Na + ratio, Ca 2+ /Na + ratio, and Mg 2+ /Na + ratio indices in the alfalfa shoots were calculated using the corresponding values. In 2019, 21 cultivars survived, and only the aboveground stem length (SL) and shoot dry mass (SDM) as supporting evidences were measured during the seedling stage to verify the results in 2018.

Comprehensive Evaluation Formula of The Saline-Alkaline-Tolerance Coefficient of Each Index
The saline-alkaline-tolerance coefficient (SATC) for every plant trait was calculated using the following formula [29,57]: SATC = value of the trait under the SA condition / value of trait under the NSA condition (1) The eigenvector of each index in Table S3 was calculated based on eigen values and factor loadings values of different indices through a principal components analysis (PCA) using the following formula: For a given alfalfa cultivar, a comprehensive index (CI) was calculated based on the eigen values and factor scores of different alfalfa cultivars through a principal components analysis (PCA) ( Table S4) [57].

CI = Factor scores values × Eigen values
(3) For a given alfalfa cultivar, a comprehensive index (CI1-CI5) was calculated based on the SATC values of 11 indices through a principal component analysis (PCA) [57]. Based on the CI1-CI5 indices, subordinate function values μ(Xj) for the different alfalfa cultivars were calculated using the following formula (Table S4): where Xj is the j comprehensive index, Xmin is the minimum value of the j comprehensive index, and Xmax denotes the maximum value of the j comprehensive index. The weight (Wj) of the comprehensive indexes (CIs) for the different alfalfa cultivars was calculated as follows (Table S4): where Wj indicates the importance of the jth comprehensive index among five comprehensive indices and Pj is the contribution rate of the jth comprehensive index of each alfalfa cultivar (obtained from the PCA). The comprehensive value (D-value) of the SA tolerance for each alfalfa cultivar was then calculated as follows (Table S4):

Statistical Analysis
PCA is a method that can transform numerous indexes into a small number of comprehensive indexes under the premise of less information loss, and effectively concentrate data and simplify the index, so as to make up for the deficiency of SA-tolerance evaluation of individual indexes. The main operational steps of PCA are summarized as follows. Firstly, the data were normalized to analyze descriptive statistical variables. Secondly, the eigen values, or variances of each principal component (PC), were extracted using dimension-reduction factor analysis. Based on standardized data, the correlation matrix was used to analyze and extract the eigen values. Thirdly, the eigenvector of each index was calculated based on the component matrix using Equation (2) and the comprehensive index (CI) was calculated based on the eigen values and factor scores using Equation (3). PCA was performed by SPSS v20.0 software (SPSS, Chicago, IL, USA). The differences in the investigated parameters of the 50 alfalfa cultivars between the NSA and SA treatments were compared using a paired t-test (also performed in SPSS v20.0). The "complete" cluster method is used for distance measuring the dataset of D values using the R statistical platform v3.5.1 (R Development Core Team, Vienna, Austria).
Stepwise forward regression was used to quantify the explanatory power of the different investigated plant growth and physiological variables for predicting the SA tolerance of alfalfa and to identify the key variables predicting SA tolerance. The SA tolerance of the 50 alfalfa cultivars was the dependent variable, with the SATCs of the plant growth and physiological variables set as the independent variables.
In this study, a PCA of SATCs was firstly conducted and then combined with subordinate function analysis to obtain the comprehensive D-value, following which the SA tolerance of the cultivars was classified according to the D-value by cluster analysis. Additionally, the stepwise regression relationship between the D-value and each growth and physiological index was established to determine the most suitable index.
To convey the SA tolerance of alfalfa, the traditional method was adopted based on the inclusion of SATCs of SL, SDM, SS, PRO, Na + , K + , and K + /Na + ratio, considering that most studies of alfalfa have frequently selected these variables as being characteristic for salinity tolerance, including SL, SDM, SS, PRO, Na + , K + , and K + /Na + ratio (literature reviewed in Table S1). For the new method proposed here, frequently neglected variables, especially Ca 2+ , Mg 2+ , Ca 2+ /Na + ratio, and Mg 2+ /Na + ratio in terms of the SATCs, were further included, based on the inclusion of the SATCs of SL, SDM, SS, PRO, Na + , K + , and K + /Na + ratio, which were adopted by the traditional method.

Genetic Variation in Plant Growth and Physiological Indices under SA Treatment vs. NSA Treatment
SA had a significant effect on both SL and SDM (P < 0.05; Table 2). The SATCs for SL and SDM ( Figure 1    SA stress had a significant effect on SS (P < 0.05) and PRO (P < 0.05; Table 2). Both SS and PRO increased in all cultivars under SA conditions, but the tolerance of the cultivars differed as did their respective accumulation amounts ( Figure 2). For SS, the cultivar WL343HQ had the highest SATC value of 4.39, whereas Algonguin had the lowest SATC value of 1.11 ( Figure 2A). PRO was highest in cultivar WL319HQ with a SATC value of 9.73 and lowest in Gannong NO.4, which had a SATC value of 0.81 ( Figure 2B). SA stress also significantly affected the contents of Na + , K + , Ca 2+ , Mg 2+ , and the ratios of Ca 2+ /Na + and Mg 2+ /Na + (P < 0.05, Table 2). The SATC of K + ranged from 0.68 in cultivar WL354HQ to 2.46 in cultivar Bingchi ( Figure 3A). For Na + content, the cultivar Vison had the highest SATC value of 8.55. By contrast, WL319 had the lowest SATC value of Na + under SA conditions, with a value of 0.86 ( Figure 3B). For Ca 2+ content, the cultivar WL903HQ had the highest SATC value of 0.82, whereas this was lowest in Eureka at 0.37 ( Figure 3C). The SATC value of Mg 2+ was lowest in cultivar WL366HQ at 0.56 and highest at 1.03 in Eureka ( Figure 3D). The SATC of the K + /Na + ratio ranged from 0.29 in WL354HQ to 3.97 in Bingchi ( Figure 3E). For the Ca 2+ /Na + ratio, the cultivar WL319HQ had the highest SATC value (0.81), whereas this was lowest in the cultivar Vison ( Figure 3F). Finally, for the Mg 2+ /Na + ratio, the SATC values were highest in cultivar WL319 at 0.97 and lowest in the cultivar Vison ( Figure 3G).

Classification of the SA-Tolerance of Alfalfa Based on the Different Methods
Following cluster analysis of the comprehensive D-value, the 50 cultivars were divided into four categories using the traditional method ( Figure 4)  In 2019, the surviving 21 cultivars were observed along with their autumn dormancy grades ≤ 3.2 (Table S2). Under the traditional method, the average SATC of SL in the sensitive SA-tolerant category was significantly lower than that in the other categories, but there was no significant difference in the SATC of SDM among the different categories ( Figures 5A,B). Under the new method, the average SATCs of SL and SDM demonstrated the same trend: strong SA-tolerant cultivars were significantly higher than weak SA-tolerant cultivars (Figures 5C,D).

Identifying the Most Important Indices for Evaluating SA Tolerance in Alfalfa Based on Different Methods.
The relationship between each index and D in the two methods was determined by stepwise regression analysis, and the index with the highest R 2 value was added first to the regression model (Table 3). Under the traditional method, the SATC of SL explained the most variation (67.90%) in the D response, and the SATC of PRO was the second-most important index with an explanatory power of 13.60% for D, followed by the SATCs of K + /Na + ratio, SDM, SS, K + , and Na + (Table 3). However, under the new method, it was the SATC of the Ca 2+ /Na + ratio that was the best predictor, explaining 62.50% of the variation in D. The SATC of SL was the second variable to be included, followed by the SATCs of Mg 2+ , Ca 2+ , SDM, PRO, Na + , SS, and K + (Table 3).

Assessing the Appropriateness of SA-Tolerance Prediction and the Identification of Suitable Cultivars for saline-alkaline (SA)-Tolerance
The traditional method explained 99.9% of the variation in alfalfa SA tolerance (Table 3), as did the newly proposed method, in which shoot Ca 2+ , Mg 2+ , and their ratios to Na + were included. The high explanatory power of both methods may be ascribed to the increased accuracy in quantifying the stress tolerance of the plant by simultaneously analyzing multiple parameters based on PCA (Table S3). Likewise, the relative values of the various plant growth and physiological parameters under SA and control conditions used in this study may reflect the tolerance of this plant more precisely when compared with their absolute values [29,30,57]. It thus seems feasible to use the plant growth and physiological variables based on the traditional salinity-characterized method as well as the new method. In fact, the D-value, as a comprehensive value of SA tolerance, changed when different parameters were included since they were generated from the reintegration of the included parameters and their assigned weights (Table S4). Specifically, the D-value of the same alfalfa cultivar for a certain plant trait obtained via the traditional method differed from that using the new method. In this respect, the most important variables identified by the PCA and linear regression analysis could only reflect the SA tolerance based on the included parameters. Further, the traditional method identified the SATC of SL as being the strongest index for SA tolerance (Table 3), while for the new method, which included shoot Ca 2+ , Mg 2+ , and their ratios to Na + , the strongest index was the SATC of shoot Ca 2+ /Na + ratio. It appears that the new method might have an improved ability to identify the most important index for predicting the SA tolerance of alfalfa.
By further analyzing the categories of the alfalfa cultivars based on the two methods (Figure 4), some known alfalfa cultivars, especially WL903HQ, WL298HQ, Victory, Gannong NO.1, and Alfaqueen, classified by the new method into the different SA-tolerance categories were observed to be more consistent with the reported categories when compared with those obtained by the traditional method ( Figure 4). According to previous studies that screened the saline or salinealkaline tolerance of alfalfa, WL903HQ [64,65], WL298HQ [66], Victory [67], and Gannong NO.1 [68] were also considered as being relatively highly tolerant to saline or saline-alkaline stress conditions, which is more consistent with the classification results of the new method. Moreover, in recent studies, Alfaqueen was identified as being weakly SA-tolerant through laboratory-controlled experiments [64,67]. These results are also more consistent with the classification results of the new method, considering that the cultivar was grouped into the category of a moderately SA-tolerant cultivar using the traditional method. We suggest that the new method is an improvement over the traditional method when aiming to classify all alfalfa cultivars.
Since the autumn dormancy grade is an important factor in cultivar selection [69], the surviving alfalfa cultivars in 2019 were used to partly verify the SA-tolerance classification using the two methods ( Figure 5). As expected in the SATC of SL, there were significant differences in the SATC of SL, with higher values in the strongly SA-tolerant category in both methods ( Figure 5A,C). However, it is noteworthy that the average SATC of SDM in the strongly SA-tolerant category was higher than that in the other categories under the new method, but there was no significant difference in the different categories of alfalfa under the traditional method ( Figure 5B,D). Many studies have demonstrated that SDM is a crucial index for expressing the SA tolerance of plants, as the more tolerant alfalfa cultivars have higher SDM under SA conditions when compared with their less tolerant counterparts [56,70]. Thus, it might be concluded that the classification of alfalfa SA tolerance based on the new method is more appropriate for discriminating cultivar properties, such as the SATC of SDM, compared to the traditional method.
Taken together, the new method holds potential for enhancing the SA-tolerance evaluation of alfalfa compared with the traditional method. However, further investigation of the suitability of the new method based on the inclusion of alfalfa shoot Ca 2+ , Mg 2+ , and their ratios to Na + will require testing under a gradient of alkaline salt-affected soil and climatic conditions.

Identification of the Best Markers for Predicting the saline-alkaline (SA)-Tolerance of Alfalfa
It is well known that the success of indirect selection for stress tolerance by using plant attributes as markers depends on the strength of the relationship of such markers with plant phenotypic responses to stress factors [27,30]. Here, we identified that the most important variable conveying SA tolerance was not the same when differing plant growth and physiological variables were included. The traditional method identified the SATC of SL as having the greatest single explanatory power at 67.90% (Table 3). Likewise, the new method also identified the SATC of SL as a strong index for saline-alkaline (SA)-tolerance (Table 3). Hence, SL likely plays an important role in governing and predicting SA-tolerance. It is worth noting that SL is the length of the longest stem on the alfalfa plant rather than the natural height. Actually, SL was once considered important in evaluating alfalfa yield [71,72], and in recent years it was reconsidered as a vital index in evaluating alfalfa yield and salinity tolerance [33,61]. From a practical and economic perspective, stem length should be considered as a chief screening criterion under SA field conditions, since its measurement is easy and fast.
When the shoot Ca 2+ , Mg 2+ , and their ratios to Na + were included in the new method, the SATC of Ca 2+ /Na + ratio was identified as the strongest index, having a single maximal explanatory power of 62.50% (Table 3). This suggests that Ca 2+ /Na + ratio could play a major role in predicting the SAtolerance of alfalfa growing under realistic SA field conditions as compared with SL when both are simultaneously evaluated. Calcium is crucial for maintaining the structural and functional integrity of plant membranes, in addition to its function in the regulation of ion transport, selectivity, ion homeostasis maintenance, signal transduction, and activation of cell wall enzymes and antioxidant enzyme activities [19,20,34,73]. According to Suarez and Grieve (1988) [74], the Ca status of the plant shoot is strongly influenced by the ionic composition of the external medium; hence, the Ca 2+ /Na + ratio may be predicted from the solution composition of the root media and shoot ion concentrations under salinity-stressed conditions. Meanwhile, some studies have shown that the maintenance of Ca acquisition and transport under salt stress conditions is a prominent determinant of salinity tolerance in plants [75]. Even under saline conditions, the Ca 2+ /Na + ratio was once proposed as a promising index for gauging plant salt tolerance [33,35]. However, the physiochemical properties of SA soils, especially the ionic composition and pH, frequently differ from those of saline soils [3,7,11], which probably leads to disparate Ca status and Ca 2+ /Na + ratios between them. In such cases, it is necessary to determine whether the shoot Ca 2+ /Na + ratio is still relevant as an index for the SA tolerance of plants, including alfalfa. Unfortunately, under alkaline salt-affected conditions, the contributions of Ca 2+ towards the maintenance of ion balance and response to SA stress are frequently ignored [36]. To the best of our knowledge, this study is the first to have quantified the contributions of shoot Ca 2+ /Na + ratio and Ca 2+ , as well as conventional plant and physiological indices in terms of relative values, to identify their roles in the SA tolerance of alfalfa in the field under alkaline salt-affected conditions. Our results highlight the importance of Ca 2+ /Na + ratio in the SA tolerance of alfalfa under alkaline salt-affected conditions.
There are practical reasons for using the shoot Ca 2+ /Na + ratio of alfalfa in terms of the relative values under SA and NSA conditions to predict SA tolerance. First, these data are easily available since the measurement of Ca 2+ and Na + in plants is relatively simple to conduct in the laboratory. Second, including the Ca 2+ /Na + ratio could improve SA-tolerance evaluations of alfalfa. We speculate that the Ca 2+ /Na + ratio could be significant for distinguishing the SA tolerance of more plant genotypes or species in similar alkaline salt-affected conditions. The utility of Ca 2+ /Na + ratio for predicting the SA tolerance of alfalfa may be further enhanced if its appropriateness in future studies is extensively evaluated by quantifying the relative contributions of Ca 2+ and Ca 2+ /Na + ratio to SA tolerance, which could be achieved by testing multiple plant species across soil types in alkaline saltaffected conditions.
It is well known that Mg 2+ is essential for protein synthesis and chlorophyll structure, being an activator for many photosynthetic and respiratory enzymes [37]. More importantly, the accumulation of cations, including Mg 2+ , as important nutrients, might also be an important mechanism in the plant response to SA stress by maintaining ion homeostasis [38,76] and contributing to signaling and enzyme activity [34]. Consistently, we found that the SATC of shoot Mg 2+ /Na + was significantly and positively related to the D-value (r = 0.74, P < 0.01) (Table S5). These results provide supporting evidence for the use of the SATC of shoot Mg 2+ /Na + as a robust predictor influencing the SA tolerance of alfalfa. However, the relationship between the SATC of the shoot Mg 2+ and the SA tolerance (D) was relatively weak (Table S5), likely due to the lower soil Mg 2+ content in contrast to soil Ca 2+ and Na + (Table 1). Additionally, the inhibition of Mg 2+ uptake by Na + and Ca 2+ occurs under SA conditions [37,39]. Furthermore, the SATC of the shoot Mg 2+ /Na + ratio was significantly positively correlated with that of shoot Ca 2+ /Na + ratio (r = 0.93, P < 0.01), indicating the importance of the SATC of the shoot Mg 2+ /Na + ratio in predicting the SA tolerance of alfalfa. This result corroborates the earlier findings of Bernstein et al. [77], who argued that the Mg 2+ /Na + ratio was a critical index of the saline tolerance of sorghum. However, because it was wholly or partly displaced by the SATC of Ca 2+ /Na + ratio and Ca 2+ in the linear regression analysis (Table 3), this variable did not emerge as an important predictor in the present study. Meanwhile, although the vital role of Mg 2+ in conveying the SA tolerance of plants is often ignored [36], we still speculate that Mg 2+ and Mg 2+ /Na + ratio have non-negligible roles in determining the SA tolerance of alfalfa in SA field conditions.
Agronomic characters represent the combined genetic and environmental effects upon plant growth and yield, as they include the physiological mechanisms conferring salinity and alkalinity tolerance [27,34,40,78]. When using the two methods, the SATC of SDM was separately identified as an important predictor and was significantly related to the SL (Table S5). These results indicate that SDM may contribute to SA tolerance, which supports previous findings [56,70]. Interestingly, except for SL, the SATC of SDM was not significantly related to any other investigated independent variable, which is similar to the results reported by Vadez et al. [79]. In our study, the soil EC before the trial was 347.22 uS cm −1 , which was equal to the EC of saturated paste extract with a value of 3.78 mS cm -1 according to Chi and Wang [4]; this is higher than the salinity threshold of alfalfa [80]. Despite the synergistic interaction between pH and salinity on the growth and physiological properties of alfalfa [32,81], stresses, including soil salinity and the pH (9.09) under the SA field conditions, are likely insufficient to satisfy the SA intensity and stress duration [30,32,82,83]. In such cases, the significant relationships between SDM and the investigated variables in the present SA conditions likely did not exist, in contrast to the stronger intensity previously reported for SA solutions-based laboratories [30,32].
Additionally, the SATCs of shoot proline and SS were both identified as predictors of SAtolerance using both methods, suggesting their involvement in SA-tolerance, though their contributions were relatively low compared with the other variables (Table 3). Debate persists on proline and SS as indicators of plant stress, such as salinity or alkalinity. Some studies deemed PRO as an important index [41], whereas others have argued against this interpretation [66,68]. Soluble sugar, as an osmotic regulator under stress, depends not only on plant species, but also on the cultivar, salt level, duration of the exposure to stress, developmental stage, and experimental conditions. It may increase [84], decrease, or even stay constant [28,85]. Thus, the functional relevance of proline and sugar contents for predicting the SA tolerance of alfalfa merits further testing and verification.
Other factors that were highly ranked included the SATCs of shoot Na + and K + as predictors of the SA tolerance of alfalfa (Table 3). This result is consistent with those obtained using the traditional method, as indicated by an adaptive mechanism wherein the absorption of Na + increases with the competing uptake of nutrient cations such as K + under SA conditions [32,36]. Interestingly, using the new method, the respective predictive ability of SA tolerance was lower than that of Ca 2+ (Table 3). This finding suggests a minor role of shoot Na + or K + in SA tolerance under SA field conditions at our site, probably because the soil available K under SA field conditions is relatively high and adequate for uptake by alfalfa roots (Table 1) [86]. Consequently, the SATC of shoot K + /Na + might not serve as a reliable marker for plant tolerance to SA under field conditions. Despite this, the K + /Na + ratio is still widely employed as a robust marker for salinity tolerance [22,[87][88][89][90] and SA tolerance (Table S1). Practically, the initial soil cation nutrients, including K, under SA field conditions were obviously different from the initial solution nutrients in the Hoagland's solution (Table S6.1) [23,30,32,37,39], as indicated by the cation contents in the alfalfa plants under nonstressed treatments between controlled solution conditions and the present field soil conditions in Tables S6.1 and 6.2. In such cases, alfalfa plants were supplied with initial contrasting contents of important cation nutrients, especially K and Ca, under SA field conditions compared with solution conditions (Tables S6.1-6.3). This could possibly account for the different ratios of important cations, such as K + /Na + ratio, Ca 2+ /Na + ratio, or Mg 2+ /Na + ratio, in the alfalfa when subjected to salinity and alkalinity stresses [21,22,24,31,32,74,77]. In this sense, our present study helps enhance our mechanistic understanding of plant-soil interactions by quantifying the relationship between alfalfa tolerance and plant growth, and physiological traits based on the soil systems in a more realistic and practical manner compared with the solution system. Future study is needed to quantify the roles of initial soil cation nutrients under SA field conditions in determining the most important factors for assessing the SA tolerance of alfalfa from the perspective of initial soil cation nutrients under SA field conditions.
Taken together, we infer that the SATCs of shoot Ca 2+ /Na + ratio, stem length, and Mg 2+ are the best indices for predicting the SA tolerance of alfalfa, as they were entered first into the regression analysis and SATCs of other investigated variables and were significantly related to the identified variables under SA conditions (Table 3). Furthermore, significant differences in SDM, SL, shoot Ca 2+, and shoot Ca 2+ /Na + ratio were detected using both methods (Table S6.3) when further analyzing the differences in the investigated variables among different categories of alfalfa under SA field conditions. This result is consistent with Wang et al. [24], who found that a salt-tolerant alfalfa cultivar had a higher SDM, shoot Ca 2+ , and Ca 2+ /Na + ratio compared with a low salt-tolerant cultivar. Coupled with the significant relationship identified between SL and SDM (Table S5), these results provide supporting evidence for the selection of the SATC of shoot Ca 2+ /Na + ratio and SATC of SL as important markers for predicting the SA tolerance of alfalfa crops.
Salt reduces plant growth via several distinct processes [22]. In the first osmotic phase there is an independent effect of shoot salt accumulation, in which the growth of shoots significantly slows due to the salt outside the roots, for which there is surprisingly little genotypic variation [21,22]. Next, in the ionic phase of salt toxicity, which requires time to progress and results from internal injury, either species or genotypes can utilize two strategies of ion efflux and tissue tolerance to resist ionic stress [19,21,22], likely explaining their divergent ability for tolerating SA soil conditions. Hence, the tolerance mechanisms and the temporal scale on which they operate may be critical for controlling growth at different time periods under exposure to saline-alkaline stress conditions; thus, these mechanisms should be considered. In our study, plant samples were collected only at the seedling stage. Several studies have shown that using this stage to evaluate the responses of alfalfa to salinity and alkalinity is important for elucidating its tolerance mechanisms and identifying important makers and suitable cultivars with different abilities to tolerate saline or SA conditions [17,20,25,26,60]. However, the differences in the levels of salt tolerance at the seedling stage did not completely reflect their salinity tolerance as adults [19,27,91]. Thus, the suitability of markers for SA tolerance, especially the Ca 2+ /Na + ratio, in distinguishing suitable alfalfa cultivars may be enhanced if the tolerance of alfalfa to other SA conditions is further studied at the whole-plant level through its ontogeny [27,91].

Conclusions
The results revealed that the inclusion of shoot Ca 2+ , Mg 2+ , and their ratios to Na + holds potential for enhancing SA-tolerance evaluations of alfalfa plants. The SATCs of shoot Ca 2+ /Na + ratio and stem length are key markers for predicting the SA tolerance of alfalfa based on the SA field conditions in this study. We speculate that Ca 2+ and the Ca 2+ /Na + ratio could play important roles in reflecting the SA tolerance of alfalfa growing in alkaline salt-affected soil conditions. The 50 alfalfa cultivars were classified into strong SA-tolerant cultivars with maximum shoot Ca 2+ /Na + ratios, stem lengths, SDM and Ca 2+ , as well as moderate SA-tolerant and sensitive SA-tolerant cultivars with minimum corresponding values under SA field conditions. These findings encourage investigations of how Ca 2+ participates in saline-alkaline tolerance and its predictions across multiple plant species, climates, and soil types in alkaline salt-affected conditions.