Impact Factor 3.677

The world's most-cited Plant Sciences journal

This article is part of the Research Topic

Biostimulants in Agriculture

Original Research ARTICLE

Front. Plant Sci., 09 October 2018 | https://doi.org/10.3389/fpls.2018.01480

The Association With Two Different Arbuscular Mycorrhizal Fungi Differently Affects Water Stress Tolerance in Tomato

  • 1Department of Life Sciences and Systems Biology, University of Turin, Turin, Italy
  • 2Council for Agricultural Research and Economics, Centre of Viticulture and Enology Research, Conegliano, Italy
  • 3National Research Council, Institute for Sustainable Plant Protection, Turin-Florence-Portici (NA) Units, Portici, Italy
  • 4CNR, Institute of Food Sciences, Avellino, Italy
  • 5Department of Health Sciences, University of Florence, Florence, Italy
  • 6Geriatric Intensive Care Unit, Experimental and Clinical Medicine Department, University of Florence, AOU Careggi, Florence, Italy

Arbuscular mycorrhizal (AM) fungi are very widespread, forming symbiotic associations with ∼80% of land plant species, including almost all crop plants. These fungi are considered of great interest for their use as biofertilizer in low-input and organic agriculture. In addition to an improvement in plant nutrition, AM fungi have been reported to enhance plant tolerance to important abiotic and biotic environmental conditions, especially to a reduced availability of resources. These features, to be exploited and applied in the field, require a thorough identification of mechanisms involved in nutrient transfer, metabolic pathways induced by single and multiple stresses, physiological and eco-physiological mechanisms resulting in improved tolerance. However, cooperation between host plants and AM fungi is often related to the specificity of symbiotic partners, the environmental conditions and the availability of resources. In this study, the impact of two AM fungal species (Funneliformis mosseae and Rhizophagus intraradices) on the water stress tolerance of a commercial tomato cultivar (San Marzano nano) has been evaluated in pots. Biometric and eco-physiological parameters have been recorded and gene expression analyses in tomato roots have been focused on plant and fungal genes involved in inorganic phosphate (Pi) uptake and transport. R. intraradices, which resulted to be more efficient than F. mosseae to improve physiological performances, was selected to assess the role of AM symbiosis on tomato plants subjected to combined stresses (moderate water stress and aphid infestation) in controlled conditions. A positive effect on the tomato indirect defense toward aphids in terms of enhanced attraction of their natural enemies was observed, in agreement with the characterization of volatile organic compound (VOC) released. In conclusion, our results offer new insights for understanding the molecular and physiological mechanisms involved in the tolerance toward water deficit as mediated by a specific AM fungus. Moreover, they open new perspectives for the exploitation of AM symbiosis to enhance crop tolerance to abiotic and biotic stresses in a scenario of global change.

Introduction

In the last decade, climatic change heavily influenced environmental conditions and a negative impact on plant development and productivity has been recorded. Plant responses to the environmental stresses include morphological, physiological, cellular and molecular changes (reviewed in Gray and Brady, 2016). One of the major abiotic threats on agriculture productivity is the progressive diffusion of water deficit in different areas of world (Dai, 2011), which induces a reduction in plant growth and development. Water deficit causes stomatal closure with a consequent decrease of CO2 adsorption followed by a reduction in photosynthetic activity and carbon partitioning (Osakabe et al., 2014; Chitarra et al., 2016). Additionally, drought has a negative impact on nutrient supply, leading to a decrease in phosphate availability (Sardans and Peñuelas, 2004). Root growth and development are important traits for plant survival in water stressed soil. A deep and extensive root system enables plants to access moisture in water-absorbed zones of soils, and is thus considered as a vital strategy for drought adaptation (Zou et al., 2017). In many crops (e.g., rice, wheat, and maize) an alteration of root system architecture has been observed under water deficit conditions, with a consistent increase in the elongation of primary roots, in the total root biomass and in the production of lateral roots, along with a change in the root growth angle and a reduction of root diameter (Comas et al., 2013). Root responses to drought stress depend on the induction of abscisic acid (ABA) pathway. This hormone is involved in different root cellular responses such as the activation of enzymes involved in cell wall remodeling (Wu et al., 1994), the regulation of cell type-specific development (Enstone et al., 2002) and the increase of lateral roots (De Smet et al., 2003). However, in agro-ecosystems, in addition to water limitation, crops have often to face the concurrent presence of other abiotic and biotic stresses, e.g., insect pests and pathogens (Bai et al., 2018). Hence, plants have to modulate a multifaceted defense by activating a number of molecular, biochemical and morphological changes, which can differ from the responses to a single stress factor (Atkinson and Urwin, 2012). Different plant defense responses lead to the release of volatile organic compounds (VOC) blends having different composition in term of presence or quantity of single compounds, depending on the stress features (Holopainen and Gershenzon, 2010). VOC released in response to pest attack can have direct and indirect defensive effect on insect performance. For example, methyl salicylate and the cis-3-hexen-1-ol, released by tomato plants in response to aphid attack (Sasso et al., 2007), have negative effects on subsequent aphid fixation (Digilio et al., 2012). These same compounds are powerful attractant of aphid natural enemies (Sasso et al., 2007, 2009). Moreover, the response of different genotypes, i.e., commercial genotypes or traditional varieties/landraces adapted to local environmental conditions, to stresses may be different. In this scenario, there is a consolidated evidence of the positive role played by “biostimulants,” e.g., organic (humic acids and seaweed extracts) and inorganic (silicium) substances, and/or microorganisms, which can be applied to improve plant nutrition/growth and to protect against stresses (du Jardin, 2015). Among them, great interest has been developed in the last years to beneficial soil microorganisms, including arbuscular mycorrhizal (AM) fungi, with a role to improve plant nutrition and tolerance to several environmental stresses. AM symbiosis is one of the most ancient interactions between the roots of more than 80% of terrestrial plants, including several crops (Balestrini and Lumini, 2018). During this interaction, in addition to the improvement of plant nutritional status, AM fungi can promote plant performance and protection from various stresses, including water deficit (Lenoir et al., 2016; Balestrini et al., 2018). In addition to an improved root system capacity to absorb nutrients and water due to the presence of external fungal hyphae, it has been also demonstrated that AM colonization influences the architecture of the host root system, leading to a better adaptation of morphology in response to water stress (Gutjahr et al., 2009; Fusconi, 2014; Zou et al., 2017). Recently, Chitarra et al. (2016) showed that the AM symbiosis positively affects the tomato tolerance to severe water deficit and how the adaptive plant response is dependent on the AM fungi species involved, underlying the importance to identify the optimal genotype/microorganism(s) combination to maximize plant resilience. Considering the role of isoprenoids in resistance to stress and plant defense, Asensio et al. (2012) tested the interaction between tomato plants and a mixed AM fungal inoculum under a drought stress conditions or jasmonic acid (JA) application, demonstrating an impact of the AM colonization in the pathways involved in isoprenoid production mainly under the considered treatments (drought and JA).

However, although the effects of AM symbiosis on the interactions with aphids have been already reported (Babikova et al., 2013, 2014), scarce information is available on these interactions under an abiotic stress condition such as water stress.

In this work, a multidisciplinary approach, involving eco-physiological, morphometric, biochemical and molecular analyses, and targeted metabolomics, has been used. Particularly, the mechanisms and the species-specificity of AM symbiosis in plant tolerance to moderate (MS) and severe (SS) water stress have been evaluated in comparison with unstressed plants. For this purpose, we tested two AM fungi, Funneliformis mosseae and Rhizophagus intraradices, in their association with a cultivated tomato cultivar (i.e., San Marzano nano). In particular, root gene expression analyses were focused on plant and fungal genes involved in phosphate (Pi) uptake and transport because of the importance of phosphorous use efficiency (PUE) in plants grown under abiotic stress (reviewed in Wang et al., 2017). It has been recently reported that, Pi transporter genes are regulated by a drought stress, both in poplar (Zhang et al., 2016) and apple (Sun et al., 2017). Here, we aimed to assess whether this applies to tomato and mainly how different AM fungi alter the expression of these genes under moderate (MS) and severe (SS) water stress conditions. In addition, considering that in the field plants are usually subjected to multiple stresses, the AM fungus (R. intraradices), resulting more efficient in increasing water stress tolerance, was selected to assess its role on tomato resistance to a combination of abiotic (water deficit) and biotic (aphid attack) stresses.

Materials and Methods

Plant, Fungal Material, and Experimental Design

Solanum lycopersicum cv. “San Marzano nano” (i.e., dwarf) seeds were surface-sterilized in sodium hypochlorite 5% (NaOCl) for 20 min, washed five times in sterile water, and germinated on wet paper. Seedlings were then moved to pots (10 × 10 × 12) containing a mixture of quartz sand (50%), sterile pumice (20%), and an inoculum (30%) of either Funneliformis mosseae (formerly Glomus mosseae; BEG 12) or Rhizophagus intraradices (FR 121), both purchased from MycAgro (Bretenière, France1). Propagules of AM fungi (i.e., a mix of spores, mycelium and mycorrhizal root pieces) are minimum equal to 10 propagules/g. For non-colonized plants (Ctrl), sterile inoculum carrier only (i.e., mix of inert mineral predominantly made with zeolite) was used instead of the specific inoculum (30%). For each pot, a quantity of about 650 g of sand plus inoculum/carrier has been used. Plants were grown in controlled conditions, at 23/21 ± 1°C (day/night), 16/8 h light/dark photoperiod, and 65 ± 10% relative humidity. From transplanting to the beginning of the experiment (after about 6/7 weeks), all the plants were watered twice a week with tap water and, once a week, with a modified Long-Ashton nutrient solution (Hewitt, 1966) containing 300 μM [Pi]. Two independent experiments (Experiment 1 and 2) were performed to determine the role of AM symbiosis under water deficit.

Experiment 1

This experiment was devoted to verify the impact of the two different AM fungi on the regulation of Pi transporter genes in roots and on the production of specific metabolites in leaves. Biometric and eco-physiological parameters have been also registered to verify the plant stressed status. Considered treatments were: (i) AM fungal colonization [non-colonized (Ctrl), F. mosseae-colonized (Fmos) and R. intraradices-colonized (Rin)] and (ii) water stress [none (NS), moderate (MS), and severe (SS)]. For each condition, 12 plants were maintained in a well-watered state (at container capacity, irrigated or unstressed, NS). The remaining 24 plants were subjected to two water stress levels. Irrigation was withheld about 6/7 weeks after fungal inoculation: moderate stress (MS) was achieved in about 2 weeks (leaf water potential of about -0.9 MPa), whilst 3 weeks were necessary to produce a severe water stress (SS) status (petiole water potential lower than -1.0 MPa). From the beginning of the water deficit, plants were moved (at least one time a week) inside the climatic chamber to avoid positional differences. Root colonization has been assessed at the beginning and at the end of the water stress imposition (not shown).

Biometric and Physiological Measurements

At the end of the water stress experiment (about 10 weeks after the beginning of the experiment), plants were harvested and plant height and shoot diameter were recorded. Parameters were taken as reported by Chitarra et al. (2016). Briefly, leaf water potential (Ψleaf) was measured on one transpiring leaf per plant, using a Scholander-type pressure chamber (Soil Moisture Equipment Corp., Santa Barbara, CA, United States). Measurements of transpiration rate (E), stomatal conductance (gs), and net photosynthetic rate (AN) were performed on adult, non-senescing, leaves at the same physiological age (in the middle part of the plant, considering the third-fourth leaf form the shoot apex). Intrinsic water use efficiency (iWUE) was calculated as the ratio between AN and gs. Measurements were taken with an infrared gas analyzer ADC-LCPro+ system (The Analytical Development Company Ltd., Hoddesdon, United Kingdom). During measurement, light intensity in the leaf chamber was set at 1200 μmol photons m-2 s-1, temperature was 25°C, and concentration of CO2 was maintained at 450–470 ppm. Measurements were taken between 10:00 and 13:00 h. The chlorophyll content index (CCI) was determined at the end of the experiment (16 or 20 days after treatment, DAT) using the portable chlorophyll meter SPAD 502 (CCM-200, Opti-Sciences, Inc., Hudson, NH, United States). Readings were collected from the second or third fully developed leaves from the top on five randomly selected tomato plants for each experimental condition.

Phosphorus Determination

Analytical grade reagents and Phosphorus (P) standard (1000 mg/L) were purchased from Sigma-Aldrich chemical. The digestion of leaf samples was performed according to Olowu et al. (2015) with some modifications. 0.10 g of each sample was weighed and 3 mL of a mixture of nitric acid: perchloric (1: 1; v: v) was added to each. The mixture was placed in a digester at 220°C for about 60 min. After 45 min, 500 μL of H2O2 (30%) were added in order to facilitate the oxidation reaction. After cooling the whole was transferred into a 10 mL volumetric flask and brought to volume with H2O. An Inductively Coupled Argon Plasma Optical Emission Spectrometers (ICP-OES iCAP 7000 Series Thermo Scientific), equipped with ASX-520 Autosampler (CETACTM, Thermo Scientific) was used. To prepare the calibration curve with the Phosphorus concentrations equal to 0.5; 1; 2; 4 and 5 mg/L, the Multi-Element Test Solution ICAP 6000 with P concentration of 10 mg/L in 1% v/v HNO3 was used. The element content was calculated by using standard curves and the final concentrations of samples were expressed as g/kg dry weight of P. The emission line for the analysis by ICP OES was chosen according to previous interference studies. The line that exhibited low interference and high analytical signal and background ratios was selected. The emission line that was employed was 185.942 nm. The ICP-OES measurements were performed in triplicate, with R2 = 0.9994 and BEC (Background Equivalent Concentration) = 0.002 ppm.

RNA Extraction and RT-qPCR

Expression changes of target transcripts were quantified on root samples (three independent biological replicates) by quantitative real-time PCR (RT-qPCR). Roots from two plants from each treatment were pooled to form a biological replicate, immediately frozen in liquid nitrogen and stored at -80°C. Total RNA was isolated from each biological replicate and cDNA synthesis was performed as described in Chitarra et al. (2016). Genomic DNA contamination was checked before proceeding with the cDNA synthesis by PCR reactions using LeEF specific primers of tomato (Supplementary Table S1). RT-qPCR experiments were carried out in a final volume of 15 μl containing 7.5 μl of Rotor-GeneTM SYBR® Green Master Mix (Qiagen), 1 μl of 3 μM specific primers and about 10 ng of cDNA. Samples were run in the Rotor Gene apparatus (Qiagen) using the following program: 10 min pre-incubation at 95°C, followed by 40 cycles of 15 s at 95°C, and 30 s at 60°C. Each amplification was followed by melting curve analysis (60–94°C) with a heating rate of 0.5°C every 15 s. All reactions were performed with three technical replicates and only Ct values with a standard deviation that did not exceed 0.3 were considered. The comparative threshold cycle method (Rasmussen, 2001) was used to calculate relative expression levels using plant LeEF and LeUBI reference genes. Oligonucleotide sequences are listed in Supplementary Table S1. In detail, genes encoding for plant (LePT1, LePT2, LePT3, LePT4, and LePT5) and fungal (RiPT and FmPT) phosphate transporters have been considered.

Targeted Analysis of Leaf Volatiles

Leaf VOC analyses (isolation, identification and quantification) were done by headspace solid-phase microextraction-gas chromatography-mass spectrometry (HS-SPME-GC-MS). For sample preparation, 0.1 g of tomato leaf stored in a 2 ml screw cap headspace glass vial at -80°C were transferred at room temperature for 5 min and 0.5 ml of 30% NaCl solution containing the internal standard (2-hexanone at 1 ng/μl) was added to the vial. After mixing, samples were pre-incubated for 30 min at room temperature. A 65 μm PDMS/DVB SPME fiber (Supelco, Sigma Aldrich, Milan, Italy) was used for extraction procedure: the fiber was introduced into the vial headspace for 5 min.

A Hewlett Packard GC-MS system composed by a HP 5890 Series II Gas Chromatograph coupled to a HP 5971A Mass Selective Detector single quadrupole mass spectrometer was used. The GC-MS was equipped with a ZB-5MS column with column guard (Phenomenex) (30 m × 0.25 μm i.d.; film thickness 0.50 μm). Oven temperature was set up with following parameters: initial temperature 45°C maintained for 3 min, then to 220°C at 12°C/min, held isothermally at 220°C for 2.5 min. The volatiles trapped on the fiber were desorbed for 2 min at 250°C in the injection port of the GC; the injection was performed in splitless mode for 1 min. To prevent cross-contamination between successive samples, the fiber was cleaned by exposure in a different GC injection port for 5 min at 250°C before a new sampling process. Ultra high purity helium was used as the carrier gas at a constant head pressure of 12 psi (corresponding to 1 ml/min at 45°C). Transfer line temperature was 280°C.

Data were recorded in full scan mode from 40 to 550 m/z, using an electron ionization (EI, 70 eV) source, with 1.53 scan/s. Data acquisition and processing were performed using HP ChemStation software (version D.02.00).

Quantitative analysis by GC-MS was performed recording the signals in selected ion monitoring (SIM) mode; the areas of the quantifier ion peak of each analyte were measured and compared to that of the internal standard 2-hexanone. In Supplementary Table S2 the retention time and principal ion of each compound are reported.

Calibration curves constructed with pure standards – trans-2-hexenal, methyl salicylate, eugenol, α–phellandrene, β–phellandrene (Sigma Aldrich) – in the ranges of the interest concentrations allowed the calculation of the compound concentrations in the samples.

Experiment 2

Experiment 2 was devoted to evaluate the impact of AM symbiosis on a combined abiotic (moderate water stress) and biotic (aphid infestation) stress condition. To this aim, aphid survival (a measure of direct plant defense), parasitoid attractiveness (a measure of indirect defense) and VOC emissions were considered in unstressed and water stressed conditions. Considered treatments were: (i) AM fungus [non-colonized (Ctrl) and R. intraradices-colonized (Rin)], (ii) water stress [none (NS) and moderate (WS)], and (iii) Aphid infestation (Aph) and their combinations. In detail, for aphid survival, 20 plants non-colonized (Ctrl) and 20 colonized with R. intraradices (Rin) were prepared for the aphid survival experiment: 10 plants for each condition were maintained under a well water conditions (CtrlNS and RinNS) and 10 were subjected to a moderate water stress (WS and WS + Rin). By contrast, for parasitoid attractiveness and VOC collection, 6 plants for the following 8 conditions were considered: (NS, Aph, Rin, WS, Rin + Aph, WS + Rin, WS + Aph, and WS + Rin + Aph), totalling 48 plants.

Insect Rearings

Both the herbivore and its parasitoid were reared at CNR-Institute for Sustainable Plant Protection (Portici, NA, Italy) as follows:

– the potato aphid Macrosiphum euphorbiae Thomas (Hemiptera: Aphididae) was reared on tomato plants cv. San Marzano nano (i.e., dwarf) from material collected in the field (Scafati, SA, Italy) on the same plant cultivar in 2001 and periodically refreshed by adding field collected aphids. Rearing conditions were: 20 ± 1°C, 18 h light/6 h dark photoperiod, 65 ± 5% relative humidity.

– the aphid parasitoid Aphidius ervi Haliday (Hymenoptera: Braconidae) was permanently reared on the pea aphid Acyrthosiphon pisum reproduced on broad bean plants (cv. Aquadulce) from material collected in the field (Battipaglia, SA, Italy) in 2001 on alfalfa and periodically refreshed by adding field collected aphids. Rearing conditions were: 20 ± 1°C, 18/6 h light/dark photoperiod, 65 ± 5% relative humidity (Guerrieri et al., 2002).

Aphid Survival Assays

For each experimental condition 10 plants were infested with one newly born nymph of M. euphorbiae, transferred with a soft brush. All plants were checked daily to assess presence of aphids, presence of exuviae (evidence of molting), presence and number of newly laid nymphs and dead aphids. When reproduction started, the offspring were counted and removed daily. The mortality of the aphids was recorded daily until the death of the last individual. The parameters of the bioassay were set as indicated above.

Wind-Tunnel Bioassay

Plants were tested in a wind-tunnel bioassay for their attractiveness toward the aphid parasitoid wasp A. ervi. For the aphid infestation fifty mixed-aged individuals of M. euphorbiae (mimicking a natural population) were gently transferred by paintbrush on tomato leaves allowing them to feed for 1 week. For the double stress thesis, the aphids were added at the onset of water deficit imposition. These treatments were performed using both inoculated with R. intraradices and non-inoculated tomato plants. For each experimental condition, a total of 6 plants was used and offered individually every day for 6 consecutive days in a random order to reduce any bias related to the time of the experiments. One hundred parasitoid females were tested singly for each target in no-choice experiments, and observed for a maximum of 5 min. The percentage of response (oriented flights, landings on the target) to each target plant was calculated. The parameters of the bioassay were set as follows: 20 ± 1°C; 65 ± 5% relative humidity; 25 ± 5 cm s -1 wind speed; 50 cm distance between releasing vial and plant target.

VOC Collection and Analysis

Immediately after the bioassay in the wind tunnel, each plant was used for the collection of head-space volatiles by an airtight entrainment system. The system consisted of a bell jar containing a single potted plant and the jar was connected to a circulation pump forcing the air in an adsorbent trap made of Tenax TA, 60–80 mesh (Sigma-Aldrich). Each collection lasted 3 h under controlled conditions (24 ± 2°C, 18/6 h light/dark photoperiod, 70 ± 10% relative humidity, 700 μmol m 2 s 1 PPFD). Collected volatiles were analyzed by an integrated system including thermal desorber (Tekmar TD-800) mounted on a gas chromatograph (column: RTX-200, 60 m, 0.25 mm ID, 0.25 μm, carrier gas: He) coupled to a mass spectrometer detector. The resulting peaks were compared with a compound database library (National Institute of Standards and Technology) and the following available authentic standards: anisole-p-allyl, camphor, 3-carene, (E)-β-caryophyllene, chlorobenzene, α-copaene, α-cubebene, p-cymene, decane, p-dichlorobenzene, dodecene, eucalyptol, eugenol, (E)-β-farnesene, α-gurjunene, hexanal, (Z)-3-hexen-1-ol, humulene (=α-caryophyllene), (R)-(+)-limonene, (S)-(-)-limonene, linalool, (+)-longifolene, menthol, 6-methyl-5-hepten-2-one, methyl salicylate, b-myrcene, (Z)-nerolidol,(E)-β-ocimene, (R)-(-)-α-phellandrene, α-pinene, skatol, α-terpinene, γ-terpinene, α-terpineol, and terpinolene.

Statistical Analysis

Analysis of variance (ANOVA) of the experimental data was performed using the SPSS software. When ANOVA indicated that either stress (stress: NS, MS, and SS) or mycorrhizal colonization (myco: Ctrl, Fmos, and Rin) or their interaction was significant, mean separation was performed using the Tukey HSD test, adopting a probability level of P < 0.05. The volatile emission patterns, measured as peak areas divided by fresh plant weight, were analyzed by Detrended Correspondence Analysis (DCA). Data were square root transformed to meet normal distribution. Analysis was conducted using the R statistics programming environment 3.3.3 (R Development Core Team, 2012) and the ‘vegan’ package ver. 2.5-2. Multivariate analysis of variance by permutation (PMANOVA) was used to test for significant differences between group centroids using vegan’s adonis functions (Oksanen et al., 2013). The number of parasitoids responding to each target was compared by a G-test for independence with William’s correction. The resulting values of G were compared with the critical values of χ2 (Sokal and Rohlf, 1995). To analyze the effect of each experimental condition on aphid survival, we used function survfit and Kaplan–Meier plot to visualize the survival curves, contained in the R package ‘survival’ version 2.41-3.

Results

Biometric Parameters and Ecophysiology

In our experiments, biometric parameters including plant height, the internodes/height ratio and shoot diameter have been affected by both water treatments and AM colonization, while no significant differences were registered for the chlorophyll content index (CCI; Supplementary Figures S1A–D). Plant height and the internodes/height ratio were influenced by stress (S) and myco (M) treatments with no significant S × M interaction, although Fmos-inoculated plants under MS and SS conditions showed higher height when compared with RinMS/SS and CtrlSS plants (Supplementary Figure S1A and Table S3). Considering the internodes/height ratio, highest values were recorded for RinMS and SS plants with respect to other treatments (Supplementary Figure S1B and Table S3). Conversely, shoot diameter did not show a different trend among the M variables, while a significant impact of the S variables and S × M interaction was observed (Supplementary Figure S1C).

The imposed water stress conditions differently affected leaf water potential (Ψleaf) levels independently from controls and the AM inoculum used. At the end of the experiment, significant differences were observed. NS plants showed in fact values lower than/equal to -0.5 MPa, while Ψleaf values of about -0.9 MPa and -1.4 MPa were observed in MS and SS plants, respectively (Figure 1A).

FIGURE 1
www.frontiersin.org

FIGURE 1. Leaf water potential and gas exchanges in AM– and AM+ plants unstressed (NS) and subjected to two different water deficit levels (moderate stress, MS; severe stress, SS), measured at the end of the experiment. (A) Leaf water potential (Ψleaf); (B) net photosynthetic rate (AN); (C) stomatal conductance (gs); and (D) iWUE levels in AM– (Ctrl) and AM+ (Rin and Fmos) plants. All data are expressed as mean ± SE (n = 5). ns, ∗, ∗∗, and ∗∗∗: non-significant or significant at P ≤ 0.05, P ≤ 0.01, and P ≤ 0.001, respectively. Different letters above the bars indicate significant differences according to Tukey HSD test (P ≤ 0.05), considering S × M interaction. Analysis of variance on the single variables is reported in Supplementary Table S3.

Considering AN rates, a decreasing trend was observed from NS to MS and SS conditions. Significant higher AN rates were recorded in NS condition, particularly in RinNS plants. Conversely, significant lower values were measured under SS condition with the exception of RinSS plants that were statistically similar to FmosMS ones (Figure 1B). A similar trend was also evident for gs, although a significant influence was observed only for S and M variables (Figure 1C and Supplementary Table S3). In addition, iWUE index, irrespectively of treatment or water status, showed similar values, with the exception of RinSS plants where significant higher iWUE was observed (Figure 1D).

Phosphorus Determination and Expression of Phosphate Transporter Genes

In no stress condition, no significant differences were detected among the three treatments in phosphorous (P) determination in leaves, i.e., CtrlNS, RinNS, and FmosNS plants, although higher values have been measured in AM-colonized plants in respect to non-colonized ones, with the highest P content (g/kg) recorded in F. mosseae-colonized plants. While no strong differences have been registered under a moderate stress condition (MS), under a severe water stress condition (SS) a higher leaf P content was detected in F. mosseae-colonized plants mainly respect to Ctrl plants (data not shown).

Additionally, the expression of plant and fungal Pi transporter genes was evaluated in roots from plants at different growth conditions (no stress and water deficit). Expression of PT genes was not influenced by M x S interaction. In detail, although statistical analysis showed an effect of the S variable only for LePT2, both LePT1 and LePT2 changed their expression under water deficit, but with an opposite trend (Figures 2A,B and Supplementary Table S3). LePT1 expression in fact increased under water deficit with the highest value in RinSS, while LePT2 transcripts progressively decreased with the water deficit. Among the mycorrhizal-inducible PT genes (LePT3, LePT4, and LePT5), LePT3 resulted to be overexpressed in RinNS roots, while decreased under water stress mainly in AM-colonized plants (Figure 2C). LePT4 resulted to be influenced by S and M, while LePT5 only by S (Figures 2D,E and Supplementary Table S3). In detail, these two last transporters showed a similar expression trend, with increased levels in SS plants (Figures 2D,E) with respect to the other treatments (NS and MS). However, LePT4 resulted overexpressed in RinSS plants, while LePT5 was in FmosSS ones, suggesting a species-specificity effect. From the fungal side, F. mosseae PT gene showed lower expression levels than R. intraradices PT in all the considered conditions (Figures 3A,B). Under water deficit, the expression of these fungal genes was not significantly affected, except for RiPT that is significantly up-regulated in MS (Figure 3B).

FIGURE 2
www.frontiersin.org

FIGURE 2. Expression changes of tomato phosphate transporter (PT) genes in roots of AM– and AM+ plants upon irrigation (NS), moderate water stress (MS), and severe water stress (SS) conditions. RT-qPCR analysis of (A) LePT1, (B) LePT2, (C) LePT3, (D) LePT4, and (E) LePT5. Lower letters denotes significant differences attested by Tukey HSD test (P < 0.05). All data are expressed as mean ± SE (n = 3). ns, ∗, ∗∗, ∗∗∗: non-significant or significant at P ≤ 0.05, P ≤ 0.01, and P ≤ 0.001, respectively. Different letters above the bars indicate significant differences according to Tukey HSD test (P ≤ 0.05), considering S × M interaction. Analysis of variance on the single variables is reported in Supplementary Table S3.

FIGURE 3
www.frontiersin.org

FIGURE 3. Expression changes of fungal phosphate transporter (PT) genes in AM+ roots upon irrigation (NS), moderate water stress (MS), and severe water stress (SS) conditions. RT-qPCR analysis of (A) F. mosseae PT (FmPT) and (B) R. intraradices PT (RiPT) genes, respectively. Lower letters denotes significant differences attested by Tukey HSD test (P < 0.05). Data are expressed as mean ± SE (n = 3).

Target Metabolomics

The target leaf VOCs present in NS and SS plants in the presence or in the absence of either F. mosseae or R. intraradices have been determined and quantified. The main VOCs, trans-2-hexenale, methyl salicylate, eugenol, α-phellandrene and β-phellandrene increased in SS, mainly in the presence of the AM fungus F. mosseae, in agreement with statistical analysis that showed a highly significant influence of both S and M variables for all the considered VOCs (Table 1). Differently, the S × M interaction resulted to be significant for all the metabolites with the exception of trans-2-hexenal (Table 1).

TABLE 1
www.frontiersin.org

TABLE 1. Quantification of VOC products in tomato leaves.

Insect Performance and VOC Emission

Root colonization by R. intraradices significantly improved the development of Macrosiphum euphorbiae on tomato in both conditions tested (NS and WS). When aphid survival curves were calculated, survival probability proved to be higher on colonized plants (Figure 4, black lines, log-rank: colonized vs. non-colonized plants, P = 0.015). Conversely, water stress (dashed lines) did not affect the development of the aphids (Figure 4, log-rank: WS vs. NS, P = 0.87).

FIGURE 4
www.frontiersin.org

FIGURE 4. Probability of aphid Macrosiphum euphorbiae survival when fed on tomato AM– and AM+ plants (inoculated with R. intraradices), upon irrigation (NS and Rin) and water-stress (WS and WS + Rin) conditions.

Root colonization by R. intraradices affected the foraging behavior of A. ervi in almost all conditions tested (Figure 5). Oriented flights significantly increased by root colonization in control, aphid infested and double stressed colonized plants. The same pattern was recorded for the landings on the source with the exception of control plants where root colonization did not affect parasitoid response. The presence of R. intraradices did not enhance the foraging behavior of A. ervi in condition of water stress.

FIGURE 5
www.frontiersin.org

FIGURE 5. Flight behavior (oriented flights and landing on source) of Aphidius ervi females (%) toward tomato plants inoculated with R. intraradices, infested by aphids, their combination and relative control. Asterisks assigned by G-test for independence (P < 0.05 and ∗∗∗P < 0.001).

Several of the collected volatiles were significantly affected by the presence of R. intraradices mainly in stressed plants (Table 2, Figure 6, and Supplementary Table S4). The DCA analysis separate plots on the basis of different treatments (Figure 6 and Supplementary Table S4). Distance between centroids (Supplementary Table S4 and Figure 6) indicated that volatile blends differ significantly in all treatments expect for (WS vs. Aph) and in almost all their combinations (double and triple, Supplementary Table S4). Whilst in non-stressed plants an increase was recorded in colonized plants in respect to non-colonized ones, in double stressed plants the AM fungus affected the VOCs emission with a significant decrease of most of volatiles, with the exception of methyl salicylate emission that increased under double stress conditions with respect to all other treatments (Table 2). Conversely, three terpenes (α-terpinene, α-phellandrene and β-phellandrene) were released at about half rate by R. intraradices colonized plants under double stress in respect to non-colonized ones.

TABLE 2
www.frontiersin.org

TABLE 2. VOCs released by tomato plants that were not stressed (NS), aphid-attacked (Aph), water-stressed (WS), and exposed to both water stress and aphid infestation (WS + Aph).

FIGURE 6
www.frontiersin.org

FIGURE 6. Detrended Correspondence Analysis (DCA) score (A) and loading (B) plots of the first two components based on volatiles emitted tomato plants that were not stressed (NS), aphid-attacked (Aph), water-stressed (WS), and exposed to both water stress and aphid infestation (WS + Aph). (A) Each experiment was carried out using uninoculated (NS, Aph, WS, and WS + Aph) and inoculated with R. intraradices (Rin, Rin + Aph, WS + Rin, and WS + Rin + Aph) plants. In brackets the percentage of variation explained is indicated. Multivariate analysis of variance by permutation was used to test for significant differences between group centroids (see Supplementary Table S4). (B) Loading plot of the DCA first two components, showing the contribution of each VOC (C1–C31, listed in Table 2) to the model.

Discussion

Although several papers were dedicated to describe the effect of AM symbiosis on tomato response to water deficit (Dell’Amico et al., 2002; Subramanian et al., 2006; Aroca et al., 2008; Wang et al., 2014; Chitarra et al., 2016; Ruiz-Lozano et al., 2016; Rivero et al., 2018), here we have focused on some specific aspects less considered so far. In this study, we have examined the impact of AM symbiosis on the tomato response to water deficit considering some belowground traits related to symbiosis, i.e., the expression of Pi transporter genes, and some aboveground responses, in terms of plant performance traits and changes in metabolites having a role in tolerance/defense. Two different AM fungal species have been used, and results here presented confirmed a species-specific impact on belowground/aboveground interactions in tomato. Previous works already showed the species-specificity of the interactions between plants and AM fungi under stressed conditions (Chitarra et al., 2016; Quiroga et al., 2017; Pollastri et al., 2018). More recently, Rivero et al. (2018) reported an untargeted metabolomic analysis in tomato roots colonized by three AM fungi of different genera, to verify their impact on tomato tolerance to drought or salt stress. Plant growth responses were also considered and, overall, results showed that whilst some responses were common to all AM fungi tested, others were specifically related to single isolates. The common effect of all AM fungal species tested was a higher biomass (both above- and belowground) in respect to non-colonized plants, correlated with the degree of the stress. Here, tomato plants inoculated with R. intraradices showed higher internodes/height ratios, calculated as a measure of plant growth, under water deficit with respect to the other conditions. This result suggests a positive influence of this isolates in regimes of water deficit, probably due to a more compact plant architecture, and thus less subject to water dispersion. As confirmed by statistical analysis (Supplementary Table S3), no significant impact on the considered growth parameters has been recorded in tomato plants colonized by F. mosseae. The discrepancy with the paper by Rivero et al. (2018) could be due to the different growth conditions and most probably to a different tomato cultivar. A different impact of the two fungi on plant performance under water deficit can be also highlighted, with R. intraradices that leads to higher iWUE, as already reported in Chitarra et al. (2016). Overall, these results confirmed the importance to produce ad hoc formulates of AM fungi in relation to a specific plant species/cultivar and environmental condition, for a better exploitation/protection in the field (Berruti et al., 2013, 2015).

Two AM Fungal Species Differently Affect Pi Uptake

Limited water resources and increasing soil salinity affect crop performance and productivity worldwide. These abiotic factors can be also accompanied by the low availability of many mineral nutrients including phosphorous (P) (Grattan and Grieve, 1999) and by the inhibition of mineral nutrient uptake and translocation (Wang et al., 2017). P fertilization can increase stress tolerance and productivity in several plant species and it appears mandatory to understand the mechanism/s of Pi acquisition and utilization under salinity and drought, to improve PUE and consequently crop productivity under abiotic stresses (Sawers et al., 2017; Wang et al., 2017). AM symbiosis promotes an improved nutritional status of the host plant (Smith and Smith, 2012), particularly under low nutrient availability. The success in nature of this symbiosis, known for more than 80% terrestrial plants, stays on the bidirectional nutritional exchanges, where the fungus supplies the plants with mineral nutrients and, in turn, receives carbon compounds (Balestrini and Lumini, 2018). Pi uptake, transfer and delivery have been largely investigated in AM roots, resulting in the characterization of a symbiotic Pi uptake pathway (Bucher, 2007; Javot et al., 2007; Smith and Smith, 2011). Particularly, PHT1 transporters, which are mainly involved in Pi uptake from soil and translocation, have been reported to be involved in the acquisition of Pi from AM (Glassop et al., 2005; Javot et al., 2007; Tamura et al., 2012; Yang et al., 2012; Duan et al., 2015; Liu et al., 2016; Tian et al., 2017; Liu et al., 2018). In tomato, eight PHT1 genes are expressed in the roots (Chen et al., 2014) and mycorrhizal-induced PT genes have been identified (Nagy et al., 2005, 2009). Different PT genes can be differentially regulated by drought. For example, several genes encoding Pi transporters from poplar and apple were found to be regulated by drought stress and Pi levels. These genes, especially those up-regulated by drought stress at low Pi level, might contribute to drought tolerance of crops in Pi-limited soils (Zhang et al., 2016).

Here, new information on the expression of five tomato PT genes in presence/absence of the AM fungus and in response to two different water deficit conditions (moderate and severe) is reported. As reported in Chitarra et al. (2016), root colonization seemed not to be strongly affected by a severe water stress in this tomato cultivar, independently of the AM fungal species considered. It is worth noting that the proportion of AM root colonization was greater for R. intraradices, while the percentage of arbuscules in the colonized portion was significantly greater for F. mosseae, both in NS and SS conditions, suggesting a different colonization strategy adopted by the two fungi. Looking at the regulation of PT genes, a significant up-regulation for the AM-inducible LePT3, LePT4, and LePT5 has been observed in well-watered conditions, in agreement with previous data (Nagy et al., 2005, 2009; Balestrini et al., 2007). Under water stress, our results showed a different regulation of the considered PT genes. Particularly, it is worth noting the opposite trend for the two genes involved in the direct Pi uptake from soil (LePT1 and LePT2), independently from the presence of the AM fungus. The poplar PHT1;2 has been already reported as induced by drought stress independently from Pi level (Zhang et al., 2016). Among the AM-inducible PT genes, LePT3 seems not involved in response to water deficit at least in non-colonized plants as indicated by the lowest transcript level recorded in NS, whilst in colonized plants its expression fell independently by the stress level. LePT4 and LePT5 transcript levels both increased under water deficit with highest values in FmosSS, in agreement with the higher P content found in the leaves of these plants in respect to control ones. Overall results suggest a major role played by LePT2 and LePT4 in promoting tolerance to water deficit, particularly under a severe condition, although this is dependent from the fungal species. Since PT gene induction is also related to Pi-starvation, it is important to note that, in our experimental system, MS plants still received some nutrients during watering, while the nutritional status of SS plants did not received any external inputs even if the P content in leaves does not show a starvation level. We can then speculate that the transcript level for LePT5 is more affected by a lower Pi content in soil than to water deficit. In fact, it has been already reported that the expression of some PT genes under drought can be Pi-dependent (Zhang et al., 2016). However, we cannot exclude a relation between our results on PT gene expression and the root architecture modification, which might happen under the different treatments. For example, it has been reported that in legumes phosphate starvation induced a higher number of lateral roots and PT4 genes, homologs of LePT4, can play a role in root tips, creating a link among Pi-perception, root branching, and Pi-signaling mechanisms (Volpe et al., 2013, 2016).

The different behavior of the two AM fungal species tested further highlights the species-specificity of these complex interactions, regulating their final outcome, depending on both fungal and plant species considered (Sawers et al., 2017).

AM Symbiosis Impacts Aphid Survivorship and VOCs Under Water Deficit and Tomato Response to Multiple Stresses

Root-associated microorganisms can alter the development of phytophagous insects in different ways. They can promote defensive responses that hamper the subsequent development of the invasive insect (Guerrieri et al., 2004) or, conversely, by improving plant quality that can enhance plant infestation (Battaglia et al., 2013). In our experiment, R. intraradices promoted the survival of M. euphorbiae in both control and water stressed plant. The null impact of water stress on aphid survival in R. intraradices colonized plants indirectly demonstrates the role of this AM species on plant tolerance to water deficit. Indeed, while water deficit should have had a deep impact on the survival of the sap feeding insects, by reducing the feeding time on sieve elements (Nachappa et al., 2016), this was not recorded for the aphid M. euphorbiae in our tests.

Plants synthetize and emit a large variety of VOC that function as important protective and signaling molecules (Guerrieri, 2016). Both biotic and abiotic stresses can induce plant volatile emission (Loreto and Schnitzler, 2010; Dong et al., 2016; Catola et al., 2018), including emissions of lipoxygenase (LOX) pathway products (various C6 aldehydes such trans-2-hexenale), shikimate pathway products (e.g., methyl salicylate), specific mono- and sesquiterpenes (such as β-phellandrene) and methanol (Niinemets, 2010; Dong et al., 2016). Additionally, the rate of emissions induced by a specific stress depends on genotype stress tolerance, timing, duration and severity of the stress (Niinemets, 2010; Copolovici et al., 2011; Niederbacher et al., 2015). Our target metabolomics analysis was focused on the impact of AM symbiosis on plant tolerance to a severe water stress, to evaluate an extreme condition that has been still less considered. The data underlined that the SS has an impact on the VOC production, demonstrating that plants, both non-colonized and F. mosseae-colonized, showed the highest concentration of the main VOC products with respect to the NS plants. Furthermore, a different effect of the AM fungal species was again evident, with R. intraradices-inoculated plants that showed similar level of all targeted metabolites in comparison with NS plants, suggesting an enhanced tolerance to this type of stress promoted by this AM fungal species. In tomato, a vast array of VOCs was found to be emitted under biotic (Digilio et al., 2010) and abiotic stresses, including water deficit (Catola et al., 2018), flooding (Copolovici and Niinemets, 2010), cold and heat (Copolovici et al., 2012; Kask et al., 2016). In our experimental study, where necessarily a MS has been used to permit the combined test with aphids, the emission of a typical tomato VOC bouquet under NS and water deficit conditions was determined. A change in VOC emission was particularly observed between non-colonized and R. intraradices-colonized plants both in well-watered condition and under a double stress (water deficit plus the aphid attack). However, the trend emission was opposite between NS and double stress plants. In NS plants colonized by R. intraradices, an increase in the release of VOC has been recorded particularly for the methyl salicylate. Significant changes in VOC emission has been registered for the single stress conditions in AM-colonized plants, while the presence of the AM fungus in double stressed plants lead to a decreasing trend in the emission of several compounds with one exception: methyl salicylate. This result strongly supported the behavioral observations given that methyl salicylate plays a major role in the attraction of A. ervi (Sasso et al., 2007), being detected by this parasitoid at antennal level at a concentration as low as 0.01 gr/ml (Sasso et al., 2009). On the other hand, the strong reduction in the release of the terpenes in double stressed plants colonized by R. intraradices in respect to non-colonized ones could be due to a shift toward the production of non-volatile compounds (essential isoprenoids such as carotenoids) in respect to VOC biosynthesis (not essential isoprenoids). This has been already reported by Asensio et al. (2012) that found that water stressed colonized plants plus exogenous jasmonic acid (JA) treatment (mimicking a wound) have a higher level of carotenoids, which have an important role in photosynthesis, with respect to non-colonized plants. The reduction of monoterpenes recorded in our experiment is also in line with the data of these same authors that explained such a reduction with the high demand of carbon by the AM fungus, particularly under combined stress (drought + JA treatment) in respect to unstressed and non-colonized plant (Asensio et al., 2012). However, a correlation between the VOC emission in colonized plants under double stress and a possible better P availability in the same plants, cannot be excluded, considering that nutritional status might have an impact on plant isoprenoid production (Asensio et al., 2012).

The effect of R. intraradices on plant resistance toward aphids resulted bifaceted. AM-colonized plants represented a better substrate for the survival aphid M. euphorbiae in respect to non-colonized ones thus reducing plant direct resistance toward insect pests. This negative effect on plant fitness was largely compensated by the significant enhancement of indirect defense as represented by the significant higher attraction of the parasitoid A. ervi recorded for colonized plants in respect to non-colonized ones. The increased level of indirect defense persisted in R. intraradices colonized plants also in conditions of double stress and was supported by the release of key attractive compounds (methyl salicylate).

Conclusion

Our results suggest that the two AM fungi can trigger different adaptation strategies against environmental stresses, with F. mosseae that seems more effective on VOC production and R. intraradices on the considered plant performance traits, e.g., leading to a significant higher water use efficiency under a severe water stress. Additionally, R. intraradices was demonstrated to be effective on response to combined abiotic and biotic stress, the latter in terms of attractiveness toward the natural enemies of aphids. However, although AM fungi enhance crop tolerance to environmental stresses, the exploitation and the practical application in the field of these symbionts require a thorough identification of the mechanisms involved in nutrient transfer, the metabolic pathways induced by single and multiple stresses and the physiological mechanisms leading to improved tolerance. Our results well fit with this aim, adding new tiles to the pieces of information already present about the mechanisms involved in the tolerance to water and biotic stresses in tomato during AM symbiosis, and in nutrition aspects such as Pi uptake. Overall, our results confirmed the importance to develop mixed AM-based products, depending on the environmental conditions to be faced, choosing the right symbiotic partner/s (i.e., plant cultivar, AM fungus) that lead to different outcomes in terms of plant tolerance to stress. This is particularly urgent in a scenario of climate change as characterized by a progressive lack of water that combines, in the field, with other abiotic (e.g., nutritional availability and salinity) and biotic (e.g., insects and pathogens) stresses. However, considering the species-specificity in affecting tolerance/resistance traits, it is relevant to test individually an AM species/isolate before to produce a mixed inoculum, to avoid a null effect on plant response.

Author Contributions

RB, EG, and BM conceived and designed the experiment. VV and WC performed the molecular analyses and the eco-physiological measurements. PC performed the multiple stress experiment and VOC collection. MV performed the phosphate determination. WC and PC performed the statistical analysis for all the data. PB participated to the metabolite extraction. GM, GP, CDS, and BM performed the target metabolomics experiment. EG participated to the multiple stress experiment and analyzed the VOC data. RB prepared the plants and performed the molecular analyses with VV. RB and VV wrote the manuscript with the contribution of EG and WC. All authors read and approved the manuscript.

Funding

This research was funded by the AQUA project (Progetto Premiale CNR).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01480/full#supplementary-material

FIGURE S1 |(A–C) Biometric parameters and (D) chlorophyll content index (CCI) and of AM- and AM+ tomato plants upon unstressed and water stressed (MS and SS) conditions. All data are expressed as mean ± SD. ns, ∗, ∗∗, ∗∗∗: non-significant or significant at P ≤ 0.05, P ≤ 0.01, and P ≤ 0.001, respectively. Different letters above the bars indicate significant differences according to Tukey HSD test (P ≤ 0.05), considering S × M interaction. Analysis of variance on the single variables is reported in Supplementary Table S3.

TABLE S1 | List of the oligonucleotides used in this study.

TABLE S2 | Retention time and principal ion of each compound.

TABLE S3 | Analysis of variance on single variables referred to Figures 1, 2 and Supplementary Figure S1.

TABLE S4 | Results of pairwise multivariate analysis of variance by permutation among VOCs released by tomato plants (see Table 2 for treatments). Significant values are in bold.

Footnotes

  1. ^ http://www.mycagrolab.com/

References

Aroca, R., Del Mar Alguacil, M., Vernieri, P., and Ruiz-Lozano, J. M. (2008). Plant responses to drought stress and exogenous ABA application are modulated differently by mycorrhization in tomato and an ABA-deficient mutant (sitiens). Microb. Ecol. 56, 704–719. doi: 10.1007/s00248-008-9390-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Asensio, D., Rapparini, F., and Penuelas, J. (2012). AM fungi root colonization increases the production of essential isoprenoids vs. nonessential isoprenoids especially under drought stress conditions or after jasmonic acid application. Phytochemistry 77, 149–161. doi: 10.1016/j.phytochem.2011.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Atkinson, N. J., and Urwin, P. E. (2012). The interaction of plant biotic and abiotic stresses: from genes to the field. J. Exp. Bot. 63, 3523–3543. doi: 10.1093/jxb/ers100

PubMed Abstract | CrossRef Full Text | Google Scholar

Babikova, Z., Gilbert, L., Bruce, T., Dewhirst, S. Y., Pickett, J. A., and Johnson, D. (2014). Arbuscular mycorrhizal fungi and aphids interact by changing host plant quality and volatile emission. Funct. Ecol. 28, 375–385. doi: 10.1111/1365-2435.12181

CrossRef Full Text | Google Scholar

Babikova, Z., Gilbert, L., Bruce, T. J. A., Birkett, M., Caulfield, J. C., Woodcock, C., et al. (2013). Underground signals carried through common mycelial networks warn neighbouring plants of aphid attack. Ecol. Lett. 16, 835–843. doi: 10.1111/ele.12115

PubMed Abstract | CrossRef Full Text | Google Scholar

Bai, Y., Kissoudis, C., Yan, Z., Visser, R. G. F., and van der Linden, G. (2018). Plant behaviour under combined stress: tomato responses to combined salinity and pathogen stress. Plant J. 93, 781–793. doi: 10.1111/tpj.13800

PubMed Abstract | CrossRef Full Text | Google Scholar

Balestrini, R., Chitarra, W., Antoniou, C., Ruocco, M., and Fotopoulos, V. (2018). Improvement of plant performance under water deficit with the employment of biological and chemical priming agents. J. Agric. Sci. doi: 10.1017/S002185961800012

CrossRef Full Text | Google Scholar

Balestrini, R., Gómez-Ariza, J., Lanfranco, L., and Bonfante, P. (2007). Laser microdissection reveals that transcripts for five plant and one fungal phosphate transporter genes are contemporaneously present in arbusculated cells. Mol. Plant Microbe Interact. 20, 1055–1062. doi: 10.1094/MPMI-20-9-1055

PubMed Abstract | CrossRef Full Text | Google Scholar

Balestrini, R., and Lumini, E. (2018). Focus on mycorrhizal symbioses. Appl. Soil Ecol. 123, 299–304. doi: 10.1016/j.apsoil.2017.09.001

CrossRef Full Text | Google Scholar

Battaglia, D., Bossi, S., Cascone, P., Digilio, M. C., Prieto, J. D., Fanti, P., et al. (2013). Tomato below ground–above ground interactions: trichoderma longibrachiatum affects the performance of Macrosiphum euphorbiae and its natural antagonists. Mol. Plant Microbe Interact. 26, 1249–1256. doi: 10.1094/MPMI-02-13-0059-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Berruti, A., Borriello, R., Della Beffa, M. T., Scariot, V., and Bianciotto, V. (2013). Application of nonspecific commercial AMF inocula results in poor mycorrhization in Camellia japonica L. Symbiosis 61, 63–76. doi: 10.1007/s13199-013-0258-7

CrossRef Full Text | Google Scholar

Berruti, A., Lumini, E., Balestrini, R., and Bianciotto, V. (2015). Arbuscular mycorrhizal fungi as natural biofertilizers: let’s benefit from past successes. Front. Microbol. 6:1559. doi: 10.3389/fmicb.2015.01559

CrossRef Full Text | Google Scholar

Bucher, M. (2007). Functional biology of plant phosphate uptake at root and mycorrhiza interfaces. Tansley review. New Phytol. 173, 11–26. doi: 10.1111/j.1469-8137.2006.01935.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Catola, S., Centritto, M., Cascone, P., Ranieri, A., Loreto, F., Calamai, L., et al. (2018). Effects of single or combined water deficit and aphid attack on tomato volatile organic compound (VOC) emission and plant-plant communication. Environ. Exp. Bot. 153, 54–62. doi: 10.1016/j.envexpbot.2018.05.001

CrossRef Full Text | Google Scholar

Chen, A., Chen, X., Wang, H., Liao, D., Gu, M., Qu, H., et al. (2014). Genome-wide investigation and expression analysis suggest diverse roles and genetic redundancy of Pht1 family genes in response to Pi deficiency in tomato. BMC Plant Biol. 14:61. doi: 10.1186/1471-2229-14-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Chitarra, W., Pagliarani, C., Maserti, B., Lumini, E., Siciliano, I., Cascone, P., et al. (2016). Insights on the impact of arbuscular mycorrhizal symbiosis on tomato tolerance to water stress. Plant Physiol. 171, 1009–1023. doi: 10.1104/pp.16.00307

PubMed Abstract | CrossRef Full Text | Google Scholar

Comas, L. H., Becker, S. R., Cruz, V. M., Byrne, P. F., and Dierig, D. A. (2013). Root traits contributing to plant productivity under drought. Front. Plant Sci. 4:442. doi: 10.3389/fpls.2013.00442

PubMed Abstract | CrossRef Full Text | Google Scholar

Copolovici, L., Kännaste, A., Pazouki, L., and Niinemets, Ü. (2012). Emissions of green leaf volatiles and terpenoids from Solanum lycopersicum are quantitatively related to the severity of cold and heat shock treatments. J. Plant Physiol. 169, 664–672. doi: 10.1016/j.jplph.2011.12.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Copolovici, L., Kännaste, A., Remmel, T., Vislap, V., and Niinemets, Ü. (2011). Volatile emissions from Alnus glutionosa induced by herbivory are quantitatively related to the extent of damage. J. Chem. Ecol. 37, 18–28. doi: 10.1007/s10886-010-9897-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Copolovici, L., and Niinemets, Ü. (2010). Flooding induced emissions of volatile signalling compounds in three tree species with differing waterlogging tolerance. Plant Cell Environ. 33, 1582–1594. doi: 10.1111/j.1365-3040.2010.02166.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, A. (2011). Drought under global warming: a review. Wires Clim. Chang. 2, 45–65. doi: 10.1002/wcc.81

CrossRef Full Text | Google Scholar

De Smet, I., Signora, L., Beeckman, T., Inzé, D., Foyer, C. H., and Zhang, H. (2003). An abscisic acid-sensitive checkpoint in lateral root development of Arabidopsis. Plant J. 33, 543–555. doi: 10.1046/j.1365-313X.2003.01652.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dell’Amico, J., Torrecillas, A., Rodríguez, P., Morte, A., and Sánchez-Blanco, M. J. (2002). Responses of tomato plants associated with the arbuscular mycorrhizal fungus Glomus clarum during drought and recovery. J. Agric. Sci. 138, 387–393. doi: 10.1017/S0021859602002101

CrossRef Full Text | Google Scholar

Digilio, M. C., Cascone, P., Iodice, L., and Guerrieri, E. (2012). Interactions between tomato volatile organic compounds and aphid behaviour. J. Plant Interact. 7, 322–325. doi: 10.1080/17429145.2012.727104

CrossRef Full Text | Google Scholar

Digilio, M. C., Corrado, G., Sasso, R., Coppola, V., Iodice, L., Pasquariello, M., et al. (2010). Molecular and chemical mechanisms involved in aphid resistance in cultivated tomato. New Phytol. 187, 1089–1101. doi: 10.1111/j.1469-8137.2010.03314.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, F., Fu, X., Watanabe, N., Su, X., and Yang, Z. (2016). Recent advances in the emission and functions of plant vegetative volatiles. Molecules 21:124. doi: 10.3390/molecules21020124

PubMed Abstract | CrossRef Full Text | Google Scholar

du Jardin, P. (2015). Plant biostimulants: definition, concept, main categories and regulation. Sci. Hortic. 196, 3–14. doi: 10.1016/j.scienta.2015.09.021

CrossRef Full Text | Google Scholar

Duan, J., Tian, H., Drijber, R. A., and Gao, Y. (2015). Systemic and local regulation of phosphate and nitrogen transporter genes by arbuscular mycorrhizal fungi in roots of winter wheat (Triticum aestivum L.). Plant Physiol. Biochem. 96, 199–208. doi: 10.1016/j.plaphy.2015.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Enstone, D. E., Peterson, C. A., and Ma, F. (2002). Root endodermis and exodermis: structure, function, and responses to the environment. J. Plant Growth Regul. 21, 335–351. doi: 10.1007/s00344-003-0002-2

CrossRef Full Text | Google Scholar

Fusconi, A. (2014). Regulation of root morphogenesis in arbuscular mycorrhizae: what role do fungal exudates, phosphate, sugars and hormones play in lateral root formation? Ann. Bot. 113, 19–33. doi: 10.1093/aob/mct258

PubMed Abstract | CrossRef Full Text | Google Scholar

Glassop, D., Smith, S. E., and Smith, F. W. (2005). Cereal phosphate transporters associated with the mycorrhizal pathway of phosphate uptake into roots. Planta 222, 688–698. doi: 10.1007/s00425-005-0015-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Grattan, S. R., and Grieve, C. M. (1999). Salinity-mineral nutrient relations in horticultural crops. Sci. Hortic. 78, 127–157. doi: 10.1016/S0304-4238(98)00192-7

CrossRef Full Text | Google Scholar

Gray, S. B., and Brady, S. M. (2016). Plant developmental responses to climate change. Dev. Biol. 419, 64–77. doi: 10.1016/j.ydbio.2016.07.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Guerrieri, E. (2016) “Who’s listening to talking plants?” in Deciphering Chemical Language of Plant Communication. Signaling and Communication in Plants, eds J. Blande and R. Glinwood (New York, NY: Springer Publishing Company).

Google Scholar

Guerrieri, E., Lingua, G., Digilio, M. C., Massa, N., and Berta, G. (2004). Do interactions between plant roots and the rhizosphere affect parasitoid behaviour? Ecol. Entomol. 29, 753–756. doi: 10.1111/j.0307-6946.2004.00644.x

CrossRef Full Text | Google Scholar

Guerrieri, E., Poppy, G. M., Powell, W., Rao, R., and Pennacchio, F. (2002). Plant-to-plant communication mediating in flight orientation of Aphidius ervi. J. Chem. Ecol. 9, 1703–1715. doi: 10.1023/A:1020553531658

PubMed Abstract | CrossRef Full Text | Google Scholar

Gutjahr, C., Casieri, L., and Paszkowski, U. (2009). Glomus intraradices induces changes in root system architecture of rice independently of common symbiosis signaling. New Phytol. 184, 829–837. doi: 10.1111/j.1469-8137.2009.02839.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hewitt, E. J. (1966). Sand and Water Culture Methods Used in the Study of Plant Nutrition. Farnham Royal: Commonwealth Agricultural Bureaux.

Google Scholar

Holopainen, J. K., and Gershenzon, J. (2010). Multiple stress factors and the emission of plant VOCs. Trends Plant Sci. 15, 176–184. doi: 10.1016/j.tplants.2010.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Javot, H., Penmetsa, R. V., Terzaghi, N., Cook, D. R., and Harrison, M. J. (2007). A medicago truncatula phosphate transporter indispensable for the arbuscular mycorrhizal symbiosis. Proc. Natl. Acad. Sci. U.S.A. 104, 1720–1725. doi: 10.1073/pnas.0608136104

PubMed Abstract | CrossRef Full Text | Google Scholar

Kask, K., Kännaste, A., Talts, E., Copolovici, L., and Niinemets, Ü. (2016). How specialized volatiles respond to chronic and short-term physiological and shock heat stress in Brassica nigra. Plant Cell Environ. 39, 2027–2042. doi: 10.1111/pce.12775

PubMed Abstract | CrossRef Full Text | Google Scholar

Lenoir, I., Fontaine, J., and Sahraoui, L.-H. (2016). Arbuscular mycorrhizal fungal responses to abiotic stresses: a review. Phytochemistry 123, 4–15. doi: 10.1016/j.phytochem.2016.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, F., Xu, Y., Han, G., Wang, W., Li, X., and Cheng, B. (2018). Identification and functional characterization of a maize phosphate transporter induced by mycorrhiza formation. Plant Cell Physiol. 59, 1683–1694. doi: 10.1093/pcp/pcy094

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, F., Xu, Y. J., Jiang, H. H., Jiang, C. S., Du, Y. B., Gong C., et al. (2016). Systematic identification, evolution and expression analysis of the Zea mays PHT1 gene family reveals several new members involved in root colonization by arbuscular mycorrhizal fungi. Int. J. Mol. Sci. 17:E930. doi: 10.3390/ijms17060930

PubMed Abstract | CrossRef Full Text | Google Scholar

Loreto, F., and Schnitzler, J. P. (2010). Abiotic stresses and induced BVOCs. Trends Plant Sci. 15, 154–166. doi: 10.1016/j.tplants.2009.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Nachappa, P., Culkin, C. T., Saya, P. M. II, Han, J., and Nalam, V. J. (2016). Water stress modulates soybean aphid performance, feeding behavior, and virus transmission in soybean. Front. Plant Sci. 7:552. doi: 10.3389/fpls.2016.00552

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagy, R., Drissner, D., Amrhein, N., Jakobsen, I., and Bucher, M. (2009). Mycorrhizal phosphate uptake pathway in tomato is phosphorus-repressible and transcriptionally regulated. New Phytol. 181, 950–959. doi: 10.1111/j.1469-8137.2008.02721.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagy, R., Karandashov, V., Chague, V., Kalinkevich, K., Tamasloukht, M., Xu, G., et al. (2005). The characterization of novel mycorrhiza-specific phosphate transporters from Lycopersicon esculentum and Solanum tuberosum uncovers functional redundancy in symbiotic phosphate transport in solanaceous species. Plant J. 42, 236–250. doi: 10.1111/j.1365-313X.2005.02364.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Niederbacher, B., Winkler, J., and Schnitzler, J.-P. (2015). Volatile organic compounds as non-invasive markers for plant phenotyping. J. Exp. Bot. 66, 5403–5416. doi: 10.1093/jxb/erv219

PubMed Abstract | CrossRef Full Text | Google Scholar

Niinemets, Ü. (2010). Mild versus severe stress and BVOCs: thresholds: priming and consequences. Trends Plant Sci. 15, 145–153. doi: 10.1016/j.tplants.2009.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’Hara, R. B., et al. (2013). Vegan: Community Ecology Package. R Package Version. 2.0.

Google Scholar

Olowu, R. A., Adewuyi, G. O., Onipede, O. J., Lawal, O. A., and Sunday, O. M. (2015). Concentration of heavy metals in root, stem and leaves of Acalypha indica and Panicum maximum jacq from three major dumpsites in Ibadan Metropolis, South West Nigeria. Am. J. Chem. 5, 40–48. doi: 10.5923/j.chemistry.20150501.06

CrossRef Full Text | Google Scholar

Osakabe, Y., Osakabe, K., Shinozaki, K., and Tran, L.-S. P. (2014). Response of plants to water stress. Front. Plant Sci. 5:86. doi: 10.3389/fpls.2014.00086

PubMed Abstract | CrossRef Full Text | Google Scholar

Pollastri, S., Savvides, A., Pesando, M., Lumini, E., Volpe, M. G., Ozudogru, E. A., et al. (2018). Impact of two arbuscular mycorrhizal fungi on Arundo donax L. response to salt stress. Planta 247, 573–585. doi: 10.1007/s00425-017-2808-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Quiroga, G., Erice, G., Aroca, R., Chaumont, F., and Ruiz-Lozano, J. M. (2017). Enhanced drought stress tolerance by the arbuscular mycorrhizal symbiosis in a drought-sensitive maize cultivar is related to a broader and differential regulation of host plant aquaporins than in a drought-tolerant cultivar. Front. Plant Sci. 8:1056. doi: 10.3389/fpls.2017.01056

PubMed Abstract | CrossRef Full Text | Google Scholar

R Development Core Team (2012). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available at: http://www.R-project.org/

Google Scholar

Rasmussen, R. (2001). “Quantification on the lightcycler,” in Rapid Cycle Real-time PCR, Methods and Applications, eds S. Meuer, C. Wittwer, and K. Nakagawara (Heidelberg: Springer Press).

Google Scholar

Rivero, J., Álvarez, D., Flors, V., Azcón-Aguilar, C., and Pozo, M. J. (2018). Root metabolic plasticity underlies functional diversity in mycorrhiza-enhanced stress tolerance in tomato. New Phytol. doi: 10.1111/nph.15295

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruiz-Lozano, J. M., Aroca, R., Zamarreño, A. M., Molina, S., Andreo-Jiménez, B., Porcel, R., et al. (2016). Arbuscular mycorrhizal symbiosis induces strigolactone biosynthesis under drought and improves drought tolerance in lettuce and tomato. Plant Cell Environ. 39, 441–452. doi: 10.1111/pce.12631

PubMed Abstract | CrossRef Full Text | Google Scholar

Sardans, J., and Peñuelas, J. (2004). Increasing drought decreases phosphorus availability in an evergreen Mediterranean forest. Plant Soil 267, 367–377. doi: 10.1007/s11104-005-0172-8

CrossRef Full Text | Google Scholar

Sasso, R., Iodice, L., Digilio, M. C., Carretta, A., Ariati, L., and Guerrieri, E. (2007). Host-locating response by the aphid parasitoid Aphidius ervi to tomato plant volatiles, J. Plant Interact. 2, 175–183. doi: 10.1080/17429140701591951

CrossRef Full Text | Google Scholar

Sasso, R., Iodice, L., Pickett, J. A., Woodcock, C. M., and Guerrieri, E. (2009). Behavioural and electrophysiological responses of Aphidius ervi (Hymenoptera: Braconidae) to tomato plant volatiles. Chemoecology 19, 195–201. doi: 10.1007/s00049-009-0023-9

CrossRef Full Text | Google Scholar

Sawers, R. J. H., Svane, S. F., Quan, C., Grønlund, M., Wozniak, B., Gebreselassie, M.-N., et al. (2017). Phosphorus acquisition efficiency in arbuscular mycorrhizal maize is correlated with the abundance of root-external hyphae and the accumulation of transcripts encoding PHT1 phosphate transporters. New Phytol. 214, 632–643. doi: 10.1111/nph.14403

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, F. A., and Smith, S. E. (2011). What is the significance of the arbuscular mycorrhizal colonisation of many economically important crop plants? Plant Soil 348, 63–79. doi: 10.1007/s11104-011-0865-0

CrossRef Full Text | Google Scholar

Smith, S. E., and Smith, F. A. (2012). Fresh perspectives on the roles of arbuscular mycorrhizal fungi in plant nutrition and growth. Mycologia 104, 1–13. doi: 10.3852/11-229

PubMed Abstract | CrossRef Full Text | Google Scholar

Sokal, R. R., and Rohlf, F. J. (1995). Biometry: The Principles and Practice of Statistics in Biological Research, 3rd Edn. New York, NY: W.H. Freeman and Co.

Google Scholar

Subramanian, K. S., Santhanakrishnan, P., and Balasubramanian, P. (2006). Responses of field grown tomato plants to arbuscular mycorrhizal fungal colonization under varying intensities of drought stress. Sci. Hortic. 107, 245–253. doi: 10.1016/j.scienta.2005.07.006

CrossRef Full Text | Google Scholar

Sun, T., Li, M., Shao, Y., Yu, L., and Ma, F. (2017). Comprehensive genomic identification and expression analysis of the phosphate transporter (PHT) gene family in Apple. Front. Plant Sci. 8:426. doi: 10.3389/fpls.2017.00426

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, Y., Kobae, Y., Mizuno, T., and Hata, S. (2012). Identification and expression analysis of arbuscular mycorrhiza-inducible phosphate transporter genes of soybean. Biosci. Biotech. Biochem. 76, 309–313. doi: 10.1271/bbb.110684

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, H., Yuan, X., Duan, J., Li, W., Zhai, B., and Gao, Y. (2017). Influence of nutrient signals and carbon allocation on the expression of phosphate and nitrogen transporter genes in winter wheat (Triticum aestivum L.) roots colonized by arbuscular mycorrhizal fungi. PLoS One 12:e0172154. doi: 10.1371/journal.pone.0172154

PubMed Abstract | CrossRef Full Text | Google Scholar

Volpe, V., Dell’Aglio, E., Giovannetti, M., Ruberti, C., Costa, A., Genre, A., et al. (2013). An AM-induced, MYB-family gene of Lotus japonicus (LjMAMI) affects root growth in an AM-independent manner. Plant J. 73, 442–455. doi: 10.1111/tpj.12045

PubMed Abstract | CrossRef Full Text | Google Scholar

Volpe, V., Giovannetti, M., Sun, X. G., Fiorilli, V., and Bonfante, P. (2016). The phosphate transporters LjPT4 and MtPT4 mediate early root responses to phosphate status in non mycorrhizal roots. Plant Cell Environ. 39, 660–671. doi: 10.1111/pce.12659

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D., Lv, S., Jiang, P., and Li, Y. (2017). Roles, regulation, and agricultural application of plant phosphate transporters. Front. Plant Sci. 8:817. doi: 10.3389/fpls.2017.00817

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, B., Yao, Z., Zhao, S., Guo, K., Sun, J., and Zhang, H. (2014). Arbuscular mycorrhizal fungal application to improve growth and tolerance of processing tomato (Lycopersicon esculentum Miller) under drought stress. J. Food Agric. Environ. 12, 452–457.

Google Scholar

Wu, Y., Spollen, W. G., Sharp, R. E., Hetherington, P. R., and Fry, S. C. (1994). Root growth maintenance at low water potentials: increased activity of xyloglucan endotransglycosylase and its possible regulation by abscisic acid. Plant Physiol. 106, 607–615. doi: 10.1104/pp.106.2.607

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, S.-Y., Grønlund, M., Jakobsen, I., Grotemeyer, M. S., Rentsch, D., Miyao, A., et al. (2012). Nonredundant regulation of rice arbuscular mycorrhizal symbiosis by two members of the phosphate transporter1 gene family. Plant Cell 24, 4236–4251. doi: 10.1105/tpc.112.104901

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Meng, S., Li, M., and Zhao, Z. (2016). Genomic identification and expression analysis of the phosphate transporter gene family in poplar. Front. Plant Sci. 7:1398. doi: 10.3389/fpls.2016.01398

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, Y.-N., Wang, P., Liu, C.-Y., Ni, Q.-D., Zhang, D.-J., and Wu, Q.-S. (2017). Mycorrhizal trifoliate orange has greater root adaptation of morphology and phytohormones in response to drought stress. Sci. Rep. 7:41134. doi: 10.1038/srep41134

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: aphid, arbuscular mycorrhizal symbiosis, phosphate transporter, plant tolerance, Solanum lycopersicum, volatile organic compound, water deficit

Citation: Volpe V, Chitarra W, Cascone P, Volpe MG, Bartolini P, Moneti G, Pieraccini G, Di Serio C, Maserti B, Guerrieri E and Balestrini R (2018) The Association With Two Different Arbuscular Mycorrhizal Fungi Differently Affects Water Stress Tolerance in Tomato. Front. Plant Sci. 9:1480. doi: 10.3389/fpls.2018.01480

Received: 05 July 2018; Accepted: 21 September 2018;
Published: 09 October 2018.

Edited by:

Youssef Rouphael, Università degli Studi di Napoli Federico II, Italy

Reviewed by:

Graziella Berta, Università degli Studi del Piemonte Orientale, Italy
Rosa Hermosa, Universidad de Salamanca, Spain

Copyright © 2018 Volpe, Chitarra, Cascone, Volpe, Bartolini, Moneti, Pieraccini, Di Serio, Maserti, Guerrieri and Balestrini. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Raffaella Balestrini, raffaella.balestrini@ipsp.cnr.it

These authors have contributed equally to this work