Study on essential oils from four species of Zhishi with gas chromatography–mass spectrometry

Abstract

Background

Citrus fruits are widely used as food and or for medicinal purposes, and they contain a host of active substances that contribute to health. The immature fruits of Citrus sinensis Osbeck and its cultivars (CS), C. junos Sieb. ex Tanaka (CJ), C. aurantium L. and its cultivars (CA) and Poncirus trifoliate Raf. (PT) are the most commonly used medicinal herbs in Traditional Chinese Medicine, called Zhishi. And their mature fruits can be used as food.

Results

In this study, the essential oils of four different Zhishi species were extracted by steam distillation and detected using gas chromatography- mass spectrometry (GC-MS). A total of 39 volatiles from the four species were tentatively identified. The limonene was the most abundant amongst the four species. Principal component analysis (PCA) of essential oils showed a clear separation of volatiles among CS, CJ and PT. However, CA could not be separated from these three species. Additionally, the volatiles accounting for the variations among the widely separated species were characterized through their corresponding loading weight.

Conclusion

Sesquiterpenes were identified as characteristic markers for PT. The content of some monoterpenes could be as taxonomic markers between CS and CJ. This work is of great importance for the evaluation and authentication of Zhishi samples through essential oils.

Background

Citrus genus is the most important fruit tree crop in the world, with an annual production of approximately 102 million tons [1]. Citrus fruits are used in the food, cosmetic and pharmaceutical industry. They are consumed fresh and processed, as juices, jam, jellies, molasses, etc.. The fresh fruits are important starting materials for juice production [2]. Especially, chemical industry extracts from Citrus bioactive compounds like flavonoids, vitamins, dietary fiber and essential oils, etc. have beneficial effects on human health. The essential oils of Citrus fruits are widely used as flavoring in foods, perfumes and pharmaceutical formulations due to their functional properties like antimicrobial, antifungal, as well as analgesics, cough suppressants, and expectorants for eliminating phlegm [3–8]. Citrus oils extraction methods include steam distillation [9], cold pressing [10] and static headspace solid-phase microextraction [11].

The immature fruits of CA, CS, CJ and PT can be used as Zhishi, one of the primary traditional Chinese medicinal plants in China [12, 13]. The herbs from different species and cultivation locations are widely used in clinical applications. Such herbs have been used to treat stuffiness, intestinal fullness sensation and distending pain, diarrhea and dysentery in addition to the retention and preservation of food [14]. However, the efficacy of these species of Zhishi for treating digestive disturbances varies due to differences in the types and quantities of the chemical substances they contain [15, 16]. Our previous studies based on four flavanone compounds demonstrate that compounds from different species of Zhishi, such as CA, CS and CJ, can be distinguished by established differential modes of acquisition [17]. In addition to flavanones, essential oils are an important type of pharmacodynamic substance in Zhishi. Therefore, it is important to elucidate the differences between volatile substances among different species. Numerous studies [18–23] on the essential oils of Zhishi have been conducted using GC and GC-MS. However, because there is still a lack of information regarding the essential oils of different species of Zhishi, a comparative study is necessary.

Therefore, in this study, 75 authentic Zhishi samples from four species of Citrus fruit (CA, CS, CJ and PT) were collected. The essential oils were obtained by using steam distillation, and then were separated and identified by GC-MS. In total, 39 volatiles were tentatively identified and quantified, and among them were 26 monoterpenes, 8 sesquiterpenes, and 5 esters, aliphatic alcohols and aliphatic aldehydes. The details were shown in the Additional file 1. Then, a multivariate statistical analysis method was employed to elucidate the variation of volatile substances among different species of Zhishi. This is the first report on comparison study of essential oils from these four Zhishi species.

Results and discussion

Quantification of the volatile components

In total, 39 compounds were tentatively identified, constituting approximately 85.14-100% of the entire volatile concentration (Table 1). Each compound was identified by matching its retention characteristics and MS fragmentation patterns against standards and library spectra. The results are expressed as relative weight percentages calculated from peak areas. Volatiles were categorized into five chemical groups consisting of monoterpenes, sesquiterpenes, esters, aliphatic alcohols and aliphatic aldehydes. For all of the four Citrus species, the total content of essential oils primarily consisted of 26 monoterpenes, which accounted for 52.17-100% of the entire volatile concentration. The other four chemical groups – sesquiterpenes, esters, aliphatic alcohols and aliphatic aldehydes – contained in CA, CS and CJ had minor contributions to the entire volatile concentration accounting for only 0–2.28%. Sesquiterpenes were present in the oils analyzed in the greatest amount (3.8-34.28%) in the PT species compared to the other three species. This result is not surprising as sesquiterpenes hydrocarbons are important in the characteristic aroma of many kinds of Citrus fruits [24]. Alpha and β-pinene are chemically unstable bicyclic terpenes due to a strained four membered ring. Thus, they are found at low levels and evenly distributed in all the four species.

Table 1 Identifications of Zhishi volatiles and their relative total ion current peak area from four different species

GC-MS analyses of the four different species of Zhishi samples

In the 35 samples of CA, monoterpenes produced the highest percentage of volatiles (84.67-99.61%). Among the volatiles, the most predominant was limonene (27–82.84%, Mean 63.47%), followed by γ-terpinene (0.25-32.31%, Mean 12.01%), β-linalool (0–26.15%, Mean 9.10%), β-pinene (0.4-17.01%, Mean 2.59%), β-cis-ocimene (0–9.99%, Mean 2.21%) and (+)-sabinene (0.1-8.26%, Mean 1.77%). Meanwhile, the mean values of volatiles such as α-pinene, β-myrcene, p-mentha-1,4(8)-diene, o-cymene, (−)-terpinene-4-ol and α-terpieol were higher than 0.5%, with an average percentage of total volatile concentration of 1.25%, 1.05%, 0.66%, 0.63%, 0.61% and 0.60%, respectively. Other chemicals, such as germacrene D, α-methyl-α-[4-methyl-3-pentenyl]-oxiranemethanol and n-octanal, belonging to sesquiterpenes, aliphatic alcohols and aliphatic aldehydes chemical groups, were minor volatile components in CA species, averaging 0.26%, 0.35% and 0.33% of the total percentage, respectively.

For CS species, monoterpenes were also the highest concentration of the volatile, averaging 85.3-100% of the total percentage. Limonene was the most abundant monoterpene in this species accounting for 72.75-89.35%, with an average value at 82.25%. Additionally, (+)-sabinene and β-linalool were recorded at a relatively lower proportion, 3.61-12.82% and 1.34-9.21%, respectively. Other compounds within the monoterpenes group, such as β-myrcene, γ-terpinene, (−)-terpinene-4-ol, β and α-pinene, also occurred at a higher proportion, the mean amounts were 1.14%, 0.91%, 0.73%, 0.62% and 0.50%, respectively. Compounds of other chemical groups had a minor amount, and the mean values were lower than 0.1%.

For CJ species, the total amounts of monoterpene compounds ranged from 89.31-99.81%. The most two abundant essentials, limonene and γ-terpinene accounted for 35.53-60.60% (Mean 46.47%) and 21.71-46.74% (Mean 34.16), respectively. Other monoterpenes, such as β-linalool, β-pinene, α-pinene, p-mentha-1,4(8)-diene, o-cymene, α-thujene, β-myrcene, (+)-4-carene, α-terpieol, p-cymen-2-ol and (+)-sabinene, each had average amounts higher than 0.5% at 3.53%, 2.80%, 2.61%, 1.63%, 1.51%, 1.12%, 0.94%, 0.88%, 0.60%, 0.55% and 0.54%, respectively. Except for germacrene D, which had an average value at 0.15% and belongs to the sesquiterpenes, compounds of other chemical groups were less than 0.1%.

Interestingly, in the PT species, sesquiterpenes (3.80-34.28%) were an important chemical group of character impact volatiles because of the relatively high percentage in the overall volatile composition compared to the other three species. In addition, it serves a characteristic aroma for many types of Citrus fruits. Among these sesquiterpenes, the most predominant was caryophyllen (Mean 5.32%), followed by γ-elemene (Mean 3.81%), β-farnesene (Mean 2.80%) and germacrene D (Mean 1.83%). Limonene and β-myrcene were the major monoterpenes in this species as the average value accounted for up to 48.14%and 18.29%, respectively.

Principal component analysis

The variability of essential oils in the four Citrus species of Zhishi shown in Table 1 was not anticipated. Principal component analysis (PCA) was employed to examine this large data set without having to assign classifications before analysis. PCA is an unsupervised clustering method that requires little prior knowledge of the data set and acts to reduce the dimensionality of multivariate data without losing important information. Par-scaled (scaled to square root of standard deviation) mathematical methods were performed to pre-treat the data set resulting from the volatiles of different Citrus species. The multivariate statistical analysis methodology helps to identify inherent patterns in the data in an unbiased manner and highlights the similarities and differences amongst samples. This methodology also helps identify those essential oils that are the most representative within the entire data set.

PCA score plot

As shown in Figure 1, each coordinate represents a sample, and the separation of the four different species of Zhishi samples was observed in the PCA scores plot. A two-component PCA model cumulatively accounted for 54.5% of the total variance. The CA species could not be separated from the other three species; however, according to our previous studies [17], CS, CJ and CA could be differentiated from one other by choosing flavanones as marker compounds. The high degree of clustering and minimal overlap of CS, CJ and PT suggest that these three species have volatile profiles that are unique to each cultivar.

Figure 1

PCA scores plot (PC1 vs. PC2) of the four species of Zhishi samples: 35 samples from CA , 15 samples from CJ, 21 samples from CS and 4 samples from PT . (CA = (red star) CJ = (green triangle), CS = (black diamond) and (PT = blue circle)).

These three species, CS, CJ, and PT, were further analyzed by PCA. We observed three widely separated clusters in Figure 2 consisting of CS, CJ and PT. The PT samples were clearly separated from CS and CJ by principal component 2 (PC2), whereas the CJ and CS samples were clearly separated by principal component 1 (PC1). The results confirmed that CS, CJ and PT were different regarding the levels and occurrence of these essential oils. However, the volatiles in CA could not be distinguished from the other three species.

Figure 2

PCA scores plot (PC1 vs. PC2) of the three species of Zhishi samples: 15 samples from CJ , 21 samples from CS and 4 samples from PT . (CJ = (red star), CS = (green triangle) and PT = (black diamond)).

PCA loading plot

The corresponding PCA loading plot was utilized to identify the differential volatiles accountable for the separation among CS, CJ and PT. Although some information of low intensity might be lost due to Par scaling, the loading plot clearly facilitated the profiling of the volatiles in our experiments. The volatiles that accounted for maximum variance in the data set are given more weight or loading. As shown in Figure 3, the variation within PC2 in the loading plot was due to the marker volatiles that accounted for the differences between PT and the others as observed in the scores plot. Analogously, the distribution in PC1 of the loadings plot demonstrated the variation of CS and CJ, respectively. For example, sesquiterpenes, such as germacrene D (0.34 on PC2), caryophyllen (0.33 on PC2), β-myrcene (0.32 on PC2), β-elemene (0.31 on PC2) and α-farnesene (0.27 on PC2) were the most differentiating volatiles, which played an important role in differentiating PT from the CS and CJ. Meanwhile, in PC1, monoterpene compounds, such as p-mentha-1, 4(8)-diene (loading of 0.29), limonene (0.28), (+)-4-carene (0.28), β-pinene (0.28), γ-terpinene (0.28), α-thujene (0.27), p-cymen-2-ol (0.26) and α-pinene (0.26), appeared to be of some relevance to the observed variability between CS and CJ.

Figure 3

Overall PCA loading plot of the GC-MS signals associated with the three species of Zhishi samples: 15 samples from CJ , 21 samples from CJ and 4 samples from PT .

Discussion

Literature reveals that limonene and γ-terpinene were the major components of the Citrus oils obtained by steam distillation and cold pressing [23]. The essential oils of four different species of Zhishi samples have been well studied and the major chemical groups are similar to the previous studies [21–23]. Limonene as a major constituent in several citrus oils (orange, lemon, mandarin, lime and grapefruit) is identified as the most predominant compound in the four species of Zhishi samples. The chemical compositions and their content of the essential oils are varied in different species. Especially, sesquiterpenes were present in the oils analyzed in the greatest amount in the PT species compared to the other three species. The contents of some monoterpene compounds such as p-mentha-1,4(8)-diene, limonene, (+)-4-carene, β-pinene, γ-terpinene, α-thujene, p-cymen-2-ol and α-pinene could be used as taxonomic markers between CS and CJ. The volatiles in CA could not be distinguished from the other three species. This may be because that CA samples distributed in a wide regions in China (Jiangxi, Hunan, Sichuan and Zhejiang). In the future, the active constituents contained in Zhishi samples from different regions need to be further investigated. Considering that all of the samples were collected in May or June under the same conditions impact of the growth stage factors on essential oil composition variability has been excluded.

Experimental

Plant material

All of the 75 Zhishi samples were collected from the main Zhishi-producing provinces of China: Jiangxi, Sichuan, Zhejiang, Fujian and Guizhou in May or June during a three-year period. For each of the 75 Zhishi samples, the species, local name, collection location, year of collection and growing environment are listed in Table 2. The diameters of collected Zhishi samples were limited to 0.5 cm ~ 2.5 cm, which is consistent with the requirements of the Chinese Pharmacopoeia (Chinese Pharmacopoeia, 2010). Professors Ge Fei, Yan Zhuyun, Zhang Yungui, Xu Jianguo and Ke Fuzhi identified them as genuine samples of C. aurantium L. and its cultivars (CA), C. sinensis Osbeck and its cultivars (CS), C. junos Sieb. ex Tanaka (CJ) and Poncirus trifoliate Raf. (PT). The dried specimens (marked as CA-1 ~ CA-35, CS-1 ~ CS-21, CJ-1 ~ CJ-15 and PT-1 ~ PT-4) were deposited at the Institute of Basic Theory, China Academy of Chinese Medical Sciences, Beijing, P. R. China. These four PT samples were wild growing.

Table 2 The origins of the 75 Zhishi samples collected

Extraction of essential oils

Fresh, immature fruits of Zhishi were cut and dried in a manner similar to that described in the Pharmacopoeia of China [14]. The dried fruits were soaked in water for 2 hours and steam-distilled in a clevenger-type apparatus for 8 hours. This procedure was repeated 6 times. The oils were then dried over anhydrous sodium sulfate, diluted by acetone, and stored at 4–6°C before analysis.

GC–MS analysis

GC–MS analysis was performed on GCMS-QP2010 Plus (Shimadzu, Kyoto) equipped with a capillary column (Rxi-5 ms, 30 m × 0.25 mm, 0.25 μm). Helium was used as the carrier gas at a flow rate of 1.0 mL/min. Oven temperature was varied from 50°C (1 min held) to 140°C(1 min held) at 5°C/min, and then from 140°C to 200°C (3 min held) at 10°C/min. The injector and interface temperatures were held at 250°C. Mass spectra in the electron impact mode (EI-MS) were generated at 70 eV. The ion source temperature was held at 250°C. The scanning time was 0.5 sec over a range of m/z 45–450. An oil sample of 1 μL was injected in the split mode injection (split ratio, 60:1).

Identification and quantitative determination

The volatile components were tentatively identified based on linear retention index (RI) and by the comparison of mass spectra with MS data of reference compounds. The linear retention indices were determined for all constituents by using a homologous series of n-alkanes (C₈–C₂₀). The components were identified by comparison of their mass spectra with those of the NIST05 and NIST05S mass spectral library and further confirmed by comparison with the literature values [19, 25]. The relative proportions of the essential oil constituents were expressed as percentages obtained by peak area normalization with all relative response factors being taken as one. The samples were run in triplicate.

Statistical analysis

The GC-MS data of different species of Zhishi samples were analyzed to identify potential discriminant variables. Multivariate statistical analyses, including unsupervised PCA and supervised partial least squares-discriminant analysis (PLS-DA), were performed using the SAS 9.1.3 statistical package (order no. 195557). PCA was used to observe the natural interrelationship among the chemical components for each of the four Citrus species. Analysis of variance was employed to identify those volatiles that would be most differentiating among the four species. Furthermore, the corresponding loading plot was used to interpret the variations among the samples. The critical p value for all analyses in this study was set to 0.05.

Conclusion

GC-MS coupled with multivariate statistical analysis of volatiles from different Citrus species, revealed a diverse volatile distribution among CS, CJ and PT. However, the volatiles contained in CA exhibited less specificity; therefore, they could not be separated from the other three species. The PCA scores plot of CS, CJ and PT appeared as widely distinct clustering. The chemical compounds accountable for the different volatile profiles of the three species are germacrene D, caryophyllen, β-myrcene, β-elemene, α-farnesene, p-mentha-1, 4(8)-diene, limonene, (+)-4-carene, β-pinene, γ-terpinene, α-thujene, p-cymen-2-ol and α-pinene, all of which are significantly up- or down- regulated in the different three species. In short, sesquiterpene compounds could serve as characteristic markers for PT. The content of some monoterpene compounds could be used as taxonomic markers between CS and CJ. This work demonstrates the variations of essential oils among different species of Citrus fruits and is of great significance in the pharmacological and clinical use of Zhishi.