# Metabolic Fingerprint of PS3-Induced Resistance of Grapevine Leaves against Plasmopara viticola Revealed Differences in Elicitor-Triggered Defenses

Marielle Adrian, Marianna Lucio, Chloé Roullier-Gall, Marie-Claire Héloir, Sophie Trouvelot, Xavier Daire, Basem Kanawati, Christelle Lemaître-Guillier, Benoît Poinssot, Régis Gougeon, et al.

► To cite this version:

Marielle Adrian, Marianna Lucio, Chloé Roullier-Gall, Marie-Claire Héloir, Sophie Trouvelot, et al. "Metabolic Fingerprint of PS3-Induced Resistance of Grapevine Leaves against Plasmopara viticola Revealed Differences in Elicitor-Triggered Defenses." Frontiers in Plant Science, 2017, 8, pp.101. 10.3389/fpls.2017.00101 hal-01534517

HAL Id: hal-01534517

https://u-bourgogne.hal.science/hal-01534517

Submitted on 25 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Check for
updates

# Metabolic Fingerprint of PS3-Induced Resistance of Grapevine Leaves against *Plasmopara viticola* Revealed Differences in Elicitor-Triggered Defenses

Marielle Adrian<fnref n="1" />*, Marianna Lucio², Chloé Roullier-Gall², Marie-Claire Héloir<fnref n="1" />, Sophie Trouvelot<fnref n="1" />, Xavier Daire<fnref n="1" />, Basem Kanawati², Christelle Lemaître-Guillier<fnref n="1" />, Benoît Poinssot<fnref n="1" />, Régis Gougeon<fnref n="3" /> and Philippe Schmitt-Kopplin²,<fnref n="4" />

**OPEN ACCESS**

**Edited by:**
Fumiya Kurosaki,
University of Toyama, Japan

**Reviewed by:**
Daryl Dryden Rowan,
Plant and Food Research,
New Zealand
Luigi Bavaresco,
Università Cattolica del Sacro Cuore,
Italy

*Correspondence:
Marielle Adrian
marielle.adrian@u-bourgogne.fr

**Specialty section:**
This article was submitted to Plant Metabolism and Chemodiversity, a section of the journal Frontiers in Plant Science

Received: 23 November 2016

Accepted: 18 January 2017

Published: 14 February 2017

**Citation:**
Adrian M, Lucio M, Roullier-Gall C, Héloir M-C, Trouvelot S, Daire X, Kanawati B, Lemaître-Guillier C, Poinssot B, Gougeon R and Schmitt-Kopplin P (2017) *Metabolic Fingerprint of PS3-Induced Resistance of Grapevine Leaves against Plasmopara viticola Revealed Differences in Elicitor-Triggered Defenses*. Front. Plant Sci. 8:101. doi: 10.3389/fpls.2017.00101

Induction of plant resistance against pathogens by defense elicitors constitutes an attractive strategy to reduce the use of fungicides in crop protection. However, all elicitors do not systematically confer protection against pathogens. Elicitor-induced resistance (IR) thus merits to be further characterized in order to understand what makes an elicitor efficient. In this study, the oligosaccharidic defense elicitors H13 and PS3, respectively, ineffective and effective to trigger resistance of grapevine leaves against downy mildew, were used to compare their effect on the global leaf metabolism. Ultra high resolution mass spectrometry (FT-ICR-MS) analysis allowed us to obtain and compare the specific metabolic fingerprint induced by each elicitor and to characterize the associated metabolic pathways. Moreover, erythritol phosphate was identified as a putative marker of elicitor-IR.

**Keywords:** induced resistance, elicitor, grapevine, downy mildew, metabolomics

## INTRODUCTION

Wine industry is of high economic, social and cultural value worldwide (Bisson et al., 2002). However, viticulture has to adapt in order to face major changes such as climate evolution and environmental constraints. A true challenge is currently the evolution towards production systems combining sustainability, economical viability, and more eco-friendly practices.

The majority of the grown grapevine cultivars (*Vitis vinifera* L. cvs) are susceptible to cryptogamic diseases such as downy mildew or powdery mildew. Repeated pesticide applications are generally used to ensure the expected yield and harvest quality. However, some fungicides are harmful for the environment and human health, and contribute to the selection of resistant strains of pathogens (Casida, 2009; Nanni et al., 2016). Besides organic farming or the development of resistant varieties by grape breeding, one strategy that is likely to allow the reduction of the use of fungicides is the induction of plant resistance to pathogens by defense elicitors (Walters et al., 2013).

<footnote n="1">¹ Agroécologie, AgroSup Dijon, CNRS, INRA, Université Bourgogne Franche-Comté, Dijon, France, ² Analytical BioGeoChemistry, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany,</footnote>

<footnote n="3">³ UMR PAM Université de Bourgogne/AgroSupDijon, Institut Universitaire de la Vigne et du Vin, Jules Guyot, Dijon, France,</footnote>

<footnote n="4">⁴ Chair of Analytical Food Chemistry, Technische Universität München, Freising-Weihenstephan, Germany</footnote>
Elicitors are compounds of natural origin (derived from plants or microbes) or synthetically generated (Walters et al., 2013). Their perception by plant cells activates a cascade of signaling events, including H₂O₂ production, that induces defense gene expression, leading to defense responses, i.e., the synthesis of Pathogenesis-related (PR-) proteins, the antimicrobial compounds phytoalexins, and cell-wall reinforcement (Garcia-Brugger et al., 2006). Phytohormones such as jasmonic acid (JA), salicylic acid, and ethylene also contribute to signaling (Shigenaga and Argueso, 2016). Defense activation thus triggers a global cell reprogramming towards defense responses (Wiesel et al., 2014) and therefore generates a cost that the plant has to fuel with primary metabolism (Bolton, 2009; Rojas et al., 2014), sometimes at the expense of growth (Huot et al., 2014).

Several compounds were reported as elicitors of grapevine defenses (Adrian et al., 2012; Delaunois et al., 2013). Among them, the β-glucans laminarin and its sulfated derivative PS3 (Ménard et al., 2004) have been studied in details as they induce resistance against *Plasmopara viticola* (Pv), the obligate biotrophic oomycete responsible of downy mildew (Aziz et al., 2003; Trouvelot et al., 2008; Gauthier et al., 2014). The perception and subsequent biological activity of oligosaccharides is highly dependent on their degree of polymerization (DP) and molecular pattern, i.e., acetylation, methylation, and/or sulfation (Trouvelot et al., 2014; Paris et al., 2015). As example, the short laminarin H13 (or Lam13, DP of 13) triggers the production of H₂O₂ in grapevine leaves without subsequent induced resistance (IR) against Pv (Allègre et al., 2009). Conversely, PS3 (DP of 25, degree of sulfation of 2.4) primes H₂O₂ production and induces protection against downy mildew (Trouvelot et al., 2008; Allègre et al., 2009; Gauthier et al., 2014). Such observations also show that plant defense activation is necessary, but not sufficient to induce resistance against pathogens. Numerous studies have reported new elicitors and mainly focused on the analysis of induced defense events and their regulation, and their efficiency to protect plants against diseases<fnref n="1" /> (Wiesel et al., 2014). However, to our knowledge, no comparative study using elicitors effective/ineffective to trigger IR was previously performed. What makes the difference between defense activation and IR thus remains unclear and merits to be investigated.

Global approaches have become powerful tools for diverse biological features. Metagenomic, transcriptomic, proteomic and metabolomic methods, now called “omics,” indeed provide useful large-scale data to understand biological systems (Buescher and Driggers, 2016). Metabolomic is suitable for different scales, from cell to organs. Among the different analytical systems available, the high sensitivity of mass spectrometry makes it a tool of choice for most of the plant metabolism analyses (Jorge et al., 2016). In plants, metabolomics approaches have been used for various studies (Jorge et al., 2016) including development, such as leaf senescence (Kim et al., 2016), and response to biotic or abiotic stresses, such as salt resistance (Janz et al., 2010). In grapevine, metabolomic was mainly used for variety characterization (Flamini et al., 2013, 2014) and to study

berry evolution (Suzuki et al., 2015), impact of environmental conditions on berry, must and wines (Roullier-Gall et al., 2014a; Matus, 2016), and response to drought stress (Hochberg et al., 2013; Griesser et al., 2015; Shiratake and Suzuki, 2016) or pathogens (Batovska et al., 2009; Becker et al., 2013; Agudelo-Romero et al., 2015).

The objective of this study was to better characterize elicitor-IR of the grapevine leaves against downy mildew. The two β-glucan elicitors H13 and PS3, respectively, ineffective and highly effective to induce resistance against downy mildew (i.e., unable and able to limit the pathogen development and confer a protection against the disease, respectively), were chosen. Ultra high resolution mass spectrometry (FT-ICR-MS) was used to perform a global metabolomic analysis in order to (i) obtain and compare the metabolic fingerprint induced by each elicitor, (ii) characterize the metabolic fingerprint and pathways associated to elicitor-IR, and (iii) find putative markers of elicitor-IR.

# MATERIALS AND METHODS

## Plant Material

The grapevine cultivar *V. vinifera* L. cv. Marselan (Cabernet sauvignon × Grenache), susceptible to *Plasmopara viticola*, was used for this study. Plants were obtained from herbaceous cuttings placed in individual pots (10 cm × 10 cm × 7 cm) containing a mixture of blond peat and perlite (3:2, v/v). They were grown in a glasshouse at a temperature of 24 and 18°C (day and night, respectively) with a photoperiod of 16 h of light and at a relative humidity (RH) of 70 ± 10 % until they developed six leaves. Plants were daily sub-irrigated with a fertilizer solution (NPK 10-10-10, Plantin, France).

## Elicitor Treatment

Two elicitors provided by Goëmar Laboratories were used for this study: H13 (or Lam13, here noted “H”), a short laminarin with a DP of 13 (Allègre et al., 2009) and the sulfated laminarin PS3 (here noted “P”; Ménard et al., 2004; Trouvelot et al., 2008).

Leaf disks (14 mm diameter) were punched from leaves (the third below the apex) of grapevine plants grown in greenhouses as described above. They were washed in distilled water and placed at the surface (upper side above) of distilled H₂O (W, as control), H or P solutions in water (2.5 g.l⁻<fnref n="1" />) in sterile Petri dishes (90 mm diameter) for 48 h (Figure 1). Afterwards, they were removed, rinsed in distilled water and deposited on a moist filter paper in a plastic box at room temperature before mock- or Pv inoculation. All disks were deposited on the same filter paper in the same plastic box to allow further tight comparison between treatments. Thirty disks were prepared per condition, from the third leaves of three plants (n = 30).

## *Plasmopara viticola* Infection

Sporangia were collected from sporulating leaves of Marselan herbaceous cuttings as previously described (Trouvelot et al., 2008) and their concentration was adjusted to 10⁴ sporangia per milliliter using an hemacytometer. At 48 hours after the

<footnote n="1">¹http://www.elicitra.org/index.phpvrs=english</footnote>
<figure><img src="image_3.png" /><figcaption>**FIGURE 1 | Outline of the experiment.** Disks were punched out from leaves of *Vitis vinifera* L. cv. Marselan herbaceous plantlets and treated by distilled H₂O (W, as control), H or P elicitor solutions (2.5 g.l⁻¹). At 48 hpt, they were rinsed in distilled water and deposited on a moist paper filter before mock- or *P. viticola* inoculation (suspension of 10⁴ sporangia per milliliter of water). Thirty disks were prepared per condition, from the third leaves (under the apex) of three plants (n = 30). Disks were collected at 48, 72, and 96 hpt for analysis. Other ones were let until sporulation (144 hpt).</figcaption></figure>

beginning of the treatment (48 hpt), leaf disks were inoculated by a droplet (70 µl) of either the Pv sporangia suspension or water (mock) and incubated at room temperature (Figure 1).

## Sample Preparation

At 48 hpt (i.e., just before Pv- or mock- inoculation), 72 and 96 hpt, the disks were reduced to 12 mm in diameter using a cork borer (to eliminate the outlying area where defense responses to injury occured) and grind in liquid nitrogen. Samples were then stored at -80°C until extraction and analysis. Three disks were frozen separately per condition. The remaining disks were kept until 7 days post inoculation (dpi) in the same conditions in order to quantify the level of sporulation as previously described (Kim Khoook et al., 2013; Figure 1).

## FT-ICR-MS Analysis

Prior to analysis, 15 mg of each sample were mixed with 1 ml methanol (LC-MS grade, Fluka Analytical; Sigma-Aldrich, St. Louis, USA) and sonicated for 30 min. After centrifugation (25,000 g, 10 min, room temperature), the supernatant was collected and re-diluted in methanol (1/50 v/v). Ultra high-resolution mass spectra were acquired using an Ion Cyclotron Resonance Fourier Transform Mass Spectrometer (FT-ICR-MS) (solariX, BrukerDaltonics GmbH, Bremen, Germany) equipped with a 12 Tesla superconducting magnet (Magnex Scientific Inc., Yarnton, GB) and an APOLO II ESI source (Bruker Daltonics GmbH, Bremen, Germany) operated in the negative ionization

mode. Samples were introduced into the micro electrospray source at a flow rate of 120 µl.h⁻¹. The MS was externally calibrated on clusters of arginine (10 mg.l⁻¹ in methanol). Spectra were acquired with a time domain of four mega words over a mass range of 100 to 1,000 Da and 300 scans were accumulated per sample.

## Data Analysis

Spectra were internally recalibrated using a list of ubiquitous fatty acids and recurrent wine compounds, allowing mass accuracies of 0.1 ppm (Gougeon et al., 2009). The *m/z* peaks with a signal-to-noise ratio (S/N) of 4 and higher were exported to peak lists. Exact masses were then subjected to Netcalc algorithm and in-house software tool to obtain chemical formulas (Tziotis et al., 2011) validated by setting sensible chemical constraints (n rule; O/C ratio ≤ 1; H/C ratio ≤ 2n+2; element counts: C ≤ 100, H ≤ 200, O ≤ 80, N ≤ 3, S ≤ 3, and P ≤ 1). In order to facilitate the interpretation of elemental formulas attributed to *m/z* values, formulas were next represented using two-dimensional van Krevelen (VK) diagrams. These diagrams display the hydrogen/carbon (H/C) versus oxygen/carbon (O/C) ratios of annotated elemental formulas and provide a commonly used qualitative description of the molecular complexity of biological samples (Tziotis et al., 2011; Roullier-Gall et al., 2014b). These plots provide additional information to function or biochemical pathways, through the representation of areas covering the various
metabolite classes of samples. As such, they allow the qualitative comparison between series of related samples, in terms of their chemical complexity. The m/z distribution according to their elemental compositions (CHO, CHOS, CHON, CHONS, CHOP, CHONP, CHONSP) was also represented. Mass annotations were obtained by KEGG, HMDB and LipidMaps databases queries. The metabolic pathways associated to selected masses were obtained by KEGG query with *Vitis vinifera* (vvi) as organism using the web server MassTRIX. Network analysis was performed from the complete sample set data as previously described (Roullier-Gall et al., 2014b). Nodes represent m/z values (metabolite candidates) and edges represent chemical reactions.

OPLS-DA (Orthogonal Projections to Latent Structures-Discriminant Analysis) was applied to datasets in order to discriminate between different sample groups. The overfitting was checked by the p-value calculated with the CV-ANOVA (Cross Validation ANOVA). From the classification models, the m/z values with the highest regression coefficient were extrapolated in different lists, as representative masses for the different clusters (called “top m/z”). These analysis were performed with SIMCA-P+13.0.3 (Umetrics, Umea, Sweden). Moreover for the analysis of the 96i and 72i datasets, it was shown a great improvement of the classification models applying before them the Relief algorithm. This chooses the features that can be most distinguished between classes (Package FSelector, Rstudio Inc., Version 0.99.896).

# RESULTS

## W48, H48, and P48 Samples Had Specific Metabolic Fingerprints

The FT-ICR-MS analysis of all the samples provided a total of 8,067 m/z. The corresponding chemical composition

histogram and VK diagram highlighted a wide chemical diversity, with a significant contribution of CHO elemental formulas (Supplementary Figure S1). A total of 5,283 m/z was present in the 48 hpt samples. Among them, 677, 831, and 1,103 ones were specific to W48, H48, and P48, respectively, whereas 1,553 ones were common to the three series (Figure 2A). The corresponding VK diagrams and chemical composition histograms obtained from W48, H48, and P48 did not show notable differences between these three sample sets (Supplementary Figure S2).

Orthogonal Projections to Latent Structures-Discriminant Analysis was next applied to discriminate W48, H48, and P48 samples (Figure 2B). W48 could be separated from P48 and H48 (p-value = 0.02), H48 from P48 and W48 (p-value = 0.05), and P48 from H48 and W48 (p-value = 0.06) (Table 1). Masses with the highest regression coefficient values in each sample group (specific and common masses with a variable importance in projection (VIP) value > 1; thereafter called “top m/z”) were next used for comparison. A total of 1,070 top m/z were isolated for W48, H48, and P48 (Table 2). Among them, 90, 34 and 81 were specific to W48, H48, and P48, respectively (Figure 2C), representing about 50, 6, and 26% of their total top m/z number. The discrimination between the three sample sets was further characterized using the chemical histogram and VK diagram built from the elemental formulas attributed to their respective top m/z list (Figure 3). W48 chemical histogram was highly different from H48 and P48 ones, with a higher number of CHOS elemental formulas and a lower number of CHO, CHOP, CHONS, CHON, CHONP, and CHONSP ones (Figures 3A–C). VK diagrams also showed noteworthy differences in the nature and the number of formulas assigned to masses between the three sample groups (Figures 3D–F), leading to specific fingerprints. W48 had a greater diversity in CHOS formulas. H48 presented a wide chemical diversity, with important CHO formulas. P48 was

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Discrimination of FT-ICR-MS data of the 48 hpt samples. Venn diagrams showing the distribution of the total (A) and top (C) m/z obtained by (–) FT-ICR-MS analysis of W48, H48, and P48 samples. Score plot (B) for the OPLS-DA analysis showing the discrimination of the W48, H48, and P48 samples. W48, H48, and P48 correspond to methanol extracts of grapevine leaf disks treated with H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g.l⁻¹) for 48 h. Top m/z correspond to m/z with the highest regression coefficient value (VIP ≥ 1).</figcaption></figure>
<table><caption>TABLE 1 | Significance of the sample discrimination performed by OPLS-DA analysis.</caption><thead><tr><th>Model</th><th>Type</th><th>CV-Anova</th><th>RY (cum)</th><th>Q (cum)</th></tr></thead><tbody><tr><td colspan="5"><strong>48</strong></td></tr><tr><td>H vs. (P, W)</td><td></td><td>0.05</td><td>0.94</td><td>0.62</td></tr><tr><td>P vs. (H, W)</td><td>OPLS-DA</td><td>0.06</td><td>0.82</td><td>0.6</td></tr><tr><td>W vs. (H, P)</td><td></td><td>0.02</td><td>0.87</td><td>0.69</td></tr><tr><td colspan="5"><strong>72i</strong></td></tr><tr><td>H vs. (P, W)</td><td></td><td>0.03</td><td>0.97</td><td>0.881</td></tr><tr><td>P vs. (H, W)</td><td>OPLS-DA</td><td>0.02</td><td>0.998</td><td>0.963</td></tr><tr><td>W vs. (H, P)</td><td></td><td>0.02</td><td>0.886</td><td>0.689</td></tr><tr><td colspan="5"><strong>96i</strong></td></tr><tr><td>H vs. (P,W)</td><td></td><td>0.04</td><td>0.977</td><td>0.866</td></tr><tr><td>P vs. (H, W)</td><td>OPLS-DA</td><td>0.02</td><td>1</td><td>0.991</td></tr><tr><td>W vs. (H, P)</td><td></td><td>0.04</td><td>0.964</td><td>0.865</td></tr></tbody></table>

W, H, and P correspond to methanol extracts of grapevine leaf disks treated with H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g.l⁻¹) for 48 h. R²Y and Q² are two values calculated after the cross validation of the model; they represent the measurements of the goodness of the fit and the predictive quality of the model.

<table><caption>TABLE 2 | Number of top m/z determined for each treatment and sampling time.</caption><thead><tr><th></th><th>W</th><th>H</th><th>P</th><th>TOTAL</th></tr></thead><tbody><tr><td>48</td><td>181</td><td>580</td><td>309</td><td>1070</td></tr><tr><td>72i</td><td>114</td><td>92</td><td>159</td><td>365</td></tr><tr><td>96i</td><td>36</td><td>79</td><td>146</td><td>261</td></tr></tbody></table>

W, H, and P correspond to methanol extracts of grapevine leaf disks treated with H₂O (W, as control), H (H13) or P (PS3) elicitor solutions (2.5 g.l⁻¹) for 48 h.

chemically less diversified than H48 and showed a CHOS specificity.

Query of KEGG database (Vitis vinifera organism) with the MassTRIX interface allowed the annotation of top m/z for W48, H48, and P48, respectively, and their assignment to metabolic pathways (Figure 4). Few pathways were associated to W48 and P48 (Figure 4). “Biosynthesis of secondary metabolites” and “Purine metabolism” pathways were attributed to P48 top m/z. Among P48 annotated compounds were AMP, ADP, ATP, UTP, tryptophan, palmitic acid, oleic acid, JA, abscisic acid, traumatic acid, and syringin (Figure 4, Supplementary Table S1). H48 showed the highest diversity of pathways, covering “Biosynthesis of secondary metabolites”, “Pentose and glucuronate interconversions”, “Ascorbate and aldarate metabolism”, “Starch and sucrose metabolism”, “Amino sugar and nucleotide sugar metabolism”, and “Flavonoid biosynthesis” (Figure 4). These pathways included compounds annotated as ribose; fructose 6P, fructose 1,6 diP, trehalose 6P, glucosamine 6P, glyceraldehyde 3P, gluconic acid, ascorbic acid, the oxidized and reduced forms of glutathion, quercetin 3-O-glucoside, polydatine and ε-viniferin (Supplementary Table S1). “Biosynthesis of secondary metabolites” was the unique pathway common to W48, H48, P48 top m/z (for a number of attribution exceeding 5) whereas “Purine metabolism” was specific to P48 (Figure 4).

## W, H, and P Fingerprints Evolved with *P. viticola* Infection and Time But Kept Specificity

At 48 hpt, leaf disks were mock- or *P. viticola* inoculated (thereafter noted “ni” and “i”, respectively) and collected at 72 and 96 hpt (i.e., 24 and 48 hpi, respectively; Figure 1). Relief algorithm combined with OPLS-DA analysis allowed the separation of W72i vs. P72i and H72i, H72i vs. P72i and W72i, and P72i vs. H72i and W72i with a significant *p*-value < 0.05 (Table 1). W96i, H96i, and P96i were similarly separated with a significant *p*-value < 0.05 (Table 1).

A total of 2,102 m/z was recorded for the 72i sample extracts. Among them, 83, 30, and 83 were specific to W72i, H72i, and P72i, respectively, whereas 402 were common to the three groups (Supplementary Figure S3A). As previously described for the 48 sample sets, the top m/z were identified and used for comparison. A total of 365 ones were obtained (Table 2), among with 70, 25, and 58 were specific to W72i, H72i, and P72i (Supplementary Figure S3B), i.e., 61, 27, and 36% of their total top m/z. The corresponding chemical composition and VK diagrams highlighted differences between W72i, H72i, and P72i (Supplementary Figures S4A-C). Compared to W72i and H72i, P72i metabolic fingerprint was characterized by a low level of CHO and CHOP formulas, and a high level of CHOS ones (Supplementary Figures 4A-C). The other formulas were few represented and distributed similarly to H72i ones. Compared to W72i, H72i had less CHOS, CHONP, and CHONSP formulas, and more CHOP ones (Supplementary Figures S4A,B). VK diagrams highlighted specific and highly different metabolic fingerprints for each sample set (Supplementary Figures S4A-C).

Query of KEGG database with the MassTRIX interface showed that “Pentose and glucuronate interconversions” and “Ascorbate and aldarate metabolism” were the main pathways associated to W72i top m/z (Figure 5A). “Biosynthesis of secondary metabolites” and “Flavonoid metabolism” were preferentially associated to H72i (Figure 5B) whereas the main pathways in P72i samples were “Purine metabolism” and “Pyrimidine metabolism” (Figure 5C). Among annotated masses were ADP, ATP, UMP, UDP for P72i and quercetin 7-O-glucoside, leucocyanidin, leucodelphinidin, 4-coumaroylshikimate, and ampelopsin for H72i (Supplementary Table S2).

A total of 1,808 m/z was recorded for the 96i sample extracts. Among them, 8, 25, and 153 were specific to W96i, H96i, and P96i, respectively, whereas 1,093 were common to the three series (Supplementary Figure S3C). A total of 261 top m/z were isolated (Table 2), among which 8, 25, and 49 were specific to W96i, H96i, and P96i, respectively (i.e., 22, 32, and 34% of their total top m/z) (Supplementary Figure S3D). P96i fingerprint was characterized by higher numbers of CHO, CHOS, and CHON formulas than W96i and H96i (Supplementary Figures S4D-F). The other chemical compositions were slightly present (CHOP, CHONS, CHONSP) or absent (CHONP). H96i displayed a majority of CHO and, to a lesser extent, of CHONSP formulas whereas all the other ones were meaningless. As few top m/z were obtained for
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

**FIGURE 3** | Detailed visualization of the 48 hpt samples analyzed by FT-ICR-MS. Chemical histograms show the distribution (number) of the elemental formulas attributed to the top *m/z* of W48 (A), H48 (B), and P48 (C) samples analyzed by (-) FT-ICR-MS, according to their elemental composition (CHO, CHOS, CHON, CHONS, CHOP, CHONP, or CHONSP). Their corresponding Van Krevelen diagrams (D-F) represent these formulas onto two axes according H/C and O/C atomic ratios. Dots are colored according to their elemental composition (CHO, CHOS, CHON, CHONS, CHOP, CHONP, CHONSP) and sized according to their relative intensity in mass spectra. W48, H48, and P48 correspond to methanol extracts of grapevine leaf disks treated with H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g·l⁻¹) for 48 h. Top *m/z* correspond to *m/z* with the highest regression coefficient value (VIP ≥ 1).

W96i, rare or no formulas were found in the different chemical compositions. Query of KEGG database with the MassTRIX interface showed that the main pathways for W96i were “pentose and glucuronate interconversions”, and “ascorbate and aldarate metabolism” (Figure 5D), as for W72i (Figure 5A). For H96i, the main pathways were more diversified and mainly encompassed “biosynthesis of secondary metabolites” “carbon metabolism”, “carbon fixation in photosynthetic organisms”, and “biosynthesis of amino acids”, “pentose phosphate pathways”, and “pentose and glucuronate interconversions” (Figure 5E). “Linoleic acid metabolism” and “purine metabolism” were the main pathways for P96i (Figure 5F). Among the annotated masses were ADP, tetradecanoic acid, palmitoleic acid and hexadecanoic acid for P96i, d-ribose 5-P, d-ribose 5-diP, sedoheptulose-1,7-diP for H96i (Supplementary Table S2).

## Erythritol Phosphate: A Putative Marker of PS3-IR

Two metabolites previously reported as markers of grapevine defenses, resveratrol and ε-viniferin, were first searched in the data files. Putative resveratrol (*m/z* 227.0714) was present in

all samples, and transiently accumulated in response to elicitor treatment and infection (Figure 6). Putative ε-viniferin (*m/z* 453.1345) was only induced in response to H treatment at 48 and 72 hpt, and in response to *P. viticola* infection at 72 hpt (for W and P treated leaves) and 96 hpt (for P treated leaves) (Figure 6).

In order to find a putative marker of elicitor-IR, annotated top *m/z* specific to Pi series and present in both P72i and P96i were next investigated. Among the 58, and 49 top *m/z* specific to P72i and P96i, only 5 and 8 could be annotated, respectively. A single one, designed as erythritol phosphate (*m/z* = 201.0169), was common to both sample sets. The assessment of the presence of this compound in all samples confirmed that it specifically accumulated in P treated leaves (Figure 6). A network analysis based on a series of biochemical transformations (Roullier-Gall et al., 2014b) was performed to obtain an additional overview of the global metabolism of grapevine leaves (Figure 7A). A focus on erythritol phosphate *m/z* (Figure 7B) highlighted correlated elemental formulas. About half of them could be annotated with Masstrix and home made database. A large number of them was annotated to pentose phosphate pathway (Supplementary Table S3).
<figure><img src="image_11.png" /></figure>

**FIGURE 4** | Metabolic pathways associated to the top *m/z* identified in the 48 hpt samples. The metabolic pathways correspond to the top *m/z* of W48 (red), H48 (green), and P48 (blue) samples. These pathways were obtained after KEGG query with the MassTRIX interface using *Vitis vinifera* organism (vi). The x axis corresponds to the number of annotated formulas obtained for each pathway. Only pathways including at least five annotated formulas are presented. W48, H48, and P48 correspond to methanol extracts of grapevine leaf disks treated with H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g.l⁻¹) for 48 h. Top *m/z* correspond to *m/z* with the highest regression coefficient value (VIP ≥ 1).

DISCUSSION

In order to better characterize elicitor-IR, a non-targeted analysis
was performed to obtain a metabolite fingerprint associated to

elicitor-IR and also to identify putative markers of elicitor-IR.
FT-ICR-MS analyses were performed from methanolic extracts of
leaf samples and in the negative mode. This method of extraction
allows access to a wide number of *m/z* (Maia et al., 2016).
<figure><img src="image_12.png" /></figure>

**FIGURE 5** | Metabolic pathways associated to the top *m/z* identified in the 72i and 96i samples. The metabolic pathways correspond to the top *m/z* of W72i (A), H72i (B), P72i (C), W96i (D), H96i (E), and P96i (F) samples. These pathways were obtained after KEGG query with the MassTRIX interface using *Vitis vinifera* organism (vvi). The x axis correspond to the number of annotated formulas obtained for each pathway. Only pathways including at least five annotated formulas are presented. W72i/96i, H72i/96i, and P72i/96i correspond to methanol extracts of grapevine leaf disks treated with distilled H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g.l⁻¹) for 48 h and collected 72 and 96 h after *P. viticola* inoculation. Top *m/z* correspond to *m/z* with the highest regression coefficient value (VIP ≥ 1).

Supervised OPLS-DA analyses allowed the discrimination of the different sample series. The first sampling time used for comparison was 48 hpt, just before inoculation. This time was long enough to allow the accumulation of defense compounds (Trouvelot et al., 2008; Gauthier et al., 2014). W48, H48, and P48 sample series could be discriminated, indicating that grapevine leaf metabolism was differently modulated in response to these treatments. W48 had few top *m/z* with a dominance of CHOS containing metabolites. The presence of sulfur in these formulas was confirmed by highly exact mass measurements, including isotopic pattern analyses. However, it must be born in mind that elemental formulas do not represent all of the detected *m/z* features, but only those top features, which weigh for the discrimination presented. The

significance of this sulfur-containing signature certainly merits to be investigated in future experiments. H48 had the highest number of top *m/z*, indicating that a large number of chemically diverse metabolites contributed to H48 series discrimination. This was confirmed by the diversity of the activated metabolic pathways and suggested a global activation of the cell metabolism. The activation of carbohydrate metabolism suggested that energy may be recruited to fuel defense reactions, as previously reported (Bolton, 2009; Rojas et al., 2014). H48 metabolic fingerprint also involved "Secondary metabolites" and "Flavonoid biosynthesis", indicating that the plant is in a defensive status. Despite the *m/z* corresponding to the grapevine phytoalexin resveratrol (*m/z* = 227,0714046) was significantly present in H48 and P48 samples series (Figure 6), it was not identified among the top
<figure><img src="image_13.png" /><figcaption>**FIGURE 6** | Occurrence of resveratrol, ε-viniferin and erythritol P in the different samples. The signal intensities of putative resveratrol (m/z = 227, 0714), ε-viniferin (m/z = 453,1345) and erythritol-P (m/z = 201.0169) are expressed in relative intensity. W, H, and P correspond to methanol extracts of grapevine leaf disks treated with distilled H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g.I⁻¹) for 48 h before *P. viticola* (i) or mock (ni) inoculation. Disks were collected at 48, 72, and 96 hpt.</figcaption></figure>

*m/z*. Conversely, its glycosylated and dimer derivatives piceid (or polydatine) and ε-viniferin, respectively, were listed among the annotated top *m/z* attributed to H48 samples **(Supplementary Table S1)**.

As for W48, a limited number of metabolic pathways was associated to P48 top *m/z*. As reminded above, PS3 is known to act by priming and little is known about its effects in plants before pathogen infection. Gauthier et al. (2014) have previously reported transcriptional changes in response to PS3 treatment but this is the first time that a direct effect of a priming elicitor is observed at the metabolome level. P48 metabolic response involved “purine metabolism”, with AMP, ADP, ATP, UTP annotations, suggesting the constitution of a pool of nucleotides and the mobilization of energy sources and/or cofactors for enzyme activities. Traumatic acid (or trans-2-dodeceno-1.12-dicarboxylic acid) and its precursor linolenic acid were among the annotated top *m/z*. Despite traumatic acid remains little studied, in plant, it is known to be involved in response to injury by stimulation of cell division (Zimmerman and Coudron, 1979). The production of traumatic acid and other 12C derivatives of the hydroperoxide lyase pathway was also recently studied in *Nicotiana attenuata* leaves in response to wounding (Kallenbach et al., 2011). Traumatic acid could also play an important role in plant adaptation to environmental stresses, as reported by Pietryczuk et al. (2014) for *Chlorella vulgaris*. It would be interesting to further analyze its putative role in elicitor-IR. Other compounds such as JA and abscisic acid were also among the annotated top *m/z*, suggesting that the alpha-linolenic acid cascade (Rojas et al., 2014) was largely

involved in response to PS3 treatment. These results confirm the involvement of JA in PS3-IR, as previously demonstrated by Trouvelot et al. (2008). Besides, Gauthier et al. (2014) have reported that grapevine leaves treated by PS3 share a common stress-responsive transcriptome profile that partly overlaps the salicylate- (SA-) and JA-dependent one, although the microarray dataset was obtained at 12 hpt. The precise role of SA and JA in PS3-IR of grapevine against downy mildew thus remains to be elucidated.

Leaf samples were mock- or *P. viticola* inoculated at 48 hpt and collected for analysis at 72 and 96 hpt (corresponding to 24 and 48 hpi, respectively). As previously reported (Allègre et al., 2009), only PS3 IR against downy mildew (with mean values of disease severity of about 25 for W and H samples, and 4 for P ones at 144 hpt; data not shown). The global metabolic profile of mock-inoculated W samples did not change during the time course of experiment (data not shown), thus providing evidence for the “stability” of the global metabolism during the time course of the experiment. *P. viticola* infection clearly impacted the leaf metabolism. As it developed into W- and H-treated leaves, metabolites could also be the result of its own biosynthesis. However, samples were collected at the beginning of infection, when the biomass of the oomycete was low (Gamm et al., 2011), thus suggesting that it was rather due to plant response to infection. Ali et al. (2009) used NMR spectroscopy and (2D)-NMR techniques to obtain a metabolic profiling of susceptible *vs* resistant grapevine leaves upon inoculation with *P. viticola* at different time points. However, they did not provide a comparative global metabolite
<figure><img src="image_14.png" /></figure>

profile at 24 and 48 hpi, making it impossible to determine if a similar response was observed in their downy mildew susceptible variety.

All infected samples could also be discriminated at 72 and 96 hpt (72i and 96i samples) and presented fingerprints different to those observed at 48 hpt. Despite a time- and

infection-dependent evolution, each sample set thus kept a specific metabolite fingerprint. Whatever the treatment, the number of top *m/z* decreased from 48 to 72 to 96 hpt, indicating that as time post treatment / post inoculation progressed, less masses marked the metabolic fingerprints. Interestingly, the weaker decrease of this number of top *m/z* was observed

<table><caption>TABLE 3 | Summary of the main metabolic pathways associated to annotated top *m/z* identified for each treatment and sampling time.</caption><thead><tr><th></th><th>W</th><th>H</th><th>P</th></tr></thead><tbody><tr><td>48</td><td>Biosynthesis of secondary metabolites</td><td><strong>Biosynthesis of secondary metabolites</strong><br><strong>Pentose and glucuronate interconversions</strong><br><strong>Ascorbate and aldarate metabolism</strong><br>Amino sugar and nucleotide sugar metabolism<br><strong>Flavonoid biosynthesis</strong><br>Starch and sucrose metabolism</td><td>Biosynthesis of secondary metabolites<br><strong>Purine metabolism</strong></td></tr><tr><td>72i</td><td><strong>Pentose and glucuronate interconversions</strong><br><strong>Ascorbate and aldarate metabolism</strong></td><td><strong>Biosynthesis of secondary metabolites</strong><br><strong>Flavonoid biosynthesis</strong></td><td>Pyrimidine metabolism<br><strong>Purine metabolism</strong></td></tr><tr><td>96i</td><td><strong>Pentose / glucuronate interconversions</strong><br><strong>Ascorbate and aldarate metabolism</strong></td><td><strong>Biosynthesis of secondary metabolites</strong><br>Carbon fixation in photosynthetic organisms<br>Carbon metabolism<br>Biosynthesis of amino acids<br>Pentose phosphate pathway<br><strong>Pentose and glucuronate interconversions</strong></td><td>Linoleic acid metabolism<br><strong>Purine metabolism</strong></td></tr></tbody></table>

In bold are highlighted the redundant ones. W, H, and P correspond to methanol extracts of grapevine leaf disks treated with H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g.l⁻¹) for 48 h.
in response to PS3 treatment, i.e., in leaves where IR was
effective.

H72i and H96i leaves did not present the same metabolic
fingerprint as W72i and W96i ones, although *P. viticola* infection
managed to develop inside them. These results show that
despite H13 failed in inducing resistance, it nonetheless durably
modulated the leaf metabolism. This is partly due to defense
activation, as indicated by the occurrence of the phytoalexin
resveratrol and derivatives, and also to the metabolic costs
probably associated to this response.

The metabolic fingerprint of the three infected sample
series changed upon time. The main pathways associated
to each treatment and sampling time are summarized in
**Table 3**. One of the main metabolic pathways associated to
H72i and H96i was "secondary metabolism", as for H48. The
production of secondary metabolites was induced by H13
treatment and probably also by *P. viticola* infection. A greater
number of metabolic pathways was associated to H48 and
H96i, and included "carbon metabolism", "carbon fixation
in photosynthetic organism", "pentose phosphate pathway",
and "biosynthesis of amino acid". These results suggest a
mobilization of the primary metabolism, with sugar and amino
acid biosynthesis, and maybe photorespiration that represents
a converging point for carbohydrate and amino acid metabolic
pathways. These different pathways were reported for playing
a role in plant defense responses (Rojas et al., 2014). They
ensure the production of secondary metabolites and the
associated energy costs. However, all these response failed to
prevent the pathogen development. The "pentose glucuronate
interconversions" pathway was mobilized in W72i, W96i and
H96i samples and might highlight the plant response to *P. viticola*
infection process. The metabolic response of P72i and P96i
involved the "purine metabolism" pathway, as for P48, suggesting
an active gene transcription and energy mobilization. The lipid
metabolism was also involved in response to PS3 treatment,
especially at 96 hpt, with masses annotated as linoleic acid,
tetradecanoic acid, hexadecanoic acid. This means a possible
membrane and cell wall rearrangement and the synthesis of fatty
acid derivative compounds such as the signaling hormone JA
which are key events associated to plant response to biotic cues
(Rojas et al., 2014).

Metabolomics thereby provided us a better understanding of
elicitor-induced resistance of grapevine leaves against *P. viticola*.
Ballester et al. (2013) also used this approach to highlight the
involvement of the phenylpropanoid pathway in the induction
of resistance of Navelate oranges. Nevertheless, the use of
metabolomics to investigate elicitor-induced resistance remains
little documented. More studies have focused on plant interaction
with microbes and pests (for review, see Balmer et al., 2013;
Shiratake and Suzuki, 2016). As example, Hong et al. (2011)
could discriminate healthy and botrytized berries of botrytized
bunches, and berries of healthy bunches. All the analytical
methods do not have the same performance and accuracy (Jorge
et al., 2016; Maia et al., 2016; Shiratake and Suzuki, 2016)
and plant/pathogen interactions have their specificity; making it
difficult to compare the different results reported in the literature.
However, a common observation revealed by metabolomic

studies is the reprogramming of the plant metabolism, with an
impact on both primary (often carbohydrates, amino acids) and
secondary metabolites during interaction. This was described
in barley, rice and *Brachypodium distachyon* in response to
*Magnaporthe grisea* (Parker et al., 2009) or in rice in response
to *Chilo suppressalis* attack (Liu et al., 2016). In grapevine, it
was reported in *V. vinifera* cv. Trincadeira berries challenged
by *Botrytis cinerea* (Agudelo-Romero et al., 2015) and leaves
infected by *Plasmopara viticola* (Ali et al., 2009; Becker et al.,
2013). Metabolomic also allows the identification of biomarkers
of defensive or infected state. As example, Chin et al. (2014)
reported biomarkers of citrus infection by the bacterium
*Candidatus Liberibacter asiaticus*.

Putative erythritol phosphate (m/z = 201,01698, C₁₁H₄O₇P)
was the unique annotated compound specifically listed among
the top m/z of PS3 treated samples. Little is known about
the biological activity of this phosphorylated polyol. It was
demonstrated that *Brucella* protobacteria can convert d-
erythrose-4-phosphate to glyceraldehyde 3-phosphate and
fructose 6-phosphate by enzymes of the pentose phosphate
pathway, thus bypassing fructose-1,6-bisphosphatase (Barbier
et al., 2014). In plants, its derivative methylerythritol phosphate
was more studied as it constitutes an important pathway for the
biosynthesis of isopentenyl diphosphate, a precursor of terpenes.
As example, the synthesis of the diterpenoid phytoalexin in
rice in response to elicitors (including JA) is mediated by this
pathway (Okada et al., 2007). In grapevine, this pathway is also
involved in the production of monoterpenes and sequiterpenes
in response to methyljasmonate (Hampel et al., 2005). In this
study, putative methylerythritol phosphate was present in all
samples. Further studies will be necessary to confirm the identity
of erythritol phosphate, to understand its regulation related to
its methylated derivative, and to determine its role in PS3-IR.
The reliability of this marker will have to be checked in other
pathosystems and in response to other plant defense elicitors.
Such a marker of elicitor-IR would be useful as a tool of selection
for further screening of putative new elicitors. It would be also
helpful to follow the plant response to elicitor treatment in field
conditions and to study the impact of environmental conditions
on elicitor-IR.

CONCLUSION

The metabolomic approach used in this study allowed us to
obtain specific and time-dependent metabolite signatures of PS3-
IR, and also to identify a putative marker of PS3-IR. Despite
*m/z* annotations and metabolic pathways could be associated to
these signatures, this constitutes only partial information since
numerous *m/z* remain "unknown" in databases. However, these
results pave the way for future studies that will contribute to
enrich our knowledge of elicitor-IR.

AUTHOR CONTRIBUTIONS

MA designed the work, conducted the experiments, contributed
to data interpretation, and drafted the manuscript. BK performed
FT-ICR-MS analysis. ML and CR-G performed data process and statistical analysis. BP, CL-G, M-CH, RG, PS-K, ST, and XD contributed to design the work and to interpret data. PS-K supervised the FT-ICR-MS analysis, and contributed to data interpretation. All authors reviewed, edited, and approved manuscript.

## ACKNOWLEDGMENTS

The authors thank Goëmar laboratories for providing the elicitors, BIVB (Bureau Interprofessionnel des Vins de Bourgogne) and the Conseil Régional de Bourgogne for financial support.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.00101/full#supplementary-material

**FIGURE S1 | Detailed visualization of all leaf samples analyzed by FT-ICR-MS.** The chemical histogram (left part) shows the distribution (number) of the elemental formulas attributed to the *m/z* obtained from all sample extracts according to their elemental composition (CHO, CHOS, CHON, CHONS, CHOP, CHONP, or CHONSP) and the Van Krevelen diagram (right part) represents these formulas onto two axes according to their H/C and O/C atomic ratios. Dots are colored according to their elemental composition (CHO, CHOS, CHON, CHONS, CHOP, CHONP, CHONSP) and sized according to their relative intensity in mass spectra.

**FIGURE S2 | Detailed visualization of the 48 hpt samples analyzed by FT-ICR-MS.** The chemical histograms show the distribution (number) of the elemental formulas attributed to the *m/z* obtained from W48, H48, and P48

samples, respectively, according to their elemental composition (CHO, CHOS, CHON, CHONS, CHOP, CHONP, or CHONSP). The Van Krevelen diagrams represent these formulas onto two axes according to their H/C and O/C atomic ratios. Dots are colored according to their elemental composition (CHO, CHOS, CHON, CHONS, CHOP, CHONP, CHONSP) and sized according to their relative intensity in mass spectra. W48, H48, and P48 correspond to methanol extracts of grapevine leaf disks treated with H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g·l⁻¹) for 48 h.

**FIGURE S3 | Discrimination of the 72 and 96i samples from FT-ICR-MS data.** Venn diagrams showing the distribution of the total (A,C) and top (B,D) *m/z* obtained by (-) FT-ICR-MS analysis of W, H, P infected samples collected at 72 and 96 hpt, respectively. W, H, and P correspond to methanol extracts of grapevine leaf disks treated with H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g·l⁻¹) for 48 h and inoculated with *P. viticola*. Top *m/z* correspond to *m/z* with the highest regression coefficient value (VIP ≥ 1).

**FIGURE S4 | Detailed visualization of the 72i and 96i samples analyzed by FT-ICR-MS.** The Van Krevelen diagrams represent the elemental formulas attributed to the top *m/z* of W72i (A), H72i (B), P72i (C), W96i (D), H96i (E), and P96i (F) samples onto two axes according to their H/C and O/C atomic ratios. Dots are colored according to their elemental composition (CHO, CHOS, CHON, CHONS, CHOP, CHONP, CHONSP) and sized according to their relative intensity in mass spectra. Diagrams (in the bottom left corner) show the distribution of these formulas according to their elemental composition (CHO, CHOS, CHON, CHONS, CHOP, CHONP, or CHONSP). W72i/96i, H72i/96i, and P72i/96i correspond to methanol extracts of grapevine leaf disks treated with distilled H₂O (W, as control), H (H13), or P (PS3) elicitor solutions (2.5 g·l⁻¹) for 48 h and collected 72 and 96 h after *P. viticola* inoculation. Top *m/z* correspond to *m/z* with the highest regression coefficient value (VIP ≥ 1).

**TABLE S1 | Annotations of W48, H48 and P48 top *m/z* after query of KEGG database (Vitis vinifera organism) with the MassTRIX interface.**

**TABLE S2 | Annotations of 72i and 96i top *m/z* for W, H, and P series after query of KEGG database (V. vinifera organism) with the MassTRIX interface.**

**TABLE S3 | list of the transformations targeted on erythritol phosphate metabolism.**

## REFERENCES

Adrian, M., Trouvelot, S., Gamm, M., Poinssot, B., Héloir, M. C., and Daire, X. (2012). "Activation of grapevine defense mechanisms theoretical and applied approaches," in *Plant Defense: Biological control, Progress in Biological control* 12, eds J. M. Merillon and K. G. Ramawat (Dordrecht: Springer).

Agudelo-Romero, P., Erban, A., Rego, C., Carbonell-Bejerano, P., Nascimento, T., Sousa, L., et al. (2015). Transcriptome and metabolome reprogramming in *Vitis vinifera* cv. Trincadeira berries upon infection with *Botrytis cinerea*. *J. Exp. Bot.* 66, 1769-1785. doi: 10.1093/jxb/eru517

Ali, K., Maltese, F., Zyprian, E., Rex, M., Choi, Y. E., and Verpoorte, R. (2009). NMR metabolic fingerprinting based identification of grapevine metabolites associated with downy mildew resistance. *J. Agric. Food Chem.* 57, 9599-9606. doi: 10.1021/jf902069f

Allègre, M., Héloir, M. C., Trouvelot, S., Daire, X., Pugin, A., Wendehenne, D., et al. (2009). Are grapevine stomata involved in the elicitor-induced protection against downy mildew? *Mol. Plant. Microbe Interact.* 22, 977-986. doi: 10.1094/MPMI-22-8-0977

Aziz, A., Poinssot, B., Daire, X., Adrian, M., Bézier, A., Lambert, B., et al. (2003). Laminarin elicits defense responses in grapevine and induces protection against *Botrytis cinerea* and *Plasmopara viticola*. *Mol. Plant. Microbe Interact.* 16, 1118-1128. doi: 10.1094/MPMI.2003.16.12.1118

Ballester, A. R., Lafuente, M. T., de Vos, C. H. R., Bovy, A. G., and González-Candelas, L. (2013). Citrus phenylpropanoids and defence against pathogens. Part I: metabolic profiling in elicited fruits. *Food Chem.* 136, 178-185. doi: 10.1016/j.foodchem.2012.07.114

Balmer, D., Flors, V., Glauser, G., and Mauch-Mani, B. (2013). Metabolomics of cereals under biotic stress: current knowledge and techniques. *Front. Plant Sci.* 4:82. doi: 10.3389/fpls.2013.00082

Barbier, T., Collard, F., Zuniga-Ripa, A., Moriyon, I., Godard, T., Becker, J., et al. (2014). Erythritol feeds the pentose phosphate pathway via three new isomerases leading to D-erythrose-4-phosphate in *Brucella*. *Proc. Natl. Acad. Sci. U.S.A* 111, 17815-17820. doi: 10.1073/pnas.1414622111

Batovska, D. I., Todorova, I., Parushev, S., Valentinova Nedelcheva, D., Stefanova Bankova, V., et al. (2009). Biomarkers for the prediction of the resistance and susceptibility of grapevine leaves to downy mildew. *J. Plant Physiol.* 166, 781-785. doi: 10.1016/j.jplph.2008.08.008

Becker, L., Poutaraud, A., Hamm, G., Muller, J. F., Merdinoglu, D., Carré, V., et al. (2013). Metabolic study of grapevine leaves infected by downy mildew using negative ion electrospray-Fourier transform ion cyclotron resonance mass spectrometry. *Anal. Chim. Acta* 795, 44-51. doi: 10.1016/j.aca.2013.07.068

Bisson, L. F., Waterhouse, A. L., Ebeler, S. E., Walker, M. A., and Lapsley, J. T. (2002). The present and future of the international wine industry. *Nature* 418, 696-699. doi: 10.1038/nature01018

Bolton, M. D. (2009). Primary metabolism and plant defense - fuel for the fire. *Mol. Plant Microbe Interact.* 22, 487-497. doi: 10.1094/MPMI-22-5-0487

Buescher, J. M., and Driggers, E. M. (2016). Integration of omics: more than the sum of its parts. *Cancer Metab.* 19, 4. doi: 10.1186/s40170-016-0143-y

Casida, J. E. (2009). Pest toxicology: the primary mechanisms of pesticide action. *Chem. Res. Toxicol.* 4, 609-619. doi: 10.1021/tx8004949

Chin, E. L., Mishchuk, D. O., Breksa, A. P., and Slupsky, C. M. (2014). Metabolite signature of *Candidatus Liberibacter asiaticus* infection in two citrus varieties. *J. Agric. Food Chem.* 62, 6585-6591. doi: 10.1021/jf5017434
Delaunois, B., Farace, G., Jeandet, P., Clément, C., Baillieul, F., Dorey, S., et al. (2013). Elicitors as alternative strategy to pesticides in grapevine? Current knowledge on their mode of action from controlled conditions to vineyard. *Environ. Sci. Pollut. Res.* 21, 4837–4846. doi: 10.1007/s11356-013-1841-4

Flamini, R., de Rosso, M., de Marchi, F., Dalla Vedova, A., Panighel, A., Gardiman, M., et al. (2013). An innovative approach to grape metabolomics: stilbene profiling by suspect screening analysis. *Metabolomics* 9, 1243–1253. doi: 10.1007/s11306-013-0530-0

Flamini, R., De Rosso, M., Panighel, A., Dalla Vedova, A., De Marchi, F., and Bavaresco, L. (2014). Profiling of grape monoterpene glycosides (aroma precursors) by ultra-high performance-liquid chromatography-high resolution mass spectrometry (UHPLC/QTOF). *J. Mass Spectrom.* 49, 1214–1222. doi: 10.1002/jms.3441

Gamm, M., Héloir, M. C., Bligny, R., Vaillant-Gaveau, N., Trouvelot, S., Alcaraz, G., et al. (2011). Changes in carbohydrate metabolism in *Plasmopara viticola* -infected grapevine leaves. *Mol. Plant Microbe Interact.* 24, 1061–1073. doi: 10.1094/MPMI-02-11-0040

Garcia-Brugger, A., Lamotte, O., Vandelle, E., Bourque, S., Lecourieux, D., Poinssot, B., et al. (2006). Early signalling events induced by elicitors of plant defenses. *Mol. Plant Microbe Interact.* 19, 711–724. doi: 10.1094/MPMI-19-0711

Gauthier, A., Trouvelot, S., Kelloniemi, J., Frettinger, P., Wendehenne, D., Daire, X., et al. (2014). The sulfated laminarin triggers a stress transcriptome before priming the SA- and ROS-dependent defenses during grapevine's induced resistance against *Plasmopara viticola*. *PLoS ONE* 9:e88145. doi: 10.1371/journal.pone.0088145

Gougeon, R. D., Lucio, M., Frommberger, M., Peyron, D., Chassagne, D., Alexandre, H., et al. (2009). The chemodiversity of wines can reveal a metabologeography expression of cooperage oak wood. *Proc. Natl. Acad. Sci. U.S.A* 106, 9174–9179. doi: 10.1073/pnas.0901100106

Griesser, M., Weingart, G., Schoedl-Hummel, K., Neumann, N., Becker, M., Varmuza, K., et al. (2015). Severe drought stress is affecting selected primary metabolites, polyphenols, and volatile metabolites in grapevine leaves (*Vitis vinifera* cv. Pinot noir). *Plant Physiol. Biochem.* 88, 17–26. doi: 10.1016/j.plaphy. 2015.01.004

Hampel, D., Mosandi, A., and Wüst, M. (2005). Induction of de novo volatile terpene biosynthesis via cytosolic and plastidial pathways by methyl jasmonate in foliage of *Vitis vinifera* L. *J. Agric. Food Chem.* 53, 2652–2657. doi: 10.1021/ jf040421q

Hochberg, U., Degu, A., Toubiana, D., Gendler, T., Nikoloski, Z., Rachmilevitch, S., et al. (2013). Metabolite profiling and network analysis reveal coordinated changes in grapevine water stress response. *BMC Plant Biol.* 20:184. doi: 10. 1186/1471-2229-13-184

Hong, Y. S., Cilindre, C., Liger-Belair, G., Jeandet, P., Hertkorn, N., and Schmitt-Kopplin, P. (2011). Metabolic influence of *Botrytis cinerea* infection in champagne base wine. *J. Agric. Food Chem.* 59, 7237–7245. doi: 10.1021/ jf200664t

Huot, B., Yao, J., Montgomery, B. L., and He, S. Y. (2014). Growth-Defense tradeoffs in plants: a balancing act to optimize fitness. *Mol. Plant* 7, 1267–1287. doi: 10.1093/mp/ssi049

Janz, D., Behnke, K., Schnitzler, J. P., Kanawati, B., Schmitt-Kopplin, P., and Polle, A. (2010). Pathway analysis of the transcriptome and metabolome of salt sensitive and tolerant poplar species reveals evolutionary adaption of stress tolerance mechanisms. *BMC Plant Biol.* 10:150. doi: 10.1186/1471-2229-10-150

Jorge, T. F., Mata, A. T., and Antonio, C. (2016). Mass spectrometry as a quantitative tool in plant metabolomics. *Philos. Trans. A Math. Phys. Eng. Sci.* 374, 20150370. doi: 10.1098/rsta.2015.0370

Kallenbach, M., Gilardoni, P. A., Allmann, S., Baldwin, I. T., and Bonaventure, G. (2011). C12 derivatives of the hydroperoxide lyase pathway are produced by product recycling through lipoxygenase-2 in *Nicotiana attenuata* leaves. *New Phytol.* 191, 1054–1068. doi: 10.1111/j.1469-8137.2011.03767.x

Kim, J., Woo, H. R., and Nam, H. G. (2016). Toward systems understanding of leaf senescence: an integrated multi-omics perspective on leaf senescence research. *Mol. Plant* 9, 813–825. doi: 10.1016/j.molp.2016.04.017

Kim Khiook, I. L., Schneider, C., Héloir, M. C., Bois, B., Daire, X., Adrian, M., et al. (2013). Image analysis methods for assessment of H2O2 production and *Plasmopara viticola* development in grapevine leaves: application to the

evaluation of resistance to downy mildew. *J. Microb. Meth.* 28, 235–244. doi: 10.1016/j.jimet.2013.08.012

Liu, Q., Wang, X., Tzin, V., Romeis, J., Peng, Y., and Li, Y. (2016). Combined transcriptome ans metabolome analyses to understand the dynamic responses of rice plants to attack by the rice stem borer *Chilo suppressalis* (Lepidoptera: Crambidae). *BMC Plant Biol.* 16:259. doi: 10.1186/s12870-016- 0946-6

Maia, M., Monteiro, F., Sebastiana, M., Marques, A. P., Ferreira, A. E. N., Freire, A. P., et al. (2016). Metabolite extraction for high-throughput FTICR-MS- based metabolomics of grapevine leaves. *EuPA Open Proteom.* 12, 4–9. doi: 10.1016/j.euprot.2016.03.002

Matus, J. T. (2016). Transcriptomic and metabolomic networks in the grape berry illustrate that it takes more than flavonoids to fight against ultraviolet radiation. *Front. Plant Sci.* 7:1337. doi: 10.3389/fpls.2016.01337

Ménard, R., Alban, S., de Ruffay, P., Jamois, F., Franz, G., Fritig, B., et al. (2004). β-1,3 glucan sulfate, but not β-1,3 glucan, induces the salicylic acid signaling pathway in tobacco and Arabidopsis. *Plant Cell* 16, 3020–3032. doi: 10.1105/ tpc.104.024968

Nanni, I. M., Pirondi, A., Mancini, D., Stammler, G., Gold, R., Ferri, I., et al. (2016). Differences in the efficacy of Carboxylic Acid Amides (CAA) fungicides against less sensitive strains of *Plasmopara viticola*. *Pest Manag. Sci.* 72, 1537–1539. doi: 10.1002/ps.4182

Okada, A., Shimizu, T., Okada, K., Kuzuyama, T., Koga, J., Shibuya, N., et al. (2007). Elicitor induced activation of the methylerythritol phosphate pathway toward phytoalexins biosynthesis is rice. *Plant Mol. Biol.* 65, 177–187.

Paris, F., Krzyżaniak, Y., Gauvrit, C., Jamois, F., Domergue, F., Joubert, J. M., et al. (2015). An ethoxylated surfactant enhances the penetration of the sulfated laminarin through leaf cuticle and stomata, leading to increased induced resistance against grapevine downy mildew. *Physiol. Plant* 156, 338–350. doi: 10.1111/ppl.12394

Parker, D., Beckmann, M., Zubair, H., Enot, D. P., Caracuel-Rios, Z., Overy, D. P., et al. (2009). Metabolic analysis reveals a common pattern of metabolic re-programming during invasion of three host species by *Magnaporthe grisea*. *Plant J.* 59, 723–737. doi: 10.1111/j.1365-313X.2009. 03912.x

Pietryczuk, A., Biziewska, I., Imierska, M., and Czerpak, R. (2014). Influence of traumatic acid on growth and metabolism of *Chlorella vulgaris* under conditions of salt stress. *Plant Growth Regul.* 73, 103–110. doi: 10.1007/s10725- 013-9872-x

Rojas, C. M., Senthil-Kumar, M., and Tzin, V. K. S. (2014). Regulation of primary plant metabolism during plant-pathogen interactions and its contribution to plant defense. *Front. Plant Sci.* 5:17. doi: 10.3389/fpls.2014.00017

Roullier-Gall, C., Boutegrabet, L., Gougeon, R. D., and Schmitt-Kopplin, P. (2014a). A grape and wine chemodiversity comparison of different appellations in Burgundy: Vintage vs terroir effects. *Food Chem.* 152, 100–107. doi: 10.1016/ j.foodchem.2013.11.056

Roullier-Gall, C., Wittig, M., Gougeon, R. D., and Schmitt-Kopplin, P. (2014b). High precision mass measurements for wine metabolomics. *Front. Chem.* 2:102. doi: 10.3389/fchem.2014.00102

Shigenaga, A. M., and Argueso, C. T. (2016). No hormone to rule them all: interactions of plant hormones during the responses of plants to pathogens. *Semin. Cell Dev. Biol.* 56, 174–189. doi: 10.1016/j.semcdb.2016. 06.005

Shiratake, K., and Suzuki, M. (2016). Omics studies of citrus, grape and rosaceae fruit trees. *Breed Sci.* 66, 122–138. doi: 10.1270/jsbbs.66.122

Suzuki, M., Nakabayashi, R., Ogata, Y., Sakurai, N., Tokimatsu, T., Goto, S., et al. (2015). Multiomics in grape berry skin revealed specific induction of the stilbene synthetic pathway by ultraviolet-C irradiation. *Plant Physiol.* 168, 47–59. doi: 10.1104/pp.114.254375

Trouvelot, S., Héloir, M. C., Poinssot, B., Gauthier, A., Paris, F., Guillier, C., et al. (2014). Carbohydrates in plant immunity and plant protection: roles and potential application as foliar sprays. *Front. Plant Sci.* 4:592. doi: 10.3389/fpls. 2014.00592

Trouvelot, S., Varnier, A. L., Allègre, M., Mercier, L., Baillieul, F., Arnould, C., et al. (2008). A β-1,3 Glucan sulfate induces resistance in grapevine against *Plasmopara viticola* through priming of defense responses, including HR-like
cell death. *Mol. Plant. Microbe Interact.* 21, 232–243. doi: 10.1094/MPMI-21-2-0232

Tziotis, D., Hertkorn, N., and Schmitt-Kopplin, P. (2011). Kendrick-analogous network visualisation of ion cyclotron resonance Fourier transform mass spectra: improved options for the assignment of elemental compositions and the classification of organic molecular complexity. *Eur. J. Mass Spectrom.* (Chichester) 17, 415–421. doi: 10.1255/ejms.113

Walters, D. R., Ratsep, J., and Havis, N. D. (2013). Controlling crop diseases using induced resistance: challenges for the future. *J. Exp. Bot.* 6, 1263–1280. doi: 10.1093/jxb/ert026

Wiesel, L., Newton, A. C., Elliott, I., Booty, D., Gilroy, E. M., Birch, P. R., et al. (2014). Molecular effects of resistance elicitors from biological origin and their potential for crop protection. *Front. Plant Sci.* 5:655. doi: 10.3389/fpls.2014.00655

Zimmerman, D. C., and Coudron, C. A. (1979). Identification of traumatin, a wound hormone, as 12-oxo-trans-10-dodecanoic acid. *Plant Physiol.* 63, 536–541. doi: 10.1104/pp.63.3.536

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

Copyright © 2017 Adrian, Lucio, Roullier-Gall, Héloir, Trouvelot, Daire, Kanawati, Lemaître-Guillier, Poinssot, Gougeon and Schmitt-Kopplin. 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.