Effect of Pisolithus tinctorious on Physiological and Hormonal Traits in Cistus Plants to Water Deficit: Relationships among Water Status, Photosynthetic Activity and Plant Quality

Cistus species can form ectomycorrhizae and arbuscular mycorrhizal fungus that can bring benefits when plants are under water stress conditions. However, the application of some ectomycorrhizae on the water uptake under drought through physiological traits and hormonal regulation is less known. The experiment was performed during three months in a growth chamber with Cistus albidus plants in which the combined effect of the ectomycorrhiza Pisolithus tinctorious inoculation and two irrigation treatments (control and water-stressed plants) were applied. Irrigation absence caused significant decrease in aerial growth and tended to decrease soil water potential at the root surface, leading to a decrease in leaf water potential. Under these conditions, the abscisic acid and salicylic acid content increased while the precursor of ethylene decreased. Although the mycorrhization percentages were not high, the inoculation of P. tinctorious improved the water status and slightly cushioned the rise in leaf temperature of water-stressed plants. The ectomycorrhiza decreased the scopoletin values in leaves of plants subjected to deficit irrigation, indicating that inoculated plants had been able to synthesize defense mechanisms. Therefore, Pisolithus tinctorious alleviated some of the harmful effects of water scarcity in Cistus plants, being its use a sustainable option in gardening or restoration projects.


Introduction
Cistus albidus L. is a typical shrub of the Mediterranean climate. They are widely cultivated both for reforestation and for its use in Mediterranean gardening, especially in Xeriscape, where water management for saving this resource is one of the best options. These plants began to be widely used in gardening and landscaping due to their rusticity and attractive blooms.
Cistus albidus responds to water deficit by developing avoidance mechanisms for the regulation of transpiration based on stomatal closure, reduction of the leaf area and epinastia [1]. Furthermore, tolerance to water stress can be explained by other functional and structural adaptations of plants, such as osmotic adjustment, changes in the elasticity of the cell wall, and mineral and hormonal balance [2].
It is widely known that mycorrhizae confer physiological characteristics to plants that help them improve their water and nutritional status, especially in adverse conditions such as water stress. However, the response of these plants will depend on both the species of the plant and the type of fungus with which it is inoculated, as well as the extension, duration, and type of stress [3,4].
Cistus species can form both ectomycorrhizae (ECM) and arbuscular mycorrhizal fungus (AMF); they are dualistic plants [5][6][7]. The dual mycorrhization state is also present in other species such as Populus, Salix, Eucalyptus [5,6]. In an earlier experiment we studied whether an AMF, Glomus iranicum var. Tenuihypharum, despite being a genus (Cistus) that generally forms ectomycorrhizae, could bring benefits, especially when plants are under water stress conditions. In that experiment, mycorrhizae helped improve the water status of the plant when irrigation was removed, as demonstrated by the less negative leaf water potential values; it also minimized the decrease in chlorophyll levels, showing plants with higher ornamental value and increased in quantum yield (Y(II)) and photochemical quenching (qP) were related to the greater photosynthetic efficiency recorded [8].
Several studies have shown beneficial effects of ectomycorrhizal symbioses on the performance of tree species under drought stress [9]. For example, seedlings of Picea colonized by ECM fungi exhibited increased stomatal conductance and photosynthesis, shoot water potential and growth compared to non-mycorrhizal plants [10]. These ameliorative effects have, at least partly, been ascribed to the increasing in the absorbing surface area and the exploration of larger soil volumes by the extramatrical mycelia, as well as to the role of aquaporins which function is related to post-translational regulation and to the coordination of the phytohormones [11,12]. Transport of water and nutrient uptake through hyphae to host plants improves plant water and nutrient supply improving the performance of plants during periods of drought stress [9]. Furthermore, drought-stressed ECM-colonized plants show increased hydraulic conductance compared with non-mycorrhizal plants [13,14]. Concretely, Pisolithus tinctorius is one of the most widespread ECM fungi and can establish associations with a broad variety of species and has been recorded in a wide range of habitats [15]. P. tinctorius has great mycorrhizal capacity and provides the plant with greater efficiency in the absorption of water and nutrients [16], so it has been used in reforestation program sand and in commercial ECM inoculum production [17].
In general, mycorrhizal fungi plays an important role in regulating phytohormones such as jasmonic acid (JA), salicylic acid (SA), and ethylene (ET), which are involved in the host defense response, and abscisic acid (ABA), gibberellin (GA), and citokinin (CK), which modulate the growth of the plants. However, while the role of hormones in AMF-inoculated plants is widely known [18,19], less is known in ECM-inoculated plants, although, fundamentally, studies of the effect of ABA and SA in ECM-inoculated plants under adverse conditions have been carry out. For example, Rincón et al. [20] evaluated whether two ectomycorrhizal fungi (Laccaria bicolorthe and Cenococcum geophilum) affected ABA production in larch during osmotic stress and Luo et al. [21] determined the phytohormone changes in plants inoculated with Paxillus involutus under salinity. However, in the case of P. tinctorius, the relationships between water status and hormonal production in plants subjected to water stress is less known.
For all the above, in this study we examine how the interaction of Cistus albidus with the ectomycorrhizal fungus Pisolithus tinctorious affects water uptake under drought through physiological traits and hormonal regulation. The objectives of the study were (i) to study physiological and hormonal response to water deficit; (ii) to evaluate the effect of P. tinctorious on physiological traits and hormonal regulation of Cistus plants; and (iii) to analyze the changes in the relationships among physiological, morphological, and ornamental parameters under water deficit and mycorrhizal inoculation. Comparative data relating these physiological processes to water stress may prove beneficial toward understanding the drought tolerance of plants.

Mycorrhization Percentage and Growth
At the end of the experiment, mycorrhizal inoculation produced a colonization of around 30% in root systems of well-irrigated Cistus plants. The level of colonization in roots of mycorrhizal plants decreased near to 20% with the suppression of the irrigation (Table 1). In both cases, the mycorrhization percentages were not high. Water stress caused significant decrease in aerial growth parameters such as height and leaf DW, being the latter the one that decreases the most (30% respect to C) ( Table 2). The inoculation of Pisolithus tinctorious increased almost all growth parameters for both irrigation conditions, as observed when we compare CM with C and WSM with WS. Indeed, inoculated plants of both irrigation treatments (CM and WSM) were taller and their leaf and root DW were higher than those corresponding to non-inoculated plants (C and WS) ( Table 2).

Leaf Mineral Ontent
Water stress had a significant effect on leaf mineral content decreasing almost the entire mineral measured, when comparing C and WS, having a highly significant effect in Ca, Na, Mg, and Cu. The effect of the mycorrhiza treatment was hardly significant. Under well-watered conditions, only the leaf content of Ca and S in CM plants increased their values both in around 30% respect C plants (Table 3). Under water stress conditions, ECM produced a slight tendency to decrease almost all the elements of WSM plants respect to WS (Table 3).

Water Relations
Soil water potential at the root surface (Ψ r ) tended to decrease in WS plants due to water stress, causing greater resistance to water absorption, which leads to a decrease in leaf water potential. No differences were observed due to the effect of mycorrhizae ( Figure 1A).
In well-watered plants, both inoculated and non-inoculated (CM and C, respectively), the values of leaf water potential (Ψ l ) remained at −1.0 MPa, showing a good water status of the plants ( Figure 1B). When stress was applied, non-inoculated plants (WS) decreased Ψ l significantly, reaching minimum values of −3.0 MPa. The use of Pisolithus tinctorius alleviated significantly the negative effect of water stress and the fall in values was not so severe (−2.1 MPa) ( Figure 1B). Figure 1C shows the contribution to dehydration to change in osmotic water potential (∆Ψ SS ). Under well-watered conditions, ∆Ψ SS was almost negligible for both inoculated and non-inoculated plants (C and CM, respectively). Water stress caused a decrease in ∆Ψ SS , reaching minimum values of around −0.27 MPa in non-inoculated plants (WS), while the ECM dampened this drop (−0.12 MPa).

Gas Exchange
Under water stress conditions, the stomatal conductance (g s ) fell drastically, almost reaching the stomatal closure in WS plants. Differences of 150 mmol m −2 s −1 between C and WS treatments were observed ( Figure 2A). Pisolithus tinctorious did not favor opening in either case, neither under conditions of good irrigation nor under water deficit (CM and WSM). Net photosynthesis (P n ) showed a behavior very similar to that of g s ( Figure 2B). Leaf temperature (T l ) of plants subjected to water stress (WS) was 1.5 • C higher than C plants ( Figure 2C).

Phytohormones
At the end of the experiment period, the abscisic acid (ABA) and salicylic acid (SA) content in WS leaves were approximately 5 and 3 times as high as those in control leaves, respectively ( Figure 3A,B). By contrast, the precursor of ethylene (ACC) content was decreased significantly by the irrigation suppression (45% of C) ( Figure 3C). Mycorrhizae only had a significant effect on scopoletin (SC) under water stress conditions (WSM), significantly decreasing their values by approximately 46% respect to non-inoculated plants (WS) ( Figure 3D).

Relationships between Physiological, Morphological and Ornamental Parameters
The relationship between g s and P n was well represented by an exponential function when all the values of the four treatments obtained during the experimental period were considered ( Figure 4A). Highly significant linear regressions between g s versus P n can be obtained when considering control treatments (C and CM) and water stress treatments (WS and WSM) individually, suggesting that g s explain differences in P n within each irrigation treatment. A strong correlation (R 2 = 0.86) between g s and leaf water potential (Ψ l ) was observed, fitting with an exponential sigmoid model. Ψ l of plants subjected to water stress (WS and WSM) declined sharply (variation between −1.0 and −3.2 MPa) in a narrow range of g s ( Figure 4B). Relationship between T l and g s fits with an exponential decay model. Leaf temperature was maintained between 23.3 and 24.3 • C when the stomatal conductance exceeded approximately 65 mmol m −2 s −1 ( Figure 4C). With increasing water stress, the dependence of Ψ l and T l on g s increased, as well as the dependence of the relative chlorophyll content (RCC) on P n . WS as compared to C promoted intrinsic water use efficiency (WUE, P n /g s ). However, P n was less dependent on shoot dry weight when the plants were subjected to irrigation suppression (Table 4). Under water stress conditions, Pisolithus tinctorious only had effect on the relation between RCC and P n , being WSM higher than WS (Table 4).

Discussion
The percentage of mycorrhization of Pisolithus tinctorius was not excessively high, despite being an ectomycorrhiza, which has special relevance in this genus with high mycorrhizal capacity [6]. In our conditions, the cistus seedlings were cultivated in pots, which prevented the extension of the fungal mycelium and of the ECM roots. Imposed water stress meant a reduction in the percentage of mycorrhization, like plants, fungi are also affected by water limitation [22]. Studies of dually colonized plant species such as members of the genus Cistus indicate that ECM may be more sensitive to drought than AMF [23,24]. In fact, in a previous study with AMF (Glomus iranicum var. tenuihypharum) in Cistus albidus plants under water stress [8], the percentage of mycorrhization was 65%, more than double that obtained with Pisolithus tinctorius.
Drought is one of the main adverse factors for seasonal plants, especially for plants grown in pots [1,25,26]. At the end of the experiment, water stress caused a decrease in leaf dry weight, which could be interpreted as a conservative strategy to reduce transpiration and maintain hydraulic conductivity under water depletion [27][28][29]. Despite the low percentage of mycorrhization, ECM improved both shoot and root growth of the plants under well-irrigated conditions. Furthermore, unlike what occurred in Cistus plants inoculated with AMF [8], Pisolithus tinctorious increased root growth under water stress. Usually, more root branching has been found in ECM species than AMF, what is related to higher ability to absorb more water and nutrient, improving plant water relations under low water conditions [30][31][32][33][34].
Leaf water potential and actual osmotic potential are good indicators of water stress [35,36]. In our study, plants subjected to water stress showed no osmotic adjustment, perhaps, the speed of the development of water stress and the low inorganic solutes accumulation, such as potassium and sodium, could not contribute to osmo-regulation [8,37]. The fact that ∆Ψ SS at the end of the experiment in control treatment were almost null indicated that dehydration was the major mechanism involved in ∆Ψ SS , while the more noteworthy ∆Ψ SS in WS treatment seems to indicate that the dehydration could be produced by the transpiration [38]. Soil water potential at the root surface (Ψ r ) reflects the accumulation of net solutes allowing reduction of Ψ l (around −3.0 MPa), in order to guarantee water transport to the leaves [37,[39][40][41]. Concerning ECM effect, it improved the plant water status under stress conditions [42]. This ameliorative effect has been ascribed to the extended external mycelia of the root systems of ECM fungi, which reach soil pores inaccessible to the roots in water [5,43,44].
As a consequence of the reduction of the substrate water potential at the root surface and the leaf water potential in plants subjected to water stress, the stomatal conductance decreased drastically, acting as a mechanism to prevent excessive loss of water [25,26,[45][46][47]. In our conditions, the relationship between g s and Ψ l was adjusted to an exponential curve, so that under water-stress plants exhibited progressively lower Ψ l , but almost no change in g s . Likewise, it has been well established that plants regulate rates of transpiration and photosynthesis in parallel, maintaining a balance between g s and P n [48]. A strong correlation between stomatal conductance with net photosynthesis observed in the current study appears to reflect the gas exchange limitation of photosynthesis [49][50][51]. At the end of the experiment, the imposed water stress induced an increase in intrinsic WUE, which is in agreement with Costa et al. [46], who reported that deficit irrigation strategies can be successfully applied to different crops, in order to improve water savings. Water deficit by closing stomata, causes increasing leaf temperature, considering the major determinant of leaf temperature is the rate of evaporation or transpiration from the leaf [52][53][54]. This behavior was observed at the end of our experiment.
ABA synthesis is one of the fastest responses of plants to abiotic stress causing stomatal closure. In addition, SA is involved in the regulation of drought responses [55], inducing the generation of reactive oxygen species (ROS) and, consequently, causing the stomatal closure [56]. In our experiment, endogenous ABA and SA levels in stressed plants increased up to of six and almost twofold than those of control plants, respectively, while endogenous 1-aminocyclopropane1-carboxylic acid (ACC) levels in leaves decreased. ACC is the precursor for the production of ethylene and under drought, a close relationship between ACC expression and ethylene synthesis has been demonstrated [57]. In these conditions, ACC levels decrease and, therefore, plant growth is reduced. However, some studies suggest that ABA and ethylene are antagonists and have described decreases in ethylene production as a consequence of an increased concentration of ABA in water stressed plants, regulating some drought responses in plants, such as root and leaf growth [58][59][60][61][62]. This is in accordance with our results, which show that in conjunction with a decrease in stomatal conductance and leaf water potential, ABA in leaves increased, while ethylene decreased, increasing root growth.
Hormonal profiles are also altered to alleviate the negative effects of water stress on mycorrhizated plants [32,[63][64][65]. Several studies confirmed that inoculation with mycorrhizal fungi results in decreasing the endogenous ABA and SA [20,66]. However, in our assay ECM treatments (CM and WSM) had not significant changes in any of both leaf hormones neither under well-watered nor water stress conditions. Since ABA is known to induce stomatal closure [67], its lower concentration in ECM plants could play a role in the stomatal conductance, allowing no increase in net CO 2 assimilation rate in ECM inoculated plants [17]. This fact leads no significant improvement in intrinsic water use efficiency (P n /g s ).
Nevertheless, ECM significantly decreased leaf scopoletin (SC). Scopoletin (6-methoxy-7-hydroxycoumarin) is a typical phytoalexin which is an important secondary metabolite synthesized in plants as a defense mechanism against various environmental stresses [68]. Its synthesis is activated once some kind of infection has occurred in plants, but it can also be triggered due to various types of abiotic stresses. Scopoletin has also been found in many other plant species (e.g., Solanaceae, such as tobacco or potato, and sunflower, among others), showing antifungal and antibacterial activity [68][69][70]. Its accumulation has been correlated with resistance to microbial attack and other stresses, as well as mechanical damage and dehydration [71]. Furthermore, it appears to be the product that most increases in concentration in infected plants compared to other coumarins and coumaric glycosides, such as scopolin, esculetin, and esculin [72][73][74]. In our study, SC decreased significantly in mychorrhizal plants indicating that ECM plants have been able to synthesize defense mechanisms against abiotic and biotic stresses against non-mycorrhizal plants, and therefore SC is lower in ECM plants.
In conclusion, the mechanism of C. albidus to avoid water stress was related to its ability to decrease aerial growth and to modify leaf gas exchange, increasing water use efficiency. Water stress positively stimulated the levels ABA and SA, who strongly enhanced drought tolerance. On the other hand, despite the low percentages of mycorrhization, Pisolithus tinctorious inoculation in addition to improving plant growth, and slightly gas exchange and quality plant parameters, it improved the water status of Cistus under water stress conditions, probably due to a decrease in SC values, suggesting that the selected ECM can alleviate some of the harmful effects of water scarcity. Therefore, the use of Cistus plants inoculated with Pisolithus tinctorious in gardening or restoration projects can be a sustainable economic and environmental option.

Plant Material and Experimental Conditions
Forty-eight plants of Cistus albidus L. from the nursery were used. On January 18, 2017, these plants were transplanted into 1.5 L pots with a commercial substrate based on Sphagnum peat, coconut fiber and perlite in an 8:7:1 ratio and were placed in a controlled growth chamber. The climatic conditions of the chamber were those necessary for optimal growth: temperature (23 • C/18 • C, day/night), photosynthetic photon flux density (350 µmol m −2 s −1 ), photoperiod (16h/8h, light/dark), and relative humidity (RH) (60%). Plants were watered at field capacity.

Treatments
After two weeks of acclimatization, half of the plants (24) were inoculated with the ectomycorrhizal fungus Pisolithus tinctorius. Two months later, once the plants were stabilized, the irrigation of half of inoculated and non-inoculated plants were removed for a month, while the remaining plants continued to be irrigated at field capacity. Thus, the following treatments were established: C, Control, well-watered plants; CM, well-watered and inoculated plants; WS, plants subjected to withholding; WSM, inoculated plants subjected to water stress. Each treatment included three replications. Each replication consisted of four plants. The experiment lasted three months.

Fungal Colonization
At the end of the experiment, roots samples with the surrounding rhizosphere soil were collected to treat to evaluate fungal development. Three root samples were used in each replication. The percentage of mycorrhizal root colonization was estimated as following: Once cleaned, the roots were immersed in KOH (100 • C for 7 min), followed by a bath in H 2 O 2 (100 • C, 5-6 min) and finally, trypan blue staining (4 min). The percentage of colonization was calculated using the methodology proposed by Kormanik and McGraw [75]. The colored roots were placed on specialized plates for counting and were observed under the magnifying glass by counting 100 fields neatly (positive and negative). The percentage of colonization was calculated using the following formula proposed by Sieverding [76]: %colonization = (number of colonized fields/total number of fields observed) × 100 (1)

Leaf Mineral Content
At the end of the experimental period, the inorganic mineral content of dry leaves was determined in three plants per replication by means of emission spectrophotometry. The nutrient concentrations were determined in a digestion extract with HNO3:HClO4 (2:1, v/v) by Inductively Coupled Plasma optical emission spectrometer (ICP-OES IRIS INTREPID II XDL, Thermo Fisher Scientific Inc., Loughborough, UK).

Biomass and Height
At the end of the experiment, three plants per replication were selected and all the substrate was gently washed from their roots. After washing the substrate from the roots, the plants were individually separated into leaves, stems and roots and the fresh weight of each organ was determined. After that, each sample was dried in an oven at 80 • C, until samples reached a constant weight, and the dry weights (DW) were obtained.
Every 15 days, the height of three plants per replication was measured. Plant height was measured from the base of the plant at the substrate surface to the most distal growth.

Water Relations
Leaf water potential (Ψ l ), leaf osmotic potential (Ψ s ) and osmotic water potential at full turgor (Ψ 100s ) were determined in three plants per replication throughout the experimental period. The leaf water potential was determined during light hours according to the technique described by Scholander et al. [77], using a pressure chamber (Model 3000; Soil Moisture Equipment Co., Santa Barbara, CA, USA), the leaves were pressurized at a rate of 0.03 MPa s -1 . Ψ s and Ψ 100s was measured using a Wescor 5520 vapor pressure osmometer (Wescor Inc., Logan, UT, USA), according to Gucci et al. (1991) [78]. For Ψ 100s the leaf samples were previously subjected to a rehydration treatment by dipping their petioles in distilled water for 24 h to achieve complete saturation.
Changes in Ψ s (∆Ψ s ) and in Ψ 100S (∆Ψ 100S ) were calculated according to Girma and Krieg [79] as the difference in Ψ s and Ψ 100S measured at the initial (i) and at the end (e) of water stress period.
The contribution of dehydration to changes of Ψ S (∆Ψ SS ) was calculated using the following equation: Soil water potential at the soil-root interface (Ψ r ) was computed according to Jones (1983) [80]: where, Ψ l WS and Ψ l C correspond to the mean value of the leaf water potential in the WS and C treatments, respectively, and the g s WS and g s C correspond to the mean value of the respective treatments. A value of Ψ r was calculated for each of the three drought pots using the value of Ψ l C and g s C . The Ψ r is assumed to be zero for control plants.

Gas Exchange and Thermography
Stomatal conductance (g s ) and net photosynthesis rate (P n ) were measured with the LI-COR 6400 portable meter (LI-COR Inc., Lincoln, NE, USA). The flow rate of circulating air within the system was approximately 200 mmol s −1 , with a leaf vapour pressure deficit to air of about 2 KPa. The CO 2 concentration was fixed at 380 ppm and the photosynthetically active radiation (PAR) at 600 µmol m −2 s −1 . The measurements were carried out in three plants per replication on the same days as the water relations.
At the same time and in the same plants that gas exchange was determined, leaf temperature (T l ) was measured with an infrared camera (FLIR-e50 System, Inc., Danderyd, Sweden) which consisted of a 240 × 180 pixels line scan imager operating in the 7.5e13 mm region, with a noise equivalent temperature difference of 0.05 • C at 30 • C and an accuracy of 2 • C or 2% of the reading. The background temperature, distance of the camera from the canopy, air temperature, emissivity, and relative humidity were used as input at the start of each series of measurements; so, the camera automatically corrects for atmospheric transmission based on these data. Background temperature was determined as the temperature of a crumpled sheet of aluminium foil in a similar position to the leaves of interest with the emissivity set at 1.0. Emissivity for leaf measurements was set at 0.96 [81] and the distance at which images were taken was 0.5 m. The images were processed with ThermaCam Researcher Professional 2.10 software (FLIR Quick Report, Danderyd, Sweden).

Relative Chlorophyll Content
The relative chlorophyll content (RCC) was determined periodically in three plants per replication with a Minolta SPAD-502 chlorophyll meter (Konica Minolta Sensing Inc., Osaka, Japan), a non-destructive method.

Hormonal Determination
The extraction and analysis of the plant hormones were developed in agreement with Albacete et al. [82] including slight modifications in the protocol, such as the measurement of scopoletin (phytoalexin). Firstly, 100 mg of fresh plant material (leaf) was homogenized in liquid nitrogen and placed in 0. The separation of plant hormones and phytoalexins was developed with slight modifications in accordance with Albacete et al. [82] by using UHPLC coupled to a 6460 triple quadrupole mass spectrometer (Agilent Technologies, Waldbronn, Germany), and a BEH C 18 column (2.1 × 50 mm, 1.7 µm) (Waters, Milford, MA, USA) with a guard column (2.1 × 5.0 mm, 1.7 µm). The column temperature was 40 • C. Water/acetic acid (99.99:0.01, v:v) (solvent A) and acetonitrile (solvent B) were used as mobile phases at the flow rate of 0.2 mL min −1 . The injection volume was 10 µL. The gradient program used was: 19% B at 0 min, 90% B at 2.5 min, 90% B at 4.5 min, 19.0% B at 6 min, and 19% B at 8.0 min for column equilibration. The electrospray interface (ESI) was set up in the negative and positive mode and the mass spectrometry analysis were run in the multiple reaction monitoring modes (MRM). The ionization and fragmentation conditions were as follows: gas temperature 325 • C, gas flow 8 L/min, nebulizer 45 psi, sheath gas temperature 375 • C, jet stream gas flow 11 L/min, capillary voltage 4000 V and 2750 V (positive and negative mode, respectively), and nozzle voltage 1000 V and 1500 V (positive and negative mode, respectively) according to the most abundant product-ions. The quantitative evaluation of plant hormones and scopoletin was carried out using authentic standards. The standards used are included in Table 5 and were diluted in water/methanol (20:80, v/v) for the quantification curves.

Statistical Analysis
In the experiment, 48 plants were randomly attributed to each treatment, with three replications for each treatment. The data were analysed by one-way ANOVA using Statistical Package for the Social Sciences (IBM SPSS Statistics 26 for Windows, CA, USA). Treatment means were separated with Duncan's Multiple Range Test (p ≤ 0.05). Prior to the statistical analysis, percentage of root colonization was subjected to an arcsine squareroot transformation to ensure the homogeneity of the variance. Relationships between parameters were fitted to different regressions using SigmaPlot v. 14.5 software (SPSS Inc., Chicago, IL, USA).

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.