# Ecotypic Variations Affected the Biological Effectiveness ofThymus daenensis Celak Essential Oil

Authors: Fatemeh Elahian; Maryam Garshasbi; Zahra Mehri Asiabar; Neda Gholamian Dehkordi; Alireza Yazdinezhad; Seyed Abbas Mirzaei
Journal: Evidence-Based Complementary and Alternative Medicine
 (2021)

---

## Abstract

Thymus (Lamiaceae) is famous for its pharmacological properties. Thymus daenensis Celak (Avishan-e-denaee in Persian) is an endemic Thymus species in Iran and is traditionally used for its digestive, carminative, antitussive, antispasmodic, and expectorant attributes in folk medicine. Ecotypic oils were extracted and analyzed with the GC-MS. Their biological properties in terms of antimicrobial, antioxidant, and antigenotoxic activities were evaluated using the minimal inhibitory concentration, minimal bactericidal concentration, and DPPH, β-carotene, and comet assays. The GC-MS results for Thymus daenensis Celak oils revealed thymol (73.86%) and carvacrol (51.89%) as the most abundant components. Due to the results, reasonable bactericidal activity values range from 0.14 to 5.00 mg/ml, and fungicidal activity ranges from 0.17 to 0.58 mg/ml. The necessary oil free radical scavenging capacity (0.41–1.79 mg/ml), bleaching inhibitory activity (0.01–1.06 mg/ml), and genoprotective potential (1.04–7.78 mg/ml) indicated the dose-dependent activity. The results suggest that Thymus daenensis is an important antibacterial and antifungal bioresource. Additionally, the antioxidant and radical scavenging capacity suggests this species has a role as a natural preservative in oxidative diseases and in the prevention of food spoilage.


---

## Body

## 1. Introduction

Essential oils are complex natural and aromatic compounds, which are derived from different parts of the plant using extraction, fermentation, expression, and steam distillation. Terpenoids, phenolics, polyphenolics, polyacetylenes, lectins, polypeptides, and alkaloids are the main plant essential oil (EO) constituents, which can exhibit potent antibacterial, antiviral, antifungal, antioxidative, and cytotoxic properties. The wide spectrum of medicinal properties of herbal EOs medicinal make them strong candidates for use in the pharmaceutical, food, and cosmetic industries [1, 2].Amongst the 928 species ofThymus or “Avishan,” which are widely used for medicinal purposes, Thymus daenensis or “Avishan-e-denaee” is one of the most important Iranian endemic species. The genus Thymus belongs to the family Lamiaceae (Labiatae), one of the best-known flowering plant families with approximately 236 genera [3]. Thymus daenensis comprises two known subspecies: 1. T. daenensis Celak subsp. daenensis and 2. T. daenensis subsp. lancifolius. Thymus daenensis is a valuable spice plant that occurs in the Zagros mountain chain and grows as a perennial dwarf shrub (6–30 cm high) with flowers and lanceolate leaves. Thymus phenolic and nonphenolic derivatives are mainly used to treat gastrointestinal disorders, nervousness, and headaches in the form of tonic or herbal tea. Thymus is an aromatic plant that contains monoterpenoid essential oils in its leaves and floral parts, similar to most species of the Lamiaceae. Previous research has also focused on the EO-mediated biomedical activities of other species of Thymus including antimicrobial, antioxidant, and also to less extent their antigenotoxic effects [4, 5].The main concern about the medicinal plants is related to the environmental conditions effect on their medicinal activities diversity. Several studies have indicated that the climate, temperature, light, height above sea level, latitude, and longitude are the most important factors affecting the physical and biological behaviors of plants including their genetics, immunology, physiology, morphology, population divergence, and evolution properties [6]. Earlier studies highlighted the climate fluctuations influence on the EOs’ composition of ecotypic Thymus kotschanus and as well examined biological variety quantitatively in their antimicrobial, antioxidative, and antigenotoxic behaviors [7]. At the first step of this study, we focused on the T. daenensis EO-mediated antimicrobial, antioxidative, and antigenotoxic behaviors. Then, climate and altitude effects on the EOs’ composition of ecotypic T. daenensis and their biological diversity were investigated. In addition to the antimicrobial activity, the antioxidant and antigenotoxic activities were reported.


## 2. Materials and Methods

### 2.1. Reagents and Media

Butylated hydroxytoluene (BHT), 2,2-diphenyl-1-picrylhydrazyl (DPPH), polyunsaturated linoleic acid, and beta-carotene (β-Carotene) were provided from Sigma-Aldrich (Sigma-Aldrich, Deisenhofen, Germany). Low-melting agarose and Ficoll-Paque (density 1.077 g/ml) were attained from Thermo Fisher Scientific (Rockford, IL, US) and GE Healthcare (Madison, WI, USA), respectively. Luria-Bertani, Mueller-Hinton, and Sabouraud dextrose Media were purchased from Merck Company (Darmstadt, Germany). Also, other pure analytical grade compounds (>99% purity) were attained from the Iranian commercial resources.


### 2.2. Standard Bacterial and Fungal Microorganisms

Three pure lyophilized Gram-negative bacteria ofPseudomonas aeruginosa (ATCC 27853), Salmonella typhimurium (ATCC 13311), and Escherichia coli (ATCC 8739); three Gram-positive bacteria of Enterococcus faecalis (ATCC 29212), Bacillus subtilis (ATCC 12711), and Staphylococcus aurous (ATCC 25923), and also two fungi Candida albicans (ATCC 10231), and Aspergillus niger (ATCC 9142) were purchased from the “Iranian Research Organization for Science and Technology” (a member of the WFCC). Methicillin-resistant Staphylococcus aurous (MRSA) was a hospital-resistant sample and was devotedly provided by Professor Moghimi. Our department confirmed the purity and character.


### 2.3. Medicinal Plant Materials and Hydrodistillation

Fresh aerial parts of ecotypicThymus daenensis Celak were collected from their habitation in Ghazvin, Isfahan, Markazi, Lorestan, and Zanjan provinces. After that was the identification of species; they were harvested from the proposed areas randomly, and the shades dried for two weeks at ambient temperature (20–25°C). As a part of an ecotypic morphology study, the fresh plant ecotypes macroscopic diameters were measured, and also images were recorded with a digital camera (COOLPIX P510, Nikon, Japan) equipped with GPS. Voucher specimens were recognized and deposited at the Herbarium of Department of Pharmacognosy, Faculty of Pharmacy, Zanjan University of Medical Sciences, Zanjan, Iran. The plants were subjected to hydrodistillation using a modified Clevenger’s apparatus at 100°C for 3 hours. During the hydrodistillation process, plant materials were placed in boiling water. Essential oils could be easily separated from the collecting tube water phase. Anhydrous sodium sulfate was used as drying agent. The oils were stored in the opaque air-tight sealed bottles at 4°C for 3 months, until the time of physicochemical and biological tests running [8].


### 2.4. Physical Characterization of the Essential Oils

The yield was defined as the amount of essential oils formed (ml) related to the 100 g of consumed plant. The oil refractive indices and densities were measured using an Atago RX-7000a refractometer, and a liquid micro-pycnometer. The oils different serial concentrations were provided in methanol (1 : 2, 1 : 4, 1 : 8, 1 : 16, 1 : 32, 1 : 64, 1 : 128, 1 : 256 mg/ml), and absorption spectrum was recorded in absorbance range (230–1000 nm) using spectrophotometer (TECAN, Infinite®-M200; Austria). The oils molar absorptivity potencies against the UV light were attained using Beer’s law (A = εbC, where A, ε, b, and C stand for absorbance, molar absorptivity, path length, and oil concentration, respectively) [7].


### 2.5. Oils Characterization by Using Gas Chromatography (GC-FID)

An HP-6890 series gas chromatography system was used for theT. daenensis oil analyses. This Agilent instrument was equipped with an HP-5 capillary column packed with 5% phenyl, and also 95% methylpolysiloxane (30 m × 0.32 mm ×  0.25 μm capillary dimensions), and a flame ionization detector (FID). Nitrogen was applied as the carrier, and this gas flow was set on 1 ml/min. The column oven temperature was accurately set for 3 minutes at 60°C, and after that progressively increased to 250°C for 65 minutes. The injection port and the flame ionization detector were heated to 255 and 290°C, respectively. The sample split ratio was 1 : 10 [9].


### 2.6. Identification of the Oils Constituents Using GC-MS

Volatile constituents were identified using a gas chromatograph (Agilent, GC-7890A, USA) connected to an Electron Ionization Mass Spectroscopy, EIMS-5975C, detector (70 eV energy), and a capillary HP-5 MS well-packed column containing 95% of the methylpolysiloxane, and 5% of phenyl (dimensions: 30 m × 0.25 mm × 0.25μm). The central oven was preset for 5 minutes at 40°C. After that, the temperature was progressively increased to 230°C for 20 minutes, and at the final step, temperature was raised up to 280°C for 2 minutes (30°C/min ramping rate). Detector and injector temperatures were fixed at 240 and 250°C, respectively. One μl of the oils was injected with a split ratio of 1 : 10 to the instrument, and also the helium flow rate was tuned to 1 ml/min. The GC-MS analysis was completed using various internal standards for components quantification. The oil components arrangement was conducted based on the retention index ratios along with the homologous normal alkane series retention times under the same running condition. Additional parameters, which were used for the composition description, were the Wiley7n.L library, NIST library, published data fragmentation model, and the authentic compounds mass spectroscopy. The results were expressed as percentages of the individual constituent whole peak area in association with the whole region of all standardized peaks in the chromatogram [10, 11].


### 2.7. Antibacterial Activities

The study of essential oils antibacterial activity was conducted using NCCLS (National Committee for Clinical Laboratory Standards) guidelines, and the broth microdilution method was applied for the reliable and reproducible results. Cultivated bacterial strains were subcultured on Mueller Hinton Broth (MHB) to an initial OD600 = 0.1, and incubated again for a few hours at 37°C to an OD600 = 0.5, and were used as the inocula after 1 : 200 dilutions in MHB (standardized at approximately 1 × 106 CFU/ml). 100 μl of the bacterial suspension was dispensed in a 96-well plate consisting of 100 μl of serially diluted oils in MHB during the range of 0.01–5.12 mg/ml. Up and about 2.0% dimethyl sulfoxide (DMSO) was used as a co-solvent, and it does not alter the test organisms’ growth. Negative controls included a standard bacterial suspension, and the DMSO maximum doses without any test components. A serial dilution (a range of 0.5–256 μg/ml) containing pure ampicillin sodium salt was used as the standard antibacterial agent. Plates were agitated and incubated for 24 hours at 37°C. Subsequently, 20 μl of MTT solution (0.2 mg/ml) was applied to each well and incubated for 30 minutes at 37°C. The lowest sample concentrations without any color changes were determined as the MIC (minimum inhibitory concentration). A 5 μl portion from the colorless wells was subcultured on Mueller Hinton Agar (MHA) and was incubated for 24 hours at 37°C. MBC (minimum bactericidal concentration) was considered as the least oil concentration that killed ≥99.5% of the inoculated bacteria [12, 13].


### 2.8. Antifungal Activities

All of the essential oils ecotypes were examined for their fungal toxic activities individually, by an adapted broth microdilution assay. The Aspergillus was cultured in Sabouraud dextrose agar at 30°C for a week and also stored at 4°C for another week in order to spore germinating. The fungal spores were suspended in sterile normal saline solution supplemented with 0.1% v/v Tween-80, and isolated using several sterile cheesecloth aseptically layers. That flow through the spore suspension was counted using a hemocytometer and adjusted at 106 spores/ml. 100 μl of the conidial spore suspensions was added to 100 μl of the essential oil serial dilutions from 0.01 to 5.12 mg/ml in Sabouraud dextrose broth. The plates were shaken for 30 seconds and after that incubated at 28 ± 0.5°C and in 80% relative humidity for 72 hours. Appropriate negative controls containing the standard fungal spore suspensions, and DMSO without the test components were produced, in order to confirm the sporangiospores viability, and any unwanted intrinsic effects absence. A fresh Sabouraud dextrose culture was diluted and adjusted to 5 × 106 CFU/ml for Candida albicans, and 100 μl of the suspension was dispensed in a 96-well plate containing 100 μl of serially diluted oils and was cultivated for 48 hours at 28°C. The pure amphotericin-B serial dilution (0.1–12.8 μg/ml) was used as the reference compound. MIC (minimum inhibitory concentration) was the least oil concentration and represented no visible fungal mycelia next to the incubation period. A 5 μl portion was subcultivated on Sabouraud dextrose agar and was from the wells with complete fungal growth inhibition and incubated for 72 hours at 28°C. MFC (minimum fungicidal concentration) was considered as the lowest oil concentration that killed ≥99.5% of the inoculated fungi [14].


### 2.9. DPPH Method and Antioxidant Activities

The oil samples radical scavenging activities were evaluated against the discoloration of DPPH free radicals purple solution. In order to achieve this goal, 2 ml of prepared methanolic solutions of DPPH (80μg/ml) was mixed with the oils’ methanolic dilutions equivalent volumes (final concentration: 0.01–1.28 mg/ml). The mixtures were incubated for 30 minutes in the darkness at ambient temperature. Absorbance values were recorded at 517 nm using a Tecan Infinite M200 spectrophotometer, and the oils free radical scavenging capacities (SC) were calculated according to the following equation:(1)SC%=ANC−ASampleANC×100.Negative controls (NCs) were composed of all the reagents excluding the antioxidant oils. In addition, reference antioxidants were examined like butylated hydroxytoluene and vitamin C (0.5 to 128μg/ml) in order to compare the compounds scavenging potencies. The relative scavenging activities were plotted against the essential oil concentrations, and the SC50 was considered as the essential oil concentration in order to result in 50% inhibition [15].


### 2.10.β-Carotene/Linoleic Acid Method and Antioxidant Activities

Thymus antioxidant potency was measured using the β-carotene bleaching inhibition in the linoleic acid presence. β-carotene/linoleic acid emulsion was prepared by mixing 1 mg β- carotene, 25 μl of linoleic acid, and 200 mg Tween-20 in 1 ml of chloroform. The solvent was vacuum-evaporated using a rotary evaporator for 30 minutes at 40°C. After that, 100 ml of oxygenated water was composed with robust shaking, in order to form a stable emulsion. 350 μl of the essential oil dilutions in ethanol (ranging 0.01–1.28 mg/ml) was mixed with 2.5 ml of the β-carotene/linoleic emulsion gently, at the final stage. All tubes were placed in a shaker incubator (50 rpm) for 2 hours, in the light and at 50°C. The absorbance difference was recorded for 2 hours, using a spectrophotometer against a negative control (NC), which contains all the reagents except oils at 470 nm. The beta-carotene oxidation inhibition percentage (BIC %) was calculated as follows:(2)BIC%=1−ASample0min−ASample120minANC0min−ANC120min×100.Appropriate routine reference substances amounts were also applied including vitamin C and butylated hydroxytoluene (0.5 to 128μg/ml) for the antioxidant capacity evaluation. The relative antioxidative activities plot versus the oil concentrations was drawn, and BIC50 was measured from the best fitted line equation in the scatter plot [16].


### 2.11. Antigenotoxicity Activity (Comet Assay)

Comet assay is a highly sensitive method for studying of the compounds antigenotoxic activity against the DNA damage in lymphocytes. Lymphocytes were isolated from the volunteer blood samples using Ficoll-Paque TM PLUS2, and a cell suspension of 2 × 105 cells/ml was prepared in the phosphate buffered saline. Then, 1 × 104 cells were treated with the sampled oils serial concentration (0.01–5 mg/ml) for 30 minutes at 4°C. The cells were centrifuged and resuspended in 50 μl of PBS and were finally mixed with the equal preheated low-melting-point agarose volume (1.5% w/v). By passing 30 minutes from incubation in ice, the cells were coated again with 100 μl of low-melting agarose (0.75% w/v) in order to form a flat sandwich. After that, cell lysis was performed using ice-cold lysis buffer for 4 hours (90 mM EDTA, 2.25 M NaCl, 0.7% w/v NaOH, 9 mM Tris base, 1% v/v Triton X-100, 10% v/v DMSO; pH = 10). Then, the slides were neutralized and electrophoresed on ice under the alkaline conditions (1 mM EDTA, 300 mM NaOH; pH ≥ 13) at 0.30 A for 45 minutes. At the next stage, the microgels were neutralized, stained, and pictured with neutral Tris-buffer, ethidium bromide, and fluorescent microscope, respectively. DNA fragmentation was scored using Tritek CometScore13TM software, and ([tail DNA/(head DNA + tail DNA)]×100) equation was applied for the tail DNA% determining. CIC50 is the sample concentration, which protects 50% genomic DNA from fragmentation using the scatter plot equation [7].


### 2.12. Statistical Analyses

All procedures were repeated three times independently, and their results were reported as the mean values ± the mean standard error. Quantitative data was statistically interpreted using IBM SPSS-22 software. Differences between experimental groups were determined using one-way ANOVA with Tukey’s post hoc test.P values less than 0.05 were considered as statistically significant.


## 2.1. Reagents and Media

Butylated hydroxytoluene (BHT), 2,2-diphenyl-1-picrylhydrazyl (DPPH), polyunsaturated linoleic acid, and beta-carotene (β-Carotene) were provided from Sigma-Aldrich (Sigma-Aldrich, Deisenhofen, Germany). Low-melting agarose and Ficoll-Paque (density 1.077 g/ml) were attained from Thermo Fisher Scientific (Rockford, IL, US) and GE Healthcare (Madison, WI, USA), respectively. Luria-Bertani, Mueller-Hinton, and Sabouraud dextrose Media were purchased from Merck Company (Darmstadt, Germany). Also, other pure analytical grade compounds (>99% purity) were attained from the Iranian commercial resources.


## 2.2. Standard Bacterial and Fungal Microorganisms

Three pure lyophilized Gram-negative bacteria ofPseudomonas aeruginosa (ATCC 27853), Salmonella typhimurium (ATCC 13311), and Escherichia coli (ATCC 8739); three Gram-positive bacteria of Enterococcus faecalis (ATCC 29212), Bacillus subtilis (ATCC 12711), and Staphylococcus aurous (ATCC 25923), and also two fungi Candida albicans (ATCC 10231), and Aspergillus niger (ATCC 9142) were purchased from the “Iranian Research Organization for Science and Technology” (a member of the WFCC). Methicillin-resistant Staphylococcus aurous (MRSA) was a hospital-resistant sample and was devotedly provided by Professor Moghimi. Our department confirmed the purity and character.


## 2.3. Medicinal Plant Materials and Hydrodistillation

Fresh aerial parts of ecotypicThymus daenensis Celak were collected from their habitation in Ghazvin, Isfahan, Markazi, Lorestan, and Zanjan provinces. After that was the identification of species; they were harvested from the proposed areas randomly, and the shades dried for two weeks at ambient temperature (20–25°C). As a part of an ecotypic morphology study, the fresh plant ecotypes macroscopic diameters were measured, and also images were recorded with a digital camera (COOLPIX P510, Nikon, Japan) equipped with GPS. Voucher specimens were recognized and deposited at the Herbarium of Department of Pharmacognosy, Faculty of Pharmacy, Zanjan University of Medical Sciences, Zanjan, Iran. The plants were subjected to hydrodistillation using a modified Clevenger’s apparatus at 100°C for 3 hours. During the hydrodistillation process, plant materials were placed in boiling water. Essential oils could be easily separated from the collecting tube water phase. Anhydrous sodium sulfate was used as drying agent. The oils were stored in the opaque air-tight sealed bottles at 4°C for 3 months, until the time of physicochemical and biological tests running [8].


## 2.4. Physical Characterization of the Essential Oils

The yield was defined as the amount of essential oils formed (ml) related to the 100 g of consumed plant. The oil refractive indices and densities were measured using an Atago RX-7000a refractometer, and a liquid micro-pycnometer. The oils different serial concentrations were provided in methanol (1 : 2, 1 : 4, 1 : 8, 1 : 16, 1 : 32, 1 : 64, 1 : 128, 1 : 256 mg/ml), and absorption spectrum was recorded in absorbance range (230–1000 nm) using spectrophotometer (TECAN, Infinite®-M200; Austria). The oils molar absorptivity potencies against the UV light were attained using Beer’s law (A = εbC, where A, ε, b, and C stand for absorbance, molar absorptivity, path length, and oil concentration, respectively) [7].


## 2.5. Oils Characterization by Using Gas Chromatography (GC-FID)

An HP-6890 series gas chromatography system was used for theT. daenensis oil analyses. This Agilent instrument was equipped with an HP-5 capillary column packed with 5% phenyl, and also 95% methylpolysiloxane (30 m × 0.32 mm ×  0.25 μm capillary dimensions), and a flame ionization detector (FID). Nitrogen was applied as the carrier, and this gas flow was set on 1 ml/min. The column oven temperature was accurately set for 3 minutes at 60°C, and after that progressively increased to 250°C for 65 minutes. The injection port and the flame ionization detector were heated to 255 and 290°C, respectively. The sample split ratio was 1 : 10 [9].


## 2.6. Identification of the Oils Constituents Using GC-MS

Volatile constituents were identified using a gas chromatograph (Agilent, GC-7890A, USA) connected to an Electron Ionization Mass Spectroscopy, EIMS-5975C, detector (70 eV energy), and a capillary HP-5 MS well-packed column containing 95% of the methylpolysiloxane, and 5% of phenyl (dimensions: 30 m × 0.25 mm × 0.25μm). The central oven was preset for 5 minutes at 40°C. After that, the temperature was progressively increased to 230°C for 20 minutes, and at the final step, temperature was raised up to 280°C for 2 minutes (30°C/min ramping rate). Detector and injector temperatures were fixed at 240 and 250°C, respectively. One μl of the oils was injected with a split ratio of 1 : 10 to the instrument, and also the helium flow rate was tuned to 1 ml/min. The GC-MS analysis was completed using various internal standards for components quantification. The oil components arrangement was conducted based on the retention index ratios along with the homologous normal alkane series retention times under the same running condition. Additional parameters, which were used for the composition description, were the Wiley7n.L library, NIST library, published data fragmentation model, and the authentic compounds mass spectroscopy. The results were expressed as percentages of the individual constituent whole peak area in association with the whole region of all standardized peaks in the chromatogram [10, 11].


## 2.7. Antibacterial Activities

The study of essential oils antibacterial activity was conducted using NCCLS (National Committee for Clinical Laboratory Standards) guidelines, and the broth microdilution method was applied for the reliable and reproducible results. Cultivated bacterial strains were subcultured on Mueller Hinton Broth (MHB) to an initial OD600 = 0.1, and incubated again for a few hours at 37°C to an OD600 = 0.5, and were used as the inocula after 1 : 200 dilutions in MHB (standardized at approximately 1 × 106 CFU/ml). 100 μl of the bacterial suspension was dispensed in a 96-well plate consisting of 100 μl of serially diluted oils in MHB during the range of 0.01–5.12 mg/ml. Up and about 2.0% dimethyl sulfoxide (DMSO) was used as a co-solvent, and it does not alter the test organisms’ growth. Negative controls included a standard bacterial suspension, and the DMSO maximum doses without any test components. A serial dilution (a range of 0.5–256 μg/ml) containing pure ampicillin sodium salt was used as the standard antibacterial agent. Plates were agitated and incubated for 24 hours at 37°C. Subsequently, 20 μl of MTT solution (0.2 mg/ml) was applied to each well and incubated for 30 minutes at 37°C. The lowest sample concentrations without any color changes were determined as the MIC (minimum inhibitory concentration). A 5 μl portion from the colorless wells was subcultured on Mueller Hinton Agar (MHA) and was incubated for 24 hours at 37°C. MBC (minimum bactericidal concentration) was considered as the least oil concentration that killed ≥99.5% of the inoculated bacteria [12, 13].


## 2.8. Antifungal Activities

All of the essential oils ecotypes were examined for their fungal toxic activities individually, by an adapted broth microdilution assay. The Aspergillus was cultured in Sabouraud dextrose agar at 30°C for a week and also stored at 4°C for another week in order to spore germinating. The fungal spores were suspended in sterile normal saline solution supplemented with 0.1% v/v Tween-80, and isolated using several sterile cheesecloth aseptically layers. That flow through the spore suspension was counted using a hemocytometer and adjusted at 106 spores/ml. 100 μl of the conidial spore suspensions was added to 100 μl of the essential oil serial dilutions from 0.01 to 5.12 mg/ml in Sabouraud dextrose broth. The plates were shaken for 30 seconds and after that incubated at 28 ± 0.5°C and in 80% relative humidity for 72 hours. Appropriate negative controls containing the standard fungal spore suspensions, and DMSO without the test components were produced, in order to confirm the sporangiospores viability, and any unwanted intrinsic effects absence. A fresh Sabouraud dextrose culture was diluted and adjusted to 5 × 106 CFU/ml for Candida albicans, and 100 μl of the suspension was dispensed in a 96-well plate containing 100 μl of serially diluted oils and was cultivated for 48 hours at 28°C. The pure amphotericin-B serial dilution (0.1–12.8 μg/ml) was used as the reference compound. MIC (minimum inhibitory concentration) was the least oil concentration and represented no visible fungal mycelia next to the incubation period. A 5 μl portion was subcultivated on Sabouraud dextrose agar and was from the wells with complete fungal growth inhibition and incubated for 72 hours at 28°C. MFC (minimum fungicidal concentration) was considered as the lowest oil concentration that killed ≥99.5% of the inoculated fungi [14].


## 2.9. DPPH Method and Antioxidant Activities

The oil samples radical scavenging activities were evaluated against the discoloration of DPPH free radicals purple solution. In order to achieve this goal, 2 ml of prepared methanolic solutions of DPPH (80μg/ml) was mixed with the oils’ methanolic dilutions equivalent volumes (final concentration: 0.01–1.28 mg/ml). The mixtures were incubated for 30 minutes in the darkness at ambient temperature. Absorbance values were recorded at 517 nm using a Tecan Infinite M200 spectrophotometer, and the oils free radical scavenging capacities (SC) were calculated according to the following equation:(1)SC%=ANC−ASampleANC×100.Negative controls (NCs) were composed of all the reagents excluding the antioxidant oils. In addition, reference antioxidants were examined like butylated hydroxytoluene and vitamin C (0.5 to 128μg/ml) in order to compare the compounds scavenging potencies. The relative scavenging activities were plotted against the essential oil concentrations, and the SC50 was considered as the essential oil concentration in order to result in 50% inhibition [15].


## 2.10.β-Carotene/Linoleic Acid Method and Antioxidant Activities

Thymus antioxidant potency was measured using the β-carotene bleaching inhibition in the linoleic acid presence. β-carotene/linoleic acid emulsion was prepared by mixing 1 mg β- carotene, 25 μl of linoleic acid, and 200 mg Tween-20 in 1 ml of chloroform. The solvent was vacuum-evaporated using a rotary evaporator for 30 minutes at 40°C. After that, 100 ml of oxygenated water was composed with robust shaking, in order to form a stable emulsion. 350 μl of the essential oil dilutions in ethanol (ranging 0.01–1.28 mg/ml) was mixed with 2.5 ml of the β-carotene/linoleic emulsion gently, at the final stage. All tubes were placed in a shaker incubator (50 rpm) for 2 hours, in the light and at 50°C. The absorbance difference was recorded for 2 hours, using a spectrophotometer against a negative control (NC), which contains all the reagents except oils at 470 nm. The beta-carotene oxidation inhibition percentage (BIC %) was calculated as follows:(2)BIC%=1−ASample0min−ASample120minANC0min−ANC120min×100.Appropriate routine reference substances amounts were also applied including vitamin C and butylated hydroxytoluene (0.5 to 128μg/ml) for the antioxidant capacity evaluation. The relative antioxidative activities plot versus the oil concentrations was drawn, and BIC50 was measured from the best fitted line equation in the scatter plot [16].


## 2.11. Antigenotoxicity Activity (Comet Assay)

Comet assay is a highly sensitive method for studying of the compounds antigenotoxic activity against the DNA damage in lymphocytes. Lymphocytes were isolated from the volunteer blood samples using Ficoll-Paque TM PLUS2, and a cell suspension of 2 × 105 cells/ml was prepared in the phosphate buffered saline. Then, 1 × 104 cells were treated with the sampled oils serial concentration (0.01–5 mg/ml) for 30 minutes at 4°C. The cells were centrifuged and resuspended in 50 μl of PBS and were finally mixed with the equal preheated low-melting-point agarose volume (1.5% w/v). By passing 30 minutes from incubation in ice, the cells were coated again with 100 μl of low-melting agarose (0.75% w/v) in order to form a flat sandwich. After that, cell lysis was performed using ice-cold lysis buffer for 4 hours (90 mM EDTA, 2.25 M NaCl, 0.7% w/v NaOH, 9 mM Tris base, 1% v/v Triton X-100, 10% v/v DMSO; pH = 10). Then, the slides were neutralized and electrophoresed on ice under the alkaline conditions (1 mM EDTA, 300 mM NaOH; pH ≥ 13) at 0.30 A for 45 minutes. At the next stage, the microgels were neutralized, stained, and pictured with neutral Tris-buffer, ethidium bromide, and fluorescent microscope, respectively. DNA fragmentation was scored using Tritek CometScore13TM software, and ([tail DNA/(head DNA + tail DNA)]×100) equation was applied for the tail DNA% determining. CIC50 is the sample concentration, which protects 50% genomic DNA from fragmentation using the scatter plot equation [7].


## 2.12. Statistical Analyses

All procedures were repeated three times independently, and their results were reported as the mean values ± the mean standard error. Quantitative data was statistically interpreted using IBM SPSS-22 software. Differences between experimental groups were determined using one-way ANOVA with Tukey’s post hoc test.P values less than 0.05 were considered as statistically significant.


## 3. Results

### 3.1. Geographical Effects on the Physicochemical Properties of the Ecotypes

Different morphological analysesof Thymus daenensis ecotypes indicated a microenvironmental conditions distinct effect including altitude, humidity, and temperature on the plant height and also the total dry matter. This flowering plant is growing up to 35 cm tall with short trichomes on its leaves. The isolation process yields the yellowish oil with different productivity (0.36–0.86% v/w), based on the Thymus geographical origins. In addition, other ecotypic oils physical properties presented diversity including density and molar extinction coefficient at 280 nm, and it seems that they can also be directly affected by environmental conditions (Table 1).Table 1
Environmental conditions and physical properties ofThymus daenensis essential oils.
EcotypeEcotype locationCity, provinceAltitude (m)Temperature (°C)Average humidity (%)Density (g/ml) at 25°Cε280 nm (l/(cm × g))Refractive index at 25°CYield (v/w %)MaxMinTD134 59 17NGharechaie highlands, Marekazi190023.96400.94 ± 0.122.83 ± 1.011.32 ± 0.010.82 ± 0.2149 19 49ETD233 15 00NKhoramabad, Lorestan183038.8−1.6570.84 ± 0.041.98 ± 0.771.45 ± 0.080.36 ± 0.2848 30 00ETD333 25 00NKhoramabad, Lorestan190038.8−1.6570.98 ± 0.057.41 ± 2.321.38 ± 0.030.56 ± 0.1348 40 00ETD439 52 42NDaran, Isfahan250023−3.1420.96 ± 0.113.22 ± 1.741.39 ± 0.060.60 ± 0.2350 23 21ETD534 05 06NKharkan, Markazi196535.5−4.4310.88 ± 0.045.38 ± 0.821.44 ± 0.030.62 ± 0.1749 24 50ETD634 11 38NArak, Markazi2404273460.93 ± 0.062.92± 0.191.46 ± 0.090.86 ± 0.2549 29 41ETD732 55 52NFereidonshahr, Isfahan261521.3−4.5500.88 ± 0.027.99 ± 1.581.51 ± 0.100.59 ± 0.2350 07 01ETD833 46 24NRobatkhomein, Markazi240725.23510.89 ± 0.0811.43 ± 1.631.47 ± 0.040.69 ± 0.1249 52 57ETD936 33 15NZanjan, Zanjan195032.1−7.5540.90 ± 0.076.10 ± 1.941.51 ± 0.030.74 ± 0.2648 22 14ETD1036 30 00NChore, Ghazvin1950182510.95 ± 0.032.66 ± 0.191.48 ± 0.060.75 ± 0.1149 43 00E


### 3.2. Essential Oil Constituents in theThymus Ecotypes

Sixty-one compounds were identified by the GC/MS analysis, in differentT. daenensis ecotypes, and were representing 95 up to more than 99% of the total oil compositions. As displayed in Table 2, thymol (0.31–73.86%), carvacrol (5.14–51.89%), linalool (0.3–22.49%), pulegone (up to 17.57%), geraniol (2.02–14.84%), α-terpineol (0.58–14.01%), p-cymene (0.12–11.19%), o-cymene (5.19–8.05%), geranyl acetate (0.16–6.11%), neomenthol (up to 5.69%), β- caryophyllene (1.81–4.97%), terpineol (0.59–4.68%), and borneol (0.01–4.04%) are the major constituents in the T. daenensis volatile oil. Some compounds were only detected in special ecotypes, and they are identified as δ-3-carene, δ-elemene, α-farnesene, α-himachalene, isoborneol, isopulegol, Z-jasmone, Z-linalool oxide, 1-menthol, neomenthol, menthone, α- muurolene, piperitone, pulegone, and thymol acetate.Table 2
Essential oil compositions for ecotypicThymus daenensis (%).
CompoundsRIRI∗ [11]TD1TD2TD3TD4TD5TD6TD7TD8TD9TD10α-Pinene9389390.860.420.260.830.10Camphene9539540.140.080.36β-Pinene9789790.580.260.291-Octen-3-ol9809790.080.230.18β-Myrcene9909910.810.70α-Terpinene101710170.060.462.800.75p-Cymene102510260.122.565.126.5311.190.239.261.640.38O-Cymene102610275.198.05δ-3-Carene102910310.121,8-Cineole103010320.101.762.761.080.400.920.181.730.130.21Z-β-Ocimene103710380.20E-β-Ocimene104710500.210.18γ-Terpinene105910600.030.441.780.661.401.683.910.26Z-Linalool oxide107410750.34α-Terpinolene108710890.130.41Linalool109810971.5718.241.481.021.600.4422.490.322.010.30Z-β-Terpineol114411440.160.600.240.590.12Camphor114811502.880.260.480.62Menthone115211532.22Isopulegol115011522.12Isoborneol116111620.62Neomenthol116711685.69Borneol116911702.851.992.443.991.441.644.040.690.01Terpinen-4-ol117511771.040.604.680.700.591.100.90Isomenthol118111833.00α-Terpineol118911901.920.581.620.582.2814.010.61Thymol methyl ether123012350.220.12Pulegone1234123717.57Carvacrol methyl ether124312453.141.90Piperitone124912531.97Geraniol125312542.0214.84Thymol1303130460.610.3158.9464.5270.9560.600.7339.4525.0073.86Carvacrol130813107.4520.9014.865.148.826.2251.8916.2926.127.89δ-Elemene133713390.08Piperitenone134013431.90α-Terpinyl acetate134713491.100.080.41Thymol acetate135413520.15Eugenol135813590.390.05α-Copaene137613770.24Geranyl acetate138113811.850.166.110.94β-Bourbonene138913880.140.440.040.040.120.120.14β-Elemene139013930.08Z-Jasmone139313950.14Β-Caryophyllene142014194.563.733.084.601.812.323.504.973.403.60Aromadendrene144014410.160.200.53α-Himachalene145114520.35α-Humulene145614550.240.120.170.100.110.241.150.20γ-Muurolene148014820.310.220.290.250.390.78Germacrene D148114850.210.750.20Viridiflorene149214970.460.260.470.620.51α-Muurolene150115000.28α -Farnesene150415060.34β-Bisabolene150715061.470.400.820.672.650.830.660.851.24Cis-α-Bisabolene150815070.660.260.233.820.111.680.201.332.684.61γ-Cadinene151515140.410.330.260.160.150.150.28δ-Cadinene152415230.410.280.520.310.440.590.66Spathulenol157915780.380.650.370.260.210.460.320.36Caryophyllene oxide158415830.270.410.381.840.710.840.932.001.691.56Cubenol164816460.100.52α -Cadinol165516560.080.723.580.06α-Bisabolol168716850.060.11Others——4.711.310.741.883.404.832.812.373.241.28Results are the average of two independent GC experiments. In order to reduce table complexity, the standard deviation values were ignored. RI represents the experimental retention index and RI∗ donates the published retention index, extracted from Adams book [11].


### 3.3. Antimicrobial Activities of Essential Oils

The MIC and MBC values for theT. daenensis oils are listed in Table 3. Generally, Gram-positive bacteria indicated high sensitivity against the T. daenensis oils antibacterial effects (P<0.001). Pseudomonas aeruginosa and Salmonella typhimurium are the most resistant and sensitive strains amongst the Gram-negative strains, respectively. Staphylococcus aurous and Bacillus subtilis indicated the most resistance and sensitivity against the oils antibacterial effects, in the Gram-positive strains. TD5, TD3, and TD10 extracts were the most potent antimicrobial oils, while TD2 was the least bioactive oil against all of the bacterial strains. Additionally, diversity in oils compositions led to differences in antifungal activity. The MIC values of T. daenensis extracts for Candida albicans and Aspergillus Niger were 0.13–0.58 mg/ml and around 0.08–0.1 mg/ml respectively. The MFC spectrum of the T. daenensis methanol extracts against the Candida albicans (0.13–0.58 mg/ml) and Aspergillus Niger (0.16–0.7 mg/ml) was determined using the adapted broth microdilution assay (Table 4).Table 3
Minimum inhibitory concentration (MIC, mg/ml) and minimum bactericidal concentration (MBC, mg/ml) values for ecotypicT. daenensis essential oils.
BacteriaAmpicillinTD1TD2TD3TD4TD5TD6TD7TD8TD9TD10Escherichia coli (Gram-negative)MIC0.002 ± 0.00a0.29 ± 0.00b,c,d0.39 ± 0.15c,d0.15 ± 0.00a,b0.15 ± 0.00a,b0.21 ± 0.08a,b,c0.44 ± 0.17d0.1 ± 0.08a,b0.14 ± 0.00a,b0.28 ± 0.00b,c,d0.30 ± 0.00b,c,dMBC0.064 ± 0.00a0.58 ± 0.00b0.79 ± 0.30b0.15 ± 0.00a0.3 ± 0.00a0.28 ± 0.00a0.58 ± 0.00b0.55 ± 0.00b0.28 ± 0.00a0.28 ± 0.00a0.59 ± 0.00bBacillus subtilis (Gram-positive)MIC<0.0005 ± 0.00a0.44 ± 0.17c,d0.53 ± 0.00d0.46 ± 0.18c,d0.3 ± 0.00b,c0.14 ± 0.00a,b0.29 ± 0.00b,c0.14 ± 0.00a,b0.14 ± 0.00a,b0.28 ± 0.00b,c0.30 ± 0.00b,c,dMBC<0.0005 ± 0.00a0.62 ± 0.08d1.05 ± 0.00e0.61 ± 0.00d0.3 ± 0.00c0.14 ± 0.00b0.29 ± 0.00c0.14 ± 0.00b0.14 ± 0.00b0.28 ± 0.00c0.59 ± 0.00dStaphylococcus aureus MRSA (Gram-positive)MIC0.14 ± 0.00a0.58 ± 0.00b,c1.05 ± 0.00d0.86 ± 0.18c,d0.45 ± 0.17a,b0.42 ± 0.16a,b0.44 ± 0.17a,b0.41 ± 0.16a,b0.42 ± 0.16a,b0.56 ± 0.00b,c0.30 ± 0.00a,bMBC>0.25 ± 0.00a0.84 ± 0.00c,d2.14 ± 0.00e1.06 ± 0.18d0.45 ± 0.17a,b0.88 ± 0.00d0.44 ± 0.17a,b0.41 ± 0.16a,b0.56 ± 0.00b,c0.56 ± 0.00b,c0.30 ± 0.00a,bSalmonella typhimurium (Gram-negative)MIC0.003 ± 0.001a0.15 ± 0.00b,c0.53 ± 0.00e0.15 ± 0.00b,c0.11 ± 0.04a,b,c0.07 ± 0.00a,b0.07 ± 0.00a,b1.09 ± 0.00f0.21 ± 0.08c,d0.28 ± 0.00d0.22 ± 0.09c,dMBC0.016 ± 0.00a0.30 ± 0.00c1.20 ± 0.00e0.28 ± 0.00c0.29 ± 0.00c0.14 ± 0.00b0.14 ± 0.00b1.56 ± 0.00f0.22 ± 0.09b,c0.56 ± 0.00d0.30 ± 0.00cEnterococcus faecalis (Gram-positive)MIC<0.0005 ± 0.00a0.15 ± 0.00b,c1.05 ± 0.00f0.15 ± 0.00b,c0.11 ± 0.04a,b0.14 ± 0.00a,b,c0.15 ± 0.00b,c0.55 ± 0.00e0.42 ± 0.16d,e0.28 ± 0.00c,d0.30 ± 0.00b,cMBC0.01 ± 0.00a0.58 ± 0.00d1.05 ± 0.00e0.15 ± 0.00b0.15 ± 0.00b0.14 ± 0.00b0.15 ± 0.00b0.55 ± 0.00d0.56 ± 0.00d0.28 ± 0.00c0.30 ± 0.00cStaphylococcus aureus (Gram-positive)MIC0.08 ± 0.01a0.58 ± 0.00a,b1.58 ± 0.61c0.91 ± 0.35b,c0.9 ± 0.35a,b,c0.56 ± 0.00a,b0.44 ± 0.17a,b0.82 ± 0.32a,b,c1.12 ± 0.00b,c0.84 ± 0.32a,b,c0.89 ± 0.34a,b,cMBC0.08 ± 0.05a0.58 ± 0.00a,b2.11 ± 0.00b,c0.91 ± 0.35a,b,c0.9 ± 0.35a,b,c0.56 ± 0.00a,b1.45 ± 1.00a,b,c2.73 ± 1.89c1.12 ± 0.00a,b,c1.12 ± 0.00a,b,c2.49 ± 0.00b,cPseudomonas aeruginosa (Gram-negative)MIC>0.25>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12MBC>0.25>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12According to the homogeneity of variances using an ANOVA analysis, values in the same row with common letters are not statistically different (P>0.05).Table 4
Minimum inhibitory concentration (MIC, mg/ml) and minimum fungicidal concentration (MFC, mg/ml) values forT. daenensis essential oils from different ecotypes.
FungiAmphotericin (μg/ml)TD1TD2TD3TD4TD5TD6TD7TD8TD9TD10Aspergillus nigerMIC0.82 ± 0.05a0.08 ± 0.00b0.08 ± 0.00b0.10 ± 0.00c0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00bMFC4.69 ± 1.80a0.36 ± 0.00c0.16 ± 0.00b0.38 ± 0.00c0.37 ± 0.00c0.34 ± 0.00c0.18 ± 0.00b0.17 ± 0.00b0.17 ± 0.00b0.70 ± 0.00d0.18 ± 0.00bCandida albicansMIC0.40 ± 0.00a0.58 ± 0.00h0.52 ± 0.00g0.30 ± 0.00f0.30 ± 0.00f0.13 ± 0.00b0.29 ± 0.00e,f0.27 ± 0.00c0.27 ± 0.00c,d0.28 ± 0.00d,e0.29 ± 0.00e,fMFC1.56 ± 0.02a0.58 ± 0.00h0.52 ± 0.00g0.30 ± 0.00f0.30 ± 0.00f0.13 ± 0.00b0.29 ± 0.00e,f0.27 ± 0.00c0.27 ± 0.00c,d0.28 ± 0.00d,e0.29 ± 0.00e,fAccording to the homogeneity of variances using an ANOVA analysis, values in the same row with common letters are not statistically different (P>0.05).


### 3.4. Antioxidative and DNA-Protecting Activities

Hydroperoxides presence can be spectrophotometrically evaluated using the beta-carotene/linoleic acid method. TD8 was the weakest bleaching inhibitor, but other capacities of this oil were remarkable and almost equal to the BHT or ascorbic acid (P>0.05). In contrast, hydrogen donor compounds are strong antioxidants and also associated with the DPPH-H formation. Table 5 represents the T. daenensis oils hydrogen donating potency. Although TD5 and TD6 presented similar and the most DPPH scavenging potency (P<0.01),TD2 was the weakest antioxidant (P<0.001) (Supplementary Figures S1 and S2). According to the comet assay analysis and also attained CIC50 (1.04 to 7.78 mg/ml), it can be concluded that the TD10 and TD9 would exhibit promising inhibitory effect on DNA damage and closer antigenotoxic activity versus BHT and ascorbic acid. The oil protecting properties on DNA rates as the following order: TD10 > TD9 > TD1 > TD3 > TD2 > TD5 > TD7 > TD8 > TD6 > TD4 (Table 5 and Supplementary Figure S3).Table 5
Thymus daenensis ecotypic oil values for antioxidant and antigenotoxic activities.
SampleVIT CBHTTD1TD2TD3TD4TD5TD6TD7TD8TD9TD10Bleaching inhibitory (BIC50 mg/ml)0.01 ± 0.00a0.03 ± 0.00a,b0.08 ± 0.01a,b0.08 ± 0.01a,b0.05 ± 0.01a,b0.02 ± 0.00a,b1.06 ± 0.03c0.01 ± 0.00a0.04 ± 0.01a,b0.18 ± 0.01b0.03 ± 0.00a,b0.04 ± 0.00a,bRadical scavenging (SC50 mg/ml)0.03 ± 0.00a0.02 ± 0.00a0.80 ± 0.05e1.79 ± 0.11g0.61 ± 0.03c,d0.52 ± 0.01b,c0.41 ± 0.01b0.44 ± 0.02b1.24 ± 0.05f0.71 ± 0.03d,e1.15 ± 0.05f0.65 ± 0.05c,dComet inhibitory concentration (CIC50 mg/ml)0.008 ± 0.00a0.02 ± 0.00a1.50 ± 0.32b>5.00e2.46 ± 1.00b,c3.34 ± 0.70c,d1.67 ± 0.47b>5.00e4.00 ± 0.50e,d4.28 ± 0.61c,d1.36 ± 0.54a,b1.04 ± 0.13a,bAccording to the homogeneity of variances using an ANOVA analysis, values in the same row with common letters are not statistically different (P>0.05). The estimated CIC50 values for TD2 and TD6 are 7.78 and 5.68 mg/ml, respectively.


## 3.1. Geographical Effects on the Physicochemical Properties of the Ecotypes

Different morphological analysesof Thymus daenensis ecotypes indicated a microenvironmental conditions distinct effect including altitude, humidity, and temperature on the plant height and also the total dry matter. This flowering plant is growing up to 35 cm tall with short trichomes on its leaves. The isolation process yields the yellowish oil with different productivity (0.36–0.86% v/w), based on the Thymus geographical origins. In addition, other ecotypic oils physical properties presented diversity including density and molar extinction coefficient at 280 nm, and it seems that they can also be directly affected by environmental conditions (Table 1).Table 1
Environmental conditions and physical properties ofThymus daenensis essential oils.
EcotypeEcotype locationCity, provinceAltitude (m)Temperature (°C)Average humidity (%)Density (g/ml) at 25°Cε280 nm (l/(cm × g))Refractive index at 25°CYield (v/w %)MaxMinTD134 59 17NGharechaie highlands, Marekazi190023.96400.94 ± 0.122.83 ± 1.011.32 ± 0.010.82 ± 0.2149 19 49ETD233 15 00NKhoramabad, Lorestan183038.8−1.6570.84 ± 0.041.98 ± 0.771.45 ± 0.080.36 ± 0.2848 30 00ETD333 25 00NKhoramabad, Lorestan190038.8−1.6570.98 ± 0.057.41 ± 2.321.38 ± 0.030.56 ± 0.1348 40 00ETD439 52 42NDaran, Isfahan250023−3.1420.96 ± 0.113.22 ± 1.741.39 ± 0.060.60 ± 0.2350 23 21ETD534 05 06NKharkan, Markazi196535.5−4.4310.88 ± 0.045.38 ± 0.821.44 ± 0.030.62 ± 0.1749 24 50ETD634 11 38NArak, Markazi2404273460.93 ± 0.062.92± 0.191.46 ± 0.090.86 ± 0.2549 29 41ETD732 55 52NFereidonshahr, Isfahan261521.3−4.5500.88 ± 0.027.99 ± 1.581.51 ± 0.100.59 ± 0.2350 07 01ETD833 46 24NRobatkhomein, Markazi240725.23510.89 ± 0.0811.43 ± 1.631.47 ± 0.040.69 ± 0.1249 52 57ETD936 33 15NZanjan, Zanjan195032.1−7.5540.90 ± 0.076.10 ± 1.941.51 ± 0.030.74 ± 0.2648 22 14ETD1036 30 00NChore, Ghazvin1950182510.95 ± 0.032.66 ± 0.191.48 ± 0.060.75 ± 0.1149 43 00E


## 3.2. Essential Oil Constituents in theThymus Ecotypes

Sixty-one compounds were identified by the GC/MS analysis, in differentT. daenensis ecotypes, and were representing 95 up to more than 99% of the total oil compositions. As displayed in Table 2, thymol (0.31–73.86%), carvacrol (5.14–51.89%), linalool (0.3–22.49%), pulegone (up to 17.57%), geraniol (2.02–14.84%), α-terpineol (0.58–14.01%), p-cymene (0.12–11.19%), o-cymene (5.19–8.05%), geranyl acetate (0.16–6.11%), neomenthol (up to 5.69%), β- caryophyllene (1.81–4.97%), terpineol (0.59–4.68%), and borneol (0.01–4.04%) are the major constituents in the T. daenensis volatile oil. Some compounds were only detected in special ecotypes, and they are identified as δ-3-carene, δ-elemene, α-farnesene, α-himachalene, isoborneol, isopulegol, Z-jasmone, Z-linalool oxide, 1-menthol, neomenthol, menthone, α- muurolene, piperitone, pulegone, and thymol acetate.Table 2
Essential oil compositions for ecotypicThymus daenensis (%).
CompoundsRIRI∗ [11]TD1TD2TD3TD4TD5TD6TD7TD8TD9TD10α-Pinene9389390.860.420.260.830.10Camphene9539540.140.080.36β-Pinene9789790.580.260.291-Octen-3-ol9809790.080.230.18β-Myrcene9909910.810.70α-Terpinene101710170.060.462.800.75p-Cymene102510260.122.565.126.5311.190.239.261.640.38O-Cymene102610275.198.05δ-3-Carene102910310.121,8-Cineole103010320.101.762.761.080.400.920.181.730.130.21Z-β-Ocimene103710380.20E-β-Ocimene104710500.210.18γ-Terpinene105910600.030.441.780.661.401.683.910.26Z-Linalool oxide107410750.34α-Terpinolene108710890.130.41Linalool109810971.5718.241.481.021.600.4422.490.322.010.30Z-β-Terpineol114411440.160.600.240.590.12Camphor114811502.880.260.480.62Menthone115211532.22Isopulegol115011522.12Isoborneol116111620.62Neomenthol116711685.69Borneol116911702.851.992.443.991.441.644.040.690.01Terpinen-4-ol117511771.040.604.680.700.591.100.90Isomenthol118111833.00α-Terpineol118911901.920.581.620.582.2814.010.61Thymol methyl ether123012350.220.12Pulegone1234123717.57Carvacrol methyl ether124312453.141.90Piperitone124912531.97Geraniol125312542.0214.84Thymol1303130460.610.3158.9464.5270.9560.600.7339.4525.0073.86Carvacrol130813107.4520.9014.865.148.826.2251.8916.2926.127.89δ-Elemene133713390.08Piperitenone134013431.90α-Terpinyl acetate134713491.100.080.41Thymol acetate135413520.15Eugenol135813590.390.05α-Copaene137613770.24Geranyl acetate138113811.850.166.110.94β-Bourbonene138913880.140.440.040.040.120.120.14β-Elemene139013930.08Z-Jasmone139313950.14Β-Caryophyllene142014194.563.733.084.601.812.323.504.973.403.60Aromadendrene144014410.160.200.53α-Himachalene145114520.35α-Humulene145614550.240.120.170.100.110.241.150.20γ-Muurolene148014820.310.220.290.250.390.78Germacrene D148114850.210.750.20Viridiflorene149214970.460.260.470.620.51α-Muurolene150115000.28α -Farnesene150415060.34β-Bisabolene150715061.470.400.820.672.650.830.660.851.24Cis-α-Bisabolene150815070.660.260.233.820.111.680.201.332.684.61γ-Cadinene151515140.410.330.260.160.150.150.28δ-Cadinene152415230.410.280.520.310.440.590.66Spathulenol157915780.380.650.370.260.210.460.320.36Caryophyllene oxide158415830.270.410.381.840.710.840.932.001.691.56Cubenol164816460.100.52α -Cadinol165516560.080.723.580.06α-Bisabolol168716850.060.11Others——4.711.310.741.883.404.832.812.373.241.28Results are the average of two independent GC experiments. In order to reduce table complexity, the standard deviation values were ignored. RI represents the experimental retention index and RI∗ donates the published retention index, extracted from Adams book [11].


## 3.3. Antimicrobial Activities of Essential Oils

The MIC and MBC values for theT. daenensis oils are listed in Table 3. Generally, Gram-positive bacteria indicated high sensitivity against the T. daenensis oils antibacterial effects (P<0.001). Pseudomonas aeruginosa and Salmonella typhimurium are the most resistant and sensitive strains amongst the Gram-negative strains, respectively. Staphylococcus aurous and Bacillus subtilis indicated the most resistance and sensitivity against the oils antibacterial effects, in the Gram-positive strains. TD5, TD3, and TD10 extracts were the most potent antimicrobial oils, while TD2 was the least bioactive oil against all of the bacterial strains. Additionally, diversity in oils compositions led to differences in antifungal activity. The MIC values of T. daenensis extracts for Candida albicans and Aspergillus Niger were 0.13–0.58 mg/ml and around 0.08–0.1 mg/ml respectively. The MFC spectrum of the T. daenensis methanol extracts against the Candida albicans (0.13–0.58 mg/ml) and Aspergillus Niger (0.16–0.7 mg/ml) was determined using the adapted broth microdilution assay (Table 4).Table 3
Minimum inhibitory concentration (MIC, mg/ml) and minimum bactericidal concentration (MBC, mg/ml) values for ecotypicT. daenensis essential oils.
BacteriaAmpicillinTD1TD2TD3TD4TD5TD6TD7TD8TD9TD10Escherichia coli (Gram-negative)MIC0.002 ± 0.00a0.29 ± 0.00b,c,d0.39 ± 0.15c,d0.15 ± 0.00a,b0.15 ± 0.00a,b0.21 ± 0.08a,b,c0.44 ± 0.17d0.1 ± 0.08a,b0.14 ± 0.00a,b0.28 ± 0.00b,c,d0.30 ± 0.00b,c,dMBC0.064 ± 0.00a0.58 ± 0.00b0.79 ± 0.30b0.15 ± 0.00a0.3 ± 0.00a0.28 ± 0.00a0.58 ± 0.00b0.55 ± 0.00b0.28 ± 0.00a0.28 ± 0.00a0.59 ± 0.00bBacillus subtilis (Gram-positive)MIC<0.0005 ± 0.00a0.44 ± 0.17c,d0.53 ± 0.00d0.46 ± 0.18c,d0.3 ± 0.00b,c0.14 ± 0.00a,b0.29 ± 0.00b,c0.14 ± 0.00a,b0.14 ± 0.00a,b0.28 ± 0.00b,c0.30 ± 0.00b,c,dMBC<0.0005 ± 0.00a0.62 ± 0.08d1.05 ± 0.00e0.61 ± 0.00d0.3 ± 0.00c0.14 ± 0.00b0.29 ± 0.00c0.14 ± 0.00b0.14 ± 0.00b0.28 ± 0.00c0.59 ± 0.00dStaphylococcus aureus MRSA (Gram-positive)MIC0.14 ± 0.00a0.58 ± 0.00b,c1.05 ± 0.00d0.86 ± 0.18c,d0.45 ± 0.17a,b0.42 ± 0.16a,b0.44 ± 0.17a,b0.41 ± 0.16a,b0.42 ± 0.16a,b0.56 ± 0.00b,c0.30 ± 0.00a,bMBC>0.25 ± 0.00a0.84 ± 0.00c,d2.14 ± 0.00e1.06 ± 0.18d0.45 ± 0.17a,b0.88 ± 0.00d0.44 ± 0.17a,b0.41 ± 0.16a,b0.56 ± 0.00b,c0.56 ± 0.00b,c0.30 ± 0.00a,bSalmonella typhimurium (Gram-negative)MIC0.003 ± 0.001a0.15 ± 0.00b,c0.53 ± 0.00e0.15 ± 0.00b,c0.11 ± 0.04a,b,c0.07 ± 0.00a,b0.07 ± 0.00a,b1.09 ± 0.00f0.21 ± 0.08c,d0.28 ± 0.00d0.22 ± 0.09c,dMBC0.016 ± 0.00a0.30 ± 0.00c1.20 ± 0.00e0.28 ± 0.00c0.29 ± 0.00c0.14 ± 0.00b0.14 ± 0.00b1.56 ± 0.00f0.22 ± 0.09b,c0.56 ± 0.00d0.30 ± 0.00cEnterococcus faecalis (Gram-positive)MIC<0.0005 ± 0.00a0.15 ± 0.00b,c1.05 ± 0.00f0.15 ± 0.00b,c0.11 ± 0.04a,b0.14 ± 0.00a,b,c0.15 ± 0.00b,c0.55 ± 0.00e0.42 ± 0.16d,e0.28 ± 0.00c,d0.30 ± 0.00b,cMBC0.01 ± 0.00a0.58 ± 0.00d1.05 ± 0.00e0.15 ± 0.00b0.15 ± 0.00b0.14 ± 0.00b0.15 ± 0.00b0.55 ± 0.00d0.56 ± 0.00d0.28 ± 0.00c0.30 ± 0.00cStaphylococcus aureus (Gram-positive)MIC0.08 ± 0.01a0.58 ± 0.00a,b1.58 ± 0.61c0.91 ± 0.35b,c0.9 ± 0.35a,b,c0.56 ± 0.00a,b0.44 ± 0.17a,b0.82 ± 0.32a,b,c1.12 ± 0.00b,c0.84 ± 0.32a,b,c0.89 ± 0.34a,b,cMBC0.08 ± 0.05a0.58 ± 0.00a,b2.11 ± 0.00b,c0.91 ± 0.35a,b,c0.9 ± 0.35a,b,c0.56 ± 0.00a,b1.45 ± 1.00a,b,c2.73 ± 1.89c1.12 ± 0.00a,b,c1.12 ± 0.00a,b,c2.49 ± 0.00b,cPseudomonas aeruginosa (Gram-negative)MIC>0.25>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12MBC>0.25>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12>5.12According to the homogeneity of variances using an ANOVA analysis, values in the same row with common letters are not statistically different (P>0.05).Table 4
Minimum inhibitory concentration (MIC, mg/ml) and minimum fungicidal concentration (MFC, mg/ml) values forT. daenensis essential oils from different ecotypes.
FungiAmphotericin (μg/ml)TD1TD2TD3TD4TD5TD6TD7TD8TD9TD10Aspergillus nigerMIC0.82 ± 0.05a0.08 ± 0.00b0.08 ± 0.00b0.10 ± 0.00c0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00b0.08 ± 0.00bMFC4.69 ± 1.80a0.36 ± 0.00c0.16 ± 0.00b0.38 ± 0.00c0.37 ± 0.00c0.34 ± 0.00c0.18 ± 0.00b0.17 ± 0.00b0.17 ± 0.00b0.70 ± 0.00d0.18 ± 0.00bCandida albicansMIC0.40 ± 0.00a0.58 ± 0.00h0.52 ± 0.00g0.30 ± 0.00f0.30 ± 0.00f0.13 ± 0.00b0.29 ± 0.00e,f0.27 ± 0.00c0.27 ± 0.00c,d0.28 ± 0.00d,e0.29 ± 0.00e,fMFC1.56 ± 0.02a0.58 ± 0.00h0.52 ± 0.00g0.30 ± 0.00f0.30 ± 0.00f0.13 ± 0.00b0.29 ± 0.00e,f0.27 ± 0.00c0.27 ± 0.00c,d0.28 ± 0.00d,e0.29 ± 0.00e,fAccording to the homogeneity of variances using an ANOVA analysis, values in the same row with common letters are not statistically different (P>0.05).


## 3.4. Antioxidative and DNA-Protecting Activities

Hydroperoxides presence can be spectrophotometrically evaluated using the beta-carotene/linoleic acid method. TD8 was the weakest bleaching inhibitor, but other capacities of this oil were remarkable and almost equal to the BHT or ascorbic acid (P>0.05). In contrast, hydrogen donor compounds are strong antioxidants and also associated with the DPPH-H formation. Table 5 represents the T. daenensis oils hydrogen donating potency. Although TD5 and TD6 presented similar and the most DPPH scavenging potency (P<0.01),TD2 was the weakest antioxidant (P<0.001) (Supplementary Figures S1 and S2). According to the comet assay analysis and also attained CIC50 (1.04 to 7.78 mg/ml), it can be concluded that the TD10 and TD9 would exhibit promising inhibitory effect on DNA damage and closer antigenotoxic activity versus BHT and ascorbic acid. The oil protecting properties on DNA rates as the following order: TD10 > TD9 > TD1 > TD3 > TD2 > TD5 > TD7 > TD8 > TD6 > TD4 (Table 5 and Supplementary Figure S3).Table 5
Thymus daenensis ecotypic oil values for antioxidant and antigenotoxic activities.
SampleVIT CBHTTD1TD2TD3TD4TD5TD6TD7TD8TD9TD10Bleaching inhibitory (BIC50 mg/ml)0.01 ± 0.00a0.03 ± 0.00a,b0.08 ± 0.01a,b0.08 ± 0.01a,b0.05 ± 0.01a,b0.02 ± 0.00a,b1.06 ± 0.03c0.01 ± 0.00a0.04 ± 0.01a,b0.18 ± 0.01b0.03 ± 0.00a,b0.04 ± 0.00a,bRadical scavenging (SC50 mg/ml)0.03 ± 0.00a0.02 ± 0.00a0.80 ± 0.05e1.79 ± 0.11g0.61 ± 0.03c,d0.52 ± 0.01b,c0.41 ± 0.01b0.44 ± 0.02b1.24 ± 0.05f0.71 ± 0.03d,e1.15 ± 0.05f0.65 ± 0.05c,dComet inhibitory concentration (CIC50 mg/ml)0.008 ± 0.00a0.02 ± 0.00a1.50 ± 0.32b>5.00e2.46 ± 1.00b,c3.34 ± 0.70c,d1.67 ± 0.47b>5.00e4.00 ± 0.50e,d4.28 ± 0.61c,d1.36 ± 0.54a,b1.04 ± 0.13a,bAccording to the homogeneity of variances using an ANOVA analysis, values in the same row with common letters are not statistically different (P>0.05). The estimated CIC50 values for TD2 and TD6 are 7.78 and 5.68 mg/ml, respectively.


## 4. Discussion

Natural plant products were the most of the medicinal agent source traditionally and were used in order to provide enthusiasm for industrial drug discovery.Thymus genera are famous in folk medicine. The ecotypes of T. daenensis essential oils compositions have been evaluated, and approximately 62 constituents were identified in order to form 95–100% of the total volatile oil spectra, generally. Thymol and carvacrol were the main constituents with average amounts of 45.5% and 16.5%, respectively. Table 2 indicates that the oxygenated compounds (>80%) are more dominant in comparison with the hydrocarbon constituents (<20%). Oxygenated monoterpenes were the most abundant (60–84%) consisting of carvacrol, linalool, terpineol, and thymol. Monoterpene hydrocarbons, oxygenated sesquiterpenes, and sesquiterpene hydrocarbons were present in low amounts amongst the ecotypes.The results demonstrated that the total phenolics contents strongly depended on the plant species, geographical origin, and genetic content. Tables1 and 2 are in accordance with the earlier research data, which stated that the high temperature with low altitude would increase the phenolic components (thymol and carvacrol) ratio in the oils [17]. In addition to these observations, it was also found that the ecological properties like high relative humidity would significantly increase the nonphenolic terpenes (linalool) and noncyclic monoterpenes (geraniol) ratios in some ecotypes. Despite these correlations, some components of T. daenensis oils variation could not be interpreted only by the environmental conditions, and other factors like genetic background may affect the type and amount of the constituents. For instance, this compound was not detected in other ecotypes, even those with similar environmental conditions, while the content of noncyclic monoterpenes pulegone has increased in TD2 as an ecotype growing in high temperature [18].Due to the other researches, while TD2 with the lowest phenolic (thymol and carvacrol) content, near 21%, indicated the weakest antibacterial effect in all strains, TD3, TD5, and TD10 were the most effective antibacterial oils with high phenols levels (74%, 80%, and 82%, respectively). The microorganism susceptibility to essential oil derivatives is also microorganism-dependent. As a conclusion, these phenolic components hydrophobic properties enable them to partition with the bacterial cell membrane lipids, and lead to disturbing the cell integrity, and also rendering them more permeable. Alkyl groups, hydroxyls, or aldehyde moieties would influence the phenols antimicrobial activity, possibly through interfering with cell wall enzymes and reducing the surface tension [19, 20].AlthoughAspergillus niger was more sensitive to the oils cytostatic effect (P<0.001), but the T. daenensis ecotypes indicated similar fungicidal effect on both fungi strains (P>0.05). Based on the MIC data, the differences between the oils’ cytostatic effects on the Aspergillus Niger were not statistically significant (P>0.05). However, MFC data indicated that TD9 had less fungicidal effects in comparison with the other ecotypes (P>0.05). TD2 and TD1 were the less bioactive agents against Candida albicans (P<0.001), while TD4 and TD5 presented the most antifungal effects against Candida albicans (P<0.001). Consistent with many publications, variations in the fungicidal effects were not consistent with the ecotypic oils phenolic contents; for example, TD1 with high thymol and carvacrol (68%) content presented a similar and even less antifungal activity compared to the TD2, TD7, TD8, and TD9 [21–23].Large natural polyphenols are ordered from simple to complex molecules according to the phenol rings and structural ingredients number and arrangement [17]. They act as antioxidant agents by donating hydrogen from their phenolic hydroxyl groups in the food industry. It is believed that the hydroxyl moieties number, position, and the degrees of glycosylation, esterification, and polymerization would determine the polyphenol components’ antioxidant activity. According to many researches, thymol and carvacrol are the most effective constituents in lipid peroxidation. A survey conducted on the different Thymus species has presented that the effective antioxidant derivatives’ order is as follows: thymol > carvacrol > terpinene > myrcene > linalool > cymene > limonene > cineole > pinene [24, 25]. DPPH radical scavenging activity and β-carotene supplementation results also revealed that the T. daenensis ecotypes exhibited dose-dependent antioxidants scavenging properties, acting as free radical terminators. TD2 represented the highest SC50 and BIC50 value, due to the fact that it contains the least effective phenolic antioxidants percent. Although some researchers reported that phenolic compounds’ high concentrations might induce DNA single/double strand breaks, this study’s results indicated that the T. daenensis oil thymol and carvacrol contents could act as oxidative DNA damage effective inhibitors [26]. In this regard, TD10 were the best antigenotoxic extracts and TD2 with the lowest phenolic compounds level represented the least DNA protection ability.


## 5. Conclusion

TheThymus essential oil chemical composition was affected by plant microenvironment variation directly. In conclusion, the most biological aspects regarding T. daenensis antimicrobial, antifungal, and antioxidant activities are associated with the oil constituents. Although direct correlations between the major constituents and biological activities have not been found in some cases, the total phenolic contents are responsible for most behaviors. It seems that the presence of some trace components in the oils may demonstrate severe antagonistic or synergistic behaviors and lead to unpredictable biological functions.



---

*Source: 6686558-2021-02-10.xml*