Optimizing Antioxidant Activity and Phytochemical Properties of Peppermint (Mentha piperita L.) by Integrative Application of Biofertilizer and Stress-Modulating Nanoparticles under Drought Stress Conditions

Drought stress (DS) negatively affects plant growth, productivity, and quality in semi-arid and arid regions. Nowadays, application of biofertilizers and stress-modulating nanoparticles (NPs) improves plant performance under stressful conditions. The study evaluated the impacts of arbuscular mycorrhizal fungi (Myco-Root) and TiO2 NPs on the nutrient uptake, dry yield, essential oil (EO) productivity, and EO quality of peppermint (Mentha piperita L.) under different irrigation regimes. The treatments included three irrigation regimes containing irrigation after 20% (I20, well-watered), 40% (I40, mild DS), and 60% (I60, severe DS) maximum allowable depletion (MAD) percentage of the soil’s available water as well as four fertilizer sources contain no fertilization (control), Myco-Root biofertilizer, TiO2 NPs, and an integrative application of Myco-Root + TiO2 NPs. The results demonstrated that the highest (195.72 g m−2) and the lowest dry yield (78.76 g m−2) of peppermint was obtained in well-watered conditions with integrative application of Myco-Root + TiO2 NPs and severe drought stress (I60) without fertilization, respectively. The dry yield of peppermint was reduced by 27.7 and 53.4% in mild (I40) and severe drought stress (I60), respectively. The maximum EO content (1.49%) and EO yield (2.30 g m−2) was recorded in mild drought stress (I40) treated with Myco-Root + TiO2 NPs. Based on the GC-MS and GC-FID analysis, 29 constituents were identified in peppermint EO, with the major constituents being menthol (38.99–52%), menthone (12.72–20.13%), 1,8-cineole (6.55–7.84%), and neo-menthol (3.14–4.52%), respectively. The maximum content of menthol, 1,8-cineole, and neo-menthol was obtained under mild drought stress (I40) fertilized with Myco-Root + TiO2 NPs. The results indicate that the integrative application of Myco-Root + TiO2 NPs could be used as an alternative method of using chemical fertilizers in sustainable agricultural systems for improving the EO quantity and quality of peppermint grown under drought stress conditions.


Introduction
Worldwide demand for the cultivation and usage of medicinal and aromatic plants (MAPs) has increased due to their pharmaceutical and therapeutic potentials [1]. Peppermint (Mentha piperita L.), as one of the MAPs belonging to the Lamiaceae family, is a natural hybrid derived from Mentha aquatica L. and Mentha spicata L. [2]. The global demand for peppermint increases day by day due to medicinal and flavoring properties and the use of this plant in cosmetic, pharmaceutical, and perfumery industries. Peppermint, as a traditional and folkloric medicine, is used for the treatment of colic, common cold, headache, sinusitis, sore throats, muscle aches, bronchitis, joint pain, gallstones, tuberculosis, and skin similation by more activation of Rubisco activase [27][28][29][30]. Earlier studies reported the positive impact of TiO 2 NPs on the growth, yield performance, and especially, productivity of biochemical components of MAPs. For example, Ahmad et al. [31] noted that the foliar spraying of TiO 2 NPs on peppermint increased the dry weight, EO content, EO yield, and main constituents such as menthol, menthone, and menthyl-acetate. Golami et al. [32] noted that foliar spraying of TiO 2 NPs affected the EO productivity and quality of rosemary (Rosmarinus officinalis L.) through increasing the major EO constituents such as α-pinene, camphene, limonene, 2-bornanone, and bornyl acetate.
Due to decreasing the efficiency of chemical fertilizers in drought stress conditions as well as the negative impacts of the synthetic fertilizers on the bioactive compounds of MAPs, the application of sustainable and eco-friendly fertilizers and new anti-stress nanoparticles is necessary. Among different methods of TiO 2 NPs application (seed coating, root application, and foliar spraying) as stress-modulating nanoparticles in plants, the foliar spraying of these nanoparticles is the most effective method for improving plant performance under stressful conditions [27]. It is worth noting that in the foliar spraying, a lower concentration of these nanoparticles is needed. In contrast, in the root application of TiO 2 NPs, higher concentrations are needed to improve plant growth, and the risk of its toxicity for plants in this condition is greater than foliar spraying [27]. The literature is scant about the integrative application of bio-fertilizers along with stress-modulating NPs. We hypothesis that i) the integrative application of Myco-Root + TiO 2 NPs enhanced the nutrient concentration, ii) integrative application of Myco-Root + TiO 2 NPs promotes photosynthesis pigments and productivity, and iii) integrative application of Myco-Root + TiO 2 NPs improves essential oil quantity and quality of peppermint. Therefore, the study evaluated the separate and integrative application effects of Myco-Root (AMF) as a bio-fertilizer and TiO 2 as stress-modulating NPs on the morphological, physiological, and phytochemical characteristics of peppermint under drought stress conditions.

Results
The analysis of variance (ANOVA) results showed that the AMF colonization, nutrient content, dry yield, essential oil (EO) content, EO yield, chlorophylls (a, b, and total), relative water content (RWC), proline, antioxidant enzymes (SOD, APX, and GPX) activity, total phenolic and flavonoid content, and net income of peppermint was impacted significantly by interaction of irrigation regimes × fertilizers application. In addition, the plant height, canopy diameter, and carotenoid content were significantly impacted by main factor effects, and the interaction of irrigation regimes × fertilizers application did not significantly impact the three mentioned traits. Additionally, the malondialdehyde (MDA) content was significantly affected by irrigation regimes (Tables 1 and 2). Table 1. The ANOVA (p-value) results of root colonization, nutrient content, agronomic traits, dry yield, essential oil content and yield, water use efficiency, and net income of peppermint affected by experimental factors.  Y: year; I: irrigation regimes; F: fertilizer sources. ns, * and ** indicated no significant difference, significant at 5% probability level, and significant at 1% probability level, respectively.

AMF Colonization
The highest root colonization (81.83%) was achieved in well-watered conditions (I 20 ) treated with Myco-Root, followed by integrative application of Myco-Root + TiO 2 NPs. The AMF inoculation rate decreased by 19.7 and 45.4% in mild (I 40 ) and severe drought stress (I 60 ), respectively ( Figure 1). Figure 1. The effects of different irrigation regimes and fertilizer sources on the root colonization of peppermint. I 20 , I 40 , and I 60 indicating well-watered, mild, and severe drought stress, respectively. Different letters indicate significant differences at the 5% level according to LSD's test.

Nutrient Content
The highest content of N, P, and K in peppermint was obtained in I 20 (well-watered conditions) fertilized with Myco-Root + TiO 2 NPs. Furthermore, the lowest content of the three mentioned elements was observed in I 60 (severe DS) without fertilization. On average, the content of N, P, and K decreased by 12.2, 13.9, and 12% in I 40 and decreased by 25.6, 22.3, and 24.2% in I 60 , respectively, when compared with well-watered conditions. Additionally, the content of the three mentioned nutrients was enhanced by 16, 19.1, and 15.81% after application of Myco-Root + TiO 2 NPs (Table 3).

Plant Height
Among different irrigation levels, the highest (51.49 cm) and the lowest (35.70 cm) plant height of peppermint was measured in I 20 and I 60 , respectively. The height of peppermint was reduced by 17.7 and 30.7% in I 40 and I 60 (Figure 2a). In addition, the application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs enhanced the height of peppermint by 11.4, 11.6, and 14.4%, respectively, when compared with control ( Figure 2b).

Figure 2.
The effects of different irrigation regimes (a) and fertilizer sources (b) on the plant height of peppermint. I 20 , I 40 , and I 60 indicating well-watered, mild, and severe drought stress, respectively. Different letters indicate significant differences at the 5% level according to LSD's test.

Canopy Diameter
Application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs enhanced the canopy diameter of peppermint by 10.6, 13.8, and 19.2%, respectively, when compared with control ( Figure 3b). Among different irrigation levels, the highest (41.36 cm) and the lowest (28.81 cm) canopy diameter of peppermint was measured in I 20 and I 60 , respectively. The canopy diameter of peppermint was reduced by 13.5 and 30.3% in I 40 and I 60 ( Figure 3a).  40 , and I 60 indicating well-watered, mild, and severe drought stress, respectively. Different letters indicate significant differences at the 5% level according to LSD's test.

Dry Yield
Integrative application of Myco-Root + TiO 2 NPs under well-watered conditions produced the highest dry yield of peppermint (195.72 g m −2 ), while the lowest dry yield (78.76 g m −2 ) was obtained in severe drought stress without fertilization. The peppermint dry yield was reduced by 27.7 and 53.4% in mild (I 40 ) and severe DS (I 60 ). Application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs increased the dry yield of peppermint by 8.7, 12.1, and 20.1%, respectively, when compared with control (Table 4).

Essential Oil Content and Yield
The maximum EO content (1.49%) and yield (2.30 g m −2 ) of peppermint was recorded in mild DS (I 40 ) treated by Myco-Root + TiO 2 NPs. The EO productivity under mild DS (I 40 ) was 22.3% higher than non-stress conditions. Also, integrative application of Myco-Root + TiO 2 NPs increased the EO content and EO yield of peppermint by 11.7 and 35.1%, respectively (Table 4).

Water Use Efficiency (WUE)
The highest WUE (177.77 g m −3 ) in peppermint was recorded in mild drought stress (I 40 ) fertilized with Myco-Root + TiO 2 NPs. Under mild and severe DS, the WUE of peppermint decreased by 3.6 and 6.8% in comparison with well-watered conditions. Additionally, application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs was enhanced WUE by 9.4, 11.6, and 18.5%, respectively (Table 4).

Chlorophylls Content
Integrative application of Myco-Root + TiO 2 NPs under I 20 irrigation regimes produced the maximum content of Ch a , Ch b , and Ch total . In contrast, the lowest content of the three mentioned traits was obtained in severe DS (I 60 ) without fertilization. The content of Ch a , Ch b , and Ch total decreased by 39, 22.2, and 35.2% in mild DS (I 40 ) and decreased by 77.6, 68.3, and 75.4% in severe DS (I 60 ), respectively. Moreover, integrative application of Myco-Root + TiO 2 NPs enhanced three photosynthesis pigments by 39.2, 55.9, and 43.3%, respectively ( Table 6). Table 6. Chlorophyll a, b, total, and relative water content of peppermint in different irrigation regimes and fertilizer sources.

Treatments
Chlorophyll a (mg g −1 Fresh Weight)

Carotenoid Content
Among different irrigation levels, the maximum and minimum carotenoid concentration of peppermint was measured in I 40 and I 60 , respectively. The carotenoid concentration was enhanced by 27% in I 40

Relative Water Content (RWC)
The RWC of peppermint was the highest (90.12%) in well-watered conditions (I 20 ) fertilized with Myco-Root + TiO 2 NPs. Furthermore, the lowest RWC content (61.13%) was obtained in severe DS (I 60 ) without fertilization. Under mild and severe DS conditions, the RWC content decreased by 10.6 and 21.1% in comparison with well-watered conditions. Additionally, application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs enhanced RWC content by 12.3, 8.5, and 13.7%, respectively ( Table 6).

Proline
The maximum proline concentration (2.07 µmol g −1 fresh weight) was obtained under severe drought stress conditions following application of Myco-Root + TiO 2 NPs. The proline concentration was enhanced by 124.6 and 188.5% in mild and severe drought stress conditions, respectively. Moreover, application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs enhanced proline concentration by 21.4, 22.3, and 40.8%, respectively (Table 7). Table 7. Antioxidant enzymes activity, proline, and total phenolic and flavonoid content of peppermint leaves in different irrigation regimes and fertilizer sources.

Malondialdehyde (MDA)
Among different irrigation regimes, the highest and lowest MDA content in peppermint leaves was observed under I 60 and I 20 irrigation regimes. The MDA content was enhanced by 16.5 and 25.8% in mild and severe drought stress conditions, respectively ( Figure 5). Different letters indicate significant differences at the 5% level according to LSD's test.

Antioxidant Enzymes Activity
The highest SOD, APX, and GPX activity was observed in mild drought stress (I 40 ) following application of Myco-Root + TiO 2 NPs. In comparison with non-stress conditions, the activity of APX, GPX, and SOD was enhanced by 52, 115, 124% under mild DS and 30, 66, and 42 under severe DS, respectively. Interestingly, integrative application of Myco-Root + TiO 2 NPs enhanced the three enzyme activities by 32, 43, and 31% when compared with control (Table 7).

Phenolic and Flavonoid Content
The maximum phenolic (34.79 mg g −1 ) and flavonoid (17.24 mg g −1 ) content was obtained in mild drought stress (I 40 ) treated with Myco-Root + TiO 2 NPs. The lowest phenolic (25.26 mg g −1 ) and flavonoid (9.77 mg g −1 ) content was also observed in nonstress conditions (I 20 ) without fertilization. In comparison with non-stress conditions, the phenolic and flavonoid content of peppermint was enhanced by 20.3, 34.4% in mild stress and enhanced by 9.4, 22.5% in severe stress, respectively. Additionally, application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs increased the phenolic content by 10.5, 9.4, and 19.2% and increased the flavonoid content by 23, 16.1, and 28.2% when compared with control, respectively (Table 7).

Net Income
Mean comparisons showed that the net income of peppermint among different irrigation regimes and fertilizer sources was higher in well-watered conditions (I 20 ) with application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs and also mild drought stress (I 40 ) fertilized with integrative application of Myco-Root + TiO 2 NPs. Under mild and severe DS conditions, the net income decreased by 19.1 and 59.5% in comparison with well-watered conditions. Additionally, application of Myco-Root, TiO 2 NPs, and integrative application of Myco-Root + TiO 2 NPs enhanced net income by 15.3, 19, and 28.5%, respectively, when compared with control ( Figure 6).

Discussion
The results of the study demonstrated the AMF inoculation rate was reduced in DS conditions. The decreasing soil moisture negatively affects spore germination and spore density, which decrease the microbial population and efficiency in soils [24]. Similarly, Ghanbarzadeh et al. [25] noted that the AMF root colonization decreased under severe DS due to reducing the availability of optimum soil moisture and insufficient supply photosynthetic carbon by the host plant for growth and spore germination of mycorrhiza fungi.
The nutrient content of peppermint decreased under DS conditions may be due to the reduction in nutrient diffusion and mass flow of nutrients under water limitation conditions [33]. In addition, the nutrient availability in the rhizosphere area depends on the activity of the soil's microbial population. In water deficit conditions, the decreasing soil microbial activity reduced the solubility and availability of nutrients [34]. Similarly, Salehi et al. [35] showed that the concentration of nitrogen, phosphorus, and potassium of chamomile (Matricaria chamomilla L.) decreased under DS condition. On the other hand, the integrative application of Myco-root + TiO 2 NPs increased the nutrient content of peppermint. The AMF-host plant symbiosis improves nutrient content in plants through enhancing surface absorbing capability as a result of the extensive hyphal network and acidification of the rhizosphere area by the release of H + , which increases nutrient solubility, especially bioavailable P-minerals and other micronutrients [19]. Similarly, Begum et al. [36] observed that the concentration of nitrogen, phosphorus, and potassium of tobacco (Nicotiana tabacum L.) increased by 48.8, 89.1, and 111.2%, respectively, with AMF inoculation. Additionally, it has been reported that the usage of TiO 2 NPs enhanced the content of N in plants by increasing enzyme activities involved in N metabolisms and also improving transformation of inorganic N to organic N [37,38].
Our results showed that the agronomic traits (plant height and canopy diameter) and dry yield of peppermint reduced in mild and severe DS conditions. The decrease in plants productivity under DS conditions was attributed to the negative effects of water scarcity on the cell division and elongation, reducing stomata conductivity and CO 2 uptake as a result of closing the stomata, and the increase of lipid peroxidation of membranes which decrease the photosynthesis rate and productivity of plants [39]. Similarly, Govahi et al. [40] indicated that the sage productivity decreased under DS conditions due to chloroplast destruction and decreaseing photosynthesis efficiency. On the other hand, the integrative application of Myco-Root + TiO 2 NPs enhanced the peppermint's agronomic traits such as plant height and canopy diameter as well as plant productivity under well-watered and DS conditions. Since there was a significant positive correlation between dry yield of peppermint and nutrient content (N, P, and K) and chlorophylls content (Figure 7), the higher productivity of peppermint under co-application of Myco-Root + TiO 2 NPs was attributed to the increasing nutrient accessibility, which leads to increasing leaf area, chlorophyll formation, and photosynthesis capacity [30,41]. Moreover, it seems that the increasing plant height through integrative application of Myco-Root + TiO 2 enhanced the photosynthesis rate by improving the absorption of sunlight, which led to increasing the agronomic traits and the plant's productivity. Results of previous studies showed that shoot dry matter of coriander (Coriandrum sativum L.) increased by 92.3% with AMF inoculation [42]. Another study on thyme (Thymus vulgaris L.) showed that the shoot dry matter yield was enhanced by 56.1%, by foliar application of TiO 2 NPs [43].
The EO content and main EO constituents of peppermint such as menthol, 1,8-cineole, and neo-menthol were enhanced in mild drought stress. In DS conditions, the closing stomata reduced absorption of CO 2 and the photosynthesis rate of plants. By closing stomata, the NADPH+H + concentration increases significantly in plant cells in comparison with NADP + , which can act as an inhibitor of photosynthesis. In this way, the biosynthesis of EO compositions conditions through usage the NADPH+H + increases plant efficiency under DS conditions [11]. Similarly, Abbaszadeh et al. [12] noted that EO yield and composition of Rosmarinus officinalis L., such as α-pinene and camphene of rosemary, increased by 31.6, 17.8, and 38.2% in mild water stress. Furthermore, integrative application of Myco-Root + TiO 2 NPs improved EO quali-quantitative of peppermint under DS conditions. The integrative application of the mentioned fertilizer can improve the EO content and its quality in two different ways: (i) Direct effects: in this case, the EO biosynthesis was enhanced due to increasing the accessibility of nutrients which play an important role in EO precursor compositions as well as intermediate compounds such as acetyl coenzyme A, NADPH, and ATP [44]. (ii) Indirect effects: in this case, the increasing photosynthesis rate as a result of higher accessibility of nutrients and modulating the negative effects of drought stress enhanced carbohydrates' productivity for the growth of cells and production of EO secreting glands, and in this way, it will be able to increase EO content in the plant [45]. Similarly, Karagiannidis et al. [46] noted that inoculation with AMF increased the EO content in mint (Mentha viridis L.) and oregano (Origanum onites L.) plants. It has been reported that the EO content, yield, and main constituents of sage such as cis-thujone, camphor, and 1,8 cineole increased by 75, 174.5, 87.5, 30.8, and 123.1%, respectively, after foliar application of TiO 2 [47].
Additionally, the EO yield of peppermint was enhanced in mild drought stress (I 40 ) treated by Myco-Root + TiO 2 NPs. The EO yield of peppermint is calculated from the dry yield and the EO productivity and has a direct relationship with the two mentioned factors (Figure 7). Therefore, the increasing EO yield of peppermint could be explained by the role of the integrative application of Myco-Root + TiO 2 NPs in increasing dry yield (by increasing a plant's yield components such as canopy diameter, etc.) and EO productivity under DS conditions.
The content of Ch a , Ch b , and Ch total of peppermint reduced under mild and severe DS due to the decomposition of chlorophylls as a result of the higher activity of ROS and lipid peroxidation under stressful conditions [48]. Jahani et al. [49] concluded that the total chlorophyll concentration of peppermint reduced by 9.1 and 31.4% in moderate and severe DS conditions due to increased photo oxidation of chlorophyll and prevention of chlorophyll biosynthesis as a result of chloroplast breakdown. Bhusal et al. [50] reported that the decreasing chlorophyll content under drought stress conditions could be due to stomata closing and also inhibiting the biosynthesis of chloroplast proteins. In contrast, the carotenoid concentration was enhanced in mild drought stress. Oxidative damage of drought stress in plant tissues is modulated by the activity of enzymatic and non-enzymatic antioxidant systems. The non-enzymatic antioxidants include carotenoids (xanthophylls and carotenes), alpha-tocopherol, and ascorbates, which improve plant tolerance in the face of stressful conditions by decreasing the activity of ROS compounds [51]. The integrative application of Myco-Root + TiO 2 NPs enhanced the chlorophylls and carotenoid concentration under DS conditions. The chlorophylls, found in the chloroplast of photosynthesizing plants, consist of two parts, including porphyrin head and a long hydrocarbon with a phytol sequence [52]. Porphyrin is composed of four nitrogen-containing pyrrole rings. The complement of the chlorophyll molecule is a Mg ion that forms a chelate with four nitrogen atoms in the center of the ring. Since the chlorophyll and carotenoid concentration have a significant and positive correlation with nutrient content (Figure 7), the enhancement of chlorophylls content could be attributed to the higher nutrient accessibility by integrative application of Myco-Root + TiO 2 NPs.
Relative water content (RWC) of plants refers to the state of water content in plant cells, especially under DS conditions [53]. The reduction in RWC under DS conditions was due to decreasing the absorption rate of water by roots and also increasing transpiration rate by leaves. Interestingly, separate application of Myco-root and integrative application with TiO 2 NPs enhanced the RWC content of peppermint under drought stress conditions. The symbiotic association between AMF and host plant roots improves the hydraulic conductivity of water into the roots, which will ultimately increase the relative water content in plant cells. Penetration of fungal hyphae into the root cortex and endoderm provides a path of low resistance across the root for the movement of water, and as a result, water will encounter less resistance across the root until it reaches the xylem [23]. Similarly, Gholinezhad et al. [54] determined that the RWC content of sesame (Sesamum indicum L.) decreased by 23 and 25% in moderate and severe DS conditions, while inoculation of the plant with two AMF species (Funneliformis mosseae and Rhizophagus irregularis) increased the RWC content by 35 and 26%, respectively.
The content of proline and MDA was enhanced in mild and severe DS conditions due to the role of ROS compounds in increasing lipid peroxidation of the cell membrane and releasing various aldehydes, such as MDA [55]. In these conditions, the increasing osmotolerant metabolites, such as proline, decrease osmotic potential and maintain turgescence pressure of plant cells, which reduces the negative effects of ROS compounds and lipid peroxidation and enhances stabilization of the membranes [56]. Similarly, Javanmard et al. [57] noted that the proline content of Lallemantia iberica was enhanced under moderate and severe DS conditions. Additionally, integrative application of Myco-Root + TiO 2 NPs enhanced the proline content, especially under drought stress conditions. Since proline is one of the amino acid compounds and nitrogen is one of the essential elements in the structure of amino acids, the higher proline concentration by integrative application of Myco-Root + TiO 2 NPs could be attributed to the availability of N through symbiotic association between AMF-host plants and also the role of TiO 2 NPs in the conversion of inorganic N to organic N [58,59].
The antioxidant enzymes' activities increased sharply under DS conditions. The higher activity of SOD, APX, and GPX enzymes under drought stress conditions play an important role in scavenging ROS compounds and decreasing lipid peroxidation, which leads to stabilization of the cell structure and an improvement in plant growth [60]. Additionally, the integrative application of Myco-Root + TiO 2 NPs enhanced antioxidant enzymes' activity. Since enzymes are protein compounds, the higher activity of enzymes by application of Myco-Root + TiO 2 NPs was due to the accessibility of nutrients, especially nitrogen [57]. Similarly, foliar spraying of TiO 2 NPs increased the antioxidant enzymes' activities, such as CAT and APX in Moldavian balm and peppermint [29,61]. Furthermore, Begum et al. [62] showed that inoculation with AMF increased the antioxidant enzymes' activities (SOD, CAT, POD, and APX) in Maize under drought stress conditions. Phenolic and flavonoid compounds play a crucial role in the formation of various biomolecules, protecting plants against stresses. Plants produce the ROS when exposed to the stress condition, and it seems that the biosynthesis of antioxidants, flavonoids, and secondary metabolites plays an important role in protecting plant cells for detoxifying ROS and also improving protein and amino acid stabilization [63]. It is worth noting that integrative application of Myco-Root + TiO 2 NPs, especially under mild and severe DS, improved phenolic and flavonoid content of the peppermint plant, which was attributed to enhancing nutrient availability that enhances the activity of enzymes involved in the biosynthesis of phenolic and flavonoid compounds [64]. Similarly, Shoarian et al. [65] noted that the application of TiO 2 NPs enhanced the phenolic and flavonoid compounds of Lallemantia iberica under drought stress conditions. In addition, Rashidi et al. [66] reported that the content of phenolic compounds in flowers of Ipomoea purpurea L. enhanced by 50%, 55.8%, and 71%, respectively, after colonization with Funneliformis mosseae, Rhizoglomus fasciculatum, and Rhizoglomus intraradices.
Integrative application of Myco-Root + TiO 2 NPs under mild DS enhanced the WUE of peppermint. These results indicated that the integrative application of the mentioned fertilizers enhanced the crop productivity per unit of water used as a result of improving nutrient availability, chlorophyll formation, and photosynthesis rate as well as decreasing the negative impacts of DS on the plant cells. In addition, the integrative application of Myco-Root + TiO 2 NPs enhanced net income and economic value of peppermint plants through increasing the dry yield, EO content, and EO yield under DS conditions. Amiri et al. [67] noted that inoculation with AMF under moderate water stress on rose geranium (Pelargonium graveolens L.) increased water use efficiency. Our previous study results showed that integrative application of AMF + TiO 2 increased water use efficiency of sage by 34.7% [13].

Study Site
A field study was conducted from 2019 to 2020 at the research farm of Maragheh University, Iran. The soil physico-chemical properties (depth of 0-30 cm) and climatic data of the experimental area are presented in Tables 8 and 9, respectively. Table 8. Physico-chemical properties of field soil (average of two years).

Experiment Design and Crop Management
The experiment was conducted with a two-factor split-plot design with three replications. The main factor was three irrigation regimes which included 20% (I 20 ), 40% (I 40 ), and 60% (I 60 ) maximum allowable depletion (MAD) percentage of the soil's available water (SAW) as normal irrigation, mild DS, and severe DS, respectively. Moreover, the sub-factor was the application of fertilizer sources including no-fertilizer as a control, TiO 2 NPs (100 mg L −1 ), arbuscular mycorrhizal fungi (Myco-Root) inoculation, and integrative application of Myco-Root and TiO 2 (Myco-Root + TiO 2 NPs). Each plot contains 5 rows of peppermint with 3 m length. Further, the distance of between rows was set to be 40 cm. Peppermint seedlings were sown by hand on the 6th and 10th day of April in 2019 and 2020, respectively. The first irrigation was performed immediately after sowing. In Myco-Root treatment, 70 g of soil including 1000 spores (Funneliformis mosseae) per 10 g soil obtained from Zist Fanavar Pishtaz Varian Company, Karaj, Iran (with the commercial name Myco-Root) was added to each hole at planting.

Implementation of Irrigation Regimes
The drought stress regimes were applied one month after sowing. The soil moisture was continuously measured with a TDR (Model German, FM-Trime) probe at depth of 0 to 30 cm. To determination of the irrigation depth used the following equations [11]: (1) in which θ fc and θ pwp are the soil field capacity (27.1%) and soil permanent wilting point (13.7%), respectively. D is the soil layer depth (cm), I d is the depth of net irrigation, P is the fraction of SAW (20%, 40%, and 60%) which must be depleted from the root zone, Ig is a gross irrigation depth, and Ea is the efficiency of irrigation that was considered as 65% [40]. The average of the final volume water used in irrigation treatments of I 20 , I 40 , and I 60 in both years was calculated as 1160, 870, and 580 L (=1.16, 0.87, and 0.58 m −3 ), respectively.

Synthesis of TiO 2 Nanoparticle
Titanium dioxide was synthesized based on the previous method reported by Gohari et al. [29]. Firstly, an amount of desirable titanium isopropoxide was stirred and hydrolyzed under cold conditions (0 • C) to form the TiO(OH) 2 white precipitate. The TiO(OH) 2 precipitate was washed with distilled water. Subsequently, TiO(OH) 2 was dissolved in nitric acid to produce a clear and homogeneous titanyl nitrate (TiO(NO 3 ) 2 ) solution. The TiO(NO 3 ) 2 was saved in a 250 mL beaker with a urea solution (molar ratio of 1:1). The obtained solution was placed into the muffle furnace for 2 h at 400 • C to produce solid titanium dioxide, and then it was kept inside the vacuum oven until usage. Before foliar spraying, TiO 2 nanoparticles were dissolved in distilled water and then mixed with a mechanical stirrer for 2 h at 35 • C.
Transmission electron microscopy (TEM) (Figure 8) and FT-IR spectrum analysis (Figure 9) was performed by a Zeiss EM-90 instrument at 80 kV tension and Perkin Elme-Spectrum Two FT-IR Spectrometer instrument, respectively. In addition, the chemical composition of the produced TiO 2 NPs was determined using the FT-IR spectrum (Figure 9a,b). A broad peak was seen between 460 and 870 cm −1 , which was consistent with Ti-O-Ti stretching vibration bonds [68] (Figure 9a). The broad bands at 2700-3600 and 1625 cm −1 were identified as the result of stretching and bending vibrations of O-H groups of adsorbed water on the surface of TiO 2 , which are crucial for photocatalytic activity, respectively. The Ti-O modes were associated with the peak at 1473 cm −1 . The stretching vibration of TiOO-H is represented by the double absorption peak at 2288 and 2324 cm −1 [69,70] (Figure 9b).  In addition, the XRD and SEM analyses were conducted at the Central Laboratory of the University of Tabriz, Iran using Tongda-TD-3700 and VWGA3 TESCAN instruments, respectively (Figures 10 and 11). X-ray diffraction (XRD) was used to confirm that TiO 2 nanoparticles were crystalline and had a certain shape ( Figure 10). The characteristic peaks of TiO 2 in the anatase phase were connected to all of the diffraction peaks in the XRD spectrum (JCPDS file   (101) [69,71,72]. The size (D hkl ) of the crystal domains was determined by using the widths at half maximum height (B) in the XRD peak's broadening. This was achieved with the Debye-Scherrer equation: where λ is the X-ray beam's wavelength (1.5406 Å), θ is the Bragg angle, and K is a constant attributed to the shape of the domain, that is about 0.9. The average size of crystallite in the TiO 2 nanoparticles that were made was estimated to be about 18 nm. The SEM was used to study the morphology of the synthesized NPs, which are depicted in Figure 11a,b for TiO 2 at two different magnifications. In the SEM image of TiO 2 , some agglomeration was also seen, which may be caused by calcination NPs. Synthesized TiO 2 NPs were discovered to be spherical in shape.

Agronomic Traits, Dry Yield, and Water User Efficiency
Before harvesting, the agronomic traits of peppermint were measured in ten plants of each plot randomly. In both years, the aerial parts of peppermint were harvested randomly at the full flowering stage from a 1.4 m −2 area. In order to determine the dry yield of peppermint, the harvesting samples were kept at room temperature without direct sunlight for two weeks. Additionally, the water use efficiency (WUE) was calculated as [73]: WUE (g m −3 ) = Dry yield/Final volume water used in irrigation treatments (5)

Essential Oil Analysis Essential Oil Extraction
The peppermint EO was extracted from 40 gr of dried aerial parts of peppermint mixed with 400 mL that hydro-distillated for 3 h by a Clevenger-type apparatus. After extracting EO, the samples were dried by adding sodium sulfate and kept at 4 • C for chemical analysis. The EO yield of peppermint was calculated based on the following equation [74]: Essential oil (%) = [EO extracted by hydro distillation/40 gr dried aerial parts] × 100 (6) EO yield (g m −2 ) = EO content × dry yield of aerial part (7) GC/MS and GC-FID Analysis The peppermint EO compounds were identified using a gas chromatograph instrument (Agilent 7990 B, USA) equipped to a 5988A mass spectrometer with HP-5MS capillary column (length 30 m, 5% phenylmethyl polysiloxane, 0.25 mm inner diameter, and 0.25 µm film thickness). Briefly, the first temperature of the oven reached 60 • C in 5 min, and was subsequently at a rate of 3 • C/min reached in 240 • C and maintained for 20 min. The carrier gas was helium, which was set up with a flow rate of 1 mL/min. The injector and detector temperatures were set up at 230 and 240 • C, respectively. The injector was adjusted based on the split ratio of 1:30. Mass scan range and ionization temperatures were 40-400 m/z and 220 • C, respectively. The mass spectra were recorded using an electron impact ionization technique (ionization energy of 70 eV). In order to compute the linear retention indices (RIs), a mixture of n-alkane series was co-injected with the essential oil into the GC instrument based on the above temperature conditions [5]. For identification of volatile ingredients, their retention times compared with those retention indices that were reported in reference libraries [75] or through comparison of their mass spectra with data in the NIST 05 database. Volatile components were separated using an Agilent 7990B gas chromatography instrument (USA) with a capillary column VF-5MS (0.50 µm f.t., 30 mL, 0.25 mm i.d., 5% phenyl methylpolysiloxane) and flame ionization detector (FID). The oven temperature program was the same as the above analytical conditions. The injection volume of essential oil was 1 µL that was diluted in an n-hexane solution (1:100). The peak area normalization method was also used for calculating the percentage compositions of essential oil [74].

Nutrient Content
The nitrogen (N) content of peppermint leaves was measured based on the Kjeldahl method [76]. In order to determine the phosphorus (P) and potassium (K) concentration of peppermint, 1 gr of dried leaves of peppermint was put into the porcelain crucible and placed in a muffle furnace. The temperature of the muffle furnace reached 550 • C gradually and was maintained for 8 h until a white ash formed. After cooling the muffle furnace, the samples were dissolved in 5 mL hydrochloric acid (HCL) 2N. After 15 min, the samples filtered through filter paper (Whatman No. 42), were brought to the volume (50 mL), and stood for 30 min. The concentration of K and P was determined by flame photometry (Jenway PFP7C) and the yellow color method (colorimetric method by Ammonium Vanadate/Molybdate) using a spectrophotometer at 470 nm absorption spectra, respectively [77].

Relative Water Content (RWC)
For the determination of RWC in peppermint leaves, firstly, the fresh weight of 10 fully developed leaves was measured (LFW). Then, the samples were kept in distilled water at 4 • C and measured leaf turgor weight after 24 h (LTW). After that, the leaves dried in the oven at 70 • C for 48 h and the leaf dry weight (LDW) was weighed [78]. Finally, RWC content was calculated by the following formula:

Chlorophyll and Carotenoid
In order to quantify the content of chlorophyll a (Ch a ), b (Ch b ) and carotenoid, 0.5 gr of peppermint fresh leaves was digested in liquid nitrogen and mixed with 10 mL of 80% acetone. Samples centrifuged at 10,000 rpm and supernatant after 10 min were transferred to a new tube, and the absorbance was read at 645, 663, and 470 nm spectrophotometrically (UV-1800, Shimadzu, Japan). After that, the content of photosynthetic pigments was calculated by following equations [79]: For determination of Malondialdehyde (MDA) concentration, firstly, 0.1 g of peppermint fresh leaves sample was extracted using 1 mL 0.1% TCA (trichloroacetic acid) and centrifuged at 12,000 rpm. After 15 min, the supernatant was mixed with 4 mL of reaction mixture including 20% TCA and 0.67% 2-thiobarbituric acid (TBA). After that, the sample was heated in a water bath at 95 • C for 15 min and then quickly cooled down in an ice bath for 10 min and centrifuged at 10,000 rpm for 5 min at 4 • C. Finally, the absorbance was recorded at 532 and 600 nm spectrophotometrically (UV-1800, Shimadzu, Japan), and the MDA concentration was represented as nmol g −1 FW [80].

Proline
Firstly, we took a 0.5 g sample of fresh peppermint leaves after being digested in liquid nitrogen, mixed with 10 mL 3% sulfosalicylic acid and then centrifuged at 12,000 rpm. After 10 min, 2 mL of the separated supernatant and the samples were put into a tube and mixed with 2 mL of acid-ninhydrin and 2 mL of glacial acetic acid. The obtained mixtures were heated in a water bath for 1h at 100 • C and then placed on an ice bath. Then, 4 mL of toluene was added to the samples and vortexed for 20 s. Finally, absorbance was read at 520 nm spectrophotometrically, and the proline content was represented as µmol g −1 of fresh weight [81].

Antioxidant Enzyme Activity Guaiacol Peroxidases (GPX) Activity
In order to prepare the reaction mixture, firstly, we mixed a hydrogen peroxide (70 mM), phosphate buffer (50 mM) and 10 mM of guaiacol and then added 20 µL of enzyme extract to the reaction mixture. The absorbance was read at 420 nm for 180 s and represented as µmol min −1 mg −1 protein [82].

Ascorbate Peroxidase (APX) Activity
The reaction mixture for APX was prepared with a mixture of 0.1 mM H 2 O 2 , potassium phosphate buffer (50 mM, pH = 7), 0.5 mM ascorbic acid, and 0.1 mM EDTA. Then, 150 µL of enzyme extract was added to the reaction mixture. After that, the absorbance was read at 290 nm (coefficient of absorbance = 2.8 mmol −1 cm −1 ) [83].

Superoxide Dismutase (SOD) Activity
The reaction mixture for the determination of SOD included L-methionine (12 mM), potassium phosphate buffer (50 mM, pH = 7.8), sodium carbonate (50 mM), 2 µM of riboflavin, 0.1 mM EDTA, and 75µM of nitro blue tetrazolium (NBT). A 100 µL amount of enzyme extract was added to the reaction mixture and kept in 25 • C under illumination. After 15 min, absorbance was read at 560 nm and represented as µmol min −1 mg −1 protein [84].

Total Phenolic and Flavonoid Content
The total phenolic and flavonoid content of peppermint was determined based on the previous method, which was reported by Singleton and Rosi [85] and Nagy and Grancai [86].

Determination of Root Colonization
To specify the AMF colonization rate, firstly, the collected root samples were washed with water for removing extra particles and then cut into pieces of 1 cm length. In order to eliminate cytoplasmic contents from cells, which facilitates staining, root pieces were cleared in hot 10% KOH for 10 min. After rinsing three times with distilled water, the root sample was acidified in 2% HCL for 15 min at room temperature and stained in a hot 20% lactophenol cotton blue solution for 15 min [87,88]. After coloring, the AMF colonization percentage was determined using the grid line intersection method, which was described by Giovannetti and Mosse [89].

Statistical Analysis
The normality and homoscedasticity of data regarding the morphological and physiological traits were verified using the Kolmogorov-Smirnov and Levene test, respectively. Combined variance analysis of data was carried out with SAS v9.4 (SAS Institute, Cary, NC, USA) and SPSS v25 software. The mean comparison between treatments was analyzed by the LSD test at probability levels of 0.05%. Irrigation levels and fertilizer sources were considered as fixed factors, and block and years were random factors. In addition, for drawing the Pearson correlation matrix (correlation plot) between main studied traits such as dry yield, EO content, EO yield, photosynthetic pigments, nitrogen, phosphorus, potassium, and major EO compounds (menthol, menthone, and 1,8-cineole), R software v3.2.4 was used.

Conclusions
The results of the study showed that the integrative application of Myco-Root + TiO 2 has great potential for improving plant performance under drought stress conditions by enhancement of the nutrient uptake, chlorophyll formation, and relative water content as well as antioxidant activity. Importantly, improving essential oil quantity and quality of peppermint demonstrated another positive impact of the integrative application of Myco-Root + TiO 2 under stressful conditions. Generally, the results of the study indicated that the integrative application of Myco-Root + TiO 2 NPs could be suggested as a sustainable strategy for improving EO quantity and quality of peppermint under DS conditions. Data Availability Statement: The datasets generated and analyzed during the current study are available from the corresponding author upon reasonable request.