Transcriptomic Analysis Reveals Salt Tolerance Mechanisms Present in Date-Plum Persimmon Rootstock ( Diospyros lotus L.)

: Agriculture needs solutions for adapting crops to increasing salinity globally. Research on physiological and molecular responses activated by salinity is needed to elucidate mechanisms of salinity tolerance. Transcriptome proﬁling (RNA-Seq) is a powerful tool to study the transcriptomic proﬁle of genotypes under stress conditions. Persimmon species have di ﬀ erent levels of tolerance to salinity, this variability may provide knowledge on persimmon species and development of salt–tolerant rootstocks. In this study, we conducted a physiological and transcriptomic proﬁling of roots and leaves in tolerant and sensitive plants of persimmon rootstock grown under saline and control conditions. Characterization of physiological responses along with gene expression changes in roots and leaves allowed the identiﬁcation of several salt tolerance mechanisms related to ion transport and thermospermine synthesis. Di ﬀ erences were observed in putative H +/ ATPases that allow transmembrane ionic transport and chloride channel protein-like genes. Furthermore, an overexpression of thermospermine synthase found in the roots of tolerant plants may indicate that alterations in root architecture could act as an additional mechanism of response to salt stress. These results indicate that Diospyros lotus L. exhibits genetically-controlled variability for salt tolerance traits which opens potential opportunities for breeding salt-tolerant persimmon rootstocks in a Mediterranean environment challenged by drought and salinity. M.R.G.-G.; resources, M.L.B.; data curation, N.D., A.Q., and F.G.-M.; writing—original draft preparation, F.G.-M. and M.L.B.; writing—review and editing, M.L.B., N.D., and M.R.G.-G.; supervision, M.L.B., N.D., and M.R.G.-G., project administration, M.L.B.; acquisition, M.L.B.


Introduction
Drought and salinity are two of the main challenges in agriculture. Reduction of water availability and increase of salinity in soils and water reservoirs can limit crop production due to severe inhibition of plant growth and development via osmotic and ionic stresses [1,2]. The negative effect of salinity in photosynthesis is caused either by a reduction in available water, by the toxic accumulation of ions in the cell, or both. Furthermore an attenuated photosynthesis leads to an increase in reactive oxygen species (ROS) and a consequent oxidative damage of the cells [3]. The cultivated surface affected by either stress increases year after year as a consequence of more frequent drought episodes associated to climate change. Among the affected areas, the Mediterranean basin registered a significant increase in salinity in soils and water reservoirs in the last 10 years [4]. Therefore, there is a pressing need to find solutions to alleviate the negative impact on plant production in an area dominated by fruit tree crops. physiological and RNA-Seq analyses of roots and leaves in persimmon rootstocks (Diospyros lotus L.) grown under saline and control conditions.

Plant Material and RNA Extraction
A full-sibling family of Diospyros lotus L. (DL) obtained from a cross between two D. lotus trees located at the Instituto Valenciano de Investigaciones Agrarias (IVIA) persimmon germplasm bank was used in this study. Seeds were stratified in plastic bags filled with perlite in a cold chamber at 4 • C at the end of March. After 30 days of stratification, seeds were transferred to peat-moss and perlite (4:1 ratio) and transferred to a greenhouse at 18-24 • C for two months (from 24 April to 26 June 2017).
One hundred and fifty seedlings were transplanted into 1L pots filled with coarse sand. Plants were acclimatized in the greenhouse until apical meristem growth was observed. During this period, the plants were watered with a nutrient solution (3% Cristaljisa 18-18-18, soluble fertilizer with micronutrients). The plants were distributed randomly in the greenhouse to avoid possible biases. The plants were acclimated before exposition to salinity treatment. After acclimation, plants were exposed to saline stress treatment for 60 days. The saline stress consisted of adding 40 mM NaCl to the nutrient solution. A total of 25 control plants remained watered with only nutrient solution.

Plant Phenotyping
Plant phenotyping was carried out using the same measurement methods from a previous work [38]. The morphological parameters total height (cm), leaves (no.), nodes (no.), internodes (cm), and defoliation (1-no. leaves/no. nodes) were measured at days 0, 30, 45, and 60. Relative growth rate (RGR) was calculated as shown in the Equation (1) Treated plants were divided into tolerant, intermediate, and sensitive phenotypes according to the strong differences shown in the phenotype after the saline stress. Further analyses were carried out in the tolerant, sensitive, and control plant groups. The salt stress tolerance of tolerant and sensitive plants was further confirmed by analysis of physiological parameters and principal component clustering.
Stem water potential (ψH, MPa) was measured in fully expanded leaves using a Model 600 Schölander Pressure Chamber (PMS Instrument Company, Albany, OR, USA). Measured leaves were kept in reflective plastic bags for 30 min (Levin, 2019). Measurements of leaf net CO 2 assimilation rate (A CO2 ), substomatal CO 2 concentration (C i ), and stomatal conductance (g s ) were measured on three single leaves of each plant. Determinations were performed in fully expanded leaves from the mid-stem zone. The measurements were made using a CIRAS-2 Portable Photosynthesis System (PP Systems, Amesbury, MA, USA) using a PLC 6 (U) universal leaf autocuvette in a closed-circuit model and kept at 25 ± 0.5 • C, with a leaf-to-air vapor deficit of about 1.7 kPa. The photosynthetically active radiation (PAR) was adjusted to 1000 µmol m −2 s −1 photon flux density and the air flow rate through the cuvette was 0.5-1.5 L min −1 . The measurements were made on a sunny day between 9:30 a.m. and 12:30 p.m. at the end of the salt treatment (day 60).
After saline treatment, three leaves from each plant were collected. Sample pre-treatments were performed as described previously [38]. Na + , Ca 2+ , K + , Mg 2+ , P, and S content was determined using a multiple-collector inductively coupled plasma mass spectrometry (MC-ICP MS, Thermo Finnigan Neptune). Chloride concentration was quantified by silver ion-titration [50] with a Corning 926 automatic chloridometer (Corning Ltd., Halstead Essex, UK).
Statistical differences calculation was done using RStudio v1. 1.447 (2018) with packages from the Comprehensive R Archive Network (CRAN). Due to the small number of samples within groups (<20), normality cannot be assumed, therefore non-parametric tests were used for average comparison. Parameters were statistically tested by Kruskal-Wallis test (p ≤ 0.05) averages between the three groups were compared with the Pairwise Wilcoxon-Mann-Whitney test at 95% confidence level (p ≤ 0.05). Principal component analyses were calculated using Statgraphics Centurion, 16.1 version (Statistical Graphics, Englewood Cliffs, NJ, USA) with centering and standardization.

RNA Extraction
After 60 days of salt treatment, young fully expanded leaves and root tip tissue were collected and immediately frozen and powdered using liquid nitrogen. Control samples from all populations were also collected and processed. Roots and leaves from three of the plants from control treatment, and six treated plants (three tolerant, three sensitive) seedlings were collected for RNA extraction.
The samples of each of the three individuals were mixed. The three biological replicates were made resampling and mixing new tissue samples. RNA was isolated according to Gambino et al. [51]. DNA was removed with the RNase-Free DNase Set (Qiagen, Valencia, CA, USA), using the RNeasy Plant Mini Kit (Qiagen). After extraction, RNA quantity and integrity were checked with electrophoresis and Qubit fluorometer (Invitrogen, Carlsbad, CA, USA) quantification. Purified RNA (500 ng) was reverse transcribed with PrimeScript RT Reagent Kit (Takara Bio, Otsu, Japan) in a total volume of 10 µL.

RNA Sequencing
RNA samples were sequenced by Novogene ® using the following procedure: mRNA was enriched using oligo (dT) beads and fragmented randomly using a fragmentation buffer. cDNA was synthetized using random hexamers and reverse transcriptase. A second strand was generated using nick-translation with RNAse H and Escherichia coli Polymerase I with Illumina second strand synthesis buffer and dNTPs. AMPure XP beads were used to purify the cDNA and the first strand was degraded with USER enzyme. Then, the resulting cDNA was purified, terminal repaired, A-tailed, ligated to the sequencing adapters, size-selected, and enriched with PCR. The resulting library was sequenced using Illumina.

Plant Physiological Responses Associated to Salt Stress Tolerance
Plants exposed to saline treatment showed significant reduction in growth, increased leaf damage and defoliation. The saline effect was very critical in these variables comparing to the control treatment, as shown by the PCA analysis performed to the morphological phenotyping at 0, 30, and 45 days ( Figure 1). The first two components of the analysis explained more than 60% of the variability and a total of 5 principal components were identified with eigenvalue >1 (Supplementary Figure S1). Among the treated plants, those that showed severe symptoms of salinity were categorized as sensitive, whereas the plants that showed less severity to the salt stress were categorized as tolerant (Supplementary Figure S2). The morphological PCA showed a separation between tolerant and sensitive groups, mainly in the Component 1. This component presents a great influence of growth and relative height parameters ( Figure 1). These differences in morphological parameters were statistically confirmed between groups (Table 1). In addition, the tolerant plants showed higher relative growth rate (RGR) than the sensitive ones. Furthermore, differences between tolerant and sensitive plants were significant for the relative growth rate (RGR) (Figure 2). The sensitive plants showed less RGR than the tolerant ones at the end of the saline treatment (day 60).
Physiological responses and morphological phenotype between the groups showed greater differences compared to the morphological traits measured, showing a clear separation between the tolerant and sensitive plants as shown in the PCA analysis ( Figure 3). The first two components of the analysis explained more than 70% of the variability and a total of 4 principal components were identified with eigenvalue >1 (Supplementary Figure S3). The main differences between groups were both found in the first and second components of the analysis.
Differences have been observed related to plants hydric potential. Sensitive plants showed significant higher water potential compared to the control and tolerant plants. In the case of gas exchange parameters, sensitive plants showed less photosynthesis than both control and tolerants, altogether with higher substomatal CO 2 concentration. In addition, sensitive plants showed reduced transpiration compared to the control. Regarding to ion contents, as compared to tolerant plants, the sensitive ones showed a significantly higher Cl − concentration. Interestingly, no significant differences were observed for Na + accumulation, neither were differences significant for Na + /K + and Na + /Ca 2+ ratios. Moreover, P accumulation was more pronounced in tolerant plants. The overall reduction in carbon fixation capacity and water potential due to salinity were significantly more pronounced in the sensitive plants. (Table 2) Agronomy 2020, 10, 1703 6 of 19 Agronomy 2020, 10, x FOR PEER REVIEW 6 of 19 Furthermore, differences between tolerant and sensitive plants were significant for the relative growth rate (RGR) (Figure 2). The sensitive plants showed less RGR than the tolerant ones at the end of the saline treatment (day 60).   Physiological responses and morphological phenotype between the groups showed greater differences compared to the morphological traits measured, showing a clear separation between the tolerant and sensitive plants as shown in the PCA analysis ( Figure 3). The first two components of the analysis explained more than 70% of the variability and a total of 4 principal components were identified with eigenvalue >1 (Supplementary Figure S3). The main differences between groups were both found in the first and second components of the analysis.
Differences have been observed related to plants hydric potential. Sensitive plants showed significant higher water potential compared to the control and tolerant plants. In the case of gas exchange parameters, sensitive plants showed less photosynthesis than both control and tolerants, altogether with higher substomatal CO2 concentration. In addition, sensitive plants showed reduced transpiration compared to the control. Regarding to ion contents, as compared to tolerant plants, the sensitive ones showed a significantly higher Cl − concentration. Interestingly, no significant differences were observed for Na + accumulation, neither were differences significant for Na + /K + and Na + /Ca 2+ ratios. Moreover, P accumulation was more pronounced in tolerant plants. The overall reduction in carbon fixation capacity and water potential due to salinity were significantly more pronounced in the sensitive plants. (Table 2) Table 2. Physiological and nutritional phenotype of the plants between groups after 0, 30, 45, and 60 days. Different letters represent significant differences (p < 0.05) between groups.

Variable
Tolerant

Gene Expression Changes in Response to Salinity
The plants chosen for transcriptomic sequencing were selected with the confirmation of the tolerance by the physiological responses after the saline stress ( Figure 3). To visualize the gene expression profiling results, a principal component analysis (PCA) was conducted using the transcriptomic data for all samples. In this analysis, a clear separation was observed between roots and leaves. Leaves showed little variability between the three categories ( Figure 4). However, the percentage of variability explained between control and treated roots was of the same order of magnitude as that between leaves and roots. Furthermore, the separation between tolerant and sensitive plants was also significant based on root transcriptomic data as shown in Figure 4. This separation was consistent with the phenotypic data presented above.

Gene Expression Changes in Response to Salinity
The plants chosen for transcriptomic sequencing were selected with the confirmation of the tolerance by the physiological responses after the saline stress ( Figure 3). To visualize the gene expression profiling results, a principal component analysis (PCA) was conducted using the transcriptomic data for all samples. In this analysis, a clear separation was observed between roots and leaves. Leaves showed little variability between the three categories ( Figure 4). However, the percentage of variability explained between control and treated roots was of the same order of magnitude as that between leaves and roots. Furthermore, the separation between tolerant and sensitive plants was also significant based on root transcriptomic data as shown in Figure 4. This separation was consistent with the phenotypic data presented above.  The comparison between expression in tolerant vs. sensitive plants revealed 2901 differentially expressed genes. These genes were further filtered using the BLAST results and annotations to look for putative genes related to salinity tolerance. Additionally, an orthogonal design analysis was conducted to be able to differentiate gene expression differences specific to roots or to leaves. The orthogonal analysis identified 1277 genes differentially expressed between tolerant and sensitive leaves (Supplementary Figure S4). In addition, 1395 genes were differentially expressed between tolerant vs. sensitive roots (Supplementary Figure S5).
Differentially expressed genes were divided according to several salt tolerance mechanisms: ion transport, photosynthesis and respiration systems, and ROS detoxification. No water transport related genes were found differentially expressed. Ion transport related genes were divided according to each ion type (Table 3). Among Ca 2+ transport related genes, a calcium ATPase-like gene was down-regulated in the leaves of sensitive plants as compared to those in tolerant plants, whereas the opposite was true for the roots. Furthermore, a calcium uniporter protein-like gene was found to be upregulated in the leaves of sensitive plants. Potassium channels SKOR and KAT3-like genes were downregulated in the leaves of sensitive plants, whereas a probable potassium transporter was upregulated in the leaves of sensitive plants. With regard to ionic channels, vacuolar cation/proton exchangers were upregulated in the leaves of sensitive plants. Other cation/proton antiporters did not show a clear response to salinity stress, as some isoforms had different expression patterns. On the other hand, chloride channel protein-like genes were upregulated in both leaves and roots of sensitive plants and downregulated in tolerant leaves. Comparing between samples showed that CLC-1-like gene showed higher expression in both tissues of sensitive plants. In the case of the CLC-2-like gene, the difference in expression was much marked in roots ( Figure 5). Phosphate, magnesium and boron transport-related genes were upregulated in tolerant roots. Phosphate, zinc, and ammonium transport-related genes were also upregulated in tolerant leaves. Regarding ATPases, plasma membrane ATPase-like genes were strongly upregulated in the roots of tolerant plants and downregulated in sensitive leaves. Furthermore, a V-type proton ATPase subunit G-like gene was strongly upregulated in leaves of tolerant plants as compared to the sensitive ones. Comparing the reads between samples, tolerant plants showed in both tissues higher expression of these genes compared to both control and sensitive groups in both tissues ( Figure 6). The comparison between expression in tolerant vs. sensitive plants revealed 2901 differentially expressed genes. These genes were further filtered using the BLAST results and annotations to look for putative genes related to salinity tolerance. Additionally, an orthogonal design analysis was conducted to be able to differentiate gene expression differences specific to roots or to leaves. The orthogonal analysis identified 1277 genes differentially expressed between tolerant and sensitive leaves (Supplementary Figure S4). In addition, 1395 genes were differentially expressed between tolerant vs. sensitive roots (Supplementary Figure S5).
Differentially expressed genes were divided according to several salt tolerance mechanisms: ion transport, photosynthesis and respiration systems, and ROS detoxification. No water transport related genes were found differentially expressed. Ion transport related genes were divided according to each ion type (Table 3). Among Ca 2+ transport related genes, a calcium ATPase-like gene was down-regulated in the leaves of sensitive plants as compared to those in tolerant plants, whereas the opposite was true for the roots. Furthermore, a calcium uniporter protein-like gene was found to be upregulated in the leaves of sensitive plants. Potassium channels SKOR and KAT3-like genes were downregulated in the leaves of sensitive plants, whereas a probable potassium transporter was upregulated in the leaves of sensitive plants. With regard to ionic channels, vacuolar cation/proton exchangers were upregulated in the leaves of sensitive plants. Other cation/proton antiporters did not show a clear response to salinity stress, as some isoforms had different expression patterns. On the other hand, chloride channel protein-like genes were upregulated in both leaves and roots of sensitive plants and downregulated in tolerant leaves. Comparing between samples showed that CLC-1-like gene showed higher expression in both tissues of sensitive plants. In the case of the CLC-2-like gene, the difference in expression was much marked in roots ( Figure 5). Phosphate, magnesium and boron transport-related genes were upregulated in tolerant roots. Phosphate, zinc, and ammonium transport-related genes were also upregulated in tolerant leaves. Regarding ATPases, plasma membrane ATPase-like genes were strongly upregulated in the roots of tolerant plants and downregulated in sensitive leaves. Furthermore, a V-type proton ATPase subunit G-like gene was strongly upregulated in leaves of tolerant plants as compared to the sensitive ones. Comparing the reads between samples, tolerant plants showed in both tissues higher expression of these genes compared to both control and sensitive groups in both tissues ( Figure 6).   Regarding photosynthesis and respiration differentially expressed genes, in leaves they were all upregulated in the tolerant plants as compared to the sensitive ones (Table 4), similarly to the differences in photosynthesis values previously measured ( Table 2).  Regarding photosynthesis and respiration differentially expressed genes, in leaves they were all upregulated in the tolerant plants as compared to the sensitive ones (Table 4), similarly to the differences in photosynthesis values previously measured ( Table 2). Regarding photosynthesis and respiration differentially expressed genes, in leaves they were all upregulated in the tolerant plants as compared to the sensitive ones (Table 4), similarly to the differences in photosynthesis values previously measured ( Table 2). Concerning the ROS detoxification related genes (Table 5), all except Thermospermine synthase ACAULIS5-like presented a higher expression in the leaves. Peroxiredoxin-like and peroxidase-like genes were upregulated in the tolerant plants except for one peroxidase-like that was upregulated in sensitive plants. Thermospermine synthase ACAULIS5-like was strongly upregulated in the roots of the tolerant plants as shown in the sample comparison ( Figure 7).   (Table 5), all except Thermospermine synthase ACAULIS5-like presented a higher expression in the leaves. Peroxiredoxin-like and peroxidase-like genes were upregulated in the tolerant plants except for one peroxidase-like that was upregulated in sensitive plants. Thermospermine synthase ACAULIS5-like was strongly upregulated in the roots of the tolerant plants as shown in the sample comparison ( Figure 7).

Ion Transport
Following salt treatment, no significant differences between salt-tolerant and sensitive plants were observed for Na + , Ca 2+ , or K + accumulation in leaves. However, a higher K + and Ca 2+ accumulation in the leaves of tolerant plants coincides with a higher leaf expression of several potassium channels and calcium transporting ATPases in the tolerant plants.
In our study, the total content of Na + and Cl − was analyzed and therefore a possible vacuolar compartmentalization could not be detected. However, we observed higher expression of all the analyzed H + -ATPases in tolerant plants. This suggest the involvement of these genes in the compartmentalization mechanism of cell protection from Na + toxicity. In spite of that the responsible proteins of the removal of sodium ions from the cytoplasm into the apoplast or vacuole are the specific plasma membrane Na + /H + antiporters, the activity of this protein depends on the electrochemical proton gradient between membranes. The H + /ATPases are the only proteins that are able to generate the adequate electrochemical proton gradient that allows these antiporters to perform the Na + extrusion [64]. Furthermore, the upregulation of this protein under salt conditions has been previously reported [65][66][67][68] and its critical role on the Na + extrusion has been confirmed using transgenic plants [69,70].
Another ion directly involved in damage caused by toxicity is Cl − , especially in the case of persimmon as previously shown by other works [44,71]. In a previous work, we detected differential chloride accumulation in a D. lotus population, however we could not link any previously measured parameters to these phenotypes [72]. In our study, we have observed a direct relationship between expression of chloride channel proteins and chlorine accumulation in leaf. Recently, chloride channels were proposed as key transporters of Cl − into the higher parts of the plant [43]. In persimmon, our results suggest that these channels might be playing a key role in avoiding high Cl − leaf accumulation in tolerant plants.

Photosynthesis, Respiration Systems, and ROS Signaling
A decreased metabolism in sensitive plants is supported both by the downregulation of photosynthetic and respiratory genes, and by a significant reduction in carbon fixation. Sensitive plants showed a reduction in photosynthesis ( Table 2) not limited to the stomatal aperture, as shown by the water use efficiency (ratio between ACO 2 and g s ). Sensitive and tolerant plants showed statistically similar WUE (Water Use Efficiency) after saline stress ( Figure 8). Therefore, the reduction of photosynthesis in sensitive plants compared to the tolerants is probably caused by photosystem damage due to the accumulation of toxic ions causing leaf necrosis and defoliation [2]. This hypothesis is also supported by the strong downregulation in photosystem genes in sensitive plants compared to the tolerant plants. ROS inactivation is also higher in both control and tolerant leaves an effect that probably leads to higher leaf ROS content in the sensitive plants. Although ROS can cause toxic effect in the leaves, it is a part of a signaling mechanism that may not be key to explaining the salinity tolerance, as an efficient ROS detoxification is essential for plant survival under natural conditions [2]. ROS are known to trigger programmed cell death [73] and tissue necrosis and defoliation, which are believed to be a mechanism for protecting the most photosynthetically active cells [74].

Root Architecture Involvement
Under stress conditions, thermospermine synthesis and K + uptake are linked to the inhibition of lateral root formation [75,76]. Although root architecture was not characterized in this study, transcriptomic data revealed thermospermine synthesis and K + transport activation by salinity in roots of tolerant plants, which suggests that root architecture may also play a role in the mechanism of tolerance to salt stress. Furthermore, the observed differences in sensitive plants regarding water potential (Table 2) cannot be linked to the transpiration as shown by the ratio between ψH and gs, as no significative differences were found between groups (Figure 9). We have neither found significative differences between expression of aquaporin genes related to water transport, as previously reported in a previous work [72]; therefore, observed water potential differences might be linked to water absorption through roots. Although previous reports have pointed out the association between reduced lateral root growth and sensitivity to saline stress trough better maintenance of Na + /K + ratio by accumulating Na + in the secondary root stele [77], other studies have pointed that in some cases, such as in seedlings, accumulation of Na + in the stele can cause lateral root quiescence or even lateral root primordia damage [78] and therefore reduced root development. In this line, D. virginiana, which is considered a salt-tolerant persimmon rootstock, the root system is characterized by a well-developed taproot [79]. D. lotus, a salt-sensitive rootstock, is characterized by a fibrous root system as compared to D. virginiana that presents a well-developed taproot [80], this morphological differences are believed to play an important role in the salinity tolerance in Diospyros rootstocks [81]. Our expression results indicate a possible association between root architecture and tolerance to salinity warrants further studies.

Root Architecture Involvement
Under stress conditions, thermospermine synthesis and K + uptake are linked to the inhibition of lateral root formation [75,76]. Although root architecture was not characterized in this study, transcriptomic data revealed thermospermine synthesis and K + transport activation by salinity in roots of tolerant plants, which suggests that root architecture may also play a role in the mechanism of tolerance to salt stress. Furthermore, the observed differences in sensitive plants regarding water potential (Table 2) cannot be linked to the transpiration as shown by the ratio between ψH and g s , as no significative differences were found between groups (Figure 9). We have neither found significative differences between expression of aquaporin genes related to water transport, as previously reported in a previous work [72]; therefore, observed water potential differences might be linked to water absorption through roots. Although previous reports have pointed out the association between reduced lateral root growth and sensitivity to saline stress trough better maintenance of Na + /K + ratio by accumulating Na + in the secondary root stele [77], other studies have pointed that in some cases, such as in seedlings, accumulation of Na + in the stele can cause lateral root quiescence or even lateral root primordia damage [78] and therefore reduced root development. In this line, D. virginiana, which is considered a salt-tolerant persimmon rootstock, the root system is characterized by a well-developed taproot [79]. D. lotus, a salt-sensitive rootstock, is characterized by a fibrous root system as compared to D. virginiana that presents a well-developed taproot [80], this morphological differences are believed to play an important role in the salinity tolerance in Diospyros rootstocks [81]. Our expression results indicate a possible association between root architecture and tolerance to salinity warrants further studies.

Root Architecture Involvement
Under stress conditions, thermospermine synthesis and K + uptake are linked to the inhibition of lateral root formation [75,76]. Although root architecture was not characterized in this study, transcriptomic data revealed thermospermine synthesis and K + transport activation by salinity in roots of tolerant plants, which suggests that root architecture may also play a role in the mechanism of tolerance to salt stress. Furthermore, the observed differences in sensitive plants regarding water potential (Table 2) cannot be linked to the transpiration as shown by the ratio between ψH and gs, as no significative differences were found between groups (Figure 9). We have neither found significative differences between expression of aquaporin genes related to water transport, as previously reported in a previous work [72]; therefore, observed water potential differences might be linked to water absorption through roots. Although previous reports have pointed out the association between reduced lateral root growth and sensitivity to saline stress trough better maintenance of Na + /K + ratio by accumulating Na + in the secondary root stele [77], other studies have pointed that in some cases, such as in seedlings, accumulation of Na + in the stele can cause lateral root quiescence or even lateral root primordia damage [78] and therefore reduced root development. In this line, D. virginiana, which is considered a salt-tolerant persimmon rootstock, the root system is characterized by a well-developed taproot [79]. D. lotus, a salt-sensitive rootstock, is characterized by a fibrous root system as compared to D. virginiana that presents a well-developed taproot [80], this morphological differences are believed to play an important role in the salinity tolerance in Diospyros rootstocks [81]. Our expression results indicate a possible association between root architecture and tolerance to salinity warrants further studies.

Conclusions
In Diospyros lotus L., several mechanisms are present to overcome the saline stress. In this paper, we could elucidate the ones that present variability within these species and therefore could be exploited for breeding new salt-tolerant rootstocks. We have observed reduced water potential in sensitive plants and we have observed differential expression in the thermospermine synthase gene that might have a role in root architecture. The role of root architecture in salt tolerance among Diospyros species has been previously pointed out by other authors [67]. On the other hand, chloride channel expression could facilitate Cl − transport inside the plant, and therefore cause higher Cl − accumulation in the higher parts of the plant. Finally, in the leaves, H + -ATPases can contribute to salt stress tolerance by generating the electrochemical membrane gradient required for Na + extrusion outside the cytoplasm.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4395/10/11/1703/s1, Figure S1: Acummulated explained variance of each principal component and eigenvalues for the principal component analysis of morphological variables. Figure S2: Example of Control, Tolerant and Sensitive phenotype of D. lotus L. plants after 60 days of salt stress. Figure S3: Acummulated explained variance of each principal component and eigenvalues for the principal component analysis of physiological and morphological variables. Figure S4: Venn diagram of the orthogonal analysis for identification of leaf-specific differential expressed genes between tolerant and sensitive plants. Figure S5: Venn diagram of the orthogonal analysis for identification of root specific differential expressed genes between tolerant and sensitive plants.