Funneliformis mosseae inoculation under water deficit stress improves the yield and phytochemical characteristics of thyme in intercropping with soybean

Intercropping of medicinal plants/legumes along with bio-fertilizer application is a relatively new sustainable practice for improving the yield and secondary metabolites production. Here, a 2-years field experiment was performed to evaluate the effects of water deficit stress and arbuscular mycorrhizal fungi (AMF) application (as bio-fertilizer) on nutrients concentration, dry matter yield, essential oil quantity and quality of thyme in intercropping with soybean. Three irrigation levels, including (i) irrigation after depletion of 20% (I20) as non-stressed, 50% (I50) as moderate water deficit and 80% (I80) available water as severe water deficit were applied as the main factor. The sub-factor was represented by different cropping patterns including thyme sole culture, replacement intercrop ratio of 50:50 and 66:34 (soybean: thyme) and the third factor was non-usage (control) and usage of AMF. According to our results, the thyme dry yield under moderate and severe water deficit stress decreased by 35 and 44% in the first year, and by 27 and 40% in the second year compared with non-stressed (I20) plants, respectively. Also, the macro- and micro-nutrients of thyme leaves increased significantly in intercropping patterns after application of AMF. The maximum essential oil percentage of thyme was achieved in 50:50 intercropping ratio treated with AMF. Under moderate and severe water deficits, the major constituents of thyme essential oil including thymol, p-cymene and γ-terpinene were increased in intercropping patterns treated with AMF. Generally, AMF application in intercropping ratio of 50:50 may be proposed to farmers as an eco-friendly approach to achieve desirable essential oil quality and quantity in thyme under water deficit stress conditions.

www.nature.com/scientificreports/ oil is also characterized by the presence of ɣ-terpinene and p-cymene which are recognized as the biogenetic precursors of thymol and carvacrol 6 . The climate-changes around the world indicate clearly that annual maximum temperature and the aridity are increasing 7 . Water deficit stress, known as one of the limiting abiotic factors, affects every aspect of the plant growth, morphology, physiology and productivity in the arid and semi-arid regions 8,9 . In medicinal and aromatic plants, increasing the essential oil content is considered as one of main defense mechanisms against drought stress conditions, because in case of stress, more metabolites are produced in the plants to prevent cell oxidization 10 .
In addition, the nutrient use efficiency (NUE) of chemical fertilizer is reduced under drought conditions due to the decrease of mobility and low rate of mineral diffusion. Thus, the lower NUE along with the negative effects given by the chemical fertilizer application (e.g., soil, air and water pollution, soil acidification, mineral depletion, decrease of soil micro-organisms and harmful effects on human health) forced the agrochemical companies to replace the chemical inputs with eco-friendly fertilizers. Recently, the interaction among plants and soil microbial population is a conventional method to limit the demand of synthetic fertilizers over sustainable agricultural systems. Arbuscular mycorrhizal fungi (AMF) are soil borne micro-organisms forming a mutualistic symbiotic association with more than 80% of plant species 11,12 . The coexistence between AMF and host plants improves the conditions for an optimal plant growth. In medicinal and aromatic plants, previous research indicated that inoculation of AMF in dill (Anethum graveolens L.), rose geranium (Pelargonium graveolens L.) and Mentha × piperita improves the nutrient uptakes and essential oil productivity 13,14 .
Conventional agricultural systems are based on cultivation and productivity using sole plant cultures with intensive application of chemical inputs including synthetic fertilizer, pesticides, herbicides and others 15 . The expansion of intensive crop production systems for long terms caused several problems, such as the reduction of biodiversity and the environmental contamination. The increasing concern regarding the negative impacts of these conventional agricultural methods on the ecosystems suggests the exploration of alternative methods endowed with improved efficiency and sustainability as well as high productivity 16 .
The intercropping system is one of the well-known sustainable and eco-friendly agricultural methods in which two or more plants are grown simultaneously. This provides interconnections between plants leads to more efficient use of environmental resources including nutrients, water, radiation, land and also reduced damages of pests and plant diseases 17,18 . Among the various intercropping systems available, intercropping based on legume species provides an added value in order to minimizing the adverse impact of chemical fertilizer (due to nitrogen transfer from legume to companion crops or N as an input for subsequent crops) 19,20 .
Iran is basically an arid country, with an average annual rainfall ranging from 224 to 275 mm/year (70% less than global average) 21 . Recently, a considerable demand for thyme and its different products has been registered around the world. In order to comply with this objective, the improvement of thyme productivity based on the application of sustainable strategies that increase nutrient uptake under water scarcity conditions, especially in the arid and semi-arid regions of the world such as Iran, is urgently needed. Thus, an experiment was aimed to compare the impacts of different water deficit levels and application of AMF on the dry matter yield, nutrients uptake, essential oil quantity and quality of thyme in different intercropping patterns with soybean.

Materials and methods
Area of study. The experiments were performed over two successful growing year (April-October 2018 and April-October 2019) at the Maragheh University research farm, East Azerbaijan Province, Iran. The filed soil had a sandy clay loam texture with pH of 7.53; the electrical conductivity (EC) and the TOC (total organic-carbon) content was 1.216 dS/m and 0.76%, respectively. The contents of macro-and micro-elements and physicochemical characteristics are listed in Table 1. Also, the monthly weather data in the research area are provided in Table 2.

Treatments.
A split plot experiment, based on a RCBD design with three replications and 18 treatments, was performed. The main factor was represented by different water deficit stress levels, including (i) irrigation after depletion of 20% available water (I 20 ) as non-stressed, (ii) 50% of available water (I 50 ) as moderate water deficit, and (iii) 80% available water (I 80 ) as severe water deficit. The sub-factor was different cropping patterns, including (i) thyme sole culture (T s ) and (ii) soybean-thyme replacement intercropping series (in proportion of 50:50 and 66:34). The experimental treatments based on proportions of soybean-thyme intercrops pattern is given in Table 3. In addition, the third factor was represented by (i) non-application (control) and (ii) application of AMF as bio-fertilizer. Thyme and soybean were sown in 4 m long rows and the distance of each row was set to 0.4 m (Fig. 1). Also, a 1.5 m distance was maintained between different treatments in order to exclude any influence of lateral water movement. For AMF treatments, the mycorrhizal fungal spores (Funneliformis mosseae; 1000 spore/10 g soil) were added to the soil at planting (Zist Fanavar Sabz Company, Iran). Based on the Table1. Physico-chemical properties of field soil (depth of 0-30 cm) an average over the 2 years. In these equations, SAW, θ fc , θ pwp , d, I d , p, E a and I g are: soil available water, soil field capacity (27.4%), permanent wilting point (13.9%), soil layer depth (cm), irrigation depth (cm), fraction of SAW (20%, 50% and 80%) that can be depleted from the root zone, irrigation efficiency (%) averagely assumed 65% 23 and the gross depth of irrigation (cm), respectively.
Measurements. Dry yield. At full flowering stage, the thyme aerial parts were harvested randomly on 15 August 2018 and 27 August 2019. The plants were harvested from a 1.6 m 2 area (middle line of each treatment). After harvesting, the samples were dried at room temperature for two weeks and the dry matter yield was recorded.

Mycorrhizal fungal colonization.
To determine the root colonization percentage, the young fresh roots of thyme were collected from the depth of 10-20 cm of the soil. The root colonization was investigated according to the method described by Phillips and Hayman 24 . Firstly, the roots were washed with water and then cut into the small pieces (1.5 cm length) and followed by in 10% potassium hydroxide at 90 °C for 10 min. Afterwards, the segments were placed in 2% HCl at 25 °C for 10 min and then stained with trypan blue (0.05%) and lactic acid (80%) for 12 h. Eventually, the roots were washed and stored in lactic acid, glycerol and distilled water (lactic acid and glycerol in proportions of 1:1:1 v/v/v) solution. Finally, the colonization percentage of thyme root was measured according to the method of Giovannetti and Mosse 25 .
Essential oil content and chemical analyses. 40 g of dried aerial parts of herbs were roughly grounded and placed in a 1 L flask where 400 mL of distilled water and conventional hydro-distillation for 3 h, using a Clevenger-type apparatus, was performed. In order to remove possible water drops, Na 2 SO 4 was added to each extracted essential oil and kept at 4 °C in the refrigerator before analysis. For calculating the essential oil content and yield of thyme, the following equations were used: Essential oil content of thyme = (distilled essential oil (g)/40 g) × 100 Essential oil yield of thyme (g m −2 ) = dry yield of thyme (g m −2 ) × essential oil content (%). www.nature.com/scientificreports/ The thyme essential oil constituents were analyzed using GC-FID and GC-MS instruments. GC-MS analysis was performed in an Agilent gas chromatograph (7990B) coupled to a 5977A mass spectrometer. The separation was done in a HP-5 MS capillary column (5% Phenyl Methylpolysiloxane, 30 m length, 0.25 mm i.d., 0.25 µm film thickness). The following oven temperature was set up for the column: 10 min at 60 °C, then the temperature raised with the rate of 3 °C/min to 250 °C, finally held for 15 min at 250 °C. The carrier gas (Helium) was used at a flow rate of 1 mL/min. Electron-impact (EI) was set to be 70 eV. The injector split ratio was 1:50 and mass range acquisition was from 40 to 400 m/z. Gas chromatographic (GC-FID) separations was performed in an Agilent technology instrument (Agilent 7990B, USA), coupled with a flame ionization detector (FID). VF-5MS column used in the GC-FID device had the same stationary phase and dimensions as the HP-5MS one. The oven temperature program was the same used for GC-MS analysis. Injector and detector temperatures were set at 230 and 240 °C, respectively. Diluted essential oils (1: 100 in n-hexane) were used for the injection. The peak areas were normalized without using correction indices.
To identify essential oil constituents, a combinational procedure was performed including determination of arithmetic retention indices using coherence of homologous series of hydrocarbons (Supelco, Bellefonte, USA), comparing retention indices with those reported in the reference literature 26 , and the computer matching of mass acquisition with the WILEY275 and NIST 05 libraries.
Plant nutrient concentration. Thyme leaf samples were analyzed for nutrients concentrations. The Kjeldahl method was used to determine the N content. The concentration of K was also measured using flame photometry 27 . P was determined by the yellow method 28 . The content of P was determined at 470 nm using a spectrophotometer. An atomic absorption spectrophotometer (AA-6300 F; Shimadzu, Japan) was used to determine the concentration of micro-nutrients including Ca, Fe and Mg of each samples.

Statistical analysis. The homoscedasticity and normality test of all data were checked with Levene test and
Anderson-Darling methods, respectively. Analysis of variance (ANOVA) and the means comparison (compared with LSD at 5% probability level) were performed using the SAS software version 9.1.

Results
Root colonization percentage. In both years, the highest root colonization percentage was measured in the 50:50 intercrop ratio with non-stressed plants (I 20 ). Also, the lowest percentage of thyme root colonization (46.81% in the first year and 56.87% in the second year) was observed in the sole culture of thyme under severe stress (I 80 ). The colonization rate in moderate (I 50 ) and severe water deficit (I 80 ) decreased by 5 and 16% in the first year, and 2 and 10% in the second year, respectively, in comparison with I 20 (Fig. 2).  www.nature.com/scientificreports/ Dry yield. The dry matter yield of thyme was significantly impacted by single effect of three factors and interaction of cropping patterns × water deficit levels. The highest and lowest dry yield (in both years) was achieved in sole culture of thyme with well-watered condition (I 20 ) and intercrop ratio of 66:34 with I 80 water deficit (Fig. 3).
In the first and second year, application of AMF enhanced the dry productivity by 17.3 and 14.5%, respectively, in comparison with control. Furthermore, the thyme dry yield in the moderate (I 50 ) and severe water deficit (I 80 ) decreased by 35 and 44% in the first year, and by 27 and 40% in the second year, respectively, in comparison with I 20 (Table 4).

Nutrient concentrations. Most of macro-and micronutrient concentrations was significantly impacted
by water deficit levels, cropping patterns and AMF application. In both years, the maximum and minimum concentration of N, K, Ca and Fe was obtained in intercrop ratio of 50:50 and thyme sole culture. In addition, the highest content of Mg was recorded in intercrop ratio of 66:34. Base on the average of both years, the concentration of N, P, K, Ca, Mg and Fe with AMF application enhanced by 10, 31, 14, 8, 17 and 16, respectively, in comparison with control. Among water deficit levels, the highest and lowest concentration of micro and macronutrients was achieved in non-stressed plants and severe water deficit conditions (Table 5).
Essential oil content. AMF application, different cropping patterns and water deficit levels and interaction of cropping patterns × AMF had a significant impact on the essential oil content of thyme. Among three levels of water deficit, the maximum essential oil productivity was measured in I 50 water scarcity that was 52% (in 2018) and 22% (in 2019) higher than those observed in I 20 , respectively (Table 4). Based on the interaction of cropping patterns × AMF application, the highest essential oil content was achieved in intercropping ratio of 50:50 after AMF application. Also, the lowest essential oil content was observed in the sole culture of thyme without usage of AMF (Fig. 4).
Essential oil yield. AMF application, different cropping patterns and water deficit levels and interaction of cropping patterns × water deficit had a significant impact on the essential oil yield of thyme. In the first and second year, application of AMF increased the essential oil yield by 59.6 and 23%, respectively, in comparison with   www.nature.com/scientificreports/ control (non-AMF). Based on the interaction of cropping patterns × water deficit, the highest essential oil yield was achieved in the non-stressed plants (I 20 ) in sole culture (Fig. 5). Finally, the essential oil yield in I 50 and I 80 water deficit stresses decreased by 3 and 35% in the first year, and by 9 and 31% in the second year, respectively, when compared with well-watered plants (I 20 ) (Table 4).  www.nature.com/scientificreports/ Essential oil composition. In total, 29 compounds were identified in the thyme essential oils, accounting for 97.45-99.45% and 97.10-99.77% of the total compositions at the first and second year, respectively (Tables 6  and 7). In the first year, the main essential oil constituents were thymol (31.12-37.26%), γ-terpinene (13.40-17.48%), p-cymene (11.74-13.92%), camphene (6.52-8.51%), 1,8-cineole (2.94-5.83%), camphor (3.74-5.61%), (E)caryophyllene (2.57-4.19%), carvacrol (2.01-3.41%) and β-myrcene (1.72-2.31%). In the second year, the main constituents of thyme essential oil were thymol (40.89-48.59%), γ-terpinene (16.93-21.55%), p-cymene (13.07-16.26%) and carvacrol (2.95-4.06%). The highest content of thymol and γ-terpinene was measured in 50:50 intercrop ratio with moderate and severe water deficit under application of AMF. However, the highest content of p-cymene was observed in 66:34 intercropping with severe water deficit (I 80 ) and application of AMF. In addition, the lowest content of the three mentioned constituents was achieved in the sole culture of thyme and non-stressed plants without AMF. Also, classification of essential oil constituents based on the chemical classes showed that oxygenated monoterpenes and monoterpene hydrocarbons represented the major fraction in the essential oils of thyme. In both years, the highest content of the two mentioned groups was obtained in 50:50 intercrop treated with AMF under moderate and severe water stress (Tables 6 and 7).

Discussion
The results demonstrated that the root colonization rate increased in intercropping patterns in comparison with sole culture due to an improvement of the soil microbial communities through an enhancement of the organic matter and accessibility of legume-fixed nitrogen into the soil which have a positive effect on root colonization 29 . De Araujo Pereira et al. 30 showed that the presence of legume species in intercropping patterns could stimulate the alkaline phosphatase activity and increase the AMF root colonization. However, the AMF root colonization of thyme decreased under higher water deficit levels which could be due to the reduction of soil microbial activities under drought conditions. Alike, Shukla et al. 31       moisture, changed the AMF root colonization. These authors noted that the maximum and minimum AMF colonization in the vigna mungo (Vigna mungo (L.) Hepper) and wheat (Triticum aestivum L.) was observed in the field capacity (FC) and half-field capacity (FC/2), respectively. Our results showed that the macro-and micro-nutrients of thyme increased in intercropping patterns. The difference in rooting depth and root expansion of plants affects the competition of intercropping components over nutrients. Therefore, the higher nutrient content in intercropping systems could be due to the increase of the environmental use efficiency as a result of different root distribution (e.g., thyme is a shallow-rooted and soybean is a deep-rooted type crop) 32,33 . Duchene et al. 18 showed that the nutrients availability in intercropping patterns improved due to the higher symbiotic nitrogen fixation by legumes, root exudation of enzymes such as phosphatases, carboxylates, and decrease of soil acidity by production of H + in comparison with sole culture of plants.
Also, the nutrient concentrations of N, P, K, Ca, Mg and Fe improved after application of AMF. The extensive network of AMF hyphae in the plant roots increases the absorption rate and area leading to the enhancement of the micro-and macro-elements uptake, especially lowly mobile nutrients such as P 11 . Alike, Rezaei-Chiyaneh et al. 34 reported that the nutrients uptake, especially N and P, in intercropping of Isabgol (Plantago ovata)/lentil (Lens culinaris) significantly increased with AMF inoculation. Also, Weisany et al. 14 showed that application of AMF (Funneliformis mosseae) significantly enhanced the concentration of P, K, Fe, Na, Zn and Mn in the dill/ common bean intercrop.
Furthermore, the results of this study exhibited that nutrient concentration decreased by raising the water deficit levels. This could be explained by the reduction of nutrient supply through mineralization, decreasing mass flow and nutrient diffusion that affect the kinetics of nutrient uptake per unit of roots 35 . In addition, Bista et al. 36 reported that the activity or concentration of major uptake-proteins for nutrients decreased under water stress and had a negative impact on the nutrient's uptake from roots. Also, these authors noted that the nutrients Table 7. Proportion of thyme essential oil constituents in different cropping patterns and irrigation levels with AMF application (second year). a RI, linear retention indices on DB-5 MS column, experimentally determined using homologue series of n-alkanes0 b Relative retention indices taken from Adams. c Ts, 50:50 and 66:34 indicating thyme sole culture and intercropping ratio of thyme with soybean; I 20 , I 50 and I 80 indicating nonstressed, moderate and severe water deficit levels; M 0 : non-usage of AMF, M 1 : usage of AMF. d Identification methods: MS, by comparison of the mass spectrum with those of the computer mass libraries Wiley, Adams and NIST 08; RI, by comparison of retention index with those reported in Adams and NIST 08; Std, by comparison of the retention time and mass spectrum of available authentic standard. Bold values show the major constituents of thyme essential oil. www.nature.com/scientificreports/ concentration is higher in the upper soil layers and drought stress reduced the nutrients absorption by decreasing the water uptake in the upper soil layers. The results showed that dry matter yield of thyme significantly increased under sole culture. In sole culture, the entire of each plot was assigned with one of the plant species. Therefore, the decrease of plants yield in intercrops could be explained by a lower number of partial plants density in these cropping patterns in comparison with sole culture. It seems that the interspecific competition in the thyme/soybean intercrops was higher than the intraspecific competition in sole culture 37 . In addition, the shading of the companion plant (soybean) over thyme seedlings in intercropping patterns could decrease the photosynthetic rate, leading to a decrease in thyme productivity 32 .
The thyme dry matter yield significantly decreased with the increase of water deficit levels. Similar to other plants, thyme plants that tolerate water or drought stresses have apparently evolved mechanisms such as early flowering and structural traits including the reduction of plant leaf area and height and also metabolic responses, such as photosynthetic alterations and proline accumulations (data not shown), to survive under water stress conditions 22 . Therefore, the observed decrease in the thyme growth parameters could be explained by the lower availability of sufficient moisture around the root zone and lower absorption of nutrients as a result of the lower proliferation of root biomass. These conditions had a negative impact on the chlorophyll pigments, the photosynthetic efficiency and finally plant yield 23 . Alike, Askary et al. 4 showed that total dry matter yield in two different species of thyme (T. daenensis Celak. and T. vulgaris) decreased by 21.6 and 38.9% in the 67 and 33% FC (field capacity) in comparison with 100% FC, respectively.
Also, our results demonstrated that application of AMF significantly increased the biomass yield of thyme. The mycorrhizal rhizosphere pH is lower due to the absorption of ammonium ions and releasing of H + ions. Reducing soil acidity and also exudation of organic acids and phosphatase enzymes increase the solubility and accessibility of nutrients to plant roots 11 . Therefore, an improvement of plant growth parameters and productivity of thyme after AMF application could be due to the enhancement of the macro-and micronutrients uptake such as N, P, K, Ca, Fe and Mg through the expansion of the mycelium and the development of the root system. Moreover, a higher production of plant growth hormones such as gibberellin (affecting the longitudinal cellular growth and primarily stem internodes), auxin and cytokinin (affecting the cell division) has been observed with AMF inoculation which has positive role in improving plant performance 11,38 . Golubkina et al. 39 showed that the AMF inoculation increased the biomass and dry yield of Artemisia dracunculus L., Hyssopus officinalis L. and Lavandula angustifolia Mill. compared with untreated plants.
Essential oils are mixtures of low molecular weight compounds that are responsible for the characteristic aroma of medicinal and aromatic plants. The ratio of one special compound or group of compounds to another affects the quality of essential oil in these herbs 3 . Although secondary metabolites such as those contained in medicinal and aromatic plant essential oils are conventionally controlled by the plant genotypes, their biosynthesis is strongly affected by the environmental factors 9 . The obtained results showed that the essential oil content and main constituents (namely thymol, γ-terpinene and p-cymene) increased in moderate and severe water deficit conditions when compared with non-stressed plants. The increase of essential oil productivity under moderate and severe water stress could be explained by a decrease in the leaf area of plants leading to a raise of the glandular trichomes density and essential oil accumulation per unit of leaf tissue. Also, CO 2 fixation rate in water stress condition reduced due to decreasing uptake of CO 2 by closing stomata. In these condition, the ratio of NADPH + H + /NADP + enhanced as a result of generating massive oversupply of NADPH + H + . It can be concluded that biosynthesis of secondary metabolites such as essential oil, alkaloids, phenols and etc. using additional NADPH + H + known as one of the plant defense mechanism for improving performance in medicinal and aromatic plants 10 . Turtola et al. 40 noted that the increase of essential oil productivity in medicinal and aromatic plants under drought stress conditions may be related to a low allocation of carbon to the growth, representing a trade-off between growth and plants defense. Govahi et al. 23 concluded that the maximum essential oil content of sage (Salvia officinalis L.) was achieved in moderate drought stress. Also, these authors noted that the essential oil productivity in the moderate and severe drought stress increased by 109 and 84%, compared with non-stressed plants, respectively.
The results showed that the essential oil content and main constituents increased in intercropping patterns with AMF Application. The nutrients accessibility, especially N, plays an important role in the essential oil gland cells size and essential oil productivity in medicinal and aromatic plants 41 . Rostaei et al. 42 noted that the accessibility of nutrients, especially N and P, play a key role in the development and division of the glandular trichomes, essential oil channels and secretory ducts. In addition, the nutrients availability improves the plant performance and photosynthetic rate in plants. Calsamiglia et al. 43 noted that availability of glucose in plant cells, which is produced in the process of photosynthesis, enhanced terpenes constituents (especially monoterpenes like thymol) and essential oil productivity in medicinal and aromatic plants. Therefore, the increment of essential oil content and constituents in intercropping patterns with AMF application may be related to the higher accessibility of nutrients by expansion of mycorrhiza hyphae as well as atmospheric nitrogen fixation by legumes and directly/indirectly transfer to component plant (thyme) due to the superior complementarity (special, chemical, and temporal) in intercropping patterns compared with pure cultures 18 . Likewise, Amani Machiani et al. 44 concluded that the essential oil content and constituents of Foeniculum vulgare Mill. and Dracocephalum moldavica L. improved in different intercropping patterns with Phaseolus vulgaris L. after organic fertilizer usage (humic acid). Rezaei-chiyaneh et al. 45 showed that different intercropping ratios of fennel/common bean enhance the essential oil content by 24% in comparison with fennel sole culture due to availability of nutrients by nitrogen fixation of legume species.
Interestingly, the thymol content increased by 29.3% in the second year, representing that the thyme essential oil quality was significantly improved in the second year compared with the first year of growing. The quantity and quality of essential oil in medicinal and aromatic plants depends on many factors such as the plant species, www.nature.com/scientificreports/ climatic conditions, harvest year and plant age. In perennial medicinal plants such as thyme, the first year of planting is considered as the establishment year. In the second year, which is usually considered as the year of economic yield, the plant performance will increase by higher nutrient uptake and photosynthetic rate. Therefore, the improvement of essential oil quantity and quality in aromatic plants during the second year could be explained by the increase of nutrients availability that have positive rule in increasing the oil glands and terpene biosynthesis 46 . Alike, Mechergui et al. 47 reported that the content of thymol in oregano (Origanum vulgare L.) increased in the second year from 31.8 to 41.7% in Tunisian Nefza population and 18.4-31.5% in Tunisian Krib population.
The essential oil yield of thyme increased in the sole culture and decreased after raising the water deficit stress levels. The essential oil yield, calculated from the herbage yield multiply with essential oil content and had a positive relationship with the two mentioned factors. These differences in thyme essential oil yields could be due to the higher dry productivity in the sole culture and application of AMF and to the decrease of the dry productivity under water deficit conditions. In addition, the higher essential oil yield of thyme after application of AMF may be attributed to the higher dry matter and essential oil productivity. Likewise, Bahreininejad et al. 22 showed that the essential oil yield in T. daenensis decreased under moderate and severe water stress due to a decrease of herbage yield in these conditions.

Conclusions
Our study showed that the water deficit levels decreased dry matter yield of thyme. In contrast, the essential oil content of thyme increased under moderate and severe water deficit conditions. Also, the macro-and micronutrients concentration, essential oil content of thyme and main constituents (thymol, γ-terpinene and p-cymene) increased in intercrop ratio of 50:50 and 66:34 with AMF inoculation. Based on the obtained results, in the semi-arid and arid regions, in which water deficit has a negative impact on the plant productivity, intercropping of thyme with soybean under application of AMF as bio-fertilizer is recommended as a feasible and eco-friendly alternative strategy to improve the thyme essential oil quality and quantity compared with the sole culture.