Abstract
Beneficial plant–microorganism interactions are widespread in nature. Among them, the symbiosis between plant roots and arbuscular mycorrhizal fungi (AMF) is of major importance, commonly improving host nutrition and tolerance against environmental and biotic challenges. Metabolic changes were observed in a well-established symbiosis between tomato and two common AMF: Rhizophagus irregularis and Funneliformis mosseae. Principal component analysis of metabolites, determined by non-targeted liquid chromatography–mass spectrometry, showed a strong metabolic rearrangement in mycorrhizal roots. There was generally a negative impact of mycorrhizal symbiosis on amino acid content, mainly on those involved in the biosynthesis of phenylpropanoids. On the other hand, many intermediaries in amino acid and sugar metabolism and the oxylipin pathway were among the compounds accumulating more in mycorrhizal roots. The metabolic reprogramming also affected other pathways in the secondary metabolism, mainly phenyl alcohols (lignins and lignans) and vitamins. The results showed that source metabolites of these pathways decreased in mycorrhizal roots, whilst the products derived from α-linolenic and amino acids presented higher concentrations in AMF-colonized roots. Mycorrhization therefore increased the flux into those pathways. Venn-diagram analysis showed that there are many induced signals shared by both mycorrhizal interactions, pointing to general mycorrhiza-associated changes in the tomato metabolome. Moreover, fungus-specific fingerprints were also found, suggesting that specific molecular alterations may underlie the reported functional diversity of the symbiosis. Since most positively regulated pathways were related to stress response mechanisms, their potential contribution to improved host stress tolerance is discussed.
Introduction
Beneficial organisms are common in the rhizosphere and they provide important ecosystem services (Philippot et al., 2013). They can greatly contribute to plant performance by improving nutrition, stress tolerance and plant phenotypic plasticity, an important advantage in heterogeneous environments where precise allocation of limited resources between growth and stress resistance is critical for survival (Pozo et al., 2015). Among these beneficial organisms, soil-borne fungi from the phylum Glomeromycota, known as arbuscular mycorrhizal fungi (AMF) deserve special attention. They are able to establish the most ancient and widespread plant–fungal symbiosis, known as arbuscular mycorrhizas (AM), with more than 80% of all terrestrial plant species (Smith and Read, 2008). AMF are obligate biotrophs, and it is assumed that the host plant allocates photosynthates to the fungus for the formation, maintenance, and function of mycorrhizal structures (). In return, the AMF improve plant acquisition of water and mineral nutrients. This symbiosis is extremely important for the uptake of inorganic phosphate, but also contributes to the uptake of nitrogen (N) and various trace elements (Smith and Smith, 2011; ). Besides plant nutrition, the symbiosis impacts the plant’s ability to overcome biotic and abiotic stresses, commonly improving host tolerance to unfavorable environmental conditions and resistance to pathogens (; ; Ruíz-Lozano et al., 2012; Selosse et al., 2014). The establishment and maintenance of the association requires a high degree of coordination between both partners, and bidirectional (plant and fungal) control assures a fair trade of resources between the symbionts (). Indeed, a precise regulation of host hormone levels has been proposed as a central mechanism in the regulation of the interaction (Pozo et al., 2015).
Although there is no strict partner specificity in AM, the outcome of AM interactions depends on the interacting partners and the environmental conditions (Walder et al., 2012; Smith and Smith, 2015). Actually, there is evidence for “functional diversity” occurring as plant and fungal genotypes determine the benefits of the interaction; some combinations being more efficient than others in terms of nutrition and/or stress resistance (; ; Mensah et al., 2015; Smith and Smith, 2015). For example, colonization of roots by Funneliformis mosseae or Rhizophagus irregularis, the two widespread AMF used in this study, resulted in different levels of bioprotection against Phytophthora root rot in tomato or Fusarium wilt in melon, and to drought stress in lettuce (Pozo et al., 2002; Marulanda et al., 2003; Porcel et al., 2006; Martínez-Medina et al., 2009).
Different studies have tried to unravel the molecular mechanisms regulating AM and their impact on plant fitness, most of them focused on differential gene expression and protein profiles (Liu et al., 2007; ; López-Ráez et al., 2010; ). A consistent output of “omic” studies on AM has been the observation that mycorrhizal colonization significantly impacts host gene expression and metabolomic profiles (Salvioli and Bonfante, 2013). Metabolomics is a valuable technology which provides comprehensive quantitative profiling of metabolites in biological systems; and liquid or gas chromatography coupled with mass spectrometry (LC–MS or GC–MS) are widely used analytical tools for such untargeted metabolomic studies (). The first metabolomic studies in AM interactions were targeted analyses focusing on a few well-characterized metabolites to monitor plant responses (Stumpe et al., 2005; Sawada et al., 2009; López-Ráez et al., 2010), however, a whole picture of the impact of AM on the general metabolic profile was missing. Nowadays, metabolomic untargeted approaches allow the separation and detection of a wide range of metabolites, such as amino acids, fatty acids, organic acids, sugar phosphates, nucleotides, and glycoside derivatives, providing a global fingerprint about the quantitative and qualitative changes in the secondary metabolism of the host plant (). The LC–MS technique is highly sensitive for the detection of key molecules in the phenotypic mechanisms underlying organism responses to abiotic or biotic interactions (Sardans et al., 2011; ). Indeed, untargeted metabolomic studies have revealed the importance of metabolic reprogramming as a determinant in other plant–microbe symbioses, e.g., legume–rhizobia associations (Zhang et al., 2012) and ectomycorrhizas in poplar (Tschaplinski et al., 2014).
Few studies have addressed the metabolome reprogramming associated with AM and knowledge concerning metabolomic transitions in mycorrhizal plants remains restricted to a few recent studies. Concerning changes in mycorrhizal roots, research has been limited to legumes, in particular to the model plant Medicago truncatula and to one AMF strain, R. irregularis (Schliemann et al., 2008; ). The studies have revealed an impact of R. irregularis colonization on primary and secondary metabolism, mainly on amino acids (glutamic acid, aspartic acid, and asparagine), fatty acids (palmitic and oleic acids), apocarotenoids (cyclohexenone and mycorradicin derivatives), and isoflavonoids (daidzein, ononin, and malonylononin). Remarkably, the impact of the symbiosis on the host metabolism extends to aboveground tissues and can vary with environmental conditions (). Moreover, recent multispecies metabolomic analysis of leaves from R. irregularis-colonized plants showed a common core of mycorrhiza-related and highly accumulated metabolites shared by dicotyledonous and monocotyledonous plants, although they also evidenced many species-specific responses (Schweiger et al., 2014).
This study aims to decipher the impact of a well-established mycorrhizal association with two widespread and well characterized AMF (F. mosseae and R. irregularis) in the root metabolic profile of a non-legume, tomato (Solanum lycopersicum), where the benefits of AM have been shown to be agronomically relevant in terms of stress resistance and fruit quality (Pozo et al., 2002; ; ; ; ; ; Zouari et al., 2014). LC–MS revealed important changes in the metabolome of mycorrhizal tomato roots and we discuss the potential relevance of these changes in host fitness.
Materials and Methods
Plant Material and AMF Inoculation
Arbuscular mycorrhizal fungi isolates of R. irregularis (BEG 121; formerly Glomus intraradices) and F. mosseae (BEG12; formerly Glomus mosseae) from the International Bank of Glomeromycota1 are continuously maintained in an open-pot culture of Trifolium repens L. mixed with Sorghum vulgare Pers. (Steud.) Millsp. & Chase plants in a greenhouse. The inocula consist of substrate (vermiculite/sepiolite, 1:1), spores, mycelia, and infected root fragments from those cultures. Tomato seeds (Solanum lycopersicum L. cv. Moneymaker) were surface disinfected by immersion in 4% NaHClO (10 min) containing 0.02% (v/v) Tween 20®, rinsed thoroughly with sterile water and incubated for 3 days in an open container with sterile vermiculite at 25°C in darkness. Plantlets were transferred to 1 L pots containing a sterile sand:soil (1:1) mixture. Pots for mycorrhizal treatments were inoculated by adding 10% (v/v) F. mosseae or R. intraradices inoculum. Uninoculated control plants received the same amount of autoclaved mycorrhizal inoculum together with a 3 ml aliquot of a filtrate (<20 μm) of both AM inocula, in order to provide the general microbial population but free of AMF propagules.
A total of ten plants were used for each treatment. Plants were randomly distributed and grown in a greenhouse at 24/16°C with a 16/8 h diurnal photoperiod and 70% humidity. Plants were watered three times a week with Long Ashton nutrient solution () containing 25% of the standard phosphorus (P) concentration, and water was supplied daily to maintain the substrate at 100% field capacity, as reported in . Plants were harvested after 8 weeks, the fresh weight of shoots and roots was determined, and the material immediately frozen in liquid N and stored at -80°C. An aliquot of each individual root system was reserved for mycorrhizal quantification.
Determination of Mycorrhizal Colonization
Mycorrhizal colonization was estimated after clearing washed roots in 10% KOH and subsequent staining of fungal structures with 5% ink in 2% acetic acid (Vierheilig et al., 2005). The extent of mycorrhizal colonization (expressed as percentage of total root length colonized by the AMF) was calculated according to the gridline intersection method () using a Nikon Eclipse 50i microscope and bright field conditions.
Phosphorus, Carbon, and Nitrogen Content
Total P, carbon (C), and N content in the roots was measured at the Ionomic Laboratory of Technical Services of the Centro de Edafología y Biología Agraria del Segura (CSIC), Murcia, Spain. Three biological replicates, each consisting of a pool of roots from three independent plants (nine plants in total), were analyzed for each treatment. Frozen roots were ground to a fine powder and lyophilized. P concentrations were analyzed after an acid digestion of the samples, by inductively coupled plasma optical emission spectrometry (ICP-OES; ICAP 6500 DUO THERMO). Total C and N contents were determined using an Elemental Analyzer (LECO TRUSPEC CN) according to standard procedures.
Analysis of Gene Expression by RT-qPCR
Total RNA from tomato roots was extracted and treated with DNase using the Direct-zol RNA MiniPrep kit (Zymo Research). Subsequently, the RNA was purified through a column using the RNA Clean & Concentrator-5 kit (Zymo Research), and stored at -80°C until use. The first-strand cDNA was synthesized with 1 μg of purified total RNA using the iScript cDNA Synthesis kit (Bio-Rad). Four independent biological replicates were analyzed per treatment. All kits were used according to the manufacturer’s suggested protocols.
Relative quantification of specific mRNA levels was performed using the comparative 2-Δ(ΔCt) method (Livak and Schmittgen, 2001). Expression values were normalized using the housekeeping gene SlEF-1α (López-Ráez et al., 2010), which encodes for the tomato elongation factor-1α. The functionality of AM symbiosis was quantified using the marker gene LePT4, which encodes a mycorrhiza-inducible phosphate transporter expressed in arbusculated cells (). Nucleotide sequences of the primers used were: SlEF-1α-F 5′-GATTGGTGGTATTGGAACTGTC-3′, SlEF-1α-R 5′-AGCTTCGTGGTGCATCTC-3′; LePT4-F 5′-GAAGGGGAGCCATTTAATGTGG-3′, LePT4-R 5′-ATCGCGGCTTGTTTAGCATTTC-3′.
Reagents and Standards
All standards, including amino acids, salicylic acid, phenols, IAA, 5-Hydroxyindole-3-acetic acid, Indole-3-acetamide, N-(3-indoley lacetyl)-L-alanine, Indole-3-carboxaldehyde, Methyl indole-acetate, jasmonic acid (JA), abscisic acid, salicylic acid glucoside ester, OPDA, carboxylic acids, and sugars were purchased from SIGMA (Barcelona, Spain). Methanol (HPLC grade) was obtained in SIGMA (Barcelona, Spain), formic acid and NaOH were obtained from J.T Baker (Deventer, Holland). Indole-3-carboxylic acid and 1,4-diaminobutane were obtained from VWR (Barcelona, Spain).
Liquid Chromatography and ESI Mass Spectrometry
LC–ESI Full Scan Mass Spectrometry (Q-TOF Instrument)
Freeze-dried roots (50 mg) were homogenized on ice in 1 ml of MeOH:H2O (10:90) containing 0.01% of HCOOH. The homogenate was centrifuged at 15000 rpm for 15 min at 4°C. The supernatant was recovered and filtered through 0.2 μm cellulose filters (Regenerated Cellulose Filter, 0.20 μm, 13 mm D. pk/100; Teknokroma). A 20 μl aliquot of was injected in the HPLC. The full metabolomic profiling was performed using an Acquity UPLC system (Waters, Mildford, MA, USA) interfaced to hybrid quadrupole time-of-flight (QTOF MS Premier). Analytes were eluted with an aqueous methanol gradient containing 0.01% HCOOH. Three biological replicates, each consisting of a pool of three independent plants (nine plants in total), were randomly injected in duplicate for every treatment. The LC separation was performed with an HPLC SunFire C18 analytical column, 5 μm particle size, 2.1 mm × 100 mm (Waters). Solvent gradients and further chromatographic conditions were performed as previously described (; ). The LC–ESI Q-TOF MS library of plant compounds was used for a straight identification in full-scan analysis. Standards for phenols, indolic compounds, amino acids, hormones and their derivatives (up to 93 compounds) were prepared (100 ppb) in a composite solution (Supplementary Table S2). The standard solution was injected through the HPLC in both positive and negative electro-spray ionization (ESI+; ESI-) to identify compounds by matching exact mass and retention time between standard and experimental samples.
Full Scan Data Analysis
Data were acquired in centroid mode and subsequently transformed into cdf files using the Databridge from MassLynx 4.1 software (MassLynx 4.1, Waters). Chromatographic signals were processed using the software R for statistical purposes2. Signals from ESI+ and ESI- were processed separately. Peak peaking, grouping and signal corrections were performed using the XCMS algorithm (Smith et al., 2006). Metabolite amounts were analyzed on the basis of normalized peak area units relative to the dry weight. The Kruskal–Wallis test (p < 0.05) was applied to analyze the metabolomic differences between treatments. To determine a global behavior of the signals, principal component analyses (PCA) plots were generated using the Multibase 2015 algorithm3. Statistical and heat map analysis were performed using the MarVis Suit 2.0 software tool for clustering and visualization of metabolic biomarkers (). Adduct, isotope correction, clustering, and color heat map visualization were also performed by using associated software packages MarVis Filter and MarVis Cluster.
Statistical Analyses
All statistical analyses (ANOVA, post hoc, and t-test) were conducted using Statgraphics Plus 3.1 (Rockville, MD, USA), “R” software version 2.9.2 (R Development Core Team)2 and the XCMS package.
Results
Root Colonization by F. mosseae and R. irregularis and Physiological Status of the Plant
Plants were harvested 8 weeks after inoculation with the mycorrhizal fungi. Root and shoot fresh weights were determined and mycorrhizal colonization, P, C, and N content in the roots were analyzed (Figure 1 and Supplementary Table S1). Staining of fungal structures within the roots showed that the mycorrhizal symbiosis was well established in both inoculation treatments, with abundant fungal colonization of the root cortex and well-formed arbuscules. Vesicles, the fungal storage structures, were more abundant in the roots with the most effective colonizer (R. irregularis), as has been described in previous studies (López-Ráez et al., 2010; Figures 1A,B). Absence of fungal structures was confirmed in roots of the non-mycorrhizal controls, and the extent of root length colonized by F. mosseae or R. irregularis differed significantly (Figure 1C, p < 0.01). The functionality of the symbiosis was assessed by analyzing the expression of the tomato gene LePT4, which encodes a phosphate transporter induced in arbuscule containing cells, where most of the nutrients exchange takes place; it is therefore used as marker of a functional symbiosis. A very strong induction of LePT4 expression was detected in mycorrhizal roots, reaching similar levels in the interaction with both fungi (Figure 1D).
FIGURE 1
The symbiosis did not have a significant effect on shoot or root biomass under our experimental conditions (Supplementary Table S1). However, both AM treatments enhanced the P content in roots compared to those from non-mycorrhizal controls (Figure 1E). Similarly, N content was higher in mycorrhizal roots (Supplementary Table S1), while total carbon content in the roots remained unaltered. Accordingly the C/N ratio showed a significant reduction in mycorrhizal roots (Figure 1F).
AMF Colonization has a Strong Impact on the Metabolic Profile of the Host Roots
We analyzed the reprogramming of the tomato root metabolism associated with well-established symbiosis with each AMF. Following the chromatographic analysis, a bioinformatic processing of the detected signals was performed, and cluster and functional pathway analyses were performed in order to obtain plausible biological information of such metabolic reprogramming.
Untargeted metabolomic analysis of root extracts via HPLC coupled with a quadrupole time-of-flight mass spectrometer revealed a total of 1407 signals in ESI- mode and 1860 signals in ESI+ mode. A supervised principal component analysis of these signals (p < 0.1, 847 and 1029 signals in ESI- and ESI+ mode, respectively) showed a clearly separated behavior between roots colonized by F. mosseae (Fm) or by R. irregularis (Ri) and non-mycorrhizal roots (Nm), (Figure 2A). According to the two main components, no overlap was observed between the mycorrhizal and Nm groups in any of the ESI modes. It is noteworthy that ESI- showed a similar behavior between signals detected for Ri and Fm. Hierarchical cluster analysis of the different groups confirmed previous observations: signals from roots colonized by both AMF clustered closely compared to those from Nm roots (Supplementary Figure S1A).
FIGURE 2
The clusters corresponding to compounds with the most contrasting accumulation patterns across different treatments were selected from the heat map analysis (Figure 2B) for detailed analysis. The clusters of selected signals were analyzed separately (Supplementary Table S2). Firstly, the number of over-accumulated compounds specific to Nm, Fm, and Ri was subtracted from the heat map. Secondly, those signals that were highly accumulated in two of the experimental conditions (such as Nm + Fm, Nm + Ri, or Fm + Ri) were also isolated for subsequent Venn-diagram and pathway analysis (Supplementary Figure S1B and Table S2). The selection, including 1876 signals, contained 300 differentially accumulated in mycorrhizal roots, as illustrated in the Venn diagram (Supplementary Figure S1B). These signals are of interest as they may correspond to compounds relevant for the known benefits of the mycorrhizal interaction, including improved host stress resistance. These metabolites can be generated as a plant response to the AMF colonization or by the AMFs themselves. Despite the core of compounds highly accumulated in both AM roots, there are many specific signals only triggered either by Fm (85 signals) or Ri (35 signals), (Supplementary Figure S1B). Interestingly, the metabolic impact of Fm is stronger than that of Ri. These results suggest that tomato plants retain common responses to different AMF, as they are conserved in both interactions, but there is, in addition, a set of responses that may be specific to particular interactions.
In order to understand the biological meaning of this metabolic transition, we classified the signals contained in the selected clusters from the heat map analysis and Venn diagram (Supplementary Figure S1B) into a pathway ontology using the MarVis Pathway 2.0 () linked to the KEGG Solanum lycopersicum database (Supplementary Table S2). We particularly focused on those signals that were strongly reduced in both mycorrhizal treatments (Cluster 1), signals highly accumulated in both mycorrhizal interactions (Cluster 2), signals exclusively accumulated in Fm colonized roots (Cluster 3), and finally, signals exclusively accumulated in Ri colonized roots (Cluster 4). Clearly, the major impact of AM on plant metabolism takes place in the primary metabolism, mainly in the amino acid and sugar metabolism (including many hits among the tricarboxylic and other carboxylic acids) but also in some specific secondary metabolites, such as phenolic alcohol derivatives, vitamins, and plant hormones, particularly oxylipins and cytokinins (Supplementary Table S2). Both fungi impacted the 13-LOX oxylipin pathway, with multiple hits in the linoleic and α-linolenic acid metabolism. Remarkably, a clearly overrepresented category in mycorrhiza-enriched compounds is that of metabolites related to ATP-binding cassette (ABC) transporters. ABC transporters are largely expressed in roots and mediate the transport of many secondary metabolites with signaling and defensive functions (Yazaki, 2006).
A closer look at the pathways containing signals with lower levels in mycorrhizal roots (Cluster 1) revealed that other compounds from the same pathways are strongly over-represented in AM. This suggests that the reduction of these compounds in mycorrhizal roots is a consequence of the metabolic flux along the pathways that reduces substrates of a given reaction accumulating the product compounds. Regarding the specific signals, Fm had a stronger impact on amino acids, sugars, and phenolics compared with Ri.
Impact of the Arbuscular Mycorrhizal Symbiosis on Amino Acid Metabolism
One of the major pathways altered in mycorrhizal roots was the metabolism of the amino acids (Figure 3 and Supplementary Table S2). Phe, Tyr, Tryp, and Leu/Ile were consistently less concentrated in mycorrhizal roots colonized by both AMF. Contrastingly, many metabolites derived from Cys, Lys, Ala, Gln, Phe, Tyr, and Trp were highly accumulated in mycorrhizal roots (Supplementary Table S2). The higher concentration of amino acid derived metabolites in AM would explain the lower concentration of free amino acids as metabolic sources (Figure 3). It is noteworthy that Phe and Tyr are the main amino acids that generate phenolic acids and their derivatives, highly accumulated in AM roots. This observation suggests a very likely circulation of the basic amino acids into more complex secondary metabolites, which are indeed found in higher concentrations in the symbiotic roots. On the other hand, glutamate (Glu) and aspartate (Asp) were found in higher levels in both mycorrhizal roots. As their active role in the incorporation of N in AM plants is well reported (; Schliemann et al., 2008), their higher levels are consistent with the increase in N observed in both mycorrhizal treatments (Figure 1B). Finally, some AMF-dependent regulation of the amino acids was also observed. For example, His and Met accumulated at higher levels in Fm, whilst they were hardly present in Nm and Ri.
FIGURE 3
Impact of the Arbuscular Mycorrhizal Symbiosis on Amino Acid Derived Compounds: Phenolic Alcohol Derivatives, Benzylisoquinolines, and Conjugated Polyamines
The phenyl-alcohol metabolism was also strongly affected in AM. This metabolism includes deamination of Phe and Tyr by the Phenylalanine ammonia lyase enzyme (), after which a set of phenolic acids are converted into aldehydes and alcohols by successive reductions. The phenolic alcohols are precursors of important cell wall components such as monolignans and lignins. In addition, coumaryl and coniferyl alcohols can be converted into more complex flavonoids with cell protective functions (Wang et al., 2013). As described above, the identified upstream compounds of this pathway, Phe and Tyr, were found in lower levels in mycorrhizal roots (Figures 3 and 4), while the content of other intermediary compounds, such as ferulic acid, cumaryl alcohol, and coniferyl alcohol, was higher in the colonized roots (Figure 4). Additionally, other tentatively identified monolignans (400.152, 362.173, 354.110, and 398.137 m/z) were also more concentrated in mycorrhizal roots (Figure 4). These observations suggest that AMF stimulates a reorganization of specific cell wall components.
FIGURE 4
Other amino acid derived compounds related to defense were also found in higher levels in mycorrhizal roots. Plants often produce alkaloids to defend themselves against pests, diseases, and other external biological stimuli. Several mass signals identified as benzylisoquinoline alkaloids (BIAs), such as 271.107; 369.126; 332.112; and 273.124 m/z, were found in high quantities in AMF colonized roots (Supplementary Figure S2). Higher accumulation of several polyamines and their conjugates were also found, some putatively identified as spermidine (145.112 m/z), tricaffeoyl spermidine (631.274 m/z), and triferuloyl spermidine (673.242 m/z), (Supplementary Figure S2).
Impact of Arbuscular Mycorrhizal Symbiosis on the Oxylipin Pathway
The untargeted metabolomic analysis revealed α-linolenic acid derivatives as major metabolic targets for mycorrhizal symbiosis. Among the 45 signals related to the oxylipin pathways altered in mycorrhizal roots (Supplementary Table S2), 11 compounds were fully identified by either exact mass or fragmentation spectrum (Figure 5), and they corresponded to the 13-LOX branch of the oxylipin pathway. This branch leads to the biosynthesis of the phytohormone JA and derivatives, known to be altered in AM in different plant species (Wasternack and Hause, 2013; ). With the exception of α-linolenic acid, the source metabolite, the compounds identified in this pathway showed higher levels in mycorrhizal roots (Figure 5). Most of them showed higher concentrations in roots colonized by both AMF, although to different levels depending on the particular compound and the colonizing fungi. Remarkably, the levels of the bioactive forms of JA methyl-JA (Me-JA) and JA-Ile conjugates were accumulated in significantly higher levels only in the Fm roots.
FIGURE 5
Discussion
Several studies have detailed the transcriptional reprogramming in the host plant during interaction with AMF, not only transiently during the establishment of the symbiosis, but also in the maintenance of the symbiosis (; Liu et al., 2007; ). While there is evidence for these mycorrhiza associated transcriptional changes in multiple plant families, the metabolic impact on host roots have been only monitored in a few plant species, mainly legumes (Schliemann et al., 2008; ). Our research presents a complete metabolomic analysis in roots of a relevant crop, tomato, in symbiotic interaction with two different AMF: F. mosseae and R. irregularis, both known to increase tomato resistance against biotic and abiotic stresses and yield (Pozo et al., 2002; ; ; ; ; Vos et al., 2012). These fungi are among the most studied and widely distributed AMF in agricultural and ecological settings. R. irregularis is the most commonly used AMF in commercial inoculants, and is widely used as a model organism in AM research, because it is readily grown using in vitro cultivation in monoxenic conditions, and its genome is now available (; Tisserant et al., 2013). In contrast, F. mosseae cannot be cultivated in monoxenic cultures but is usually very efficient in increasing host resistance to pests and pathogens ().
Increases in P levels in the host following root colonization by AMF is one of the major and most reported benefits of mycorrhizal interactions, although the increase depends on the partners involved and the experimental conditions (Smith and Smith, 2015). However, AM are not always associated to increased vegetative biomass (Smith and Smith, 2011), and improved stress tolerance has been proposed as another major benefit of the symbiosis (; Selosse et al., 2014). In our experimental system, both F. mosseae and R. irregularis colonization increased the total P levels in tomato roots (cv. Moneymaker), whilst shoot and root biomass, root/shoot ratio and total carbon content were not significantly altered. Interestingly, total N was significantly higher in mycorrhizal tomato roots, and consequently the C/N ratio was reduced.
Our untargeted metabolome analysis confirmed, following restrictive statistical analyses, that the metabolome of mycorrhizal tomato roots is significantly different from that of non-mycorrhizal tomatoes. Pioneering work by Schliemann et al. (2008), and later showed that there are clear differences in the development and symbiosis-dependent primary and secondary metabolism of M. truncatula roots colonized by the AMF R. irregularis.
A large number of signals related to sugar and carboxylic acid metabolism showed elevated levels in mycorrhizal roots. This suggests that primary sugar metabolism was activated by the symbiosis as has been shown in multiple mycorrhizal systems (; Zouari et al., 2014) probably related to an increase in the host photosynthesis to increase C-fixation (). Total C content, however, remained unaltered in mycorrhizal roots, probably due to the fact that part of the host-derived C is taken by the fungal partner to maintain the mycelial network (). Indeed, transcriptional regulation of carbohydrate related genes and activation of sugar transporters are reported in AM (; ). Another notable target of mycorrhizal metabolic reprogramming corresponds to amino acid metabolism, which is one of the pathways with more hits amongst the AM-related differential signals (Supplementary Table S2). A strong reduction in the accumulation of several amino acids (Trp, Tyr, Phe, Ala, Leu) was observed, probably because of their function as source compounds of amino acid-derived secondary metabolites (Figure 3). For example, Phe and Tyr are precursors of the phenylpropanoid pathway, and several intermediaries of this pathway were found to be highly accumulated in the mycorrhizal roots. Both amino acids were found in low levels in mycorrhizal tomato roots, and the data suggest that most phenolic derivatives may have been redirected to the formation of lignans and lignins. An increase in lignans has been previously reported in other root-beneficial fungus interaction involving Piriformospora indica (). Regarding the lignins, it has been shown that mycorrhizal colonization can increase the lignin content of the root cell walls (Ziedan et al., 2011), and cell-wall lignification is one of the proposed mechanisms restricting penetration by phytopathogenic fungi in mycorrhizal roots (). Remarkably, a similar reduction of amino acid content was observed in aboveground tissues of mycorrhizal L. japonica (), and Arabidopsis plants treated with the defense-priming agent β-amino butyric acid showed a lower content of all amino acids except Glu (Pastor et al., 2014). Thus, it is tempting to speculate that the reduction of basic and aromatic amino acids is a common response to defense-priming stimuli. However, this putative relationship requires further experimental confirmation.
Despite an overall reduction in most amino acids, our study revealed a higher accumulation in the mycorrhizal roots of Glu and Asp (Figure 3), important amino acids in uptake of N by AMF extra-radical mycelium (). The elevated N levels observed in AM, together with the elevated levels of these amino acids, also reported in other mycorrhizal systems (Schliemann et al., 2008) suggest that N uptake and assimilation is stimulated in mycorrhizal tomato roots. An impact of AMF colonization on enzymes catalyzing the biosynthesis of N rich compounds such as alkaloids has been also described (Zeng et al., 2013). Our metabolomic analysis showed that both BIAs and conjugated polyamines, all with defense-related functions in plants, are also over accumulated in mycorrhizal compared to non-mycorrhizal roots.
One of the clear targets of root reprogramming in mycorrhizal roots was the oxylipin pathway. All the metabolites of the 13-LOX branch identified either through the fragmentation spectrum, exact mass or using standards, were more concentrated in mycorrhizal roots, including several bioactive forms of the phytohormone JA. This may explain the reduced levels of the source metabolite α-linolenic acid in AM. The induction of most of the intermediates of the pathway is consistently reproduced in both mycorrhizal root systems, although the relative levels of the different compounds differ according to the AMF. The alteration of multiple metabolites in the pathway, known to be precisely regulated by inter conversion among them (Wasternack and Hause, 2013) and differential accumulation in response to the particular AMF, support their involvement in fine-tuning of metabolic reprogramming in response to particular growth conditions, the partners involved and the symbiotic stage reached (Wasternack and Hause, 2013; ). Elevated levels of the JA-related volatile compounds Me-JA and CIS-jasmone, both with known roles in defense against biotic stresses, are reported here for the first time in mycorrhizal roots. Moreover, Me-JA is also known to be involved in the plant response to abiotic stresses like drought or salinity (). The exact role of JA and its derivatives in the control of AM remains controversial, as exogenous application of the hormone provides contrasting results, and JA deficient mutants have relatively subtle mycorrhizal phenotypes that seem to depend largely on the host plant species (Wasternack and Hause, 2013). In some plant systems LOX-silencing does not significantly affect AMF colonization, therefore it has been proposed that activation of JA-signaling is a downstream event triggered by this symbiosis (Wasternack and Hause, 2013). In fact, PvLOX2-silencing in common bean roots is reported to have no effect on mycorrhiza establishment, whilst it does impair mycorrhiza-induced resistance (Mora-Romero et al., 2014). Thus, the results suggest that AM functioning implies a precise regulation of the oxylipin pathway that may contribute to improving stress resistance in mycorrhizal plants. It is noteworthy that the genes coding for alkaloid biosynthetic enzymes are JA-inducible (Mishra et al., 2013; Wasternack and Hause, 2013). Additionally, phenylpropanoid-polyamine conjugated (PPCs), other N rich compounds related to defense responses, were found in elevated levels in the mycorrhizal roots and are also described to be under JA regulation in several plants, including tomato ().
Although the metabolic pathways altered by the mycorrhizal symbiosis were common to both F. mosseae and R. irregularis interactions, many compounds showed specific responses to one of the interactions. These results suggest, therefore, a fine-tuned developmental regulation of these pathways in an AMF-dependent manner. In terms of fitness costs, we have observed no reduction in plant biomass despite the extensive metabolic changes encountered in AM roots. AMF-induced resistance against soil pathogens in tomato roots seems to be related to cost-efficient defense regulation mechanisms (Pozo et al., 2002; Steinkellner et al., 2012; Vos et al., 2012; Selosse et al., 2014). There is experimental evidence of the role of the JA-signaling pathway in priming of plant defenses by mycorrhizas (Song et al., 2013). The contribution of other particular metabolic pathways altered in mycorrhizal roots to the enhanced host stress resistance–tolerance remains to be experimentally demonstrated.
To sum up, following a non-targeted full metabolomic approach, the metabolic transition in roots colonized by two widely distributed AMF, F. mosseae and R. irregularis, was successfully characterized. An important reprogramming of some major metabolic pathways in both mycorrhizal interactions was observed, pointing to common responses associated to AM, although there were also some AMF-specific responses. Many of the changes are related to plant defense mechanisms, and may underlay the well-known effects of the interaction on plant-stress tolerance. The identification of differentially regulated pathways in this study is instrumental to functional studies aiming to reveal the mechanistic basis of AM benefits, such as the improved resistance–tolerance to biotic and abiotic stresses. Moreover, these studies can pave the way to improving the biotechnological applications of AMF in agricultural settings and in the production of plant secondary metabolites with medicinal or nutritional properties (; Pedone-Bonfim et al., 2015).
Statements
Acknowledgments
We thank the SCIC of the Universitat Jaume I for technical support. This study was financially supported by grants AGL2012-39923 from the Spanish MINECO and P10-AGR-6516 from Junta de Andalucía. J. Rivero is supported by Ph.D. fellowship BES-2013-062638. We acknowledge J. A. López-Ráez for helpful discussions.
Conflict of interest
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: http://journal.frontiersin.org/article/10.3389/fmicb.2015.00598
References
1
AbdallahC.ValotB.GuillierC.MounierA.BalliauT.ZivyM.et al (2014). The membrane proteome of Medicago truncatula roots displays qualitative and quantitative changes in response to arbuscular mycorrhizal symbiosis.J. Proteomics108354–368. 10.1016/j.jprot.2014.05.028
2
AgutB.GamirJ.JacasJ. A.HurtadoM.FlorsV. (2014). Different metabolic and genetic responses in citrus may explain relative susceptibility to Tetranychus urticae.Pest Manag. Sci.701728–1741. 10.1002/ps.3718
3
ArocaR.Del Mar AlguacilM.VernieriP.Ruiz-LozanoJ. 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.56704–719. 10.1007/s00248-008-9390-y
4
BagoB.PfefferP. E.Shachar-HillY. (2000). Carbon metabolism and transport in arbuscular mycorrhizas.Plant Physiol.124949–958. 10.1104/pp.124.3.949
5
BaldiA.FarkyaS.JainA.GuptaN.MehraR.DattaV.et al (2010). Enhanced production of podophyllotoxins by co-culture of transformed Linum album cells with plant growth-promoting fungi.Pure Appl. Chem.82227–241. 10.1351/PAC-CON-09-02-09
6
BalestriniR.Gómez-ArizaJ.LanfrancoL.BonfanteP. (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.201055–1062. 10.1094/mpmi-20-9-1055
7
BarzanaG.ArocaR.PazJ. A.ChaumontF.Martínez-BallestaM. C.CarvajalM.et al (2012). Arbuscular mycorrhizal symbiosis increases relative apoplastic water flow in roots of the host plant under both well-watered and drought stress conditions.Ann. Bot.1091009–1017. 10.1093/aob/mcs007
8
CochraneF. C.DavinL. B.LewisN. G. (2004). The Arabidopsis phenylalanine ammonia lyase gene family: kinetic characterization of the four PAL isoforms.Phytochemistry651557–1564. 10.1016/j.phytochem.2004.05.006
9
DeclerckS.StrulluD. G.FortinJ. A. (2005). In Vitro Culture of Mycorrhizas.New York, NY: Springer-Verlag, Berlin Heidelberg. 10.1007/b138925
10
De VosR. C.MocoS.LommenA.KeurentjesJ. J.BinoR. J.HallR. D. (2007). Untargeted large-scale plant metabolomics using liquid chromatography coupled to mass spectrometry.Nat. Protoc.2778–791. 10.1038/nprot.2007.95
11
DoidyJ.Van TuinenD.LamotteO.CorneillatM.AlcarazG.WipfD. (2012). The Medicago truncatula sucrose transporter family: characterization and implication of key members in carbon partitioning towards arbuscular mycorrhizal fungi.Mol. Plant51346–1358. 10.1093/mp/sss079
12
El-MesbahiM. N.AzcónR.Ruíz-LozanoJ. M.ArocaR. (2012). Plant potassium content modifies the effects of arbuscular mycorrhizal symbiosis on root hydraulic properties in maize plants.Mycorrhiza22555–564. 10.1007/s00572-012-0433-3
13
FahadS.HussainS.MatloobA.KhanF. A.KhaliqA.SaudS.et al (2015). Phytohormones and plant responses to salinity stress: a review.Plant Growth Regul.75391–404. 10.1007/s10725-014-0013-y
14
FeddermannN.FinlayR.BollerT.ElfstrandM. (2010). Functional diversity in arbuscular mycorrhiza – the role of gene expression, phosphorous nutrition and symbiotic efficiency.Fungal Ecol.31–8. 10.1016/j.funeco.2009.07.003
15
FernándezI.MerlosM.López-RáezJ. A.Martínez-MedinaA.FerrolN.AzcónC.et al (2014). Defense related phytohormones regulation in arbuscular mycorrhizal symbioses depends on the partner genotypes.J. Chem. Ecol.40791–803. 10.1007/s10886-014-0473-6
16
FesterT.FetzerI.BuchertS.LucasR.RilligM. C.HartigC. (2011). Towards a systemic metabolic signature of the arbuscular mycorrhizal interaction.Oecologia167913–924. 10.1007/s00442-011-2037-6
17
FritzM.JakobsenI.LyngkjærM. F.Thordal-ChristensenH.Pons-KühnemannJ. (2006). Arbuscular mycorrhiza reduces susceptibility of tomato to Alternaria solani.Mycorrhiza16413–419. 10.1007/s00572-006-0051-z
18
GamirJ.PastorV.CerezoM.FlorsV. (2012). Identification of indole-3-carboxylic acid as mediator of priming against Plectosphaerella cucumerina.Plant Physiol. Biochem.61169–179. 10.1016/j.plaphy.2012.10.004
19
GianinazziS.GollotteA.GollotteA.BinetM. N.Van TuinenD.RedeckerD.et al (2010). Agroecology: the key role of arbuscular mycorrhizas in ecosystem services.Mycorrhiza20519–530. 10.1007/s00572-010-0333-3
20
GiovannettiM.AvioL.BaraleR.CeccarelliN.CristofaniR.IezziA.et al (2012). Nutraceutical value and safety of tomato fruits produced by mycorrhizal plants.Br. J. Nutr.107242–251. 10.1017/S000711451100290X
21
GiovannettiM.MosseB. (1980). An evaluation of techniques for measuring vesicular arbuscular mycorrhizal infection in roots.New Phytol.84489–500. 10.1111/j.1469-8137.1980.tb04556.x.
22
GovindarajuluM.PfefferP. E.JinH.AbubakerJ.DoudsD. D.AllenJ. W.et al (2005). Nitrogen transfer in the arbuscular mycorrhizal symbiosis.Nature435819–823. 10.1038/nature03610
23
GuetherM.BalestriniR.HannahM.HeJ.UdvardiM. K.BonfanteP. (2009). Genome-wide reprogramming of regulatory networks, transport, cell wall and membrane biogenesis during arbuscular mycorrhizal symbiosis in Lotus japonicus.New Phytol.182200–212. 10.1111/j.1469-8137.2008.02725.x
24
HewittE. J. (1966). Sand and Water Culture Methods Used in the Study of Plant Nutrition. Technical Communication No. 22.London: Commonwealth Agriculture Bureau.
25
HodgeA.StorerK. (2015). Arbuscular mycorrhiza and nitrogen: implications for individual plants through to ecosystems.Plant Soil3861–19. 10.1007/s11104-014-2162-1
26
HohnjecN.ViewegM. E.PuhlerA.BeckerA.KusterH. (2005). Overlaps in the transcriptional profiles of Medicago truncatula roots inoculated with two different Glomus fungi provide insights into the genetic program activated during arbuscular mycorrhiza.Plant Physiol.1371283–1301. 10.1104/pp.104.056572
27
JungS. C.Martínez-MedinaA.López-RáezJ. A.PozoM. J. (2012). Mycorrhiza-induced resistance and priming of plant fefenses.J. Chem. Ecol.38651–664. 10.1007/s10886-012-0134-6
28
KaeverA.LandesfeindM.FeussnerK.MosblechA.HeilmannI.MorgensternB.et al (2014). MarVis-pathway: integrative and exploratory pathway analysis of non-targeted metabolomics data.Metabolomics11764–777. 10.1007/s11306-014-0734-y
29
KaschukG.KuyperT. W.LeffelaarP. A.HungriaM.GillerK. E. (2009). Are the rates of photosynthesis stimulated by the carbon sink strength of rhizobial and arbuscular mycorrhizal symbioses?Soil Biol. Biochem.411233–1244. 10.1016/j.soilbio.2009.03.005
30
KaurH.HeinzelN.SchottnerM.BaldwinI. T.GalisI. (2010). R2R3-NaMYB8 regulates the accumulation of phenylpropanoid-polyamine conjugates, which are essential for local and systemic defense against insect herbivores in Nicotiana attenuata.Plant Physiol.1521731–1747. 10.1104/pp.109.151738
31
KiersE. T.DuhamelM.BeesettyY.MensahJ. A.FrankenO.VerbruggenE.et al (2011). Reciprocal rewards stabilize cooperation in the mycorrhizal symbiosis.Science333880–882. 10.1126/science.1208473
32
LaparreJ.MalbreilM.LetisseF.PortaisJ. C.RouxC.BecardG.et al (2014). Combining metabolomics and gene expression analysis reveals that propionyl- and butyryl-carnitines are involved in late stages of arbuscular mycorrhizal symbiosis.Mol. Plant7554–566. 10.1093/mp/sst136
33
LiuJ.Maldonado-MendozaI.López-MeyerM.CheungF.TownC. D.HarrisonM. J. (2007). Arbuscular mycorrhizal symbiosis is accompanied by local and systemic alterations in gene expression and an increase in disease resistance in the shoots.Plant J.50529–544. 10.1111/j.1365-313X.2007.03069.x
34
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method.Methods25402–408. 10.1006/meth.2001.1262
35
López-RáezJ. A.VerhageA.FernándezI.GarcíaJ. M.Azcón-AguilarC.FlorsV.et al (2010). Hormonal and transcriptional profiles highlight common and differential host responses to arbuscular mycorrhizal fungi and the regulation of the oxylipin pathway.J. Exp. Bot.612589–2601. 10.1093/jxb/erq089
36
Martínez-MedinaA.PascualJ. A.LloretE.RoldánA. (2009). Interactions between arbuscular mycorrhizal fungi and Trichoderma harzianum and their effects on Fusarium wilt in melon plants grown in seedling nurseries.J. Sci. Food Agric.891843–1850. 10.1002/jsfa.3660
37
MarulandaA.AzcónR.Ruíz-LozanoJ. M. (2003). Contribution of six arbuscular mycorrhizal fungal isolates to water uptake by Lactuca sativa plants under drought stress.Physiol. Plant.119526–533. 10.1046/j.1399-3054.2003.00196.x
38
MensahJ. A.KochA. M.AntunesP. M.KiersE. T.HartM.BuckingH. (2015). High functional diversity within species of arbuscular mycorrhizal fungi is associated with differences in phosphate and nitrogen uptake and fungal phosphate metabolism.Mycorrhiza10.1007/s00572-015-0631-x[Epub ahead of print].
39
MishraS.TriptahiV.SinghS.PhukanU. J.GuptaM. M.ShankerK.et al (2013). Wound induced tanscriptional regulation of benzylisoquinoline pathway and characterization of wound inducible PsWRKY transcription factor from Papaver somniferum.PLoS ONE8:e52784. 10.1371/journal.pone.0052784
40
Mora-RomeroG. A.González-OrtízM. A.Quiroz-FigueroaF.Calderón-VázquezC. L.Medina-GodoyS.Maldonado-MendozaI.et al (2014). PvLOX2 silencing in common bean roots impairs arbuscular mycorrhiza-induced resistance without affecting symbiosis establishment.Funct. Plant Biol.4218–30. 10.1071/FP14101
41
PastorV.BalmerA.GamirJ.FlorsV.Mauch-ManiB. (2014). Preparing to fight back: generation and storage of priming compounds.Front. Plant Sci.5:295. 10.3389/fpls.2014.00295
42
Pedone-BonfimM.Da SilvaF.MaiaL. (2015). Production of secondary metabolites by mycorrhizal plants with medicinal or nutritional potential.Acta Physiol. Plant.371–12. 10.1007/s11738-015-1781-3
43
PhilippotL.RaaijmakersJ. M.LemanceauP.Van Der PuttenW. H. (2013). Going back to the roots: the microbial ecology of the rhizosphere.Nat. Rev. Microbiol.11789–799. 10.1038/nrmicro3109
44
PorcelR.ArocaR.AzcónR.Ruíz-LozanoJ. M. (2006). PIP aquaporin gene expression in arbuscular mycorrhizal Glycine max and Lactuca sativa plants in relation to drought stress tolerance.Plant Mol. Biol.60389–404. 10.1007/s11103-005-4210-y
45
PozoM. J.CordierC.Dumas-GaudotE.GianinazziS.BareaJ. M.Azcón-AguilarC. (2002). Localized versus systemic effect of arbuscular mycorrhizal fungi on defence responses to Phytophthora infection in tomato plants.J. Exp. Bot.53525–534. 10.1093/jexbot/53.368.525
46
PozoM. J.López-RáezJ. A.Azcón-AguilarC.García-GarridoJ. M. (2015). Phytohormones as integrators of environmental signals in the regulation of mycorrhizal symbioses.New Phytol.2051431–1436. 10.1111/nph.13252
47
Ruíz-LozanoJ. M.PorcelR.AzcónC.ArocaR. (2012). Regulation by arbuscular mycorrhizae of the integrated physiological response to salinity in plants: new challenges in physiological and molecular studies.J. Exp. Bot.634033–4044. 10.1093/jxb/ers126
48
SalvioliA.BonfanteP. (2013). Systems biology and “omics” tools: a cooperation for next-generation mycorrhizal studies.Plant Sci.203–204107–114. 10.1016/j.plantsci.2013.01.001
49
SardansJ.PeñuelasJ.Rivas-UbachA. (2011). Ecological metabolomics: overview of current developments and future challenges.Chemoecology21191–225. 10.1007/s00049-011-0083-5
50
SawadaY.AkiyamaK.SakataA.KuwaharaA.OtsukiH.SakuraiT.et al (2009). Widely targeted metabolomics based on large-scale MS/MS data for elucidating metabolite accumulation patterns in plants.Plant Cell Physiol.5037–47. 10.1093/pcp/pcn183
51
SchliemannW.AmmerC.StrackD. (2008). Metabolite profiling of mycorrhizal roots of Medicago truncatula.Phytochemistry69112–146. 10.1016/j.phytochem.2007.06.032
52
SchweigerR.BaierM. C.PersickeM.MullerC. (2014). High specificity in plant leaf metabolic responses to arbuscular mycorrhiza.Nat. Commun.53886. 10.1038/ncomms4886
53
SelosseM. A.BessisA.PozoM. J. (2014). Microbial priming of plant and animal immunity: symbionts as developmental signals.Trends Microbiol.22607–613. 10.1016/j.tim.2014.07.003
54
SmithC. A.WantE. J.O’mailleG.AbagyanR.SiuzdakG. (2006). XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification.Anal. Chem.78779–787. 10.1021/ac051437y
55
SmithF. A.SmithS. E. (2015). How harmonious are arbuscular mycorrhizal symbioses? Inconsistent concepts reflect different mindsets as well as results.New Phytol.2051381–1384. 10.1111/nph.13202
56
SmithS. E.ReadD. (2008). Introduction, in Mycorrhizal Symbiosis3rd Edn. London: Academic Press1–9. 10.1016/B978-012370526-6.50002-7
57
SmithS. E.SmithF. A. (2011). Roles of arbuscular mycorrhizas in Plant nutrition and growth: new paradigms from cellular to ecosystem scales.Annu. Rev. Plant Biol.62227–250. 10.1146/annurev-arplant-042110-103846
58
SongY. Y.YeM.LiC. Y.WangR. L.WeiX. C.LuoS. M.et al (2013). Priming of anti-herbivore defense in tomato by arbuscular mycorrhizal fungus and involvement of the jasmonate pathway.J. Chem. Ecol.391036–1044. 10.1007/s10886-013-0312-1
59
SteinkellnerS.Hage-AhmedK.Garcia-GarridoJ. M.IllanaA.OcampoJ. A.VierheiligH. (2012). A comparison of wild-type, old and modern tomato cultivars in the interaction with the arbuscular mycorrhizal fungus Glomus mosseae and the tomato pathogen Fusarium oxysporum f. sp. lycopersici.Mycorrhiza22189–194. 10.1007/s00572-011-0393-z
60
StumpeM.CarsjensJ. G.StenzelI.GobelC.LangI.PawlowskiK.et al (2005). Lipid metabolism in arbuscular mycorrhizal roots of Medicago truncatula.Phytochemistry66781–791. 10.1016/j.phytochem.2005.01.020
61
TisserantE.MalbreilM.KuoA.KohlerA.SymeonidiA.BalestriniR.et al (2013). Genome of an arbuscular mycorrhizal fungus provides insight into the oldest plant symbiosis.Proc. Natl. Acad. Sci. U.S.A.11020117–20122. 10.1073/pnas.1313452110
62
TschaplinskiT. J.PlettJ. M.EngleN. L.DeveauA.CushmanK. C.MartinM. Z.et al (2014). Populus trichocarpa and Populus deltoides exhibit different metabolomic responses to colonization by the symbiotic fungus Laccaria bicolor.Mol. Plant Microbe Interact.27546–556. 10.1094/mpmi-09-13-0286-r
63
VierheiligH.SchweigerP.BrundrettM. (2005). An overview of methods for the detection and observation of arbuscular mycorrhizal fungi in roots.Physiol. Plant.125393–404. 10.1111/j.1399-3054.2005.00564.x
64
VosC. M.TesfahunA. N.PanisB.De WaeleD.ElsenA. (2012). Arbuscular mycorrhizal fungi induce systemic resistance in tomato against the sedentary nematode Meloidogyne incognita and the migratory nematode Pratylenchus penetrans.Appl. Soil Ecol.611–6. 10.1016/j.apsoil.2012.04.007
65
WalderF.NiemannH.NatarajanM.LehmannM. F.BollerT.WiemkenA. (2012). Mycorrhizal networks: common goods of plants shared under unequal terms of trade.Plant Physiol.159789–797. 10.1104/pp.112.195727
66
WangY.ChantreauM.SiboutR.HawkinsS. (2013). Plant cell wall lignification and monolignol metabolism.Front. Plant Sci.4:220. 10.3389/fpls.2013.00220
67
WasternackC.HauseB. (2013). Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in Annals of Botany.Ann. Bot.1111021–1058. 10.1093/aob/mct067
68
YazakiK. (2006). ABC transporters involved in the transport of plant secondary metabolites.FEBS Lett.5801183–1191. 10.1016/j.febslet.2005.12.009
69
ZengY.GuoL. P.ChenB. D.HaoZ. P.WangJ. Y.HuangL. Q.et al (2013). Arbuscular mycorrhizal symbiosis and active ingredients of medicinal plants: current research status and prospectives.Mycorrhiza23253–265. 10.1007/s00572-013-0484-0
70
ZhangN.VenkateshwaranM.BoersmaM.HarmsA.Howes-PodollM.Den OsD.et al (2012). Metabolomic profiling reveals suppression of oxylipin biosynthesis during the early stages of legume-rhizobia symbiosis.FEBS Lett.5863150–3158. 10.1016/j.febslet.2012.06.046
71
ZiedanE.-S.ElewaI.MostafaM.SahabA. (2011). Application of Mycorrhizae for controlling root diseases of sesame.J. Plant Prot. Res.51355–361. 10.2478/v10045-011-0058-0
72
ZouariI.SalvioliA.ChialvaM.NoveroM.MiozziL.TenoreG. C.et al (2014). From root to fruit: RNA-Seq analysis shows that arbuscular mycorrhizal symbiosis may affect tomato fruit metabolism.BMC Genomics15:221. 10.1186/1471-2164-15-221
Summary
Keywords
arbuscular mycorrhiza, metabolomics, Funneliformis mosseae, Rhizophagus irregularis, oxylipins
Citation
Rivero J, Gamir J, Aroca R, Pozo MJ and Flors V (2015) Metabolic transition in mycorrhizal tomato roots. Front. Microbiol. 6:598. doi: 10.3389/fmicb.2015.00598
Received
02 April 2015
Accepted
01 June 2015
Published
23 June 2015
Volume
6 - 2015
Edited by
Essaid Ait Barka, Reims University, France
Reviewed by
Richard S. Winder, Natural Resources Canada, Canada; Prasun Ray, The Samuel Roberts Noble Foundation, USA
Copyright
© 2015 Rivero, Gamir, Aroca, Pozo and Flors.
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) or licensor 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: María J. Pozo, Department of Soil Microbiology and Symbiotic Systems, Estación Experimental del Zaidín – Consejo Superior de Investigaciones Científicas, Prof Albareda 1, 18008 Granada, Spain, mjpozo@eez.csic.es; Víctor Flors, Metabolic Integration and Cell Signaling Laboratory, Associated Unit UJI-CSIC, Plant Physiology Section, Department of Ciencias Agrarias y del Medio Natural, Universitat Jaume I, Avenida Vicente Sos Baynat, 12071 Castellón, Spain, flors@uji.es
This article was submitted to Plant-Microbe Interaction, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.