Effect of rhizosphere bacteria on antioxidant enzymes and some biochemical characteristics of Medicago sativa L. subjected to herbicide stress

Abiotic stresses such as herbicides can affect plant growth and yield. Using herbicide-resistant plant growth-promoting bacteria is a new approach to mitigate these side effects. This study was conducted to evaluate the effect of three native plant growth-promoting bacteria isolated from the Medicago sativa rhizosphere, including Serratia rubidaea (A) and Pseudomonas putida (B), Serratia sp. (C) plus Synorhizobium meliloti (R) and their combinations (AB, AC, BC, ABC, AR, BR, CR, ABR, ACR, BCR, and ABCR) on microbial population, plant biomass, antioxidant enzymes (CAT, APX, and GPX) activities, and hydrogen peroxide and malondialdehyde contents at the presence and absence of imazethapyr herbicide. The results indicated that herbicide application decreased plant biomass but increased microbial population, antioxidant enzymes activities, and the concentrations of hydrogen peroxide and malondialdehyde of all inoculated and non-inoculated plants. Bacterial inoculation in most cases increased microbial population, plant biomass, and antioxidant activities. These increases were more evident under herbicide application. The highest increase in these attributes was achieved by AB, AR, and ABR inoculums in the presence and absence of the herbicide. The microbial population, plant biomass and antioxidant activities were decreased under BC, CR, BCR, and ABCR inoculations. It can be concluded that in addition to growth promotion, these bacteria increase resistance against herbicide stress by controlling free-radical induced oxidative stress and lipid peroxidation through antioxidant enzymes. These findings create new visions in biofertilizer preparation for reducing environmental stresses.


Introduction
Medicago sativa L. (M. sativa) is high quality and valuable cattle forage and a critical element in crop rotation in sustainable agricultural systems (Yanes et al. 2012;Rad et al. 2020;Aghajanlou et al. 2021;Kumar et al. 2020;Moradi et al. 2021). The high quantitative productivity of this forage crop, with low fiber, high protein content, enriched with vitamins and minerals, high tolerance against drought and other environmental stresses, etc., are the most essential factors that grant special status to this forage crop (Putnam and Orloff 2014;Mastinu et al. 2021). The optimal productivity of this crop is very diverse due to M. sativa cultivars and environmental conditions like as low soil fertility, affliction with diseases, presence of weeds, etc., which decrease the yield (Hancock 2015;Ghorbani et al. 2008;Yousefi et al. 2020;Zangani et al. 2021). Pesticide application is the primary method for weed and disease control in plants. The over consumption of these chemical substances and their soil penetration property significantly reduces plant productivity in addition to their harmful environmental impacts, like as underground water pollution and extermination of beneficial soil microorganisms (Volenec et al. 2002;Lamb et al. 2006). Plant damage caused by herbicide application is the result of increased reactive oxygen species. These compounds lead to lipid peroxidation, cell membrane damage, and reactive oxygen production (Silva et al. 2013;Karimmojeni et al. 2021). The cell contents transude due to cell membrane damage and hence specific physiological processes like as photosynthesis will be disrupted (Kruse et al. 2006;Naservafaei et al. 2021). Plants activate their antioxidant defense systems to reduce damages caused by reactive oxygen species. This system involves low-weight molecular compounds like as ascorbic acid, glutathione and antioxidant enzymes, including glutathione peroxidase (GPX), ascorbate peroxidase (APX), catalase (CAT), etc. (Hussain et al. 2016;Gupta et al. 2020Gupta et al. , 2021. The herbicide imazethapyr also known as Pursuit, is a selective systemic herbicide that has been and is being applied to control M. sativa weeds. Studies have revealed that plant-growth-promoting bacteria, colonized in M. sativa rhizosphere, improve productivity in different environmental conditions specially in stressed situations (Carrillo-Castañeda et al. 2002;Ramakrishna et al. 2019). These bacteria promote plant growth by providing nutrients, secretion of plant growth hormones, and many other mechanisms (Chiappero et al. 2019). Another function of plant-growthpromoting bacteria is decomposition of herbicide residues in soil (Akbar and Sultan 2016). Weeds and microbes destroyed by herbicides act as nutritional carbon and nitrogen sources for other soil microorganisms, leading to bacterial recovery and increasing their strength for crop-growth restoration (Lupwayi et al. 2010). Studies in recent years revealed that plant-growthpromoting bacteria are highly contributive in plant resistance against environmental stresses, like as herbicide application, by activating antioxidant mechanism (Khan et al. 2019;Vurukonda et al. 2016). Due to the emphasized role of biological fertilizers in sustainable agriculture in the past few decades, the aim of this study was to investigate the potential of native herbicide-resistant bacteria from M. sativa in prevention of yield reduction in alfalfa through promoting the activity of antioxidant enzymes under herbicide stress. Indeed, the majority of studies on PGPR have been focused on growth promotion potential of them under normal condition while rare studies have been concentrated on their role in mitigation of abiotic stresses in plants.

Materials and methods
The cultivated M. sativa samples in Isfahan University of Technology (located at Lavark, Najaf-Abad, south-west of Isfahan Province, Iran) were harvested in winter. Samples were randomly collected from five points of field at three replicates from each point. The soil surrounding the roots was removed, and the bacterial isolation was done through preparing serial tenfold dilutions from the rhizosphere soil and subsequent culturing on Nutrient agar. The pure isolates were identified based on phenotypic and genotypic (16S rRNA sequencing) identification methods. These isolates were screened based on their plant growth promotion capabilities including nitrogen fixation, auxin production, phosphate, and potassium solubilization, lipase, protease, and cellulose enzymes production, and herbicide resistance. The activity level of these bacteria differed concerning the different characteristics and based on the biological activity of them the best isolates in terms of growth-promoting attributes were A, B, and C, respectively.

Pot experiment
This experiment was conducted in Chah-Anari field in Isfahan University of Technology (40 km south-west of Isfahan, 32º 32´ N 51º 23´ E) as a pot experiment. Hamedani M. sativa cultivar was cultured and subjected to 16 microbial treatments including non-inoculated (control), inoculation with A, B, C, and R (Synorhizobium meliloti) isolates, dual microbial inoculation, triple microbial inoculation, inoculation with four bacterial species, and two levels of herbicide including zero (control) and one liter per hectare as the usual herbicide dose in M. sativa fields. Imazethapyr (Bayer) was applied post-emergence at the 5-6 leaf stage. The herbicide was sprayed in such a way that all parts of the plants were washed and the soil surface was wet. This two-factor factorial experiment was run based on a randomized complete block design with five replications. Twenty seeds were sown in all pots; thinning was performed to have the same bush count in each pot before herbicide treatment and to have uniform germination with optimum density. M. sativa was planted on April 3, 2019, and harvested on June 11, 2019.

Inoculum preparation
To prepare bacterial inoculum, 100 μl of bacterial suspension (10 9 CFU/ml) was added to 20 ml nutrient broth, and the seeds' primary inoculation took place just before cultivation. The secondary inoculation took place one month after cultivation, during the plants initial establishment. For this purpose, approximately 300 ml water was consumed for irrigation of 5-l pot; a proper measure to prevent water discharge. Then, according to the concentration of inoculum in the seed inoculation stage, 1500 μl of bacterial suspension (10 9 CFU/ml) was administered for 300 ml of irrigation water. For rapid bacterial penetration in the roots and preventing drought stress, the pots were irrigated before inoculation.

Biochemical parameters determination
Leaf samples were harvested 48 h after herbicide application and preserved at −80 °C before physiological measurement.

Hydrogen peroxide (H 2 O 2 )
The concentration of H 2 O 2 was measured according to Velikova et al. (2000). Leaf tissue (200) mg was homogenized in 2 ml TCA solution (0.1% w:v) and the homogenate was centrifuged at 4 °C and 10,000 g for 10 min; 0.5 ml of the supernatant was mixed with 0.5 ml potassium phosphate buffer (0.1 M, pH 7.0) and 1 ml KI (1 M) and kept for 1 h at room temperature in a dark place. The adsorption was recorded at 390 nm (Chaichi et al. 2022).

Malondialdehyde (MDA)
Lipid peroxidation was measured according to the method described by Heath and Packer (1968) for malondialdehyde, and Meir et al. (1992) for aldehydes, with minor modifications (Yousefvand et al. 2022). Leaf tissue (200 mg) was homogenized in 2 ml TCA solution (0.1% w:v), and the homogenate was centrifuged at 10,000 rcf and 4 °C for 10 min; 0.5 ml of the supernatant was added to 1 ml of 20% TCA containing 0.5% (w:v) TBA. The reaction mixture was incubated at 95 °C for 30 min. Then the tubes were submerged in an ice bath to stop the reaction. The samples were centrifuged again at 10,000 rcf for 15 min. The lipid peroxidation was measured for MDA and the non-specific absorption samples at 532 and 600 nm, respectively. The absorption of the colored reaction was recorded at 455 nm for aldehydes. Red pigment (The MDA-TBA complex volume) was presented as μmol/g FW and determined by the 155 mM −1 cm 1 extinction coefficient obtained through the following equation:

Determination of Antioxidant enzymes activity
Catalase (CAT), glutathione peroxidase (GPx), and ascorbate peroxidase (APX) activities of M. sativa were measured according to the methods described by (Chance and Maehly 1955;Herzog and Fahimi 1973;Kadkhodaie et al. 2013) with slight modifications. Fresh leaf tissue (0.1 g) was powdered by applying liquid nitrogen and (2) homogenized in 1 ml Na-phosphate buffer (50 mM, pH 7) containing 2 mM a-dithiothreitol, 2 mM EDTA, 0.2% Triton X-100 and 50 mM Tris-HCl in a cold mortar. The sample mixture was centrifuged at 12,000 rpm and 4 °C for 30 min, and the supernatant was used for determining CAT, GPx, and APX enzymes activities.

Catalase enzyme activity (CAT)
This activity was assayed in a total volume of 3 ml of 50 mM Na-phosphate buffer (pH7.0) containing 4.51 μl of H 2 O 2 (30%) and 50 μl enzyme extract (pH 7.8). The absorbance was recorded at 240 nm every 30 s for 2 min to measure H 2 O 2 reduction. Specific catalase activity (enzyme unit) was expressed as the enzyme decomposed per milligram of protein (U mg −1 protein). The volume of the particular catalase activity was obtained through the following equation: U: A unit of catalase activity is equal to the volume of the enzyme that catalysis the H 2 O 2 to O 2 and H 2 O 2 in 1 min; ΔA: Absorbance differences in 240 nm in 1 min; TV: Total bulk of buffer and extract (3 ml); EV: Extract volume (0.05 ml); ε: Extinction coefficient for catalase (39.4 mM −1 cm −1 ); D: Dilution coefficient.

Ascorbate peroxide enzyme activity (APX)
The APX was measured in a total volume of 3 ml of 50 mM Na-phosphate buffer (pH 7.0), containing 4.51 μl of H 2 O 2 (30%), 100 μl of 5 mM ascorbate and 50 μl of enzyme extract (pH 7.8). To evaluate APX, the absorbance , was recorded through a spectrophotometer at 290 nm every 30 s for 2 min. The amount of specific ascorbate peroxide activity was calculated through the same equation applied for specific catalase activity, at ε = 2.8.
The GPX was measured through recording the absorbance at 470 nm every 30 s for 2 min. The volume of specific peroxidase activity was calculated through the same equation applied for specific catalase activity, at ε = 26.6.

Protein content
Protein content was measured by consuming 3 ml of Bradford reagent (consisting of 50 ml ethanol (95%), 100 ml orthophosphoric acid, and 100 mg Coomassie brilliant blue in 1000 ml) and 100 μl of protein extract, mixed and incubated at room temperature for 30 min (Lazzari et al. 2011;Manca et al. 2013). The absorbance was recorded at 595 nm. Bovine serum albumin (BSA) was used as a standard reagent (Bradford 1976) (Lazzari et al. 2012;Sanna et al. 2009).

Microbial population
The M. sativa plants were harvested from the soil 20 days after herbicide application, and the microbial population in the plant rhizosphere was counted. In this process, the soil surrounding the plant root was removed by gentle shaking, and the rhizosphere soil on the root was suspended in saline buffer. The obtained suspension was used to prepare 10 -1 till 10 -5 dilutions, and 0.1 ml of each dilution was streak cultured on nutrient agar. After 24 h incubation at 28 °C, the formed colonies were counted. Shoot dry matter.
To measure the plants' dry weight in each pot, the shoots were dried at 72 °C oven for 48 h and their weights were measured.

Statistical analysis
The obtained data were subjected to statistical analysis (ANOVA) by SAS statistical software (version 9.4), and the means were compared by applying the least significant difference (LSD) at a 0.05 probability level.

Shoot dry matter
The dry shoot matter was affected (at 5% level) by the main and the interaction effects of herbicide and microbial treatments (Table 1). Biomass production decreased by herbicide application in both inoculated and non-inoculated plants (Fig. 1). Microbial inoculation in most cases increased the dry shoot matter under with and without herbicide treatments. The highest increases were found in case of ABR, Ab and AR, respectively (Fig. 2).

Hydrogen peroxide
H 2 O 2 content at 1% level was significantly affected by the main effects and the interaction effects of herbicide and microbial treatments (Table 1). Herbicide application increased the H 2 O 2 content in all inoculated and non-inoculated plants ( Table 2). The H 2 O 2 content in most cases was decreased by microbial treatments. The highest decreases in H 2 O 2 content were observed by AB and ABR, respectively (Table 3).

Malondialdehyde
MDA content at 1% level was significantly affected by the main effects and the interaction effects of herbicide and microbial treatments ( Table 1). The herbicide application increased malondialdehyde content in the sprayed plants, indicating a degree of damage to the plant tissues caused ( Table 2). The MDA content in most cases was decreased by microbial treatments. The highest decreases in MDA content were observed by ABR and AB, respectively (Table 3).

Antioxidant enzyme activity
The results obtained from ANOVA revealed that the main effects of the herbicide and microbial treatments were (at a 1% level) on the activities of CAT, APX, and GPX enzymes. The interaction effects of herbicide and microbial treatment were only significant (at 1% level) in CAT and APX activities (Table1). Herbicide application, increased the activities of CAT, APX, and GPX ( Table 2). The microbial inoculation, in most cases, increased the activities of all three antioxidant enzymes under both the applied and non-applied herbicide conditions (Table 3). The highest increases in CAT were observed by AB, ABR, and AR, respectively. Also, the  highest increases in GPX and APX were observed by AB, ABR, AR, ACR, and AC, respectively (Table 3).

Microbial population
The main and the interaction effects of herbicide and microbial treatments were significant (at 5% level) on microbial population (Table 1). Herbicide application increased the microbial population under all inoculated and non-inoculated treatments except under ABCR inoculation (Fig. 3). Microbial population in most cases increased following inoculation, under both herbicide treatments. The maximum increases in microbial population were observed by ABR and AB, respectively (Fig. 4).

Discussion
Biomass production was decreased as a consequent of herbicide application in both inoculated and non-inoculated plants. According to Bastiaans et al., field application of herbicide parallel to removing the weeds slows the crop growth (Bastiaans et al. 2008). Some herbicides have similar composition to plastoquinone, which can interfere with electron transfer process of photosynthesis and imposing oxidative stress on the plant (Zimdahl 2007  leads to cell membrane lysis through membrane lipids peroxidation Mahdavi et al. 2020); therefore, the plant has to consume part of its energy to produce antioxidant compounds to alleviate the oxidative stress related to herbicide. These consequences reduce plant growth and yield (Zimdahl 2007).
Herbicide application increased the hydrogen peroxide content and consequently led to increase in malondialdehyde in treated plants which indicates a degree of damage to the plant tissues. In the study conducted on the treatment of Salvia officinalis by five herbicides, the concentrations of hydrogen peroxide and malondialdehyde were increased in the herbicide-treated plants (Teimouri Jervekani et al. 2018). More assessments run on the changes of fatty acid profile and antioxidant systems in a Nostoc muscorum exposed to bentazon indicated that bentazon led to increase in reactive oxygen radicals, oxidative stress, and programmed cell death, which ultimately increased the malondialdehyde content of the plant (Galhano et al. 2011) (Figs. 5 and 6).
Application of herbicide was led to increase in microbial population of root rhizosphere in all conditions except under the ABCR treatment. Pozo et al. found that consuming 2-20 kg per hectare of alachlor increased the total population of bacteria and fungi (Pozo et al. 1994). Similarly, study of short-term effects of napropamide on the activity and structure of soil microbial population revealed that the microbial population increased following herbicide application as previously reported (Cycoń et al. 2013). It is assumed that pesticide-resistant microorganisms can metabolize herbicide as carbon and nitrogen sources and form dominant population in their habitat. Moreover, the biomass of sensitive and dead microorganisms act as nutrient resources for these population. Similar to our research, atrazine degrading bacteria were isolated from previously treated field with this herbicide and were able decompose atrazine as carbon and nitrogen sources, leading to an increase in their ability to recover the microbial population in the presence of herbicides (Radosevich et al. 1995).
Bacterial treatments in most cases increased plant biomass and microbial population in root rhizosphere; while, less reductions were found consequent to herbicide application in inoculated plants. The maximum increase in microbial population and plant biomass, in both herbicide and non-herbicide treatments, were obtained following inoculation with A bacterium accompanied with B and R but not C bacterium. Auxin accelerates symbiotic interaction between rhizobium and legume through signaling and leading to improvement in nitrogen fixation, which can eventually stimulate the plant growth (Lin et al. 2020). The plant biomass and microbial population was decreased following BC, CR, BCR, and ABCR bacterial treatments. This may indicate that inoculum C had antagonistic effects, especially when accompanied with inoculums B and R. Bacterium C can produce hydrogen cyanide, a gas that negatively affects root metabolism and is a very strong inhibitor for many metalloenzymes. The microbial population of rhizosphere can be negatively affected due to increase in hydrogen cyanide production.
Some researchers have found that PGPB strains can increase growth and modulate negative effects of herbicides on plant growth. PGPB-related enhanced growth and increased secondary metabolites and antioxidant potential have been reported in Cannabis sativa (Pagnani et al. 2018;Kumar et al. 2019). According to Ahemad and Khan, although consuming Quizalafop-p-Ethyl and Clodinafop herbicides led to decrease in chickpea growth, the plants inoculated with Mesorhizobium sp. MRC4 had more yields (Ahemad and Khan 2009). Ma et al. assessed the improvement in Brassica juncea growth and nickel tolerance by nickel resistant-plant-growth-promoting bacteria (Ma et al. 2009). Their results revealed that Psychrobacter sp. SRA2 significantly increased the fresh and dry biomass of subjected plants. Shahid and Khan found that due to inoculating chickpea with Burkholderia cepacia, a phosphate solubilizing and glyphosate-resistant bacterium, the concentrations of chlorophyll b, a, and carotenoid pigments were increased, as well as the negative effects of glyphosate on chickpea were decreased (Shahid and Khan 2018). In Rego and colleagues study, two rhizobacteria, including Pseudomonas fluorescens BRM-32111 and Burkholderia pyrrocinia BRM-32113, were isolated from rice rhizosphere, and their role on mitigation of toxic allelochemicals and crop yield were assessed in field experiment in the presence of allelopathic rice residues (Rêgo et al. 2018). They found that rice inoculation with these bacteria can significantly reduce the yield drop rate related to toxic allelochemicals. Pandey et al. assessed the potential of metal toxicity amelioration and growthpromoting abilities of Ochrobactrum sp., Bacillus sp. and Bacillus sp.; cadmium, lead and arsenic resistant species, respectively, while applied on rice. Following inoculation of metal-treated rice with these bacterial species, the toxic effect of metal was reduced, and the overall biomass and root/shoot ratio increased (Pandey et al. 2013). Rizvi and Khan reported that the presence of a metal tolerant nitrogen fixing Azotobacter chroococcum in the soil polluted with Cu and Pb, caused growth and yield enhancement of maize (Rizvi and Khan 2018).
Another finding was that in most cases bacterial treatments increased the activity of antioxidant enzymes under both herbicide-treated and -untreated conditions. The maximum increases with regard to plant biomass production and antioxidant enzyme activities were observed in the plants inoculated with AB, AR, and ABR, which highlights the contribution of antioxidant enzyme activation in promoting plant growth by bacterial inoculation. This correlation became evident through the simultaneous reductions in plant biomass and enzyme activities under BC, CR, BCR, and ABCR inoculums. The significant correlations between plant biomass and antioxidant enzyme, lipid peroxidation and microbial population under both herbicide treatments suggest the higher incidence of oxidative stress under herbicide application which greatly stimulated antioxidant enzyme activation to reduce herbicide stress in plants. The higher increase in CAT, APX, and GPX activity induced by the most beneficial bacterial inoculations (AB, AR, and ABR) may suggest high contribution of these enzymes in plant biomass improvement through inoculation under herbicide application. This is confirmed with high negative correlations between antioxidant enzymes activity and free radical production and lipid peroxidation (Tables 4 and 5).
In the study of Shaid and Khanon, the productivity of glyphosate exposed chickpea, the chickpea was inoculated with glyphosate-tolerant phosphate solubilizing Burkholderia cepacia and the antioxidant enzymes of the herbicidetreated plant was measured (Shahid and Khan 2018). In the study of Pramanik et al the role of cadmium-resistant PGPR strain named as Enterobacter aerogenes MCC 3092 in reduction of cadmium toxicity in rice was evaluated and it was found that superoxide dismutase and catalase activity in inoculated plants with Enterobacter aerogenes was higher than the non-inoculated seedlings (Pramanik et al. 2018). Islam et al. assessed the effectiveness of resistant Pseudomonas aeruginosa on wheat productivity in soils with high concentrations of zinc. Their results revealed that total biomass, antioxidant enzyme activity (SOD, POD, and CAT), and antioxidant compounds, including ascorbic acid and phenolic compounds, increased in inoculated plants (Islam et al. 2014). Under herbicide stress, the plant with the more potent antioxidant system will be subjected to less lipid peroxidation with free radicals which leading to low MDA production. The high antioxidant activity in the pesticideresistant plants can indicate that these plants can maintain their growth and productivity at a normal level by neutralizing the harmful oxidants and preventing severe damage to membrane. The results of this study are consistent with several similar studies. Assessing the effect of herbicidetreated soil on the antioxidant system of sweet corn vegetative organs, revealed that the herbicide application increases the activities of catalase, superoxide dismutase, and peroxidase in the leaves and roots (Grigoryuk et al. 2016). Mohsin et al (2020) pretreated wheat seedlings with 2,4-D to assess the effect of this toxin on the improvement of the plant tolerance against oxidative stress caused by salinity and tolerance to methylglyoxal toxicity by activating the antioxidant defense system, ion homeostasis, and the glyoxalase system. They found that the activity level of antioxidant enzymes, ascorbic acid, and glutathione content, were increased in the 2,4-D-treated seedlings (Mohsin et al. 2020).
This fact that microbial population has a positive correlation with plant biomass and antioxidant enzyme activities and has negative correlation with free radical production and lipid peroxidation (Tables 4 and 5), indicates the contribution of rhizosphere microbial population in alleviating oxidative stress through potentiating antioxidant enzyme activities in plants. The highest decrease in hydrogen peroxide and MDA contents and the highest increase in antioxidant enzymes' activities were observed in the plants inoculated with AB, AR, and ABR, especially under herbicide application. The hydrogen peroxide and MDA contents were high under BC, CR, BCR, and ABCR inoculations, while antioxidant enzymes activities were low under the same bacterial treatments. The results revealed enzymatic antioxidant defense system's activation in inoculated plants to control lipid peroxidation and oxidative stress induced by herbicide application. In general, increases in microbial population, plant biomass, and antioxidant enzymes activities, and consequently decreases in hydrogen peroxide and MDA concentrations due to bacterial inoculation, were higher under herbicide conditions. Pandey et al. used three different bacterial isolates (a cadmium-resistant Ochrobactrum sp., a lead resistant Bacillus sp., and an arsenic resistant Bacillus sp.) as inoculum for metal-treated rice (Pandey et al. 2013). They found that the increased of MDA level in rice roots, due to metal stress, was decreased following bacterial inoculation. Inoculation of metal tolerant Azotobacter chroococcum in Cu and Pb contaminated corn field, caused decreasing in MDA levels (Rizvi and Khan 2018).
Excessive consumption of herbicides is one of the most important problems in conventional agriculture, affecting soil fertility negatively and cause considerable reduction in the crop yield. Most of the microbial inoculations obtained from the M. sativa rhizosphere can increase the crop dry weight, enhance the antioxidant system against the herbicide stresses and decrease the damage caused by herbicide. It can be concluded that there exists a direct relation between the population of the bacterial treatment in the soil and the antioxidant activity and the defensive power of the plant against the herbicide. This phenomenon can prevent the drastic drop in the crop yield by decreasing the hydrogen peroxide contents; thus, an inhibitory function in the peroxidation of membrane lipids. As these plant growth-promoting bacteria introduced and used in this study, their ability to enhance plant stability against abiotic stresses may make them as appropriate candidates for biofertilizers in farmlands.
Author contribution statement MM, MZ, and HK conceived the study; MZ, HK, and HM performed biochemical, anatomical and morphological analysis, and helped interpret results. MM, MZ, and HK provided plant care, performed histological analysis and helped interpret results. MZ, HK, and AM interpreted results and wrote original draft of the manuscript; HK and AM performed statistical analysis, review and editing of the manuscript. included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.