Unbiased profiling of volatile organic compounds in the headspace of Allium plants using an in-tube extraction device.

BACKGROUND: Plants produce and emit important volatile organic compounds (VOCs), which have an essential role in biotic and abiotic stress responses and in plant-plant and plant-insect interactions. In order to study the bouquets from plants qualitatively and quantitatively, a comprehensive, analytical method yielding reproducible results is required.
RESULTS: We applied in-tube extraction (ITEX) and solid-phase microextraction (SPME) for studying the emissions of Allium plants. The collected HS samples were analyzed by gas chromatography-time-of-flight-mass spectrometry (GC-TOF-MS), and the results were subjected to multivariate analysis. In case of ITEX-method Allium cultivars released more than 300 VOCs, out of which we provisionally identified 50 volatiles. We also used the VOC profiles of Allium samples to discriminate among groups of A. fistulosum, A. chinense (rakkyo), and A. tuberosum (Oriental garlic). As we found 12 metabolite peaks including dipropyl disulphide with significant changes in A. chinense and A. tuberosum when compared to the control cultivar, these metabolite peaks can be used for chemotaxonomic classification of A. chinense, tuberosum, and A. fistulosum.
CONCLUSIONS: Compared to SPME-method our ITEX-based VOC profiling technique contributes to automatic and reproducible analyses. Hence, it can be applied to high-throughput analyses such as metabolite profiling.

Background

Plants produce various kinds of volatile organic compounds (VOCs) that are a part of the metabolome. By today the total number of identified VOCs is about 1700, and they account for 1 % of secondary metabolites [1, 2]. The major chemical classes of VOCs emitted from plants are terpenoids, phenylpropanoids/benzenoids, and derivatives of fatty acids and amino acids [3]. The genus Allium, is comprised of onions, leeks, and garlic, the total number of species is to 600–750 [4]. Allium plants can produce sulfur-containing VOCs through enzymatic reaction of sulfur-storage compounds [4]. For example, primary “aroma” compounds are thiosulfinates including allicin that are produced from aliphatic cysteine sulfoxides as “aroma” precursors in Genus Allium. Dithiins, ajoenes, and sulfides are known to as secondary “aroma” compounds [5].

We chose Allium fistulosum (Japanese bunching onions), A. chinense (rakkyo), and A. tuberosum (Oriental garlic), because these plants have been cultivated in Japan since the 8th century and are favorites of the Japanese. Allium plants emit VOCs that result in strong odors. The odoriferous compounds whose moiety contains sulfur in their moieties function not only as a defense against pathogens [6] and insects [7], but they also attract special herbivores and insect-eating insects such as moths [8, 9] and bees [10]. The chemical composition of such metabolites is diverse [11]. Since sulfur-containing VOCs produced by Allium plants exhibit anticancer- [12, 13], antithrombotic- [14], and antibacterial activity [15, 16], they are thought to be beneficial to human health.

There are several methods to collect VOCs in various matrices. A traditional way is steam distillation by which oils produced by plants can be collected. Meanwhile headspace (HS) sampling is a non-destructive solvent-free method for collecting VOCs emitted from plants [17, 18] including vegetables [19], humans [20], and microbes [21]. Moreover, the investigation of HS composition is much more meaningful than volatile analysis of samples collected by distillation or extraction methods. In case of high concentration capacity HS (HCC–HS) sampling methods [17, 22] such as solid-phase microextraction (SPME) [23-25], in-tube extraction (ITEX) [26-28], and stir bar sportive extraction (SBSE) [29, 30], VOCs can be easily concentrated. However, there is still a need for developing a comprehensive, reproducible, and high-throughput analysis for detection and quantification of VOCs in biological samples of various cultivars. Less than approximately 20 samples can be analyzed as one batch with one SPME fiber due to capacity of sorbent materials of SPME. Trapping of VOCs depends on SPME fibers’ properties [31]. To date, many types of sorbent materials are commercially available for ITEX-based method. Choosing the appropriate sorbent material of ITEX is important to trap non-polar and/or polar VOCs. Compared to SPME-method sampling according to ITEX procedure is fully-automated at the four steps, i.e., sample conditioning, analyte extraction/sorption, desorption/injection, and trap conditioning. Plus more samples can be analyzed by using ITEX- than SPME-method [32]. As high-throughput analysis is required for VOC profiling, we applied ITEX method in this study.

After HCC-HS sampling, VOCs are directly analyzed by gas chromatography (GC)-based techniques, because target analytes are easily released by heating sorbent materials. Of these methods, GC combined with electron ionization–time-of-flight–mass spectrometry (EI-TOF–MS) may help to identify and estimate the structure of VOCs, because EI-TOF–MS yields comprehensive information on molecular fragments in terms of mass-to-charge ratios [33], and because of well-documented libraries such as NIST/EPA/NIH mass spectral library (NIST-L) [34], Adams library (Ad-L) [35], the terpenoids library (Te-L; http://www.massfinder.com/wiki/Terpenoids_Library), and VocBinbase (Vo) [29], which contain mass spectral and retention index (RI) information of compounds that can be analyzed by GC–MS. Furthermore, several alignment tools such as AMDIS [36], ChromA [37], H-MCR [38], metalign [39], Tagfinder [40], and XCMS [41] have been developed and are freely available for GC–MS data interpretation.

The goal of the study was to develop a comprehensive, reproducible, and high-throughput profiling method for VOC collection from many samples by using fully-automated ITEX procedure and to then provisionally identify the detected VOCs in the HS of plants using the summarized mass spectral libraries. By applying our pipeline, we performed comprehensive HS-VOC profiling of the sheaths and basal plates of 12 Allium cultivars with ITEX-method in this study.

Results and discussion

Optimization of HCC-HS sampling and comparison of HS-VOC profiles in the HS of Allium fistulosum using ITEX and SPME techniques

To achieve the best performance for HCC-HS sampling in HS-GC vials, a method is needed that suits the goal of comprehensive, reproducible analyses. To this end, we modified the method of Tikunov et al. [23] and Kusano et al. [31]. Allium plants produce sulfides such as dipropyl disulfide as the main VOC component [7, 42]. We used ITEX and SPME to conduct HCC-HS collection from the Allium plants and evaluated statistically. The choice of the internal standards (ISs) is also critical for non-targeted metabolite profiling [43, 44]. Several ISs with different physicochemical properties (i.e., RI and chemical structure) are required for a comprehensive VOC analysis to evaluate whether the analytes participate in cross-contribution [44] and whether the RI of each IS peak is reproducible. Therefore we carefully examined dissolving agents for the ISs based on the value of the partition coefficient and the solubility of each IS [45] and selected methanol as the solvent.

We conducted HCC-HS sampling using ITEX and SPME to compare their performance for peak detection and to assess their comprehensiveness and the reproducibility of the results obtained with each technique. First we estimated the lower limit of quantification (LLOQ) of dipropyl disulfide, the major disulfide in A. fistulosum [7] and A. cepa [46], using ITEX- and SPME-GC-TOF–MS (see Additional file 2). The LLOQ of the peak detected by ITEX-GC-TOF–MS analysis was 250 pmol; it was 25 pmol by SPME-GC-TOF–MS analysis (data not shown). Then, using both methods, we analyzed the sheath and the basal part of A. fistulosum (brand name, Mikata spring onion; class01 in Table 1, Fig. 1, Additional file 1). The total ion chromatogram (TIC) of each analyte showed that peak detection was more sensitive with SPME device (Additional file 1). The score scatter plot of samples analyzed with the ITEX device and the SPME fiber showed clear separation of the first principal component (Additional file 1). It may be due to the use of different resins (TGR/CSIII for ITEX and PDMS/DVB for SPME).

Table 1

Allium species used in this study

Class in PCA	Binomial name	Species name	Bland name	Harvested field in Japan
08	Allium chinense	Rakkyo	Young rakkyo	Namegata, Ibaraki
01	Allium fistulosum	Spring onion	Mikata spring onion	Hamamatsu, Shizuoka
02	Allium fistulosum	Green spring onion	Aoi-chan green spring onion	Akitakata, Hiroshima
03	Allium fistulosum	Scallion	Hakata scallion	Hakata, Fukuoka
04	Allium fistulosum	Green spring onion	Green spring onion from Nagareyama	Nagareyama, Chiba
05	Allium fistulosum	White spring onion	White spring onion from Nagano	Nagano
06	Allium fistulosum	Leek	Shimonita leek	Gunma
07	Allium fistulosum	Leek	Shirakami leek	Noshiro, Akita
09	Allium fistulosum	Scallion	Kujo scallion	Nagahama, Shiga
11	Allium fistulosum	Spring onion	Goudo spring onion	Anpachi, Gifu
10	Allium tuberosum	Oriental garlic	Oriental garlic	Nagahama, Shiga
12	Allium fistulosum	Red spring onion	Red spring onion	Tsuruoka, Yamagata

All samples were harvested in October 2012

We used the same class names in the PCA score scatter plot (see Fig. 4)

Fig. 1

Visual phenotypes of the Allium samples used in this study. a A. chinense (rakkyo, class08 in Table 1 and Fig. 4). b A. fistulosum (spring onion, class01), c A. fistulosum (green spring onion, class 02), d A. fistulosum (scallion, class03), e A. fistulosum (green spring onion, class04), f A. fistulosum (white spring onion, class05), g A. fistulosum (leek, class06), h A. fistulosum (leek, class07), i A. fistulosum (scallion, class09), j A. fistulosum (spring onion, class11), k A. tuberosum (Oriental garlic, class10), l A. fistulosum (red spring onion, class12). The red and white areas of the scale bar are 5 cm

Fig. 4

Score (a) and loading (b) plots of PCA of the VOC profiles of the Allium samples. Principal components one and two (PC1, PC2) represent the first two principal components that account for a total of 53.2 % of the variance. Each plot represents an independent plant. In the loading plot, black dots and white triangles represent tentatively identified- and unknown peaks, respectively. All compound names and IDs are listed in the Additional file 3. ID011 2-butenal, 2-ethyl; ID026 3,4-dimethylthiophene; ID091 diallyl disulphide; ID154 2,5-thiophenedicarboxaldehyde

An advantage of the use of ITEX lies in its use of a stainless-steel needle and a special purge-, trap-, and trap-cleaning system [27, 47, 48]. This makes it possible to run more samples while maintaining the high reproducibility of data. On the other hand, the SPME method is appropriate for semi-targeted analysis because it features a wide variety of fibers rather than the ITEX sorbent materials. Although the SPME showed more sensitivity than ITEX to detect VOCs, less than approximately 20 samples can be analyzed as one batch due to capacity of sorbent materials of SPME. Thus, we applied the ITEX method for non-targeted HS-VOC profiling of the Allium samples.

Comparison of the libraries for the tentative identification of HS-VOCs

We estimated how many volatile peaks in the mass spectra overlapped in NIST-L, Ad-L, Te-L, and in Vo before provisional identification of the detected peaks. Non-processed MS data from the HS-ITEX-GC-TOF–MS analysis can be exported and then processed using our method for metabolite profiling (Fig. 2). However, the putative identification of the detected VOCs is limited because few libraries show the EI mass spectra and RI and because it is very difficult to obtain authentic standards for VOCs. Despite this limitation, we estimated how many mass spectra overlapped among Vo and the three commercially-available libraries for volatiles (Ad-L, Te-L, NIST05). The estimation procedure is clarified in details in the Materials and Methods section. Instead of complete matching of the compounds using CAS numbers and/or compound names, we used the similarity of each mass spectrum and the RI difference of the corresponding peak in the query library (Ad-L, Te-L, Vo) and the NIST05 reference library (Tables 2 and 3). Approximately 35 % of the mass spectra in Ad-L (3rd edition, 555/1607; 4th edition, 765/2205) exhibited high similarity against NIST05. On the other hand, only four compounds (β-maaliene, methyl tridecanoate, methyl undecanoate, and methyleugenol) showed a similarity value greater than 900 in Te-L; the SD of the RI differences of the four compounds was 5.6 (Table 3). Using their chemical structures in NIST05 we compared these compounds and found that they were identical in Ad-L and NIST05. The mass spectra in Te-L tend to be unique. Consequently, the difference shown in Table 3 may increase the number of compounds that can be annotated.

Fig. 2

Schema of the workflow for data processing and peak annotation to obtain the data matrix. Non-processed data for GC-TOF–MS analysis of each sample were exported as NetCDF files. These files were imported in MATLAB for baseline correction, peak alignment, and deconvolution by the H-MCR method. Libraries were prepared for the provisional identification of the extracted mass spectra of the VOC peaks (gray box). After merging the information into a data matrix, we obtained a data matrix comprised of the compound name, sample name, and the sum of the peak area of each extracted mass

Table 2

Libraries used for the provisional identification of VOCs

Library name	RI	Phase composition of the GC column(s)
Adams library (3rd ed.)	Available	5 % diphenyl, 95 % dimethyl polysiloxane
Adams library (4th ed.)	Available	5 % diphenyl, 95 % dimethyl polysiloxane
Terpenoids library	Available	100 % dimethyl polysiloxane
VocBinBase	Available	5 % diphenyl, 95 % dimethyl polysiloxane
NIST05	Available	Various types (polar and non-polar)

Table 3

Estimation of the number of similar compounds in the Adams (Ad-L) and the Terpenoids library (Te-L), and in VocBinBase (Vo) against NIST05

Library name	Number of EI spectra	Number of identified compounds	≥850a	SD of RI diff (≥850)	≥900a	SD of RI diff (≥900)
Ad-L (3rd ed.)	1607	1607	794	8.36	555	8.29
Ad-L (4th ed.)	2205	2205	1077	8.73	765	8.69
Te-L	1982	1982	91	8.48	4	5.56
Vob	1632	212	258	8.67	143	8.09
NIST05	190,825	163,198	–	–	–	–

The RI difference (diff) was calculated by subtracting the RI of a compound peak in the query library from that in the reference library (NIST05). The values were transformed into absolute values

SD standard deviation; RI diff absolute RI difference

aThe value represents similarity defined as described in “Methods”

bVocBinBase contains 1420 unidentified EI spectra

HS-VOC profiling of the 12 Allium plants using the ITEX technique

We conducted VOC profiling in the HS of 10 A. fistulosum-, one A. chinense-, and one A. tuberosum cultivars with the ITEX technique. The visual phenotypes of each Allium plant are presented in Fig. 1. We focused on the sheaths and basal plates to analyze the VOCs. The entire aerial parts of the other Allium cultivars used in this study are eaten in Japan. We obtained VOC profile data on 35 samples [three biological replicates except for A. fistulosum (class05, n = 2)] and 354 extracted mass spectral peaks as a data matrix. The detected peaks were identified or provisionally identified using our fully-automated annotation pipeline. Of these, 52 peaks, including two artifacts (Si-containing peaks derived from column breeding) were tentatively identified by comparing their mass spectra and the RI corresponding to those in the four libraries (Table 2), or identified using authentic standards (Additional file 3). The molecular formula of each annotated peak was investigated and the proportion of sulfur-containing peaks in the annotated peaks was calculated (Fig. 3). Approximately half of the annotated peaks contained sulfur atom(s) in their moieties. According to Pino et al. [49], sulfur-containing compounds account for approximately 90 % of the total volatile content in diethyl ether extracts of A. chinense and A. tuberosum. Our findings suggest that ITEX-based VOC profiling could detect not only sulfur-containing peaks but also other types of VOCs.

Fig. 3

Proportion of sulfur-containing peaks in the 52 annotated peaks, including two artifacts (Si-containing peaks), in the HS of Allium plants. The proportion was calculated by counting the number of annotated compounds that consisted of CHOS, CHO, CH, or CHOSi

We conducted principal component analysis (PCA) to visualize the similarities/differences in the VOC composition of each Allium cultivar (Fig. 4). The score scatter plot of the VOC profile data showed subspecies-dependent separations among A. chinense, A. fistulosum, and A. tuberosum (Fig. 4a). Next, we investigated the distribution of tentatively identified peaks in the profiles of the Allium cultivars. The PCA loading plot showed that, some peaks tended to be abundant in A. fistulosum cultivars [e.g. 3,4-dimethylthiophene (ID026)], while the levels of the two sulfur-containing compounds [2,5-thiophenedicarboxaldehyde (ID154) and diallyl disulphide (ID091)] were more abundant in A. tuberosum than in A. fistulosum cultivars (Fig. 4b).

Discriminative VOCs among the Allium cultivars

We compared the VOC profiles of each Allium cultivar to determine whether the VOC composition in the HS can be used in their differentiation. VOC changes in the HS of Allium samples were recorded by subtracting the average of the normalized responses of the annotated peaks (log2-transformed value) in each Allium cultivar from those of the control, Mikata spring onion (class01, Fig. 1b). The extent of the VOC changes tended to be similar to that shown by PCA (Fig. 4, Additional file 3). For example, the visual phenotype of the control cultivar Mikata spring onion (class01), and of Aoi-chan green spring onion (class02) was very similar (Fig. 1b, c). There was no significant difference in the level of the annotated VOCs between these cultivars (Additional file 3). In the VOC profiles of other cultivars of A. fistulosum, there were a few differences in the VOC levels when compared to the control (data not shown). Thus, we focused on the subspecies-dependent differences.

We compared changes in the level of the 50 annotated VOC peaks in the profiles of A. chinense and A. tuberosum against the control (class01) (false discovery rate, FDR < 0.05). Of these, The 15 compound peaks showed significant changes in the profiles of A. chinense, while the level of 36 peaks was changed in A. tuberosum (Fig. 5).

Fig. 5

The log2-fold changes in the VOCs of 50 annotated peaks in the VOC profiles of sulfur-containing peaks (a), peaks consisting of CHO (b), and peaks consisting of CH (c). The log2-fold changes (log2FC) in the normalized response of peaks of each cultivar (A. chinense or A. tuberosum) against that of the control cultivar class01 (A. fistulosum) are shown in the Additional file 3. We analyzed three biological replicates of each Allium plant except the white spring onion (class05, n = 2) cultivar. Asterisks on the top of the bars present that the level of VOCs in A. chinense and/or A. tuberosum was significantly changed using the LIMMA package (see the “Methods” section) when compared to the control cultivar (class01 in Table 1). The significance level was set at FDR <0.05 (asterisk *). Sulfur-containing peaks were: ID020 1-propene, 3,3’-thiobis, ID021 methanesulfonic acid, methyl ester, ID024 thiophene, 2,5-dimethyl, ID026 3,4-dimethylthiophene, ID028 thiophene, 2-ethenyl, ID032 thiophene, 2,4-dimethyl, ID054 dimethyl trisulfide, ID065 3-thiophenecarboxaldehyde, ID086 S-methyl methanethiosulphonate, ID091 diallyl disulphide, ID092 5-methyl-2-thiophenecarboxaldehyde, ID099 1,2,4-trithiolane, ID101 dipropyl disulfide, ID119 prop-1-enyl dithiopropanonate, ID124 thiophene, 2-propyl, ID134 3-vinyl-1,2-dithiacyclohex-5-ene, ID154 2,5-thiophenedicarboxaldehyde, ID164 trisulfide, di-2-propenyl, ID172 trisulfide, dipropyl, ID176 1,2,4-trithiolane, 3,5-diethyl, ID200 thieno[2,3-b]thiophene, 2-methyl, ID358 cyclic octaatomic sulfur, ID377 disulfide, methyl 1-propenyl, Peaks consisting of CHO: ID002 hexanal, ID011 2-butenal, 2-ethyl, ID012 2-pentenal, 2-methyl, ID050 heptenal, (2E), ID063 2-furanone, 2,5-dihydro-3,5-dimethyl; ID083 octen-1-al, (2E); ID088 2-octen-1-ol, (E); ID141 nonanoic acid; ID143 decenal, (2E); ID167 decadienal, (2E,4E); ID179 2-dodecenal, (E); ID201 3(2H)-furanone, 2-hexyl-5-methyl; ID218 2-tridecanone; ID259 3(2H)-furanone, 5-methyl-2-octyl; ID328 hexadecanoic acid methyl ester, n; ID334 1,2-benzenedicarboxylic acid, butyl octyl ester; ID354 9,12-octadecadienoic acid, methyl ester, (E,E); ID365 [1,1′,3′,1′-terphenyl]-2′-ol; ID371 2,5-cyclohexadiene-1,4-dione, 2,5-diphenyl; ID383 benzaldehyde; Peaks consisting of CH: ID139 benzene, 1,3-bis(1,1-dimethylethyl); ID266 8-heptadecene; ID273 heptadecane, n; ID279 1,1′-biphenyl, 2,2′,5,5′-tetramethyl; ID317 5-octadecene, (E); ID322 nonadecane; ID368 tricosane

Among 23 sulfur-containing peaks, 10 peaks showed significant changes in A. chinense, while 19 peaks were significantly changed in A. tuberosum (Fig. 5a). Of these, there were nine discriminant peaks in both subspecies. Thiosulfinates are the initial compounds in the HS of Allium species when their tissues are chopped or homogenated [50]; they decompose immediately and then sulfides are emitted as major aroma compounds [11, 42]. HS-ITEX VOC profiling detected a monosulfide (ID020), two disulfides (ID091 and ID101), and two trisulfides (ID054 and ID164). Their level was higher in A. tuberosum than in the control except for dipropyl disulfide (ID101) that is the major disulfide in Allium plants. The level of this compound was significantly lower in A. chinense and A. tuberosum than in the control (Fig. 5a). Interestingly, the level of this compound differed among the cultivars of A. fistulosum (Additional file 3). It suggests that this compound can be used for the chemotaxonomic classification of A.fistulosum cultivars. For bunching onions like A.fistulosum, DNA markers such as simple sequence repeats (SSRs), amplified fragment length- and single-nucleotide polymorphisms (AFLPs, SNPs) are available (http://www.vegmarks.nivot.affrc.go.jp/VegMarks/jsp/index.jsp). However, its high cost hampers the data collection of many cultivars. As a first step, our VOC profiling is useful for choosing representative cultivars in Allium plants for further analyses.

We detected nine peaks that are considered to be thiophene compounds [thiophene, 2,5-dimethyl, 3,4-dimethylthiophene, thiophene, 2,4-dimethyl, thiophene, 2,4-dimethyl, 3-thiophenecarboxaldehyde, thiophene, 2-propyl, 2,5-thiophenedicarboxaldehyde, 1,2,4-trithiolane, 3,5-diethyl, and thieno(2,3-b)thiophene, 2-methyl] in the Allium samples. Of these, two dimethylthiophenes (thiophene, 2,5-dimethyl, and 3,4-dimethylthiophene) were found to be thermal decomposition products from dialkyl disulfides in the distilled oils of Allium species [42, 51]. In the HS of Allium cultivars next to sulfur-containing volatiles there were VOCs, which consisted of CHO atoms. Out of those VOCs five compounds (hexanal, 2-pentenal, 2-methyl, heptenal, (2E), 2-tridecanone, and 3(2H)-furanone, 5-methyl-2-octyl) are previously found in distilled oils of A. fistulosum cultivars [52]. The level of four compound peaks such as hexanal, 2-butenal, 2-ethyl-, 2-tridecanone, and 3(2H)-furanone, 5-methyl-2-octyl was significantly different in the profiles of A. chinense, while 11 VOCs with asterisk (*) on the top of the bars were significantly changed in the profiles of A. tuberosum (Fig. 5b).

Annotated or identified compounds whose moieties included only CH atoms were categorized as alkanes or alkenes (Fig. 5c), out of which odd-numbered alkanes (heptadecane and nonadecane) are previously found in the methanol extract of garlic (A. sativum) [53] and in the HS of flowers of heliotrope and mandarin [54, 55]. The function(s) and biosynthetic pathway(s) of such compounds remain largely limited, except for dipropyl trisulphide in A. fistulosum and diallyl disulphide in A. tuberosum as described in [31]

Among 50 annotated peaks, 12 metabolite peaks showed significant changes in A. chinense and A. tuberosum when compared to the control cultivar, A. fistulosum (class 01). Dipropyl disulphide (ID101) known as the main VOC component in Allium plants was included in the 12 metabolite peaks. The 12 compounds as well as previously reported compounds were listed in Table 4 including thiosulfinates produced from S-alk(en)yl cysteine sulfoxides which are sulfur-storage compounds. These peaks are like to be used as discriminative compounds in VOC profiles of A. chinense, A. tuberosum, and A. fistulosum.

Table 4

List of previously reported compounds in Allium species and 12 VOCs with significant changes in A. chinense and A. tuberosum against A. fistulosum in this study

Name	ID	Detected Allim subspecies	Reference
Sulfur-containing peaks
Thiophene, 2,4-dimethyl	ID032	A. chinense, A. tuberosum and A. fistulosum
Dimethyl trisulfide	ID054	A. chinense, A. tuberosum and A. fistulosum	[42]
3-Thiophenecarboxaldehyde	ID065	A. chinense, A. tuberosum and A. fistulosum
Diallyl disulphide	ID091	A. chinense, A. tuberosum and A. fistulosum
1,2,4-Trithiolane	ID099	A. chinense, A. tuberosum and A. fistulosum
Dipropyl disulfide	ID101	A. chinense, A. tuberosum and A. fistulosum
3-Vinyl-1,2-dithiacyclohex-5-ene	ID134	A. chinense, A. tuberosum and A. fistulosum
2,5-Thiophenedicarboxaldehyde	ID154	A. chinense, A. tuberosum and A. fistulosum
Thieno[2,3-b]thiophene, 2-methyl	ID200	A. chinense, A. tuberosum and A. fistulosum
Fully saturated thiosulfinates	ND	A. cepa, A. sativum, A. ursinum, A. porrum, A. fistulosum, A. ascalonicum, A. ampeloprasum, A. schoenoprasum and A. tuberosum	[10]
Mono-S-b,c-unsaturated thiosulfinates	ND	A. cepa, A. sativum, A. ursinum, A. porrum, A. fi A. m, A, A. ascalonicum, A. ampeloprasum, A. schoenoprasum and A. tuberosum	[10]
Di-S-b, c-unsaturated thiosulfinates	ND	A. cepa, A. sativum, A. ursinum, A. porrum, A. fistulosum, A. ascalonicum, A. ampeloprasum, A. schoenoprasum and A. tuberosum	[10]
Mono-a, b-unsaturated thiosulfinates	ND	A. cepa, A. sativum, A. ursinum, A. porrum, A. fistulosum, A. ascalonicum, A. ampeloprasum, A. schoenoprasum and A. tuberosum	[10]
Mixed a, b- and c-unsaturated thiosulfinates.	ND	A. cepa, A. sativum, A. ursinum, A. porrum, A. fistulosum, A. ascalonicum, A. ampeloprasum, A. schoenoprasum and A. tuberosum	[10]
Allicin (diallylthiosulphinate)	ND	A. sativum	[5]
Alliin (S-allyl-l-cysteine sulphoxide)	ND	A. sativum,A. ursinum, A. ampeloprasum and A. longicuspis	[5, 10]
Dipropyl disulphide	ID101	A. fistulosum and A. tuberosum	[6, 42]
Dipropyl trisulphide	ND	A. fistulosum and A. tuberosum	[6]
1-Propenyl propyl disulphide	ND	A. fistulosum and A. tuberosum	[6]
Methiin (S-methyl-l-cysteine sulfoxide)	ND	A. cepa, A. sativum,A. chinens and A. longicuspis	[10]
Propiin (S-propyl-l-cysteine sulfoxide)	ND	A. cepa, A. porrum, A. porrum, A. altaicum and A. fistulosum	[10]
Isoalliin (S-propenyl-l-cysteine sulfoxide)	ND	A. cepa, A. nutans,A. ascalonicum and A. schoenoprasum	[10]
Ethiin (S-ethyl-l-cysteine sulfoxide)	ND	A. aflatunens, A. ampeloprasum, A. ochotense and A. victorialis	[10]
Butiin (S–n-butyl-l-cysteine sulfoxide)	ND	A. siculum	[10]
1-Propenyl-containing disulfides	ND	A. uictorialis	[13]
Thiopropanal S-oxide	ND	A. cepa	[42]
Propenyl propyl disulphide	ND	A. cepa	[42]
1-Propenyl propyl disulphide	ND	A. cepa	[42]
Di-1-propenyl disulphide	ND	A. cepa	[42]
Methyl propyl trisulphide	ND	A. cepa	[42]
Propenyl propyl trisulphide	ND	A. cepa	[42]
Peaks consisting of CHO
2-Butenal, 2-ethyl	ID011	A. chinense, A. tuberosum and A. fistulosum	[13]
Decenal, (2E)	ID143	A. chinense, A. tuberosum and A. fistulosum	[13]
2-Methyl-2-pentenal	ND	A. uictorialis and A. cepa	[13, 42]
Prop(en)yl aldehydes	ND	A. cepa	[42]
2-Methyl-2-pentenal	ND	A. cepa
Peaks consisting of CH
Nonadecane	ID322	A. chinense, A. tuberosum and A. fistulosum	[53–55]

ND not detected in this study

Conclusions

Since Allium plants emit various types of sulfur-containing compounds and other VOCs, comprehensive profiling techniques are needed. We developed a VOC profiling method for the HS of Allium samples that is based on an ITEX technique and our metabolomics pipeline by using GC-TOF–MS [56, 57]. The amount of sample material needed for HS collection was much lower with our ITEX-method than the traditional methods such as solvent extraction and steam distillation. Our findings suggest, that ITEX-based VOC profiling yields good reproducibility for the detection of various types of VOC in Allium plants. As ITEX-based VOC profiling captures differences in the composition of VOCs in the HS of Allium plants, it is probably appropriate for chemotaxonomic classification of these plants. For odor analysis of samples with strong odors, for example Allium plants, GC–olfactometry (GC–O) coupled with MS is useful because it facilitates the evaluation of odor compounds and yields MS spectral information. However, as the concentration of odor compound is often very low and odor-related VOCs can interact synergistically or additively, the identification of actual “odor” peaks remains difficult. Taken together, we think that the HS sampling- and the ITEX-based VOC profiling methods presented here help to improve the detection of odor compounds in Allium plants.

Methods

Chemicals

All chemicals and reagents used for this study were of spectrometric grade. The n-alkane standard solution C8–C20 for determination of RI was purchased from Fluka Chemical (Tokyo, Japan), deuterium-labeled alkanes used to distinguish natural alkanes collected from Allium samples were obtained from Cambridge Isotope Laboratories (Andver, USA), and dipropyl disulfide (98 %) and surrogate standard mixture (EPA524.2) from Sigma-Aldrich Japan (Tokyo, Japan). The other chemicals were purchased from Nacalai Tesque (Kyoto, Japan) or Wako Pure Chemical Industries (Osaka, Japan).

Plant material and sample preparation procedure

Metadata for this study are provided in Additional file 2.

Ten Allium (A.) fistulosum species, six spring onion cultivars, two scallions, and two Japanese-leek cultivars, rakkyo (A. chinense) and Oriental garlic (A. tuberosum), were purchased from a grocer in Kawasaki, Japan or harvested in a Japanese field (see Table 1 and Additional file 2). After removing the roots, a 10-cm length of the sheath and the basal plate of each plant sample were collected and chopped with stainless steel surgical blades (Feather, Tokyo, Japan). Out of the A. fistulosum cultivars, four were grown by applying a method (hilling) similar to that used for growing the leek A. ampeloprasum var. porrum to obtain longer white stems for consumption in Japan (Fig. 1f, g, h, l). Each sample was immediately frozen in liquid nitrogen and kept at −80 °C until use. As the group of samples of Mikata spring onion (class01) was gathered center of the PCA score scatter plot (Fig. 4), this cultivar was chosen as the control.

The samples were crushed into powder (2 min at 4 °C) in a Mixer Mill MM 311 instrument featuring a grinding jar with a stainless steel screw cap (Restech, Tokyo, Japan) and the frozen powder from each sample (flesh weight, 1 g) was weighed in a 20-ml HS vial (Supelco, MO, USA). For VOC profiling of Allium plants we used a modified method of Tikunov et al. [23] and Kusano et al. [31]. Briefly, the 20-ml HS-GC vial (Supelco) containing the frozen powder was closed with a magnetic screw cap (AMR, Tokyo, Japan) for ITEX- and SPME-analysis. Then, 1 ml of 100 mM 2,2′,2′’,2′’’-(ethane-1,2-diyldinitrilo) tetraacetic acid (EDTA)- NaOH water solution (pH 7.5) was added to each vial; the water derived from an Allium sample was considered to be equal to 1 ml. After vortexing, 10 μl of solution containing n-decane (d, 99 %; 50 μM), n-pentadecane (d, 98 %; 50 μM), n-eicosane (d, 98 %; 50 μM) for definition of RI and EPA524.2 fortification solution (20 μg/ml of fluorobenzene, 4-bromofluorobenzene, and 1,2-dichlorobenzene-d) as ISs was mixed in methanol, then solution was added to each vial as IS. Solid CaCl2 was added to obtain a final concentration of 5 M and the samples were stored overnight at 22 °C.

HS collection using the SPME fiber

The SPME device for a CTC CombiPAL auto-sampler (CTC Analytics, Zwingen, Switzerland) was purchased from AMR (Tokyo, Japan). We used an SPME fiber comprised of a 65-μm-thick layer of polydimethylsiloxane (PDMS)/divinylbenzene (DVB)-fused silica (FS) fiber/stainless-steel (SS) tube. Before analysis, the fiber was conditioned at 250 °C for 30 s in the injection port of an Agilent 6890 N gas chromatograph (Agilent Technologies, Wilmington, USA) equipped with a 30 m × 0.25 mm inner diameter fused-silica capillary column with a chemically bound 0.25-μl film Rtx-5 Sil MS stationary phase (RESTEK, Bellefonte, USA). Collection of volatiles was carried out by inserting the SPME-fiber to the vial and by trapping the VOCs for 20 min at 80 °C under continuous agitation. After HS collection it was placed in the injection port of the gas chromatograph that was coupled to a Pegasus III TOF mass spectrometer (LECO, St. Joseph, USA). The thermodesorption of VOCs occurred for 15 s at 250 °C.

HS collection using the ITEX device

We used a CTC CombiPAL auto-sampler (PAL COMBI-xt) featuring the ITEX device PAL ITEX-2 option (CTC Analytics). The ITEX procedure was controlled with a PAL Cycle Composer (CTC analytics). We conducted preliminary experiments to choose an appropriate sorbent material from the four materials, Tenax TA, Tenax GR (TGR), Carbosieve SIII (CSIII) and mixed TGR and CSIII (TGR/CSIII), that are commercially available (data not shown). Then, we chose that the sorbent material for the ITEX-2 portion was TGR (80/100 mesh)/CSIII (60/80 mesh). The parameters for HS collection were as described in the Additional file 2. After HS collection, 500 μl of the HS sample were injected into the injection port of the gas chromatograph coupled to the mass spectrometer used for HS collection by ITEX.

GC-TOF–MS analysis

GC-TOF–MS conditions were as described in the Additional file 2. Data acquisition was on a Pegasus III TOF mass spectrometer (LECO); the acquisition rate was 30 spectra/s in the mass range of a mass-to-charge ratio of m/z = 30–550. Five ISs were used for data normalization.

Data analysis

Raw data were exported in the network common data form (NetCDF) file format using LECO ChromaTOF software (version 2.32) and then processed with the hierarchical multi-curve resolution (H-MCR) method [38]. We obtained the normalized response for calculating the signal intensity of each metabolite from the mass-detector response by using the cross-contribution compensating multiple standard normalization (CCMN) method [44]. The resolved mass spectra were matched against reference mass spectra in the NIST-L (version NIST05) using NIST MS search program (version 2.0, http://www.chemdata.nist.gov/dokuwiki/doku.php?id=chemdata:ms-search). Peaks were tentatively identified according to the guidelines for metabolite identification [58]. When mass spectra exhibited a match value greater than 799 and the corresponding peaks had RIs with small differences upon comparison of their resolved mass spectra and RIs against those in the reference libraries (Ad-L, 3rd and 4th edition, and Te-L) and against Vo and NIST-L (see Table 2 and “Results and discussion” section), the peaks were considered to be putatively annotated compounds. We compared the RIs of sulfur-containing metabolites and compounds we detected with those reported in the literature [51, 52, 59].

To estimate the number of compounds that overlapped with each reference library, we first exported the mass spectral information, including the compound name, RI, synonyms, and m/z value, with relative peak intensity (maximum, 999; minimum, one) from each library in ASCII text format (.MSP) automatically. Then we compared the similarity of the mass spectra in each library using MassBank [60]. The similarity (≥850 or 900) and the RI difference (<|30 unit|) were used to extract the same or very similar compounds from the query library and NIST05. It should be noted that the standard deviation (SD) of the absolute RI difference of these compounds is less than 8.8 when we applied similarity of ≥850 (Table 3).

The lower limit of quantification (LLOQ) and the limit of detection (LOD) of dipropyl disulfide obtained from ITEX-GC-TOF–MS- and SPME-GC-TOF–MS analyses were estimated as described in the Additional file 2.

Statistical analysis

Multivariate analysis was performed with SIMCA-P + 12.0 software (Umetrics AB, Umeå, Sweden). For our analysis, profile data were log10-transformed, centered, and scaled to unit variance.

Log2-transformed profile data were statistically analyzed using the LIMMA package [61]. It includes FDR correction for multiple testing [62] in the R environment for statistical computing (version 2.14.1 for 64-bit).