Molecular Interactions between the Specialist Herbivore Manduca sexta (Lepidoptera, Sphigidae) and Its Natural Host Nicotiana attenuata . VIII. An Unbiased GCxGC-ToFMS Analysis of the Plant’s Elicited Volatile Emissions [W][OA]

Treating wounds in Nicotiana attenuata leaves with Manduca sexta oral secretions (W 1 OS) mimics most changes elicited by M. sexta herbivory, but an unbiased analysis of the effect of the different OS constituents on volatile emissions is lacking. We used two-dimensional gas chromatography/time-of-ﬂight (GCxGC-ToF) mass spectrometry combined with multivariate statistics to parse volatiles into regulatory patterns. Volatiles released by wounding alone and by the alkalinity of OS were assessed by applying a buffer known to mimic the pH-mediated changes of OS elicitation (pectin methyl esterase activation and meth- anol release). The activities of fatty acid amino acid conjugates, well-known elicitors of antiherbivore defenses, and of 2-hydroxyoctadecatrienoic acid, a newly discovered signal in OS, were determined. Approximately 400 analytes were detected after deconvolution and alignment of GCxGC data; 35 volatiles were signiﬁcantly regulated upon W 1 OS. Two-thirds of these were speciﬁcally regulated by OS, being either ampliﬁed (most terpenoids and certain hexenylesters) or strongly repressed (many short-chain alcohols and some aromatic and hexenylester derivatives). Fatty acid amino acid conjugates played a central role in this pattern of regulation, since they induced the emission of half of OS-elicited volatiles and inhibited the production of almost all OS-repressed volatiles; 2-hydroxyoctadecatrienoic acid inﬂuenced emission of trans- a -bergamotene, while other unknown OS constituents ampliﬁed hexenylester production. We conclude that the complex bouquet of herbivory-elicited volatiles results from the complex modulations of the wound response by diverse cues found in OS. This work also underscores the value of ultra-high-resolution GCxGC-ToF analysis combined with the nontargeted mining of the resulting data. in intensity across univariate ANOVA ( P , 0.05) on log 2 -transformed normalized peak areas. Hierarchical clustering analyses done autoscale using the Euclidean distance as a clustering metric and the complete linkage aggregation method. The statistical robustness of the dendogram structures was tested by the bootstrap resampling procedure. PCA was performed on log 2 -transformed and mean-centered values using an add-in group Projection plots were obtained from the coordinates for PCs

Plants employ a complex arsenal of defense mechanisms to protect themselves from insect herbivores. In addition to constitutive defenses such as trichomes or thick secondary walls, plants also use induced defenses, which are deployed specifically when herbivores start to chew on a leaf (Karban and Baldwin, 1997). These conditional responses allow plants to forgo the costs of defense production (Zavala et al., 2004) and optimize their ability to compete with neighboring plants or to tolerate the damage of resistant herbivores (Stowe et al., 2000;Schwachtje et al., 2006). Induced defenses function either directly via the production of amino acid-catabolizing enzymes, antidigestive proteinase inhibitors, or toxic and repelling chemicals that render plant tissues less suitable as food for herbivores (Duffey and Stout, 1996) or indirectly through traits that increase the apparency of infested plants to the natural enemies of their herbivores (Karban and Baldwin, 1997;Pare and Tumlinson, 1999).
An intensively studied example of an inducible indirect defense is the production and emission of volatile organic compounds (VOCs) by plants attacked by insects (for review, see Turlings and Wäckers, 2004). The first complete evidence that plants use odors to attract natural enemies of their predators was obtained from studies on plant-mite interactions (Sabelis and van de Baan, 1983;Dicke et al., 1988;Dicke et al., 1990aDicke et al., , 1990b. These first studies and many since, performed under optimized laboratory conditions (Dicke, 1994;Mattiacci et al., 1995Alborn et al., 1997Turlings et al., 1990Turlings et al., , 1991aTurlings et al., , 1991bTurlings et al., , 1998 or under field conditions (Shimoda et al., 1997;De Moraes et al., 1998;Bernasconi Ockroy et al., 2001;Kessler and Baldwin, 2001), have demonstrated that parasitoids and predators make effective use of herbivory-induced plant VOCs (HIPV) to locate their prey. Certain of these volatiles have notably been shown to decrease oviposition rates and increase egg predation on the emitting plant (Feeny et al., 1989;Kessler and Baldwin, 2001;Colazza et al., 2004;Johne et al., 2006). Moreover, some of these plant odors can also serve as direct repellents against insect herbivores (Quiroz et al., 1997;Bernasconi et al., 1998;Landolt et al., 1999). In addition, the volatiles can be perceived by neighboring plants and increase or potentiate their defensive abilities (Arimura et al., 2000a(Arimura et al., , 2000bBaldwin et al., 2006;Kessler et al., 2006;Kost and Heil, 2006). On the other hand, volatile signaling is also used by some plant species to overcome the vascular constraints on systemic signaling to elicit defense responses in adjacent leaves with little or no vascular connection with the attacked leaves (Frost et al., 2007).
Despite their common feature of being volatile at ambient temperatures, these major players of plant signaling are chemically highly diverse (Holopainen, 2004). HIPV bursts have been partially characterized in several plant species (for overview, see van Poecke and Dicke, 2004;Van Dam and Poppy, 2008), including, for instance, maize (Zea mays), cotton (Gossypium hirsutum), lima bean (Phaseolus lunatus), and native tobacco (Nicotiana attenuata). The major constituents of these VOC bouquets consist of products of the shikimic acid pathway, terpenes and fatty acid (FA) derivatives. The best characterized metabolites from this latter subpopulation are the green leaf volatiles (GLVs), originating from the degradation of C 18 FAs (linolenic and linoleic acids) when hydroperoxidated by a lipoxygenase and cleaved into C 12 and C 6 units by a hydroperoxide lyase Matsui, 2006).
Fatty acid amino acid conjugates (FACs), produced in the insect gut by conjugation of host-derived FAs to amino acids (Spiteller et al., 2000) and contained in the oral secretions (OS) of many caterpillar species, were the first nonenzymatic elicitors of HIPV characterized (Alborn et al., 1997;Halitschke et al., 2001). More generally, recent ''omic'' approaches have demonstrated that FACs are responsible for eliciting a large portion of the hundreds of genes regulated during the plant-herbivore interaction (Halitschke et al., 2003;Roda et al., 2004) as well as the reconfiguration of the proteome (Giri et al., 2006). However, such untargeted analyses have not been conducted for HIPV, and until now the regulation of only a few components of plants' volatile blends (e.g. trans-a-bergamotene in N. attenuata [Halitschke et al., 2001] and hexenylacetate or transb-farnesene in maize [Alborn et al., 1997]) has been associated with FAC signaling.
Notably, FACs are not always active: in lima bean and cowpea (Vigna unguiculata), for instance, FAC treatment of wounds does not elicit VOCs (Spiteller et al., 2001). Rather, Mithö fer et al. (2005) demonstrated that a continuous wounding treatment that resembled insect feeding was sufficient for HIPV emission. In cowpea, inceptin, an elicitor produced from the postingestive digestion of host plant ATP synthase, triggers changes in phytohormone and VOC production similar to those detected during herbivory (Schmelz et al., 2006;Carroll et al., 2008). In N. attenuata, the large methanol burst, which is many times larger than the GLV burst elicited by Manduca sexta herbivory, originates from the demethylation of pectin cell walls and in turn is elicited by the alkalinity of OS . b-Glucosidase hydrolytic enzymes from Pieris brassicae OS also trigger HIPV emissions when applied to cabbage (Brassica capitata) plants (Mattiacci et al., 1995). 2-Hydroxyoctadecatrienoic acid (2-HOT) is a newly identified component of M. sexta OS (E. Gaquerel, A. Steppuhn, and I.T. Baldwin, unpublished data). This oxylipin is produced from linolenic acid through the action of the a-dioxygenase (a-DOX; Hamberg et al., 2002Hamberg et al., , 2003. In N. attenuata, the production of 2-HOT during herbivory and its accumulation in OS allow the plant to monitor the progression of the insect's attack and to sustain its production of jasmonic acid, the central hormone that coordinates antiherbivore defenses. The striking postingestive stability and high activity of N. attenuata's a-DOX1 proteins in the alkaline and proteolytic milieu of the M. sexta midgut mediates 2-HOT's particular signaling role (E. Gaquerel, A. Steppuhn, and I.T. Baldwin, unpublished data).
Given the rapid advances in characterizing the active elicitors of M. sexta OS, a full unbiased analysis of N. attenuata's HIPV blend, an analysis that has not been conducted for any plant, is overdue. To examine the extent to which FACs, 2-HOT, and alkalinity contribute to the OS-elicited HIPV bouquet, we used comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry (GCxGC-ToFMS) in combination with multivariate statistics. This nontargeted approach was applied to a suite of elicitation treatments in which the amount of mechanical wounding was held constant.

A Stepwise Preprocessing of Volatile Blends
OS collected from M. sexta caterpillars fed on wildtype N. attenuata plants was purified by ion-exchange chromatography (IEX) to remove, among others, FACs and 2-HOT (IOS; Fig. 1A). This latter extract, as well as solutions of synthetic FACs and 2-HOT at concentrations that mimicked those found in OS (Halitschke et al., 2001;E. Gaquerel, A. Steppuhn, and I.T. Baldwin, unpublished data), were prepared as summarized in Figure 1B. We analyzed VOCs emitted locally during three consecutive 12-h periods (Fig. 1, B and C) after treating wounded (W) leaves with these solutions and compared these blends with those elicited by the complete OS or an alkaline buffer (B). We started VOC collections 1 h after treatment, when the effect of the mechanical disturbance had attenuated, to emphasize the exact contributions of each of the solutions tested.
Unbiased comparative analyses require exceptionally high resolution given the high structural diversity of HIPV constituents (Holopainen, 2004) and the range of concentrations at which they occur, which in this study was provided by GCxGC-ToFMS ( Fig. 2A). Once the 108 samples were analyzed, the raw chromatograms were deconvoluted to find analyte peaks using the LECO ChromaToF software version 2.21 (Fig. 2B). Routinely, approximately 600 peaks were detected. The peak lists were then checked prior to further processing for known artifact peaks and contaminants (e.g. column bleeding and plasticizers) identified from the analysis of ambient air samples trapped in the glasshouse.
Inconsistencies in the size of the peak tables rendered it difficult to place the information into a suitable matrix format in which rows represented an individual peak, columns represented an individual sample, and the values are peak intensities obtained by integration of deconvoluted peaks. As discussed previously by Shellie et al. (2005), the comparison feature embedded in the LECO ChromaToF software helped in aligning the different peak lists. This required the comparison of each individual chromatogram with a unique reference matrix in order to obtain output tables with the same numbers of peaks. To that end, for each of the three volatile-trapping periods, we selected the largest peak table obtained from the deconvolution of OS-induced samples as a reference matrix; this matrix contained the largest range of HIPV released during each of the sampling intervals. Samples from each of the three trapping periods were processed separately. A mass spectrum similarity threshold of 500/1,000 was used, as it provided sufficient selectivity for peak matching. Retention time (RT) window parameters for peak comparison were set as one modulation period for column 1 (C1) and 0.2 s for column 2 (C2); an evaluation of RT precision with these parameters revealed the parame-ters to be within the range of typical peak widths of nontailing peaks (Supplemental Table S2). Only the peaks providing matches with metabolites from the reference list (i.e. constituents of the OS-elicited HIPV profiles) were retained. These procedures shortened matrices to approximately 400 peaks, which were subsequently normalized by the tetralin internal standard @ RT (846; 3.76).

W1OS Triggers Pronounced Changes in Volatile Emission
In order to evaluate the performance and value of data deconvolution and pretreatment for statistical analysis, we compared chromatographic peak areas in W1OS and control (CTRL) samples. The intensity and statistical significance of the changes elicited after W1OS were examined using Volcano plots, a procedure commonly used to assess the quality of microarray data sets (Cui and Churchill, 2003). To that end, ANOVAs (W1OS versus CTRL) were performed on normalized log 2 -transformed intensities to screen for differentially regulated analytes. Volcano plots were obtained by plotting the negative log 10 -transformed P value issued from these ANOVA tests against the log 2 ratio of mean intensities (W1OS/CTRL). We selected a P value of 0.05 for statistical significance threshold and a 1.5-fold change ratio as cutoff threshold to discriminate between significantly up-and down-regulated peaks. sexta's OS were collected from third to fourth instar larvae reared on wild-type N. attenuata plants. OS purified by IEX (IOS), which are free of FACs, 2-HOT, and other FA metabolites, were prepared as described previously (Halitschke et al., 2001) using Amberlite IRA-400 resin. OS and IOS extracts were flushed with argon and stored at 220°C until use. B, The first fully expanded leaf (L11) of rosette-stage N. attenuata plants (randomized; n 5 6 biological replicates per treatment) was wounded (W) with a fabric pattern wheel, and 20 mL of one of the different eliciting solutions to test was directly applied to the leaf surface. Volatile emissions strictly induced by the mechanical wounding and/or the alkalinity of M. sexta's OS were assessed by applying a 0.1 M Tris, pH 9, buffer solution (B) containing 0.1% (v/v) Triton X-100 to wounded leaves. This nonionic surfactant was added to the B solution to evaluate its potential eliciting effect, as it was used for the preparation of the FAC solution. Synthetic N-linolenoyl-L-Gln (C18:3-Gln) and N-linolenoyl-L-Glu (C18:3-Glu), the two major FACs, and 2-HOT were diluted in the B solution at concentrations similar to those found in M. sexta's OS. Control plants (CTRL) were left untreated. C, Leaves were enclosed, 1 h after elicitation, in food-quality plastic containers secured with miniature clips. Plant volatiles were collected, using self-packed Super Q traps connected to a manifold vacuum pump, as described by Halitschke et al. (2000), during three consecutive periods of 12 h, the second one occurring during the glasshouse dark time. Background contaminants present in ambient air were also tested using empty trapping containers.
A total of 61 peaks were significantly more intense in chromatograms obtained from leaves induced by W1OS than from untreated leaves during the first trapping period. Concomitantly, no plant volatiles were significantly decreased (Fig. 3A). The number of significantly regulated peaks decreased to 23 and finally to 17 in the two next data points, which under-scores that W1OS treatment elicited rapidly waxing and waning changes in VOC releases. Prior to being ''blasted'' against libraries, lists of significantly regulated analytes were checked for false positives potentially created during peak detection and alignment. We estimated the false discovery rate to be 7% of the total of differentially regulated peaks. These errors, which were mainly produced during the alignment of peak matrices, were manually corrected, and contaminants remaining from the preprocessing were removed. Finally, significantly emitted plant-derived compounds were identified from comparison with standards, based on their mass spectral and retention index homologies shared with entries from the NIST and homemade VOC libraries or reports from the literature, or they were putatively assigned to VOC classes (Supplemental Table S1).

OS Reconfigures the Plant Wound Response Partly through the Action of FACs
To explore the extent to which the introduction of OS to wounds specifically regulated volatile emission, we compared the total changes elicited after W1OS with those detected after W1B (Fig. 4A). We called the ensemble of volatiles found differentially regulated in such comparisons the OS response (Supplemental Table S1, column C). More precisely, up-regulated volatiles were interpreted as being amplified by OS application, while down-regulated volatiles were considered repressed by this treatment. Similarly, volatile substances regulated by the wound and/or the alkalinity of OS were identified by statistically comparing Figure 2. Preprocessing of GCxGC-ToFMS data. A, Two-dimensional chromatograms (C1 3 C2) obtained after analysis of the volatile extracts by GCxGC-ToFMS were processed in two steps. B, Raw data were first deconvoluted to find analyte peaks using the LECO ChromaToF software version 2.21, and peak lists were then checked prior to further processing for known artifact peaks and contaminants. C, Peak matrices from each of three trapping periods were then aligned separately using the comparison feature embedded in the LECO ChromaToF software. For each data point, the largest peak table obtained from the deconvolution of samples collected from leaves wounded and treated with M. sexta's W1OS was used as a reference matrix. Peaks providing matches with metabolites from the reference list (i.e. W1OS-induced plant volatiles) were retained for statistical analysis. Processing parameters and peak list alignment procedures are explained in more detail in the text. Table S1, column B).

W1B and CTRL blends (Supplemental
During the first period of the experiment, the emission of 12 volatiles (58% of the W1OS response; Fig. 4B; Supplemental Table S1, column A): one acid (pentanoic acid, 3-methyl-, ANOVA W1OS versus W1B, P 5 0.043), one aldehyde (octanal, P 5 0.012), one aromatic (unknown @ RT [954; 4.13], P 5 0.012), five hexenylesters (cis-3-hexenylpropionate, P 5 0.005; cis-3-hexenyl,trans-2-butenoic ester, P 5 0.023; cis-3hexenyl-3-methylbutyrate, P 5 0.012; cis-3-hexenyltigilate, P 5 0.038; and cis-3-hexenyl-methylpentanoate, P 5 0.033), and four terpenoids (b-pinene, P 5 0.026; unknowns @ RT [714; 2.71], P 5 0.045; @ RT [720; 2.34], P 5 0.034; and @ RT [1,176; 3.27], P 5 0.039) was amplified by OS factors. By contrast, we observed that the emission of certain wound-inducible hexenylesters, such as cis-3-hexenylacetate (ANOVA W1B ver- Figure 3. W1OS triggers pronounced changes in volatile emission. A, Volcano plot analysis of the orientation, intensity, and statistical significance of the changes in chromatographic peak areas between volatile samples collected 1 to 13 h after N. attenuata plants were wounded with a fabric pattern wheel and treated with M. sexta's W1OS and those collected from untreated plants (CTRL). The log 2 ratio of mean intensities (W1OS/CTRL; n 5 6 biological replicates per treatment) is shown on the x axis, and the negative log 10 -transformed P value, from ANOVA on normalized log-transformed peak intensities, is shown on the y axis. Each dot represents an analyte from the common set of 463 normalized peaks between samples. The vertical dashed lines represent 1.5-fold change ratios, and the horizontal dashed line represents the statistical significance level at P 5 0.05. Peaks with large fold change values lie outside the vertical threshold lines. Up-and down-regulated peaks are localized in the top left and top right parts of the plot, respectively. Gray dots represent significantly regulated contaminants or false positives created by the peak alignment procedure. Plant volatiles found to be up-regulated when W1OS elicited were compared with CTRL samples and are highlighted in red. No plant volatiles were found to be significantly down-regulated. B, Classes of plant volatiles significantly regulated after plants were elicited by W1OS. Nineteen volatiles were emitted at significantly higher rates by plants elicited by W1OS than by untreated plants at 1 to 13 h after elicitation (Supplemental Table S1). Seven volatiles were upregulated at 13 to 25 h and 25 to 37 h after elicitation; none was down-regulated. Volatiles detected were classified as acids (e.g. pentanoic acid, 3-methyl-), aldehydes (e.g. octanal), aromatics (e.g. hexenylesters [e.g. cis-3-hexenylacetate]), methylesters (e.g. pentanoic, 3-methyl-, methylester), and terpenoids (e.g. b-pinene).
To determine the extent to which FACs, 2-HOT, or other OS components not removed by IEX could account for the OS response, we examined peak matrices for FAC-responsive (Supplemental Table S1, column E), 2-HOT-responsive (column F), and IOSresponsive (column D) volatiles. Following the same logic, we defined these responses as the subsets of volatiles being differentially emitted in W1FAC-, W12-HOT-, or W1IOS-treated samples when compared with W1B-treated samples. During the first trapping period, the FAC response accounted for 75% (three of four) of the amplified terpenoid emissions detected after OS was supplied to wounds (Supple- Figure 4. Application of OS reconfigures the plant wound response. A, Volcano plot analysis of the orientation, intensity, and statistical significance of the changes in chromatographic peak areas between volatile samples collected at 1 to 13 h after plants were wounded with a fabric pattern wheel and treated with M. sexta's W1OS and those collected from wounded plants treated with a buffer at the same pH as M. sexta's OS (W1B). The log 2 ratio of mean intensities (W1OS/W1B; n 5 6 biological replicates per treatment) is shown on the x axis, and the negative log 10 -transformed P value, from ANOVAs on normalized log-transformed peak intensities, is shown on the y axis. Each dot represents an analyte from the common set of 463 normalized peaks between samples. The vertical dashed lines represent 1.5-fold change ratios, and the horizontal dashed line represents the statistical significance level at P 5 0.05. Peaks with large fold change values lie outside vertical threshold lines. Up-and down-regulated peaks are localized in the top left and top right parts of the plot, respectively. Gray dots represent significantly regulated contaminants or false positives created by the peak alignment procedure. Plant volatiles found to be up-regulated and downregulated in W1OS-and W1B-elicited samples are highlighted in red and green. B, Two types of responses were considered from samples collected at 1 to 13 h, 13 to 25 h, and 25 to 37 h after elicitation: the W1OS response was detected by comparing W1OS to control (CTRL; 1; dashed circle; see also Fig. 3), and the OS response was detected by comparing W1OS versus W1B (2; black circle). Venn diagrams show the number of overlapping and nonoverlapping volatiles induced (red arrows) or repressed (green arrows) between those two sets of comparisons (Supplemental Table S1). Numbers specific to the W1OS response (1) represent W-induced/repressed volatiles for which no significant effect of the application of OS was observed. Numbers in the common area represent volatiles for which the application of OS amplifies the W regulation. Numbers specific to the OS response (2) represent W-repressed volatiles being induced by OS or W-induced volatiles being repressed by OS application. mental Table S1, column E): unknown @ RT (714; 2.71), ANOVA, W1FACs versus W1B, P 5 0.037; unknown @ RT (720; 2.34), P 5 0.044; and unknown @ RT (1,176; 3.27), P 5 0.032. Surprisingly, none of the OS-induced hexenylesters was significantly up-regulated after FAC application. Similar distinctions in the regulation of these two classes of volatiles were observed for the two next sampling times. Remarkably, FAC treatments inhibited the production of 2-pentanol,4-methyl (P 5 0.013), 2-pentanol,3-methyl (P 5 0.021), cis-3-hexenol (P 5 0.016), and 1-hexanol (P 5 0.020) in the same range of statistical significance as had been observed after application of OS. FAC signaling also accounted for the decreases, compared with wounded leaves, in the release of benzothiazole (P 5 0.047) and terpineol (P 5 0.004) at 13 to 25 h after treatment. As expected, removing FACs from OS by IEX abolished most of the FAC-mediated amplification and inhibitory patterns detected above (Supplemental Table S1, column D). Only two volatile substances, cis-3-hexenylbutyrate (P 5 0.047) and trans-a-bergamotene (P 5 0.039 and 0.007), were significant amplified by the 2-HOT treatment compared with the W1B response (Supplemental Table S1, column F).

The OS and FAC Responses Are Most Similar
The above analysis was conducted by comparing one volatile at a time. Although this univariate procedure reveals interesting patterns, it does not provide insights into the relationships among classes of treatments. Numerous multivariate procedures exist, and of these, hierarchical clustering analyses (HCA) and principal component analyses (PCA) are the most commonly used procedures in transcriptomic, proteomic, and metabolomic analyses of plant samples. Prior to performing such analyses, a filtering procedure is usually required to emphasize information linked to group separation. Thus, as recommended by Boccard et al. (2007), multiple univariate ANOVAs (P , 0.05, n 5 6 groups of treatments) were performed in an effort to eliminate, in a supervised way, analytes shared by all samples and therefore containing a common pattern of noninformative signals. Using this filter, only approximately 15% of the initial information was retained for the multivariate analyses. Again, these filtered lists of variables were checked for contaminants and false positives produced during peak matrix alignment.
To determine if these prefiltered variables were sufficient to discriminate classes of treatments, we performed HCA. Filtered data were autoscaled, and HCA was performed using the Euclidean distance as a clustering metric and the complete linkage aggregation method. For each of the three data points, preselected variables with equivalent abilities to classify samples were used to group treatments during HCA. Control samples collected 1 to 13 h after treatment were, as expected, clearly separated by the first branch of the dendogram (W response, Fig. 5). The second branch subdivided profiles collected from W1OSelicited samples from the other samples (OS response). The first branch of the dendogram obtained from the clustering of samples collected 13 to 25 h after the start of the experiment separated OS-and FAC-treated samples from those untreated or elicited with FACfree solutions (OS and FAC responses). For the last sampling time, 2-HOT-treated samples clustered with OS and FAC-elicited samples (OS and FAC and 2-HOT responses).

PCA Identified Multiple Regulatory Patterns of Volatile Production
PCA was used as an unsupervised method to produce interpretable projections of samples in a reduced dimensionality (score plots) and to highlight biomarkers responsible for group separation (loading plots). Prior to performing these analyses, different types of transformations were applied to correct for a nonconstant signal variance. Since for analytical chemical measurements the total uncertainty is often proportional to the expected value of the signal, a logarithmic transformation is often appropriate. Therefore, we log 2 transformed and mean-center normalized the values. The two first principal components (PCs) extracted accounted for approximately 60% of the total variance existing in the three sample populations (Supplemental Fig. S1). In addition, as expected from the HCA analyses, CTRL and W1OS samples were the treatments best separated by these two PCs. Examining other PCs that accounted for lower amounts of variance, we were unable to obtain the additional power needed to discriminate among treatments. To address this problem, we performed two new series of PCA per data point after dividing the data sets into two subparts.
The first category of PCA aimed at classifying CTRL, W1B, and W1OS samples (Fig. 6). The two first PCs extracted using those samples as inputs captured approximately 45% to 60% of the total variance (Fig.  6A) and segregated samples according to two biological responses. The first one represented the W response and the second was associated with the OS response, as described previously. In accordance with HCA results, group separation was less efficient for volatile profiles collected in the dark phase of the experiment (13-25 h). To simplify the interpretation, metabolites contributing to these projections were categorized as follows (Figs. 6B and 8 [below]; Table  I). OS-responsive volatiles, high-ranking loadings on PC1, were divided into those up-regulated (type I, W induced and OS amplified; 28 volatiles) and those down-regulated (type II, W induced and OS repressed; 11 volatiles). We classified as type III, high-ranking loadings on PC2, metabolites being W inducible but not OS responsive (seven volatiles). Globally, the clustering of the volatiles identified with this procedure was largely commensurate with that detected from the univariate statistical comparisons.
The contributions of B, IOS, 2-HOT, and FACs to the OS response were analyzed in a second series of PCA (Fig. 7). The two first PCs extracted using W1B, W1IOS, W12-HOT, W1FAC, and W1OS samples as inputs explained approximately 35% to 40% of the total variance (Fig. 6A). Even though none of the OS factors tested could completely restore the OS response, FAC-treated samples were the ones that most closely mimicked OS-elicited samples. As already revealed by HCA, this observation was more pronounced for samples collected at later times during the experiment. In addition, samples treated with IOS, which is free of FACs, were clearly separated from OStreated samples. An efficient segregation of W1B-and W12-HOT-treated samples was apparent only 25 to 37 h after elicitation. Following the same logic as presented above, we discovered that three main biological patterns structured those projections (Figs. 7B and 8; Table I). Type I.a volatiles, high-ranked loadings on PC1 and PC2, were amplified by OS through the action of FACs (13 volatiles; W induced and OS and FAC amplified). Terpenoid metabolites were the predominant volatiles regulated in this manner (Table I). Surprisingly, we also detected that the application of 2-HOT amplified the production of trans-a-bergamotene (W12-HOT versus W1B; ANOVA, 13225 h, P 5 0.039; 25-37 h, P 5 0.007; Supplemental Table S1, column F; Fig. 7B; Supplemental Fig. S2), a well-known FAC-regulated sesquiterpene (Halitschke et al., 2001). We called the volatiles with this pattern of regulation type I.b (W induced and OS, FAC, and 2-HOT amplified). Metabolites that were amplified by unknown OS factors were classified as type I.c (12 volatiles). Most of these were identified as hexenylesters and aromatic derivatives ( Table I). Volatiles that were repressed by FACs corresponded mainly to those classified in the aforementioned type II (11 volatiles; W induced and OS repressed). Altogether, the number of differentially regulated volatiles largely decreased over the experimental time, while the relative proportion of FACregulated substances increased (Fig. 8B).

DISCUSSION
We present here a conceptual approach for the untargeted comparative processing via GCxGC-ToFMS of volatile bouquets emitted after insect herbivory. When combined with univariate and multivariate statistics, this approach (1) identifies wound-dependent metabolites and those being further modulated by the application of M. sexta's OS and (2) demonstrates that FACs act as major orchestrators of the OS-elicited response, eliciting and repressing suites of volatiles (Fig. 8).

Nontargeted Functional Analysis of GCxGC-ToFMS Data
The number of biological questions in plant science to which metabolomics have been applied is exponentially growing (for review, see Shulaev et al., 2008). However, a majority of these studies have focused on the regulation of endogenous metabolites rather than on what is sometimes called the plant's exometabolome (e.g. root exudates and floral nectars). As dis- Figure 5. Hierarchical clustering analyses of sample specimens. N. attenuata volatiles were collected at 1 to 13 h, 13 to 25 h, and 25 to 37 h after elicitation. For each data point, an ANOVA-based filter (P , 0.05) was applied on normalized peak intensities in an effort to eliminate redundant signals from GCxGC-ToFMS matrices (see text). This shortened the peak lists by approximately 70% to 80%. The resulting data sets were checked manually for contaminants and false positives from the peak alignment procedure, autoscaled, and classified by hierarchical clustering analysis with The Institute for Genomic Research MultiExperiment Viewer software using the Euclidean distance as a grouping parameter and the complete linkage aggregation method. Symbols indicate different groups of treatments: OS (black squares), IOS (white squares), FACs (black triangles), 2-HOT (white triangles), and buffer (B; white diamonds) solutions were applied onto leaves wounded with a fabric pattern wheel. Control plants (CTRL; white circles) were left untreated. One outlier sample (midvein disrupted during handling) was deleted from the W1FACs, CTRL, and W1B classes of treatments. cussed in the recent publication of Van Dam and Poppy (2008), this also applies to metabolites emitted in the gas phase by plants, for which few studies have been carried out using unsupervised approaches.
GCxGC-ToFMS technology has been shown to provide high-quality mass spectra with great sensitivity, largely as a result of the enhanced resolution and zone compression obtained from the orthogonal separation (Shellie et al., 2001). These qualities were exemplified by the 5-to 10-fold increase in the signal-to-noise ratio of a series of commonly emitted volatiles that we compared with the results from one-dimensional analyses performed with the same instrument (data not shown). The approach was found to be reliable and accurate for the automated processing of a large number of samples using parameters based on the recommendations described by Shellie et al. (2005) and preliminary analyses performed with standard compounds (Fig. 2; Supplemental Table S2). For each collection period, each individual chromatogram was compared with the OSelicited volatile profile reference list. This procedure reduced the size of the output tables presenting the same numbers of peaks and enriched in OS-specific volatiles. These peak matrices, as already documented by several authors (Shellie et al., 2005;Hagan et al., 2007;Kusano et al., 2007;Tu et al., 2007), were suitable for univariate and multivariate statistical analyses (Figs. 3-7), after correction for false positives resulting from the alignment procedure. With the caveat that these aligning and statistical procedures might exclude metabolites up-regulated only by FACs or 2-HOT and not by the OS treatment, they allowed regulatory patterns orchestrating OS induction to be efficiently discriminated (Fig. 8). In that sense, this study constitutes one of the first applications of GCxGC-ToFMS for the nontargeted processing of plant volatile bouquets. Using the approach developed here, we also report the use of this procedure in the analyses of polymorphism in volatile emission among individual plants collected from a native N. attenuata population (M.C. Schumann, N. Heinzel, E. Gaquerel, A. Svatos, and I.T. Baldwin, unpublished data).  (Table I). To simplify the analysis, we categorized volatile accumulation patterns into three main types: type I contains W-induced and OS-amplified volatiles (e.g. cis-3-hexenylmethylpentanoate between 1 and 13 h); type II contains those being W induced and OS repressed (e.g. cis-3-hexenol between 1 and 13 h); and type III contains those detected as W induced but not OS responsive (e.g. hexylisobutyrate between 1 and 13 h). Normalized intensities and magnified areas of representative two-dimensional chromatograms obtained from W1OS, W1B, and CTRL samples illustrate for one member each of the three types of regulation.

FACs Are the Major Determinants of HIPV Emissions
Applying FACs to wounded leaves has been shown to trigger the emission of particular volatiles in amounts similar to those detected during caterpillar feeding (Alborn et al., 1997;Halitschke et al., 2001). Here, we rediscovered through an unbiased approach the dominant role of this class of elicitors in coordinating N. attenuata's volatile emission. Precisely, we detected that the FAC response accounted for approximately half (sample time 1, seven of 17; sample time 2, seven of eight; sample time 3, three of five) of the significant changes in volatile emission elicited by OS Table I. Summary of differentially emitted volatiles and their corresponding regulatory patterns N. attenuata leaf volatiles were collected at 1 to 13 h, 13 to 25 h, and 25 to 37 h after elicitation, and those differentially emitted (ANOVA, P , 0.05, n 5 6) between the groups of treatments were analyzed by PCA (Figs. 6 and 7). Types I.a to I.c and type II contain leaf volatiles that were induced after W and induced or repressed after application of M. sexta OS. Amplification of the wound (W) response was mediated by the action of FACs (type I.a), the action of FACs or 2-HOT (type I.b), or the action of unknown factors (type I.c). Type II represents volatiles induced after W that were repressed after the application of M. sexta OS mainly through the action of FACs. Type III represents W-induced volatiles for which no significant effect of the application of OS was observed. Volatiles are listed by chemical classes with their RTs on the first and second dimensions. Numbers indicate the type of identification (Id.): 1, comparison with standards; 2, putative identification. Retention indices were calculated according to Kováts (1958) and compared with literature (Lit.) reports: a (Ruther, 2000), b (Adams, 2004), c (Engel et al., 2002), d (Triqui and Reineccius, 1995), e (Guth and Grosch, 1991), and f (Jirovetz et al., 2002). treatments. Terpenoids defined the largest proportion of FAC-inducible volatiles (Table I; Supplemental Table S1). In addition to trans-a-bergamotene, whose regulation through FAC signaling has already been described (Halitschke et al., 2001), a suite of new FAC targets was discovered, for instance unknown @ RT (720; 2.34) and unknown @ RT (1,176; 3.27). The identification of these two metabolites will require further work, as will the determination of their potential activity in recruiting predators. In contrast, the diversity of hexenylesters emitted after W1OS treatment appeared to be either strictly associated with the mechanical stress or regulated through the action of unknown OS factors (Table I; Supplemental Table S1).
Even though FACs could not fully account for all of the qualitative changes in volatile emission observed after OS treatment, FAC-elicited profiles were the most closely related to OS-elicited ones, as was clearly evident from the HCA dendogram and PCA projections. These results extend our appreciation of the importance of FAC signaling during attack by M. sexta larvae. Wu et al. (2007) demonstrated that the same FACs that we used were sufficient to induce similar levels of salicylic acid-induced protein kinase (SIPK) and wound-induced protein kinase (WIPK) activity, which, in turn, are responsible for eliciting most of the defensive responses elicited in N. attenuata leaves when attacked by caterpillars or OS elicited. Additionally, FACs are also known to be responsible for the majority of the dramatic changes observed in the transcriptome (65%-86% of the regulated transcripts [Halitschke et al., 2001;Roda et al., 2004]) and the Figure 7. PCA of OS-, FAC-, and 2-HOT-responsive volatiles. A, Filtered data sets used for HCA (Fig. 5) were mean centered and subjected to PCA using an Excel add-in developed by the Bristol Chemometrics group. Projection plots were obtained from the coordinates calculated for the two first PCs extracted. The percentage of the total variance explained by each PC is reported in parentheses. Symbols indicate different groups of treatments: OS (black squares), IOS (white squares), FACs (black triangles), 2-HOT (white triangles), and buffer (B; white diamonds) solutions were applied to leaves wounded with a fabric pattern wheel (W). Centers of each group of treatment were calculated and are indicated on projections by crosses. B, For each data point, analytes contributing to the PCA projection were examined and screened for plant volatiles (Table I). To simplify the analysis, we categorized volatile accumulation patterns into two main types and three subtypes: type I volatiles, whose emission was W induced and OS amplified, with type I.a for FAC induced (e.g. unknown sesquiterpene @ RT [720; 2.34] between 1 and 13 h), type I.b for OS, FAC, and 2-HOT induced (e.g. trans-a-bergamotene between 25 and 37 h), and type I.c containing volatiles that were OS induced through the action of unknown factors (e.g. cis-3-hexenylmethylpentanoate); type II (see Fig. 6) represents OS-and FAC-repressed volatiles (e.g. cis-3-hexenol between 1 and 13 h). Normalized intensities and magnified areas of representative twodimensional chromatograms obtained from W1OS, W1B, and control (CTRL) samples illustrate for one member each of the four types of regulation.
Many plants change the blend of volatiles released when attacked by different herbivore species (for review, see Dicke, 1999), and FACs may play a central role in tailoring the VOC blend. Whether plants integrate the whole FAC signature or only respond to dominant components remains unknown. Wu et al. (2007) demonstrated that the four most abundant FACs found in M. sexta OS are functionally redundant in terms of kinase activations. Recent analytical developments have demonstrated that the spectrum of FACs contained in the OS of many caterpillar species is more complex than previously thought (Yoshinaga et al., 2007). Clearly, additional unbiased analyses of the specific and common responses elicited by each of these FACs are needed.

FACs Suppress Particular Volatiles
One of the tantalizing discoveries from this unbiased analysis was that M. sexta OS strongly repressed the emission of particular components of the volatile bouquet (Fig. 4). Repressed metabolites included shortchain alcohol derivatives, certain hexenylesters, as well as some aromatic and terpenoid derivatives ( Fig. 4; Supplemental Table S1). The analysis demonstrated that this OS-mediated repression (type II) could be largely attributed to the action of FACs (Supplemental Table S1). Treating wounded leaves with FACs reproduced, in the same range of intensity and statistical significance, approximately 80% of the alterations in volatile emission observed during the first trapping period after OS treatment. This effect was particularly well exemplified by the repression of cis-3-hexenol production detected in the first hour after OS or FAC treatment. This constitutes, to our knowledge, the first observation of a dual regulatory role for FACs on a plant's volatile release.
Factors that suppress the deployment of woundinduced defenses have been characterized in the OS of different caterpillar species. Notably, Musser et al. (2002Musser et al. ( , 2005 demonstrated that Glc oxidase (GOX), an enzyme isolated from the salivary glands of Helicoverpa zea and subsequently found in other caterpillar species, suppresses the wound-inducible accumulation of foliar nicotine in cultivated tobacco (Nicotiana tabacum) and is associated with increased survival and growth of larvae that fed on GOX-treated plants compared with those that fed on leaves treated with water. GOX enzymatic activity also decreased the emission of nicotine into the head space as well as other commonly emitted OS-induced volatiles (Delphia et al., 2007). GOX represses the transcription of genes involved in terpenoid biosynthesis in alfalfa (Medicago truncatula; Bede et al., 2006). C. Diezel and I.T. Baldwin (unpublished data) found that, when the antagonistic action of ethylene signaling is repressed, the GOX activity contained in Spodoptera exigua OS triggers the accumulation of salicylic acid in N. attenuata, which, in turn, antagonizes jasmonic acid-inducible defenses, such as the aforementioned production of nicotine. The role of ethylene signaling in the repressing effect of FACs on the emission of certain VOCs, for instance, will be particularly interesting to investigate further. Alternatively, it should be recalled that the repressing effect of FACs might not necessarily be the result of signaling cross talk. For example, FACs might directly interact with wound-released volatiles at the treatment site. This hypothesis could be tested by comparing FAC-elicited blends emitted locally and systemically.
The ecophysiological consequences of this novel repression function of FACs will be exciting to explore sexta's OS during W1OS elicitations. Types I.a to I.c and type II contain leaf volatiles that were induced after W and significantly induced or repressed after application of M. sexta's OS. Amplification of the wound response was mediated by the action of FACs (type I.a; dark square), the action of FACs or 2-HOT (type I.b; dark gray square), or the action of unknown factors (type I.c; gray square). Type II (light gray square) represents volatiles induced after W that were repressed after application of M. sexta's OS through the action of FACs. Type III (white square) represents W-induced volatiles for which no significant effect of the application of OS was observed. B, Number of volatiles significantly regulated following types I to III (see Table I) at 1 to 13 h, 13 to 25 h, and 25 to 37 h after elicitation.
in future work. The central question will be to determine whether this repression is advantageous for the plant or the attacking insect. GLV emissions have been shown to increase the apparency of N. attenuata plants to herbivores Meldau et al., 2009), and getting ''off the herbivore's radar'' can provide a level of protection for plants that is comparable to that afforded by direct defenses. Meldau et al. (2009) observed that although most of the chemical defenses were reduced in plants rendered silenced in SIPK and WIPK signaling, these plants did not suffer more damage from herbivores when planted into a native habitat. This lack of herbivore attack was associated with a large attenuation of GLV emissions from SIPK-and WIPK-silenced plants, and the authors demonstrated that applying GLVs, particularly cis-3hexenol, restored, under laboratory conditions, M. sexta tissue consumption and mass gain. Alternatively, some of the same VOCs that are repressed by FACs may be involved in attracting Geocoris pallens predators and thereby interfere with the defensive function of the volatile release (Kessler and Baldwin, 2001). To further complicate the matter, neighboring unattacked plants might react to these modifications of the VOC blend composition produced by attacked plants. In that direction, Paschold et al. (2006) demonstrated that VOC blends deficient in GLVs or trans-a-bergamotene regulated numerous defense-related genes in neighboring plants; this regulation did not occur when these plants were exposed to the complete VOC blends. Clearly, much more work will be required to understand the function of FAC-mediated VOC signaling and to determine the net fitness consequences for the plant of these VOCs as cues for herbivores, predators, or neighboring plants.

Regulation of trans-a-Bergamotene Emission
Consistent with the results of Halitschke et al. (2001), we found that FACs trigger trans-a-bergamotene emission and that their removal from M. sexta OS by IEX decreased the release of this sesquiterpene. We also found that 2-HOT, when applied to wounded leaves in the same range of concentration as that calculated in M. sexta OS, was able to elicit the release of the same amounts of trans-a-bergamotene as are elicited by OS and FAC treatments. We recently discovered that both 2-HOT and FACs are removed from OS by IEX purification (E. Gaquerel, A. Steppuhn, and I.T. Baldwin, unpublished data), which explains why IOS-treated plants have severely reduced emissions of trans-a-bergamotene. 2-HOT is produced from linolenic acid through the action of the a-DOX pathway and accumulates in planta during pathogen inoculation (Hamberg et al., 2002(Hamberg et al., , 2003. In N. attenuata, the production of 2-HOT in the larval OS during the herbivory allows the plant to monitor the progress of the insect attack (E. Gaquerel, A. Steppuhn, and I.T. Baldwin, unpublished data). This appeared to be controlled by the striking stability of N. attenuata's a-DOX1 proteins in the alkaline and protease-filled gut of M. sexta larvae. Moreover, M. sexta larvae feeding on N. attenuata's a-dox1-silenced lines grew better than on wild-type ones and accumulated less 2-HOT in their OS when feeding on these plants. The elicitation of trans-a-bergamotene release as an indirect defense represents a new testable hypothesis for the function of 2-HOT, a function that could be readily addressed with a-dox1-silenced plants.

Plant Material and Treatments
We used an isogenic line, obtained after 22 generations of inbreeding, of Nicotiana attenuata (synonymous with Nicotiana torreyana; Solanaceae) derived from field-collected seeds. Seeds were germinated on agar plates containing Gamborg B5 medium (Duchefa; http://www.duchefa.com) as described by Krü gel et al. (2002). Plants were grown in the glasshouse in 1-L individual pots at 26°C to 28°C under 16 h of light supplied by Philips Sun-T Agro 400-or 600-W sodium lamps.
All of the eliciting treatments were performed with plants at the rosette stage. Six replicate plants for each treatment were used. Plants were divided into six groups of equal size, and each of six treatments was randomized inside each of these groups. The second fully elongated (11 position) leaf was used. Plants were either left untreated (CTRL) or wounded using a pattern wheel to punch three rows of holes on each side of the midvein. Wounded leaves were immediately treated with 20 mL of the eliciting solution pipetted directly onto the wounded leaf and gently dispersed across the surface with a gloved finger, changing gloves between treatments.
The compositions of the different eliciting solutions are summarized in Figure 1. Volatile emissions strictly induced by the mechanical wounding and/or the alkalinity of Manduca sexta's OS were assessed by applying a 0.1 M Tris, pH 9, buffer solution (B) containing 0.1% (v/v) Triton X-100. This nonionic surfactant was added to evaluate its potential eliciting effect, since it was used for the preparation of the FACs solution. OS were collected from third to fourth instar M. sexta larvae reared on wild-type plants, flushed with argon, stored at 220°C, and diluted 1:1 (v/v) before being used with a 23 concentrated buffer solution. OS purified by IEX (IOS), which are FACs and 2-HOT free, were prepared as described previously (Halitschke et al., 2001) using Amberlite IRA-400 resin (Sigma; http://www.sigmaaldrich.com/), stored at 220°C, and diluted 1:1 (v/v) before use with a 23 concentrated B solution. Synthetic N-linolenoyl-L-Gln (C18:3-Gln) and N-linolenoyl-L-Glu (C18:3-Glu), the two dominant FACs, and 2-HOT (Larodan Fine Chemicals; http://www.larodan.se/) were diluted in the buffer solution at 0.2, 0.4, and 0.1 nmol mL 21 , respectively, concentrations similar to those found in M. sexta's OS (Halitschke et al., 2001;E. Gaquerel, A. Steppuhn, and I.T. Baldwin, unpublished data).

Collection of Plant Volatiles
Treated leaves were enclosed, 1 h after elicitation (the time needed to complete all of the treatments), in two 50-mL food-quality plastic containers (Huhtamaki; http://www.huhtamaki.com/) secured with miniature clawstyle hair clips. Ambient air flowed into the cage primarily through a clippedoff P1000 pipette tip inserted into the bottom container and was pulled out through a self-packed glass tube (ARS, Inc.; http://www.ars-fla.com/_fpclass/ fp_contact.html) containing glass wool and 20 mg of Super Q (Alltech; http:// www.discoverysciences.com) and secured in a second clipped-off P1000 pipette tip inserted into the top container. Air flow was created by a manifold vacuum pump (model DAA-V114-GB; Gast Manufacturing; http:// www.gastmfg.com/) as described by Halitschke et al. (2000). Plant volatiles were collected during three periods of 12 h, the second one occurring during the dark period. Background contaminants present in ambient air were also collected using empty trapping containers. Immediately after collection, traps were eluted by spiking each with 400 ng of tetralin as an internal standard and flushing the trap with 250 mL of dichloromethane into a GC vial containing a glass insert.

Instrumental Parameters
An Agilent 6890N gas chromatograph equipped with an Agilent 7683 autoinjector (Agilent Technologies; http://www.agilent.com/) coupled with a LECO Pegasus III time-of-flight mass spectrometer with a 4D thermal modulator upgrade (LECO; http://www.leco.com/) was used to collect the three-dimensional GCxGC-TOFMS data. The GC inlet and transfer line were held constant at 250°C. Splitless injections of 1 mL were made onto an RTX-5MS column (i.e. column 1 [C1] of the GCxGC system [20 m 3 250 mm i.d. 3 0.5 mm; Restek, Bellefonte; http://www.restek.com/]). The collected C1 effluent was transferred to a DB-17 column 2 (C2) of the GC3GC system (0.890 m 3 100 mm i.d. 3 0.1 mm; Agilent Technologies) every 6 s (modulation time). C1 was held at 40°C for 5 min and then increased at 5°C min 21 to 190°C and finally increased at 25°C min 21 to 250°C, where it was held for 5 min. C2 was initially set at 45°C and followed the same temperature program as C1, giving a total run of 45 min. The modulator was maintained at 30°C higher temperature than C1. The ion source was maintained at 250°C. Data were collected, after a solvent peak delay of 120 s, in the mass-to-charge ratio range 50 to 300, at a rate of 200 spectra s 21 .

Data Preprocessing and Statistical Analysis
The LECO ChromaToF software version 2.21 (LECO) was used throughout to control the instruments as well as to acquire and process the data (including automatic peak deconvolution). Sample populations from each of three trapping periods were processed separately. We selected for each of the three trapping periods the largest peak table obtained from the deconvolution of OS-induced samples as a reference matrix for each of the volatile-trapping periods (Fig. 2). Processed peaks were reported at a signal-to-noise ratio of 10, since we estimated this value as the minimum required for integration and accurate identification during blasts with the NIST and homemade libraries (data not shown). Normalization (object-wise standardization) of analytical profiles is an important step in the preprocessing of metabolic profiling data. Here, we applied linear normalization by dividing all peak heights by the tetralin internal standard @ RT (846; 3.76). To minimize the effect of zero substitution, 1 was added to all areas prior to normalization.
The variations in intensity across treatments were analyzed by univariate ANOVA (P , 0.05) on log 2 -transformed normalized peak areas. Hierarchical clustering analyses were done with The Institute for Genomic Research MultiExperiment Viewer 3.1 software (http://www.tm4.org/mev.html) on autoscale using the Euclidean distance as a clustering metric and the complete linkage aggregation method. The statistical robustness of the dendogram structures was tested by the bootstrap resampling procedure. PCA was performed on log 2 -transformed and mean-centered values using an Excel add-in developed by the Bristol Chemometrics group (http://www.chm.bris. ac.uk/org/chemometrics/). Projection plots were obtained from the coordinates calculated for the PCs extracted.
VOCs markers underlying treatment-based discriminations were searched against a custom spectrum library constructed from authentic standards and the NIST98 standard and identified based on retention index and spectrum similarity match. For determination of the retention index, a C8 to C24 n-alkanes series was used. Only the first-dimension RT was used to calculate the retention index according to Kováts (1958), given that the second RT value was negligible.

Supplemental Data
The following materials are available in the online version of this article.
Supplemental Table S1. Induced changes in volatile emission.
Supplemental Table S2. Evaluation of injection, RT, and calibration precisions.