Adaptation to Water and Salt Stresses of Solanum pimpinellifolium and Solanum lycopersicum var. cerasiforme

: Solanum pimpinellifolium and Solanum lycopersicum var. cerasiforme represent a valuable tool for tomato breeding, particularly for tolerance to abiotic stresses. Water stress and salinity are major constraints to tomato’s cultivation, and for which limited genetic variability has been reported within the cultivated species. We evaluated four accessions of S. pimpinellifolium and four of S. l. var. cerasiforme for their adaptation to water deﬁcit and salinity. The CO 2 assimilation rate, stomatal conductance, substomatal CO 2 concentration, transpiration rate, and leaf chlorophyll concentration were evaluated, as well as morphological and agronomic traits. The accessions showed a remarkable inter- and intra-species response variability to both stresses. Two S. pimpinellifolium accessions and one S. l. var. cerasiforme showed unaltered physiological parameters, thus indicating a good adaptation to water deﬁcit. Two S. l. var. cerasiforme accessions showed an interesting performance under salt stress, one of which showing also good adaptation to water stress. In general, both stresses showed a negative impact on leaf size and fruit fresh weight, especially in the big-sized fruits. However, ﬂowering, fruit setting and earliness remained unaltered or even improved when compared to control conditions. Stressed plants yielded fruits with higher ◦ Brix. Response to stresses seemed to be linked to origin environmental conditions, notwithstanding, variability was observed among accessions of the same region.


Introduction
Water deficit and salinity are two of the most important abiotic stresses that limit both productivity and quality of crops [1]. Plant responses to these stresses involve adaptive changes and, frequently, deleterious effects [2]. Water deficit prevents plant development by decreasing the plant's relative water content and water potential in their tissues, and consequently the closure of the stomatal complex. This leads to osmotic stress, limited nutrient uptake, reduced photosynthetic activity, oxidative stress, and growth inhibition [3]. Likewise, net photosynthesis can be lowered by mesophyll conductance restriction or by impairing CO 2 fixation reactions [4]. The effect of salinity on plants is more complex and involves two phases: osmotic stress, similar to the one induced by water deficit, and toxic effect produced by the ionic accumulation on leaves that gives rise to imbalances between nutrients and induces metabolic damage [5,6]. Plants have developed different adaptive strategies to face these flowering and fruits traits in an experiment comparing three treatments: water deficit, salinity, and a control. The relation between the responses and the origin of their climatic conditions and accession domestication level is discussed herein.

Plant Material
Plant material consisted of four accessions of S. pimpinellifolium and four accessions of S. l. var. cerasiforme ( Table 1). The range of distribution of the accessions of S. pimpinellifolium encompasses very humid environments as the ones of Manabí province (Ecuador) and the Amazonian region in the north of Peru, and dry areas of Northern Peru (Piura) and the coastal desert of Peru. In these areas, the range of temperatures and pluviometry is very high (Figure 1). The four accessions of S. l. var. cerasiforme come from very diverse areas, from the Sinaloa desert in Mexico, where S. l. cerasiforme grows as a wild, to the highly humid area of San Martín (Peru), where a richness of morphological variability exists. Plant material was provided by the COMAV's genebank (Universitat Politècnica de València, Valencia, Spain), except LA2251 which was provided by the Tomato Genetic Resources Centre (University of California, Davis, USA) and the accession PI 487625 from the US Department of Agriculture (USDA, Washington D.C., USA) ( Table 1).
After 15 days in pots (at six true-leaves stage), plants were divided into three groups for the control, water and saline stress treatments. Plants were irrigated by drip irrigation system and water doses in control plants was calculated weekly based on estimations of crop evapotranspiration (ETc) [41], and fractionated in two daily applications (morning and afternoon). The nutrient saline solution was allowed to drain freely from pots and the control drainage was controlled from 10% to 20% depending on solar radiation. Water stress treatment began by reducing the volume of irrigation water to 40% of the control. Salinity treatment began by adding NaCl (70 mM) to the irrigation solution to reach an EC of 8 dS m -1 and pH 6.1. Plants were pruned at one stem. Both stresses were conducted simultaneously.
Physiological and phenotypic measurements were taken at 75 and 120 days after the water deficit and salinity treatments (DAT) began. The layout was completely randomized with four replications per treatment and two plants per replication. A total of eight plants per accession were used in each treatment.

Experimental Design and Greenhouse Conditions
Tomato seeds were germinated in moistened perlite at 28 • C under greenhouse conditions. Seedlings were transferred to 15 [40]. The electrical conductivity (EC) of the nutrient solution was 1.4 dS m −1 and pH 6.1. After 15 days in pots (at six true-leaves stage), plants were divided into three groups for the control, water and saline stress treatments. Plants were irrigated by drip irrigation system and water doses in control plants was calculated weekly based on estimations of crop evapotranspiration (ETc) [41], and fractionated in two daily applications (morning and afternoon). The nutrient saline solution was allowed to drain freely from pots and the control drainage was controlled from 10% to 20% depending on solar radiation. Water stress treatment began by reducing the volume of irrigation water to 40% of the control. Salinity treatment began by adding NaCl (70 mM) to the irrigation solution to reach an EC of 8 dS m -1 and pH 6.1. Plants were pruned at one stem. Both stresses were conducted simultaneously.
Physiological and phenotypic measurements were taken at 75 and 120 days after the water deficit and salinity treatments (DAT) began. The layout was completely randomized with four replications per treatment and two plants per replication. A total of eight plants per accession were used in each treatment.

Gas Exchange Measurements
The CO 2 assimilation rate (A N , µmol CO 2 m −2 s −1 ), stomatal conductance to water vapor (g s , mol H 2 O m −2 s −1 ), substomatal CO 2 concentration (Ci, µmol CO 2 mol −1 air) and transpiration rate (E, mmol H 2 O m −2 s −1 ) were measured with a portable LI-COR 6400 infrared gas analyzer (Li-Cor Inc., Lincoln, NE, USA). In addition, parameters A N /Ci and A N /E were calculated as instantaneous carboxylation and water use efficiency, respectively, at the end of the experiment. Measurements were taken under saturating light conditions (1000 µmol quanta m −2 s −1 ), with reference CO 2 (400 µmol CO 2 mol −1 ) at 24 • C (24 • C ± 2) and 75% relative humidity (75% ± 10). To guarantee their accuracy, readings were kept after all gas exchange parameters were stabilized and their corresponding variation coefficients were lower than 1%. The terminal lobe of fully expanded and non-detached leaves (3rd-4th leaf position from the apex) were used for the measurements taken from 09:00 h to 11:00 h (UT + 01:00 h). Eight different plants were used (n = 8) per accession and treatment (control, water, and salt stresses).

Leaf Chlorophyll Concentration Measurement by Soil Plant Analysis Development (SPAD) Chlorophyll Meter
SPAD was estimated using a porTable SPAD-502 meter (Konica Minolta, Tokyo, Japan). Measurements were done in the same leaves used for gas exchange, taking 3 measurements in each leaf and calculating the mean value. Eight different plants were used (n = 8) per accession and per treatment (control, water, and salt stress) at the end of the experiment.

Agronomic and Phenotypic Traits
Plants were phenotyped for vegetative, flower and fruit traits in order to determine the effect of salinity and water deficit treatments on plant development. Vegetative traits recorded were length (cm) and width (cm) of the whole leaf and of the terminal lobe at 75 and 120 days after treatment started. The following traits were recorded in each plant for the first four trusses to measure stress effects on flowering, fruit set and earliness: number of flowers per truss, number of fruits per truss, fruit set (measured as the number of fruits divided by the number of flowers for the four first trusses), flowering and maturity earliness (measured as the number of days between the transplant and the opening of the first flower, and between the transplant and the first ripe fruit, respectively). Fruit fresh weight (g) was measured for each plant by calculating the average of ten ripe fruits individually weighted and randomly collected from the first four trusses at 120 days after the treatment began. Total soluble solids content ( • Brix) was measured using the juice of the ten weighted fruits using a manual refractometer (HI96801, Hanna Instruments S.R.L, Woonsocket, Rhode Island, USA).

Statistical Analyses
The results of all parameters were analyzed by two-way ANOVA (accessions × treatments; StatPoint Technologies, Warrenton, VA, USA) to estimate the main effects and their interactions considering all accessions (Supplemental material, Tables S1-S7). Given that plants showed a great phenotypic and origin variability, it has been also done one-way ANOVA to analyze the behavior of each accession individually after water and salinity stresses were applied and in comparison to their respective controls. There were not significant differences among plants from the same accession and treatment.
The analysis was performed at 75 and 120 after treatment began. The mean comparisons were made using Fisher's least significance difference (LSD) test at p < 0.05.

Results
S. l. var. cerasiforme accession Ceras4 showed a physiological disorder after the transplant to the greenhouse. This disorder consisted on the formation of intumescences on the abaxial part of the leaves that progressed distorting the leaves and giving a generalized necrosis that killed the plants. In some cases, the disorder affected also stems and petioles. This disorder is well known in some tomato varieties when cultivated under air-controlled conditions like greenhouses, and it has been also observed in some wild relatives like Solanum habrochaites [42]. It is clearly dependent on the genotype and the environmental conditions. Hence, this accession was excluded from the experiment.

Photosynthetic Rate and Stomatal Conductance
The two-way ANOVA revealed highly significant differences for the main effects (accession and treatment), as well as for their interaction, for both parameters (A N and g s ) at 75 and 120 DAT. Salt treatment effect showed a significantly higher negative impact on the genotypes' performances compared to water stress and control treatments regarding both parameters (Table S1). At the genotype's level, the photosynthetic rate (A N ) was not significantly affected by water deficit at 75 DAT for any accession (Figure 2A) when compared to control conditions. At the end of the experiment (120 DAT), A N decreased in Pim4 and Ceras3 water-stressed plants (59% and 50%, respectively) when compared with control plants, while it increased by 64% in Pim2 accession. The effect of salt stress on A N at 75 DAT was higher compared with water stress, reducing A N in three Pim accessions (Pim2, Pim3, and Pim4) (Figure 2A). At 120 DAT ( Figure 2B), salt stress provoked a significant reduction in A N ranging from 47% to 84% in all of the studied accessions, and particularly in Pim1 salt-treated plants (84%). Only A N in Pim2 accession under salinity did not show significant differences at 120 DAT when compared with control plants ( Figure 2B).  Table S1. AN: Net CO2 assimilation rate; gs: stomatal conductance; SE: Standard error.

Instantaneous Carboxylation Efficiency
Two-say ANOVA revealed highly significant differences for the factors "accession" and "treatment", as well as for their interaction for the instantaneous carboxylation efficiency parameter, estimated as the AN/Ci ratio,. AN/Ci ratio were significantly lower under saline treatment (Table S2).
Under water stress conditions, AN/Ci ratio was reduced in Pim1, Pim4 and Ceras3 accessions (ranging from 37% to 43%), compared to their control ( Figure 3A). Pim2 and Pim3 water-stressed  Table S1. A N : Net CO 2 assimilation rate; g s : stomatal conductance; SE: Standard error.
Stomatal conductance was not affected by water stress in Pim1, Pim4, and Ceras2 accessions at 75 DAT ( Figure 2C), while it was significantly increased in Pim2 (80% with respect to control). At 120 DAT only Pim2 and Ceras1 remained unaffected, while the other genotypes were significantly affected ( Figure 2D). All plants grown under salt stress significantly reduced g s when compared with controls (ranging from 53% to 86%) at 75 DAT. Finally, at 120 DAT, g s was again diminished in all accessions under salinity conditions ( Figure 2D).

Instantaneous Carboxylation Efficiency
Two-say ANOVA revealed highly significant differences for the factors "accession" and "treatment", as well as for their interaction for the instantaneous carboxylation efficiency parameter, estimated as the A N /Ci ratio. A N /Ci ratio were significantly lower under saline treatment (Table S2).
Under water stress conditions, A N /Ci ratio was reduced in Pim1, Pim4 and Ceras3 accessions (ranging from 37% to 43%), compared to their control ( Figure 3A). Pim2 and Pim3 water-stressed plants showed higher A N /Ci ratio than controls (increase of 51% and 18%, respectively) ( Figure 3A).

Leaf Chlorophyll Concentration Measurement by SPAD
Two-way ANOVA showed highly significant effects regarding the factors "accession" and "treatment", whereas their interaction (accession × treatment) was not statistically significant (Table  S3).
SPAD values, measured in terminal lobe of fully expanded and non-detached leaves, of salttreated plants showed the lowest values, while water-stressed plants tended to show higher values than the other two treatments, although statistically not different from the values observed under control conditions (Table S3 and Figure 4).
Two accessions (Pim3 and Ceras2) did not show significant differences either under salt stress conditions or under water stress when compared with respective controls.  Table S2. A N : Net CO 2 assimilation rate; Ci: internal CO 2 ; E: transpiration; SE: Standard error.
Under salt-stress conditions, A N /C i was reduced in all accessions when compared with their respective control or even water-stressed plants ( Figure 3A). This reduction was notably important in Pim1 (87%). Only Pim2 accession did not show significant differences when compared with control ( Figure 3A).

Instantaneous Water Use Efficiency
Regarding the instantaneous water use efficiency (A N /E), estimated as the ratio between CO 2 assimilation (A N ) and transpiration (E), the factors "accession" and "treatment" were statistically significant, as well as their interaction, as revealed by the conducted ANOVA. Salt-treated plants showed a significantly higher values than water-stressed plants and control (Table S2).
At the genotypes level, A N /E in Ceras3 showed a significant increase (41% increase) under water stress when compared with its control plants ( Figure 3B).
Pim3, Ceras2, and Ceras3 under salt-stressed plants registered higher A N /E values than control plants (65%, 74%, and 67% increase, respectively), while this ratio was not significantly affected in Pim2, Pim4, and Ceras1 salt-treated plants. On the contrary, Pim1 salt-stressed plants registered the highest decrease (55%) when compared with its control ( Figure 3B).

Leaf Chlorophyll Concentration Measurement by SPAD
Two-way ANOVA showed highly significant effects regarding the factors "accession" and "treatment", whereas their interaction (accession × treatment) was not statistically significant (Table S3).
SPAD values, measured in terminal lobe of fully expanded and non-detached leaves, of salt-treated plants showed the lowest values, while water-stressed plants tended to show higher values than the other two treatments, although statistically not different from the values observed under control conditions (Table S3 and Figure 4). internal CO2; E: transpiration; SE: Standard error.

Leaf Chlorophyll Concentration Measurement by SPAD
Two-way ANOVA showed highly significant effects regarding the factors "accession" and "treatment", whereas their interaction (accession × treatment) was not statistically significant (Table  S3).
SPAD values, measured in terminal lobe of fully expanded and non-detached leaves, of salttreated plants showed the lowest values, while water-stressed plants tended to show higher values than the other two treatments, although statistically not different from the values observed under control conditions (Table S3 and Figure 4).
Two accessions (Pim3 and Ceras2) did not show significant differences either under salt stress conditions or under water stress when compared with respective controls.  Table S3. SE: Standard error.
Two accessions (Pim3 and Ceras2) did not show significant differences either under salt stress conditions or under water stress when compared with respective controls.

Leaf Morphological Parameters
The two-way ANOVA showed that both accession and treatment effects were highly significant regarding leaf and terminal lobe's dimensions (length and width).Significant interactions were observed (accession × treatment) for the next parameters: leaf length both at 75 and 120 DAT, leaf lobe length at 75 DAT, leaf lobe width at 75 DAT. By contrast, there were not significant interactions for leaf width both at 75 and 120 DAT, leaf lobe length at 120 DAT, leaf lobe width at 120 DAT (Table S4).
In general terms, water stress led to a significant decrease of both leaf and terminal lobe dimensions in all genotypes at 75 DAT, except in Pim3, which was not significantly affected for both leaf length and width, and Ceras3, which showed no differences regarding leaf and terminal lobe width ( Figure 5A). Water stress treatment effect was considerably decreased at 120 DAT ( Figure 5B). Thus, only Ceras3 was negatively affected for both leaf and terminal lobe length and width, whereas Pim3 showed shorter terminal lobes, Ceras1 showed significantly shorter leaves, and Ceras2 thinner terminal lobes ( Figure 5B). Agronomy 2020, 10, x FOR PEER REVIEW 11 of 21  Table S4. SE: Standard error.

Flowers and Fruits per Truss Parameters
Both accession and treatment effect showed a highly significant impact on the parameters used herein. Their interaction (accession × treatment) was as well statistically significant (Table S5).
At the accessions level, S. pimpinellifolium accessions showed a complex response to the different treatments. Thus, significant differences between control and water-stressed plants were only Salt stress showed an overall significant negative effect in most accessions; however, accessions behavior was more variable than under water stress conditions ( Figure 5A). At 75 DAT, S. pimpinellifolium accession Pim1 and S. l. var. cerasiforme accessions Ceras1 and Ceras2 showed a significant decrease of both leaf and terminal lobe average dimensions, compared to control plants. Pim2 showed a significant decrease of its leaf's length, Pim3 significantly increased its terminal lobe average width, Pim4 leaf length and terminal lobe width were significantly affected, and finally accession Ceras3 decreased its average leaf length while increasing its terminal lobe width ( Figure 5A). On the other side, at 120 DAT the response was more stable across genotypes. Thus, significant differences between salt-stressed and control plants were only observed for Pim1 and Ceras3, which decreased their leaf length ( Figure 5B).

Flowers and Fruits per Truss Parameters
Both accession and treatment effect showed a highly significant impact on the parameters used herein. Their interaction (accession × treatment) was as well statistically significant (Table S5).
At the accessions level, S. pimpinellifolium accessions showed a complex response to the different treatments. Thus, significant differences between control and water-stressed plants were only observed for Pim3, which decreased the number of flowers per truss, and Pim2, which increased both the number of fruits per truss and fruit set ( Figure 6A-C). S. l. var. cerasiforme accessions did not show significant differences between water stress and control treatments regarding both the number of flowers per truss and the number of fruits per truss ( Figure 6A,B). Accession Ceras1 significantly increased its fruit set under water stress treatment, while the remaining accessions showed no differences ( Figure 6C).
Agronomy 2020, 10, x FOR PEER REVIEW 12 of 21 observed for Pim3, which decreased the number of flowers per truss, and Pim2, which increased both the number of fruits per truss and fruit set ( Figures 6A, 6B and 6C). S. l. var. cerasiforme accessions did not show significant differences between water stress and control treatments regarding both the number of flowers per truss and the number of fruits per truss (Figures 6A and 6B). Accession Ceras1 significantly increased its fruit set under water stress treatment, while the remaining accessions showed no differences ( Figure 6C). Regarding salt-treated S. pimpinellifolium, all genotypes showed fewer flowers per truss compared with control plants ( Figure 7A). In addition, the number of fruits per truss increased in the case of Pim2, decreased in Pim4 and remained unalterable in the rest of accessions ( Figure 6B). Finally, fruit set was significantly increased for Pim2 ( Figure 6C). Regarding S. l. var. cerasiforme accessions, there were no significant differences between salt and control treatments in any parameters ( Figures 6A, 6B and 6C).  Table S5. SE: Standard error.

Flowering and Maturity Earliness
Regarding flowering and maturity earliness, two-way ANOVA showed highly significant effects for both accession and treatment factors. In this case, the interaction was not statistically  Table S5. SE: Standard error.
Regarding salt-treated S. pimpinellifolium, all genotypes showed fewer flowers per truss compared with control plants ( Figure 7A). In addition, the number of fruits per truss increased in the case of Pim2, decreased in Pim4 and remained unalterable in the rest of accessions ( Figure 6B). Finally, fruit set was significantly increased for Pim2 ( Figure 6C). Regarding S. l. var. cerasiforme accessions, there were no significant differences between salt and control treatments in any parameters ( Figure 6A-C).
Agronomy 2020, 10, x FOR PEER REVIEW 13 of 21 significant for the trait days to flowering, whereas for days to fructification the interaction between both factors was highly significant (Table S6). Most accessions started to flower around 27 days after the beginning of the experiment ( Figure  7A). Regarding treatments, there was not a significant effect of water and salinity stress on this trait. Only Ceras1 plants grown under water stress presented a statistically significant delay of 3 days on the days to flowering when compared to control plants ( Figure 8A). Similar behavior was observed in the days to fructification ( Figure 7B). Only Pim3 plants grown under water stress presented a statistically significant delay of 5 days on the days to fructification when compared to their corresponding control plants. Under salt stress only Ceras1 plants presented a statistically significant delay of nine days to fructification when compared to their corresponding control plants and water stress ( Figure 7B).  Table S6. SE: Standard error.

Fruit Parameters
Regarding °Brix and fruit fresh weight, the two-way ANOVA showed statistically significant effects for accession and treatment. Despite that, significant interaction (accession × treatment) was only observed for fruit fresh weight (Table S7).
Total soluble solids content was positively affected by both stress treatments in most accessions, compared to control conditions (between 21% and 28% average increase, for water and salt stresses, respectively). In fact, only Pim4 was not significantly affected by any of the treatments ( Figure 8A). This accession showed the highest average ° Brix under the three treatments. Salt treatment showed the highest °Brix values, followed by water stress and then control. In addition, accessions Ceras2 and Ceras3 showed no differences between both stress treatments, accession Pim2 was not significantly affected by the water stress treatment ( Figure 8A).
Under stress conditions, plants tended to decrease their fruit size, especially when grown under salt stress (averaging 21% and 31% under water and salt stresses, respectively). Regarding waterstressed plants, statistically significant differences were observed in S. pimpinellifolium accession Pim4 and S. l. var. cerasiforme accessions Ceras2 and Ceras3, which showed fruit fresh weight decrease (around 15%, 17% and 30%, respectively) when compared with control plants. The marked difference between the fruit weight of Ceras1 with respect to Ceras2 and Ceras3 is due to the wild condition of Ceras1, while Ceras2 and Ceras3 beard fruits remarkably bigger than the typical small rounded fruits  Table S6. SE: Standard error.

Flowering and Maturity Earliness
Regarding flowering and maturity earliness, two-way ANOVA showed highly significant effects for both accession and treatment factors. In this case, the interaction was not statistically significant for the trait days to flowering, whereas for days to fructification the interaction between both factors was highly significant (Table S6).
Most accessions started to flower around 27 days after the beginning of the experiment ( Figure 7A). Regarding treatments, there was not a significant effect of water and salinity stress on this trait. Only Ceras1 plants grown under water stress presented a statistically significant delay of 3 days on the days to flowering when compared to control plants ( Figure 8A). Similar behavior was observed in the days to fructification ( Figure 7B). Only Pim3 plants grown under water stress presented a statistically significant delay of 5 days on the days to fructification when compared to their corresponding control plants. Under salt stress only Ceras1 plants presented a statistically significant delay of nine days to fructification when compared to their corresponding control plants and water stress ( Figure 7B).

Discussion
Water stress and salinity are two of the major abiotic constraints to tomato cultivation. Despite tomato being considered moderately sensitive to salinity compared to other Solanaceae's species, yield loss may be high and depends on both the salt concentration and the duration of the stress [43][44][45]. In addition, water stress can have important consequences for tomato production, as it might  Table S7. SE: Standard error.

Fruit Parameters
Regarding • Brix and fruit fresh weight, the two-way ANOVA showed statistically significant effects for accession and treatment. Despite that, significant interaction (accession × treatment) was only observed for fruit fresh weight (Table S7).
Total soluble solids content was positively affected by both stress treatments in most accessions, compared to control conditions (between 21% and 28% average increase, for water and salt stresses, respectively). In fact, only Pim4 was not significantly affected by any of the treatments ( Figure 8A). This accession showed the highest average • Brix under the three treatments. Salt treatment showed the highest • Brix values, followed by water stress and then control. In addition, accessions Ceras2 and Ceras3 showed no differences between both stress treatments, accession Pim2 was not significantly affected by the water stress treatment ( Figure 8A).
Under stress conditions, plants tended to decrease their fruit size, especially when grown under salt stress (averaging 21% and 31% under water and salt stresses, respectively). Regarding water-stressed plants, statistically significant differences were observed in S. pimpinellifolium accession Pim4 and S. l. var. cerasiforme accessions Ceras2 and Ceras3, which showed fruit fresh weight decrease (around 15%, 17% and 30%, respectively) when compared with control plants. The marked difference between the fruit weight of Ceras1 with respect to Ceras2 and Ceras3 is due to the wild condition of Ceras1, while Ceras2 and Ceras3 beard fruits remarkably bigger than the typical small rounded fruits of this variety. Finally, salt-stressed plants significantly decreased fruit fresh weight in Pim1, Pim2, Pim4, Ceras2 and Ceras3 (around 22%, 32%, 24%, 29% and 36%, respectively) ( Figure 8B).

Discussion
Water stress and salinity are two of the major abiotic constraints to tomato cultivation. Despite tomato being considered moderately sensitive to salinity compared to other Solanaceae's species, yield loss may be high and depends on both the salt concentration and the duration of the stress [43][44][45]. In addition, water stress can have important consequences for tomato production, as it might result in yield reduction of up to 50% in the case of an equivalent reduction in irrigation [46]. In this work, we explore the adaptation to water deficit and salinity of seven parents of a MAGIC population, four belonging to S. pimpinellifolium and three to S. l. var. cerasiforme.
Overall, there was an evident high variability in the response of the tested genotypes under the experimental conditions, demonstrating the existence of genetic variability susceptible of being exploited by breeders. The conducted ANOVA revealed the existence of significant differences between the tested accessions for all recorded traits. However, when looking at the results of the LSD test it can be concluded that there is not a clear effect of the "species" factor on the results. In our opinion, this result shouldn't be taken as a general rule, but as the consequence of the specific group of accessions tested in this experiment. These accessions were selected to maximize their genetic diversity, morphological traits, geographic origin, and degree of domestication, as they are the founders of a MAGIC population. As an example, Ceras1 grows as a wild and resembles phenotypically to S. pimpinellifolium. In the same way, the four accessions of S. pimpinellifolium were collected covering a wide range of distribution of this species, showing remarkable different morphology in some plant, inflorescence and fruit traits, as well as adaptation to different environmental conditions. Hence, the maximization of the diversity between the accessions belonging to each species probably masked the differences in the response between species. Regarding the applied treatments, the effect of the water stress treatment was less harmful than the one produced by salinity. These results have also been described in other crops, such as pepper [47], as a consequence of the specific toxicity produced by salt ion accumulation, mainly Na + and Cl − , when plants are subjected to salinity stress.
In this experiment, S. pimpinellifolium accessions did not show a decrease in A N at 75 DAT under water stress. At 120 DAT only Pim4 decreased A N , remaining this parameter unaltered in Pim1 and Pim3 and increasing in Pim2. However, the stomatal conductance diminished at 120 DAT in three accessions except in Pim2. In spite of the stomatal closure detected in Pim1 and Pim3, it did not interfere with the photosynthetic rate (A N ) indicating that stomatal closure was more sensitive to water stress than A N [48]. This could be explained by the fact that only critical levels of g s , described as low as 0.1 mol H 2 O/m·s are able to affect photosynthesis [49]. At the end of the productive period, the instantaneous carboxylation efficiency increased in Pim2 and Pim3, diminishing in Pim1 and Pim4. The decrease in A N /Ci indicated that water stress affected the photosynthesis by metabolic limitations in Pim1 and Pim4 [50,51]. Water use efficiency (A N /E) remained unaltered in the four genotypes. Considering these data as a whole, we consider the accessions Pim2 and Pim3 as tolerant to the water stress treatment applied in this experiment. Regarding the recorded agronomic and morphological traits, most of them remained unaffected, decreasing only the number of flowers per truss in Pim3 and increasing the number of fruits per truss and the fruit set in Pim2. The vegetative vigor, measured as the leaves size, decreased at 75 DAT, but this effect disappeared as the treatment progressed. The fruit weight was only reduced in Pim4. The • Brix increased in almost all Pim accessions as a consequence of the water deficit treatment. This response is in agreement with the one found by other authors [52,53]. According to Albert et al. [54], this can be due to a reduction in fruit water content and not to increased synthesis of sugars, although Ripoll et al. [55] found higher fructose and glucose synthesis in tomato fruits submitted to water deficit stress at different stages of fruit development, indicating that both dilution effect and higher sugar synthesis are responsible of fruit quality enhancement in tomato under water deficit conditions.
The four tested accessions in this experiment come from areas with very different climatic conditions, ranging from the desertic climate of the coastal northern Peru (Pim2) and the coastal south Peru (Pim4), to the highly humid and hot areas or Manabí (Pim1) in Ecuador and Amazonas province in the north of Peru (Pim3). Nakazato et al. [56] conducted an interesting study using S. pimpinellifolium accessions from different western slopes of the Andes and the adjacent coastal regions of Peru, to demonstrate that environmental factors drive the phenotypic differentiation and adaptation to environmental stresses, and showing evidence of trait-environment associations. They found a negative correlation between tolerance to water stress, measured as the number of days to wilting, and the annual precipitation, indicating that individual populations that occur in arid environments are more likely drought tolerant [56]. Interestingly, Pim2, one of the most tolerant accessions to water stress in our experiment due to photosynthesis preservation, comes from desert areas with scarce precipitation. This behavior fits with the results found by Nakazato et al. [56]. Nevertheless, in a study conducted by Albert et al. [54], in which some S. pimpinellifolium accessions from different origins were tested, no correlation was found between the climate and the adaptation to water deficit, showing different adaptation to water stress of accessions coming from the same province. Herein, Pim3, coming from Amazonas province (Peru) with climatology clearly different to Piura (Pim2) was also found tolerant to the water stress conditions applied in this study. Zuriaga et al. [15] conducted a comprehensive experiment to study the molecular variability of S. pimpinellifolium, including 247 accessions covering all its range of distribution. They found a clear genetic differentiation among the studied accessions, related to their geographical origin [15]. The accessions used herein belong to groups genetically and environmentally differentiated by Zuriaga et al. [15], demonstrating that the differences found have a genetic basis and can be exploited in breeding.
In contrast to the response to water stress observed in S. pimpinellifolium, salinity treatment severely affected the photosynthetic rate, as well as the stomatal conductance, particularly at the end of experiment. The ion imbalance, ion toxicity and osmotic stress produced by the salinity treatment interfered with both parameters, as described by other authors [7,57]. Herein, salt treatment lead to a decrease in chlorophyll levels, as demonstrated by the significant differences shown by SPAD measurements in Pim1 and Pim2, not observed under water stress conditions. As a consequence of the decrease in photosynthesis efficiency, the vegetative development also decreased, as shown by the smaller leaf and leaf terminal lobe sizes, compared to control. Despite that, fruit set (except for Pim4) and flowering earliness were not affected. Finally, a negative effect was observed on fruit fresh weight, which has been described as one of the most important effects of sensitivity to salinity, being more notable in bigger fruits [30,35,58]. The decrease in fruits as small as the ones of S. pimpinellifolium, shows the negative effects of salinity in the tested accessions.
Regarding S. l. var. cerasiforme, the photosynthetic rate did not vary in water stressed compared to control plants at 75 DAT, although Ceras1 and Ceras3 reduced their g s , compared with control. At the end of the experiment (120 DAT), Ceras1 and Ceras2 showed no alteration of CO 2 fixation, whereas Ceras3 showed a decrease of g s values. This allowed for dynamic modulation of photosynthesis under water stress, with conservative values of A N /Ci indicating no photosynthesis metabolic limitations in both Ceras1 and Ceras2 [47], however in Ceras3 decrease in A N /Ci and in g s indicating that A N reduced by stomatal closure and by inhibition of mesophyll conductance and/or photochemical efficiency at 120DAT, similar results were obtained in [47]. The water use efficiency remained unaltered in Ceras1 and Ceras2 and increased in Ceras3, caused by a considerable stomata closure associated to a decrease in the transpiration rate [51]. Regarding the phenotypic traits, only leaf size, at 75 DAT, and fruit fresh weight (not in Ceras1) decreased as consequence of the treatment. Decrease in plant vigor, as a result of the water stress treatment, was also reported for a group of 55 accessions of S. l. var. cerasiforme and S. pimpinellifolium tested by Albert et al. [54] whereas the accessions' yield remained unaffected. The authors suggested that tomato plants buffer the negative effects of water deficit by limiting their vegetative growth and by reallocating the photo-assimilates to the fruits. According to our results, the most tolerant S. l. var. cerasiforme accessions were Ceras1 and Ceras2 for the absence of alteration of CO 2 fixation, unaltered water-use efficiency, and maintenance of most of the vegetative and reproductive traits unaffected. Interestingly, Ceras1 was also tested for tolerance to water deficit by Albert et al. [54] with similar conclusions to ours, where yield was not decreased and plant vigor was slightly reduced. Finally, Ceras3, the accession with bigger fruits, showed the most marked reduction in fruit fresh weight. The greater weight loss in larger fruited accessions has already been reported by several authors [35,54].
Under salinity stress, at 75 DAT, no limitation of CO 2 occurred in S. l. var. cerasiforme accessions. At 120 DAT stomatal conductance and photosynthesis were negatively affected by salinity in all accessions. Herein, A N /Ci decreased in all S. l. var. cerasiforme accessions, suggesting stomatal constraints. Previous studies have demonstrated a positive relationship between photosynthetic capacity and growth in the plants grown under salinity [47,[59][60][61]. As a consequence of the considerable reduction of carbon assimilation rate, vegetative growth, leaf and leaf terminal lobe sizes decreased in all tested genotypes.
S. l. var. cerasiforme, unlike S. pimpinellifolium grows normally in the proximity of cultivated fields or in backyards in close contact with humans [37]. In these conditions, plants can compete with fast growing species or genotypes such as those in cultivation or in backyards, but do not need to develop special adaptive traits to abiotic stresses [56]. Notwithstanding, S. l. var. cerasiforme can also be found in wild conditions. This is the case of Mexico where this species is widely distributed and can be found as a wild, weedy and even as partially cultivated varieties in tropical and subtropical areas with semiarid and humid regimes [62]. When growing in semiarid and hot climates this species is found in association with larger plants that provide them with shade. Accession Ceras1 comes from the Sinaloa desert, in Mexico, in the Pacific slopes (ca. 300-1100 m) of Sierra Madre Occidental [63], while Ceras2 and Ceras3 come from humid and hot areas. Hence, the three S. l. var. cerasiforme accessions come from very different conditions. In our experiment, Ceras1 has been the most tolerant accession under both treatments while Ceras2 showed an interesting behavior under water stress deficit. Given the fact that Ceras1 grows as a wild variety in semiarid conditions, it may indicate that this accession could have a natural adaptation to similar stresses conditions to those applied here. Morphologically, there are marked differences between the three accessions. Ceras1 is similar to S. pimpinellifolium with thin stem, small leaves, and fruits of similar weight, shape and size (between 1 and 1.5 cm of diameter). In contrast, the other S. l. var. cerasiforme accessions (Ceras2 and Ceras3) are closer to the cultivated tomato, with stronger stem and bigger leaves and fruits of more than 3 cm in diameter. The genetic and morphological variability of S. l. var. cerasiforme was already shown by Blanca et al. [37]. In the three accessions used in this work, strategically selected by their genetics, environmental conditions of origin and morphology, we have found different adaptation to abiotic stresses, showing that this species constitutes a valuable source of variability for tomato's breeding.

Conclusions
Water stress and salinity are major constraints to tomato cultivation and genetic variability has been found limited within the cultivated species. Tomato wild relatives have been reported as important genetic sources of abiotic tolerance traits. The studied accessions showed a remarkable amount of variability regarding water and salt stresses adaptation, corroborating the richness of these materials in genetic variability. Accessions Pim2, Pim3, and Ceras1 showed promising results regarding water stress, while Ceras1 and Ceras2 showed an interesting performance under salt stress. The aptitude of some genotypes to improve under a specific treatment seemed to be linked to the environmental conditions of their region of origin, although variability can be found for adaptation to abiotic stresses among accessions of the same area. The study of the MAGIC population will enable the elucidation of the genetics underlying the tolerance found in some of these accessions and their exploitation in future breeding programs.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4395/10/8/1169/s1, Table S1: Two-way ANOVA analysis of net CO 2 assimilation rate (A N ) and leaf stomatal conductance (g s ) in 4 Solanum pimpinellifolium accessions and 3 Solanum lycopersicum var. cerasiforme accessions under control, water stress and salt stress. Table S2: Two-way ANOVA analysis of instantaneous carboxylation efficiency (A N /C i ) and water use efficiency (A N /E) in 4 Solanum pimpinellifolium accessions and 3 Solanum lycopersicum var. cerasiforme accessions under control, water stress and salt stress. Table S3: Two-way ANOVA analysis of SPAD values in 4 Solanum pimpinellifolium accessions and 3 Solanum lycopersicum var. cerasiforme accessions under control, water stress and salt stress. Table S4: Two-way ANOVA analysis of leaf and terminal lobe length and width in 4 Solanum pimpinellifolium accessions and 3 Solanum lycopersicum var. cerasiforme accessions under control, water stress and salt stress. Table S5: Two-way ANOVA analysis of the number of flowers per truss, the number of fruits per truss and the fruit set in 4 Solanum pimpinellifolium accessions and 3 Solanum lycopersicum var. cerasiforme accessions under control, water stress and salt stress. Table S6: Two-way ANOVA analysis of the days to flowering and the days to fructification in 4 Solanum pimpinellifolium accessions and 3 Solanum lycopersicum var. cerasiforme accessions under control, water stress and salt stress. Table S7: Two-way ANOVA analysis of the total soluble solids content ( • Brix) and the fruit fresh weight in 4 Solanum pimpinellifolium accessions and 3 Solanum lycopersicum var. cerasiforme accessions under control, water stress and salt stress. Funding: This research was funded by the Spanish Ministerio de Economía y Competitividad and the Fondo Europeo de Desarrollo Regional/European Regional Development Fund, grant number AGL2015-71011-R. Authors also thank the G2P-SOL (Linking genetic resources, genomes, and phenotypes of Solanaceous crops) and BRESOV (Breeding for resilient, efficient, and sustainable organic vegetable production) projects for support. G2P-SOL and BRESOV projects have received funding from the European Union's Horizon 2020 research and innovation program under grant agreements 677379 (G2PSOL), and 774244 (BRESOV). David Alonso is grateful to Universitat Politècnica de València for a predoctoral (PAID-01-16) contract under the Programa de Ayudas de Investigación y Desarrollo.

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