Co-application of ACC-deaminase producing PGPR and timber-waste biochar improves pigments formation, growth and yield of wheat under drought stress

Besides other deleterious effects, drought elevates ethylene level too in plants. Increased ethylene concentration reduces root elongation and development that consequently retard plant growth and yield. There are certain PGPR which produce ACC-deaminase. The ACC-deaminase converts ACC (an immediate precursor of ethylene biosynthesis in methionine pathway in higher plants) into ammonia and α-ketobutyrate instead of ethylene. Regularization of ethylene level in plants mitigate the effects of drought. On the other hand, biochar has been reported to be rich in nutrients and exhibiting higher water holding capacity. So, a pot study was conducted with the hypothesis that the combined application of ACC-deaminase producing PGPR and biochar would minimize the drought effects on wheat growth. The ACC-deaminase producing PGPR were applied on wheat seeds in combination with two biochar doses. Three moisture levels were maintained throughout the trial. The data obtained revealed that B. amyloliquefaciens + 2BC improved the chlorophyll a, chlorophyll b, photosynthetic rate, transpiration rate, 100-grain weight, and grain N, P and K up to 114%, 123%, 118%, 73%, 59%, 58%, 18% and 23%, respectively, under drought conditions. It is concluded that co-application of PGPR and biochar is an effective technique to mitigate the drought effects.

Various biotic (pests, pathogens) and abiotic (soil compaction, drought, salinity, waterlogging, heavy metals, poor nutrition etc.) stresses are a big cause of low crops productivity around the globe 1 . Drought stress is very common in worldwide arid and semi-arid areas. Moreover, climate change is going to create the worst situation in this regard [2][3][4][5][6] . The demand for irrigation water is expected to increase by 10% up to 2050 7 . Under drought stress, growth and yield of crops are usually decreased due to less intake of nutrients, poor photosynthesis 8 and limited supply of water 9 . In addition, drought accelerates the biosynthesis of ethylene 10,11 which retards the roots elongation and development [12][13][14] .
Although, traditional breeding, water management and genetic engineering are thought to be useful tools to alleviate drought stress but high technicalities are involved to adopt and implement these approaches 15 . However, the use of plant growth promoting rhizobacteria (PGPR) is an alternative technique for mitigation of drought effects 15 . A large number of rhizospheric bacteria are well documented that show growth promotion in plants under stressful conditions 14 . As far as regularization of ethylene biosynthesis under drought stress is concerned, using 1-aminocyclopropane-1-carboxylate deaminase (ACC-deaminase) producing PGPR is found to be quite effective [16][17][18][19][20] . The ACC-deaminase cleaves the ACC (1-aminocyclopropane-1-carboxylic acid, an immediate precursor of ethylene biosynthesis through methionine pathway in higher plants) into ammonia and α-ketobutyrate instead of ethylene 11,21,22 . Besides ethylene regularization, the PGPR also help in better root development 23 , secretion of growth hormones (auxins or cytokinins) 22 and solubilization of immobile nutrients (phosphorus, potassium etc.) 24,25 .
On the other hand, the imperative role of organic amendments in mitigation of drought stress by improving soil water holding capacity and availability of nutrients cannot be denied 26,27 . Biochar (BC), is a black carbon compound which is a good source of nutrients. It is produced through pyrolysis at high temperature under low or no supply of oxygen [27][28][29] . The physio-chemical properties of BC depend on the nature of waste material used and temperature of the pyrolysis 30,31 . High surface area and pore spaces of BC-structure improve soil water and nutrients holding capacity [32][33][34][35] .
Wheat (Triticum aestivum L.) is an important cereal crop and staple food in most parts of the world. It contains 55% carbohydrates and 8-12% proteins 36 . It is an important crop due to its worldwide trade too 37 . Cultivation of wheat under a limited supply of water significantly decreases the yield 38 while its demand is increasing at the rate of 1.6%/annum 39 . The need of time is to enhance wheat yield even in the areas under drought stress.
In recent past, the researchers focused on the application of either ACC-deaminase containing PGPR or BC in separate to mitigate the drought stress. The novelty and aim of the present study are to examine the combined effect of ACC-deaminase producing PGPR and timber-waste BC for the alleviation of drought effects. Keeping in mind the importance of wheat, the current study was conducted with the hypothesis that the co-application of drought tolerant ACC deaminase producing PGPR and timber waste BC could be very effective to alleviate drought effects.

Results
Gas exchange attributes. Main effects of treatments (T) and various levels of drought (D) were significantly different, while their interaction (T × D) remained similar for the rate of photosynthesis and transpiration. For stomatal conductance, both main and interactive effects of T and D were significantly different. The photosynthetic rate was significantly improved as compared to the control where P. aeruginosa and B. amyloliquefaciens were applied ( Table 1). The BC application without PGPR significantly enhanced photosynthetic rate as compared to the control but the application rate 2BC was more effective than 1BC for the improvement in photosynthetic rate. However, among all the treatments P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC remained the best to significantly increase the photosynthetic rate. For transpiration rate, a significant improvement was noted as compared to the control in all the treatments. Application of 1BC and 2BC without PGPR gave statistically similar results regarding transpiration rate. However, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC remained significantly better for transpiration rate as compared to L. adecarboxylata, A. fabrum, P. aeruginosa, B. amyloliquefaciens, 1BC and the control. Maximum increase in photosynthetic rate (118%) and transpiration (73%) rate was noted as compared to control where B. amyloliquefaciens + 2BC was applied. Application of 1BC and 2BC remained significantly better at severe drought (SD) level as compared to the control for stomatal conductance. However, at mild drought (MD) level, application of 2BC remained significantly better than 1BC and control regarding stomatal conductance. The B. amyloliquefaciens + 2BC at SD level while 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC at MD level performed significantly better among other treatments for stomatal enoids and proline contents. In case of carotenoids, the treatments, L. adecarboxylata + 2BC, A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC were statistically alike but differed significantly as compared to the control for carotenoids at SD (Fig. 2). The 2BC proved significantly better as compared to 1BC and the control for carotenoids at SD. The 1BC, L. adecarboxylata + 1BC, A. fabrum + 1BC, P. aeruginosa + 1BC, B. amyloliquefaciens + 1BC, 2BC, L. adecarboxylata + 2BC, A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC were statistically alike but significantly different as compared to the control for carotenoids at MD (Fig. 2a). Maximum increase i.e., 42%, 48% and 220% in carotenoids was noted at NM, MD and SD respectively as compared to the control where B. amyloliquefaciens + 2BC was applied. For proline reduction, A. fabrum, P. aeruginosa and B. amyloliquefaciens differed significantly as compared to the control at SD. The 1BC significantly decreased proline too as compared to control, but the 2BC was statistically far better for proline reduction as compared to 1BC at SD. The PGPR strains remained significantly better with 2BC for proline reduction as compared www.nature.com/scientificreports www.nature.com/scientificreports/ to with 1BC at SD (Fig. 2b). However, the B. amyloliquefaciens + 2BC showed maximum decrease in proline i.e. 34% at SD as compared to the control.
Photosynthetic pigments. Main effects of T and D were significantly different but their interactions remained similar for chlorophyll a, chlorophyll b and total chlorophyll contents in wheat leaves. The treatments L. adecarboxylata, A. fabrum, P. aeruginosa and B. amyloliquefaciens were statistically alike but differed significantly as compared to the control for chlorophyll a content. With and without PGPR, the application of 1BC and 2BC gave statistically similar results but remained significantly different as compared to the control for chlorophyll a ( Table 2). For chlorophyll b contents, with and without 1BC, all the PGPR treatments (except B. amyloliquefaciens + 1BC) were statistically alike but differed significantly as compared to the control. Without PGPR, the 2BC was significantly better than 1BC as compared to the control for chlorophyll b contents. However, A. fabrum + 2BC and B. amyloliquefaciens + 2BC proved significantly better as compared to other treatments regarding chlorophyll b contents. In case of total chlorophyll contents, the treatments L. adecarboxylata, A. fabrum, P. aeruginosa and B. amyloliquefaciens differed significantly as compared to the control. Statistically, 2BC was significantly better than 1BC for total chlorophyll. However, the 1BC and 2BC significantly increased total chlorophyll as compared to the control. The 2BC, L. adecarboxylata + 2BC, A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC were statistically alike and proved better among the other treatments for total chlorophyll contents. Maximum increase in chlorophyll a (114%), chlorophyll b (123%) and total chlorophyll (115%) was noted through B. amyloliquefaciens + 2BC as compared to the control. N, P and K in shoot. Main effects of T and D were significantly different but their interaction remained similar for nitrogen (N) and phosphorus (P) concentration in shoot. For potassium (K) concentration in shoot, both main and interactive effects of T and D were significantly different. The treatments, L. adecarboxylata, A. fabrum, P. aeruginosa and B. amyloliquefaciens differed significantly as compared to the control for shoot N concentration. The 2BC remained significantly better as compared to 1BC for N concentration in shoot ( Table 3). The PGPR (L. adecarboxylata, A. fabrum, P. aeruginosa and B. amyloliquefaciens) with 2BC performed significantly better as compared to 1BC (except B. amyloliquefaciens + 1BC) for N concentration in shoot. For phosphorus concentration in shoot, 1BC and 2BC significantly differed as compared to the control. In case of P concentration in shoot, L. adecarboxylata, A. fabrum and B. amyloliquefaciens were statistically alike but significantly different as compared to the control. The 2BC significantly increased (i.e., by 28%) P concentration in shoot as compared to 1BC. The treatments A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC proved significantly better for P concentration in shoot. Maximum increase in N (73%) and P (150%) concentration in shoot was noted as www.nature.com/scientificreports www.nature.com/scientificreports/ compared to the control where B. amyloliquefaciens + 2BC was applied. In case of K concentration in shoot, all the treatments remained significantly better as compared to control at SD. The treatments 1BC and 2BC remained statistically alike without PGPR for K concentration in shoot at MD and SD. The B. amyloliquefaciens + 1BC, A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC were significantly better for K concentration Treatments Chlorophyll a (mg g −1 ) Chlorophyll b (mg g −1 ) Total Chlorophyll (mg g −1 ) Various levels of drought (D)    www.nature.com/scientificreports www.nature.com/scientificreports/ in shoot at MD. Maximum increase of 38%, 85% and 138% in K concentration in shoot was noted where B. amyloliquefaciens + 2BC was applied as compared to the control at NM, MD and SD, respectively. Yield attributes. Both main and interactive effects of T and D were significantly different for grain yield pot −1 and straw yield pot −1 . In case of 100-grain weight, main effects of T and D were significantly different but their interaction remained statistically similar. The PGPR with and without 1BC, as well as, 1BC and 2BC were statistically at par with the control at SD for grain yield pot −1 . However, the L. adecarboxylata + 2BC, A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC differed significantly as compared to the control for grain yield pot −1 at MD and SD ( Table 4). The treatments B. amyloliquefaciens + 1BC also remained significantly better as compared to the control and other PGPR with 1BC at MD for grain yield pot −1 . At NM, the B. amyloliquefaciens + 1BC, 2BC, L. adecarboxylata + 2BC, remained significantly different as compared to the control for grain yield pot −1 . Maximum increase, 40%, 155% and 215% in grain yield pot −1 was noted at NM, MD and SD as compared to the control where L. adecarboxylata + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC were applied, respectively. For 100-grain weight, all the treatments (except A. fabrum) remained significantly better as compared to the control. The 1BC and 2BC remained statistically alike for 100-grain weight. The 1BC, A. fabrum + 1BC, B. amyloliquefaciens + 1BC, 2BC, L. adecarboxylata + 2BC, A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC remained significantly better for 100-grain weight. Maximum increase (59%) in 100-grain weight was noted as compared to the control where B. amyloliquefaciens + 2BC was applied . In case of straw yield pot −1 , the PGPR with and without BC remained significantly better as compared to the control at NM, MD and SD. However, 1BC and 2BC remained statistically alike for straw yield pot −1 . The L. adecarboxylata + 1BC, A. fabrum + 1BC, P. aeruginosa + 1BC, B. amyloliquefaciens + 1BC, 2BC, L. adecarboxylata + 2BC, A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC remained significantly better for straw yield pot −1 . Maximum increase i.e., 181% and 178% in straw yield pot −1 was noted as compared to the control where B. amyloliquefaciens + 2BC and A. fabrum + 2BC were applied at SD and MD respectively. N, P and K in grains. Both main and interactive effects of T and D were significantly different for N, P and K concentration in grain. At SD, the L. adecarboxylata, A. fabrum, P. aeruginosa and B. amyloliquefaciens with and without 1BC and 2BC differed significantly as compared to the control for grain nitrogen ( Table 5). The B. amyloliquefaciens, 1BC, L. adecarboxylata + 1BC, A. fabrum + 1BC, P. aeruginosa + 1BC and B. amyloliquefaciens + 1BC, 2BC, L. adecarboxylata + 2BC, A. fabrum + 2BC, P. aeruginosa + 2BC and B. amyloliquefaciens + 2BC remained significantly better as compared to the control for grain nitrogen at MD. Application of A. fabrum + 2BC proved significantly better as compared to the control for grain nitrogen at NM. For grain phosphorus, all the treatments (except P. aeruginosa) were significantly different as compared to the control at SD. However, at MD and NM all the treatments were statistically similar to the control for grain P. In case of grain K, all the treatments were significantly different as compared to the control at SD. The A. fabrum and P. aeruginosa at 2BC remained significantly better than at 1BC for grain K at SD. At MD, inoculation of PGPR with 1BC and 2BC  www.nature.com/scientificreports www.nature.com/scientificreports/ proved significantly better than the control for grain K. The performance of A. fabrum was significantly better at 2BC than at1BC for grain K at MD. Among all the treatments A. fabrum + 1BC, B. amyloliquefaciens + 1BC, 2BC, L. adecarboxylata + 2BC, A. fabrum + 2BC, P. aeruginosa + 2BC, B. amyloliquefaciens + 2BC remained significantly better at NM. The 2BC proved significantly better than 1BC for grain K at NM, MD and SD. Maximum increase in grain N (58%), P (18%) and K (23%) was noted through B. amyloliquefaciens + 2BC as compared to the control at SD.

IE (T × D) (Means of 3 replicates) ME (T) IE (T × D) (Means of 3 replicates) ME (T) IE (T × D) (Means of 3 replicates) ME (T) NM
Discussion. In the current study, reduction in shoot length of wheat at MD and SD in the control might be due to competition for water and nutrients between root and shoot. Gargallo-Garriga et al. 40 stated that drought stress deactivated shoot metabolic activity to conserve water and food which would have facilitated roots elongation. The reduction in water and nutrients movement in shoot, enhanced the up-regulation of ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) 41 while significant amount of ACC generated ethylene restricted the root elongation 12 . Glick et al. 11 suggested the mechanism for reduction in stress ethylene through the activity of ACC-deaminase 11 . According to Glick et al. 11 , the synthesis of indole acetic acid (IAA) by PGPR stimulates the elongation of plants cells and activates ACC synthase which converts S-adenosyl methionine to ACC. A significant amount of ACC is exuded by plants roots and seeds in rhizosphere which is hydrolyzed by ACC deaminase into NH 3 and α-ketobutyrate, resulting in better roots elongation. Secretion of roots exudates (phytosiderophores, sugars, organic acids, amino acids, vitamins, nucleosides and mucilage) also attracts PGPR that colonize roots and facilitates better uptake of water and nutrients 42,43 . Besides imperative role of PGPR, high surface area, ion exchange capacity, nutrients and water holding capacity of BC make it an effective amendment for better intake of nutrients and water in plants [44][45][46] . A significant improvement in photosynthetic rate, transpiration rate and stomatal conductance under MD and SD, signified the effectiveness of co-application of L. adecarboxylata, A. fabrum, P. aeruginosa, B. amyloliquefaciens with 2BC, comparative to their sole application and the control. The increase in gas exchange attributes might be due to the better uptake of water and nutrients, improvement in soil water holding capacity (WHC), PGPR colonization and reduction in ethylene through the co-application of 2BC 47,48 and ACC deaminase containing PGPR 49,50 . According to Zheng et al. 51 and Borch et al. 52 the deficiency of nitrogen and phosphorus significantly decrease the growth of crops. Siddique et al. 53 suggested the less stomatal conductance as main cause of reduction in rate of transpiration under drought. Tholen et al. 54 argued that drought stress decreased the intake of CO 2 due to less stomatal conductance which restricted the carboxylation. Hence reduced the rate of photosynthesis 55 . However, A significant improvement in chlorophyll a, chlorophyll b, total chlorophyll, shoot and grain nutrients concentration validated the efficacious functioning of co-application of ACC deaminase producing PGPR and 2BC that also significantly increased the yield attributes (100-grains weight, straw yield and grains yield pot −1 ) of wheat plants at MD and SD. Richardson et al. 56 and Zahir et al. 57 , suggested better roots elongation and secretion of organic acids by PGPR for P and K solubilization as key factors, responsible for better nutrients uptake, improvement in dry weight and yield of crops. Moreover, the biochar ability to sorb nutrients also reduced the losses of N and improved its uptake in plants 58,59 . Chan et al. 60 suggested   www.nature.com/scientificreports www.nature.com/scientificreports/ that the high surface area of biochar is a basic reason for improved cation exchange sites in soil which resulted in better bioavailability of nutrients. Results of the current study also showed that electrolyte leakage in the wheat plants leaves was decreased where L. adecarboxylata, A. fabrum, P. aeruginosa, B. amyloliquefaciens with 2BC were applied comparative to the control under MD and SD. The reduction in electrolyte leakage might be due to the activity of ACC deaminase, better availability of water and nutrients by co-application of PGPR and 2BC. Senaratna and McKersie 61 observed a significant increase in electrolyte leakage because of cell membrane damage by drought stress which made it more permeable. According to Matile et al. 62 , cell usually lost its membrane integrity as a result of lipid degradation by ethylene. When the lipids in cell membrane become degraded, the ethylene comes in contact with the chloroplast and activates the chlorophyllase (chlase) gene which causes severe damage to the chlorophyll 62 . However, the addition of 2BC and PGPR in combination significantly decreased the electrolyte leakage which might be due to activity of ACC deaminase, better availability of water and nutrients uptake.

Conclusion.
It is concluded that the co-application of drought tolerant ACC deaminase producing PGPR and 2BC is comparatively better approach than their sole application to mitigate drought stress in wheat. Though Leclercia adecarboxylata and Pseudomonas aeruginosa were also effective enough but Agrobacterium fabrum and Bacillus amyloliquefaciens with 2BC gave maximum increase in the gas exchange attributes, nutrients concentration in shoot and grain, photosynthetic pigments and yield of wheat. More investigations are needed at field level to introduce A. fabrum and B. amyloliquefaciens with 2BC to improve growth and yield of wheat under drought stress.
www.nature.com/scientificreports www.nature.com/scientificreports/ For the production of biochar, timber waste was collected from local timber market. The timber waste was initially sun-dried and then pyrolyzed at 389 °C for 80 min in an especially designed pyrolyzer as described by Qayyum et al. 27 . All the pyrolyzed material (biochar) was then crushed in a grinder and passed through 2 mm sieve. Finally, the fine powder of timber waste biochar (BC) was stored in air tight plastic jars 27 .
The pH and ECe of BC were determined by mixing the BC and water with the ration, 1:20 (w/v) as described by Qayyum et al. 27 . Di-acid (HNO 3 : HClO 4 ) digestion 71 of biochar was done for the analysis of total phosphorus following yellow color method on spectrophotometer 72 , and those of potassium and sodium on flame photometer 73 . For the determination of nitrogen, H 2 SO 4 digestion 72 was done followed by distillation on Kjeldahl's distillation apparatus 74 . The volatile matter and ash content of biochar were analyzed according to Qayyum et al. 75 by heating the biochar in muffle furnace at 450 °C and 550 °C respectively. The fixed carbon in biochar was assessed (Table 6) The plastic bag (30 cm deep × 20 cm in diameter) was used as a pot, having capacity to carry 8 kg soil. The soil was collected from the plough layer of bank of the Chenab River, Multan, Punjab, Pakistan. The soil of selected area was previously characterized as dark yellowish brown, moderately calcareous, weakly structured and well drained with Cambic subsurface horizon and an Ochric epipedon 77 . The soil texture was determined by hydrometer method 78 which was sandy loam (USDA triangle) with mixed hyperthermic Haplocambids. The organic matter in soil was determined by Walkley 79 . The organic nitrogen in soil was determined using the equation: = Organic N (%) Soil Organic Matter/20 (2) For extractable soil P determination, Olsen and Sommers 80 method was used. Similarly, the extractable K in soil was determined according to the method described by Nadeem et al. 73 (Table 6). In each plastic pot, 8 kg soil was filled. To fulfil macro nutrients requirement nitrogen (N), phosphorus (P) and potassium (K) fertilizers were added at the rate of 120: 90 and 60 kg ha −1 respectively, as recommended dose keeping in mind the nutrients concentration of biochar where it was applied 81 . The urea was added in three split doses. As far as diammonium phosphate (DAP) and muriate of potash (MOP) fertilizers are concerned, the recommended rates of fertilizers were applied in a single dose at the time of sowing. Timber waster biochar was added at three different rates including: control i.e., no biochar (0BC), 0.75% of soil (60 g biochar per 8 Kg soil) biochar (1BC) and 1.50% of soil (120 g biochar per 8 Kg soil) biochar (2BC).
The seeds of wheat (Galaxy-2013) were obtained from the certified seed dealer of the Government of Punjab, Pakistan. Healthy seeds were separated from broken and weak seeds. The seeds were surface-sterilized with  www.nature.com/scientificreports www.nature.com/scientificreports/ sodium hypochlorite (5%) followed by 3 washes with ethanol (95%). Finally, all the seeds were washed three times with sterilized deionized water 82 . For PGPR inoculation, 10 ml of inoculum (0.5 optical density at 535 nm wavelength) 83 was added along 10% sugar (glucose) in 100 g sterilized seeds. After proper mixing of seeds, inoculum and sugar solution, top dressing of seeds was done with a mixture of peat and clay (3:1 ratio) as described by Ahmad et al. 84 . Before inoculation of seeds, the peat and clay mixture was sterilized at 121 °C for 20 min in an autoclave 83 . All the control treatment seeds were also top dressed with peat and clay mixture along with 10% sugar solution without inoculum 85 .
In each pot, 10 seeds of wheat were initially sown. In control, the soil normal moisture (NM) was maintained at the level of 70% of field capacity (FC 70 ) throughout the experiment on weight basis. However, to introduce mild drought (MD) and severe drought (SD) stress as per treatment plan, the soil moisture was maintained at the level of 50% and 30% of field capacity (FC 50 and FC 30 ), respectively, throughout the trial as suggested by Boutraa et al. 86 . After germination of seeds, five healthy seedlings were kept in each pot by thinning.
Leaf gas exchange parameters (net photosynthetic rate, net transpiration rate and stomatal conductance) were determined with the help of Infra-Red Gas Analyzer (CI-340 Photosynthesis system, CID, Inc. USA) by joining 4 leaves of wheat. On a sunny day, the readings were taken between 10:30 and 11:30 AM at saturating intensity of light 87 .
After 50 days of sowing, the seedlings were harvested from each pot for the measurement of shoot length and determination of electrolyte leakage, proline contents, photosynthetic pigments level and nutrients concentration in the shoot.
The electrolyte leakage (EL) was determined following the procedure given by Lutts et al. 88 . The leaves were washed with deionized water and then cut using a steel cylinder having diameter 1 cm. One gram of uniform sized leaf pieces were immersed in a test tube containing deionized water (20 ml) and incubated at 25 °C for 24 h. The electrical conductivity (EC1) was determined using pre-calibrated EC meter. The second EC (EC2) was noted heating the test tubes in a water bath at 120 °C for 20 min. The final value of EL was calculated using the equation as follows; = × Electrolyte Leakage (EL) EC1/EC2 100 (3) For proline assessment in wheat leaves, methodology stated by Bates et al. 89 was followed. The proline was extracted from fresh (0.1 g) leaves in 2 ml of 40% methanol. After extraction, the 1 ml mixture of glacial acetic acid and 6 M orthophosphoric acid (3:2 v/v) was mixed in 1 ml extract along with 25 mg ninhydrin. Then the solution was incubated at 100 °C for 60 min. After cooling down, 5 ml Toluene was added. For the estimation of proline contents, absorbance was noted on spectrophotometer at 520 nm wavelength.
The chlorophyll a, chlorophyll b and total chlorophyll contents were determined in the fresh leaves of wheat according to the protocol given by Arnon 90 . The extract was taken from the leaves using acetone (80%) solution. For the estimation of chlorophyll a and chlorophyll b, the absorbance was taken at 663 and 645 nm wavelength, respectively on spectrophotometer. Final calculations were made using the following relations; = . .
= + Total Chlorophyll (mg/g) Chlorophyll a Chlorophyll b (6) where, OD = Optical density (wavelength). V = Final volume made. W = Fresh leaf made (g). The samples were digested with sulfuric acid 72 followed by distillation on Kjeldahl's distillation apparatus 74 . The yellow colour method was used for the determination of phosphorus concentration noting absorbance at 420 nm on spectrophotometer 72 . As far as the K concentration in wheat shoot and grain is concerned, the samples were digested and then run on flame photometer as described by Nadeem et al. 73 .
The wheat plants were harvested after 125 days of sowing for the determination of grains yield pot −1 , straw yield pot −1 and 100-grain weight. Weight of 100 grains, straw and grains yield pot −1 were assessed on top weight balance. For straw yield, plants were harvested at 4 inches above the ground surface. Sun dried 100 grains of wheat were counted randomly and manually and then weighed on top weight balance. Total wheat grains collected from a single pot were weighed and considered as grain yield pot −1 .
Statistical analyses of the data were carried out using standard statistical procedures 91 . All the treatments were compared using Tukey's test at p ≤ 0.05.

Data Availability
No datasets were generated or analyzed during the current study. All the analyzed data can be accessed after publication by requesting the corresponding author.