Direct Contact – Sorptive Tape Extraction coupled with Gas Chromatography – Mass Spectrometry to reveal volatile topographical dynamics of lima bean (Phaseolus lunatus L.) upon herbivory by Spodoptera littoralis Boisd.

The dynamics of plant volatile (PV) emission, and the relationship between damaged area and biosynthesis of bioactive molecules in plant-insect interactions, remain open questions. Direct Contact-Sorptive Tape Extraction (DC-STE) is a sorption sampling technique employing non adhesive polydimethylsiloxane tapes, which are placed in direct contact with a biologically-active surface. DC-STE coupled to Gas Chromatography – Mass Spectrometry (GC-MS) is a non-destructive, high concentration-capacity sampling technique able to detect and allow identification of PVs involved in plant responses to biotic and abiotic stresses. Here we investigated the leaf topographical dynamics of herbivory-induced PV (HIPV) produced by Phaseolus lunatus L. (lima bean) in response to herbivory by larvae of the Mediterranean climbing cutworm (Spodoptera littoralis Boisd.) and mechanical wounding by DC-STE-GC-MS. Time-course experiments on herbivory wounding caused by larvae (HW), mechanical damage by a pattern wheel (MD), and MD combined with the larvae oral secretions (OS) showed that green leaf volatiles (GLVs) [(E)-2-hexenal, (Z)-3-hexen-1-ol, 1-octen-3-ol, (Z)-3-hexenyl acetate, (Z)-3-hexenyl butyrate] were associated with both MD and HW, whereas monoterpenoids [(E)-β-ocimene], sesquiterpenoids [(E)-nerolidol] and homoterpenes (DMNT and TMTT) were specifically associated with HW. Up-regulation of genes coding for HIPV-related enzymes (Farnesyl Pyrophosphate Synthase, Lipoxygenase, Ocimene Synthase and Terpene Synthase 2) was consistent with HIPV results. GLVs and sesquiterpenoids were produced locally and found to influence their own gene expression in distant tissues, whereas (E)-β-ocimene, TMTT, and DMNT gene expression was limited to wounded areas. DC-STE-GC-MS was found to be a reliable method for the topographical evaluation of plant responses to biotic and abiotic stresses, by revealing the differential distribution of different classes of HIPVs. The main advantages of this technique include: a) in vivo sampling; b) reproducible sampling; c) ease of execution; d) simultaneous assays of different leaf portions, and e) preservation of plant material for further “omic” studies. DC-STE-GC-MS is also a low-impact innovative method for in situ PV detection that finds potential applications in sustainable crop management.


Background
In the past ten years, the study of the interaction between larvae of the Mediterranean climbing cutworm (Spodoptera littoralis Boisd.) and leaves of the lima bean (Phaseolus lunatus L.) has provided evidence of both early and late events, and has been used as a model system to decipher plant-insect interactions [1][2][3][4][5]. Upon herbivory by S. littoralis, the lima bean responds, as do many other plants, with a cascade of events that lead to the activation of defense mechanisms. These mechanisms include the perception of molecular patterns or effectors of defense [6,7], mitogen-activated protein kinase (MAPK) activation, and protein phosphorylation [8,9], production of ethylene and jasmonates [10], expression of late defense response genes [11], and emission of herbivory-induced plant volatiles (HIPVs) [12,13].
Even if robotic mechanical wounding can simulate plant response similar to HIPV [4], the simple mechanical damage (MD) is not fully satisfactory to induce the same responses if not supported by the application of insect's oral secretions (OS) [14]. Despite the presence of several elicitors in S. littoralis OS (e.g., fatty acid conjugates) [7,15], it is not clear whether these factors originate from the salivary glands or other feeding-related organs, such as the ventral eversible gland [14,16]. However, the plant volatile (PV) blends emitted in response to herbivores differ markedly with different feeding modes [17][18][19][20].
In plant defensive strategies, the release of PVs plays multiple roles: direct deterrents against herbivores [21,22], attraction of natural enemies of the attacking herbivores [23][24][25][26], damage and disease long-distance signaling [27][28][29][30], and pathogen resistance priming [29][30][31][32]. Since volatiles are produced from several biosynthetic pathways, their qualitative and quantitative composition is the result of the concerted action of different pathways, triggered by multiple factors. To date, studies of the emission of PVs in response to herbivory have been limited to single organs or to the whole plant, either by destructive methods or by head-space analysis [33,34], and only one study analyzed PV gradients within a single leaf [35].
Direct Contact-Sorptive Tape Extraction (DC-STE) is a fast and easy-to-use sampling technique, developed to study the effect of cosmetic treatment on sebum composition, through in vivo sampling at the human skin surface [36,37]. The technique employs a thin flexible non-adhesive polydimethylsiloxane (PDMS) tape, which is placed directly in contact with a (biological) surface for a fixed time (Figure 1). Bicchi et al. [38] showed that this technique can also be applied to plants to monitor PVs, in both surface-static headspace and direct-contact (DC) modes. In DC-STE, volatiles produced at the biological surface are concentrated in the apolar PDMS layer by sorption (a sampling approach based on the partition of a compound between the sample and the bulk of a polymeric retaining phase) in amounts depending on the compound polarity and volatility. While in headspace sampling (e.g. static and dynamic headspace, high concentration-capacity solid phase microextraction) sorption is applied to the plant surrounding air space [33], DC-STE interacts directly with leaf surfaces. In DC-STE, plant-air interaction equilibrium is eliminated thus limiting the number of phases involved with sampling to two (plant and PDMS) instead of three (plant, air and PDMS). In this study, a glass coverslip was placed just above the DC-STE tape in order to exclude PDMSair interaction.
Compound recovery from PDMS is achieved either by thermal desorption and on-line transferred to the injector of a Gas Chromatography-Mass Spectrometry (GC-MS) system, or by liquid extraction with polar solvents. DC-STE can be used successfully for both qualitative and quantitative analyses [37] making DC-STE coupled with GC-MS an efficient approach to characterize the profile and dynamics of PV production in response to both biotic and abiotic stresses.
In this study, the use of DC-STE combined with GC-MS was applied in vivo for the first time to evaluate the dynamics of HIPV release, upon abiotic (MD) and biotic (herbivory wounding, HW) stresses, by using the model system S. littoralis/P. lunatus. Furthermore, MD was used in combination with S. littoralis OS (MDOS). Here we show that HIPVs are differentially produced in different parts of the wounded leaf, depending on the biotic or abiotic stress applied. The analytical method was compared to the expression of genes involved in HIPV biosynthesis, which showed the same HIPV topographical pattern.

Results
In response to herbivory, plants produce PVs, which can serve as direct deterrents [21] or to attract the herbivore's predators and parasitoids [23][24][25][26]39]. The dynamics of HIPV emission, and the relationship between damaged area and biosynthesis of bioactive molecules, remain open questions. An innovative in vivo strategy was here used to identify compounds actively related to plant-insect interactions, employing a non-destructive high concentration-capacity sampling technique to capture volatiles from lima bean leaves after abiotic and biotic wounding.

DC-STE-GC-MS analysis discriminates herbivory from mechanical wounding
To analyze the topographical distribution of HIPVs, leaves from plants grown in a growth chamber treated with HW, MD and MDOS as well as control intact leaves were sampled with PDMS rectangular tapes (4 × 15 × 0.2 mm) placed in direct contact with leaves at specific distances from the damaged areas (0 cm, 1.5 cm, 3 cm) for different sampling times (2, 6, 24 h). Adaxial and abaxial leaf laminae were sampled in three different leaf portions: a) close to the damaged area (referred as the wounding zone, 0 cm); b) in the central portion (referred as the middle zone, 1.5 cm); and c) in the basal portion of the leaf (referred as the basal zone, 3 cm) ( Figure 1). Preliminary trials showed no significant differences in PV results between adaxial and abaxial epidermises (data not shown). Analysis of camphor variation supports the repeatability of the method, accounting for 18.3% as relative standard deviation throughout the whole dataset.
Several PVs were identified by GC-MS analyses including green-leaf volatiles (GLVs, including aldehydes, alcohols and acetates), alkyl aldehydes, homoterpenes, monoand sesquiterpenoids (Additional file 1). Because of the large number of samples (337), several Principal Component Analyses (PCA) were carried out; the best results were those obtained with logarithmic scaling as data pre-treatment [40].  Figure 2B). These HIPVs were therefore used as variables for the subsequent PCA to explore the internal differences in the damaged leaf dataset. A better discrimination (about 71% of total variance explained) was obtained between HW and MD treatments, whereas MDOS samples showed a scattered pattern ( Figure 3).

DC-STE-GC-MS determines and quantifies the topography of leaf HIPV production
The ability to discriminate between MD and HW highlights the potential of DC-STE-GC-MS as a reliable technique for in vivo HIPV monitoring. This ability was used to study the dynamics of volatile production as a function of topography in lima bean leaf responses to HW, MD and MDOS.
To visualize HIPV distribution, the damaged leaf dataset was divided into three different matrices, depending on the type of damage, each including a smaller but still    In MD treated leaves, there was a clear distinction between the wounded area and the rest of the leaf ( Figure 4E). However, only GLVs and 1-octen-3-ol were produced in wounded areas, while (E)-β-ocimene and DMNT were not discriminant for the different leaf portions ( Figure 4F). The observation of the temporal differences showed that in all treatments GLVs were always produced early, whereas production of terpenoids and homoterpenes occurred later. In particular, PCAs highlighted some interesting differences in the temporal patterns between HW and other damages, with MDOS again showing intermediate values (Additional file 2).
A quantitative evaluation of the main damage-related compounds were carried out by combining in-tape camphor standardization with an external calibration by Gas Chromatography -Selected Ion Monitoring -Mass Spectrometry (GC-SIM-MS) for all types of damage, reaching a good linearity for every quantified HIPV (for quantitation parameters see Additional file 3).
In general, GLVs were the most abundant compounds in the damaged area (Table 1). (Z)-3-hexen-1-ol, (E)-2hexenal, and 1-octen-3-ol reach rates of up to 100 ng/cm 2 . (E)-β-ocimene and (E)-nerolidol were generally produced in smaller amounts far from the wounded zone; however, they were found to exceed 100 ng/cm 2 in the damage area. The homoterpenes, DMNT and TMTT, were mostly found in low quantities in HW-damaged leaves (Table 1).

Topographical gene expression analysis and DC-STE-GC-MS HIPV mapping
Because of the non-destructive DC-STE method of PV sampling, the different leaf sampled portions producing HIPVs could be used for gene expression analyses. Farnesyl Pyrophosphate Synthase (FPS) [41], P. lunatus Ocimene Synthase (PlOS) [10] and P. lunatus Terpene Synthase 2 (PlTPS2) [42] gene expressions were analyzed and compared to the results obtained by DC-STE for the related compounds . In addition, Lipoxygenase (LOX) [41] gene expression was analyzed, to assess any similarity with the observed high formation of GLVs.
Significantly higher expression of PlOS ( Figure 5A) was in all cases coherent with the measured amount of the related compound (E)-β-ocimene ( Figure 5B), with fold change values > 10 in the wounded zones of leaves treated by HW, MD or MDOS. Production of the homoterpene TMTT was associated with the gene expression pattern of PlTPS2 only for HW and MDOS treatments, whereas regulation of the gene was not comparable to the amount of the homoterpene upon MD treatment ( Figure 5: C and D). Upregulation of FPS gene expression ( Figure 5E) was consistent with (E)-nerolidol amount in HW and MDOS treatments ( Figure 5F). Finally, the total GLV -production ( Figure 5H) was in all cases higher in wounded zones, and consistent with LOX upregulation, in particular when referred to HW ( Figure 5G). These results are fully supported by the Kruskal-Wallis significance test (with Bonferroni adjustment, p < 0.017), as shown in Figure 5.

Discussion
One of the most challenging tasks in multitrophic interaction studies is the adoption of advanced analytical platforms that enable different analyses to be run simultaneously using different "omic" methodologies. DC-STE-GC-MS enabled to characterize the qualitative and quantitative topographical profile of leaf volatile emission upon herbivory, while evaluating at the same time the gene expression of the same sampled tissues.
The present results highlight the key role of the damaged area in HIPV production [35], with GLVs associated with both mechanical damage and herbivory, and monoterpenoids, sesquiterpenoids, and homoterpenes specifically associated with herbivory. In particular, MD treatment appears to be sufficient to induce higher amount of GLVs, including (E)-2-hexenal, (Z)-3-hexen-1ol, 1-octen-3-ol, (Z)-3-hexenyl acetate and (Z)-3-hexenyl butyrate [20,25,32,[45][46][47][48]. The DC-STE-GC-MS technique enabled GLVs to be determined qualitatively and to be quantified for further comparisons. Furthermore, the analysis revealed that some GLVs [(E)-2-hexenal, (Z)-3hexen-1-ol and 1-octen-3-ol] are more intensively produced during MD than they are during other stresses. Conversely, (Z)-3-hexenyl acetate and (Z)-3-hexenyl butyrate are produced in higher amount in HW leaves, supporting their herbivory-induced production pattern [23,25,26,49]. GLVs are synthesized via the LOX pathway from C 18 polyunsaturated fatty acids [50], which are cleaved to C 12 and C 6 compounds by hydroperoxide lyases (HPL) [28]. Most plants have several isoforms of LOX [51], and a specific LOX that is essential to GLV formation has been identified in a few plant species [52]. In the present study, upregulation of LOX expression was evenly distributed throughout the leaf, although GLVs were mostly found in the wounded area. This discrepancy between gene expression and GLV production may be due, on the one hand, to the wide variety of roles played by LOX [53], and, on the other hand, to the effect of the GLVs on leaf tissues [30]. For instance, (Z)-3-hexenal in the vapor phase was taken up by Arabidopsis and converted into its alcohol and acetate in the cells. This scenario was further confirmed by the fact that the isotope ratios of alcohol and acetate were almost identical to that of (Z)-3-hexenal when 13 C-labeled (Z)-3-hexenal of a given isotope ratio was used for the exposure [54]. GLVs produced in the wounded zone may therefore influence expression of genes in unwounded tissues of the same leaf  [30], DC-STE sampling also highlights that GLVs play a central role in the so-called plant damage associated molecular pattern (DAMP). Indeed they seem to be essential to trigger gene expression required to prepare an adequate damage reaction in the surrounding tissues and organs. The MD related high GLV production could be explained with their well-known anti-microbial activity [32,55]. This is a resistance trait that is required during pathogen infection, which could occur after wounding [30], even without the herbivore interaction. Among monoterpenes, (E)-β-ocimene, a well-known damage-related HIPV [5,10,44] is a significant example of HIPV distribution. Its amount is limited to the damaged area in HW, while in MDOS (E)-β-ocimene also occurs distant from the wounded tissues. This different distribution agrees with the pattern of PlOS expression, demonstrating to produce almost exclusively (E)-β-ocimene when activated [56]; production is mainly located in the wounding area [45]. Transgenic Arabidopsis, transformed with the PlOS promoter GUS fusion constructs, shows that the activity is restricted to the wounded sites [10]. Lepidopteran caterpillars continuously remove leaf tissue after every bite, even if in a time longer than that one needed for the induction [57]. Conversely, application of OS to MD enables the elicitor to remain on the leaf longer, at least throughout the sampling time. This might explain why, in MDOS treated leaves, PlOS upregulation was observed in leaf areas distant from the damage.
Homoterpenes and sesquiterpenoids, such as DMNT, TMTT and (E)-nerolidol, are often associated with damage-related emission [5,22,26,58]; they have been studied as indirect defense mediators [25,39]. DMNT distribution is comparable to that of TMTT, and shows a general distribution from the damage zone throughout the leaf. However, their amount is higher in the wounded zone after both HW and MDOS treatment. The TPS enzymes have been found to be involved in DMNT and TMTT precursor production [42,58,59] and their products have been related to herbivory events [10,35,58]. The PlTPS2 gene analyzed here showed a distribution comparable to that of homoterpene amount, in particular in leaves undergoing HW and MDOS.
Production of (E)-nerolidol is limited to the wounded zone, in particular in HW and MDOS damage, while MD does not seem to induce it. The lower amount of this compound in MDOS compared to HW is of interest because it shows the inability of OS alone to trigger the same HW-related leaf emission. Expression of FPS was found to be upregulated not only in damaged areas but also in leaf tissues distant from the wounding zone. FPS plays a key role in HIPV emission since its product, farnesyl pyrophosphate, is a basic precursor for sesquiterpenoid biosynthesis [13,60]. FPS is considered an important HW-related enzyme [43] and its inducibility by HIPVs has also been discussed and confirmed [61,62]. FPS upregulation was marked in HW leaves, underlining the relationship between herbivory and FPS activation [43].

Conclusions
The use of DC-STE-GC-MS provides a clearer picture of DAMP distribution in lima bean, by showing differential release of HIPV classes after different kinds of wounding. DAMPs, which are essential for airborne damage-signals, were found to be mainly related to disrupted tissues. The results confirm the role of HIPVs as DAMP signals and show their role as signals able to quickly spread in the surrounding environment of wounded areas. Upon herbivory a fast V m depolarization is known to affect the whole damaged leaf, whereas calcium, potassium, ROS and NO responses are limited to the wounded zones. DC-STE-GC-MS results show that GLVs are released almost immediately and their emission is topographically in concomitance with early events such as V m depolarization and calcium signaling, as previous data suggested [32,[63][64][65][66][67][68][69]. The

Plant and animal material
Feeding experiments were carried out using the lima bean (Phaseolus lunatus L. cv Ferry Morse var. Jackson Wonder Bush). Individual plants were grown from seed in plastic pots with quartz sand at 23°C and 60% humidity, using daylight fluorescent tubes at approximately 270 μE m −2 s −1 with a photophase of 16 h. Experiments were conducted with 12-to 16-day-old seedlings showing two fully-developed primary leaves, which were found to be the most responsive [1].
Spodoptera littoralis Boisd. (Lepidoptera, Noctuidae) larvae were kindly provided by R. Reist from Syngenta Crop. Protection Münchwilen AG, Switzerland, and were fed on an artificial diet comprising 125 g bean flour, 2.25 g ascorbic acid, 2.25 g ethyl 4-hydroxybenzoate, 750 μL formaldehyde, 300 mL distilled water and 20 g agar, previously solubilized in 300 mL of distilled water. The ingredients (Sigma-Aldrich, St. Louis, MO, USA) were mixed with a blender and stored at 4°C for not more than one week. With the exception of plant volatile (PV) collection (see below), plants were exposed for 2 h to third instar larvae reared from egg clutches in Petri dishes (9 cm diameter) in a growth chamber with 16 h photoperiod at 25°C and 60-70% humidity. The amount of herbivore damage was limited to 30% of leaf surface, as detected by ImageJ image analysis [4]. Feeding experiments were always performed between 1 and 3 p.m.

Collection of oral secretions
In order to evaluate the effect of S. littoralis oral secretions (OS), 5-day-old larvae were allowed to feed on lima bean leaves for 24 h. Regurgitation was caused by gently squeezing the larva with a forceps behind the head. OS was collected in glass capillaries connected to an evacuated sterile vial (peristaltic pump).

PV sampling setup
Biotic stress was caused by S. littoralis (HW); whereas abiotic stress was performed by mechanically damaging leaf tissues with a pattern wheel (MD). Furthermore abiotic and biotic stresses were connected by combining MD with S. littoralis oral secretions (MDOS). A large number of samples were analyzed (337) and multivariate methods were used to define discriminant variables (i.e., HIPVs) and to plot chemical and molecular topographical maps of leaf areas producing HIPVs in response to biotic and abiotic stress. In particular, the experiments were carried out in nine sampling steps, each representing a specific combination of type of damage (HW, MD, and MDOS) and sampling duration (2, 6, 24 h). For each sampling step, three biological replicates were analyzed, with 12 tapes for each. A control using two tapes was also sampled. HW was caused by S. littoralis caterpillars; the damaged area for each plant was as near as possible equal. MD was done by piercing the leaves manually with a pattern wheel. The damaged leaf area and the duration of time of the damaging mechanism were kept constant. The damage process in MDOS was similar to that in MD, with the addition on the wounded area of 10 μL of a solution 1:1 of S. littoralis OS and 5 mM MES (2-(N-morpholino)-ethane-sulphonic acid) buffer (pH 6.0). The OS quantity was assessed after several trials (from 0.5 to 10 μL) and was found the most appropriate to obtain reproducible experiments [43].
At the end of the sampling time, the tapes were removed and stored at −20°C. Leaves were cut into 3 parts (wounded area, middle, base) and stored at −80°C for further analyses.

Direct Contact-Sorptive Tape Extraction of PVs
Polydimethylsiloxane (PDMS) tapes (4 × 15 × 0.2 mm, ca. 33 mg) were placed on different areas of the adaxial and abaxial leaf lamina of S. littoralis-attacked and of control leaves. A glass coverslip was placed just above the DC-STE tape in order to exclude PDMSair interaction. The quantitation of the collected PVs was obtained by an external standard at known concentration levels, being difficult to calculate an analyte recovery rate with DC-STE applied to in vivo plant matrices (unlike it was done in [37] with standards). Sampling was carried out in triplicate in the positions on the leaf shown in Figure 1, for the times reported above (2, 6, 24 h). Camphor (Sigma-Aldrich, Milan, Italy) was used as internal standard (I.S.) and was sorbed onto the tapes as proposed by Wang et al. [70] for Solid Phase Micro Extraction. Preliminary analysis with tapes with and without camphor I.S. were carried out to verify any possible interference of camphor with lima bean PV production (Additional file 4). After sampling, the PDMS tapes were placed in thermal desorption tubes, stored in sealed vials, and submitted to automatic thermal desorption (see below). Sorption tapes were provided by the Research Institute for Chromatography (Kortrijk-Belgium).

GC-MS analysis
PDMS tape thermal desorption was carried out with a Thermal Desorption Unit (TDU) from Gerstel (Mülheima/ d Ruhr, Germany). Analyses were driven automatically by an MPS-2 multipurpose sampler installed on an Agilent 7890 GC unit coupled to an Agilent 5975C MSD (Agilent, Little Falls, DE, USA). The TDU thermal desorption program was: from 30°C to 250°C (5 min) at 60°C/min in splitless flow mode, and transfer line at 300°C. A Gerstel CIS-4 PTV injector was used to cryofocus compounds thermally desorbed from the PDMS tapes, and inject them into the injector GC port. The PTV was cooled to −40°C using liquid CO 2 ; injection temperature: from −40°C to 250°C (5 min) at 12°C/s. The inlet was operating in the splitless mode. Helium was used as carrier gas at a flow rate of

GC-MS data processing
Data were processed with Agilent MSD ChemStation ver. D.03.00.611 (Agilent Technologies). Components were identified by comparing their linear retention indices (I T s) (calculated versus a C 9 -C 25 hydrocarbon mixture) and their mass spectra to those of authentic samples, or by comparison with those present in commercially-available mass spectrum libraries (Wiley, Adams).

RNA extraction from lima bean leaves after HW, MD and MDOS
After each experiment, leaves were collected and immediately frozen in liquid nitrogen. Samples from time-course experiments were pooled so as to have a single pool of replicates for each stress condition (HW, MDOS, MD, undamaged leaves). Fifty mg of frozen leaf material were ground in liquid nitrogen with mortar and pestle. Total RNA was isolated using the Agilent Plant RNA Isolation Mini Kit (Agilent Technologies, Santa Clara, CA, US) and RNase-Free DNase set (Qiagen, Hilden, Germany). Sample quality and quantity were checked using the RNA 6000 Nano kit and the Agilent 2100 Bioanalyzer (Agilent Technologies), following the manufacturer's instructions. Quantification of RNA was also confirmed spectrophotometrically, using the NanoDrop ND-1000 (Thermo Fisher Scientific, Waltham, MA, US).

Quantitative real time-PCR (qPCR) reaction conditions and primers
First strand cDNA synthesis was run with 1 μg of total RNA and random primers, using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, US), and following the manufacturer's recommendations. Reactions were prepared by adding 1 μg of total RNA, 2 μL of 10X RT Buffer, 0.8 μL of 25X dNTPs mix (100 mM), 2 μL 10X RT random primer, 1 μL of Multiscribe™ Reverse Transcriptase, and nuclease-free sterile water to 20 μL. Reaction mixtures were incubated at 25°C for 10 min, 37°C for 2 h, and 85°C for 5 s.
The qPCR experiments were run on a Stratagene Mx3000P Real-Time System (La Jolla, CA, USA) using SYBR green I with ROX as an internal loading standard. The reaction mixture was 10 μL, comprising 5 μL of 2X Maxima™ SYBR Green qPCR Master Mix (Fermentas International, Inc, Burlington, ON, Canada), 0.5 μL of cDNA and 100 nM primers (Integrated DNA Technologies, Coralville, IA, US). Controls included non-RT controls (using total RNA without reverse transcription to monitor for genomic DNA contamination) and nontemplate controls (water template). Specifically, PCR conditions were the following: P. lunatus Actin1 (PlACT1), Farnesyl Pyrophosphate Synthase (FPS), Lipoxygenase (LOX) [41], P. lunatus Ocimene Synthase (PlOS) [10], P. lunatus Terpene Synthase 2 (PlTPS2) [42]: 10 min at 95°C, 45 cycles of 15 s at 95°C, 30 s at 55°C, and 30 s at 72°C, 1 min at 95°C, 30 s at 55°C, 30 s at 95°C. Fluorescence was read following each annealing and extension phase. All runs were followed by a melting curve analysis from 55°C to 95°C. The linear range of template concentration to threshold cycle value (Ct value) was determined by preparing a dilution series, using cDNA from three independent RNA extractions analyzed in three technical replicates. Primer efficiencies for all primer pairs were calculated using the standard curve method [71]. Two different reference genes (Actin1 (PlACT1) and the 18S ribosomal RNA) were used to normalize the results of the qPCR. The best of the two genes was selected using the Normfinder software [72]; the most stable gene was P1ACT1. Primers used for qPCR were as described elsewhere [3,41,42] and are reported in Additional file 5.
All amplification plots were analyzed with the Mx3000P™ software to obtain Ct values. Relative RNA levels were calibrated and normalized with the level of PlACT1 mRNA.

Statistical analyses
Analysis of variance (ANOVA) and the Tukey test were used to assess difference between treatments and control. For all other experiments, at least five samples per treatment group entered the statistical data analysis. PV chemical data are expressed as mean values ± standard error of the mean (SEM).
Principal Component Analysis (PCA) was used in three different steps, each targeting different discrimination (control-damage, different damage, different leaf areas). A log-transformation was used as GC-MS data pre-treatment [40]. Each PCA step was followed by a significance test for the discriminant compounds. To compare the different leaf areas, the Kruskal-Wallis test was applied to both chemical and gene expression data. Bonferroni adjustment (p/k; k = number of comparisons) was applied to protect against Type I Error [73,74].
All statistical data analyses were done using SPSS software for Windows.

Additional files
Additional file 1: PV DC-STE-GC-MS profile. The typical GC-MS profile of PVs captured by DC-STE on Phaseolus lunatus leaves wounded by Spodoptera littoralis obtained after thermal desorption of DC-STE. A compound table is also provided.
Additional file 2: PCA analysis of time-course experiments on different damage dataset. Score and loading plots presented in Figure 4 are here displayed taking into account sampling time of every sample (instead considering the distance from the wounded area).

Competing interests
The authors declare that they have no competing interests.