- Methodology article
- Open Access
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.
BMC Plant Biology volume 15, Article number: 102 (2015)
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.
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-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 , expression of late defense response genes , and emission of herbivory-induced plant volatiles (HIPVs) [12,13].
Even if robotic mechanical wounding can simulate plant response similar to HIPV , 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) . 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-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-26], damage and disease long-distance signaling [27-30], and pathogen resistance priming [29-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 .
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.  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 , 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 PDMS – air 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  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.
In response to herbivory, plants produce PVs, which can serve as direct deterrents  or to attract the herbivore’s predators and parasitoids [23-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, mono- and 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 .
Figure 2A reports the PCA (42% of explained variance) on the total dataset of samples, discriminating undamaged (controls) from damaged leaves. The damaged sample distribution in Figure 2A showed that HW and MD seemed divided into two different subsets, while application of OS to MD leaves produced intermediate results between them.
The resulting damage-related discriminant compounds included GLVs [(E)-2-hexenal, (Z)-3-hexen-1-ol, (Z)-3-hexenyl acetate, (Z)-3-hexenyl butyrate], a linoleic acid breakdown product (1-octen-3-ol), a monoterpene [(E)-β-ocimene], two homoterpenes [4,8-dimethyl-1,3,7-nonatriene (DMNT) and 4,8,12-trimethyl-1,3,7,11-tridecatetraene (TMTT)] and a sesquiterpenoid [(E)-nerolidol] (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 considerable number of samples (HW: 54 samples; MD: 53 samples; MDOS: 52 samples). PCA data processing was performed by using the discriminating variables identified above (GLVs, homoterpenes, mono- and sesquiterpenoids) with the aim of establishing a relationship between sampling time and leaf portion. A distinctive distribution of volatiles as a function of the damaged area was found for both HW and MDOS (Figure 4: A and C). Compared to controls, HW treated leaves showed a significantly higher production of the GLVs (E)-2-hexenal, (Z)-3-hexen-1-ol, (Z)-3-hexenyl acetate and of 1-octen-3-ol close to the HW damaged leaf portion (Figure 4B). (Z)-3-Hexenyl butyrate, DMNT, TMTT, (E)-β-ocimene and (E)-nerolidol were produced in the same area, close to the HW zone, but also in distant leaf portions (Figure 4B). A similar pattern was found when MD plants were treated with OS (Figure 4: C and D). 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)-2-hexenal, and 1-octen-3-ol reach rates of up to 100 ng/cm2. (E)-β-ocimene and (E)-nerolidol were generally produced in smaller amounts far from the wounded zone; however, they were found to exceed 100 ng/cm2 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) , P. lunatus Ocimene Synthase (PlOS)  and P. lunatus Terpene Synthase 2 (PlTPS2)  gene expressions were analyzed and compared to the results obtained by DC-STE for the related compounds . In addition, Lipoxygenase (LOX)  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.
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 , 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-1-ol, 1-octen-3-ol, (Z)-3-hexenyl acetate and (Z)-3-hexenyl butyrate [20,25,32,45-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)-3-hexen-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 C18 polyunsaturated fatty acids , which are cleaved to C12 and C6 compounds by hydroperoxide lyases (HPL) . Most plants have several isoforms of LOX , and a specific LOX that is essential to GLV formation has been identified in a few plant species . 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 , and, on the other hand, to the effect of the GLVs on leaf tissues . 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 13C-labeled (Z)-3-hexenal of a given isotope ratio was used for the exposure . GLVs produced in the wounded zone may therefore influence expression of genes in unwounded tissues of the same leaf because of GLV diffusion. In line with what Heil and Land have been recently reported , 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 , 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 ; production is mainly located in the wounding area . Transgenic Arabidopsis, transformed with the PlOS promoter GUS fusion constructs, shows that the activity is restricted to the wounded sites . Lepidopteran caterpillars continuously remove leaf tissue after every bite, even if in a time longer than that one needed for the induction . 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  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 .
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-69].
The DC-STE-GC-MS results are in agreement with the present body of knowledge of plant damage recognition and reaction, and provide a better understanding of the dynamics of plant responses to damage. The main advantages of this technique compared to classical PV sampling methods are: a) in vivo sampling; b) ease of execution; c) simultaneous assays of different leaf portions, and d) preservation of plant material for further omic studies.
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 .
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 . 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 .
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 PDMS – air 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  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.  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).
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 CO2; 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 1 mL/min. Column: HP5MS (30 m × 0.25 mm i.d. × 0.25 μm; Agilent Technologies). Temperature program: from −30°C (1 min) to 50°C at 50°C/min, then to 165°C at 3°C/min, then to 250°C (5 min) at 25°C/min. MS operated in EI mode at 70 eV with a mass range from 35 to 350 amu in full scan mode.
Quantitative Gas Chromatography – Selected Ion Monitoring – Mass Spectrometry analysis (GC-SIM-MS): appropriate amounts of 2-hexenal, 3-hexenol, 1-octen-3-ol, 3-hexenyl acetate, (Z)-3-hexenyl butyrate, 1-octen-3-ol, (E)-β-ocimene, 4,8-dimethyl-1,3,7-nonatriene and 4,8,12-trimethyl-1,3,7,11-tridecatetraene (Sigma-Aldrich, Milan, Italy) were diluted with cyclohexane (Sigma-Aldrich, Milan, Italy) to obtain nine different concentrations in the range 1 to 1000 μg/mL for each component. Calibration curves were constructed by analyzing the resulting standard solutions three times, by GC-MS in SIM mode, under the conditions reported above.
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 Ts) (calculated versus a C9-C25 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 non-template controls (water template). Specifically, PCR conditions were the following: P. lunatus Actin1 (PlACT1), Farnesyl Pyrophosphate Synthase (FPS), Lipoxygenase (LOX) , P. lunatus Ocimene Synthase (PlOS) , P. lunatus Terpene Synthase 2 (PlTPS2) : 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 . 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 ; 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.
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 . 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.
Damage associated molecular pattern
Direct contact-sorptive tape extraction
Green leaf volatile
Herbivory-induced plant volatile
Mechanical damage with oral secretions
Maffei ME, Bossi S, Spiteller D, Mithöfer A, Boland W. Effects of feeding Spodoptera littoralis on lima bean leaves. I. Membrane potentials, intracellular calcium variations, oral secretions, and regurgitate components. Plant Physiol. 2004;134(April):1752–62.
Maffei ME, Mithöfer A, Boland W. Before gene expression: early events in plant-insect interaction. Trends Plant Sci. 2007;12:310–6.
Maffei ME, Mithöfer A, Arimura G, Uchtenhagen H, Bossi S, Bertea CM, et al. Effects of feeding Spodoptera littoralis on lima bean leaves. III. Membrane depolarization and involvement of hydrogen peroxide. Plant Physiol. 2006;140:1022–35.
Bricchi I, Leitner M, Foti M, Mithöfer A, Boland W, Maffei ME. Robotic mechanical wounding (MecWorm) versus herbivore-induced responses: early signaling and volatile emission in Lima bean (Phaseolus lunatus L.). Planta. 2010;232:719–29.
Mithöfer A, Wanner G, Boland W. Effects of Feeding Spodoptera littoralis on Lima Bean Leaves II. Continuous Mechanical Wounding Resembling Insect Feeding Is Sufficient to Elicit Herbivory-Related Volatile Emission. Plant Physiol. 2005;137(March):1160–8.
Bonaventure G, VanDoorn A, Baldwin IT. Herbivore-associated elicitors: FAC signaling and metabolism. Trends Plant Sci. 2011;16:294–9.
Maffei ME, Arimura G-I, Mithöfer A. Natural elicitors, effectors and modulators of plant responses. Nat Prod Rep. 2012;29:1288–303.
Arimura G-I, Maffei ME. Calcium and secondary CPK signaling in plants in response to herbivore attack. Biochem Biophys Res Commun. 2010;400:455–60.
Arimura G-I, Ozawa R, Maffei ME. Recent advances in plant early signaling in response to herbivory. Int J Mol Sci. 2011;12:3723–39.
Arimura G, Köpke S, Kunert M, Volpe V, David A, Brand P, et al. Effects of feeding Spodoptera littoralis on lima bean leaves: IV. Diurnal and nocturnal damage differentially initiate plant volatile emission. Plant Physiol. 2008;146:965–73.
Wu J, Baldwin IT. New insights into plant responses to the attack from insect herbivores. Annu Rev Genet. 2010;44:1–24.
Baldwin IT. Plant volatiles. Curr Biol. 2010;20:R392–7.
Maffei ME, Gertsch J, Appendino G. Plant volatiles: production, function and pharmacology. Nat Prod Rep. 2011;28:1359–80.
Zebelo S, Piorkowski J, Disi J, Fadamiro H. Secretions from the ventral eversible gland of Spodoptera exigua caterpillars activate defense-related genes and induce emission of volatile organic compounds in tomato Solanum lycopersicum. BMC Plant Biol. 2014;14:140.
Mori N, Yoshinaga N. Function and evolutionary diversity of fatty acid amino acid conjugates in insects. J Plant Interact. 2011;6:103–7.
Zebelo SA, Maffei ME. The ventral eversible gland (VEG) of Spodoptera littoralis triggers early responses to herbivory in Arabidopsis thaliana. Arthropod Plant Interact. 2012;6:543–51.
Van Poecke RMP, Dicke M. Indirect defence of plants against herbivores: using Arabidopsis thaliana as a model plant. Plant Biol (Stuttg). 2004;6:387–401.
Leitner M, Boland W, Mithöfer A. Direct and indirect defences induced by piercing-sucking and chewing herbivores in Medicago truncatula. New Phytol. 2005;167:597–606.
Turlings TCJ, Ton J. Exploiting scents of distress: the prospect of manipulating herbivore-induced plant odours to enhance the control of agricultural pests. Curr Opin Plant Biol. 2006;9:421–7.
Maffei ME. Sites of synthesis, biochemistry and functional role of plant volatiles. South African J Bot. 2010;76:612–31.
Zebelo SA, Bertea CM, Bossi S, Occhipinti A, Gnavi G, Maffei ME. Chrysolina herbacea modulates terpenoid biosynthesis of Mentha aquatica L. PLoS One. 2011;6, e17195.
Mithöfer A, Boland W. Plant defense against herbivores: chemical aspects. Annu Rev Plant Biol. 2012;63:431–50.
De Moraes CM, Lewis WJ, Paré PW, Alborn HT, Tumlinson JH. Herbivore-infested plants selectively attract parasitoids. Nature. 1998;393(June):570–3.
Ponzio C, Gols R, Pieterse CMJ, Dicke M. Ecological and phytohormonal aspects of plant volatile emission in response to single and dual infestations with herbivores and phytopathogens. Funct Ecol. 2013;27:587–98.
Mäntylä E, Alessio GA, Blande JD, Heijari J, Holopainen JK. From plants to birds: higher avian predation rates in trees responding to insect herbivory. PLoS One. 2008;3, e2832.
Clavijo McCormick A, Irmisch S, Reinecke A, Boeckler GA, Veit D, Reichelt M, et al. Herbivore-induced volatile emission in black poplar: regulation and role in attracting herbivore enemies. Plant Cell Environ. 2014;37(8):1909–23.
Dicke M, Baldwin IT. The evolutionary context for herbivore-induced plant volatiles: beyond the “cry for help”. Trends Plant Sci. 2010;15:167–75.
Engelberth J, Alborn HT, Schmelz EA, Tumlinson JH. Airborne signals prime plants against insect herbivore attack. Proc Natl Acad Sci U S A. 2004;101:1781–5.
Yi H-S, Heil M, Adame-Alvarez RM, Ballhorn DJ, Ryu C-M. Airborne induction and priming of plant defenses against a bacterial pathogen. Plant Physiol. 2009;151:2152–61.
Heil M, Land WG. Danger signals - damaged-self recognition across the tree of life. Front Plant Sci. 2014;5(October):578.
Ballhorn DJ, Kautz S, Schädler M. Induced plant defense via volatile production is dependent on rhizobial symbiosis. Oecologia. 2013;172:833–46.
Scala A, Allmann S, Mirabella R, Haring MA, Schuurink RC. Green leaf volatiles: a plant’s multifunctional weapon against herbivores and pathogens. Int J Mol Sci. 2013;14:17781–811.
Tholl D, Boland W, Hansel A, Loreto F, Röse USR, Schnitzler J-P. Practical approaches to plant volatile analysis. Plant J. 2006;45:540–60.
Bicchi C, Maffei M. The Plant Volatilome: Methods of Analysis. Methods Mol Biol. 2012;918:289–310.
Köllner TG, Lenk C, Schnee C, Köpke S, Lindemann P, Gershenzon J, et al. Localization of sesquiterpene formation and emission in maize leaves after herbivore damage. BMC Plant Biol. 2013;13:15.
Sandra P, Sisalli S, Adao A, Lebel M, Le Fur I. Sorptive Tape Extraction — A Novel Sampling Method for the in vivo Study of Skin. LCGC Eur. 2006;19:33–9.
Sgorbini B, Ruosi MR, Cordero C, Liberto E, Rubiolo P, Bicchi C. Quantitative determination of some volatile suspected allergens in cosmetic creams spread on skin by direct contact sorptive tape extraction-gas chromatography–mass spectrometry. J Chromatogr A. 2010;1217:2599–605.
Bicchi C, Cordero C, Liberto E, Rubiolo P, Sgorbini B, Sandra P. Sorptive tape extraction in the analysis of the volatile fraction emitted from biological solid matrices. J Chromatogr A. 2007;1148:137–44.
Heil M. Indirect defence via tritrophic interactions. New Phytol. 2008;178:41–61.
Brereton RG. Chemometrics for Pattern Recognition. Chichester, UK: John Wiley & Sons, Ltd; 2009.
Ozawa R, Bertea CM, Foti M, Narayana R, Arimura G-I, Muroi A, et al. Exogenous polyamines elicit herbivore-induced volatiles in lima bean leaves: involvement of calcium, H2O2 and Jasmonic acid. Plant Cell Physiol. 2009;50:2183–99.
Brillada C, Nishihara M, Shimoda T, Garms S, Boland W, Maffei ME, et al. Metabolic engineering of the C16 homoterpene TMTT in Lotus japonicus through overexpression of (E, E)-geranyllinalool synthase attracts generalist and specialist predators in different manners. New Phytol. 2013;200:1200–11.
Bricchi I, Occhipinti A, Bertea CM, Zebelo SA, Brillada C, Verrillo F, et al. Separation of early and late responses to herbivory in Arabidopsis by changing plasmodesmal function. Plant J. 2013;73:14–25.
Zhang P-J, Zheng S-J, Van Loon JJ, Boland W, David A, Mumm R, et al. Whiteflies interfere with indirect plant defense against spider mites in Lima bean. Proc Natl Acad Sci U S A. 2009;106:21202–7.
Arimura G-I, Matsui K, Takabayashi J. Chemical and molecular ecology of herbivore-induced plant volatiles: proximate factors and their ultimate functions. Plant Cell Physiol. 2009;50:911–23.
Heil M, Lion U, Boland W. Defense-inducing volatiles: in search of the active motif. J Chem Ecol. 2008;34:601–4.
Holopainen JK, Gershenzon J. Multiple stress factors and the emission of plant VOCs. Trends Plant Sci. 2010;15:176–84.
De Boer J, Snoeren T, Dicke M. Predatory mites learn to discriminate between plant volatiles induced by prey and nonprey herbivores. Anim Behav. 2005;69:869–79.
D’Auria JC, Pichersky E, Schaub A, Hansel A, Gershenzon J. Characterization of a BAHD acyltransferase responsible for producing the green leaf volatile (Z)-3-hexen-1-yl acetate in Arabidopsis thaliana. Plant J. 2007;49:194–207.
Grechkin A. Recent developments in biochemistry of the plant lipoxygenase pathway. Prog Lipid Res. 1998;37:317–52.
Joo Y-C, Oh D-K. Lipoxygenases: potential starting biocatalysts for the synthesis of signaling compounds. Biotechnol Adv. 2012;30:1524–32.
Christensen SA, Nemchenko A, Borrego E, Murray I, Sobhy IS, Bosak L, et al. The maize lipoxygenase, ZmLOX10, mediates green leaf volatile, jasmonate and herbivore-induced plant volatile production for defense against insect attack. Plant J. 2013;74:59–73.
Liavonchanka A, Feussner I. Lipoxygenases: occurrence, functions and catalysis. J Plant Physiol. 2006;163:348–57.
Matsui K, Sugimoto K, Mano J, Ozawa R, Takabayashi J. Differential metabolisms of green leaf volatiles in injured and intact parts of a wounded leaf meet distinct ecophysiological requirements. PLoS One. 2012;7, e36433.
Heil M. Herbivore-induced plant volatiles: targets, perception and unanswered questions. New Phytol. 2014;204:297–306.
Degenhardt J, Köllner TG, Gershenzon J. Monoterpene and sesquiterpene synthases and the origin of terpene skeletal diversity in plants. Phytochemistry. 2009;70:1621–37.
Schittko U, Preston CA, Baldwin IT. Eating the evidence? Manduca sexta larvae can not disrupt specific jasmonate induction in Nicotiana attenuata by rapid consumption. Planta. 2000;210:343–6.
Tholl D, Sohrabi R, Huh J-H, Lee S. The biochemistry of homoterpenes. Common constituents of floral and herbivore-induced plant volatile bouquets. Phytochemistry. 2011;72:1635–46.
Arimura G, Garms S, Maffei M, Bossi S, Schulze B, Leitner M, et al. Herbivore-induced terpenoid emission in Medicago truncatula: concerted action of jasmonate, ethylene and calcium signaling. Planta. 2008;227:453–64.
Dewick PM. The biosynthesis of C5-C25 terpenoid compounds. Nat Prod Rep. 1997;14:111–44.
Ueda H, Kikuta Y, Matsuda K. Plant communication. Mediated by individual or blended VOCs? Plant Signal Behav. 2012;7:222–6.
Arimura G, Ozawa R, Shimoda T, Nishioka T, Boland W, Takabayashi J. Herbivory-induced volatiles elicit defence genes in lima bean leaves. Nature. 2000;406:512–5.
Arimura G, Ozawa R, Horiuchi J, Nishioka T, Takabayashi J. Plant–plant interactions mediated by volatiles emitted from plants infested by spider mites. Biochem Syst Ecol. 2001;29:1049–61.
Kost C, Heil M. Herbivore-induced plant volatiles induce an indirect defence in neighbouring plants. J Ecol. 2006;94:619–28.
Allmann S, Halitschke R, Schuurink RC, Baldwin IT. Oxylipin channelling in Nicotiana attenuata: lipoxygenase 2 supplies substrates for green leaf volatile production. Plant Cell Environ. 2010;33:2028–40.
Hatanaka A. The biogeneration of green odour by green leaves. Phytochemistry. 1993;34:1201–18.
Farag MA, Paré PW. C6-green leaf volatiles trigger local and systemic VOC emissions in tomato. Phytochemistry. 2002;61:545–54.
Matsui K. Green leaf volatiles: hydroperoxide lyase pathway of oxylipin metabolism. Curr Opin Plant Biol. 2006;9:274–80.
Zebelo SA, Matsui K, Ozawa R, Maffei ME. Plasma membrane potential depolarization and cytosolic calcium flux are early events involved in tomato (Solanum lycopersicon) plant-to-plant communication. Plant Sci. 2012;196:93–100.
Wang Y, O’Reilly J, Chen Y, Pawliszyn J. Equilibrium in-fibre standardisation technique for solid-phase microextraction. J Chromatogr A. 2005;1072:13–7.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29, e45.
Andersen CL, Jensen JL, Ørntoft TF. Normalization of Real-Time Quantitative Reverse Transcription-PCR Data: A Model-Based Variance Estimation Approach to Identify Genes Suited for Normalization, Applied to Bladder and Colon Cancer Data Sets. Cancer Res. 2004;64:5245–50.
Bland JM, Altman DG. Multiple significance tests: the Bonferroni method. Br Med J. 1995;310:170.
Strassburger K, Bretz F. Compatible simultaneous lower confidence bounds for the Holm procedure and other Bonferroni-based closed tests. Stat Med. 2008;27:4914–27.
This work was partly supported by the Doctorate School of Pharmaceutical and Biomolecular Sciences of the University of Turin. The study was also carried out in the framework of the project “Studio di metaboliti secondari biologicamente attivi da matrici di origine vegetale” financially supported by the Ricerca Locale (Ex 60% 2013) of the University of Turin, Turin (Italy).
The authors declare that they have no competing interests.
LB and BS carried out the PV sampling and analyses. LB performed the statistical analysis. CMB and LB carried out the molecular genetic studies. CB, MEM and PR conceived the study, and participated in its design and coordination and helped to draft the manuscript. CC helped to draft the manuscript. All authors read and approved the final manuscript.
LB is a graduate student in the Doctorate School of Pharmaceutical and Biomolecular Sciences of the University of Turin. BS and CC are Assistant Professors of Pharmaceutical Biology in the Dept. of Drug Science and Technology, University of Turin. CB and PR are Professors of Pharmaceutical Biology in the Dept. of Drug Science and Technology, University of Turin. CMB is Associate Professor of Plant Physiology in the Dept. Life Sciences and Systems Biology, University of Turin. MEM is Professor of Plant Physiology in the Dept. Life Sciences and Systems Biology, University of Turin.
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.
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).
Parameters for the volatiles’ quantitation. A table reports HIPV quantitation parameters obtained by Gas Chromatography – Selected Ion Monitoring – Mass Spectrometry (GC-SIM-MS).
Camphor effects on Phaseolus lunatus . A table and a barchart report a comparison between PVs in plants analyzed with and without camphor preloading on tapes.
Gene-specific primers used for quantitative real-time PCR. A table collecting GenBank accession number and sequences of every primer used.
About this article
Cite this article
Boggia, L., Sgorbini, B., Bertea, C.M. et al. 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.. BMC Plant Biol 15, 102 (2015). https://doi.org/10.1186/s12870-015-0487-4
- Direct Contact-Sorptive Tape Extraction (DC-STE)
- Gas Chromatography coupled with Mass Spectrometry (GC-MS)
- Herbivory-induced plant volatile (HIPV)
- Phaseolus lunatus L.
- Spodoptera littoralis Boisd.
- Plant-insect interactions
- Green leaf volatiles (GLVs)