Skip to main content

Variation of chemical composition of essential oils in wild populations of Thymus algeriensis Boiss. et Reut., a North African endemic Species

Abstract

Background

Thymus algeriensis is an endemic aromatic plant to Tunisia largely used in folk medicine and as a culinary herb. The bulks aromatic plants come from wild populations whose essential oils compositions as well as their biological properties are severely affected by the geographical location and the phase of the plant development. Therefore, the aim of the present work is to provide more information on the variation of essential oil composition of T. algeriensis collected during the vegetative and the flowering phases and from eight different geographical regions. Besides, influence of population location and phenological stage on yield and metal chelating activity of essential oils is also assessed.

Methods

The essential oil composition of Thymus algeriensis was determined mainly by GC/FID and GC/MS. The chemical differentiation among populations performed on all compounds was assessed by linear discriminate analysis and cluster analysis based on Euclidean distance.

Results

A total of 71 compounds, representing 88.99 to 99.76% of the total oil, were identified. A significant effect of the population location on the chemical composition variability of T. algeriensis oil was observed. Only 18 out of 71 compounds showed a statistically significant variation among population locations and phenological stages. Chemical differentiation among populations was high. Minor compounds play an important role to distinguish between chemical groups. Five chemotypes according to the major compounds have been distinguished. Chemotypes distribution is linked to the population location and not to bioclimate, indicating that local selective environmental factors acted on the chemotype diversity.

Conclusions

The major compounds at the species level were α-pinene (7.41-13.94%), 1,8-cineole (7.55-22.07%), cis-sabinene hydrate (0.10-12.95%), camphor (6.8-19.93%), 4-terpineol (1.55-11.86%), terpenyl acetate (0-14.92%) and viridiflorol (0-11.49%). Based on major compounds, the populations were represented by (α-pinene/1,8-cineole/cis-sabinene hydrate/camphor/viridiflorol), (1,8-cineole/camphor/terpenyl acetate), (α-pinene/1,8-cineole/camphor), (1,8-cineole/camphor/4-terpineol) and (α-pinene/1,8-cineole/cis-sabinene hydrate/camphor/4-terpineol) chemotypes. Variation of phenological stage did not have a statistically significant effect on the yield and metal chelating activity of the essential oil. These results can be used to investigate the geographical location and the harvesting time of this plant for relevant industries.

Background

In the last few years, there has been an increasing concern regarding the safety and potentially adverse effects of synthetic chemicals used for food preservation or in medicine. Therefore, the commercial development of medicinal plants as new sources of bioactive products to enhance human health and food preservation is of prime importance. Essential oils extracted by hydrodistillation from aromatic plants are needed for their various biological and pharmacological properties. However, several factors, namely climatic, geographic conditions and growth stage of collected plants may severely affect essential oil yield, their composition and their biological properties. Thus, studies of chemical variability of essential oil in relation to environmental factors might provide information on what determines its chemical polymorphism. In addition, knowledge of the chemical composition of essential oils is a very important quality criterion for their marketing and contributes to their valorization.

Thymus (Lamiaceae) is a large genus divided in eight sections, comprising more than 250 species particularly prevalent in the Mediterranean area. Thymus algeriensis Boiss. et Reut., which is endemic to Tunisia and Algeria, is an herbaceous fragrant plant largely used, fresh or dried, as a culinary herb [1]. Furthermore, this plant is also widely used in folk medicine against illnesses of the digestive tube and antiabortion [2]. Recently, the T. algeriensis essential oil was found to possess an interesting inhibitory activity towards angiotensin I-converting enzyme suggesting the potential of this plant as an antihypertensive agent [3]. In Tunisia, T. algeriensis populations are distributed from the sub-humid to the lower arid bioclimates at altitudes ranging from 120 to 1100 m. The species grows on poor fertile calcareous soils and occurs in scattered and small populations. T. algeriensis is a short lived, diploid (2n = 2x = 30) and gynodioecious shrub. It reproduces by seeds and can reach 20-50 cm in height. The leaves are opposite and linear/lanceolate (6-12 mm). The flowers, with ovate bracts and pink purplish or whitish purple corolla, are small (5-7 mm). Flowering takes place between April and June.

Previous works on T. algeriensis showed important intraspecific chemical variability of the essential oils among samples according to the geographical regions [36]. However, there are no researches in assessment of essential oil variations at the vegetative stage of this plant and in different geographical locations. Therefore, the aim of the present work is to provide more information on the variation of volatiles of Tunisian T. algeriensis collected during the vegetative and the flowering phases and from eight different localities and to determine in which way this would affect the corresponding oils yields and their metal chelating activities.

Results and discussion

Identified essential oil compounds

Eight wild populations of T. algeriensis from different regions were collected during the vegetative (S1) and the flowering (S2) stages. They belonged to 3 bioclimatic zones. Populations 1, 2 and 3 were located at the South West of Tunisia (Gafsa region) with an inferior arid climate characterized by a mean rainfall of 100-200 mm/year, while populations 5, 6, 7 and 8 were localized at the North West of Tunisia (Le Kef region) characterized by a superior semi-arid climate (rainfall: 400-500 mm/year). The intermediate population 4 was from the inferior semi-arid bioclimate and characterized by a mean rainfall of 300-400 mm/year. The altitudes ranged from 192 m (population 3) to 800 m (population 7) (Table 1 and Figure 1).

Table 1 Location and main ecological factors of the 8 T. algeriensis populations analyzed
Figure 1
figure 1

Goegraphical localization of the 8 sites of Tunisian T. algeriensis populations. For the detailed description of the locations: see Table 1.

The chemical composition of all the oil samples was mainly investigated using both GC/FID and GC/MS techniques. The percentages and the retention indices of the identified compounds of these essential oils were listed in Table 2 in the order of their elution on the HP-5MS column. Seventy-one compounds, representing 88.99 to 99.76% of the total essential oil, were identified and separated on the basis of their chemical structures into 5 classes (Table 2). Whatever the phenological stage, all these essential oils were characterized by very high percentage of monoterpenes (49.91-90.33%) and especially the oxygenated ones (32.01-62.18%) which constituted the predominant class as was found previously for T. algeriensis[3, 5, 6]. The sesquiterpenes were also represented mainly by oxygenated sesquiterpenes (2.92-21.84%) in contrast to what has been observed by Ben El Hadj Ali et al. [6] where the amount of oxygenated sesquiterpenes did not exceed 4.6% of the total essential oil of T. algeriensis. Although all the studied samples could be classified as oxygenated monoterpene-rich oils, they have shown wide range of variations in their compounds. The essential oils chemotypes detected in Tunisian T. algeriensis populations based on major compounds as well as the clusters among populations based on all the oil compounds are shown below.

Table 2 Mean percentage of compounds (%) in essential oils of the 8 T. algeriensis populations during the vegetative (S1) and the flowering (S2) stages

Chemical variation according to population locations and phenological stages

Our results showed that Fisher test was not applicable for 21 essential oil compounds (non-Normal distributed variables) (P1, Table 2). The analysis of variance showed that the means for the majority (47 over 50) of the oil compounds differed significantly among populations (P1, Table 2) and not between phenological stages (p > 0.05). In fact, only 11 compounds (sabinene 6, β-myrcene 8, trans-β-ocimene 14, cis-sabinene hydrate 16, camphenilone 18, p-cymen-8-ol 31, thymyl methyl ether 37, γ-cadinene 59, elemol 61, spathulenol 64 and caryophyllene oxide 65) differed significantly (p < 0.05) between vegetative and flowering stages (data not shown). According to Fisher test, 9 over 50 compounds differed significantly among populations and between phenological stages, which are: sabinene 6, cis-sabinene hydrate 16, camphenilone 18, p-cymen-8-ol 31, thymyl methyl ether 37, γ-cadinene 59, elemol 61, spathulenol 64 and caryophyllene oxide 65.

Among the 21 essential oil compounds (non-Normal distributed variables) non-parametric statistical tests such as Friedman (two-way analysis of variance) and Kruskal-Wallis (one-way analysis of variance) tests were applied. The resulting p-values of Friedman test, which are global p-values reflecting the combined effect of population location and phenological stage, were found to be significant at p < 0.1 for 8 over 21 compounds (linalyl acetate 40, cuminol 42, p-mentha-1,4-dien-7-ol 45, terpenyl acetate 46, α-gurjunene 51, α-humulene 53, alloaromadendrene 56 and palustrol 62) (P2, Table 2). Moreover, according to the Kruskal- Wallis test, the analysis of variance showed that means of 17 over 21 compounds differed significantly at p < 0.05 among populations (P3, Table 2). Besides, according to the same test, only 5 over 21 compounds (trans-piperitol 34, cuminol 42, carvacrol 44, α-gurjunene 51 and γ-eudesmol 68) differed significantly at p < 0.05 between vegetative and flowering stages (data not shown). Therefore, according to non-parametric statistical tests, 9 over 21 compounds (linalyl acetate 40, cuminol 42, carvacrol 44, p-mentha-1,4-dien-7-ol 45, terpenyl acetate 46, α-gurjunene 51, α-humulene 53, alloaromadendrene 56 and palustrol 62) differed significantly among populations and between phenological stages.

By combining all tests, the amounts of only 18 over 71 compounds (sabinene 6, cis-sabinene hydrate 16, camphenilone 18, p-cymen-8-ol 31, thymyl methyl ether 37, linalyl acetate 40, cuminol 42, carvacrol 44, p-mentha-1,4-dien-7-ol 45, terpenyl acetate 46, α-gurjunene 51, α- humulene 53, alloaromadendrene 56, γ-cadinene 59, elemol 61, palustrol 62, spathulenol 64 and caryophyllene oxide 65) showed a statistically significant variation among population locations and phenological stages (Table 2). Furthermore, taking into account all the oil identified compounds, a general model (Wilks' Lambda, two way analysis) was applied. In fact, a significant effect of the population location was also observed at p < 0.001. Nevertheless, neither the phenological stage, nor the interaction between population location and phenological stage were found to be statistically significant on the chemical composition of T. algeriensis essential oil (p > 0.05) (data not shown).

Table 2 showed that α-pinene 3, 1,8-cineole 13 and camphor 26 were the major compounds present in most populations belonging to different bioclimates. α-Pinene 3 was highly (7.41-13.94%) represented in most populations. Percentages of 1,8-cineole 13 ranged from 7.55% (population 2) to 22.07% (population 7). The amounts of camphor 26 ranged from 6.8% (population 2) to 19.93% (population 6), whereas Ben ElHadj Ali et al. [6] showed that camphor characterized few populations from the semi-arid zone of Tunisia. In addition to the compounds already described, 4-terpineol 30 mainly characterized the inferior arid zone population 3 (8.34-11.86%) (Table 2). In our recent work [3], 4-terpineol was also found at relatively high rate (7.36%) in a population from the inferior arid zone. Nevertheless, the amounts of 4-terpineol were very low (< 0.1%) in other described Tunisian populations from the same bioclimatic zone [6]. Our results also showed that cis-sabinene hydrate 16 characterized only population 2 (9.86%) and population 3 (12.95%) from the inferior arid zone at the vegetative stage (Table 2). Recently, cis-sabinene hydrate was also found at relatively high rate (5.29%) in a population from the same region at the flowering stage [3]. However, this compound was not detected in Tunisian populations described by Ben ElHadj Ali et al. [6] and it was reported at low amounts (0.15-2.30%) in Algerian populations [4, 5]. The amounts of terpenyl acetate 46 were very low (< 3.22%), except for population 4 (inferior semi-arid bioclimate), which was distinguished by a high proportion of this constituent (8.88-14.92%) (Table 2). Interestingly, terpenyl acetate 46, whose identification was confirmed by 13C-NMR spectroscopy, was described for the first time as a main compound in the essential oil of Tunisian T. algeriensis. Viridiflorol 66 has the highest percentage in the inferior arid population 2 (5.69-11.49%) (Table 2). A similar result was previously described by Ben ElHadj Ali et al. [6] in populations from the same bioclimatic region. Nevertheless, this compound was not found in the essential oils of Algerian populations of T. algeriensis[5]. As can be seen in Table 2, in populations from the superior semi-arid bioclimate, borneol 29 was detected with low amounts (4.60-6.21%) as compared to what has been obtained by Ben ElHadj Ali et al. [6] where this compound ranged between 18.4 and 24.3%. Thymol 43 was absent in the most of populations (Table 2). By contrast, this compound was found to be highly represented (54.9%) in a population from the superior-arid zone of Tunisia [6] and the percentage of thymol in Algerian T. algeriensis ranged from 0.2 to 29.5% [5].

Chemical clusters among populations

To identify possible relationships between volatile compounds and geographical locations, linear discriminate analysis (LDA) was applied. The LDA, performed on average contents of all compounds for each population regardless the phenological stage, showed that the first two principal axes represented 77.90% of the total variation. The first axis (54% of the total variation) was mainly correlated with linalyl acetate 40, α-gurjunene 51, α-copaene 48, cuminol 42, thymol 43, eremophilene 58, elemol 61, thymyl methyl ether 37 and bicyclogermacrene 57. The second axis represented 24% of the total variation, and linalyl acetate 40, viridiflorol 66, terpenyl acetate 46, trans-piperitol 34, eremophilene 58 and alloaromadendrene 56 were the main compounds contributing to its definition. The plot of the projection of the average values of all the compounds onto the first two principal axes, revealed a high chemical dispersion among populations (Figure 2). Furthermore, when using only the 18 compounds which show a statistically significant variation among populations and phenological stages (Table 2), the plot according to axes 1 and 2 (70.70% of the total variation) showed a similar chemical population groups (data not shown). Therefore, according to the linear discriminate analysis, four population groups in relation to the geographic location could be distinguished. The first, the second and the third group represented by populations 2, 4 and 1 respectively, situated at the periphery of the plot. Population 2 (inferior arid zone), situated at the negative sides of axes 1 and 2, constituted the first group. Population 4 (inferior semi-arid zone), situated at the positive side of axis 1 and at the negative side of axis 2, formed the second group. Population 1 from the inferior arid zone is situated at the positive sides of axes 1 and 2 formed the third group. The fourth group, situated in the centre of axis 1 and 2, is represented by the populations 3 from the inferior arid bioclimate and populations 5, 6, 7 and 8 from the superior semi-arid bioclimate.

Figure 2
figure 2

Linear descriminant analysis (LDA) for the essential oil compounds of the 8 T. algeriensis populations. Projection of the average contents of the essential oil compounds onto the first two principal axes (+ and - indicate positive and negative correlations with the axes, respectively). Coding numbers of populations' locations and for the detailed description of the bioclimatic zones: see Table 1.

In addition to the linear discriminate analysis and to better characterize populations groups, cluster analysis (dendrogram) was applied to a matrix linking essential oil composition to sample location and phenological stages. In fact, the dendrogram generated from the Euclidean distances (Figure 3) performed on the essential oils compounds of T. algeriensis populations at each phenological stage, showed population groupings globally similar to those observed by the LDA clustering. In fact, the general structure of the dendrogram showed the existance of three main clusters. The first group included populations 1 and 3 from the inferior arid bioclimate which could be divided into two subgroups represented by population 1 and population 3, respectively. The dendrogram showed that population 3 had a greater affinity with the population 1. Moreover, LDA clustering also showed that the axis 2 divided populations into two major groups where population 3 trended with population 1 (Figure 2). The second group also contained two subgroups represented by populations (5, 6, 7S1 and 8) and populations (4 and 7S2), respectively, while the third cluster was represented by population 2.

Figure 3
figure 3

Dendrogram obtained by cluster analysis based on Euclidean distance performed on the essential oil compounds of the 8 T. algeriensis populations during two phenological stages. PiSj represent population i (i = 1-8) at the vegetative stage (S1) or at the flowering stage (S2). Coding numbers of populations' locations: see Table 1.

Principal axes of Figure 2 showed that essentially minor compounds played an important role to distinguish between the chemical groups. However, conventional essential oil chemotypes were determined only on the basis of major compounds and therefore, five T. algeriensis chemotypes could be distinguished. The first group was represented by population 2 characterized by oils rich in α-pinene 3/ 1,8-cineole 13/ camphor 26 and it is distinguished from other populations by the presence of cis-sabinene hydrate 16 at the vegetative stage and viridiflorol 66 at the flowering stage as major compounds. By contrast to other populations, population 4 was characterized by the highest percentage of terpenyl acetate 46 and other compounds (14.82-22.88) that did not belong to the monoterpenes and the sesquiterpenes (Table 2). Population 4 can be defined as 1,8-cineole 13/ camphor 26/terpenyl acetate 46 chemotype and formed the second group. Populations 1, 5, 6, 7/S1 and 8 can be defined as α-pinene 3/ 1,8-cineole 13/ camphor 26 chemotype and formed the third group. The forth group represented by populations 3 and 7 at the flowering stage (3/S2 and 7/S2) which corresponded to a chemotype rich in 1,8-cineole 13/camphor 26/4-terpineol 30. The fifth group represented by population 3 at the vegetative stage (3/S1) presented an essential oil rich in α-pinene 3/ 1,8-cineole 13/ cis- sabinene hydrate 16/ camphor 26/4-terpineol 30.

In our study, Tunisian T. algeriensis showed a high chemical diversity among populations from the same region and bioclimate. In fact, populations 1, 2 and 3 from the inferior arid bioclimate and which were geographically close populations, clustered separately into different chemotypes. Nevertheless, the northern populations 5, 6, 7 and 8 from the superior semi-arid bioclimate and which were geographically near each other constituted an homogeneous group (Figure 2 and 3). Ben ElHadj Ali et al. [6] also showed a high chemical polymorphism among T. algeriensis populations. They showed that distribution of essential oil chemotypes was not always concordant with the bioclimatic zones and seemed rather to be linked with the geographic location and local selective forces acting on the chemotype diversity. In fact, local abiotic (topography, moisture, temperature and edaphic factors) and/or biotic selective factors (associated fauna and flora) act on loci terpene biosynthesis pathways and contribute to the emergence of different chemical profiles [7].

Influence of population location and phenological stage on yield and metal chelating activity of essential oil

The essential oils extracted by hydrodistillation from the dried aerial parts of T. algeriensis, collected from diverse locations during the vegetative and the flowering stages, ranged from 1.03 to 3.66% (v/w) (Table 3). These yields were higher in South West of Tunisia (populations 1, 2 and 3 from the inferior arid bioclimate) than in North West of Tunisia (populations 5, 6, 7 and 8 from the superior semi-arid bioclimate) with a maximum obtained in the population 2. It was reported that climatic conditions, soil types of collected regions and different phases of the plant development induce high variations in essential oil yield and their compounds [8]. Besides, for populations 5 and 7 these yields significantly decreased (p < 0.01) from the vegetative to the flowering stage (Table 3). Similar results were previously obtained for Malva aegyptiaca[9] where the yields of volatiles decreased from vegetative stage, full-flowering plants to seed-bearing plants. Nevertheless, the essential oil yields were found to be higher during the flowering phase than in the vegetative stage of Thymus capitatus[10] or Thymus caramanicus[11].

Table 3 Yield and chelating activity of T. algeriensis essential oils during the vegetative and the flowering stages. Coding numbers of populations' locations: see Table 1

Metal chelating activity was known as one of antioxidant mechanisms, since it reduced the concentration of the catalyzing transition metal in lipid peroxidation. Among the transition metals, Fe2+ ion was known as the most important lipid oxidation prooxidant due to its high reactivity [12]. Chelating activity is presented by IC50 value, defined as the concentration of the essential oil needed to chelate 50% of Fe2+ present in the test solution and calculated from the graph of chelating percentage against extract concentration. Lower IC50 value reflected better chelating activity. Essential oils of T. algeriensis collected from diverse locations during the vegetative and the flowering stages were subjected to screening for their chelating activities (Table 3). Our results showed that statistically significant differences of chelating activity were mainly observed when they were compared by the population location criteria. In fact, the variation of phenological stages did not have a statistically significant effect on the oil chelating activity for the most of populations (Table 3). Chelating activity was found to be very interesting (from 68 to 86 μg/ml) for populations 4 and 7 and which were comparable to the chemical EDTA (IC50 value = 40 μg/ml). Nevertheless, in the work of Bounatirou et al. [10], the antioxidant activity (DPPH assay) of essential oils obtained from the aerial parts of T. capitatus varied by the period of vegetation (vegetative, flowering or post-flowering) but no major differences were found between the antioxidant activity of the oils collected at different locations.

Conclusions

Analysis by GC/FID, GC/MS and 13C-NMR of Tunisian T. algeriensis essential oils allowed the identification of 71 compounds. The major compounds at the species level were α-pinene (7.41-13.94%), 1,8-cineole (7.55-22.07%), camphor (6.8-19.93%), 4-terpineol (1.55-11.86%), cis- sabinene hydrate (0.10-12.95%), terpenyl acetate (0-14.92%) and viridiflorol (0-11.49%). A high variation among populations was revealed for the majority of oil compounds. Nevertheless, neither the phenological stage, nor the interaction between population location and phenological stage were found to be statistically significant on the chemical composition of T. algeriensis essential oil. The spatial distribution of the populations was not concordant with the bioclimatic zones and seemed rather to be liked to local selective forces acting on the chemotype diversity. It is worthy to note that in the linear discriminant analysis, essentially minor compounds play an important role to distinguish between the chemical groups. Based on major compounds, the populations were represented by (α-pinene/1,8-cineole/cis-sabinene hydrate/camphor/viridiflorol), (1,8-cineole/camphor/terpenylacetate), (α-pinene/1,8-cineole/camphor), (1,8-cineole/camphor/4-terpineol) and (α-pinene/1,8-cineole/cis-sabinene hydrate/camphor/4-terpineol) chemotypes. The metal chelating activity of the essential oils was assessed and compared to synthetic EDTA. A variation of metal chelating activity of the oil was revealed according to population locations rather than to bioclimates or phenological stages. These results can be used to investigate the geographical location and the harvesting time of this plant for relevant industries.

Methods

Populations analyzed and sampling

The 8 populations of T. algeriensis collected from different bioclimatic and geographical zones and reported in Table 1 and Figure 1 were analyzed separately. A number of three to five individuals from each population were sampled over the entire population area at the vegetative (December 2009) and at the flowering (April 2010) stages. The distance between individuals exceeded 20 m, to avoid collection from close parents. The harvested samples size does not exceed 20 cm. After that, the fresh vegetable matter was first weighted and then dried on the shadow, until constancy of the weight (20 days). Separated from stems, aerials parts were subjected for essential oil extraction.

Essential oil extraction

The dry matter was submitted to hydrodistillation for 4 h, using a Clevenger-type apparatus. Each essential oil was dried over anhydrous sodium sulphate and stored in sealed vials protected from light at -20°C until analysis.

Essential oil analyses

Gas chromatography (GC)

A Hewlett-Packard 5890 series II gas chromatograph equipped with HP-5MS capillary column 30 m × 0.25 mm i.d., film thickness 0.25 μm; Hewlett-Packard) and connected to a flame ionization detector (FID) was used. The column temperature was programmed at 50°C for 1 min, then 7°C/min to 250°C, and then left at 250°C for 5 min. The injection port temperature was 240°C and that of the detector 250°C (split ratio: 1/60). The carrier gas was helium (99.995% purity) with a flow rate of 1.2 ml/min and the analysed sample volume was 2 μl. Percentages of the constituents were calculated by electronic integration of FID peak areas, without the use of response factor correction. Mean percentage of compounds in T. algeriensis essential oil represents the average calculated on three to five individuals. Retention indices (RI) were calculated for separate compounds relative to (C8-C25) n-alkanes mixture (Aldrich Library of Chemicals Standards) [13].

Gas chromatography/mass spectrometry (GC/MS)

The isolated volatile compounds were analysed by GC/MS, using a Hewlett-Packard 5890 series II gas chromatograph. The fused HP-5MS capillary column (the same as that used in the GC analysis) was coupled to a HP 5972A masse-selective detector (Hewlett-Packard, Palo Alto, CA, USA). The oven temperature was programmed at 50°C for 1 min, then 7°C/min to 250°C, and then left at 250°C for 5 min. The injection port temperature was 250°C and that of the detector was 280°C (split ratio: 1/100). The carrier gas was helium (99.995% purity) with a flow rate of 1.2 ml/min and the analysed sample volume was 2 μl. The mass spectrometer conditions were as follow: ionization voltage, 70 eV; ion source temperature, 150°C; electron ionization mass spectra were acquired over the mass range 50-550 m/z.

Volatile compounds identification

The essential oil compounds of T. algeriensis were identified by comparing the mass spectra data with spectra available from the Wiley 275 mass spectra libraries (software, D.03.00). Further identification confirmations were made referring to retention indices (RI) data generated from a series of known standards of n-alkanes mixture (C8-C25) [13] and to those previously reported in the literature [12, 1422].

13C-NMR analysis

NMR spectra were recorded on a Bruker AVANCE 400 Fourier Transform spectrometer operating at 100.13 MHz for 13C-NMR, equipped with a 5 mm probe, in CDCl3, with all shifts referred to internal TMS. 13C-NMR spectra of the oil samples were recorded with the following parameters: pulse width = 4 μs (flip angle 45°); acquisition time = 2.7 s for 128 K data table with a spectral width of 25 000 Hz (250 ppm); CPD mode decoupling; digital resolution = 0.183 Hz/pt. The number of accumulated scans was 2000-3000 for each sample depending of the available amount of oil (when available, 40 mg of oil in 0.5 ml of CDCl3). Identification of some compounds such as terpenyl acetate was assessed by the method developed and computerized in the laboratory of the team "Chimie et Biomasse", using home-made software, by comparison with spectral data of reference compounds compiled in a laboratory-built library [23].

Metal (Fe2+) chelating activity

Chelating activity of the essential oils was assessed by the ferrozine assay as described by Dinis et al. [24]. Ferrozine can quantitatively form complexes with Fe2+. In the presence of other chelating agents, the complex formation is disrupted with the result that the red color of the complex is decreased. Therefore, measurement of the rate of color reduction allows estimation of the chelating activity of the coexisting chelator. To 0.5 ml of essential oil solution prepared in methanol, 1.6 ml of deionised water and 0.05 ml of FeCl2 4H2O solution (2 mM) were added and left for incubation at room temperature for 5 min. Then, the reaction was initiated by adding 0.1 ml of ferrozine (5 mM), shaken vigorously and left standing at room temperature for 10 min. Absorbance of the solution was then measured at 562 nm. The chelating antioxidant activity for Fe2+ was calculated according to the following formula:

Metal chelating rate (% ) = Ac - As Ac ×100

where Ac is the absorbance of the control reaction and As is the absorbance of the tested sample. Essential oil concentration (μg/ml) corresponding to 50% ferrous iron chelating (IC50) was calculated from the graph plotting Fe2+ chelating activity against oil concentration. EDTA was used as a positive control and all determinations were carried out in triplicate.

Statistical analyses

The distribution of the 71 compounds identified from the essential oil was checked by a descriptive statistical analysis using the SPSS software for Windows™ (version 17, SPSS Inc., Chicago, IL, USA). The percentages of compounds were transformed using the arcsine transformation and all classic ones in order to improve the distribution property. However, these transformations did not yield satisfactory results. To assess the variation of the percentages of compounds (having Normal distribution) among the populations and/or phenological stages, a two-ways ANalysis Of VAriance (ANOVA) was performed using the SPSS software. For compounds having skewed distributions we have performed a non parametric two-way analysis of variance (Friedman test) implemented in the statistical programming environment R version 2.7.1 (R Development Core Team R: A Language and Environment for Statistical Computing, Vienna, Austria; 2008). We also computed Kruskal-Wallis test (SPSS) to compare stages and to test difference between populations. Statistical significance was assumed at p < 0.05. The chemical population structure and the relationship among populations were determined by Linear Discriminate analysis (LDA) performed on the percentages of all identified compounds for all populations using the SPSS program. The divergence between populations was also estimated by the Euclidean distances calculated among population pairs using the Numerical Taxonomy and Multivariate Analysis System (NTSys) software version 2.1 [25]. A dendrogram representing all populations and phenological stages was constructed using the Euclidian distance calculated among populations/phenological stages pairs and the Unweighted Pair Group Method with Arithmetic Averaging (UPGMA) algorithm [26]. The resulting tree was plotted using the treeview software version 1.6.6 [27]. Duncan's multiple range test (p < 0.05) and Student test (p < 0.05) were used to compare averages of essential oil yield and chelating activity among populations and phenological stages, respectively.

References

  1. Chaieb M, Boukhriss M: Flore suscinte et illustrée des zones arides et sahariennes de Tunisie. Association de la Protection de la Nature et de l'Environnement, l'Or du Temps (Ed.), Sfax, Tunisia. 1998, 42-

    Google Scholar 

  2. Le Floc'h E: Contribution à une étude ethnobotanique de la flore tunisienne. Ministère de l'Enseignement Supérieur et de la Recherche Scientifique (Ed.), Imprimerie Officielle de la République Tunisienne. 1983, 402-

    Google Scholar 

  3. Zouari N, Fakhfakh N, Zouari S, Bougatef A, Neffati M, Ayadi MA: Chemical composition, angiotensin I-converting enzyme inhibitory, antioxidant and antimicrobial activities of essential oil of Tunisian Thymus algeriensis Boiss. et Reut. (Lamiaceae). Food Bioprod Process. 2011, 89: 257-265. 10.1016/j.fbp.2010.11.006.

    Article  CAS  Google Scholar 

  4. Giordiani R, Hadef Y, Kaloustian J: Compositions and antifungal activities of essential oils of some Algerian aromatic plants. Fitoterapia. 2008, 79: 199-203. 10.1016/j.fitote.2007.11.004

    Article  Google Scholar 

  5. Hazzit M, Baaliouamer A, Verissimo AR, Faleiro ML, Miguel MG: Chemical composition and biological activities of Algerian Thymus oils. Food Chem. 2009, 116: 714-721. 10.1016/j.foodchem.2009.03.018.

    Article  CAS  Google Scholar 

  6. Ben El Hadj Ali I, Zaouali Y, Bejaoui A, Boussaid M: Variation of the chemical composition of essential oils in Tunisian populations of Thymus algeriensi Boiss. et Reut. (Lamiaceae) and implication for conservation. Chem Biodivers. 2010, 7: 1276-1289. 10.1002/cbdv.200900248

    Article  CAS  Google Scholar 

  7. Croteau R, Gershenzon J: Genetic control of monoterpene biosynthesis in mints. Genetic Engineering of Plant Secondary Metabolism. Recent Advances in Phytochemistry. Edited by: Ellis BE, Kuroki GW, Stafford HA. 1994, 183-New York: Plenum Press

    Google Scholar 

  8. Rahimmalek M, Tabatabaei BES, Etemadi N, Goli SAH, Arzani A, Zeinali H: Essential oil variation among and within six Achille species transferred from different ecological regions in Iran to the field conditions. Indust Crop Prod. 2009, 29: 348-355. 10.1016/j.indcrop.2008.07.001.

    Article  CAS  Google Scholar 

  9. Zouari N, Fakhfakh N, Zouari S, Sellami M, Abid M, Ayadi MA, Zaidi S, Neffati M: Volatile and lipid analyses by gas chromatography/mass spectrometry and nutraceutical potential of edible wild Malva aegyptiaca L. (Malvaceae). Int J Food Sci Nutr. 2011, 62: 600-608. 10.3109/09637486.2011.564157

    Article  CAS  PubMed  Google Scholar 

  10. Bounatirou S, Smit S, Miguel MG, Faleiro L, Rejeb MN, Neffati M, Costa MM, Figueiredo AC, Barroso JG, Pedro LG: Chemical composition, antioxidant and antibacterial activities of the essential oils isolated from Tunisian Thymus capitatus Hoff. et Link. Food Chem. 2007, 105: 146-155. 10.1016/j.foodchem.2007.03.059.

    Article  CAS  Google Scholar 

  11. Nejad Ebrahimi S, Hadian J, Mirjalili MH, Sonboli A, Yousefzadi M: Essential oil composition and antibacterial activity of Thymus caramanicu s at different phenological stages. Food Chem. 2008, 110: 927-931. 10.1016/j.foodchem.2008.02.083.

    Article  CAS  PubMed  Google Scholar 

  12. Liu L, Song G, Hu Y: GC/MS Analysis of the essential oils of Piper nigrum L. and Piper longum L. Chromatographia. 2007, 66: 785-790. 10.1365/s10337-007-0408-2.

    Article  CAS  Google Scholar 

  13. Kovàts E: Characterization of organic compounds by gas chromatography. Part 1. Retention indices of aliphatic halides, alcohols, aldehydes and ketones. Helv Chim Acta. 1958, 41: 1915-1932. 10.1002/hlca.19580410703.

    Article  Google Scholar 

  14. Adams RP: Identification of essential oil components by gas chromatography/quadrupole mass spectrometry. 2001, Carol Stream USA: Allured Publishing Corporation

    Google Scholar 

  15. Baranauskiene R, Venskutonis PR, Viskelis P, Dambrauskiene E: Influence of nitrogen fertilizers on the yield and composition of thyme (Thymus vulgaris). J Agric Food Chem. 2003, 51: 7751-7758. 10.1021/jf0303316

    Article  CAS  PubMed  Google Scholar 

  16. Asuming WA, Beauchamp PS, Descalzo JT, Dev BC, Dev V, Frost S, Ma CW: Essential oil composition of four Lomatium Raf. species and their chemotaxonomy. Biochem Syst Ecol. 2005, 33: 17-26. 10.1016/j.bse.2004.06.005.

    Article  CAS  Google Scholar 

  17. Hazzit M, Baaliouamer A, Faleiro ML, Miguel MG: Composition of the essential oils of Thymus and Origanum species from Algeria and their antioxidant and antimicrobial activities. J Agric Food Chem. 2006, 54: 6314-6321. 10.1021/jf0606104

    Article  CAS  PubMed  Google Scholar 

  18. Jalali-Heravi M, Zekavat B, Sereshti H: Characterization of essential oil components of Iranian geranium oil using gas chromatography-mass spectrometry combined with chemometric resolution techniques. J Chromatogr A. 2006, 1114: 154-163. 10.1016/j.chroma.2006.02.034

    Article  CAS  PubMed  Google Scholar 

  19. Kim MR, Abd El-Aty AM, Kim IS, Shim JH: Determination of volatile flavor components in danggui cultivars by solvent free injection and hydrodistillation followed by gas chromatographic-mass spectrometric analysis. J Chromatogr A. 2006, 1116: 259-264. 10.1016/j.chroma.2006.03.060

    Article  CAS  PubMed  Google Scholar 

  20. Ferhat MA, Meklati BY, Chemat F: Comparison of different isolation methods of essential oil from Citrus fruits: cold pressing, hydrodistillation and microwave "dry" distillation. Flavour Frag J. 2007, 22: 494-504. 10.1002/ffj.1829.

    Article  CAS  Google Scholar 

  21. Vagionas K, Graikou K, Ngassapa O, Runyoro D, Chinou I: Composition and antimicrobial activity of the essential oils of three Satureja species growing in Tanzania. Food Chem. 2007, 103: 319-324. 10.1016/j.foodchem.2006.07.051.

    Article  CAS  Google Scholar 

  22. Zouari S, Zouari N, Fakhfakh N, Bougatef A, Ayadi MA, Neffati M: Chemical composition and biological activities of a new essential oil chemotype of Tunisian Artemisia herba alb Asso. J Med Plants Res. 2010, 4: 871-880.

    CAS  Google Scholar 

  23. Gonny M, Bradesi P, Casanova J: Identification of the components of the essential oil from wild Corsican Daucus carota L. using 13C-NMR spectroscopy. Flavour Frag J. 2004, 19: 424-433. 10.1002/ffj.1330.

    Article  CAS  Google Scholar 

  24. Dinis TC, Maderia VM, Almeida LM: Action of phenolic derivatives (acetaminophen, salicylate, and 5-aminosalicylate) as inhibitors of membrane lipid peroxidation and as peroxyl radical scavengers. Arch Biochem Biophys. 1994, 315: 161-169. 10.1006/abbi.1994.1485

    Article  CAS  PubMed  Google Scholar 

  25. Rohlf FJ: Ntsys: Numerical taxonomy and multivariate analysis system, version 2.1. 2000, Setauket: Exeter software

    Google Scholar 

  26. Michener CD, Sokal RR: A quantitative approach to a problem of classification. Evolution. 1957, 11: 490-499.

    Article  Google Scholar 

  27. Page RD: Treeview: An application to display phylogenetic treeson personal computers. Comput Appl Biosci. 1996, 12: 357-358.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We are very grateful to Dr Fakher Frikha (ENIS) for his fruitful discussion and Dr Khaled Jebahi (Higher Institute of Applied Studies in Humanities, Medenine, Tunisia) for his kind help with English. Our thanks are due to Mr Béchir bouzbida (IRA) for his kind help with collecting plant material. We acknowledge Mr Slah Zaidi (IRA) for his technical assistance during essential oil analysis.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Nacim Zouari.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

SZ (chemist) carried out chemical composition analysis of essential oils. IA and AR (statisticians) contribute to the statistical analyses. NF and NZ (biochemists) realized the essential oil extraction, metal chelating activity and the redaction of the manuscript. Authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Zouari, N., Ayadi, I., Fakhfakh, N. et al. Variation of chemical composition of essential oils in wild populations of Thymus algeriensis Boiss. et Reut., a North African endemic Species. Lipids Health Dis 11, 28 (2012). https://doi.org/10.1186/1476-511X-11-28

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1476-511X-11-28

Keywords