Physiological Beneficial Effect of Rhizophagus intraradices Inoculation on Tomato Plant Yield under Water Deficit Conditions

Increasing drought, under current climate change scenarios, will reduce the sustainability of tomato cultivation in the Mediterranean region. The present study evaluates the effect of Rhizophagus intraradices inoculation on tomato plant physiology and yield in response to progressive water deficit conditions. Two commercial products (Prod1 and Prod2) containing only R. intraradices were tested at two different concentrations (1% and 5% of the substrate volume) using three methods of inoculation: (a) mixed to substrate, (b) dissolved in water, (c) spread on seedlings root blocks before transplant. The highest mycorrhization of root fragments (F%) was found with Prod2 at 1% w/w at 40 days after sowing (DAS); this product was therefore used in a second experiment to inoculate tomato plants and test their physiological response to progressive water deficit induced withholding irrigation. Phenology, plant height, stem diameter, chlorophyll content and fluorescence, whole canopy gas exchange, biomass production and partitioning and phosphorus content were investigated in inoculated and not inoculated tomato plants under well-watered and water stressed conditions. Vegetative period and plant height were shorter in inoculated than in control plants; moreover, inoculation with R. intraradices increased fruit production by enhancing chlorophyll content under water stress condition, PS2 efficiency, ETR, Fv/Fm, net photosynthetic rate and whole canopy WUE.


Introduction
According to FAO [1], tomato (Solanum lycopersicum) is the second most important vegetable crop, after potato, considering global production (about 164 million tons of fresh fruit harvested) and cultivated surfaces (4.7 million ha). For its high antioxidant content, it is recommended that an adult human eat from 25 to 32 kg of tomatoes per year [2]. In the Mediterranean region, where tomato cultivation requires high water and fertilization inputs, agro-techniques should be improved to increase the sustainability of tomato production while maintaining high yield levels to satisfy the needs of a growing population [3]. Under the current climate change scenario, the increasing temperature and the variation in rainfall intensity and distribution will increase drought and temperature-related stresses leading to reduction in crop productivity [4][5][6].
Research to increase resilience to drought and to mitigate the effects of climate change on crop yield [7] exploited conventional breeding techniques and genetic engineering to create stress-tolerant

Plant Material and Growth Conditions
The effects of methods of distribution and concentration of two commercial formulations of Rhizophagus intraradices were evaluated in a preliminary experiment (Exp1). Tomato (Solanum lycopersicum cv 'Zelig') seeds were surface sterilized in sodium hypochlorite 3.5% for 10 min, washed in sterile water, and then put on wet paper in Petri dish to germinate. Plastic plateau (S40C, http://www.plastisavio.it/contenitori) with 83 cm 3 of soil volume per plant and a final density of 245 plants m −2 were filled with Sphagnum peat (Gartenmarkt, https://www.tuttogiardino.it/it/giardino/terriccio-e-concimi/34-3121.html) having the following characteristics: bulk density 0.27 g cm −3 , organic matter content 40%, 5.5-6.5 pH, and 68% v/v of water content. 1200 g of sterilized peat (30 g per plant) were homogenously mixed to the inoculum of R. intraradices.
Three different modes of distribution of the two commercial products were tested: (1) The commercial products were homogeneously mixed to the substrate (2) The commercial products were dissolved in water and provided as irrigation on the germinating seedlings (3) The commercial products were spread on seedling root blocks before the transplant in the field.
No inoculum was added to the substrate for non-inoculated plants (M-), whilst two concentrations of inoculum (1% and 5% w/w, resulting in 12 g and 60 g of inoculum per 1200 g of topsoil) were used to induce root colonization. Twenty plants per treatment were transplanted in plastic plateau after germination and were grown in a greenhouse, with a temperature of 25 °C/21 °C day/night, 500 μmol m −2 s −1 PPFD and 16/8-h light/dark photoperiod, and relative humidity of 60% ± 10%. From germination to the first sampling time (40 days after sowing, DAS), all plants were watered twice per week with tap water and, once per week, with a modified Long-Ashton nutrient solution [39] containing 300 μM inorganic phosphate (Table S1) until plateau water capacity (~ 600 ml per plateau). A total of 40 plants per inoculum (Prod1 and Prod2) were used to assess the symbiosis development at 40, 50, 60 and 70 DAS (10 plants per sampling day). Considering the amount and the cost of each formulation and the average number of plants per square meter in an industrial tomato crop, the cost of mycorrhizal tomato plantlets production was estimated by multiplying the cost of the inoculum for the amount of inoculum used for the average number of plants per square meter (Table 1). For Experiment 2 (Exp2) tomato seeds were treated as described in Exp1 and put on wet paper in a Petri dish for germination. Seedlings were then transplanted in closed plastic pots (3.4 L capacity) containing Sphagnum peat (Gartenmarkt) and Prod2 at 1%, thereafter M+ plants. For non-inoculated plants (M-), no inoculum was added to the Sphagnum peat (Gartenmarkt). The plants were grown in greenhouse under the same growing conditions (temperature, light intensity and duration, relative humidity and nutrient solution) used in Exp1.
Six plants per treatment (M-and M+) were randomly assigned to a well-watered (WW) or to a water-stressed (WS) treatment. On 95 DAS, all plants were irrigated until field capacity and the top of the pots was covered with PVC bags. A little slit was cut at the bottom of the plastic bag to allow tomato plant to grow through it. The slit was sealed with adhesive packing tape around plant collars to minimize evaporation losses from the soil surface. The onset of water stress was monitored gravimetrically according to [40]. The Fraction of Transpirable Soil Water (FTSW) was calculated as follows: (1) where ASWC is the Available Soil Water Content, calculated subtracting the soil water content at Wilting Point (WP) from the actual Soil Water Content SWC, and TTSW is the Total Transpirable Soil Water, calculated as difference between the soil water content at Field Capacity (FC) and that at WP. FC and WP were 68 and 15% (w/w) respectively and were determined gravimetrically in a preliminary experiment according to [41]. The addition of 1% of Prod2 increased FC and WP of about 1%.
The WW M-and M+ plants were continuously irrigated to maintain FTSW above 0.8. For the remaining 6 plants, irrigation was withheld for 13 days until 0.28 FTSW and 0.35 FTSW were reached for M-and M+ plants, respectively. This level of stress was then kept for ten days, from 113DAS to 123 DAS (Figure 1).

Assessment of AM Symbiosis Development
Two hundred randomly chosen 1 cm-long root segments, 20 root segments from 10 plants representative of each experimental condition, were stained with 0.1% Cotton Blue in lactic acid at each sampling date in Exp1. Fungal colonization was then quantified according with the Trouvelot method [42] using the MYCOCALC software freely available at: www.dijon.inra.fr/mychintec/Mycocalc-prg/download.html. Intensity of the mycorrhizal colonization in the root system (M%) and arbuscular abundance in the root system (A%) were calculated in this experiment. In Exp2 the fungal colonization was determined at the end of the experiment (123 DAS) following the same protocol described for Exp1. In Exp2 the root system was split into two subsamples: one part was used to determine root symbiosis, whilst the remaining roots were weighted, dried at 65 °C for 4 days and weighted again to determine dry weight and P concentration. Plant height, stem diameter and phenology were measured once a week for the entire duration of Exp2 (0-123 DAS). Whole canopy leaf area was determined once a week using the software ImageJ (NIH, http://rsb.info.nih.gov/ij/, Bethesda, MD, USA) on scanned leaves, and the results were used to normalize gas exchange data. Chlorophyll content was also non-destructively determined with SPAD 502 Chlorophyll Meter (Minolta Co., Ltd., Osaka, Japan) once a week from 100 to 123 DAS.

Chlorophyll Fluorescence
Chlorophyll fluorescence parameters were measured by a Fluorescence Monitoring System (Hansatech Instruments, Norfolk, UK) on the last fully opened leaf according to the method described by Zhu et al. three times from 103 to 120 DAS. Measurements on dark adapted leaves were carried out before LED lamps switching on, therefore on leaves dark adapted for 8 hours. A saturating pulse of irradiation of 2000 μmol m −2 s −1 was applied for 3 second to measure the maximal fluorescence (Fm) in the dark-adapted state. The leaves were then placed under actinic light of 3000 μmol m −2 s −1 to determine the maximal fluorescence (Fmʹ) and the minimal fluorescence in the light-adapted state (F0ʹ), and the steady-state fluorescence (Fs). Photochemical quenching coefficient (qP), nonphotochemical quenching coefficient (NPQ), the maximal quantum yield of PSII in the dark-adapted state (Fv/Fm), effective photochemical efficiency of PSII (ФPSII) and electron transport rates (ETR) were calculated according to the following Equations [43]:

Gas Exchange Measurements
From 100 DAS to 120 DAS, tomato plants were enclosed in an open automated self-assembled multi-chamber gas exchange apparatus, as described in [44]. The system enables continuous recording of CO2 and H2O differentials for the entire duration of the experiment. The air flow entering in each sampling chamber is measured upstream each chamber in order to calculate net photosynthetic and transpiration rate and regulated according to plant leaf area in order to prevent CO2 and H2O stagnation inside the sampling chamber. In Exp2 the air flow rate in the chambers was set at 0.0015 m 3 s −1 . The polyethylene sampling chambers had a volume of 0.068 m 3 ± 0.008, so a complete volume air change occurred approximately every 50 second. Diurnal direct and diffuse light was always constant from 7 a.m. to 7 p.m. for the entire duration of the experiment. Net photosynthetic (Pn, μmol s −1 m −2 ) and transpiration rates (E, mmol s −1 m −2 ) were calculated from flow rates and CO2 and water vapor differentials using the formula provided in [45] and normalized for the whole canopy leaf area previously measured. Whole-canopy instantaneous WUE (wcWUEi) was calculated as Pn/E ratio (μmol CO2/ mmol H2O) according to [44]. The daily average of whole canopy transpiration rate and wcWUEi was plotted against daily FTSW values. A two-slope linear model was fitted [40] in order to calculate the FTSW thresholds (FTSWt) per each treatment. This index represents the FTSW value below which transpiration rate starts to decrease [46].

Plant Biomass Production, Partitioning, Growth and P Content
Plant biomass production and partitioning was determined on 123 DAS. The aboveground plant biomass was divided into leaves, stems, and fruits. Tomato leaves were scanned (Scanjet 3970) and the software ImageJ was used for leaf area determination. Leaves, stems and fruits were put in an oven at 65 °C for 4 days to determine dry matter and P content.
Leaf area and dry weight of tomato leaves were used to calculate Specific Leaf Area (SLA, cm 2 g −1 ) and Leaf Area Ratio (LAR, cm 2 g −1 ), according to the Equations: Half of the belowground plant biomass was dried at 65 °C for 4 days, weighted and then used for P content determination. Total P content (ppm) of leaves, stems, root, and fruit was determined through acid digestion (1:100) using nitric acid method and analysed by inductively coupled plasma emission ICP-OES (mod. 5100, Agilent Technologies, Santa Clara, CA, USA).

Statistical Analysis
Analysis of variance was performed using the Agricolae R package. Three-way ANOVA, with product type, concentrations and sampling date as main factors, was carried out in Exp1 and, when the F-test was significant, mean separation was performed by Tukey test at p < 0.05 and p <0.01. Degree of variation around means is given as standard error (SE). In Exp2 two-way ANOVA, with inoculation and drought stress as main factors, followed by post-hoc Tukey test was performed.

Assessment of AMF Symbiosis Development
Out of the three mode of AMF distribution (data not shown), only mixing commercial inoculum to peat resulted in significant level of root mycorrhization. At the highest product concentration (5%), frequency of mycorrhization of root fragments (F%) was highest with Prod2 in the first 2 sampling times but no significant difference of F% was found between Prod 1 and Prod2 at 60 and 70 DAS (days after sowing). The effect of product concentration was minor and rather inconsistent in time and at the last sampling time ( Figure 2A).
The intensity of root cortex colonization (M%) was significantly affected both by product and concentration ( Figure 2B). With Prod2 M% increased linearly from the first (40 DAS) to the last sampling time (70 DAS) without significant differences between product concentrations. With Prod1 M% started later, at DAS 50 at the highest concentration (5%) and at 60 DAS at 1%, but then progressed linearly at a rate similar to that measured with Prod2. At 70 DAS, M% was higher using Prod2 (at both concentrations) than using Prod1 at 1% (M% = 28) and Prod1 at 5% (M% = 47).
Inoculation costs were estimated considering the purchase of the AMF inoculum only, as applications cost were considered negligible. Inoculation of 1 ha of a processing tomatoes crop (40,000 plants ha −1 ) ranged from 104 € with Prod2 at 1% concentration to 920 € with Prod1 at 5% concentration (Table 1).

Figure 2.
Frequency of mycorrhization F% (A) and intensity of root cortex colonization M% (B) measured at 40, 50, 60 and 70 DAS (days after sowing) mixing the commercial products (circles for prod1 and diamonds for prod2) to the substrate (modality "a" of Exp1). Circles represent prod1: solid line for 5% and dotted line for 1% concentration respectively; diamonds represent prod2: solid line for 5% and dotted line for 1% concentration respectively. Asterisks indicate statistically significant differences between the two commercial products (p < 0.05).
We considered the amount of product used in order to obtain a frequency of mycorrhization (F%) above 30% at 40 DAS, that is an average lifetime of tomato seedlings when they are sold to the farmers. Considering the cost of Prod1 and Prod2 per seedling and an average number of plant m −2 , the use of mycorrhizal tomato seedlings with Prod2 at 1% will increase by 12% the actual cost per hectare of tomato seedlings (Table 1).
The results obtained in the preliminary experiment (Exp1), modality "a", were instrumental to determine the commercial product and the relative concentration to be used in Exp2. In Exp2, Prod2 at 1% concentration was used for a deeper characterization of plants responses to AMF colonization and drought stress. At the end of the experiment (123 DAS), frequency of mycorrhization of root fragments was F = 80%, intensity of root cortex colonization was M = 70%, average presence of arbuscules within the infected areas was a = 95% and arbuscules abundance in the root system was A = 70%.

Characterization of Drought Dynamics
The average time course of soil water deficit of M+ (mycorrhizal) and M-(non-mycorrhizal) tomato plants investigated under WW (well-watered) and WS (water stressed) conditions is expressed as Fraction of Transpirable Soil Water (FTSW) and was measured since the beginning of the experiment. In Figure 1

Chlorophyll Content and Fluorescence
Leaf chlorophyll content, expressed as SPAD values, was higher in WW M+ plants (46.9) than in WS M+ plants (41.1) and M-plants (WW and WS, 36.7 and 34.7, respectively) considering the average of all the SPAD values recorded during the experiment ( Figure 3A). During the measurements inside the gas exchange apparatus (100-120 DAS), SPAD values increased for M+ plants and slightly decreased for M-ones ( Figure 3A). All WW plants showed no significant change of SPAD values during the experiment ( Figure 3A). Maximum quantum yield (Fv/Fm), photosystem II efficiency (ΦPSII), electron transport rate (ETR), photochemical quenching (qP) and non-photochemical quenching (NPQ) were monitored once a week between 100 and 120 DAS ( Figure 3B-3F). Differences in Fv/Fm between the two water treatments increased progressively during the stress period, reaching statistical significance on 120 DAS (p < 0.05, Figure 3B). Differences in Fv/Fm between WS M+ and Mplants became significant on 114 DAS (p < 0.05) and kept increasing until 120 DAS (p < 0.01; Figure  3B). On 114 DAS, ΦPSII was on average higher for M+ plants than for M-(p < 0.01, Figure 3C). ΦPSII was reduced by drought stress to a greater extent in M-plants than in M+ ones (p < 0.01; Figure 3C). Similarly, ETR was higher in M+ than in M-plants both in WW and WS conditions ( Figure 3D). Drought stress affected ETR significantly in M-plants but not in M+ plants ( Figure 3D). ETR of WS M-plants was significantly lower than ETR of WW M-plants on 120 DAS (p < 0.01, Figure  3D). Differences in qP between water treatments and between M-and M+ plants became significant only on 120 DAS (p < 0.05, Figure 3E). On 120 DAS, qP was higher in WW M+ plants than in WS M-(p < 0.01), WW M-(p < 0.05) and WS M+ ones (p < 0.05, Figure 3E). NPQ was statistically different between M+ and M-plants since 103 DAS ( Figure 3F). On 114 DAS, NPQ started to be significantly affected by drought in M+ plants, resulting in an increase of NPQ of WS M+ plants (p < 0.05, Figure  3F). On 120 DAS, also NPQ of M-plants were significantly affected by drought stress (p<0.05, Figure  3F).

Gas Exchange Measurements
The multi-chamber system ran continuously for 20 days, from 100 DAS to 120 DAS. The system recorded the air temperature entering (inlet) and leaving (outlet) the chambers. Mean daily temperature of inlet air was 22.3 °C; mean daily temperature of outlet air was 24.  A similar trend was observed for mean daily transpiration rate (E; Figure 4B). Under WW conditions, E was lower for M+ than for M-plants (3.5 and 5.3 mmol m 2 s −1 for M+ and M-, respectively, p < 0.05). Under WS condition, E started to decrease on 106 DAS in M-plants and two days later, on 108 DAS, in M+ plants ( Figure 4B). Thereafter and until harvesting, E decreased linearly by about 38% (3.2 mmol m 2 s −1 ) in M-plants and by about 30% (2.3 mmol m 2 s −1 ) in M+ plants.
Whole-canopy instantaneous WUE (wcWUEi), expressed as Pn/E ratio (μmol CO2/ mmol H2O), was higher in M+ plants (2.8 μmol CO2/mmol H2O) than in M-ones (1.3 μmol CO2/mmol H2O) throughout the experiment ( Figure 4C). The difference in wcWUEi became statistically different between water treatments on 110 DAS for M+ plants and on 112 DAS for M-ones ( Figure 4C, p < 0.05). On average wcWUEi under WS condition was 1.6 and 3 μmol CO2/mmol H2O for M-and M+ plants, respectively. Figure 5 reports the diurnal variation of Pn, E and wcWUEi of four representative moments of the dry-down experiment: before the water stress occurred (0.77 FTSW for WS plants, 101 DAS, Figure  5A  Before the onset of drought stress (0.77 FTSW, 101 DAS), Pn was constant throughout the day for M+ (11.1 μmol m 2 s −1 ) and M-plants (7.8 μmol m 2 s −1 , Figure 5A). Daily E was also constant throughout the day for M+ and M-plants (3.1 mmol m 2 s −1 and 4.5 mmol m 2 s −1 , respectively) ( Figure  5B). Therefore, wcWUEi was constant throughout the day without significant differences between water treatments but with significant differences between the mycorrhization treatments (p < 0.05, Figure 5C). At 0.56 FTSW, on 106 DAS, water availability affected WS M-plants, whose Pn was 28% lower than that of WW M-plants (p < 0.05, Figure 5D). At 0.56 FTSW, Pn of WS M+ plant was not different from that of WW M+ ones and it was stable during the day ( Figure 5D). E of WS M-plants decreased by about 27% compared to WW M-plants ( Figure 5E); the decrease became significant after 10 a.m. and lasted until 6 p.m. (p < 0.01). On the same DAS, E of WS M+ plant was lower but not statistically different from E of the WW M+ plants ( Figure 5E). wcWUEi at 0.56 FTSW was similar to that recorded at 0.77 FTSW for M-plants, while it was slightly higher in WS M+ plant ( Figure 5F).
At 0.34 FTSW (113 DAS), Pn of WS M-plants was 44% lower than that of WW M-plants, whilst Pn of WS M+ plants were 27% lower than that of WW M+ plants ( Figure 5G). The daily trend of Pn of M+WS plants was constant throughout the day, whilst Pn of M-WS plants significantly decreased between 11a.m. and 1p.m. (Figure 5G). E of WS M-and M+ plants was 47% and 37% lower than that of WW plants, respectively ( Figure 5H). wcWUEi of WS plants at 0.34 FTSW increased by 10% and 17% in the M-and M+ plants, respectively ( Figure 5I).
At 0.28 FTSW (120 DAS), on average Pn in WS plants was 47% lower than in WW plants both in M-and M+ plants (p < 0.01). The Pn reduction started after 2 p.m. (Figure 5J) when VPD was high (data not shown). E was affected more in M-plants (reduction by 51%) than in M+ ones (reduction by 37%) ( Figure 5K). Daily E started to decrease in M-plants after 10 a.m., whilst it started to decrease after 2 p.m. in the M+ ones ( Figure 5K). wcWUEi of WS M+ plants were statistically higher than that of WS M-ones ( Figure 5L).
The effect of drought stress was different for M+ and M-plants. FTSWt was, in fact, 0.75 and 0.70 for M-and M+ plants, respectively ( Figure 6A, full symbols). Differences between M-and M+ plants were observed also in the trend of wcWUEi in response to increasing water stress ( Figure 6B). Under WW conditions M-plants had lower value of wcWUEi compared to M+ ones. wcWUEi increased linearly with increasing water stress, with a rate that was higher for M+ plants than for M-ones ( Figure 6B). Similarly, the threshold after which wcWUEi started to increase was 0.69 and 0.75 FTSW for M-and M+ plants, respectively ( Figure 6B).  Table 2). Drought stress did not affect plant height in M+ or M-plants.
Plant stem diameter of M-and M+ plants increased in response to drought stress (p < 0.05, Table  2). Mean flowering date (defined as the date at which 50% of the plants had produced the first flower) occurred 12 days before in M+ plants than in M-ones. The number of flowering branches per plant was greater in M+ than in M-plants (p < 0.05, Table 2), whilst the number of fruit per branch was lower in M+ plants than in M-ones. In addition, fresh weight of tomato fruits was higher in M+ WW plants (25.3 g) than in M-WW ones (12.4 g). Statistically significant inoculation x drought stress interaction was observed (p < 0.05, Table 2) where drought stress increased tomato fruit production of M+ plants and reduced tomato fruit production of M-ones (Table 2). The data (mean of 3) were analyzed by two-way ANOVA (α = 0.05) with significant F-values and Pvalues highlighted in bold. In case of significant interaction, values followed by the same small letter are not significantly different (Tukey HSD). Different capital letters indicate significant differences between water treatment (WW and WS) and asterisks (second column) indicate whether inoculation caused a significant effect (Tukey HSD).
Aboveground dry matter (g/pot) was higher in M-plants than in M+ ones (p < 0.05, Table 3). Drought stress did not decrease aboveground dry matter in M-whilst reduced aboveground dry matter in WS M+ plants even though not significantly (Table 3). Belowground dry matter of M-was statistically higher than in M+ plants (p < 0.05, Table 3). No differences were induced by drought stress. Total dry matter (leaves, stems, fruits and roots) was on average higher in M-plants than in M+ plants (p < 0.05, Table 3), without any significant effect induced by drought stress. SLA was higher in M+ than in M-plants and was not affected by drought stress (Table 3). LAR was higher in M+ than in M-plants (p < 0.05, Table 3) and drought stress significantly reduced LAR both in M+ and in Mplants (p < 0.05, Table 3). Table 3. Aboveground (ABG) dry matter (g plant −1 ), belowground (BG) dry matter (g plant −1 ), root to shoot ratio, total dry matter (g plant −1 ), Specific Leaf Area (cm 2 g −1 ), and Leaf Area Ratio (cm 2 g −1 ) recorded in Exp2 for M-and M+ plants subjected to WW and WS conditions.

Variable Inoculation
Water The data (mean of 3) were analyzed by two-way ANOVA (α =0.05) with significant F-values and Pvalues highlighted in bold. In case of significant interaction, values followed by the same small letter are not significantly different (Tukey HSD). Different capital letters indicate significant differences between water treatment (WW and WS) and asterisks (second column) indicate whether inoculation caused a significant effect (Tukey HSD).
Leaf area measured during the experiment resulted significantly affected by AMF colonization (p < 0.05, Figure 7A) and by water treatment (p < 0.05, Figure 7A), but no statistically significant interaction mycorrhization x water treatment was found. M+ plants had leaf area higher than Mplants through all the experiment, both under WW and WS conditions; on the last day of measurement, water stress resulted in a reduction of 32% and 40% for M+ and M-plants respectively ( Figure 7A).
Biomass partitioning (expressed as % of total dry biomass) to leaves and stems was higher in Mplants than in M+ ones (p < 0.05, Figure 7B). On the contrary, biomass partitioning to fruits and roots was higher in M+ than in M-plants (p < 0.05, Figure 7B).

Assessment of AM Symbiosis Development
Obtaining homogeneously AMF inoculated plants is limited by several constraints such as the physical (e.i. substrate used) [47] or chemical (e.i. salinity) [48] properties of the substrate used, or the nutrients availability, especially P and N [49,50]. In our experiment the major limitation to an effective homogeneous inoculation through "a" mode of distribution was the physical structure of the products used. In fact, being Prod1 a granular inoculum, it was difficult to mix it homogeneously with peat and to distribute in the plateau. Prod2, on the other hand, being a powder formulation, was much more suitable for the homogenization with peat and the distribution in plateau. The poor symbiosis development obtained with "b" mode of distribution was a consequence of the insolubility in water of Prod1 and of the impossibility to keep Prod2 in suspension. The failure of the "c" mode of products distribution is unclear and needs further investigation; probably the competition with soil microbial community, or adverse weather conditions at transplant in open field, or roots structure, prevented AMF symbiosis.
Once the modality of inoculation was assessed, we investigated the economic feasibility of tomato seedlings production with a frequency of mycorrhization (F%) above 30% at 40 DAS, that is an average lifetime of tomato seedlings when they are sold also to the farmers and transplanted in the field. Considering the cost of Prod1 and Prod2 per seedling and an average number of tomato plant m −2 in the field, the use of mycorrhizal tomato seedlings with Prod2 at 1% will increase by 12% the actual cost per hectare of tomato seedlings (Table 1). As proved by Exp2, AMF symbiosis affected tomato physiology decreasing the vegetative period and increasing of 8.5% the tomato fruit production. All these results benefit all the stakeholders of tomato industry from production to cultivation and processing. In fact, the nurseryman with the minimum increase of 12% of the cost would easily sell on the market a certified AMF tomato seedling to the farmer; the farmer, on the other hand, will experiment the benefit of colonization and the improved water use efficiency as well as the shortage of the vegetative period and the increase in tomato fresh fruit production. The reduction of vegetative stage in AMF tomato led to an extension of harvesting period and consequently more time for tomato processing plants. Of course, further investigations in open field conditions are needed to definitively assess the hypothesized increased in tomato yield.

Ecophysiological Assessment of AM Symbiosis in Tomato
Beneficial effects of AM symbiosis on tomato plants are reported in literature [15,27]. In addition, to the expected effects on roots [48,51,52], AMF colonization contributes to modify gene expression in tomato fruits [53,54] and leaves, in particular genes involved in defense process, transport and hormonal metabolisms [55,56] as well as in lignin synthesis [15]. As previously mentioned, the beneficial effects of AM symbiosis observed during this experiment are the basis of an economic feasibility study aiming at the promotion of production and commercialization of mycorrhizal tomato plantlets.
The Eco physiological beneficial effects were all related to the photosynthetic performances of AM tomato plants compared to the non-inoculated ones during a progressive drought stress. In particular, the main effects were related to the efficiency of PS2, chlorophyll content and gas exchange that determined higher WUE under drought stress.
ETR values were calculated for M+ and M-plants using the mean absorbance ratio of 0.84. According to the literature, this value might vary across succulent and desert species characterized by great variation in color and high reflectiveness [57], but it is not the case of S. lycopersicum leaves. In our study, the high level of chlorophyll pigments (SPAD) showed by M+ plants, which confirms results of previous research [47,58], is positively correlated with PS2 efficiency (ΦPSII, r = 0.68, p < 0.05) and with ETR measured in M+ plants (r = 0.85, p < 0.05), whilst it is negatively correlated with NPQ (r = -0.78, p < 0.05). According to [59] and [60], in fact, M+ plants are less sensitive to light stress being NPQ inhibited by mycorrhization, as confirmed by data shown in Figure 3. High SPAD values, together with high PS2 efficiency, high ETR, high Fv/Fm and low NPQ, confirm for tomato the beneficial effect of the AM symbiosis already observed in another crop, as in Phaseolus lunatus [61]. Especially under conditions of high light intensity, AMF, acting as a "sink" of C, may allocate and dispose excess products of photosynthesis in an efficient way [62], dissipating light in excess avoiding permanent damages to chloroplast biochemistry and resulting in high PS2 efficiency and ETR, but low NPQ, as observed in our study (Figure 3).
In this experiment, the photosynthetic performances of M+ plants were enhanced as demonstrated by high SPAD and high Pn values. In fact, SPAD and ETR are positively correlated with Pn (r = 0.63 and r = 0.75 respectively), that overall was higher for M+ plants than for M-plants [63][64][65][66] (Figure 4). [27] reported that high Pn values in M+ plants are due to increased stomatal density induced by AM symbiosis but we did not have this experimental evidence in our experiment. Surprisingly, high Pn values of M+ plants did not lead to high total dry biomass production ( Table  3). Further investigations are therefore necessary to understand why M+ plants, that had high photosynthetic potential (high SPAD, high ΦPS2, high ETR, low NPQ) and high Pn, were shorter ( Table 2) and had lower total dry biomass than M-plants (Table 3). Despite the low aboveground biomass production that M+ plant had as consequence of the shorter vegetative period, higher total tomato fruit production than M-under drought stress conditions was observed ( Table 2).
The effect of AM symbiosis on drought tolerance is reflected on stomatal conductance and transpiration rate [67]. In our study, the effect of R. intraradices was reflected in a number of physiological responses: (i) modification of leaves structure by increasing SLA (Table 3), (ii) reduction of transpiration (E) both under WW and WS conditions ( Figure 4B), (iii) decreasing the threshold of FTSW for transpiration and increasing the threshold of FTSW for wcWUEi ( Figure 6). The reduction in transpiration rate (E) of M+ plants was double checked through the gas exchange apparatus and through gravimetric method according to [44].
The decreasing threshold of FTSW for transpiration is a consequence of the exploring ability of hyphae of AMF extraradical mycelium that provide an extensive surface area for water and nutrient absorption [68]. Among nutrients, P plays a crucial role in tomato seedling growth [69], pollen quality and quantity [70], and thereby tomato fruit production and soluble sugar concentration. In our study M+ plants, subjected to P deficiency (300 μM), had higher total tomato fruit prodution ( Table 2) and higher P content in fruits, stems and roots [37,71,72] than M-, which is in agreement with result from [53] (Figure 8). This result further proves the ability of AMF extraradical mycelium in uptaking P beyond the P depletion zone generated in the soil solution by fast root P absorption [73]. Due to the positive relationship between AMF-mediated P absorption and shoot P content enhanced by the mycorrhizosphere bacteria [74], the use of AMF in tomato cultivation would have significant implication on tomato production and quality.
AM symbiosis induced also morphological modifications to plants architecture increasing the number of flower branches per plants in line with the short vegetative period that M+ plants had. Therefore, longer reproductive period led to higher number of floral branches in M+ plants. Surprisingly, the expected number of fruits per branches of M+ plants was lower than of M-; it could be due to the adverse conditions that tomato M+ plants had during the experiment where flowering was limited by lack of wind inside the greenhouse. Further investigations are needed to better understand if AM symbiosis affect the flower differentiation or if this effect was due only to the environmental conditions during the experiment. Anyway, M+ plants produced bigger fruits than M-, leading in higher tomato fresh production especially under water stress. This is the results of the higher wcWUEi that M+ plants had during the dry-down experiment.

Conclusions
The Eco physiological analysis carried out in this study highlighted that the AM symbiosis can be established and provides beneficial effect to tomato plants using low concentration of inoculum (1%) homogeneously mixed to the substrate.
AM symbiosis induced several effects on plant physiology and phenology: (i) colonized plants (M+) reduced their vegetative period of 12 days compared to the not colonized ones (M-), (ii) were shorter than M-plants, (iii) and had higher fruit production than non-mycorrhizal plants. The high fruit production was a consequence of (iv) high Pn supported by (v) high chlorophyll content (SPAD), high PS2 efficiency, high ETR, high Fv/Fm, and low NPQ observed in M+ plants; higher wcWUEi was observed in M+ than M-as result of (vi) low transpiration rate and (vii) low threshold of FTSW of the colonized plants. Drought stress decreased shoot dry matter but did not affect tomato fruit production (viii) neither tomato fruit quality, in terms of P content (ix).
Tomato plantlets inoculated with Prod2 at 1% were 12% more expensive than non-inoculated plantlets, but their fruit production increased by about 8.5% in response to water stress and by about 27% if compared to drought stressed not inoculated tomato plants. Considering also the positive physiological effect of R. intraradices in alleviating drought stress, the use of AMF seems a promising option to stabilize tomato production and to increase resilience of tomato cultivations, as well as widening the cultivation and harvesting time of processing tomato crop in Italy and pushing farmers to use tomato M+ plants rather than the conventional ones in order to increase their earnings and to save on the cost of water and irrigation.