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

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 [3][4][5][6]. 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 semiarid 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).
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.
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   Table 1. S2 S1 S2 S1 S2 S1 S2 S1 S2 S1 S2 S1 S2 S1 S2   The numbering refers to elution order of compounds from a HP-5MS column and their percentages were obtained by FID peak-area normalization. The percentage for each population represents the average calculated on n individuals (3 <n < 5). b RI, retention indices calculated against C 8 -C 25 n-alkanes mixture on the HP 5MS column. c For the detailed description of the populations (1-8) locations, see Table 1 and 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 13 C-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 semiarid 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 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.

trans-Piperitol
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 cissabinene 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  Table 1. 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].
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, Fe 2+ ion was known as the most important lipid oxidation prooxidant due to its high reactivity [12]. Chelating activity is presented by IC 50 value, defined as the concentration of the essential oil needed to chelate 50% of Fe 2+ present in the test solution and calculated from the graph of chelating percentage against extract concentration. Lower IC 50 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 (IC 50 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 13 C-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/cissabinene 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.

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 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 Populations Values represent mean ± standard deviation. Values followed by the same letter under the same row, are not significantly different (p > 0.05). For the same population, Student test was used to compare averages of essential oil yield (Y) and chelating activity (C) between phenological stages and is considered highly significant (**) at 0.001 ≤ P ≤ 0.01, significant (*) at 0.01 ≤ P ≤ 0.05 and not significant (ns) at P > 0.05.
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 (C 8 -C 25 ) 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 (C 8 -C 25 ) [13] and to those previously reported in the literature [12,[14][15][16][17][18][19][20][21][22]. 13 C-NMR analysis NMR spectra were recorded on a Bruker AVANCE 400 Fourier Transform spectrometer operating at 100. 13 MHz for 13 C-NMR, equipped with a 5 mm probe, in CDCl 3 , with all shifts referred to internal TMS. 13 [23].

Metal (Fe 2+ ) 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 Fe 2+ . 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 FeCl 2 4H 2 O 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 Fe 2+ 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 (IC 50 ) was calculated from the graph plotting Fe 2+ 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.