Chemical soil analysis
The rhizospheric soil pH of halophytes in both zones was alkaline, with similar values between plants and zones (Table 1). The lowest pH and EC were observed in the rhizosphere of B. salicifolia in Zone 1. The average EC of the rhizospheric soil of all Zone 1 plants significantly differed from those of Zone 2, except for B. salicifolia. The highest EC was in D. spicata (29.3 dS m− 1) in Zone 1, while the lowest was in B. salicifolia (5.2 dS m− 1 and 4.4 dS m− 1 in Zone 1 and 2, respectively). D. spicata of Zone 1 also had the highest sodium adsorption ratio (SAR), while in Zone 2, it was S. torreyana. The soil SAR of all Zone 1 plants differed from that of Zone 2, except for K. scoparia.
Table 1
Chemical properties of the rhizosphere soil of halophytes
| Halophyte | Zone 1 | | Zone 2 |
| | pH | EC (dS m− 1) | SAR (mmol L− 1) | | pH | EC (dS m− 1) | SAR (mmol L− 1) |
| Distichlis spicata | 9.8 ± 0.28a | 29.3 ± 2.58a* | 60.4 ± 2.34a* | | 9.3 ± 0.31A | 13.2 ± 2.33A | 15.0 ± 1.58B |
| Kochia scoparia | 9.6 ± 0.24a | 18.5 ± 2.69b* | 9.5 ± 2.67e | | 9.4 ± 0.28A | 9.4 ± 2.58A | 10.8 ± 2.06C |
| Eragrostis obtusiflora | 9.2 ± 0.27a | 22.4 ± 2.23b* | 24.3 ± 2.14c* | | 9.6 ± 0.24A | 8.5 ± 2.63B | 7.6 ± 1.06CD |
| Cynodon dactylon | 9.4 ± 0.27a | 14.6 ± 2.06c* | 22.7 ± 2.08c* | | 9.2 ± 0.22A | 5.2 ± 1.85B | 15.3 ± 2.12B |
| Suaeda torreyana | 9.8 ± 0.24a | 14.7 ± 2.10c* | 51.2 ± 2.82b* | | 9.4 ± 0.26A | 7.6 ± 2.13B | 20.9 ± 2.35A |
| Baccharis salicifolia | 8.7 ± 0.26b | 5.2 ± 2.20d | 14.1 ± 2.56d* | | 9.1 ± 0.20A | 4.4 ± 2.02B | 7.5 ± 1.02D |
| Solanum fructu-tecto | 9.4 ± 0.23a | 14.3 ± 2.38c | 14.2 ± 2.53d | | | | |
| Means and standard deviation are shown, n = 8. Identical lowercase letters indicate no statistical difference in pH, EC, and SAR between the plants in Zone 1 (Tukey, α = 0.05). Identical capital letters present no statistical difference in pH, EC, and SAR between plants in Zone 2 (Tukey, α = 0.05). (*) Shows a significant difference in pH, EC, and SAR in the same plant when comparing Zone 1 and Zone 2 (Tukey, α = 0.05). EC = Electrical conductivity, and SAR = Sodium adsorption ratio. |
In Zone 1, the soluble P soil concentration ranged from 24 to 14 mg kg− 1, while in Zone 2 from 10 to 11 mg kg− 1 of soil (data not shown). No differences in P soil concentration were found between zones of the same halophyte. The soil Na+ concentration range in Zone 1 was 124–524 mM, whereas, in Zone 2, the range was 94–447 mM (Fig. 1a). Soil Na+ concentration was similar between zones of the same halophyte species (Fig. 1a). The highest soil Na+ concentration was observed in S. torreyana. The average soil K concentration was 34.21 mM and 16.32 mM for Zones 1 and 2, respectively. In Zone 1, the highest K concentration was observed in the rhizosphere of D. spicata. In Zone 2, all plants had similar soil K concentration (Fig. 1b). The average rhizospheric Ca concentration was 3.19 mM in Zone 1 and 3.04 mM in Zone 2. In both zones, the concentration was similar between plants. When comparing the Ca concentration between zones in the same plant, only C. dactylon showed a difference (Fig. 1c). In Zone 1, the range of Mg concentration in the soil was 0.60–1.02 mM, while in Zone 2, the range was 0.28–0.59 mM. The average concentration was similar between zones, 0.80 mM in Zone 1 and 0.59 mM in Zone 2 (Fig. 1d).
Foliar physiological characteristics of halophytes
Foliar RWC values were between 46% and 92% (Fig. 2a). The average RWC of Zone 1 plants was 80%, and that of Zone 2 was 66%. The highest RWC from Zone 1 was observed in leaves of S. torreyana (92%), whereas B. salicifolia had the highest RWC (85%) in Zone 2 (Fig. 2a). D. spicata and K. scoparia had higher RWC in Zone 1 than in Zone 2. The other plants showed no statistical difference in RWC between the two zones.
The LSI of halophytes in this study was between 0.04% and 1.87%. Similar to what was observed in RWC, S. torreyana presented the highest LSI, 1.87% and 1.29% in Zones 1 and 2, respectively, while C. dactylon had the lowest LSI. No differences were observed when comparing LSIs of the same plant between zones (Fig. 2b).
In Zone 1, S. fructu-tecto had the highest leaf proline concentration, although it was similar in E. obtusiflora, D. spicata, and C. dactylon. For Zone 2, the proline concentration in leaves of D. spicata (1,312 µg g− 1 DW) and C. dactylon (1,045 µg g− 1 DW) was different from that in the other plants (Fig. 2c). The proline concentration was similar when comparing the same plant between zones with different EC.
The highest chlorophyll a concentration was observed in D. spicata (476 µg cm− 2) from Zone 1. Regardless of the zone, the other plants had low chlorophyll a concentrations; B. salicifolia (47 µg cm− 2) and K. scoparia (48 µg cm− 2) from Zone 2 had the lowest concentration (Fig. S2a). Only D. spicata and B. salicifolia showed differences in chlorophyll a between zones.
The chlorophyll b concentrations varied from 6 to 247 µg cm− 2, values lower than the chlorophyll a concentration. Like chlorophyll a concentration, D. spicata had the highest chlorophyll b concentration (247 µg cm− 2) in Zone 1 (Fig. S2b). Chlorophyll b concentration only differed between zones in K. scoparia and B. salicifolia.
The concentration range of total carotenoids fluctuated between 3,909 µg cm− 2 and 37,640 µg cm− 2. D. spicata and C. dactylon presented the highest concentration of total carotenoids (Fig. S2c). No differences in carotenoid concentration were observed between zones.
Foliar concentration and nutrient ratios in halophytes
The foliar Na concentration in halophytes ranged from 4 to 39 g kg− 1 (Fig. 3a). S. torreyana had the highest foliar concentration in both zones, 39 g kg− 1 in Zone 1 and 29 g kg− 1 in Zone 2. In contrast, the lowest Na concentration was observed in Zone 2, in C. dactylon and B. salicifolia, with 3 g kg− 1 and 3.75 g kg− 1, respectively. The foliar K concentration ranged from 3 to 66 g kg− 1. In Zone 2, all plants had a similarly low K concentration, approximately 11 g kg− 1. S. torreyana had the highest foliar concentration of K when established in Zone 1 (Fig. 3b). The foliar Ca concentration in the halophytes from Zone 1 ranged between 2 and10 g kg− 1. B. salicifolia and Solanum fructo-tecto had the highest Ca concentration in Zone 1 (10 g kg− 1). In Zone 2, plants had lower Ca concentrations (2–5 g kg− 1), and no difference was observed between plants (Fig. 3c). Only K. scoparia had different K concentrations when comparing Zone 1 and Zone 2. In Zone 1, E. obtusiflora, C. dactylon, and Solanum fructu-tecto had the highest Mg concentrations, with 6.53, 6.39, and 7.19 g kg− 1, respectively. In Zone 2, the range of foliar Mg concentration was between 1.53–7.19 g kg− 1, and C. dactylon had the highest Mg concentration (Fig. 3d).
The foliar P concentration in the halophytes was from 6 to 12 g kg− 1 (Fig. S3). The highest concentration was observed in B. salicifolia (12.2 g kg− 1) and S. fructo-tecto (12.4 g kg− 1) from Zone 1, while the lowest was in C. dactylon (6.8 g kg− 1) also from Zone 1. All plants from Zone 2 had similar foliar P concentrations (average 10.3 g kg− 1). Comparing zones, only K. scoparia and C. dactylon had different foliar P concentrations.
In Zone 1, only S. torreyana, D. spicata, and K. scoparia had a ratio Na/K higher than 1. In Zone 2, all halophytes had this ratio higher than 1, except C. dactylon (Fig. 4a). S. torreyana, E. obtusiflora, and B. salicifolia had a higher ratio in Zone 2 than in Zone 1. The ratio Na/Ca was in the range of 0.5 to 9.6. Only S. torreyana in both zones and K. scoparia from Zone 2 (Fig. 4b) had a high ratio Na/Ca (> 9). A difference in the ratio between zones was observed in K. scoparia. S. torreyana and K. scoparia established in both zones had a ratio higher than 5, and different from the ratio observed in the rest of the halophytes. S. torreyana had the highest Na/K and Na/Mg ratios in both zones (Fig. 4c and 4d).
Colonization frequency by AMF and DSE
All plant roots simultaneously had colonization by AMF and DSE structures. AMF formed in roots, coiled hyphae, vesicles, and arbuscules (Fig. S4). Total mycorrhizal colonization (14% − 48%) was higher in Zone 1 than in Zone 2 for all plants. D. spicata had the highest percentage of AMF colonization in both zones. S. torreyana had the lowest percentage of AMF in both zones (Fig. S5a). Colonization by arbuscules was between 6% and 16% (Fig. S5b), and by vesicles was from 2–9% (Fig. S5c). Roots of D. spicata and E. obtusiflora had the highest arbuscular colonization in Zone 1. In Zone 1, the highest colonization by vesicles was in the roots of D. spicata (8.4%), K. scoparia (8.9%), and B. salicifolia (9.8%). In Zone 2, all halophytes had similar vesicle colonization percentages.
DSE structures were observed in the roots of all halophytes as dark-colored septate hyphae and microsclerotia with intracellular growth in the root cortex (Fig. 5Sd). Average colonization by DSE in halophytes from Zone 1 was between 16–48%, and in Zone 2, 14–38%. S. torreyana had the highest percentage of DSE colonization in both zones (Fig. S5d). In D. spicata, colonization by DSE was lower than by AMF in both zones. The roots of S. torreyana, K. scoparia, and B. salicifolia had a higher percent colonization by DSE when established in Zone 1 than in Zone 2 (Fig. S5d).
Concentration of total glomalin in soil (T-GRSP) and roots (TR-GRP)
The range of T-GRSP concentration in Zone 1 was 0.08–0.12 mg g− 1, whereas Zone 2 had 0.02–0.15 mg g− 1 (Fig. 5a). In Zone 1, the highest T-GRSP concentration was observed in the rhizosphere of D. spicata (0.12 mg g− 1) and S. fructo-tecto (0.13 mg g− 1). In contrast, in Zone 2, the rhizosphere of K. scoparia had the highest concentration (0.15 mg g− 1). A higher concentration of T-GRSP at major ECs had C. dactylon and E. obtusiflora.
The average TR-GRP concentration was higher than T-GRSP. In Zone 1, the average TR-GRP was 0.25 mg g− 1, and in Zone 2, it was 0.18 mg g− 1 (Fig. 5b). The C. dactylon roots had the highest TR-GRP concentration in both zones. E. obtusiflora, D. spicata, and C. dactylon had higher TR-GRP in Zone 1 than in Zone 2.
Principal component and correlation analysis
The first two components of the PCA explained 58.6% of the total accumulated variance (Fig. 6a). PC1 accumulated 39.3% of the variance and the most significant variables in the upper right quadrant were colonization by AMF, Na/K, Na/Ca, Na/Mg ratios, leaf proline concentration, chlorophyll a and b, carotenoid concentration, EC, soil Ca and K concentration, and leaf Mg concentration. The most significant variables in the lower right quadrant were colonization by DSE, the LSI and RWC, leaf K and Na concentration, and rhizospheric soil Na concentration (Fig. 6a). PC2 explained 19.3% of the variance, and the most influential variables were leaf Mg, P, and Ca concentration, rhizospheric soil P, Mg, and Ca concentration, T-GRSP, and TR-GRP (Fig. 6a).
The PCA also showed that D. spicata, C. dactylon, and S. torreyana plants from the two zones clustered in PC1, whereas K. scoparia and E. obtusiflora clustered in PC2 (Fig. 6b). B. salicifolia from the two zones were part of PC, but grouped in different quadrats. Several significant correlations were observed, these are presented in Table S2.