Differential Root Exudation and Architecture for Improved Growth of Wheat Mediated by Phosphate Solubilizing Bacteria

Phosphorous (P) deficiency is a major challenge faced by global agriculture. Phosphate-solubilizing bacteria (PSB) provide a sustainable approach to supply available phosphates to plants with improved crop productivity through synergistic interaction with plant roots. The present study demonstrates an insight into this synergistic P-solubilizing mechanism of PSB isolated from rhizosphere soils of major wheat-growing agro-ecological zones of Pakistan. Seven isolates were the efficient P solubilizers based on in vitro P-solubilizing activity (233-365 μg ml–1) with a concomitant decrease in pH (up to 3.5) by the production of organic acids, predominantly acetic acid (∼182 μg ml–1) and gluconic acid (∼117 μg ml–1). Amplification and phylogenetic analysis of gcd, pqqE, and phy genes of Enterobacter sp. ZW32, Ochrobactrum sp. SSR, and Pantoea sp. S1 showed the potential of these PSB to release orthophosphate from recalcitrant forms of phosphorus. Principal component analysis indicates the inoculation response of PSB consortia on the differential composition of root exudation (amino acids, sugars, and organic acids) with subsequently modified root architecture of three wheat varieties grown hydroponically. Rhizoscanning showed a significant increase in root parameters, i.e., root tips, diameter, and surface area of PSB-inoculated plants as compared to uninoculated controls. Efficiency of PSB consortia was validated by significant increase in plant P and oxidative stress management under P-deficient conditions. Reactive oxygen species (ROS)-induced oxidative damages mainly indicated by elevated levels of malondialdehyde (MDA) and H2O2 contents were significantly reduced in inoculated plants by the production of antioxidant enzymes, i.e., superoxide dismutase, catalase, and peroxidase. Furthermore, the inoculation response of these PSB on respective wheat varieties grown in native soils under greenhouse conditions was positively correlated with improved plant growth and soil P contents. Additionally, grain yield (8%) and seed P (14%) were significantly increased in inoculated wheat plants with 20% reduced application of diammonium phosphate (DAP) fertilizer under net house conditions. Thus, PSB capable of such synergistic strategies can confer P biofortification in wheat by modulating root morphophysiology and root exudation and can alleviate oxidative stress under P deficit conditions.

Phosphorous (P) deficiency is a major challenge faced by global agriculture. Phosphatesolubilizing bacteria (PSB) provide a sustainable approach to supply available phosphates to plants with improved crop productivity through synergistic interaction with plant roots. The present study demonstrates an insight into this synergistic P-solubilizing mechanism of PSB isolated from rhizosphere soils of major wheat-growing agro-ecological zones of Pakistan. Seven isolates were the efficient P solubilizers based on in vitro P-solubilizing activity (233-365 µg ml −1 ) with a concomitant decrease in pH (up to 3.5) by the production of organic acids, predominantly acetic acid (∼182 µg ml −1 ) and gluconic acid (∼117 µg ml −1 ). Amplification and phylogenetic analysis of gcd, pqqE, and phy genes of Enterobacter sp. ZW32, Ochrobactrum sp. SSR, and Pantoea sp. S1 showed the potential of these PSB to release orthophosphate from recalcitrant forms of phosphorus. Principal component analysis indicates the inoculation response of PSB consortia on the differential composition of root exudation (amino acids, sugars, and organic acids) with subsequently modified root architecture of three wheat varieties grown hydroponically. Rhizoscanning showed a significant increase in root parameters, i.e., root tips, diameter, and surface area of PSB-inoculated plants as compared to uninoculated controls. Efficiency of PSB consortia was validated by significant increase in plant P and oxidative stress management under P-deficient conditions. Reactive oxygen species (ROS)-induced oxidative damages mainly indicated by elevated levels of malondialdehyde (MDA) and H 2 O 2 contents were significantly reduced in inoculated plants by the production of antioxidant enzymes, i.e., superoxide dismutase, catalase, and peroxidase. Furthermore, the inoculation response of these PSB on respective wheat varieties grown in native soils under greenhouse conditions was positively correlated with improved plant growth and soil P contents. Additionally,

INTRODUCTION
Phosphorus (P) is the second most essential macronutrient after nitrogen (N) for plant growth and development (Shafi et al., 2020). P-deficiency is one of the yield-limiting factors faced by over 40% of agricultural soils worldwide (Elhaissoufi et al., 2020;Brownlie et al., 2021), although P is quite abundant in the soil but mostly present in insoluble forms, rendering a very low concentration (0.1%) of P for plant uptake (Bindraban et al., 2020).
Chemical P fertilizers are produced at a cost of 4 billion USD per year to fulfill the global demand Tewari et al., 2020). However, P use efficacy of applied chemical P fertilizers is as low as 10% throughout the world due to rapid P immobilization by fixation, adsorption, or precipitation reactions with Al 3+ and Fe 3+ in acidic soils and with Ca 2+ in alkaline to normal soils (Penn and Camberato, 2019;Mitter et al., 2021). The excessive and irrational application of chemical fertilizers is not only expensive but also leads to environmental pollution and adverse effects on human health (Manisalidis et al., 2020). Therefore, there is an urgent need to exploit such strategies that can increase P availability for plant uptake from the agricultural soils (Alewell et al., 2020).
Rhizosphere-dwelling phosphate-solubilizing bacteria (PSB) convert insoluble phosphates into bioavailable forms facilitating P uptake by plants via secretion of organic acids (Chungopast et al., 2021). The P solubilizing efficacy of PSB depends on the nature and the amount of organic acid production by these bacteria (Liang et al., 2020). Among different acids produced by PSB, gluconic acid is produced predominately via a direct oxidation pathway of glucose catalyzed by periplasmic glucose dehydrogenase (GDH). This GDH enzyme is encoded by the gcd gene and a redox coenzyme pyrroquinolinequinone (PQQ) encoded by pqq operon comprising six core genes, pqq A, B, C, D, F, and G (An and Moe, 2016). Another enzyme produced by PSB involved in organic P mineralization is phytase encoded by phy genes. This enzyme is responsible for P release in soil from organic materials stored in the form of phytate (Richardson and Simpson, 2011).
Furthermore, plant growth-promoting rhizobacteria (PGPR) may indirectly stimulate plant growth and help the plants to alleviate oxidative stress by enhanced production of antioxidant enzymes [e.g., catalase (CAT), peroxidase (POD), and superoxide dismutase (SOD)] in plant tissues (Batool et al., 2020;Bhat et al., 2020;Ha-Tran et al., 2021). Abiotic factors can cause oxidative stress due to an overproduction of reactive oxygen species (ROS). Oxidative product levels like H 2 O 2 and malondialdehyde (MDA) are the primary marker of ROS and lipid peroxidation responses (Thirupathi et al., 2020). Antioxidant enzymes are of paramount importance to counteract harmful ROS in plant tissues, maintain homeostasis within the cell, and enhance the tolerance of the plant (Dumanović et al., 2020).
With the growing importance as a promising strategy for sustainable agriculture, the mechanism of P solubilization of these PSB-based biostimulants needs to be deeply explored for optimal utilization of these microorganisms under variable field conditions (Alori et al., 2017;Bargaz et al., 2018;Elhaissoufi et al., 2020). It is still difficult to point out a single mechanism responsible for P solubilization by PSB (Sharma et al., 2013). Plant traits, i.e., extensive rooting system and root exudates, can contribute to a greater extent for P uptake. Root architecture is the spatial configuration of the root system in the soil, essential for plant P acquisition (Wissuwa et al., 2020), but root architecture is highly plastic in its response to low P conditions of the soil (Peris et al., 2020). Soil microbes play a key role to build specific patterns of metabolites around the roots by preferentially stimulating root exudation of primary metabolites, i.e., sugars, amino acids, and organic acids. Therefore, soil-dwelling microbes exert a strong influence on plant nutrient sensing (Canarini et al., 2019). Most of the previous studies discussed the direct effects of PSB on plant growth and crop yield (Singh et al., 2018;Liu et al., 2019). Other studies reported beneficial effects of PSB on root parameters (Suleman et al., 2018;Rasul et al., 2019;Rezakhani et al., 2019) showing no stronger linage of root exudation, root architectural variation, and antioxidant enzymes to supplement the plant in response to PSB inoculation. Unfortunately, the relative contribution of PSB to P acquisition remains largely unknown in this symbiotic relationship. Hence, it is imperative to better understand the mechanism of plant-microbe interaction in plant P acquisition under P-limited conditions to reduce reliance on chemical fertilizers. Therefore, the current study was focused to evaluate the role of native PSB from the wheat rhizosphere of unexplored agroecological zones in the establishment of a mutualistic relationship with plants for P acquisition. The experiments were conducted in hydroponics (without the interference of soil) as well as in soil under P-deficient conditions. We hypothesized that (H1) P-solubilizing bacteria enhance nutrient (P) availability, with an increase in root exudations in response to mutualistic association with a plant. We assumed (H2) that PSB inoculation may enhance the activity of antioxidant enzymes and improve root architecture that might help in the adaptation of plants under P-deficiency, and (H3) such PSB-mediated morphological and physiological changes in the plant may improve the grain yield of wheat.

Sample Collection and Soil Physiochemical Analysis
Rhizosphere soil samples were collected from 31 different sites of eight major wheat-growing agro-ecological zones of Pakistan Figure S1). At each sampling site, soils were collected randomly from five different locations at the depth of 0-20 cm, homogenized and stored at 4 • C prior to use. A sub-sample of soil (0.5 kg) was used for the determination of physicochemical characteristics (Supplementary Table S1). Soil pH was determined using pH meter (PHS-3C, REX, Shanghai, China; Thomas, 1996). The electric conductivity (EC) of soil was measured using a conductivity meter (DDS-307A, REX; Rhoades, 1993). Soil organic matter was determined using the wet oxidation method (Nelson and Sommers, 1996). Total N was measured using the Kjeldal method (Bremner and Mulvaney, 1982). Soil available P was determined by the sodium bicarbonate (NaHCO 3 ) method (Olsen et al., 1954). Exchangeable K and Na were extracted in 1 N ammonium acetate (NH 4 CH 3 CO 2 ) solution and estimated by a flame photometer (Model 410, Corning, Halstead, United Kingdom;Simard, 1993).

Isolation and Identification of P-Solubilizing Bacteria
Phosphate solubilizing bacteria were isolated on NBRIP (National Botanical Research Institute's Phosphate) (Nautiyal, 1999) and Pikovskaya's (Pikovskaya, 1948) agar media containing TCP (tricalcium phosphate) as the sole P source. The plates inoculated with soil suspension were incubated at 28 • C ± 2 • C for 7 days. The colonies with clear halo zone were considered as phosphate-solubilizing isolates and further purified by streaking on NBRIP and Pikovskaya's agar plates to obtain single colonies. Solubilization index was determined using the formula described by Paul and Sinha (2017).
A biosafety assessment of P-solubilizing isolates was carried out on commercially available blood agar medium. Single colony of bacteria was streaked on blood agar plates and kept at 37 • C ± 2 • C for 24 h (Russell et al., 2006). Cell and colony morphology, cell motility, and Gram's staining of PSB were studied using light microscopy (Gopala Krishnamurthy et al., 1967).
Total genomic DNA of each bacterium was obtained using the CTAB method (Wilson, 2001) and was used as a template for amplification of 16S rRNA gene using universal primers (forward PA 5 -AGACTTTGATCCTGCTCAG-3 and reverse PH 5 -AGGAGGTGATCCAGCCGCA-3 ) according to Ayyaz et al. (2016). Amplified 16S rRNA PCR products were analyzed on 1% agarose gel and further purified by PCR purification kit (QIAGEN Sciences, Germantown, MD, United States) prior to sequencing (Macrogen Inc., Seoul, South Korea). The derived sequences were compared to 16S rRNA gene sequences available at NCBI Gen Bank database using the BLAST algorithm for bacterial identification. The nucleotide sequences were aligned, and a phylogenetic tree was constructed using the maximum likelihood method (Kumar et al., 2016). Sequences were deposited to NCBI GenBank 1 .

Quantitative P Solubilization and Production of Organic Acids
Quantification of soluble P and organic acids released by PSB was carried out in an in vitro assay at 7 DPI (days post-inoculation). A single purified colony of each isolate was inoculated in NBRIP and Pikovskaya's broth media. Inoculated broth cultures were incubated at 28 • C ± 2 • C and 180 rpm on the shaker. Cultures were subjected to centrifugation at 4,000 rpm for 10 min at 4 • C to get cell-free supernatants. The available P in culturesupernatant was spectrophotometrically estimated at 880 nm using the molybdenum blue method (Murphy and Riley, 1962).
For analysis of organic acids, NBRIP filtrates were quantified on high-performance liquid chromatography (HPLC; Agilent 1,200 Series; Agilent Technologies, Santa Clara, CA, United States) equipped with C-18 column using mobile phase methanol: phosphate buffer (90: 10 v/v; pH 2.7) at flow rate of 1 ml min −1 and monitored at 210 nm (Tahir et al., 2013). Organic acids, i.e., acetic, citric, gluconic, malic, and succinic acid, were analyzed. Peak area and retention time were compared to standards (Sigma) for the quantification of organic acids (Vyas and Gulati, 2009).

Plant Growth-Promoting Traits of PSB
Plant growth-promoting traits of PSB were evaluated by using standard protocols. Indole acetic acid (IAA) production by PSB was quantified by using Salkowski's method (Gordon and Weber, 1951) on HPLC (Tien et al., 1979). Gibberellic acid produced by PSB was quantified according to protocol by Tien et al. (1979).
Siderophore production by PSB was detected using chrome auzurol S (CAS) agar medium (Schwyn and Neilands, 1987). The development of pink coloration around the bacterial colonies is an indication of siderophore production. Zinc solubilization by PSB was determined on Tris-minimal salts medium supplemented with insoluble zinc oxide (ZnO: 14 mM). Formation of the halo zone around bacterial colonies showed zinc solubilization. Solubilization index was measured as described by Fasim et al. (2002).

Effect of PSB Inoculation on Wheat Germination and Seedling Vigor
Based on in vitro P solubilizing efficacy, seven P-solubilizing isolates (ZW9, ZW32, SSR, D1, S1, TJA, and TAYB) were selected to study their effect on seed vigor index, germination, and root morphological parameters. Seeds of wheat variety "Faisalabad-08" were surface-sterilized for 5 min with 1.5% sodium hypochlorite solution followed by successive washings with sterile water to remove disinfectant traces. Sterilized seeds were inoculated with bacterial cultures (1 × 10 9 CFU ml −1 ) separately for 30 min. Uninoculated seeds soaked in LB broth were used as controls. Seeds were kept on water agar (0.25%) containing petri dishes (9 cm diameter × 2 cm depth). The plates were incubated in dark for 72 h at 28 • C ± 2 • C followed by 12-h dark and light cycles for 7 days in a growth room. Six biological replicates and 15 seeds per plate for each treatment were arranged in completely randomized design (CRD). The experiment was conducted twice, and different plant growth parameters were studied at 7 days post inoculation (DPI).

Development of PSB Consortia
Seven efficient PSB (ZW9, ZW32, SSR, D1, S1, TJA, and TAYB) selected for consortia development were studied for in vitro compatibility test (Irabor and Mmbaga, 2017). Bacterial inoculum was spotted on nutrient agar plate spread with the tested strain. The plates were kept at 28 • C ± 2 • C and observed up to 72 h. No inhibition zone was observed around any bacterial colony, indicating that they were compatible with each other.
A single colony of each bacterium grown separately in LB medium at 28 • C ± 2 • C for 24-48 h. An equal volume of each bacterium grown in LB was mixed to develop a consortium (1 × 10 9 CFU ml −1 ). Seeds were surface-sterilized with sodium hypochlorite solution (1.5%) for 5 min followed by successive washings (six times) with sterile water to remove disinfectant traces. Sterilized seeds of wheat variety 1, 2, and 3 were soaked in consortium-1, consortium-2, and consortium-3 suspensions (1 × 10 9 CFU ml −1 ); respectively for 30 min. Uninoculated seeds soaked in sterilized LB broth were used as controls.

Hydroponic Experiment to Study the Effect of PSB on Root Exudation and Anti-oxidant Enzymes
A hydroponic experiment was conducted to study the inoculation effects on wheat root exudation and anti-oxidant production in P-deficient conditions under greenhouse at NIBGE, Faisalabad (31 • 25 0 N 73 • 5 28 E). Surface-sterilized seeds were inoculated as described in the previous section. Inoculated seeds of variety 1, 2, and 3 and their respective uninoculated seeds were germinated in moistened sand with relative humidity 60-70% and day/night temperature of 25 • C /23 • C in the green house (Figure 1). After 1 week, seedlings were transferred to experiment buckets (24 cm × 37 cm) each having a volume of 4.5 L filled with half-strength Hoagland solution (Hoagland and Arnon, 1950) for proper plant growth and development (Islam et al., 2008).
Subsequently, 2-week-old plants were transferred to an aerated hydroponic system comprised of the opaque plastic box (13 cm × 13 cm) each having a volume of 1.5 L, filled with treatment-wise Hoagland solutions (Hoagland and Arnon, 1950). Three biological treatments, i.e., T1: Consortia-inoculated plants supplemented with modified hoagland containing 1 g L −1 tricalcium phosphate (TCP) as P-deficient solution, T2: Positive control (uninoculated plants supplemented with Hoagland containing all nutrients including P), and T3: Negative control (uninoculated plants supplemented with modified Hoagland containing 1 g L −1 TCP as P-deficient solution) were arranged in a two-factorial completely randomized design (CRD). Each box was aerated with an adjustable air pump that operates continuously to ensure a sufficient supply of air. Five seedlings per box were grown under a daily temperature regime of 25 • C ± 2 • C per day and 18 • C ± 2 • C per night with a photoperiod of 12-h light and darkness. Hoagland solution from each treatment was changed after every seventh day (Somasegaran and Hoben, 2012).
The plants were harvested 45 days after cultivation in a hydroponic system, and data were collected to study growth parameters of the plants, i.e., fresh and dry weight, root-shoot length, and plant P. Root growth parameters including root length (cm), root volume (cm 3 ), area of projection (cm 2 ), length per volume (cm m −3 ), root surface area (cm 2 ), root diameter (cm 2 ), tips, forks, and crossings were studied using rhizoscanner (Epson photo scanner V700) equipped with software WinRHIZO (Regent Int. Co., Ltd.).

Collection and Analysis of Root Exudates
Wheat plants were harvested after 45 days of cultivation in hydroponics, and root exudates were extracted in an Erlenmeyer glass flask by immersing the plant roots in sterile 2 mM CaCl 2 solution (Vranova et al., 2013). Flasks were wrapped with aluminum foil and placed in sunlight for 4 h. Root exudates collected in CaCl 2 solution were filtered using a 0.2-µm spinfilter (Thermo Fisher Scientific, Waltham, MA, United States) and stored at −20 • C for composition analysis (Yang et al., 2006).

Determination of Amino Acids
Aliquots (10 ml) of root exudates were lyophilized, cleaned with cold acetone, and kept overnight at −20 • C to remove protein interference. The supernatant was transferred to new tubes, vacuum-dried, and filtered through 0.2-µm spin-filter (Thermo Fisher Scientific). To derivatize amino acids, 5 µl OPA-3-MPA (O-pthaladehyde 3-mercaptopropionic acid) and 7 ml borate buffer (0.4 M) were added to a 1-ml sample and kept at room temperature for 1 min prior to run. Amino acids were separated in reversed-phase HPLC (Agilent 1200 Series) supplemented with a C-18 guard column at 1 ml min −1 flow rate and detected (λexc 250 nm; λem 395 nm) with a fluorescence detector (Frank and Powers, 2007).

Determination of Sugars
Aliquots (10 ml) of root exudates were lyophilized, re-dissolved in 40 µl H 2 O, and filtered through a 0.2-µm spin-filter (Thermo Fisher Scientific). Sugars were separated by HPLC under isocratic conditions (acetonitrile and water in 75: 25 v/v) at 1 ml min −1 flow rate (Lugtenberg et al., 1999) and detected at 487 nm using a UV detector.

Determination of Organic Acids
Organic acids in root exudates were quantified by following the procedure as described by Kamilova et al. (2006). Organic acids were separated using mobile phase 10 mM H 3 PO 4 at 0.8 ml min −1 flow rate and detected at 210 nm using HPLC with a UV detector (Agilent 1200 Series).

Determination of Antioxidant Enzyme Activity
Fresh leaves (1 g) of wheat plants harvested from the hydroponic system were homogenized in ice-cold 10 ml of potassium phosphate buffer (100 Mm; pH 7) for extraction of enzymatic anti-oxidants. The homogenate was centrifuged (10,000 × g) for 20 min at 4 • C and stored at −80 • C for analysis of antioxidant enzymes.

Superoxide Dismutase (SOD)
SOD activity was determined by following the process described by Dhindsa and Matowe (1981). The reaction mixture contained 50 µl Nitro blue tetrazolium chloride (NBT), 100 µl L-methionine, 50 µl riboflavin, 250 µl sodium phosphate buffer, and 50 µl enzyme extract. The reaction mixture was placed under a fluorescent lamp (3 W) for 15 min followed by exposure to dark for 15 min and measured the absorbance at 560 nm using spectrophotometer (CamSpec M550 double beam UV-Vis; Spectronic Camspek Ltd, Leeds, United Kingdom). One unit of SOD activity is the amount of enzyme required for 50% inhibition of photo-chemical reduction of NBT.

Catalase (CAT)
CAT activity was measured as described by Aebi (1984). The reaction mixture contained 1.9 ml phosphate buffer (50 mM; pH 7) and the enzyme-extract (100 µl). To initiate the reaction, 1 ml H 2 O 2 was added to the reaction mixture the absorbance was measured at 240 nm after every 30 s until 3 min using spectrophotometer. One unit of CAT activity is defined as absorbance change of 0.01 U min −1 .

Peroxidase (POD)
POD activity was measured by mixing 50 mM phosphate buffer (pH 7), 20 mM of guaiacol, and 100 µl of enzyme extract (Rao et al., 1996). To initiate the reaction, 40 mM H 2 O 2 was added to the reaction mixture and incubated at 25 • C for 5 min. Change in the absorbance at 436 nm was monitored every 2 min. One unit of POD activity is an absorbance change of 0.01 U min −1 . All antioxidants are expressed as unit mg −1 fresh protein.

Malondialdehyde (MDA) and H 2 O 2
Change in lipid peroxidation was studied by calculating MDA production in leaves of wheat plants harvested from the hydroponic system to determine the level of oxidative damage caused by P deficiency. Fresh leaves (1 g) were homogenized in 20 ml trichloroacetic acid (TCA; 0.1 percent) and centrifuged (12,000 × g) for 20 min at 4 • C for MDA and H 2 O 2 contents estimation.
For MDA extraction aliquot of supernatant (1 ml) was mixed with 4 ml TCA containing thiobarbituric acid (TBA; 5%). The reaction mixture was incubated at 95 • C for 15 min and subsequently cooled on ice-bath to stop the reaction. Samples were centrifuged at 12000 × g for 10 min, and absorbance was recorded at 532 nm and 660 nm (Demiral and Türkan, 2005).
H 2 O 2 content was measured according to Velikova et al. (2000). The reaction mixture (3 ml) contains 0.5 ml supernatant, 10 mM potassium phosphate buffer (pH 7), and 1 M KI (potassium iodide) solution. The mixture was vortexed, and absorbance was monitored at 390 nm. Amount of H 2 O 2 content was calculated using a standard curve prepared with known H 2 O 2 concentrations.

Pot Experiment Using Soils of Different Agro-Ecological Zones
Effects of P-solubilizing consortium-1, consortium-2, and consortium-3 on recommended wheat varieties of Province 1, 2, and 3 were further assessed in a pot experiment at NIBGE (31 • 25 0 N 73 • 5 28 E) greenhouse using native soils collected from eight different agro-ecological zones of Pakistan. Among 31 sites of isolation (Supplementary Figure S1 and Supplementary Table S1), soils from 16 major wheat-growing sites were selected for evaluation of PSB consortia (Supplementary Table S5).
Consortium-1 was used for wheat variety 1 (Faisalabad-08) recommended for sites (site no. 1-9, Supplementary Table S5) of Province 1. Consortia-2 was used for wheat variety 2 (Fakhr-e-Sarhad) recommended for sites (site no.  Table S5). Seeds were surface-sterilized with sodium hypochlorite solution (1.5%) as described in previous sections. Sterilized seeds were inoculated with respective PSB consortia (1 × 10 9 CFU ml −1 ) for 30 min. Uninoculated seeds dipped in LB medium were used as controls. Sowing was done in pots (10 cm diameter) containing 300 g soil per pot (supplemented with TCP @ 1 g/pot) from 16 different wheat-growing areas of three provinces. There were 32 treatments (16 different soils with and without inoculation) and six biological replicates arranged in a two-factor completely randomized design (CRD; Supplementary Table S5). For proper seed germination, all pots were watered before sowing. Each pot has two plants. After seed germination, pots were kept moist as per requirement by providing them with water and nutrient solution (Hoagland and Arnon, 1950) without P source, alternatively.
Plants were uprooted at 45 DAS and evaluated for the growth parameters, i.e., plant fresh weight, dry weight, root length, and shoot length. Ground plant material (1 g) from each treatment was used for plant P analysis by using the tri-acid digestion method as described by Tandon (1993). Rhizospheric soils from each treatment were analyzed for viable count (Somasegaran and Hoben, 2012). Soil available P was measured by the sodium bicarbonate (NaHCO 3 ) method (Olsen et al., 1954), and soil phosphatase activity was measured by the p-nitophenyl method as described by Tabatabai and Bremner (1969).

Evaluation of PSB Consortium for Wheat Yield Parameters
Based on the higher P-solubilizing potential of consortium-1 under controlled conditions, a pot experiment was carried out to evaluate the effects of consortium-1 on yield and P content of wheat under net house conditions at NIBGE, Faisalabad (31 • 23 45.1 N, 73 • 01 3.4 E) during the wheat season November 2018-April 2019.
Seeds were pelleted by mixing press mud as carrier material (@ 50 kg seed per kg carrier material) with consortium-1 (1 × 10 9 CFU ml −1 ) and left for 1 h. Uninoculated pelleted seeds were used as a control. Press mud-based carrier material had N 4.3%, P 2%, organic-matter 65%, pH 7.5, EC 3 ms cm −1 , 40% moisture content, and 2 mm particle size. Sowing was done in earthen pots (diameter: 30 cm) filled with 12 kg soil (Loam texture, organic-matter 0.57%, pH 8.2, and available P: 1.9 mg kg −1 ) and watered for proper seed germination before sowing. There were three treatments including consortium-1 inoculated seeds supplemented with 80% recommended dose of DAP (i.e., 20% reduced DAP) and two uninoculated controls supplemented with 80 and 100% DAP (i.e., recommended dose of DAP), respectively. Six biological replicates and six plants per pot were arranged in two-factor CRD.
Recommended doses of fertilizers were applied (N: P, 15:100 kg ha −1 ). P was added in the form of DAP at the time of sowing, while N (urea) was added in split doses, i.e., initially at the time of sowing and subsequently with first and second irrigations. Plants were uprooted at 35 DAS and the growth parameters, i.e., root length, shoot length, and plant dry weight were recorded. At maturity, plants were harvested and data for grain yield, plant biomass, plant height, number of tillers, and plant P were recorded (Tandon, 1993). Rhizosphere soil from each treatment was analyzed for available P (Olsen et al., 1954) and phosphatase activity (Tabatabai and Bremner, 1969).

Detection of Inoculated PSB
The rhizosphere soil of wheat plants grown in pots was also evaluated for the survival of inoculated PSB using viable count (Somasegaran and Hoben, 2012) and BOX-PCR (Basheer et al., 2016). Re-isolated PSB was identified by comparing the morphological characteristic of inoculated bacteria and other plant growth-promoting attributes like P solubilization, IAA production, and zinc solubilization . Strain-specific fingerprints of reisolated PSB were compared to those of pure colonies using BOX-A1R primer 5 CTACGGCAAGGCGACGCTGACG 3 (Basheer et al., 2016).

Statistical Analysis
All data from in vitro studies, HPLC results, and pot experiments were analyzed statistically by ANOVA. The variation between treatments was compared by least significant difference (LSD) at 1 and 5% level of confidence for lab and net house experiments using Statistix 10 software (Analytical Software, Tallahassee, FL, United States). Principal component analysis (PCA) for various root and yield parameters and regression analysis were carried out using SPSS 23.0 software (SPSS Inc., Chicago, IL, United States). Box plots were applied for analysis of soil parameters in a pot experiment using different soils according to ANOVA and Tukey's HSD test at p < 0.05) using Origin Software Package Version 2020b (OrginLab Corporation, Northampton, MA, United States).
Among 25 selected PSB, Ochrobactrum sp. SSR (SI: 5.8), Enterobacter sp. D1 (SI: 3.8), Pantoea sp. S1 (SI: 4), Pseudomonas sp. TJA (SI: 3.6), and Bacillus sp. TAYB showed a higher P FIGURE 2 | Mineral phosphate solubilization by PSB isolated from different wheat-growing agro-ecological zones of Pakistan. P solubilizing index was calculated from halo zone formation on NBRIP and Pikovskaya's agar medium (A). Quantification of P solubilized by PSB in liquid NBRIP and Pikovskaya's medium (B). pH of medium (initial pH 7) at seventh day of pot inoculation (C). There were six biological replicates, and standard error of the means is represented as bars. The significant difference (p < 0.01) among treatments is denoted as means and represented by different letters (a, b, c).*PSB showing significantly higher P solubilizing index and P solubilized, while significantly lowering the pH of the NBRIP and Pikovskaya's medium. ZW9, ZW32, and D1 are Enterobacter spp.; SSR: Ochrobactrum sp.; S1: Pantoea sp.; TJA: Pseudomonas sp.; and TAYB: Bacillus sp.

Quantitative P Solubilization and Production of Organic Acids
All PSB strains showed P solubilization activity in both NBRIP and Pikovskaya's broth along with a subsequent decrease in pH ( Figure 2C). The solubilized P ranged from 69-365 µg ml −1 in NBRIP medium with decline in pH up to 3.6, and from 43 to 278 µg ml −1 P in Pikovskaya's medium with pH decline up to 4, respectively. Maximum release of P (365 µg ml −1 ) was observed for Ochrobactrum sp. SSR in NBRIP medium ( Figure 2B). Primary selection of seven PSB isolates (Enterobacter sp. ZW9, Enterobacter sp. ZW32, Ochrobactrum sp. SSR, Enterobacter sp. D1, Pantoea sp. S1, Pseudomonas sp. TJA, and Bacillus sp. TAYB) was carried out based on their P-solubilizing potential as indicated by in vitro qualitative and quantitative tests.

Hydroponic Experiment to Study the Effect of PSB on Root Exudation and Anti-oxidant Enzymes
The effect of PSB on wheat root exudation and anti-oxidant production was studied in a hydroponic experiment in the greenhouse. Plants inoculated with selected PSB consortia had a positive impact on growth parameters of tested wheat varieties (variety 1, variety 2, and variety 3) compared to an uninoculated positive control (plants supplemented with Hoagland containing all nutrients including P) and uninoculated negative control Effect of PSB inoculation on root morphology parameters. Seeds were inoculated with PSB in water agar medium in a plate assay and incubated at room temperature for 7 days. All values are an average of six biological replicates. ± represents standard deviation. Means followed by the same letter differ non-significantly at p = 0.05 according to LSD.
(supplemented with modified Hoagland containing 1 g L −1 TCP as P-deficient solution; Figure 4). Inoculations (supplemented with modified Hoagland containing 1 g L −1 TCP as P-deficient solution) with respective bacterial consortia resulted in improved plant growth (shoot length and plant fresh and dry weight) and plant P content (1.3-1.6%) of three tested wheat varieties, compared to uninoculated controls (negative and positive; Supplementary Table S4).

Root Exudate Analysis
Effect of PSB inoculations on wheat root exudation was analyzed hydroponically in greenhouse. Major root exudate components characterized and quantified by HPLC showed that root samples from PSB-inoculated plants of three wheat varieties produced  a higher amount of amino acids, sugars, and organic acids as compared to positive and negative controls. Thirteen amino acids were released from the roots of each wheat variety with variable compositions. Ornithine was the most abundant amino acid (356 nmol g −1 DW 4 h −1 ) in inoculated roots of wheat variety 3 ( Figure 5G) followed by 320 nmol g −1 DW 4 h −1 in inoculated roots of wheat variety 1 ( Figure 5A) and 306 nmol g −1 DW 4 h −1 in inoculated roots of wheat variety 2 ( Figure 5D) as compared to uninoculated negative and positive controls. Root samples from positive controls produced a higher amount of amino acids as compared to negative controls. Serine (63 nmol g −1 DW 4 h −1 ) and isoleucine (48 nmol g −1 DW 4 h −1 ) were the next most abundant amino acids, while citrulline, tyrosine, and valine were the least abundant amino acids released from roots of wheat variety 3. In the case of wheat variety 2, cysteine (59 nmol g −1 DW 4 h −1 ) was the second most abundant amino acid produced after ornithine, while asparagine and valine (2 nmol g −1 DW 4 h −1 ) were the least abundant amino acids. Leucine (110 nmol g −1 DW 4 h −1 ) was the most abundant, while valine and tyrosine were the least abundant amino acids produced by wheat variety 1. Six sugars (glucose, sucrose, fructose, xylose, mannose, and galactose) were released in variable quantities from roots of all three wheat varieties. Glucose was the most abundantly produced sugar in inoculated plants of both variety 1 (250 nmol g −1 DW 4 h −1 ; Figure 5B) and variety 3 (180 nmol g −1 DW 4 h −1 ; Figure 5H), respectively, in comparison to negative and positive controls, while fructose was the most abundant (280 nmol g −1 DW 4 h −1 ) and mannose was the least abundant (30 nmol g −1 DW 4 h −1 ) sugar produced in inoculated treatments of variety 2 as compared to uninoculated controls ( Figure 5E).
Principal component analysis (PCA) showed clustering of plant P content, root morphological traits (root surface area, root diameter, root length, and root tips), and root exudations (sugars, amino acids, and organic acids) of PSB-inoculated plants in one group indicating a positive correlation among these parameters (Figure 6).

Antioxidant Enzyme Analysis
Accumulated ROS-induced oxidative damage and MDA content in plant roots were reduced by PSB inoculations. All wheat varieties showed higher H 2 O 2 contents in negative controls as a result of induced P deficiency. The production of super dismutase (SOD), catalase (CAT), and peroxidase (POD) in PSB-inoculated treatments of all wheat varieties were higher as compared to negative and positive controls. In addition, CAT, SOD, and POD were the highest produced enzymes in wheat variety 2 followed by variety 1 and variety 3 (Figure 7).

Pot Experiment Using Soils of Different Agro-Ecological Zones
Effect of soil-specific consortia on recommended wheat varieties of three major wheat-growing provinces was studied in a pot experiment, nine soils (site nos. 1-9, Supplementary  Table 4).
The plant growth of wheat variety 1 inoculated with consortium-1 was significantly increased. Significant increase in plant shoot length (32 cm), root length (29.5 cm), and plant dry weight (1.05 g) of consortium-1-inoculated wheat (variety-1) were observed in soils of site 7 (Gujranwala) with 40% increase in root length as compared to uninoculated controls. Increased shoot length (28 cm), root length (27 cm), and dry weight (0.91 g) were observed in inoculated (consortium-2) wheat variety 2 grown in site 1 (Dir) soil of Province 2 with 21% increased root length as compared to uninoculated controls; whereas inoculation of plants (variety 3) with consortium-3 improved shoot length (29 cm), root length (33 cm), and plant dry weight (0.80 g) with a 48% increase in root length as compared to uninoculated controls in site 3 soil (Sanghar) of Province 3 ( Table 4).
Consortium-1, 2, and 3 significantly increased (up to 4%) plant P content in respective wheat varieties grown in soils of site 7 of province-1, site 1 of province-2, and site 3 of province -3 as compared to uninoculated plants grown in soils of respective sites (Figure 8). Soils of province 1 and province 3 are majorly comprised of irrigated plains and soils of province 2 prevalent with dry mountains. The inoculated plants with consortium-1 increased soil available P (5.2-9.8 µg g −1 soil) and phosphatase activity (24.3-29.6 µmol g −1 soil h −1 ) in inoculated soils of province-1. Consortium-2 increased soil available P up to 5-9 µg g −1 soil and phosphatase activity 25-30 µmol g −1 soil h −1 . Consortium-3 also increased available P up to 4-6.5 µg g −1 soil and phosphatase activity (17-20 µmol g −1 soil h −1 ; Figure 8 and Supplementary Table S6). A positive correlation was found as a result of inoculation of three PSB consortia on plant growth (plant fresh and dry weight, root-shoot length, and plant P), soil available P, and soil phosphatase activity for recommended wheat varieties grown in native soils of Province 1, Province 2, and Province 3 (Supplementary Figure S2).

Evaluation of PSB Consortium for Wheat Yield Parameters
P-solubilizing consortia significantly enhanced wheat growth in a pot experiment under net house conditions. PSB consortium-1 inoculated wheat variety1 showed increased grain yield (5.95 g plant −1 ) in comparison to 80% uninoculated control (80% of recommended DAP dose). A significant increase (4.1%) in seed P was observed in inoculated plants compared to 80% uninoculated control. An increase in soil available P (6.4 µg g −1 soil) and phosphatase activity (24 µmol g −1 soil h −1 ) was observed with PSB inoculation ( Table 5). A positive correlation was found between wheat yield parameters, i.e., seed P content, wheat grain yield, soil available P, and gluconic acid produced by P-solubilizing bacteria in an in vitro quantification assay (Supplementary Figure S3).

Detection of Inoculated PSB
To confirm the presence of inoculated PSB in the rhizosphere of wheat variety-1 grown in a pot experiment under net house conditions, bacterial population was recorded using viable count at 35 DAS. The re-isolated colonies of Enterobacter spp. ZW9, ZW32, and Ochrobactrum sp. SSR were identified based on the morphological characteristics and other plant growthpromoting attributes, i.e., P solubilization (224-345 µg ml −1 , IAA production (4-47 µg ml −1 ), and zinc solubilization (SI: 2.9-3.9) in comparison to pure cultures. BOX-PCR confirmed the inoculated P-solubilizing bacteria. BOX fingerprinting of re-isolated colonies were found identical to those of pure cultures (Figure 9).

DISCUSSION
Due to the depleting reservoirs of available P and high costs of phosphatic fertilizers, PSB-based biostimulants serve as an alternative and sustainable strategy in modern agriculture for food and environmental security. Therefore, the present study was designed to get comprehensive data of wheat rhizosphereassociated PSB from different soils belonging to major wheatgrowing agro-ecological zones of Pakistan (situated between 24 and 37 north and longitude of 61 to 75 east; Supplementary Table S1 and Supplementary Figure S1). We determined the beneficial effects of native PSB inoculum on wheat growth, yield, and P uptake. Previous studies reported the presence of PSB in the soils of selected areas of Pakistan (Suleman et al., 2018), but major wheat growing agro-ecological zones remained unexplored.
Ochrobactrum (SSR), Pantoea (S1), and Pseudomonas (TJA). Increase in soluble P concentration was observed in quantitative assay with a concomitant decrease in pH up to 3.6 from initial pH 7 in liquid medium (Figure 2). We observed that pH of the medium decreased as soluble P increased. This might be because P-solubilizing bacteria can secrete organic acids in a culture medium to degrade TCP (insoluble P source) in order to increase soluble P (Penn and Camberato, 2019;Chouyia et al., 2020). Among organic acids detected by HPLC acetic acid, gluconic and citric acids were predominantly produced by selected P-solubilizing bacteria ( Table 1). Three most efficient PSB strains, Ochrobactrum sp. SSR and Enterobacter spp. (ZW9, ZW32) produced higher gluconic acid and acetic acid. This indicates that apart from gluconic acid, other acids like acetic acid and citric acid are produced in significant levels by bacteria. Suleman et al. (2018) reported gluconic acid and acetic acid as dominantly produced acids by Pseudomonas sp. and Enterobacter sp., respectively. Differential production of these organic acids by PSB seems a key factor responsible for phosphate solubilization (Azaroual et al., 2020;Tian et al., 2021).
In order to link bacterial genetic potential with P-solubilizing activity, glucose dehydrogenase (gcd), pqqE, and phytase (phy) genes (Perez et al., 2007;Han et al., 2008;Chen et al., 2016) related to the organic and inorganic P solubilization were studied in selected PSB. Gcd was amplified in Pantoea sp. S1, while phy genes were amplified in both S1and Enterobacter sp. ZW32. The current study reported amplification of PqqE gene in ZW32 and SSR (Figure 3 and Supplementary  Table S3). In our previous study (Rasul et al., 2021), expression analysis of multiple genes related to P-release pathways of PSB was carried out using qPCR. This study revealed that Ochrobactrum sp. SSR having gcd and phy genes can solubilize different organic and inorganic P sources (phytate, rock phosphate, and K 2 HPO 4 ). SSR was found with multifaceted nature; hence, it was selected as a common P-solubilizing strain in consortia development. Efficient and compatible PSB (Enterobacter spp. ZW9, ZW32 and D1, Ochrobactrum sp. SSR, Pantoea sp. S1, Pseudomonas sp. TJA, and Bacillus sp. TAYB) with multiple plant growth-promoting (PGP) traits (i.e., IAA, siderophores, and zinc solubilization) was used to develop three consortia for their native soils and respective recommended wheat varieties. This is due to indigenous microbes that are most likely to promote plant growth as they are more adaptive with higher persistence in native soil (Glick, 2012;Souza et al., 2015).
Root exudation is an important process that regulates the interaction between root and soil microbes (Mommer et al., 2016). These metabolites are difficult to collect from real soil environments due to their chemical interaction with soil matrix and background microbial exudation (Olanrewaju et al., 2019;Kumar and Dubey, 2020;Pang et al., 2021). Therefore, in the present study, the hydroponic system was used to evaluate the effect of PSB consortia on root exudation with consequent changes in root architecture. Characterization and HPLC quantification of major root exudates from hydroponically grown wheat varieties showed that PSB-inoculated roots produced higher amounts of amino acids, organic acids, and sugars under P-deficient conditions as compared to uninoculated controls ( Figure 5). It has been known that soil microbes can release certain compounds to stimulate the root exudation of primary metabolites (i.e., amino acids, organic acids, and sugars) that are passively lost from roots and utilized by rhizosphere-dwelling microbes (Canarini et al., 2019).
In the present study, 13 amino acids and six sugars were released from roots of three wheat varieties with variable compositions. Among amino acids, ornithine was abundantly found in root exudation of PSB-inoculated plants of three wheat varieties as compared to uninoculated plants; cumulative release of sugars in root exudates was higher in variety 1 and variety 2 compared to variety 3 (Figure 5). Amino acids in root exudates besides sugars tailor communities because they represent carbon and nitrogen substrates for microbial growth (Pongrac et al., 2020). Higher release of these exudates contributes to increase in bacterial persistence in the root rhizosphere (Nwachukwu et al., 2021). In return, bacteria can benefit plants by enhanced the nutrient absorbing capability of roots especially to facilitate phosphorus uptake (Grover et al., 2021). Furthermore, HPLC analysis revealed a significant FIGURE 8 | Effect of PSB inoculation on soil and plant P parameters of wheat varieties grown in soils from different agro-ecological zones compared to uninoculated controls. Evaluation of consortium-1 for soil available P (A), soil phosphatase activity (B), and plant P content (C) of wheat variety-1 grown in soils of province-1. Evaluation of PSB consortium-2 for soil available P (D), soil phosphatase activity (E), and plant P content (F) of wheat variety-2 grown in soils of province-2. Evaluation of consortium-3 for soil available P (G), soil phosphatase activity (H), and plant P content (I) of wheat variety-3 in soils of province-3. Different small (a, b, c) letters above each box indicate significantly differences between the treatments according at p < 0.05. I: Inoculated, C: Uninoculated control. For Province 1: increase of organic acids specifically gluconic acid (20-39%), malic acid (30%), and citric acid (21%-41%) released in exudates collected from PSB-inoculated three wheat varieties as compared to their respective uninoculated controls (Figure 5). Increased concentrations of oxalic, malic and citric acids in root exudates act as chemo-attractants for soil microbes during beneficial interactions (Pascale et al., 2020;Tsai et al., 2020).
Root response to P plays a major role in the modulation of root architecture as compared to other nutrients (Peris et al., 2020). Therefore, in addition, to study the root exudates, rhizoscanning of hydroponically grown plants was also carried out to explore the effect of PSB on root architecture. PSB inoculation modified root architecture, as indicated by an increase in total root length (24-27%), root tips (46-53%), root surface area (29-35%), root volume (36-48%), and area of projection (28-47%) of three wheat varieties provided with an insoluble source of P (TCP) as compared to uninoculated control. Significant increases in root diameter, root tips, and root surface area were key parameters modified in response to PSB inoculation ( Table 3). Effect of bacterial inoculation on wheat growth and yield parameters, P uptake and soil nutrients at 35 days after sowing (DAS) and at 120 DAS in a pot experiment at Faisalabad grown wheat during winter season 2018-19. Control 80% and Control 100% represent non-inoculated controls supplemented with 80% of the recommended dose of P fertilizer and full recommended dose of P fertilizer, respectively. PSB (1 × 10 9 CFU ml −1 ) were seed-inoculated before sowing. Inoculated pots were supplemented with 80% of recommended dose of DAP fertilizer. Means are the average of six biological replicates arranged in CRD. Means followed by the same letter differ non-significantly at p = 0.05 according to LSD. ± Represents the standard deviations (SD). Six plants were grown in earthen pots (diameter: 12 in, height: 14 in) containing 12 kg soil.
The root tip is the first part of the plant to explore a new soil environment and plays a crucial role in response to environmental stimiuli. A growing number of studies showed root tip functions as a control center for sensing externalnutrient concentration and translating it into an alteration of root architecture (Huang and Zhang, 2020;Grover et al., 2021). Increased root diameter promotes plant growth and development since thicker roots easily penetrate the soil, thereby anchoring the plant more effectively leading to efficient nutrient uptake (Freschet et al., 2020). In addition to improved architecture, a significant increase in growth and plant P contents of three wheat varieties was observed in plants inoculated with PSB consortia as compared to the uninoculated controls (Supplementary Table S4). Principal component analysis showed clustering of P content, root morphological traits, and root exudations of PSBinoculated plants in one group, indicating a positive correlation between these parameters (Figure 6). PSB-stimulated primary metabolites may serve important functions in the soils as a key nutrient source for rhizosphere microbes and influence rootmicrobe relations (Pang et al., 2021).
Besides the effect of PSB on root architecture and exudation in hydroponics, these bacteria sustained the plant growth under P stress as indicated by a significant increase in plant antioxidant enzyme (SOD, POD, MDA, and CAT) activity with a concomitant decrease in MDA and H 2 O 2 contents (Figure 7). There is increasing evidence that PGPR confers tolerance in plants to various abiotic stresses (Mokrani et al., 2020;Bhat et al., 2020;Khan et al., 2020;Camaille et al., 2021) by activating enzymes like CAT, POD, and SOD activity which act as primary reactive oxygen species (ROS) scavengers (Grover et al., 2021). These enzymes catalyze the O 2− to H 2 O 2 , and H 2 O 2 is further converted into H 2 O and O 2 by peroxidases, thereby enhancing tolerance and eliciting resistance in plants against oxidative stress (Li et al., 2020;Eljebbawi et al., 2021). P-solubilizing efficacy of three developed PSB consortia was further subjected to in planta evaluation under natural soil conditions using soils collected from different agro-ecological zones of three provinces with their recommended wheat varieties. Results showed a significant increase (up to 4%) in plant P contents and plant root length (up to 48%) in soils of province-1 and province-3 majorly comprising irrigated plains and soils of province 2 prevalent in dry mountains in response to inoculation of soil-specific consortia (Supplementary Table S1 and Supplementary Figure S1). A significant increase of soil available P and phosphatase activity was observed in PSBinoculated soils as compared to uninoculated controls (Figure 8). This is in line with Chen et al., 2021, who reported an elevated amount of available P in wheat-grown soil inoculated with PSB as compared to uninoculated control. Boubekri et al. (2021) reported the highest level of available P and greater stimulation of wheat plant growth in soils co-inoculated with phosphate-solubilizing consortium. Lack of region/soil-specific strains is one of the major constraints in the inconsistent performance of bio-inoculants. In this study, we tried to address this problem by developing region/soil-specific inoculum with expected outcome across wheat-cultivated agricultural lands. The use of region-specific strains ensures the optimal use FIGURE 9 | Re-isolation of inoculated PSB colonies. Photograph of gel indicates BOX PCR patterns of pure culture of PSB from wheat rhizosphere, 1: 1Kb DNA ladder; 2: SSR, 3: ZW32, 4: ZW9 (A). Photograph of gel indicates re-isolated colonies of PSB that are similar morphologically to the inoculated consortium-1 PSB. 1: 1Kb DNA ladder, 2: Re-isolated SSR colony, 3: Non-specific colony, 4: Re-isolated ZW32 colony, 5: Re-isolated ZW9, 6: Non-specific colony, 7: -ve Control; 8: 100 bp DNA ladder (B). Plate showing re-isolated colonies of SSR, ZW32, and ZW9 (C). NBRIP Plate showing solubilization zones formed by re-isolated colonies of SSR, ZW32, and ZW9 (D). Plant growth-promoting attributes of re-isolated PSB strains for confirmation of inoculated PSB: P solubilization and IAA production were detected quantitatively by HPLC, and Zinc solubilization was confirmed by formation of halo zone on agar plate added with ZnO (E). Values are an average of six biological replicates. Enterobacter spp. ZW9, ZW32, Ochrobactrum sp. SSR, -ve control represents negative control.
of biofertilizers with expected outcomes over vast areas of agricultural lands (Mitter et al., 2021). PCA showed a significant positive correlation between all plant growth parameters (root length, plant P, soil available P, and soil phosphatase activity) in PSB-inoculated wheat varieties as compared to their respective uninoculated plants (Supplementary Figure S2). The increased P availability in the rhizosphere soil and its better aboveground translocation occurred either directly or indirectly in response to PSB inoculation results in better plant growth (Elhaissoufi et al., 2020).
Furthermore, to study the contribution of efficient PSB consortium-1 to increase the crop yield, a pot experiment was carried out using wheat variety 1 under net house conditions at NIBGE. Inoculated plants provided with 20% reduced DAP application showed a significant increase (7%) in grain yield, plant biomass (20%), seed P content with an increase in soil phosphatase activity, and subsequent availability of soil P ( Table 5). Regression analysis confirmed a positive correlation between wheat yield parameters, i.e., seed P content, wheat grain yield, soil available P, and gluconic acid produced by P-solubilizing bacteria in an in vitro quantification assay (Supplementary Figure S3). Higher seed P content might be due to translocation of P to seed as previously reported by Feng et al. (2021). Enhanced soil phosphatase activity in PSB-inoculated treatments showed the ability of these PSB to convert organic P complexes in soil into an available form for plant uptake (Chawngthu et al., 2020;Chen et al., 2021). Previous studies demonstrated cost-effective use of PSB inoculations to increase grain yield and plant P in wheat, rice, and sugarcane along with reduced DAP dose (Suleman et al., 2018;Rasul et al., 2019;Rosa et al., 2020). These studies indicated that P-solubilizing microorganisms showed the best effect with reduced application of DAP fertilizers. The survival of inoculated strains of consortium-1 in root rhizosphere as studied by viable count and BOX-PCR revealed that these were rhizospherecompetent phosphobacteria. Moreover, plant growth-promoting traits of re-isolated bacteria in comparison with those of the pure culture of inoculated strain indicated the persistence of inoculated PSB (Figure 9).
The ecological approach used in the study allowed us to isolate and select multifaceted phosphate-solubilizing bacteria from different unexplored agro-ecological zones. These promising strains (ZW9, ZW32, and SSR), deposited to the DSMZ German culture collection (accession no. 109592-93 and 109610), serve as a global valuable pool. The phosphate solubilization activity of these potential PSB was highly promising in different types of soils and wheat varieties as statistically validated by improved plant growth, yield, and soil fertility-related traits like soil available P and soil phosphatase activity. This has lain a solid foundation for further utilization of these effective PSB candidates as soil/genotype-specific biofertilizers for wheat under a wide range of agro-ecological conditions intended to foster sustainable agriculture. The study is of prime importance because 40% of global P-deficient soil direly needs an eco-efficient solution for better cropping on a larger scale.

CONCLUSION
This study provides new evidence that P-solubilizing bacteria employed beneficial impact on morpho-physiological attributes of inoculated plants leading to alleviation of P stress through induced sequential production of root exudates, modification of root architecture, and mitigation of oxidative damage by induced activities of antioxidant enzymes. These PSB-induced structural and physiological alterations in the plant ultimately led to improved wheat growth and P acquisition.

AUTHOR CONTRIBUTIONS
MY analyzed the data and wrote the manuscript. EI designed the setup for hydroponics. MR and NM performed data analysis and reviewed the manuscript. LR provided soils from Sindh and reviewed the manuscript. IA helped in sequencing of both 16S rRNA and gene specific amplifications. MY, MR, LR, and IF executed statistical analysis. AT analyzed HPLC results of root exudation. MI helped in the analysis of physiochemical soil characteristics. SY conceived and supervised the whole study and edited the manuscript.