Improvement of Growth, Yield, Seed Production and Phytochemical Properties of Satureja khuzistanica Jamzad by Foliar Application of Boron and Zinc

Satureja khuzistanica Jamzad is a valuable and endemic medicinal plant. Boron and zinc are essential elements for the vegetative and reproductive growth of plants and have significant effects on yield, essential oil composition and the seed production of plants. To investigate the effects of the foliar application of zinc and boron on the growth, yield, seed production and phytochemical properties of S. khuzistanica, a study was conducted in a factorial experiment with three replicates in two consecutive years based on a randomized complete block design. The foliar application of boron (B) at three concentrations (control or distilled water, 0.4% and 0.8% as H3BO3) and zinc (Zn) at three concentrations (control or distilled water, 0.3% and 0.6% as ZnSO4) was carried out. Our results showed that the foliar application of B resulted in a significant increase in the fresh and dry weights of plants, the dry weight of stems, drug yield, seed yield, seed germination and 1000-seed weight. At the same time, the application of B resulted in a significant decrease in seed emptiness. The fresh and dry weights of plants, drug yield, seed yield, 1000-seed weight and seed germination were also significantly improved by Zn foliar spraying compared to the control. Application of 0.8% B resulted in a significant decrease in seed emptiness by 14.16% and 22.37%, as compared to the control. The foliar spraying of B and Zn improved the total phenolic content, the essential oil content and the yield and antioxidant activity of S. khuzistanica. Moreover, B application generally concentrated more carvacrol in the essential oil (in the first experimental year). In contrast, no significant differences were observed between Zn treatments in carvacrol content and total flavonoids. The use of several microelements, such as B and Zn, could improve both the quantity and quality of S. khuzistanica. Additionally, improvement of seed set and seed quality by the foliar spraying of Zn and B may be useful for growing plants in arid and semi-arid areas.


Introduction
Satureja khuzistanica Jamzad, which belongs to the Lamiaceae family, is one of the most important medicinal plants endemic to Iran. Here, the plant grows wild in dry areas, sunny areas and the calcareous rocky soils of the southwestern region (Lorestan, Ilam and yield in Arachis hypogaea L. [23]. The essential oil composition of Litsea cubeba has also been significantly affected by the foliar application of Zn and B [24].
Nowadays, the herbal drugs of S. khuzistanica are collected from wild niches to meet the demands of the pharmaceutical and food industries. Due to the high demand of the pharmaceutical and food industries for S. khuzistanica raw materials and the challenges facing wild collections such as the degradation of genetic diversity, the extinction of ecotypes and the unreliable supply of herbaceous materials, there is an urgent need to cultivate the species in agroecosystems [3]. The main problem associated with the cultivation of the species is low seed production. The growth and essential oil composition of plants in agroecosystems may be affected by various complex factors including genotype, geographical and climatic conditions, soil properties, plant density and nutrient levels [3,25,26]. Regarding the influence of Zn and B on plant growth and development, the aim of the present study was to investigate the effects of the foliar application of Zn and B on the growth, yield, seed production, essential oil content and composition of S. khuzistanica.

Plant Materials and Field Experiment
Healthy cuttings of S. khuzistanica were taken from a mother plantation grown at the Medicinal and Aromatic Plants Research Garden of Lorestan University, Lorestan province, Iran (altitude 1170 m, longitude 48 • 26' E, latitude 33 • 44' N), in September 2018. The species was identified by the Department of Botany of the Research Institute of Forests and Rangelands (TARI), Tehran, Iran. A voucher specimen (herbarium No. 58416) has been deposited at the Herbarium of TARI. Cuttings were planted in plastic pots and grown in the greenhouses of the faculty of agriculture, Lorestan University.
A factorial experiment based on the Randomized Complete Block Design with three replications was conducted for two consecutive years (2018 and 2019). The results of farm soil analysis and climatic conditions of the field experiment for both cropping years are presented in Table 1 and Supplementary Table S1, respectively. Small plants of S. khuzistanica were hand-planted in the field on 4 April 2018. Each plot (2 × 3 m) contained 4 rows with 50 cm row distance and 50 cm plant distance. All the experimental plots received 50-ton cow manure, 120 kg P as triple superphosphate and 70 kg P as potassium sulfate in March 2018. In both years, 90 kg N was applied as urea in two equal splits, the first dose on 4 May and the second dose 2 months later, in July. During the field experiment, the plants were watered once a week, weeds were controlled by hand hoeing and no pesticide was used. Drip irrigations were provided to the plant at 7-day intervals. Foliar application of B at three concentrations (control or distilled water, 0.4% and 0.8% applied as H 3 BO 3 ) and Zn at three concentrations (control or distilled water, 0.3% and 0.6% applied as ZnSO 4 ) were carried out. The treatments were applied at three different times-one month after transplanting, three months after transplanting and at the floral budding stage.

Sampling and Measuring Traits
At full flowering stage, five plants were randomly selected from each plot, and yield attributes such as plant height, canopy diameter, number of main branches and number of sub-branches were measured. The aerial parts of selected plants, consisting of the leaves and upper portions of shoots and inflorescences, were harvested in each plot in early October in 2018 and 2019, then air-dried in the shade at room temperature. Plant dry weight was also measured for each plot. The oil-bearing parts of the plant (leaves and flowering shoots) were separated from the woody stems, weighed and recorded as the drug yield of S. khuzistanica [27]. The remaining plants in each plot were kept in the field until late November for seed set, then harvested to measure the seed yield. Seed emptiness and 1000-seed weight were also calculated.

Concentration of Zn and B in Leaves
Leaf samples were taken from each plot in early October 2018 and 2019. Samples were rinsed in demineralized water and then air-dried at 70 • C for three days in a forcedair oven [28]. Dried leaves were ground to pass a 30-mesh screen and digested in a ternary solution (HNO 3 :H 2 SO 4 :HCLO 4 = 10:1:4 with volume). The Zn and B content was determined according to an atomic absorption spectrophotometric method [29]. With the optimized spectroscopic conditions, a steady baseline was recorded. Standard and sample solutions of zinc and boron were aspirated. Spectrograms were recorded by measuring the absorbance at 213.9 nm and 249.8 nm, using boron and zinc lamps, respectively. An air-acetylene flame with fuel flow rate of 1600 mL/min was used. Zn and B concentrations were calculated using the following equation: element mg kg = titre value from machine × volume used molar mass

Seed Germination
For the measurement of the seed germination percentage, 75 hand-cleaned seeds were used for each treatment and divided into three replications of 25 seeds. First, seeds were disinfected with 5% sodium hypochlorite for 30 s, then washed with distilled water. Seeds were placed in Petri dishes with filter paper moistened by distilled water. The Petri dishes were kept in a germinator set for a temperature of 25 • C and photoperiod of 16 h. Constant monitoring ensured sufficient moisture throughout the experiment period. The number of seeds germinated as indicated by the radicle emergence was determined every two days. The germination experiment was continued until germinated number became constant. Finally, the seed germination percentage was calculated for each treatment [30].

Essential Oil Isolation
Essential oils from dried samples of S. khuzistanica were isolated by hydro-distillation for 3 h (20 g samples were mixed into 400 mL distilled water). The isolated oils were dried over anhydrous sodium sulfate and kept in amber glass vials at 4 • C before analysis. The essential oil content (w/w) of each sample was determined from three consecutive replicates and expressed as mg per 100 g sample [31].

GC-FID Analysis and Carvacrol Measurement
The GC-FID analysis of the oils was conducted using a Shimadzu (GC-17A, Kyoto, Japan) gas chromatograph coupled with a flame ionization detector (FID). A DB-5 fusedsilica capillary column (30 m × 0.25 mm i.d., film thickness 0.25 mm) was used to separate the essential oil compounds. The oven temperature was increased from 60 to 280 • C at a rate of 4 • C min −1 and held isothermal at 250 • C for 10 min. Ion source and transfer-line temperature were 250 • C. Ultra-pure helium was used as the carrier gas with a linear flow rate of 1.1 mL min −1 . Injector and detector (FID) temperatures were 285 • C. The split ratio was 1:60. All samples and standards were injected three times [31]. The identification of carvacrol was performed by co-injection of samples with authentic standards (Merck, Germany). Quantification was done by GC-FID peak area normalization considering an equal response factor [32].

Extraction Process
Extraction of plant extract was done by maceration method [33]. A total of 15 g of the dried samples was powdered and dissolved in 150 mL of aqueous methanol 80%. Extraction was done by soaking the mixture on a shaker at room temperature for 72 h. The mixture was filtered through filter paper (GF/A, 110 mm; Whatman, Maidstone, UK). The obtained extract was concentrated in a vacuum rotary evaporator at 50 • C, then dried in a vacuum oven at 450 • C for 3 h. The dried extract was used to measure total phenol, total flavonoids and the antioxidant activity assay [25].

Total Phenolic Content
To measure the total phenol content, the Folin Ciocalteu reagent was used [34]. First, 10 µL of each sample extract was mixed with 490 µL distilled water and 500 µL Folin-Ciocalteu reagent. After adding 500 µL saturated sodium carbonate (1%), the mixture was incubated for 2 h at room temperature. Finally, the absorbance of the samples was measured by a UV-vis spectrophotometer (UV-1800) at 765 nm. All samples were analyzed in three replications and the total phenol content was calculated by the standard curve of gallic acid at concentrations of 0, 50, 100, 150, 250 and 500 mg·L −1 . Total phenolic results were expressed as gallic acid equivalents (mg GAE·g −1 plant dry weight).

Total Flavonoid Content
Total flavonoid content was measured by aluminum chloride colorimetric method, according to Quetter-Deleu [35] but with a slight modification. Amounts of 1 mL of the standard extract or solution, 2 mL of 2% aluminum chloride hexahydrate methanolic solution and 6 mL of 5% potassium acetate solution were mixed. After 40 min incubation at 37 • C, the samples were read at 415 nm with a UV-vis spectrophotometer (UV-1800). The flavonoid content of the samples was calculated by plotting standard curve of rutin at concentrations of 0, 50, 100, 150, 250 and 500 mg·L −1 . The results were expressed in mg of rutin equivalents (RE)·g −1 plant dry weight. Samples and standards were analyzed in three replications. Antioxidant properties of methanolic extract of S. khuzistanica were evaluated by DPPH reagent (2,2-dipheny-l-picrylhydrazyl) according to Choi et al. [36]. Different concentrations of each extract (range from 0.01-1.0 mg·mL −1 ) were made in methanol, then 2.5 mL of the extract (methanol for the control) and 1.0 mL of a 3.0 × 10 −4 M DPPH solution in methanol were mixed and left for 30 min in the dark at 25 • C. The absorbance was recorded at 517 nm (As) using a UV-vis spectrophotometer (UV-1800). Subsequently, the absorbance of a blank solution containing 1.0 mL methanol and 2.5 mL plant extract solution was also recorded. The DPPH solution, plus methanol, were used as the control and ascorbic acid as a positive control. Finally, the antioxidant activity of the extract was expressed based on IC50, indicating a concentration of the extract that results in 50% free radical scavenging.

Ferric-Reducing Antioxidant Power (FRAP) Assay
The FRAP method is based on the reduction of Fe 3+ -TPTZ (yellow) to Fe 2+ -TPTZ (blue) at low pH values. First, fresh FRAP reagent was prepared by mixing 300 mm acetate buffer (pH 3.6), 10 mm TPTZ and 20 mm FeCl 3 ·6H 2 O in a ratio of 10:1:1 at 37 • C. Subsequently, 190 µL of FRAP solution was added to 10 µl of methanol extract and kept at 37 • C for eight minutes. The absorbance of the solutions was read at a wavelength of 593 nm using the Bio-Rad (Hercules, CA, USA) microtiter plate reader. A blank reading was taken, and aqueous solutions of ferrous sulphate heptahydrate (25-1000 µM) were used to generate a calibration curve. The antioxidant power of the extract was expressed as µM Fe (II). All tests were performed in triplicate [37].

Data Analysis
Analysis of variance of data was performed based on the studied experimental design with SAS statistical software. The Mean comparison of the treatments were done with the least significant difference experiment at the level of 0.05%.

Biomass and Yield Attributes
Boron (B) treatment had a significant effect on the fresh and dry weight of plants, the dry stem weight and the drug yield in both years. There was a significant difference in plant height, canopy diameter and the number of main and secondary branches between B treatments in the second year (Supplementary Tables S2 and S3). Foliar application of B significantly increased the fresh and dry weight of plants, the dry stem weight and the drug yield in both years, and the plant height, canopy diameter and the number of main and secondary branches in the second year. The highest fresh weight and dry weight of plants, dry stem weight and drug yield (40.4952 and 85.62 g·plant −1 in the first year and second year, respectively) of S. khuzistanica were obtained by applying 0.8% B in both years. Application of 0.8% B resulted in 33.32% and 31.94% increase in drug yield compared to control in the first and second year, respectively. Moreover, the highest plant height, canopy diameter, number of main branches and number of secondary branches were obtained by the highest B application in the second year of the experiment (Tables 2 and 3). Table 2. Mean comparison of the simple effects of boron and zinc foliar application on biomass and yield attributes of S. khuzistanica in the first year of the experiment.

Treatments
Number  The results of the analysis of variance (Supplementary Tables S2 and S3) showed that the effect of the Zn foliar application on the fresh and dry weight of plants and drug yield was significant in both years. The plant height and number of sub-branches were significantly affected by Zn application in the first year. Moreover, there was a significant difference between the various Zn treatments in the number of main branches and the dry weight of stems in the second year (Supplementary Tables S2 and S3). The highest values for the fresh and dry weight of plants and the drug yield were observed in both years with 0.6% Zn foliar spray. Zn foliar spraying increased the drug yield from the untreated values of 34.095 and 64.05 g·plant −1 to 37.808 and 92.61 g·plant −1 (10.89% and 44.59%) with the highest dose (0.6% Zn) in the first year and second year, respectively. Foliar fertilization with 0.6% Zn resulted in the highest plant heights, dry stem weights and the number of main branches in the second year (Tables 2 and 3).
The interaction of Zn * B was particularly significant for plant height in the first year of the experiment (Supplementary Table S2). The highest plant height (51.86 cm) was obtained in a simultaneous foliar application of 0.3% Zn and 0.8% B (Figure 1).

Seed Traits
The results from the analysis of variance showed significant effects of the B application on seed yield, 1000-seed weight, seed emptiness and seed germination in both years (Supplementary Table 4 , 5). Seed yield, seed germination and 1000-seed weight were significantly increased by the foliar application of B in both years. The highest values of seed yield, seed germination and 1000-seed weight were found at 0.8% B. Application of 0.8% B resulted in a significant decrease in seed emptiness by 14.16% and 22.37%, as compared to control in the first year and second year, respectively ( Table 4, 5).
The seed yield, 1000-seed weight and seed germination were significantly affected by the foliar application of Zn in both years (Supplementary Tables 4 and 5). Application of 0.6% Zn significantly promoted the seed yield (20.85% and 14.96%), 1000-seed weight (11.59% and 10.69%) and seed germination (13.29% and 11.79%), compared to the control in the first year and second year, respectively (Table 4, 5).

B and Zn Concentrations
The results from the analysis of variance showed that the foliar application of B had a positive effect on B content in the leaves of S. khuzistanica in both years (Supplementary Table 4, 5). The highest B concentration (19.19 and 19.96 mg.kg-1 in the first year and second year, respectively) was obtained by 0.8% B spray ( Table 4, 5). A significant difference was also observed between Zn treatments in terms of the Zn content of leaves (Supplementary Table 4, 5). The application of 0.6% Zn recorded the highest Zn content in the leaf (25.22 and 27.55 mg.kg-1 in the first and second years, respectively), while the lowest content was observed in the control in both experimental years (

Seed Traits
The results from the analysis of variance showed significant effects of the B application on seed yield, 1000-seed weight, seed emptiness and seed germination in both years (Supplementary Tables S4 and S5). Seed yield, seed germination and 1000-seed weight were significantly increased by the foliar application of B in both years. The highest values of seed yield, seed germination and 1000-seed weight were found at 0.8% B. Application of 0.8% B resulted in a significant decrease in seed emptiness by 14.16% and 22.37%, as compared to control in the first year and second year, respectively (Tables 4 and 5).
The seed yield, 1000-seed weight and seed germination were significantly affected by the foliar application of Zn in both years (Supplementary Tables S4 and S5). Application of 0.6% Zn significantly promoted the seed yield (20.85% and 14.96%), 1000-seed weight (11.59% and 10.69%) and seed germination (13.29% and 11.79%), compared to the control in the first year and second year, respectively (Tables 4 and 5).

B and Zn Concentrations
The results from the analysis of variance showed that the foliar application of B had a positive effect on B content in the leaves of S. khuzistanica in both years (Supplementary Tables S4 and S5). The highest B concentration (19.19 and 19.96 mg·kg −1 in the first year and second year, respectively) was obtained by 0.8% B spray (Tables 4 and 5). A significant difference was also observed between Zn treatments in terms of the Zn content of leaves (Supplementary Tables S4 and S5). The application of 0.6% Zn recorded the highest Zn content in the leaf (25.22 and 27.55 mg·kg −1 in the first and second years, respectively), while the lowest content was observed in the control in both experimental years (Tables 4 and 5). On the other hand, the interaction effect of the Zn and B application was not significant for the Zn and B leaf content in any of the experimental years (Supplementary Tables S4 and S5).

Essential Oil Content and Yield and Carvacrol Content
The results from the analysis of variance showed that the foliar application of B had a significant effect on the essential oil content and yield of S. khuzistanica in both years (Supplementary Tables S6 and S7). The foliar spray of 0.8% B increased the essential oil content in the leaves of S. khuzistanica by 18.6% and 21.58%, compared to the untreated plants in the first year and second year, respectively. Essential oil yield was more affected by the B application than essential oil content. It was 58.05% and 58.33% higher in the 0.8% B treatment than in the first year and second year control, respectively. The 0.8% B treatment resulted in a higher concentration of carvacrol in the essential oil in the first year, while the percentage of carvacrol as the main component of essential oil was significantly higher in 0.8% B plants than in untreated plants. In contrast, in the second year, no significant differences were observed among B treatments in terms of carvacrol content in essential oil (Tables 6 and 7).  Essential oil content and essential oil yield were significantly affected by zinc application in both years. However, the effect of Zn application on carvacrol content was not significant in any of the experimental years ( Supplementary Tables S6 and S7). The highest essential oil content (2.657% and 2.8461%) and yield (5.051 g·m −2 and 13.318 g·m −2 ) were obtained with the 0.6% Zn foliar application (Tables 6 and 7).

Total Phenol and Flavonoid Contents
The results from the analysis of variance showed that the total phenolic content of the extract of S. khuzistanica was significantly affected by B application in both years. Significant differences were observed in the total flavonoid content between the different B contents in the second year (Supplementary Tables S6 and S7). The total phenolic content was increased by 34.25% in the first year and 27.5% in the second year by the foliar application of 0.8% B compared to the control. Similarly, total flavonoid content was increased by 22.95% by the application of 0.8% B compared to the control in the second year (Tables 6 and 7).
The total phenolic content was significantly affected by the Zn application. No significant differences were observed between Zn treatments in the total flavonoid content and the flavonoid to phenolic ratio in both years and the total phenolic content in the second year (Supplementary Tables S6 and S7). The highest content of total phenolics (58.12 mg GAE·g −1 dry weight) was obtained under the 0.6% Zn spray in the first year (Table 6).

Antioxidant Capacity
The results indicated that the B application had a significant effect on the antioxidant activity of the S. khuzistanica extract calculated by DPPH and FRAP assays in both years (Supplementary Tables S6 and S7). The IC50 (50% reduction in DPPH concentration) value was declined by the foliar application of B in both years, indicating the highest free radicalscavenging activity of the extract in plants treated by 0.8% B. Similarly, the highest FRAP value (742.81 and 831.7 µmol Fe +2 ·g −1 dry weight) was observed for application of 0.8% B in the first year and second year, respectively (Tables 6 and 7).
Various Zn treatments showed significant differences in term of antioxidant capacity in both assays and both years (Supplementary Tables S6 and S7). Application of Zn significantly promoted antioxidant activity of the plant extracts in DPPH assay (decreasing IC50 by the Zn foliar application). The highest (IC50 = 0.0722 and 0.0688 mg·mL −1 in the first year and second year, respectively) and the lowest (IC50 = 0.0888 and 0.0788 mg·mL −1 in the first year and second year, respectively) antioxidant activity were observed by 0.6% Zn and the control treatment, respectively. Similarly, the highest antioxidant capacity in the FRAP assay (731.25 and 820.14 µmol Fe +2 ·g −1 dry weight in the first year and second year, respectively) was obtained by a 0.6% Zn application in both years (Tables 6 and 7). However, the interaction effect of Zn * B on the antioxidant activity of the extract was not significant in any of the applied assays (Supplementary Tables S6 and S7).

Biomass and Yield Attributes
Zinc plays an important role in the basic processes of cellular mechanism and carbon metabolism, which ultimately affects the growth and development of crops [38]. The results of the present study showed that Zn foliar application significantly increased the growth and yield of S. khuzistanica in both years of the experiment. This improvement in the growth and yield of S. khuzistanica is possibly due to the involvement of Zn in carbohydrate metabolism, indole acetic acid, RNA and ribosomal function. Our result is similar to the observations made by Wasaya et al. [39] in a study with maize, where they reported that Zn supply increased the growth and yield of crops by increasing the rate of photosynthesis, chlorophyll synthesis and Crop Growth Rate (CGR). Earlier studies demonstrated that Zn foliar application is effective in improving crop yield and quality, while its deficiency reduces crop yield and quality [40,41]. Zn is one of the most important factors affecting the activity of tryptophan synthase. Given the fact that tryptophan acts as a precursor for auxin synthesis, increasing auxin by Zn supply leads to the intensification of apex dominance and consequently increases the longitudinal growth of shoots [42]. In a study on Petunia hybrida, it was shown that the application of Zn resulted in increasing plant height, the number of leaves on a plant, the leaf area and leaf chlorophyll content. They concluded that a rise in plant growth and yield could be a result of the Zn-positive role in RNA metabolism and ribosomal content in plant cells, leading to the stimulation of carbohydrates, proteins and DNA synthesis. Zn also contributes to the synthesis of tryptophan as a growth booster [43].
Boron is of the most important essential elements for the higher plants and known to be involved in photosynthesis, N-fixation, meristematic tissue development, respiration, transpiration and other biochemical activities [10,44]. According to our results, the B foliar application increased most of the yield attribute traits of S. khuzistanica in both years. A similar result was found for sunflower, where B supply significantly improved plant growth and biomass [45], leading to the conclusion that sunflower growth components, leaf length, crop growth rate (CGR), relative growth rate (RGR), net assimilation rate (NAR) and chlorophyll content were increased by B foliar application. Previous findings have shown that there is a potential interaction between B and plant growth regulators, especially cytokinins. For instance, B deficiency decreased the amount of active cytokinin in Brassica napus seedlings [46]. Furthermore, B is one of the structural components of cell wall ingredients, i.e., rhamnogalacturonans, even though its deficiency disrupts cell wall structure and function and reduces the expression of several enzymes located in the cell wall involved in cell expansion, including pectin methylesterase, polygalacturonase and xyloglucan endotransglucosylase. Therefore, it stimulates cells division and elongation, especially in meristematic tissues, and increases plant growth [47].

Seed Traits
Zn acts in flower bud formation and flowering, playing roles in the synthesis of tryptophan as an auxin precursor, and in the translocation of carbohydrates to the site of bud development or to the bud itself [7]. The application of Zn in plants also increases fertility factors such as the formation of spore, fertility of pollen and secretion of stigma exudate, eventually stimulating seed production and fruit set [7,9]. A constant and continuous supply of Zn during plant development could improve the seed yield of the plant possibly due to better plant growth, Zn accumulation in the plant, the photosynthesis rate and enhanced production of phytohormones such as indole acetic acid (IAA), cytokinins and gibberellins [48]. In the present study, foliar application of Zn improved the 1000-seed weight, seed germination and seed yield. A similar effect of Zn application on the grain yield of wheat was also reported by Meena et al [49]. The results of a study by Farooq et al. [50] showed that the seed treatment of bean with Zn improved seed germination, which was attributed to the increase of α-amylase activity. Stimulating effects of the Zn supply on seed yield and germination have also been reported in other plants. For example, examining the impact of different Zn treatments on sunflowers revealed that the highest number of filled seeds and seed yield was obtained in the concentration of 0.5% [45]. Furthermore, it has been reported that the maximum yield and quality of okra seeds were obtained using Zn treatment [44].
According to our results, B had a significant positive effect on seed traits, increased seed yield, 1000-seed weight and seed germination and decreased seed emptiness. Improvement of seed yield and quality was a possible outcome of the improved reproductive growth of plants, carbohydrate metabolism, cell division, germination of pollen grains and the elongation of pollen tubes [7,11]. Azeem et al. [51] reported that the B foliar spray increased seed numbers and the seed weight of cotton through stimulating flowering, facilitating the formation of pollen tubes and increasing seed growth. Likewise, the B application induced the seed germination of Parthenium argentatum Gray [52]. B deficiency decreased esterase activity, which is known to be involved in a rupture of stigmatic cuticles for proper pollen tube penetration through stylar tubes for successful fertilization. It has been also observed that the stigmatic exudation that protects stigma from drying, resulting in proper pollen adhesion and germination, is reduced under B deficiency [53]. It could be concluded that the decrease in seed yield under B deficiency is likely to reduce the stigmatic receptivity at the time of anthesis, ultimately causing inhibition in pollen grain viability and resulting in poor germination and seed set/yield [53,54].

Phytochemical Evaluation and Antioxidant Activity
Although Zn application has been shown to affect plant physiology and biochemistry, debates and ambiguities remain about its effects on plant secondary metabolism [55]. In the present study, the effect of Zn application on carvacrol content as the main constituent of the essential oil and total flavonoid was not significant in any years of research. Interestingly, essential oil content and the yield and antioxidant activity of S. khuzistanica were significantly increased by the Zn application in both years. On the other hand, the total phenol was significantly risen due to the Zn application in the first year. It seems that Zn could protect the membrane proteins and lipids against free radicals. This may explain the increase in the antioxidant activity of the plant by the foliar application of zinc sulfate [56]. The foliar application of Zn increased phenol, flavonoids and the antioxidant activity of cowpea [57]. The essential oil yield of essential oil-bearing plants depends on the drug yield and essential oil content. Therefore, an increase in essential oil yield due to Zn application was caused by improved drug yield and essential oil content. The increased essential oil content and yield in S. khuzistanica as the result of Zn foliar spray are also in line with findings reported by Said-Al Ahl and Mahmoud [58], who found a positive and significant effect of Zn on essential oil content and the yield of Ocimum basilicum. Furthermore, Shahhoseini et al. [59] reported the positive effects of Zn on the biomass and terpenoid compounds of Tanacetum parthenium essential oil. The foliar application of Zn has improved the quantity and quality of the S. farinacea essential oil [19]. In general, Zn is involved in carbon assimilation, saccharide accumulation and scavenging of free radicals, and is overall imposed on photosynthetic carbon metabolism and translocation in relation to secondary metabolite biosynthesis [19,42]. Accordingly, the role of Zn in the carbon allocation to produce phenolic compounds in shikimic acid and acetate pathways can be attributed to this issue [42].
In the current study with S. khuzistanica, the total phenol, total flavonoids and antioxidant activity assayed by DPPH and FRAP were increased in response to B. These results were consistent with the findings reported by Thurzo et al. [60] in Prunus avium L. and Merino-Gergichevich et al. [61] in Vaccinium corymbosum L., who also reported a positive effect of B on the total phenol and flavonoids of studied plants. Given that B has been shown to play an important role in the metabolism of phenolic compounds, the increase in phenol and flavonoids as a result of B application is not surprising [62]. On the other hand, plants with high phenol and flavonoids exhibit high antioxidant activity [27]. In line with these findings, improved antioxidant activities by B application have been reported in other studies by Tavallali et al. [63] and Sarafi et al. [64]. The results of the present study showed the positive impact of B in the content, yield and composition of the S. khuzistanica essential oil, which is in accordance with the results reported by Sugier et al. [65] in Arnica montana L. and Arnica chamissonis. Application of B increased the essential oil content and yield of both A. montana and A.chamissonis. In addition, their essential oil constituents were altered by the B treatment. Choudhary et al. [66] also reported that the foliar application of B increased the essential oil content of Mentha arvensis and Cymbopogon flexuosus. They also mentioned that menthol and methyl acetate percentages in the M. arvensis essential oil and the citral and geraniol percentages in the C. flexuosus essential oil were increased by B application. In general, B fertilization indirectly affects the biosynthesis of primary and secondary metabolites in plants-though different results have been observed in earlier reports. However, in this study, foliar application of B led to an increase in carvacrol, as the predominant component of the S. khuzistanica essential oil only in the first year.

Conclusions
The results of our study suggest that the foliar application of B can be considered a suitable method to increase the quantity and quality of the crop and seed of S. khuzistanica. Application of Zn could be a promising method to improve the growth and yield of S. khuzistanica as it improves yield attributes, drug yield and phytochemical properties. Zn supply also leads to a higher seed set and better seed germination. The results of the current study indicate that the use of several microelements, such as B and Zn, can improve both the quantity and quality of S. khuzistanica. In addition, improving the seed set and seed quality by foliar spraying Zn and B may be useful for growing plants in arid and semi-arid areas. Overall, this study provides some useful information on the effectiveness of the foliar application of Zn and B on S. khuzistanica plant production.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/plants10112469/s1. Table S1: Climatic conditions of S. khuzistanica cultivation area for both cropping years. Table S2: Analysis of variance of the effect of boron and zinc foliar application on biomass and yield attributes of S. khuzistanica in the first year of the experiment. Table S3: Analysis of variance of the effect of boron and zinc foliar application on biomass and yield attributes of S. khuzistanica in the second year of the experiment. Table S4: Analysis of variance of the effect of boron and zinc foliar application on seed traits and B and Zn contents of leaf of S. khuzistanica in the first year of the experiment. Table S5: Analysis of variance of the effect of boron and zinc foliar application on seed traits and B and Zn contents of leaf of S. khuzistanica in the second year of the experiment. Table S6: Analysis of variance of the effect of boron and zinc foliar application on phytochemical traits of S. khuzistanica in the first year of the experiment. Table S7: Analysis of variance of the effect of boron and zinc foliar application on phytochemical traits of S. khuzistanica in the second year of the experiment.  Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.