Chemical Profiles and Insecticidal Potential of Essential Oils Isolated from Four Thymus Species against Rhyzopertha dominica (F.)

Although chemical pesticides have been efficiently used to manage insect pest, their overuse has led to environmental contamination and threats to human health, enticing researchers to introduce eco-friendly and effective agents. In this study, the insecticidal effectiveness of essential oils isolated from Thymus species, including T. eriocalyx, T. kotschyanus, T. fallax, and T. vulgaris, was evaluated against the adults of Rhyzopertha dominica. The terpenes p-cymene, 1,8-cineole, linalool, α-terpineol, and carvacrol were the prominent compounds in the hydrodistilled essential oils. All essential oils produced significant fumigant at 24, 48, and 72-exposure times. The energy reserves protein by all essential oils, glycogen by T. kotschyanus and T. vulgaris, and lipid by T. fallax and T. vulgaris were significantly decreased compared to control. All essential oils except T. vulgaris affected the amylolytic and proteolytic activity of the pest. The pest increased the α- and β-esterase enzyme activity in response to the essential oils. Nutritional indices of adults were also affected by essential oils, in which feeding deterrence index was calculated from 20.41% to 61.11%. Accordingly, based on lethal and extensive sub-lethal insecticidal activities, T. eriocalyx, T. kotschyanus, T. fallax, and T. vulgaris essential oils can be considered as efficient agents for R. dominica management.


Introduction
The lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), is one of the damaging primary insect pests of stored grains especially wheat and rice throughout most of the world [1,2]. The adults have high viability potential and are strong flyers that can find the hosts from up to 1 km distance [3]. Rhyzopertha dominica is an internal feeder, and the whole life of larvae and most of the adult life, as harmful stages of the pest, are spend inside a single kernel [4]. As fumigants can penetrate areas where insects hide, fumigation of an active substance is the main strategy in the management of stored-product insect pests [5]. Although synthetic fumigants such as phosphine (PH 3 ) have a very good penetrant feature and quick lethality actions, frequent application of these chemicals has resulted in several side effects including environmental contamination, threat to human health and non-target organisms, and developing insect pest resistance [6,7].
The promising insecticidal effects of plant-derived essential oils fumigated on several stored product insect pests have been reported [8][9][10]. Recently, the insecticidal activities of essential oils against R. dominica were reported: The significant fumigant toxicity of essential oil extracted from Satureja intermedia C. A. Mey was revealed against the adults

Fumigant Toxicity
Data of the mortality of R. dominica adults treated by the essential oils of T. eriocalyx Probit analysis of data obtained from the fumigant toxicity of Thymus essential oils against R. dominica adults revealed that a concentration of 180.9, 157.0, 174.5, and 174.7 µL/L, respectively, from T. eriocalyx, T. fallax, T. kotschyanus, and T. vulgaris essential oils was adequate to achieve 90% mortality of tested insects (LC 90 values) after 72 h. According to the high r 2 values accessible in Table 2, there is a positive direct relation between tested concentrations of Thymus essential oils and the mortality of R. dominica. Although the calculated LC 50 values are close to each other, the essential oil of T. kotschyanus with the lowest LC 50 and the highest relative potency at all exposure times showed higher toxicity than others (Table 2). On the other hand, the adults of R. dominica were more susceptible to the T. kotschyanus essential oil than the other species. Calculated 24-h LC 30 values, including 84. 16, 85.36, 75.57, and 103.23 µL/L for T. eriocalyx, T. fallax, T. kotschyanus, and T. vulgaris essential oils, respectively, were selected for sublethal experiments.

Effect of Essential Oils on Energy Reserves
Variations in protein, glycogen and lipid content of R. dominica adults treated by LC 30 and LC 50 of T. eriocalyx, T. fallax, T. kotschyanus, and T. vulgaris essential oils are shown in Table 3. The LC 30 values had no significant effect on the protein, glycogen and lipid content, but with increasing concentration to LC 50 , the content of these macromolecules decreased. The treatment of adults with LC 50 values of all tested essential oils significantly reduced their total protein content compared to the control (p < 0.05). Although the glycogen content of adults decreased by all essential oils, only the reduction affected by T. kotschyanus and T. vulgaris essential oils was significant compared to the control group (p < 0.05). Regarding lipid content, the more decrease delivered by T. fallax essential oil, followed by T. eriocalyx and T. vulgaris essential oils (Table 3).

Effect of Essential Oils on Esterase Activity
Esterase enzyme activity of R. dominica adults treated with LC 30 and LC 50 of T. eriocalyx, T. fallax, T. kotschyanus, and T. vulgaris essential oils was significantly increased compared to the control group (p < 0.05) ( Table 4). The activity of α-esterase enzyme by T. kotschyanus and T. vulgaris treatment was higher than T. eriocalyx and T. fallax. The β-esterase activity in the adults treated by the essential oils of T. fallax and T. vulgaris had greater increases than T. eriocalyx and T. kotschyanus (Table 4).

Effect of Essential Oils on Amylolytic and Proteolytic Activity
Although the LC 30 of essential oils had no considerable effect on the amylase and protease activities, treatment of R. dominica adults with LC 50 values of T. eriocalyx, T. kotschyanus and T. fallax essential oils significantly decreased the activity of these digestive enzymes (p < 0.05). The effect of T. vulgaris essential oil on the activity of amylase and protease enzymes was not significantly different compared to the control (p < 0.05). Essential oils of T. kotschyanus and T. fallax caused the most decrease in amylase enzyme activity. T. kotschyanus essential oil also reduced the protease enzyme activity more than other essential oils (Table 5).
Mean values in a column followed by different lowercase letter are significantly different on the basis of ANOVA with LSD test (p < 0.05).

Effect of Essential Oils on Nutritional Indices
The effects of LC 30 and LC 50 values of Thymus essential oils on nutritional indices of R. dominica adults are shown in Table 6. The consumption index (CI) of adults treated by both LC 30 and LC 50 of T. eriocalyx, T. kotschyanus, and T. vulgaris essential oils was significantly reduced compared to the untreated insects (p < 0.05). Both LC 30 and LC 50 of T. fallax essential oil had no significant effect on CI of adults. The efficiency of conversion of ingested food (ECI) was significantly decreased by LC 30 and LC 50 of T. fallax and T. kotschyanus (p < 0.05). Treatment by LC 30 and LC 50 values of T. eriocalyx, T. kotschyanus, and T. vulgaris essential oils resulted in a significant reduction in the relative consumption rate (RCR) of R. dominica adults (p < 0.05). While the LC 30 of T. fallax essential oil had no significant effect on the RCR, it was significantly reduced by the LC 50 value compared to the control (p < 0.05). Although tested concentrations of all essential oils significantly decreased the relative growth rate (RGR) (p < 0.05), the lowest RGR was observed in the adults treated with T. kotschyanus essential oil. The essential oils of T. vulgaris and T. kotschyanus showed the highest feeding deterrence index (FDI), with both LC 30 and LC 50 , followed by T. eriocalyx and T. fallax oils (Table 6).

Discussion
Chemical composition of T. eriocalyx, T. fallax, T. kotschyanus, and T. vulgaris essential oils was investigated in some recent studies. An extensive examination of T. eriocalyx essential oils from different locations in Iran showed a high degree of variation in the essential oil compositions [23,24], as observed in other countries on Thymus and other genera of Lamiaceae family [25]. They described three chemotypes, (1) a linalool chemotype, (2) a thymol chemotype, and (3) a geraniol chemotype. The T. eriocalyx essential oil in this work represents a chemotype that may be the best described as a thymol/carvacrol mixed chemotype, more closely resembling sample T2 from Sefidkon et al. [24] and a T. eriocalyx sample analyzed by Amiri [26]. Based on previous reports on T. fallax essential oil compositions, there are two general chemotypes: (1) a thymol-rich chemotype and (2) a carvacrol-rich chemotype. However, The T. fallax essential oil sample in this study does not fit with either of those two chemotypes; both thymol (12.1%), and carvacrol (2.4%) concentrations were low, and this essential oil sample may represent a new chemotype, rich in linalool [27][28][29][30]. Based on the essential oil composition, agglomerative hierarchical cluster analysis has revealed that three chemotypes of T. kotschyanus can be defined: (1) a carvacrolrich chemotype (>80% similarity within samples), (2) a thymol-rich chemotype (around 90% similarity within samples), and (3) a mixed chemical type dominated by α-terpinyl acetate, α-terpineol, linalyl acetate, and linalool, which includes the essential oil in this work (see Figure 1) and is completely dissimilar to the carvacrol and thymol chemotypes. An analysis of numerous T. vulgaris essential oils reported in the literature had revealed four major chemotypes, (1) a thymol-rich chemotype, (2) a carvacrol-rich chemotype, (3) a linalool-rich chemotype, and (4) a geraniol-rich chemotype [31]. Most T. vulgaris samples from Iran fall into the thymol-rich chemotype [32][33][34][35]. However, the T. vulgaris essential oil in this study belongs to a geraniol-rich chemotype. It can also be noted that differences in exogenous and endogenous factors affecting the chemical profiles of Thymus essential oils, including geographical location, genetic make-up and growth stages of plants, and extraction method [24,31,36,37], have resulted in the above-mentioned differences.
in this work (see Figure 1) and is completely dissimilar to the carvacrol and thymol chemotypes. An analysis of numerous T. vulgaris essential oils reported in the literature had revealed four major chemotypes, (1) a thymol-rich chemotype, (2) a carvacrol-rich chemotype, (3) a linalool-rich chemotype, and (4) a geraniol-rich chemotype [31]. Most T. vulgaris samples from Iran fall into the thymol-rich chemotype [32][33][34][35]. However, the T. vulgaris essential oil in this study belongs to a geraniol-rich chemotype. It can also be noted that differences in exogenous and endogenous factors affecting the chemical profiles of Thymus essential oils, including geographical location, genetic make-up and growth stages of plants, and extraction method [24,31,36,37], have resulted in the above-mentioned differences.  The insecticidal efficiency of essential oils was related to their compounds [9,46]. Rozman et al. [47] showed that the monoterpenes 1,8-cineole, carvacrol, and linalool occurring in T. vulgaris essential oil, and recognized as main compounds of all essential oil studied in our study, had promising fumigant toxicity against the adults of R. dominica. Fumigant toxicity of p-cymene and α-terpineol, the other main compounds in the present study, against the adults of R. dominica was also demonstrated by Shaaya et al. [48] and Filomeno et al. [49]. Rozman et al. [47] also concluded that linalool was highly effective than 1,8-Cineole and carvacrol. Although linalool was identified in T. eriocalyx (3.2%) and T. vulgaris (4.1%) essential oils, its level was higher in the essential oils of T. fallax (15.4%), T. kotschyanus (9.5%). Furthermore, according to the study of Rozman et al. [47], linalyl acetate had significant fumigant toxicity against R. dominica adults, which was recognized in both T. fallax (0.6%) and T. kotschyanus (7.9%) essential oils, and these can justify the toxicity of T. fallax and T. kotschyanus essential oils, which were higher than T. eriocalyx and T. vulgaris.
Reducing the amount of energy resources the proteins, carbohydrates, and lipids considered one of the main strategies in the insect pest management because of their important roles in insect biochemical pathways, growth, metamorphosis, reproduction, and diapause [50][51][52]. On the other hand, changes in energy reserves specify the susceptibility of the insect to tested materials [53]. In this study, though reductions in glycogen content by T. eriocalyx and T. fallax essential oils and in lipid content by T. kotschyanus essential oil were not significantly different than controls (p < 0.05), the energy content of R. dominica adults in other treatments was reduced. These results agreed with Yazdani et al. [54], Jyoti et al. [55], and Moutassem et al. [18] who discovered the reduction in protein, carbohydrate, and lipid contents of insect pests treated with Thymus essential oils. Insects normally are able to convert carbohydrates into lipids [53], which could explain the relatively lower reduction in lipid content due to the essential oil of T. kotschyanus compared to other essential oils.
Exposure to lethal and sublethal concentrations considerably affected the enzymatic activities of an organism, reflecting the biochemical disturbances [56]. Hence, to clarify the mechanism responsible for the toxicity of Thymus essential oils, effects on digestive amylase and protease enzymes and detoxifying αand β-esterase activities were also studied. In this study, the activity of digestive enzymes amylase and protease was reduced in R. dominica adults treated by T. eriocalyx, T. kotschyanus, and T. fallax essential oils. Similar supportive results have been reported in the case of reducing the digestive enzyme activities in insect pests treated by essential oils [57][58][59], showing strong insecticidal action. It may be due to cytotoxic effects on epithelial cells of alimentary canal which reduced the secretion rate of their enzymes [52,59].
Treatments with the LC 30 and LC 50 of T. eriocalyx, T. fallax, T. kotschyanus, and T. vulgaris essential oils significantly increased the activity of the αand β-esterases compared to control. The esterases can hydrolyze the ester bonds in a toxic substance and detoxify them in the insect's body [60]. The high activity of esterase enzymes indicates that the insect's defense mechanisms are active against the treating agents and can be consider as development of resistance [61,62]. In general, regardless of the type of influencing exogenous factor, the insect activates its immune or defense mechanisms, such as increasing detoxifying enzymes, which will require significant energy or expense [53]. In this study, the decrease in protein, lipid, and glycogen content of R. dominica treated with Thymus essential oils was probably due to the pest's greater need for these energy sources to synthesize esterase enzymes. Indeed, based on the multiple modes of action of essential oils, including the inhibition of acetylcholinesterase, adenosine triphosphatases (ATPases) and glutathione-S-transferase activity and disruption in gamma-aminobutyric acid (GABArs) and octopamine receptors [63], along with reduction in energy reserves and digestive enzyme activity and antinutritional effects established in this study, the chances of the pest resistance will be low. In other words, the co-evolution of herbivorous insects and plants continues, and according to this study, Thymus essential oils can overcome the herbivorous insect R. dominica.
The nutritional indices, including consumption index (CI), efficiency of conversion of ingested food (ECI), relative consumption rate (RCR), and relative growth rate (RGR), of R. dominica adults treated with Thymus essential oils particularly T. kotschyanus essential oil were changed. The feeding deterrence index (FDI) of the both LC 30 and LC 50 of essential oils from 20.41% to 66.11% was also considerable. The nutritional indices describe the utilization of diets and or effectiveness of digestion that how the insect can convert foods or their own biomass [54]. Decreases in nutritional indices of insects can lead to growth retardation due to reduce in energy resources, fertility, and fecundity, and even life stage longevity [64]. In agreement with our findings, ECI, RCR, and RGR of R. dominica were decreased by the essential oils of Eucalyptus dundasii Maiden [65], E. floribendi Hugel ex Endi [66], Ferula assafoetida L., Juglans regia L., and Pelargonium hortorum Bailey [67]. The antinutritional effects of Thymus essential oils were also investigated by the some recent study. For example, a reduction in ECI, RCR, and RGR along with 57.83% FDI by 2 µL/mg food of the larvae of Indian meal moth (Plodia interpunctella (Hübner)) treated by T. daenensis Celak were reported [68]. In another study, ECI, RCR, and RGR of the larvae of Colorado potato beetle (Leptinotarsa decemlineata (Say)) treated with T. daenensis essential oil, rich in thymol (72.3%), carvacrol (7.1%), and γ-terpinene (4.8%), were significantly decreased compared to the control group [69]. Furthermore, 14 ppm of essential oil caused about 40% FDI for treated larvae. Reduction of nutritional indices of the mentioned pests by T. daenensis essential oil is in line with the results of this study. However, differences in FDI values may be due to differences in insects tested, essential oils used, and their chemical components. Strong antifeedant effects of thymol-rich chemotype essential oils of T. vulgaris and T. zygis Loefl. ex L. (with FDI of 74.9% and 72.3%, respectively) against the larvae of Egyptian cotton leaf worm Spodoptera littoralis (Boisduval) were reported by Valcárcel et al. [70]. They also indicated that the terpenes α-pinene, limonene, linalool, 1,8-cineole, camphor, p-cymene, carvacrol, thymol, and menthone had significant antifeedant effects, in which α-pinene, carvacrol, and thymol with high FDI (67.3%, 55.8%, and 52.4%, respectively) were more effective. Almost all above-mentioned terpenes were recognized in Thymus essential oils investigated in the present study, and accordingly, it can be said that their antinutritional effects may be related with such terpenes. Menthone, despite in low percentage, existed only in T. vulgaris and T. kotschyanus essential oils (0.4% and 0.5%, respectively), which can have an effect on the higher FDI of these essential oils than others. Furthermore, the more FDI of T. vulgaris than T. kotschyanus may be related to the high percentage of thymol in this oil which was not identified in T. kotschyanus.

Plant Materials and Extraction of Essential Oils
The Thymus eriocalyx, T. fallas, T. kotschyanus, and T. vulgaris plants were collected from Malayer (34 • 17 50" N 48 • 49 23" E Hamedan province), Saqqez (36 • 12 20" N 46 • 13 19" E Kurdistan province), Sardabeh (38 • 17 01" N 48 • 02 07" E Ardabil province), and Ardabil (38 • 14 46" N 48 • 14 26" E Ardabil province) of Iran, respectively (Geographic Coordinates from Google Maps). Thymus species were recognized according to the keys reported by Jamzad [71]. The voucher specimens were deposited with their scientific names at the Department of Plant Sciences, Moghan College of Agriculture and Natural Resources, University of Mohaghegh Ardabili, Ardabil, Iran. Five cm aerial parts including stems, leaves and flowers from the top of each separate plant were air dried in the room temperature. The specimens were pulverized with an electric grinder. The extraction of essential oils was carried out using a Clevenger apparatus with a 2000-mL flask, 100 g of the each pulverized specimen, and 1200 mL distilled water over a 3-h period. Anhydrous sodium sulphate was used to remove water from obtained essential oils, which were reserved in glass containers covered by aluminum foil and stored in refrigerator at 4 • C.

Insect Rearing
A colony of R. dominica was obtained from the Department of Plant Protection, University of Mohaghegh Ardabili, Iran. Insect rearing was performed in cylindrical plastic containers, in which the openings were covered by a mesh cloth for aeriation. The adults were released on crushed wheat (200 g, Aftab cultivar) into the containers which were kept in a growth chamber at 28 ± 1 • C, 60 ± 5% relative humidity, and 14:10 h light: dark conditions. The synchronized 1-d-old adults were selected for the experiments.

Fumigant Toxicity
Twenty 1-d-old adults were transferred to glass containers (140 mL) as fumigant chambers. The calculated concentrations (78.54, 86.61, 99.46, 114.24, 131.16, and 149 24, 125.24, 140.52, and 157.08 µL/L for T. vulgaris) were poured on filter paper discs (3 cm in diameter) using a micropipette. Based on the preliminary tests, the lower and upper concentrations causing about 25 to 75% mortality of treated insects, respectively, were selected, and the others were calculated via logarithmic intervals [73]. The filter papers treated with each essential oil concentration were affixed to the inner surface of the glass container lids using adhesive tape. The lids of glass containers were rendered impermeable to air using parafilm. Mortality of treated adult insects was recorded after 24, 48, and 72 h exposure times. All steps were considered for control groups without adding essential oil concentrations or any solvent, and the experiments were repeated three times. Bovine serum albumin and Fast Blue RR salt were purchased from Roche Co. (Germany) and Merck Co. (Germany), respectively. Other biochemical materials from reagents to enzyme substrates were purchased from Sigma Chemical Co. (USA). Treatments by the essential oil concentrations and the control were repeated three times in all biochemical assays.

Effect of Essential Oils on Energy Reserves
The whole bodies of 100 1-d-old adults treated with LC 30 and LC 50 of essential oils for 24 h were homogenized using a hand-held glass homogenizer in in 250 µL of sodium phosphate buffer (0.04 M, pH 7.0) at 4 • C, and the resulting homogeneous mixtures were centrifuged at 10,000 rpm for 10 min (Sigma 1-14 K refrigerated centrifuge, USA). The lipid content of adult insects was measured using vanillin reagent according to the method described by Van Handel [74]. The protein content was estimated using bovine serum albumin as a standard, according to the Bradford method [75]. The enthrone reagent was used to estimate the glycogen content [76], and the absorption was recorded spectrophotometrically at 626 nm (Unico, UV/Vis 2100, Dayton, NJ, USA).

Effect of Essential Oils on Detoxifying Esterase Enzymes
The whole bodies of fifty alive adults treated with essential oil concentrations were homogenized in 250 µL of sodium phosphate buffer (0.04 M, pH 7.0) on ice. The homogeneous mixtures were centrifuged at 10,000 rpm for 15 min, and the resulting supernatant was separated and refrigerated at −20 • C as an enzymatic extract. The activity of αand β-esterases was measured by the following procedures [77]: The enzyme extract (12.5 µL) was mixed with 112.5 µL of sodium phosphate buffer and 50 µL of α-naphthyl acetate (α-NA) and β-naphthyl acetate (β-NA), respectively, for αand β-esterase. The mixtures were incubated at 30 • C for 15 min, and then50 mg of Fast Blue RR salt in 50 mL of sodium phosphate buffer (0.075 M, pH 7.0) was added. The absorbance was recorded for α-NA and β-NA during 7 min intervals at 450 and 540 nm, respectively, using a microplate reader (ELIZA-Reader, Anthos 2020, UK).

Effect of Essential Oils on Digestive Amylase and Protease Enzymes
To measure the amylase enzyme activity, 20 µL of the enzyme extract, 500 µL of acetate buffer (50 mM, pH 6.0), and 40 µL of 1% soluble starch were reacted at 37 • C within 30 min. Then, 100 µL of dinitro salicylic acid reagent (DNSA: Sigma Chemical Co., St. Louis, USA) was added, boiled in a water bath for 10 min, and centrifuged at 15,000 rpm at 4 • C for Plants 2022, 11, 1567 13 of 17 5 min. Three replications were considered for each treatment and blank and the absorbance was recorded at 540 nm. The amount of enzyme required to produce 1 mg of maltose at 37 • C was considered as amylase activity unit [78].
To measure the protease activity, azocasein substrate was used based on the method of Elpidina et al. [79]. Twenty µL of enzyme extract was mixed with 80 µL of azocasein solution (1.5%) in 2-morpholinoethanesulfonic acid buffer (MES: 50 mM, pH 6.0), and incubated at 37 • C for 50 min. The enzymatic reaction was blocked by adding 100 µL of 30% trichloroacetic acid (TCA). The unhydrolyzed azocasein was then precipitated by refrigeration at 4 • C within 0.5 h, and the mixture was centrifuged at 15,000 rpm for 10 min. One hundred µL of the supernatant was added to 100 µL of 2 M sodium hydroxide and the absorbance was recorded at 440 nm. Three replications were considered for each treatment and blank. In the blank, the enzyme extract was added to the reaction mixture after adding 30% TCA. A change in optical density per milligram protein per minute was defined as the protease activity unit.

Effect of Essential Oils on Nutritional Indices
To evaluate the antinutritional effects of Thymus essential oils against R. dominica adults, 200 1-d-old adults were separately treated with LC 30 and LC 50 of the essential oils. For each concentration, the alive insects were divided into 7 groups (10 insects for each). The insects were transferred to 6-cm Petri dishes containing 3 g of crushed wheat grains (Aftab cultivar). After two weeks, the food consumption, the adults' weight, and the weight gain of insect pest were determined. The initial and remaining food and the fed and unfed adults were separately weighed after two weeks. To determine the percentage of dry weight of food and R. dominica, twenty samples were weighed, dried in an oven (60 • C for 48 h) and then re-weighed (Sartorius AG Germany GCA803S, d = 0.001 ct). Nutritional indices, including Consumption Index (CI), Relative Consumption Rate (RCR), Relative Growth Rate (RGR), and Efficiency of Conversion of Ingested food (ECI) were calculated based on the formulae described by Waldbauer [80]: CI = F/A, RCR = F/TA, RGR = G/TA, and ECI= G/F, Where, A = mean dry weight of insect during feeding period (mg), G = dry weight gain of insect during feeding period (mg), F = dry weight of food eaten (mg), and T = feeding period (day). Feeding deterrence index (FDI) was calculated as (FDI) (%) = [(C-T)/C] × 100, where C = the mean weight of food eaten in control and T = the mean weight of food eaten in treatment [81].

Statistical Analysis
For the agglomerative hierarchical cluster (AHC) analysis, the chemical compositions of 20 T. kotschyanus essential oils from Iran were treated as operational taxonomic units (OTUs). The percentage compositions of 11 major components (p-cymene, 1,8-cineole, linalool, borneol, α-terpineol, geraniol, geranial, thymol, carvacrol, α-terpinyl acetate, and geranyl acetate) were used to evaluate the phytochemical relationship between the Iranian T. kotschyanus essential oils by AHC analysis using the XLSTAT software, version 2018.1.1.62926 (Addinsoft, Paris, France). Pearson correlation was used as a measure of similarity, and the unweighted pair-group method with arithmetic average (UPGMA) was used to define the clusters.
Statistical software SPSS version 16.0 (IBM, Chicago, IL, USA) was used to analyze the insecticidal experiments. The normality of mortality data was checked by Kolmogorov-Smirnov test. Probit analysis was used to calculate the Lethal Concentrations (LC) and their 95% confidence limits along with concentration-mortality regression lines details. The data heterogeneity was also evaluated by χ 2 test. Data obtained from all insecticidal bioassays were subjected to analyses of variance, and the comparison of means were performed by the least significant difference (LSD) test.

Conclusions
The use of plant-derived essential oils to manage stored-product insect pests has been highly considered due to their rapid degradation in the environment and lower impact on non-target organisms than conventional chemicals. According to this study, the cosmopolitan stored-grain insect pest R. dominica is susceptible to the essential oils isolated from aerial parts of wild-grown Thymus species, including T. eriocalyx, T. fallax, T. kotschyanus, and T. vulgaris. Although the activity of detoxifying enzymes αand β-esterase was augmented in the adults treated by the essential oils, along with significant fumigant toxicity, the terpene-rich essential oils showed promising sublethal effects from biochemical disturbances to antinutritional activities. The energy resources was affected by the essential oils: The protein content decreased by all essential oils, the glycogen content decreased by T. kotschyanus and T. vulgaris essential oils, and the lipid content significantly reduced by T. fallax, T. eriocalyx and T. vulgaris essential oils. The amylase and protease digestive enzyme activities were also affected by the essential oils except T. vulgaris. Furthermore, nutritional indices of insect pest especially RGR and RCR were affected by the essential oils so that they caused considerable FDI from 20.4% to 61.1%. It was also found that Thymus essential oils are complex mixture of terpenic compound that are probably responsible agents in the multiple modes of action of essential oils. However, further research is recommended to investigate the side effects of these essential oils on stored products and beneficial non-target insects and to evaluate the pesticidal effects of their compounds to other insect pests along with effects on human health. Furthermore, the introduction of durable formulations against environmental factors may be effective in their application.

Data Availability Statement:
The data that support the findings of this study are available upon reasonable request.