Seed Metabolism and Pathogen Resistance Enhancement in Pisum sativum During Colonization of Arbuscular Mycorrhizal Fungi: An Integrative Metabolomics-Proteomics Approach

Pulses are one of the most important categories of food plants, and Pea (Pisum sativum L.) as a member of pulses is considered a key crop for food and feed and sustainable agriculture. Integrative multi-omics and microsymbiont impact studies on the plant's immune system are important steps toward more productive and tolerant food plants and thus will help to find solutions against food poverty. Didymella pinodes is a main fungal pathogen of pea plants. Arbuscular mycorrhizal fungi (AMF) promote plant growth and alleviate various stresses. However, it remained unclear as to how the AMF effect on seed metabolism and how this influences resistance against the pathogen. This study assesses the AMF impacts on yield components and seed quality upon D. pinodes infection on two different P. sativum cultivars, susceptible versus tolerant, grown in pots through phenotypic and seed molecular analyses. We found that AMF symbiosis affects the majority of all tested yield components as well as a reduction of disease severity in both cultivars. Seeds of mycorrhizal pea plants showed strong responses of secondary metabolites with nutritional, medicinal, and pharmaceutical attributes, also involved in pathogen response. This is further supported by proteomic data, functionally determining those primary and secondary metabolic pathways, involved in pathogen response and induced upon AMF-colonization. The data also revealed cultivar specific effects of AMF symbiosis that increase understanding of genotype related differences. Additionally, a suite of proteins and secondary metabolites are presented, induced in seeds of P. sativum upon AMF-colonization and pathogen attack, and possibly involved in induced systemic resistance against D. pinodes, useful for modern breeding strategies implementing microsymbionts toward increased pathogen resistance.


INTRODUCTION
A sustainable supply of crops to feed human societies could be the most important humanitarian action and its impact on declining of current and future food restrictions on this planet should be more sensed especially in regions with the highest food insecurity. Hence, the stable food supply is completely dependent on the sustainable production of crops by increasing yield and quality improvement through the development of natural potentials locally. To supply the food sustainably, large-scale production of pulses or legumes is inevitable (Fernandez-Aparicio et al., 2010). Pea (Pisum sativum L.) is one of the major legumes in the world (FAOSTAT 1 , updated 2019). Moreover, the use of microsymboints is in line with sustainable agriculture (Sayeed Akhtar et al., 2011).
Considering the importance of pulses as a source of vegetable protein in food basket and global need to improve their quantity and quality sustainably should be the most urgent need particularly in regions where people are suffering from the highest poverty rate. Thus, the study of microsymbionts which affect yield components of legumes especially on seed and upon stress conditions not only due to better use of natural growth promoters in agricultural systems but also most probably reduces the losses caused by abiotic and biotic stresses. Arbuscular mycorrhizal fungi (AMF) have interactions with a broad range of plants in nature Read, 2008, as cited in Wehner et al., 2010). AMF inoculation of crops due to an increase in yield and growth (Hamel and Plenchette, 2007, as cited in Horii and Ishii, 2014). Promoting tolerance against water deficit and growth has been observed in papaya plants associated with AMF (Cruz et al., 2000, as cited in Horii andIshii, 2014).
Chemical properties, physiology, and seed development of soybean could be influenced by AMF root colonization (Bethlenfalvay et al., 1997). Bona et al. (2016) noted to positive effects of AMF root colonization in prior studies on quality and yield in strawberry (Castellanos-Morales et al., 2010;Castellanos-Morales et al., 2012;Lingua et al., 2013;Bona et al., 2015) and saffron (Aimo et al., 2010). Also, the inoculation of Allium sativum with AMF has enhanced the yield and growth (Borde et al., 2009, as cited in Bona et al., 2016. The formed symbiosis between AMF and maize (Zea mays) affected its seed proteome through modulation of enzymes related to functional categories such as stress, nucleotide metabolism, energy, storage, and development (Bona et al., 2016). Colonization of Amorpha fruticosa roots by Glomus mosseae has been modified the root proteome (Song et al., 2015). The G. mosseae by promoting cell integrity and osmotic stress depletion in wheat protected the root system against drought stress (Bernardo et al., 2017). The protein synthesis, metal handling, RNA metabolism, and reactive oxygen species (ROS) are influenced in leaf proteome of pea plants inoculated with AMF . In leaves of AMF colonized Medicago truncatula plants, not only genes related to synthesis of jasmonic acid (JA), flavonoid, abscisic acid (ABA), and terpenoids were up-regulated but also flavonoids and anthocyanins were intensified (Adolfsson et al., 2017).
Biotic stresses such as plant diseases are serious threats for food plants and consequently for food sources (Strange and Scott, 2005). One of the main diseases known in P. sativum is ascochyta blight caused by D. pinodes (Moussart et al., 1998). D. pinodes as the most abundant pathogen among seeds of P. sativum (Deneufbourg et al., 1994, as cited in Moussart et al., 1998) causes the poor quality and yield reduction (Khan et al., 2013). Under the field conditions, disease management of this pathogen is difficult in comparison with the control of other diseases upon P. sativum (Khan et al., 2013).
Chemical treatment of crop seeds is considered one of the major strategies in integrated control of ascochyta blight (Davidson and Kimber, 2007). On the other hand, using chemical compounds is not only a sustainable solution for the control of plant diseases but also acts as a destroyer of the environment or against sustainable agriculture systems. Therefore, symbiosis formation between food plants and microsymbionts with bio-control potential could be used as stable and safe disease management. Promoting of resistance against biotic stress in plants by AMF is proved (Azcoń-Aguilar and Barea, 1996, as cited in Rebollo Couto et al., 2013. Dehne (1982) reviewed the former studies on the role of mycorrhizae for effective protection of varied plant species against fungal pathogens such as Pythium ultimum (Stewart and Pfleger, 1977), Phytophthora cinnamomi (Bäertschi et al., 1982), P. megasperma (Chou and Schmitthenner, 1974), P. parasitica (Schenck et al., 1977;Davis and Menge, 1980), F. oxysporum cucumerinum , F. oxysporum lycopersici , Rhizoctonia solani (Stewart and Pfleger, 1977), Cylindrocladium scoparium (Barnard, 1977), Phoma terrestris (Becker, 1976), Pyrenochaeta terrestris (Safir, 1968), Olpidium brassicae Schoenbeck and Dehne, 1981), and Thielaviopsis basicola (Baltruschat and Schonbeck, 1972;Schoenbeck and Dehne, 1977).
Replacing common chemical fertilizers and toxicant with AMF has been focused in line with sustainable agriculture (Gange et al., 1999;Harrier and Watson, 2004;Wang et al., 2008;Lecomte et al., 2011, as cited in Rebollo Couto et al., 2013. Additionally, AMF could be applied as a bio-control agent and growth promoter (Sayeed Akhtar et al., 2011). It was found that AMF, directly and indirectly, inhibit the prevalence of fungal root pathogens (Wehner et al., 2010). The role of AMF in strengthening the host plant upon pathogen infection by upregulation of defensive pathways and signals has been described (Haneef Khan et al., 2010).
Some studies have investigated the AMF impacts on the plant's growth at the maturity stage (Bethlenfalvay et al., 1994;Bethlenfalvay et al., 1997;Al-Karaki and Clark, 1999). Prior researches have mainly investigated the role of AMF on soilborne pathogens in non-legumes and few investigations have attended to AMF impacts as a bio-control agent on aboveground fungal pathogens including ascochyta blight of P. sativum caused by D. pinodes. Further, most studies have analyzed the growth parameters of host plants upon mycorrhizal colonization rather than yield components. For instance, recent studies that have investigated the D. pinodes infection of P. sativum upon different treatments including mycorrhizal plants, were focusing on the leaves where no clear disease reduction in mycorrhizal treatments compared with non-mycorrhizal was reported Turetschek et al., 2017). In general, only a few studies exist, using multilevel-omics approaches by integrating plant morphological and molecular phenotyping (Zivy et al., 2015).
Here, we assess the integrative effects of AMF on several yield components, growth parameters, seed secondary metabolome, and -proteome upon tolerant and susceptible pea cultivars against D. pinodes as an aboveground pathogen in pots. This study answers the following research questions: I) Does AMF affect above-belowground growth and yield parameters in pea plants? II) How does AMF promote the productivity and seed quality of P. sativum from the phenotypic-omics perspective? III) Has AMF-symbiosis a protecting effect on seed quality of tolerant and susceptible pea cultivars upon pathogen attack? IV) Do AMF influence the seed secondary metabolome and proteome? and V) How do AMF affect the response of the seed metabolism upon pathogen infection?

Performing Experiments and Planting
In line with processes carried out before (Ranjbar Sistani et al., 2017) and according to Begum et al. (2001), the disinfected pea seeds (ethanol 95%) were washed by ultrapure water and then, were kept (20 min) in bleach 5% and rinsed by ultrapure water several times. Finally, the soaked seeds (4 h) were considered for planting (Begum et al., 2001). The pre-germinated pea seeds during 3 days on the sterile perlite-vermiculate substrate were planted in pots containing sterilized soil with the following chemical properties: K 2 PO 120 mg/L, N 7 mg/L, pH 5.6 and P 13 mg/L . The pots were kept under the following conditions: 600 μmol m −2 s −1 lighting (14-h day/10-h night), humidity 60%-70%, and 22°C day/16°C night (Ludidi et al., 2007;Larrainzar et al., 2014). Treatments consisting of tolerant [cultivar Protecta (cv. Pr)] and susceptible [cultivar Messire (cv. Me)] genotypes, symbiont (M = mycorrhizal, NM = non˗mycorrhizal) and pathogen (infected/diseased = I and uninfected/healthy = U) were arranged in pots. Four plants per biological replicate (three biol. replicates) and two technical replicates per each treatment were considered. For metabolomics, proteomics, and phenotypic evaluations of seed, a hundred seeds per each biological replicate were considered. A modified recipe of B&D (Broughton and Dilworth, 1970)   . An overview of the experimental setup can be found in Figure 1.

Pathogen Inoculation and Disease Severity Estimation
According to our previous study (Ranjbar Sistani et al., 2017) and in detail: the prepared pathogen (D. pinodes) by Rubiales Lab in Spain was multiplied on Potato Dextrose Agar (PDA) and incubated upon 12-h photoperiod and 22°C (Davidson et al., 2012). Incubated colonies of pathogen after a week was used to make fungal suspension by adding Milli-Q (ultrapure) water and scratching of colonies upon sterile conditions (Zimmer and Sabourin, 1986;Carrillo et al., 2013). The autoclaved cheesecloth was applied to filter the prepared suspension and the filtered pathogen suspension with concentration 3 × 10 5 spores/ml was selected for inoculation (Zimmer and Sabourin, 1986;Carrillo et al., 2013). The fungal suspension plus TWEEN 20 (120 ml/100 ml) was handled to inoculate the pea seedlings after 20 days of planting separately from healthy (uninfected) treatments that were sprayed by a mixture of water and TWEEN 20 (Carrillo et al., 2013). Then, inoculated seedlings were kept under transparent covers during 48 h (Carrillo et al., 2013;Okorska et al., 2014) and after this period, infected seedlings were set apart for one week (12-h photoperiod, 21 ± 2°C) (Garry et al., 1998). Afterward, all seedlings were arranged as described in Performing Experiments and Planting, the aerial parts of infected pea plants by D. pinodes were regularly and periodic investigated by using a USB digital microscope (25X-400X, BMSCI, Japan) to image and computing the lesions (Ranjbar Sistani et al., 2017). The captured dimensions of lesions were applied to calculate the area of lesions and finally disease severity assay (Hwang et al., 2006;Carrillo et al., 2013). To complete the disease severity assessments, the infection rate of seeds in diseased pea plants was specified through the paper towel and PDA methods (Xue et al., 1996;Gaurilˇcikien˙e et al., 2012;Mahmoud et al., 2013) besides the area of lesions.

Sample Preparation
The pea plants were harvested when pods and seeds were ripe (BBCH 81-88). The root samples were prepared for investigation of mycorrhizal root colonization when the seeds were matured (Al-Karaki and Clark, 1999;Horii and Ishii, 2014). According to the modified method of Jin et al., 2013, soil particles were removed from the root samples using tap water, washed again with Milli-Q (ultrapure) water, and then dried with a paper towel. For each pot, subsamples of about 0.1-0.2 g were separated from the root system, freshly weighted (FW1) and immediately the fresh weight of remaining root was recorded (FW2) besides determining of the dry weight of that larger sample (DW2) [ (Giovannetti and Mosse, 1980;Plant & Mycorrhiza Root Lengths 2 )]. The small subsamples of pea roots were boiled (10% KOH, 15 min) and the boiled root samples were washed by using water (Vierheilig et al., 1998). For staining, the cleaned root samples were put in boiling inkacetic acid (5%, 3 min), and then, the stained root samples were washed (20 min) with acidified water by acetic acid (Vierheilig et al., 1998).

Mycorrhizal Colonization Assay
The mycorrhizal colonization of P. sativum roots was studied based on the gridline intersect method (Giovannetti and Mosse, 1980). On average, 20 stained root samples per pot were investigated. The count of intersections between lines, NM roots (R1), and M roots (R2) were recorded (INVAM 2 ). The percentage of root colonization by AMF was calculated per each root sample by the following formula (Giovannetti and Mosse, 1980;Jin et al., 2013;INVAM 2

Investigation of Physical Properties and Yield of Pea Seeds
In compliance with prior work (Ranjbar Sistani et al., 2017), the total weight of pea seeds per each pot was recorded and the seed yield (kg/ha) formula that was mentioned in Sajid et al. (2012) work was adapted to determine the seed yield per treatment (kg pot −1 ). The physical quality of pea seeds was evaluated through quantification of absorption and hardness rates of pea seeds per treatment that were kept in water (16 h) and then were assessed by hydration coefficient (%) and non-soakers (%) formulas (Abdelgani et al., 1999). The vigor index of pea seeds was quantified based on seed germination rate and seedling length measurements by using rolled paper towels (Abdul Baki and Anderson, 1973;Farrag and Moharam, 2012).

Extraction of Seed Secondary Metabolites and Seed Metabolomics by NanoESI LC-MS/MS
A modified protocol of De Vos et al. (2007) was applied for seed metabolomics as described before (Ranjbar Sistani et al., 2017): the mixed ground pea seed (100 mg) with 1 ml methanol (80%) were arranged in a filled ultrasonic bath with ice and water (10 min) followed by centrifugation (10 min, 21,000 × g). The resulting supernatant was dried in a vacuum centrifuge concentrator and then each dried sample was mixed with 50 ml [methanol (50%) + FA (0.1%)] and was centrifuged. After that, the reserpine (3ml, Sigma-Aldrich) was added to each sample which was diluted (1:10) by methanol (5%) in FA (0.1%) and were centrifuged again. For mass spectrometry, 20 ml of each supernatant was used. The NanoESI LC-MS/MS was carried out by using LTQ-Orbitrap XL Hybrid Ion Trap-Orbitrap Mass Spectrometer (Thermo Fisher Scientific, Germany) with HPLC column and same parameters setting to proteomics analysis except for required signal (min 50,000) and scan range (130-1,800 m/z).

Profiling, Quantification, and Identification of Seed Metabolites
Annotation, quantification, and identification of the extracted metabolites were carried out as previously described (Ranjbar Sistani et al., 2017), with some changes: for identification and quantification of metabolites, Xcalibur (version 2.3.26, Thermo Fisher Scientific Inc.) was used for conversion of mass spectrometry output (RAW files) to CDF files. The resulting CDF files were used as input of MET-COFEA (Zhang et al., 2014). Then, the output files of MET-COFEA were loaded into the MET-XAlign (Zhang et al., 2015) tool. In parallel, the created mzXML files from the conversion of RAW files by MassMatrix (version 3.9, Case Western Reserve University 3 ) were processed by using ProtMAX 2012_rev.2.14 (Egelhofer et al., 2013) as well as a spectral survey (Wang et al., 2016) and peak integration (Xcalibur 2.2, Thermo Scientific). The resulting data from MET-XAlign and ProtMAX outputs were used to complete the metabolite identification by using HMDB 4 (Wishart et al., 2018), METLIN 5 (Smith et al., 2005;Guijas et al., 2018), KNApSAcK family 6 (Afendi et al., 2012). The metabolite candidates were finalized through literature reviews.

Statistical Data Analysis
The normality of data was checked by using the Shapiro-Wilk test and to find the significant differences in treatments, ANOVA with Tukey HSD (P-value < 0.05) as a post hoc test, and the Kruskal-Wallis test were applied by using STATGRAPHICS Centurion 18 7 . To study the major impacts of factors (microsymbiont, cultivar, and pathogen) and their interactions for phenotypic data, multifactor (three-way) ANOVA plus Tukey HSD procedure (P-value < 0.05) was performed in STATGRAPHICS Centurion 18 7 . To analyze the correlations among all phenotypic variables, Pearson correlation coefficients through STATGRAPHICS Centurion 18 7 were calculated. Only proteins and metabolites that were reproducibly detected in all three biol. replicates (averages of two technical replicates per biological replicate) of at last one treatment were considered for quantification (Ranjbar Sistani et al., 2017). Filling of missed values base on a prior distribution followed by z-transformation was performed for independent component analysis (ICA) and cluster analysis by using COVAIN (Sun and Weckwerth, 2012) (Updated 2019; MATLAB R2017a, The MathWorks, Inc., Natick, Massachusetts, United States.). The q-values were calculated using an online Shiny application 8 of q-value R package 9 (Storey, 2002;Storey, 2003;Storey and Tibshirani, 2003;Storey et al., 2004;R Core Team, 2013).

Mycorrhizal Colonization of P. sativum Roots
The mean mycorrhizal colonization in roots of inoculated pea plants was ranged between 64.78% and 80.13% (Table S1). Although the mean AMF colonization between two cultivars was not significantly different, cultivar Protecta (cv. Pr) indicated higher colonization (11%) compared to cultivar Messire (cv. Me) ( Table S1). Similarly, mycorrhizal colonization was not significantly affected by pathogen infection (Table S1), however, uninfected pea plants were more colonized (12%) by AMF than infected plants (Table S1). Although, TMRL and TPRL were not significantly different between the two cultivars, cv. Pr had higher root lengths (mycorrhizal-colonized and not) in comparison with cv. Me (Table S1). Root (FW and DW) were strongly and positively correlated with TMRL and TPRL (Table  S2). Also the symbiont factor showed a significant impact on root fresh and dry weights ( Table 1).

AMF Impact on Yield Components, Growth Parameters, and Seed Physical Properties
The following results are according to multifactor ANOVA to determine the variability of phenotypic parameters under the influence of three factors (symbiont, cultivar, and pathogen) besides comparison of treatments (M vs. NM) per each phenotypical variable. Although symbiont (AMF) factor had a significant effect on pod size, pod weight, seed (FW and DW) per seed, seed (FW) per plant, thousand seed weight TSW (FW and DW), seed yield, root (FW and DW), and shoot (FW) ( Table 1), the yield components showed no significant changes between M and NM treatments (Tables 2 and 3). Cultivar and pathogen factors separately indicated significant effects on pod number per plant, pod weight, shoot length, number of flowers, number of nodes, shoot (FW and DW), and root (FW and DW) ( Table 1). Cultivar and symbiont factors individually showed a significant influence on pod size, seed (FW) per seed, and TSW (FW) ( Table  1). Also, the pathogen factor had a significant impact on seed (FW) per plant, seed yield, seed (DW) per plant, and seed number per plant ( Table 1). Among seed physical properties, cultivar and pathogen factors separately demonstrated a significant effect on non-soaker and vigor indexes respectively ( Table 1). Among yield components, symbiont, cultivar, and pathogen factors had significant effects separately on pod weight while the triple interaction effect of these factors on seed (FW) per plant and seed yield were significant (Table 1). Additionally, a significant and positive correlation was observed between seed number per plant and mycorrhizal colonization percentage while the none-soaker index was significantly and negatively correlated with AMF colonization, TMRL, and TPRL (Table S2). Although the majority of yield components showed no significant changes between M and NM treatments (Tables 2 and 3) but the seed (FW and DW) per seed and TSW (FW and DW) in cv. Pr were significantly enhanced in M plants vs. NM under healthy and diseased conditions ( Table 2). Furthermore, root (FW and DW) were significantly enhanced in mycorrhizal (M) vs. non˗mycorrhizal (NM) plants upon healthy and diseased across both cultivars ( Table 3).

Mycorrhizal Colonization Effect on Disease Severity
The disease severity of P. sativum plants infected by D. pinodes was determined by the evaluation of the seed infection level and area of lesions. Although the mean area of lesions between the two cultivars was not significantly different, cv. Pr had lower lesions area (Table 4) and significantly lower levels of seed infection (Tukey HSD test, p < 0.05, nearly 2 folds) compared with cv. Me (Table 4). We found that AMF treatment of cv. Me resulted in a significantly lower level (Tukey HSD test, p < 0.05, nearly 2 folds) of seed infection (Table 4). Also, the inoculation of AMF significantly decreased the mean area of lesions in both infected cultivars ( Table 4).

The Effect of AMF on Pea Seed Metabolites Under Infection Caused by D. pinodes
In total, 47 metabolites from 12 chemical classes and families were identified. The identified metabolites mostly belonged to glycerophospholipids, flavonoids, isoflavonoids, prenol lipids, carboxylic acids and derivatives, and organooxygen compounds (Table S3A). Glycerophosphoinositol phosphates (PIPs) as well as 2',4',5-Trihydroxy-7,8-[2-(1-methylethenyl) furo] isoflavone and lupinisoflavone A could not be assigned separately (due to identical chemical formula, retention time and m/z) and were therefore assigned as compound groups of PIPs and putative metabolites (*2',4',5-Tri_Lup_Isoflavon), respectively (Table S3A). Compared to our previous study, we increased our metabolite identification due to the additional use o f t h e KN A p S A c K f a m i l y (h t t p : / / k a n a y a . n a i s t . j p / KNApSAcK_Family/) database. Eighty-seven percent of identified metabolites were significantly changed (ANOVA, Kruskal-Wallis, Tukey HSD tests, P-value < 0.05, and ≥ twofold change) among treatments (M vs. NM of diseased plants across both cultivars) ( Table S3B). Most of the metabolites that significantly increased in M vs. NM treatments were induced upon pathogen infection and mainly in cv. Pr (Table S3B). For example, (Z)-3-Oxo-2-(2-pentenyl)-1-cyclopenteneacetic acid, cyanidin 3-sophoroside 5-glucoside, luteone 7-glucoside, and stachyose were significantly enhanced in the seed metabolome of M vs. NM of cv. Pr under pathogen infection (Table S3B).
In contrast, L-2-Amino-3-(oxalylamino) propanoic acid (or L -3 -A m i n o -2 -( o x a l y l a m i n o ) p r o p a n o i c a c i d ) , N -(Carbethoxyacetyl)-4-chloro-L-tryptophan, and vignatic acid A were only significantly accumulated in M vs. NM treatments of cv. Me against pathogen infection ( Table S3B).

Influence of AMF-Root Colonization on Pea Seed Proteome
Only proteins without missing values (LFQ intensities) across replicates of the same treatment were analyzed statistically. The resulting 1,325 proteins were functionally categorized by applying the Mercator pipeline (Lohse et al., 2014) and MapMan tool (Thimm et al., 2004) (Table S4). Table S4 illustrates the comprehensive seed proteomics data including categorized proteins, protein descriptions, mapped functional classes, statistical analyses, and fold change (FC) ratio for all treatments. Overall, remarkable changes (ANOVA, Kruskal-   Wallis, Tukey HSD test, and P-value < 0.05 and ≥ two-fold change) were found in 25% of mapped proteins by comparing M vs. NM under diseased and healthy conditions in both cultivars ( Table S4).
Among mapped functional categories that were significantly changed (P-value < 0.05 and ≥ two-fold change) in M vs. NM, the category of protein degradation had the highest number of proteins ( Table S4). Most of the identified proteins were involved in protein degradation and thus several of those were found with significantly increased levels in M vs. NM treatments ( Table S4). Among M vs. NM treatments under healthy conditions, protein degradation, protein synthesis, and stress categories in cv. Pr and plastid (PS) in cv. Me had the largest number of proteins changed upon M (Table S4). Protein functional categories such as amino acid metabolism (Figure 2; frv2_55445, frv2_85107, frv2_60349, frv2_112647, frv2_48897 and frv2_83949), miscellaneous (misc) (Figure 2; frv2_61709, frv2_111740, Q9FN08, frv2_46536, frv2_82443, P93479, frv2_103111, frv2_103598, frv2_78262, and frv2_83016) and redox (Table S4; frv2_125364, frv2_76239, frv2_47501, frv2_76239, frv2_77791, frv2_81620, and Q9FF55) showed a significant increase in M vs. NM upon disease and healthy conditions in both cultivars.
Under pathogen infection, several categories including DNA (Table S4; frv2_98925 and frv2_50256), mitochondrial electron transport (Table S4; frv2_77583, frv2_103383, frv2_81222, and frv2_53151), nucleotide metabolism (Table S4; frv2_87000, frv2_81447, frv2_112313, and frv2_48788), oxidative pentose phosphate (OPP) ( Table S4; Q43848 and Q9SZE1), protein amino acid (aa) activation (Table S4; Q9ZPI1 and frv2_121810) and signaling (Figure 2; frv2_113758 and Q96453) were remarkably enhanced in M vs. NM of both cultivars. In contrast, the protein synthesis category (Figure 2; frv2_75059, frv2_46668, frv2_75112, frv2_54389, and frv2_76505) indicated a notable accumulation in M vs. NM exclusively under healthy conditions. In M vs. NM, the categories development (Figure 2; frv2_47399, frv2_88457, and frv2_94799) and PS (Table S4; Q43848, frv2_80442, frv2_77556, frv2_84713, and Q9SHE8) categories only in cv. Me besides C1-metabolism (Table S4; frv2_90362 and frv2_86521) and glycolysis (Table S4; frv2_114679, frv2_112199, frv2_47599, and Q88C93) uniquely in cv. Pr were significantly pronounced upon disease and health conditions. Comparing of I vs. U upon M and NM treatments in both cultivars indicated that levels of several proteins engaged in amino acid metabolism, cell, DNA, lipid metabolism, metal handling, misc, mitochondrial electron transport, nucleotide metabolism, photosynthesis (PS), redox and signaling were significantly reduced while abundances of some proteins involved in proteins synthesis were remarkably enhanced (Table S4). Interestingly, some levels of proteins involved in protein degradation and stress categories were either notably accumulated or reduced in I vs. U (Table S4). Upon pathogen infection, the largest functional categories were protein synthesis (increased levels) and protein degradation (decreased levels) ( Table S4). Some proteins assigned to fermentation, major CHO metabolism, and protein targeting were remarkably accumulated only in the seed proteome of NM treatments of both cultivars against pathogen infection. In contrast, some proteins of C1-metabolism, development, protein posttranslational modification, TCA cycle, biodegradation of xenobiotics, and transport categories showed a significant reduction ( Table S4).
Comparison of cultivars across all treatments including M and NM under disease and healthy conditions demonstrated that several proteins of N-metabolism, secondary metabolism, and signaling in cv. Me besides some proteins of amino acid metabolism and misc in cv. Pr were significantly (p < 0.05) and distinctively enhanced (≥2-fold) ( Table S4). Moreover, several levels of proteins involved in cofactor and vitamin metabolism were significantly intensified in both cultivars of NM treatments against pathogen infection ( Table S4). Numbers of identified proteins from OPP and protein post-translational modification classes in cv. Me in addition to some TCA cyclerelated proteins in cv. Pr showed a significant enhancement in NM treatments upon pathogen infection ( Table S4).
The proteins from main functional classes (including >2 proteins) with significantly increased fold change ratio (p < 0.05, FC ≥ 2) in M vs. NM under healthy and disease conditions of both cultivars are presented in Figure S1. Amino acid metabolism, misc and redox categories were only the main protein functional groups with enhanced FC significantly upon all paired comparisons of M vs. NM ( Figure S1A). The most increased FC in M vs. NM belonged to protein synthesis in cv. Pr under healthy conditions ( Figure S1A). Hormone metabolism and OPP had the highest increased FC in M treatments of cv. Me and cv. Pr respectively against pathogen infection ( Figure S1A). Under pathogen infection, proteins associated with mitochondrial electron transport, nucleotide metabolism, OPP, protein amino acid (aa) activation and DNA were exclusively and largely accumulated in M vs. NM of both cultivars while, protein synthesis was uniquely promoted in M treatments of cultivars under healthy conditions ( Figure S1A). When comparing M vs. NM, a significant increase of C1-metabolism, glycolysis, major CHO metabolism, and minor CHO metabolism were observed only in cv. Pr besides the notable accumulation of proteins related to development, PS, and TCA cycle individually in cv. Me were observed ( Figure S1A). Protein degradation, protein synthesis, and stress functional classes were strongly increased in seed proteome of all I vs. U treatments in both cultivars, and on average among these categories, protein synthesis had the highest summed FC ratio ( Figure S1B). Under disease conditions, engaged proteins in cell and biodegradation of xenobiotics illustrated noteworthy enhancement exclusively in cv. Pr while development, protein post-translational modification, and redox were highly intensified only in cv. Me ( Figure S1B). The proteins included in signaling, stress, and TCA cycle were significantly promoted in cv. Pr compared with cv. Me upon M and NM treatments of both healthy and diseased ( Figure S1C). Among these protein functional categories, stress-related proteins were highly increased in Pr vs. Me than proteins mapped in signaling and TCA cycle ( Figure S1C).

Integrated Analysis of Seed Metabolomics-Proteomics Data
To determine the metabolites and proteins with the highest impact on discrimination of treatments under stress conditions caused by pathogen, an ICA was performed (Figure 3). Only intensities of identified seed metabolites and proteins with significant alteration (ANOVA, Kruskal-Wallis, Tukey HSD tests, P-value < 0.05, and ≥ two-fold change) in M vs. NM upon diseased conditions of both cultivars were analyzed. The integrated seed metabolomics and proteomics data displayed distinctive separations of cultivars (Me vs. Pr) on IC2 and mycorrhizal symbiont (M vs. NM) on IC3 (Figure 3). Figure 3 shows the ten highest positive and negative loadings responsible for these discriminations are depicted in separate bar plots. Kudzusaponin SA1 and a Nodulin-like protein (frv2_125506) had considerable effects on the separation of Me vs. Pr in disease treatments (Figure 3). Some proteins like phosphoenolpyruvate carboxylase, allergenic isoflavone reductase-like protein, and cysteine synthase were highly involved in discrimination of M vs. NM besides Me vs. Pr upon diseased treatments ( Figure 3 and Table S5). An uncharacterized protein (frv2_55445) belonged to amino acid metabolism (degradation, arginine) had a high impact on the separation between both cultivars (Me vs. Pr) and M vs. NM in diseased treatments (Figure 3 and Table S5). Two proteins (small nuclear ribonucleoprotein and mediator of RNA polymerase II transcription subunit 36a-like) mapped in RNA category were recognized with high influences on the separation of M vs. NM among diseased treatments ( Figure 3 and Table S5).
Cofactor-independent phosphoglycerate mutase (glycolysis) and protein disulfide isomerase-like (redox) had the highest impact on the separation of Me vs. Pr under diseased conditions ( Figure 3 and Table S5). Furthermore, an uncharacterized protein (frv2_75059) from protein synthesis class as well as flavonoid glucosyltransferase (secondary metabolism) indicated a strong effect on discrimination of M vs. NM upon diseased treatments (Figure 3 and Table S5). Among disease treatments, two metabolites including 2'hydroxydihydrodaidzein (isoflavonoid) and mesquitol (flavonoid) showed a remarkable impact on separating M vs. NM (Figure 3 and Table S5).
The data were then further reduced to focus on the most relevant protein groups potentially involved in Induced Systemic Resistance (ISR). For this, proteins that showed increased levels in M vs. NM of healthy plants overlapping with significantly increased levels upon pathogen infection of NM plants were selected for MapMan visualization (Figure 4). Names of the protein groups which were visualized, can be drawn from supplemental Table S6 as they were in the same order. Figure  4A schematically shows the priming effect where protein levels of AMF symbiotic and healthy plants were induced compared to NM plants. During pathogen attack these protein levels also were increased in NM plants while they were not changed in infected M plants as their levels were already induced (primed) ( Figure  4). Altogether, cv. Protecta shows a strong overlap of induced proteins ( Figure 4B and Table S6) between healthy M treated and stressed non-symbiotic (PrI) plants while only one protein was also slightly further increased upon infection in M plants (MI). In contrast, only a few proteins (four proteins) were induced in healthy M plants of cv. Messire and overlapping with a stress response of non-symbiotic plants ( Figure 4C and Table S6). The strongest responses of the seed metabolism of cv. Protecta (in terms of numbers) were proteins related to primary metabolism such as gluconeogenesis and amino acid synthesis and secondary metabolism-related proteins, including lipid and hormone (jasmonate) regulating proteins. In terms of fold FIGURE 3 | ICA of integrated seed metabolomics-proteomics data including remarkably changed (Kruskal-Wallis; ANOVA, Tukey HSD test; p < 0.05 and ≥ two-fold change) seed metabolites and proteins in M vs. NM upon diseased conditions of both cultivars. Numeric data were z-transformed. The loadings graphs illustrate the components with the top ten of the highest loadings (> 0.01 and < −0.01) on IC2 and IC3. Me, cv. Messire; Pr, cv. Protecta; M, mycorrhizal; I, infected (diseased). changes, two proteins of yet unknown function revealed the highest induction (>100 fold) upon AMF symbiosis of healthy and after infection of NM Protecta plants ( Table S6). The four proteins possibly involved in a priming effect in cv. Messire were a major intrinsic protein, a Polyvinylalcohol dehydrogenase-like protein (signaling), a Cytochrome b5-like protein, and a storage protein (legumin) ( Table S6).

AMF Efficiently Colonize P. sativum Roots and Influence Root Biomass
The successful colonization of plants by AMF is influenced by temperature, nutritional properties of soil, type of mycorrhizal symbiont, plant species, and light conditions (Smith and Smith, 1996;Johnson et al., 1997, as cited in Jansa et al., 2008. In this study, the high average AMF colonization (70%) of pea roots indicated a stable symbiosis between AMF species and the root system of P. sativum cultivars similar to our previous observations . Moreover, the considerable enhancement of root biomass in mycorrhiza symbiotic (M) plants compared with non-mycorrhizal (NM) plants in both cultivars and under disease stress prove the high impact of AMF on the root system as a plant growth-promoting symbiont in line with reported observations (Kavitha and Nelson, 2014;Oruru et al., 2017;Wang et al., 2018). Here, the remarkable increase of root biomass in M vs. NM plants besides significant positive correlation among root biomass (FW and DW) as well as between total mycorrhizal and plant root length (TMRL and TPRL) confirm the mutual connection between mycorrhizal root colonization and architecture of root system that was found before by Sinclair et al. (2014) in colonized roots of strawberry. However, to reach optimal AMF promoted crop productivity, the choice of AMF species that match with the various crops is determinative (Van Geel et al., 2016). Chen et al. (2017) have concluded that the application of AMF in a combined form of several mycorrhizal species is more effective  Table S6). LFQ, Label-Free Quantification. compared with single inoculation on cucumber. This is important, as we can say that the AMF mixture, we used in this study, is effective even if we cannot distinguish between the putative differences in efficiency of the different species. Hence, although our approach is more close to nature, a differential analysis of the various AMF species remains a future task. Nevertheless, the mutual impact of AMF on root growth and vice versa, as indicated previously (Atkinson et al., 1994, as cited in Sinclair et al., 2014, corroborates our results for both pea cultivars. Further studies on AMF species and inoculum composition may further increase growth-promoting efficiency. Interestingly, the lower TMRL and TPRL in cv. Messire than cv. Protecta related not only to different root structures but also indicate the dissimilar response of the two cultivars to mycorrhizal root colonization by the same AMF species. Similarly, we previously observed an individual cultivar rootmicrosymbiont colonization pattern between cv. Messire and Protecta when inoculated with Rhizobium (Turetschek et al., 2017). Different AMF root colonization and root architecture were also observed in olive cultivars (Chatzistathis et al., 2013). Oruru et al., 2017 not only found the dissimilarities in AMF root colonization between modern and wild cultivars of cowpea, they also reviewed the differences among cultivars of tomato (Steinkellner et al., 2012), wheat (Tarawaya, 2003), maize (Njeru et al., 2013), and barley (Zhu et al., 2003). However, not much is known about the reason for and consequences of these root architectural differences in relation to symbiotic efficiency and hence, further root phenotyping (incl. molecular) studies will be needed. Finally, it is noteworthy that no significant negative effect of the pathogen on AMF colonization compared to healthy plants could be found regardless of cultivars. In contrast, our prior study (Ranjbar Sistani et al., 2017) showed significant reduction of Rhizobium density isolated from nodulated roots caused by D. pinodes. These results indicate that AMF root colonization is more stable under pathogen attack compared to Rhizobium symbiosis.

Positive Effects of AMF Symbiosis on Growth, Yield, and Seed Physical Properties
In general, our results show that AMF as belowground symbionts also promote the growth of the above-ground parts of pea plants. Like in our previous study (Ranjbar Sistani et al., 2017), which focused on effects of bacterial (Rhizobium) symbiosis on pea seeds, here, we found that AMF had a significant impact on seed yield, seed (FW and DW) per seed, TSW (FW and DW), seed (FW) per plant, pod weight, and pod size. In earlier studies, the notable influence of AMF root colonization on seed yield enhancement (Cely et al., 2016;Adeyemi et al., 2017), seed weight per plant (Abdel-Fattah et al., 2016), pod weight (Abdel-Fattah et al., 2016;Adeyemi et al., 2017), and hundredseeds weight (Kavitha and Nelson, 2014;Abdel-Fattah et al., 2016) has been observed. Specific proteins such as glutamine synthetase (N-metabolism), enolase (glycolysis), the protein disulfide isomerase (redox), malate dehydrogenase (TCA cycle), peptidyl-prolyl cis-trans isomerase (cell) were remarkably enhanced in seed proteome of M vs. NM pea plants. These proteins have been recognized previously as essential proteins during Brassica campestris seed development (Li et al., 2012). The data indicate that AMF influence seed growth and development of P. sativum cultivars. Interestingly, Indole-3-acetamide (IAM), crucial for the synthesis of the auxin indole-3-acetic acid (IAA) (Mano et al., 2010), was reduced in our M treated plants. IAA is included in AMF symbiosis relation (Ludwig-Müller and Güther, 2007) and plant growth-regulation (Kögl and Kostermans, 1934;Went and Thimann, 1937). The role of AMF in promoting the synthesis of auxin and consequently enhancement of root-hair growth in Poncirus trifoliate was reported (Liu et al., 2018). The depletion of IAM upon AMF treatment might indicate that it was required for the formation of IAA, not detected in our study. Nevertheless, enhanced IAA production of seeds is not only of advantage for seed germination and development but also for microsymbiont attraction and interaction (Fu et al., 2015).
Strong genotype related impacts of AMF on pod weight in okra (Nwangburuka et al., 2012) have been reported earlier.
Proteome analysis also revealed that secondary metabolism, general stress, and redox response and also proteins of the primary metabolism of amino acid synthesis and glycolysis are more enhanced in cv. Protecta compared to cv. Messire. Interestingly, most of the AMF responsive proteins of those metabolic pathways are also involved in pathogen defense.

Belowground Mycorrhizal Symbiont Controls Seed Infection of the Aboveground Pathogen
The majority of studies associated with biocontrol and impact of mycorrhiza against plant pathogens are restricted to soil-borne and root fungal phytopathogens (Vierheilig et al., 2008). Here, the pathogen factor had significantly negative effects on aboveground parameters. The intensity of several proteins mapped for many functional classes including C1-metabolism, cell, development, DNA, metal handling, nucleotide metabolism, protein post-translational modification, redox, biodegradation of xenobiotics, signaling, TCA cycle, and transport were considerably reduced under stress caused by the pathogen. These negative effects were much more pronounced in cv. Messire compared to Protecta. These findings show the broad metabolic dampening effects of D. pinodes on P. sativum seeds, particularly for cv. Messire especially in the absence of AMF. Protecta, in contrast, shows a stronger increase in several proteins of the signaling, stress response, development, amino acid, glycolysis, and of the secondary metabolism, which might be related to the enhanced protection of cv. Protecta against the pathogen. Here, several of our previously identified metabolites and proteins involved in general pathogen resistance can be confirmed (Ranjbar Sistani et al., 2017).
For instance, we previously found proteins of the Late Embryo Abundant (LEA)-related family among the most significantly involved in pathogen resistance response (Ranjbar Sistani et al., 2017). Several groups of those (frv2_129189; frv2_25524, frv2_11113; frv2_54342, frv2_79176, P28639, frv2_128555) again showed a significant increase upon pathogen attack, however, independent of M treatment. LEA proteins are mainly known to be induced during seed development and in resistance to dehydration (Battaglia and Covarrubias, 2013). Our data confirm the findings of our previous work that LEA proteins are important for enhanced seed pathogen resistance and are not only involved in drought stress tolerance (Ranjbar Sistani et al., 2017).
Likewise, vicilin storage proteins (frv2_74601 and frv2_80935) were again found with a stronger increase under infection in cv. Protecta, however, slightly but not significantly stronger in M vs. NM treated plants. Nevertheless, this effect was more pronounced in rhizobia symbiotic plants (Ranjbar Sistani et al., 2017).
Also, the phytohormone abscisic acid (ABA) is a known component in plant environmental stress reactions (Luo et al., 2009;Song et al., 2014;as cited in Porcel et al., 2012;Adolfsson et al., 2017). A higher ABA content has been detected in leaves of Medicago truncatula colonized by AMF (Adolfsson et al., 2017). In Arabidopsis seeds, it has been reported that salicylic acid application induced ABA signaling leading to increased synthesis of ABA-regulated proteins, such as LEA proteins, dehydrins, and heat shock proteins, (Rajjou et al., 2006). Although here, ABA was notably enhanced in seeds upon pathogen infection, this accumulation was not enhanced in AMF symbiotic plants. Thus, our data support a role of ABA during pathogen defense and enhanced LEA protein levels independent on AMF symbiosis.
Interestingly, we found several plastidic proteins in pea seeds strongly depleted in both cultivars upon pathogen defense, independent of AMF symbiosis. Plastids are important components for embryos development in P. sativum (Smith et al., 1990). The role of chloroplast in the seeds biotic stress response, however, has not been described before. We suggest that the huge breakdown of plastidic proteins, we observed upon infection, might serve as a source for the enhanced secondary metabolite synthesis. The findings are in line with Fondevilla et al. (2011) who found that genes assigned to the primary metabolism were mostly down-regulated in accession P665, a resistant genotype, compared to Messire as susceptible pea cultivar against D. pinodes infection. Altogether, a common pathogen response of both cultivars is evident through enhanced levels of proteins involved in (iso) flavonoid production, of which the most significant and common compound identified was Sativan. Unfortunately, not much is known about the function of this isoflavone. Our data, however, support that it is involved in pathogen defense, especially of NM treated plants.
We found a significant reduction in seed infection levels, as well as lower levels of leaf lesion areas, when comparing M vs. NM treated plants of both cultivars. Remarkably, with an almost 2-fold reduction of seed infection levels the bio-control potential of mycorrhizal symbiosis against D. pinodes was very similar to the potential of rhizobial symbiosis demonstrated in our previous study (Ranjbar Sistani et al., 2017). This is particularly interesting because our earlier studies on AMF symbiosis showed no reduction in leaflets disease severity upon D. pinodes attack Turetschek et al., 2017). Nevertheless, it has been reported that AMF protect plants against pathogens through the dampening of phytopathogenic effects (Borowicz, 2001). Also, AMF inhibiting impact on disease caused by Xanthomonas translucens on leaves of wheat was observed (Fiorilli et al., 2018). Altogether, the significant decrease of seed infection besides smaller lesion areas in AMF treated plants highlights the bio-control potential of AMF as a belowground microsymbiont reducing disease severity caused by an aboveground pathogen, in line with the previously reviewed findings by Jung et al. (2012), reporting the role of AMF in the reduction of infections caused by Botrytis cinerea (Møller et al., 2009;Pozo et al., 2010), Alternaria solani (Fritz et al., 2006;de la Noval et al., 2007), and Magnaporthe grisea (Campos-Soriano et al., 2012). Additionally, the high potential of AMF to control B. cinerea in tomato has been identified (Fiorilli et al., 2011). Nevertheless, compared to Rhizobium inoculation , AMF seems to have less significant effects on pea cultivar growth promotion and leaf protection but a similar impact on seed pathogen resistance. This study together with our earlier work (Ranjbar Sistani et al., 2017) revealed that both mycorrhizal and rhizobial below-ground symbionts are dominant against above-ground pathogen invasion. Besides, the below-ground symbiotic activity by AMF or Rhizobium reduced the systematic infection caused by the above-ground pathogen. Both studies indicate that the efficiency of these fungal and bacterial microsymbionts could be influenced by cultivar type and pathogen infection, however, their general bio-control function and infection inhibition is cultivar independent.
Levels of Vignatic acid A were increased in seeds of M treatments of cv. Messire upon pathogen infection. Vignatic acid A is known as a preventing compound with a promoting impact on resistance against bruchid beetle pest (Sugawara et al., 1996).
This cultivar specific effect becomes even more complex when comparing metabolic pathways that are induced by AMF symbiosis of healthy plants and are potentially involved in pathogen response. This specific mechanism (priming) of the plant metabolism by microsymbionts such as AMF, before infection by biotrophic and hemibiotrophic pathogens, is known as ISR (Cameron et al., 2013). Primed plants have already developed a certain metabolic defense apparatus, which reduces or even prevents pathogen attack. For instance, phytohormone JA plays a role in inducing systemic resistance against stresses (Luo et al., 2009;Song et al., 2014as cited in Staudinger et al., 2016Adolfsson et al., 2017;Bernardo et al., 2017) and also during plant development (Bernardo et al., 2017). Also, JA has a high impact on mycorrhizal colonization especially in higher plants (Gutjahr et al., 2015, as cited in Bernardo et al., 2017. According to our results, a protein related to the synthesis of JA (lipoxygenase frv2_95615 and frv2_84258) was highly increased in M plants (Table S4). Furthermore, data indicate that the JA-responsive protein (frv2_84258) is involved in ISR (Figure 4 and Table S6). Additionally, we identified N-jasmonoyl isoleucine, and it was significantly enhanced in cv. Protecta, similar to the levels of JA in our previous study (Ranjbar Sistani et al., 2017), but this time strongly induced by AMF instead of Rhizobium symbiosis. The JA has a role in regulating of pea reaction against D. pinodes (Fondevilla et al., 2011). These findings support that cv. Protecta is not only a more resistant genotype but is also more receptive to ISR mechanisms by AMF. Hence, our data support that the effect of AMF treatment on the seed metabolism and in line with our previous findings of Rhizobium treatments is not due to a difference in P-regime but rather depending on the cultivar and due to a global change in metabolic homeostasis of the whole plant induced by microsymbiosis. Although pathogenesisrelated (PR) proteins are known to be commonly involved in pathogen defense (van Loon et al., 2006), we found only a few of them involved in pea leaf and seed protection Ranjbar Sistani et al., 2017). Here again, we found possible isoforms (frv2_110573; frv2_74661; frv2_41448; frv2_88778) induced upon pathogen infection in cv. Protecta seeds only. Interestingly, they were already induced by AMF symbiosis and further accumulated upon stress. Besides PR proteins, the detection of proteins involved in cell wall fortification, such as proteins of the monolignol pathway, seemed limited. Nevertheless, we found a UDP-glucosyltransferase as well as UDP-glucose 4-epimerase and UDP-glucose 6-dehydrogenase possibly involved in cell wall synthesis. These three enzymes were induced upon pathogen attack of cv. Protecta, but only UDP-glucosyltransferase (frv2_103598) was also enhanced in seeds of AMF symbiotic plants (cv. Protecta). Furthermore, Cyanidin 3-sophoroside 5-glucoside seemed involved in better seed pathogen resistance upon AMF symbiosis of cv. Protecta, being only induced upon infection of M treated plants. These findings suggest that AMF does not only prime.
The 2,3-Dihydroxy-2,4-cyclopentadien-1-one was also accumulated upon pathogen infection in both cultivars with a more significant effect on cv. Protecta. Very striking were also the high levels of several proteins of amino acid metabolism and glycolysis in cv. Protecta. Protein levels showed a strong increase under pathogen infection of non-mycorrhizal cv. Protecta, but some were also induced in healthy M treated plants and did not much further increase in Protecta seeds of M plants after pathogen attack, indicative for a primed primary metabolism as previously described by Schwachtje et al. (2018). By the production of pyruvate, the entry metabolite of the TCA cycle, glycolysis plays a role in energy supply and thus the production of defense metabolites (Rojas et al., 2014). Hence, the data suggest that the cultivar related difference in the primary metabolism seems another reason for the better performance of cv. Protecta.
Nevertheless, some defense mechanisms were also found to be more pronounced in cv. Messire. For instance, flavonoid Quercetin 3-(6''-acetylgalactoside)-7-rhamnoside strongly accumulated upon pathogen infection, especially in cv. Messire and most strongly in seeds of M treated plants. Similarly, Pisatin as specific phytoalexin of P. sativum has been induced against pathogen attack (De Wit-Elshove, 1969). Its synthesis pathway has previously been observed to be induced in the leaf proteome of Rhizobium inoculated pea plants under infection caused by D. pinodes . In contrast, pisatin involvement in rhizobial induced seed protection could not be detected in our previous study (Ranjbar Sistani et al., 2017). Interestingly, in the present study, a considerable AMF induced enhancement of proteins and metabolites of the secondary metabolism such as proteins of the flavonoid biosynthesis, especially the accumulation of isoflavone reductase and chalcone isomerase, along with increased levels of flavonoids and isoflavonoids, mainly in seeds of cv. Messire under disease stress was found. This finding supports that pisatin biosynthesis is induced in seeds of M plants upon pathogen infection.
Although, AMF showed cultivar specific effects, particularly for cv. Protecta, our results also indicated a general enhancement of pathogen resistance through AMF symbiosis ( Figure 5). For example, G-proteins, known as signal transduction regulators (Smith et al., 1998;van Heusden, 2009;Smith et al., 2011;Parua andYoung, 2014, as cited in Sun et al., 2018) have a role in plant responses to various stress conditions (Roberts et al., 2002;Lozano-Duran and Robatzek, 2015;Li et al., 2016, as cited in Sun et al., 2018. Mycorrhizal induced accumulation of Gproteins in wheat roots against drought stress was also reported (Bernardo et al., 2017). In our study, the group of proteins (frv2_86036; frv2_118099; Q9C5W6; frv2_94136) was particularly increased in seeds of AMF-associated plants under pathogen disease, indicating a possible role in ISR signaling.
In general, we found a similar molecular defense pattern induced by either AMF like Rhizobium against D. pinodes, described in our prior study (Ranjbar Sistani et al., 2017). This is not surprising, as the host builds fortification and munitions. Hence, Flavonoids and Soysapogenol C were important components in secondary metabolism adjustment in both Rhizobium-and AMF-induced responses. Interestingly, however, L-DOPA 3'-glucoside, discussed above, was significantly induced by AMF and in response to stress and not found responsive to Rhizobium-symbiosis. Also, the main proteins involved in signal transduction, mentioned above, were differentially induced. While, AMF enhanced levels of G-proteins and phosphatase 2C seems Rhizobium specific (Ranjbar Sistani et al., 2017). The analysis of secondary metabolites is still in its infancy and identification remains a bottleneck. Thus, major symbiont and pathogen involved differences and further key secondary metabolites regulation pathogen defense, remain to be identified in the future. The difference might also lie in response rate and/or induction level that might be different not only between cultivars but also depending on symbiosis, terms that need further evaluation.
Taken together, among common stress responses, several proteins and metabolites have been found AMF and cultivar specifically regulated, demonstrating different defense strategies. While genotype-specific resistance pattern is a common base in plant breeding, not much is known about genotype-specific microsymbiont induced priming effects.

CONCLUSIONS
In this study, we demonstrated that AMF is not only promoting growth and yield of pea plants but also protects seed yield upon pathogen attack. The seed metabolic response of two cultivars with varying susceptibility to the pathogen D. pinodes was analyzed. Besides, with the growth-promoting effect of AMF microsymbiont, we demonstrate its positive impacts on the resistance of the plants against the pathogen. We found that the genotypic effects are strong for both the general and the AMF induced pathogen response.
Altogether, cv. Protecta, the less susceptible genotype, accumulates more proteins of the secondary metabolism involved in lipid, flavonoid, and phytohormone production as well as of the primary amino acid and glycolysis metabolism. Besides, cv. Protecta showed a stronger induction (priming) of seed metabolic pathways upon AMF symbiosis, which led to a dampened and less severe but more effective response to pathogen attack. Nevertheless, pathogen infection was also reduced in cv. Messire when interacting with AMF. However, infection levels were similar compared to cv. Protecta without AMF symbiosis and thus AMF impact less effective. Thus, the data demonstrate AMF and genotype-specific pathogen defense strategies. This study demonstrates that AMF as below-ground microsymbiont not only promote the growth above-below ground parts of pea plants but enhance the resistance of P. sativum plants and hence protects seed quantity. Also, AMF have a notable influence on the seed metabolome and proteome, influencing seed metabolism and thus nutritional and medicinal pea seeds quality.
Although, mycorrhizal and rhizobial symbionts in singular applications could be prescribed for improvement of pea seed protection, little is known about impacts of co-inoculation, often occurring in nature, and therefore investigations are demanded.
We believe that sustainable development of food resources for food poverty elimination can be strongly supported through the enhancement of productivity and quality of crops by using microbial symbiotic potentials. Focusing on plant species FIGURE 5 | Schematic overview of mycorrizal induced response of pea seed (grown in pots) against D. pinodes infection. Seed metabolites and proteins significantly (Tukey HSD, p < 0.05; n = 3; Supplementary Table S3B and S4) enhanced with the highest ICA loadings (> 0.01 and < −0.01; Supplementary  Table S5) are shown. MI, Mycorrhizae inoculated, and pathogen-infected. supplying food and microsymbiont interaction studies on improving production-and seed quality-related strategies under stress conditions besides promoting current ecological properties should be a priority in plant and agricultural studies. We recommend more studies related to efficiency assessment of soil microorganisms on food plants productivity upon abiotic and biotic stresses under field conditions with considering the environmental viabilities.

DATA AVAILABILITY STATEMENT
The datasets generated for this study can be found in the PRIDE PXD006617.

AUTHOR CONTRIBUTIONS
Experiments, including phenotyping, proteomics, and secondary metabolites analyses were carried out by NRS. Data mining was performed by NRS and SW. NRS and SW wrote the manuscript. SW, GD, and H-PK conceived the work. The manuscript was revised by GD and H-PK.

FUNDING
We appreciated the support from COST Action FA1306.

ACKNOWLEDGMENTS
The authors thank the gardeners of the division of Molecular Systems Biology (Mosys, Faculty of Life Sciences, University of Vienna), Thomas Joch, and Andreas Schröfl. Thanks also to Sonja Tischler for running and maintaining the mass spectrometry facility of Mosys division. We thank Core Facility Cell Imaging and Ultrastructure Research (CIUS, Faculty of Life Sciences, University of Vienna) for supplying microscopy equipment, Dr. Irene Lichtscheidl-Schultz, and Dr. Ingeborg Lang for their instructions about microscopy equipment. Also, we warmly thank Dr. Wenchao Zhang for his supportive and effectual instructions about MET-COFEA and MET-XAlign tools.