Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2023.1135693 Microbiology Original Research Cold-tolerant phosphate-solubilizing Pseudomonas strains promote wheat growth and yield by improving soil phosphorous (P) nutrition status Dasila Hemant 1 2 Sah V. K. 3 Jaggi Vandana 4 Kumar Arun 5 Tewari Lakshmi 2 Taj Gohar 6 Chaturvedi Sumit 3 Perveen Kahkashan 7 Bukhari Najat A. 7 Siang Tan Ching 8 * Sahgal Manvika 2 * 1Department of Microbiology, Akal College of Basic Sciences, Eternal University, Rajgarh, Himachal Pradesh, India 2Department of Microbiology, College of Basic Sciences and Humanities, G.B. Pant University of Agriculture and Technology, Pantnagar, India 3Department of Agronomy, College of Agriculture, G.B. Pant University of Agriculture and Technology, Pantnagar, India 4Department of Pediatrics and Human Development, Michigan State University, Grand Rapids, MI, United States 5Department of Agronomy, Dr. Khem Singh Gill, Akal College of Agriculture, Eternal University, Rajgarh, Himachal Pradesh, India 6Department of Molecular Biology and Biotechnology, College of Basic Sciences and Humanities, G.B. Pant University of Agriculture and Technology, Pantnagar, India 7Department of Botany and Microbiology, College of Science, King Saud University, Riyadh, Saudi Arabia 8School of Pharmacy, KPJ Healthcare University College, Nilai, Malaysia

Edited by: Deep Chandra Suyal, Independent Researcher, Bhopal, India

Reviewed by: Saurabh Kumar, ICAR-Research Complex for Eastern Region, India; Shailesh Pandey, Forest Research Institute (FRI), India; Damini Maithani, IFTM University, India

 *Correspondence: Manvika Sahgal sahgal.manvika@gmail.com Tan Ching Siang tcsiang@kpjuc.edu.my

This article was submitted to Terrestrial Microbiology, a section of the journal Frontiers in Microbiology

 13 03 2023 2023 14 1135693 01 01 2023 09 02 2023 Copyright © 2023 Dasila, Sah, Jaggi, Kumar, Tewari, Taj, Chaturvedi, Perveen, Bukhari, Siang and Sahgal. 2023 Dasila, Sah, Jaggi, Kumar, Tewari, Taj, Chaturvedi, Perveen, Bukhari, Siang and Sahgal

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

It is well-known that phosphate-solubilizing bacteria (PSB) promote crop growth and yield. The information regarding characterization of PSB isolated from agroforestry systems and their impact on wheat crops under field conditions is rarely known. In the present study, we aim to develop psychrotroph-based P biofertilizers, and for that, four PSB strains (Pseudomonas sp. L3, Pseudomonas sp. P2, Streptomyces sp. T3, and Streptococcus sp. T4) previously isolated from three different agroforestry zones and already screened for wheat growth under pot trial conditions were evaluated on wheat crop under field conditions. Two field experiments were employed; set 1 includes PSB + recommended dose of fertilizers (RDF) and set 2 includes PSB – RDF. In both field experiments, the response of the PSB-treated wheat crop was significantly higher compared to the uninoculated control. In field set 1, an increase of 22% in grain yield (GY), 16% in biological yield (BY), and 10% in grain per spike (GPS) was observed in consortia (CNS, L3 + P2) treatment, followed by L3 and P2 treatments. Inoculation of PSB mitigates soil P deficiency as it positively influences soil alkaline phosphatase (AP) and soil acid phosphatase (AcP) activity which positively correlated with grain NPK %. The highest grain NPK % was reported in CNS-treated wheat with RDF (N–0.26%, P–0.18%, and K-1.66%) and without RDF (N-0.27, P-0.26, and K-1.46%), respectively. All parameters, including soil enzyme activities, plant agronomic data, and yield data were analyzed by principal component analysis (PCA), resulting in the selection of two PSB strains. The conditions for optimal P solubilization, in L3 (temperature-18.46, pH–5.2, and glucose concentration–0.8%) and P2 (temperature-17°C, pH–5.0, and glucose concentration–0.89%), were obtained through response surface methodology (RSM) modeling. The P solubilizing potential of selected strains at <20°C makes them a suitable candidate for the development of psychrotroph-based P biofertilizers. Low-temperature P solubilization of the PSB strains from agroforestry systems makes them potential biofertilizers for winter crops.

 PSB fertilizer psychrotroph principal component analysis response surface methodology 1. Introduction

Phosphorus (P) is not only an essential nutrient for plant growth but is also involved in diverse biochemical processes like lipid metabolism and the biosynthesis of cell membranes and nucleic acids (Ha and Tran, 2014). Despite its importance, P emerges as a limiting nutrient in global agricultural ecosystems (Lin et al., 2016). In plants, P plays an important role in regulating crucial physiological processes like photosynthesis and nutrient stress mitigation (Xiong et al., 2018). These properties of P as a nutrient make it an indispensable element that neither is substituted nor is replaced in agriculture. As per current research updates, the agriculture P requirement is approximately 0.1–0.5% of the total P concentration (Sharma et al., 2013). The soil P becomes unavailable to plants as it forms complexes with Al3+, Ca2+, and Fe2+ as phosphates (Vance et al., 2003; Siedliska et al., 2021). A large quantity of phosphatic fertilizers (P2O5) is added to restore soil for the available P deficit. This negatively affects the soil quality and puts an economic burden. Therefore, alternative methods for using P-based fertilizers economically without sacrificing soil quality are required to support sustainable agriculture. Since PSBs mobilize P to plants, which would otherwise remain deposited in the soil, the application of phosphate-solubilizing bacteria (PSB) becomes imperative. A century's worth of agricultural P would be supplied by the bacterial-mediated mobilization of free P from complex ore (Sharma et al., 2013). PSB comes under a broad category of potential plant growth-promoting rhizobacteria (PGPR) which provide nutritional value to plants via interaction with the rhizospheric part of the plant making them a possible candidate for biofertilizers. PGPR such as Azospirillum, Bacillus, Enterobacter, Rhizobium, Serratia, Flavobacterium, and Pseudomonas could act as possible biofertilizers. Bacteria can be utilized alone or in combination, for example, the application of consortium (combination of Azotobacter and PSB) can increase wheat crop output (Jain et al., 2021). Wheat is a winter crop and temperatures as low as 15°C are prevalent. The capacity of PSB to bind P in the soil microenvironment is affected by temperatures as low as 15°C. Screening PSBs that can promote P solubilization and plant growth and development even at low temperatures is therefore essential. As a result, utilizing plant-beneficial psychrotolerant bacteria offers an alternative strategy for increasing agriculture output by lessening the effects of low temperatures (Akhtar et al., 2016; Patni et al., 2018; Adhikari et al., 2021).

The most significant advantage of using these psychrotrophs is their ability to colonize cold habitats which are sometimes referred to as cold active microorganisms (CAMs). CAMs can be both psychrophilic and psychrotrophic. Pseudomonas, Staphylococcus, Rahnella, and Stenotrophomonas are a few significant agricultural bacterial groups that are CAMs (Vyas et al., 2010; Dolkar et al., 2018; Kadioglu et al., 2018; Araya et al., 2021). This has sparked research on new cold-tolerant bacterial species that are powerful PGPR and can assist reduce the cold effect on a crop in an agroforestry system to maximize crop yield (Kadioglu et al., 2018). To ameliorate soil P deficiency, inoculating the soil with psychrotroph PSB is a potential technique. Pseudomonas sp. has been used as a biofertilizer in early winter wheat varieties since it is a powerful psychrotroph (Selvakumar et al., 2013). PSB inoculation also minimizes reliance on chemical fertilizers (Alori et al., 2017). The production of organic acids (OAs) is one of the mechanisms involved in phosphorus (P) solubilization. OAs release P by solubilizing inaccessible organic and inorganic phosphorus and the most significant are gluconic acid, citric acid, lactic acid, and malic acid. Pseudomonas sp. is well-known to produce OA and is linked to its P solubilization ability (Kushwaha et al., 2020). As the agricultural output is positively correlated with soil health, therefore, recording variability in soil-health indicators after inoculation of biofertilizer is crucial for the assessment of the overall impact of this strategy. Soil microbial dynamics positively correlate with soil health (Jacoby et al., 2017), and one indicator of soil health is the activity of its enzymes (Dasila et al., 2022). The enzyme pool in the soil is directly influenced by the microbial population, which also directly impacts the availability of micronutrients- and macronutrients (Tahat et al., 2020). Since wheat is one of the most responsive crops toward fertilizers input, therefore, it is the model crop for PSB and plant interaction. The ability of low-temperature P solubilization of these PSB strains makes them a suitable candidate to test them as psychrotroph-based biofertilizers in field conditions. The potential of these PSB strains for plant growth-promoting activity has already been tested in pot conditions (Dasila et al., 2022). Low-temperature P solubilization adaptability of these PSB strains from agroforestry system makes them potential biofertilizers for winter crops. The present investigation also focused on the need of developing biofertilizers from agroforestry system as microorganisms from agroforestry system can be very useful tend in regulating the nutrient cycle due to rich soil organic matter.

 2. Materials and methods 2.1. PSB collection

Four PSB strains, namely Pseudomonas sp.-L3 (Accession No. MG966341), Pseudomonas sp.-P2 (Accession No. MG966347), Streptomyces sp.-T3 (Accession No. MG966352), and Streptococcus sp.-T4 (Accession No. MG966353), were isolated from three Dalbergia sissoo Roxb. provenances in North region of Uttarakhand, India, located at Lachhiwala (30.2099°N latitude, 78.1342°E longitude), Pantnagar (29.0369°N latitude and 79.4472°E longitude), and Tanakpur (29.0722°N latitude, 80.1066°E longitude), India (Joshi et al., 2019).

 2.2. Phosphate and zinc solubilization index

In vitro phosphate and zinc solubilizing ability of four bacterial strains was determined using Pikovskaya's agar (Pikovskaya, 1948) and zinc (Zn) solubilizing agar medium (Sharma et al., 2012). In the Pikovskaya agar, the inorganic P source is tri-calcium phosphate (TCP) whereas in Zn solubilizing medium ZnO serves as an inorganic Zn source. A total of 20 μL of active culture from nutrient broth (NB) was spotted onto Pikovskaya and Zn agar plates and incubated for 5 days at 15°C at a low-temperature incubator. The halo zone appears around the bacterial colony which indicates a positive test for P and Zn solubilizing potential under in vitro conditions. The phosphate-solubilizing index (PSI) and zinc solubilizing index (ZSI) of four PSB strains were calculated (Saravanan et al., 2007).

 2.3. Quantitative estimation of phosphorous and Fourier Transform Infrared Spectroscopy analysis

The quantitative inorganic phosphate-solubilizing potential of four PSB strains was assessed in the National Botanical Research Institute's Phosphate (NBRIP) medium containing glucose (10 g), tri-calcium phosphate-TCP (5 g), MgCl2 (5 g), MgSO4.7H2O (0.25 g), KCl (0.20 g), (NH4)SO4 (0.10 g), distilled water 1 (L), and the pH (7.0 ± 0.2). The log phase culture of four PSB strains was inoculated in NBRIP broth medium and incubated at 30 ± 2°C, 120 rpm in an incubator shaker for 5 days. Thereafter, cell-free supernatant was used to measure soluble P via spectroscopic analysis of reduced phosphomolybdic acid (blue color) by adding α-amino naptho sulfonic acid in the treated sample (Fiske and Subbarow, 1925) and comparing with P standard solution (100 mgL−1). At 5 DPI, prior to estimating P, the pH of bacterial cultures was monitored.

For organic acid detection, FTIR analysis was performed. The log phase culture of selected PSB strains was grown in NB medium and kept at 28 ± 2°C for 72 h in an incubator shaker. Then, the log phase culture of PSB strains was subjected to centrifugation (at 5,000 rpm, for 5 min), followed by filtration via a bacteriological filter (0.2 μM). An equal proportion of filtrate and ethyl acetate were mixed thoroughly for 20 min in a separating funnel, and after vigorous shaking, two layers were developed. The upper (ethyl acetate) layer containing the metabolites was separated, and the extract was concentrated by a rotator evaporator for analysis with FTIR (Sharma et al., 2016).

 2.4. Consortium development

The compatibility of four PSB strains was tested via permutation and combination (Roshani et al., 2020). Freshly grown colony of each of the four PSB strains was inoculated in NB medium and incubated at 28 ± 2°C and 120 rpm for 24 h. The absorbance of bacterial suspension was taken at 600 nm. Finally, an equal culture volume of compatible PSB strains {A600-0.6} was dispensed to 100 mL of NB and thoroughly mixed to develop a PSB-based bacterial consortium (Prasad and Babu, 2017).

 2.5. Seed bacterization and trial establishment

A superior wheat genotype UP 262 was obtained from the seed processing Center (SPC), G.B. Pant University of Agriculture and Technology, Pantnagar, India. The seeds were surface-sterilized with 1% sodium hypochlorite (NaOCl) solution followed by two/three times washing with sterile distilled water. The surface-sterilized seeds were then dried for 1 h in laminar followed by bacterization via immersing them in a log phase (107 cells mL−1) liquid broth culture of four test PSBs and consortia (Tariq et al., 2007; Sharma et al., 2014). Thereafter, three seedlings hill−1 were manually transplanted in the field with a planting geometry of 15–20 cm followed by proper irrigation and monitoring under standard operating procedure (SOP). Field preparation for both experiments was started 1 week before the field setting. The field was plowed, and bunds of 15 cm in height were prepared. Each plot was 2 m long and 2 m wide having a net area of 4 m2. Approximately six lanes per plot were made and leveled carefully after seeds were set.

2.5.1. Experiment field setup

Two sets of the field were prepared. The set 1 (RDF + Treatments) and set 2 (RDF – Treatments) detailed information about the treatments is mentioned in Table 1. RDF used in field experiment 1 was N-80, P-40, and K-40 Kgha−1 with field design given in Supplementary Table 1. Field experiments were set up in randomized block design (Figure 1).

Details of field trial with gross plot size.

S.No. Particular Details
1. Design Randomized block design
2. Crop Wheat (Triticum aestivum.)
3. Variety UP 262
4. Number of replications per treatment 3
6. Total number of plots 18 plots for field set I (treatment + RDF) and 18 plots for field set II (treatment-RDF).
7. Gross plot size 2 × 2 m
8. Net plot size 6 m2
9. Spacing 20 × 20 cm

The detail of treatment in the field lay out showing two experiment sets, Experiment set 1 (Treatment + RDF) and Experiment set 2 (Treatment – RDF).

 2.5.2. Sampling of plants

Plant samples were carefully uprooted at 45 and 90 days post-inoculation (DPI) and analyzed for plant vigor parameters. Soil samples were also taken and analyzed for soil physico-chemical properties and soil enzyme parameters.

 2.5.3. Biometric assessment and yield growth parameters

Wheat samples (45 and 90 DPI) were subjected to biometric observation such as plant height, plant fresh weight, and number of tillers per plant. The grain yield (GY), biological yield (BY), 1,000 g wt, and harvest index (HI) were taken at the final harvesting. At final harvesting, the plants were dried for 72 h followed by thrashing. HI of each plot was calculated with the help of the following formulae:

HI = Economic yieldBiological yieldX 100 2.6. Soil enzyme analysis

For soil health assessment, three soil enzyme activities, e.g., alkaline phosphatase (AP), acid phosphatase (AcP), and fluorescein diacetate (FDA) hydrolysis, were determined to evaluate the impact of PSB inoculation on soil health.

2.6.1. Alkaline phosphatase activity

For AP activity, one g of soil, 250 μl of toluene, 4 mL of modified universal buffer (MUB), and 1 mL of 25 mM p-nitrophenyl phosphate (PNP) were added in a test tube and incubated at 37 ± 1°C for 2 h under gentle shaking condition. Calcium chloride (CaCl2) was added to stop the reaction and Tris buffer was added. The sample was then subjected to centrifugation followed by filtration through Whatman filter paper (0.2 μm) (Agri et al., 2022).

 2.6.2. Acid phosphatase activity

All the steps are the same as described in the determination of AP activity except pH which is adjusted to 5.4 for AcP and 7.4 for AP enzyme activity (Tabatabai and Bremner, 1969).

 2.6.3. Fluorescein diacetate hydrolysis

For FDA estimation, a tri-phenyl tetrazolium chloride (TTC) solution was used. Tris buffer of strength 0.1 M (pH −7.4) and TTC solution of 5 mL were added in a 100 mL flask containing 5 g of soil and incubated for 8 h. To stop the reaction, 25 mL of toluene was added followed by centrifugation at 4,000 rpm for 10 min. The supernatant collected from this was filtered, and the absorbance was read at 485 nm (Inbar et al., 1991).

 2.7. Microbial diversity analysis on different media

The microbial population was measured in different media such as plate count agar (PCA) for total bacteria count, Ashby media for total nitrogen fixers, Aleksandrow media for potassium solubilizers, and Pikovskaya agar for P solubilizers, respectively. The soil samples were serial-diluted followed by pour plating, and then, plates were incubated for 2 days at 28 ± 30° C. Thereafter, bacterial colonies were counted (Messer and Johnson, 2000; Chai et al., 2015).

 2.8. Soil physico-chemical properties

Soil physico-chemical properties (pH, EC, and available NPK) were analyzed before and after the final harvesting of the wheat crop. Soil pH and EC were determined according to the method of Bower and Wilcox (1965) and Jackson (1967), respectively. Organic carbon (OC) of soil was estimated by titration method using the standard protocol of Walkley and Black (1934). For soil, macronutrient estimation of nitrogen (N) was done according to the method of Hanway and Heidel (1952). For phosphorous (P) estimation, Hanway and Heidel (1952) protocol was used, and for potassium (K), Subbiah and Asija (1956) method was used.

 2.9. N, P, and K status of grains

Wheat plants of both the field experiments, set 1 and set 2, were harvested, and their seeds were analyzed for N, P, and K content according to the method of Hanway and Heidel (1952) and Subbiah and Asija (1956), respectively.

 2.10. Electrospray ionization mass spectroscopy analysis for organic acid detection

The organic phase of PSB strains (L3 and P2) was extracted as mentioned earlier in Section 2.3. The organic phase was further subjected to the direct ESI-MS using a Water Q-ToF micro mass spectrometer.

 2.11. Wheat plant chlorophyll and carotenoid estimation

Chlorophyll a, chlorophyll b, and total chlorophyll (TC) of wheat plants were estimated via dimethyl sulfoxide (DMSO) method (Hiscox and Israelstam, 1979) where chlorophyll a, chlorophyll b, and TC were estimated by the following formulae:

       Chlorophyll a (mgg−1) = (12.7 × A663 + 8.02 × A645)1000 × W× V       Chlorophyll b (mgg−1) = (22.9 × A645+ 8.02 × A663)1000 × W× VTotal chlorophyll (mgg−1) = (20.2 × A645+ 8.02 × A663)1000 × W× V

The same procedure was followed for the estimation of carotenoid except that the absorbance of leaf extract was taken at 480 nm (Kirk and Allen, 1965). Carotenoid content was calculated by the following formulae:

Carotenoid (mgg −1) =  [A480 + 0.11 × A663 − 0.638 × A645]1000 × W×V

where W is the weight of leaf sample, and V is the volume of leaf extract.

 2.12. Effect of carbon sources on P solubilizing activity

To analyze the impact of different C sources on the P solubilizing potential of two Pseudomonas strains L3 and P2 which were selected on the basis of plant vigor parameters upon vigorous statistical analysis. To study this, different C sources were used in the NBRIP medium. Log phase culture of L3 and P2 PSB strains were inoculated into NBRIP broth medium amended with 1% of glucose, fructose, maltose, and sucrose separately, uninoculated NBRIP medium serves as a control, and each treatment was done in three replicates. The quantity of solubilized P was measured using the ammonium molybdate method as discussed earlier.

 2.13. Optimization of P solubilization potential on PSB strains using response surface methodology

P solubilization potential of PSB strain L3 and P2 was optimized using the RSM model in Design Expert software (Trial Version 10.0.6), via standardizing the effect of pH, temperature, and C source (glucose) concentration. To run this model, a central composite rotatable design (CCRD) including three variables with four central points was used. RSM model includes a base design that includes both upper and lower limits of the variable along with central points. 3D RSM plots were generated using a mathematical model (Firdous et al., 2017).

 2.14. Principle component analysis

Principle component analysis was done to evaluate the overall response of in vitro and in vivo field results and further mapped them according to the response shared by them and finally resulted in the selection of potential PSB strains.

 3. Results 3.1. PSI and ZSI potential of PSB strains

All four PSB strains tested positive for inorganic P solubilization. The maximum PSI response was observed in P2 (3.6 cm) followed by L3 (3.4 cm), T3 (2.7 cm), and T4 (2.4 cm). The four PSB strains also tested positive for Zn solubilizing potential, and maximum ZSI was observed in T3 (3 cm) followed by P2 (2.8 cm), T4 (2.6 cm), and L3 (1.8 cm) (Figure 2).

PSI (A) and ZSI (B) of L3, P2, T3 and T4 PSB, forming a clear halozone around the bacterial colony.

 3.2. Quantitative estimation of P and FTIR analysis

The amount of soluble P released in all four PSB strains was estimated in the NBRIP broth medium. Maximum soluble P was released in PSB strain L3 (62.83 ± 0.16cd μgmL−1h−1) followed by P2 (Figure 3). In FTIR analysis, a peak of carbonyl bond was detected in all four PSB strains and it corresponded to the carboxylic group (Figure 3). Soluble P released for selected PSB strain (L3 and P2) at temperature 15 ± 2°C was 46.68 and 47.74 μgmL−1h−1, respectively. A significant amount of P has been solubilized even at this low temperature suggesting the cold adaptive solubilization potential of L3 and P2 PSB strain.

(A) Quantitative soluble P release along with pH drop (r = −0.52) in L3, P2, T3, and T4 PSB with correlation value. (B) FTIR analysis of PSB strains showing carboxylic group peak.

 3.3. Consortium development

Out of four PSB strains tested, two strains L3 and P2 were compatible and used for the development of consortia.

 3.4. Plant biometric and vigor parameters

Plant samples uprooted carefully during sampling (45 and 90 DPI) were used for recording agronomic traits.

3.4.1. Field experiment set 1 (PSB + 100% RDF)

In the field experiment set 1, maximum agronomic response in terms of shoot length (SL), shoot fresh weight (SFW), root length (RL), and root fresh weight (RFW) was observed in CNS treatment followed by P2 treatment and agronomic response shared by them is significantly higher than the uninoculated control (C) (Table 2). Overall, there was a positive impact of inoculation of PSB in wheat plants as compared to uninoculated control (C) (Figure 4). The values of plant vigor parameters in the treated plant were significantly higher than in uninoculated control plants. The maximum response was observed in CNS treatment (Figure 6A) with GY (42.07 ± 0.6e), BY (105.84 ± 0.88f), HI (44.30 ± 0.08d), ET (12.16 ± 0.12de), 1,000 g wt (24.07 ± 0.21d), DMA (34.18 ± 0.18d), and GPS (36 ± 0.14e) being followed by P2 treatment (Table 3). In PCA loading plot, plant vigor parameters like BY and GY were found to be positively influenced upon PSB inoculation (Figures 5A, B). During cluster analysis, there were two out-groups, the first being CNS, the most effective, and the second out-group being C treatment, the least effective. Significant variation has been observed in GPS count among the PSB treatment (Figure 5E). There was a clear difference among the treatment in field conditions too (Supplementary Table 3).

Agronomic parameters of PSB treatments with two different sets Set I (PSB + RDF) and Set II (PSB − RDF) at 90 days post-inoculation (DPI), respectively, where SL, shoot length; RL, root length; SFW, shoot fresh weight; RFW, root fresh weight.

Set I (PSB + RDF) Set II (PSB −RDF)
Treatments SL SFW RL RFW SL SFW RL RFW
L3 78.64 ± 1.24de 8.25 ± 0.22d 7.3 ± 0.08cd 0.72 ± 0.01de 60.93 ± 0.41cd 4 ± 0.16d 5.1 ± 0.07c 0.44 ± 0.01b
P2 80.33 ± 1.63e 10.88 ± 0.35e 7.46 ± 0.11d 0.76 ± 0.01de 59.26 ± 0.65cd 3.53 ± 0.05c 5.25 ± 0.04de 0.45 ± 0.03bc
T3 72 ± 1.2c 6.5 ± 0.07c 5.87 ± 0.03b 0.5 ± 0.02c 53.4 ± 0.48bc 3.11 ± 0.08bc 4.73 ± 0.01ab 0.41 ± 0.01a
T4 64 ± 0.57b 5.93 ± 0.05bc 5.42 ± 0.02b 0.47 ± 0.08ab 54.73 ± 1.06bc 2.91 ± 0.07bc 4.75 ± 0.05ab 0.46 ± 0.08bc
CNS 83.33 ± 0.95f 9.38 ± 0.27de 8.4 ± 0.14e 0.82 ± 0.02e 62.66 ± 1.08e 5.42 ± 0.01ef 5.52 ± 0.03e 0.49 ± 0.06c
C 62.46 ± 0.86ab 5.4 ± 0.2ab 5.08 ± 0.2a 0.46 ± 0.06ab 45 ± 0.74a 2.01 ± 0.+ 4.74 ± 0.1ab 0.41 ± 0.01a

Superscript alphabets in the values signify overlap of significance.

Variation among the treatment during sampling of field set 1 (A) and field set 2 (B).

Plant vigor parameters of wheat UP 262 genotype in all PSB treatments in FS 1 (field set 1) and FS 2 (field set 2), where FS 1 = PSB + RDF and FS 2 = PSB – RDF, GY, grain yield; BY, biological yield; HI, harvest index; ET, effective tiller; DMA, dry matter accumulation; GPS, grain per spike.

Treatment GY BY HI ET 1,000 gwt. DMA GPS
FS 1 FS 2 FS 1 FS 2 FS 1 FS 2 FS 1 FS 2 FS 1 FS 2 FS 1 FS 2 FS 1 FS 2
CNTRL 31.00 ± 0.46a 31.77 ± 0.21a 92.85 ± 0.13ab 64.64 ± 0.17a 38.84 ± 0.14a 30.64 ± 0.15ab 8.80 ± 0.16ab 6.06 ± 0.08a 19 ± 0.08a 15.03 ± 0.14ab 18 ± 0.14a 12.24 ± 0.07ab 22.00 ± 0.12a 10 ± 0.01ab
L3 38.84 ± 0.13bc 33.55 ± 0.14bc 98.90 ± 0.07de 82.26 ± 0.22d 41.06 ± 0.58ab 33.06 ± 0.18cd 11.31 ± 0.07d 8.72 ± 0.09bc 23.12 ± 0.4bc 18.60 ± 0.11cd 33.16 ± 0.58cd 13.86 ± 0.06ab 30.00 ± 0.18bc 11 ± 0.02bc
P2 39.05 ± 0.49bc 34.08 ± 0.26bcd 100.31 ± 0.95e 84.86 ± 0.28cd 42.59 ± 0.26bc 32.02 ± 0.34bc 11.01 ± 0.34 8.16 ± 0.07de 22.35 ± 0.21bc 19.30 ± 0.08cd 32.08 ± 0.22bc 14.36 ± 0.04bc 32.00 ± 0.49cde 12 ± 0.01c
T3 37.51 ± 0.09c 32.00 ± 0.19b 93.83 ± 0.14ab 70 ± 0.14b 39.8 ± 0.32bc 31.52 ± 0.26ab 9.77 ± 0.06b 7.14 ± 0.14bcd 22.08 ± 0.06bc 16.08 ± 0.12b 30.14 ± 0.18bc 14 ± 0.02bc 30.00 ± 0.09cde 10 ± 0.02ab
T4 36.72 ± 0.16cd 31.90 ± 0.28a 94.18 ± 0.35cd 72.12 ± 0.31bc 40.36 ± 0.45c 32.18 ± 0.16bc 10.21 ± 0.04cd 7.32 ± 0.16bcd 21.07 ± 0.45ab 16.74 ± 0.14b 30.00 ± 0.16bc 14.24 ± 0.04bc 30.00 ± 0.16cde 12 ± 0.02c
CNS 42.07 ± 0.60e 38.52 ± 0.27d 105.84 ± 0.88f 86.83 ± 0.32e 44.30 ± 0.08de 34.39 ± 0.25cd 12.16 ± 0.12de 9.60 ± 0.12e 24.07 ± 0.21d 19.96 ± 0.16e 34.18 ± 0.18d 16.64 ± 0.05cd 36.00 ± 0.18e 15 ± 0.03cd

Superscript alphabets in the values signify overlap of significance.

PCA (A) and loading plot (B) of plant vigor PSB treatments of type field set 1, PCA plot (C) and loading plot (D) of plant vigor PSB treatments in the field set 2, and hierarchal cluster analysis of PSB treatment in the field set 1 (E) and in the field set 2 (F).

 3.4.2. Field experiment set 2 (PSB – 100% RDF)

The main purpose of PSB inoculation without RDF is to evaluate the potential of PSB strains on wheat plants with persistent inorganic fertilizer which was earlier applied in the soil. In the field experiment set 2 at 90 DPI, maximum agronomic response was observed in CNS treatment followed by L3 treatment. The response shared by them was significantly higher than the uninoculated control (C) (Table 2). Overall, the potential of PSB strains when applied alone was also found to be positive on wheat plants, and their response was quite significant from uninoculated control (C) (Figure 4). The influence on plant vigor parameters was maximum in CNS-treated plants with GY (38.52 ± 0.2d), BY (86.83 ± 0.14e), HI (34.39 ± 0.25cd), ET (9.6 ± 0.12de), 1,000 g wt (19.96 ± 0.16e), DMA (16.64 ± 0.05cd), and GPS (15 ± 0.03cd) being followed by P2 treatment (Table 3). The PCA of the plant vigor data also proved that the highest influence on plant variable was given by CNS followed by L3 and P2 in the order and GY was influenced the most (Figures 5C, D). Cluster analysis of plant vigor data of field experiment set 2 followed a similar pattern (Figure 5F). However, the overall response was low as compared to that in Field Experiment set 1, and GPS count was the most influenced parameter (Figure 6B).

Grain per spike (GPS) variation in different PSB treatments in the field set 1 (PSB + RDF) (A) and field set 2 (PSB – RDF) (B).

 3.5. Soil enzyme results

Alkaline phosphatase and soil acid phosphatase soil enzyme activity represent the organic P conversion to soluble P. In the case of the field experiment set 1, maximum AP activity response was observed in CNS treatment (668.84 ± 0.23f) followed by T4 treatment (666 ± 0.18ef) which is quite significant from uninoculated control (C) soil (466 ± 0.12a μgmL−1h−1). In soil AcP activity, the maximum response was observed in CNS treatment (662.62 ± 0.32e μgmL−1h−1) activity. Soil AcP activity of uninoculated control-C (466 μgmL−1h−1) was found to be significant on the lower side. FDA represents the total soil microbial activity, and among treatments, the maximum response was observed in CNS treatment (48.06 ± 0.014e) followed by P2 (46.64 ± 0.08de μgmL−1h−1) which were significant on the higher side as compared to uninoculated control-C (34.46 ± 0.16a) (Figure 7A). In a field experiment set 2, soil AP activity was maximum in CNS-treated soil (646.64 ± 0.24f μgmL−1h−1) followed by P2 (632.24 ± 0.26de μgmL−1h−1), and control (C) soil sample was found to be 412.68 ± 0.21a μgmL−1h−1, respectively. In soil AcP activity, the maximum response was observed in CNS-treated soil (645.6 ± 0.18f μgmL−1h−1) and T4 (610 ± 0.22de μgmL−1h−1). In FDA activity in soil field experiment set 2, the maximum response was shown by CNS treatment (38.87 ± 0.13cd) and P2 (38.84 ± 0.18cd μgmL−1h−1) (Figure 7B).

Soil enzyme variation includes FDA, AP, and AcP in response to PSB inoculation under field conditions in (A) Set I and (B) Set II after final harvesting.

 3.6. Microbial diversity analysis

The microbial count was taken out in different time intervals, and there was a clear significant difference among the population count, not only in two different field soil types but also within a treatment. Microbial CFU count on different media was taken at 104 cells mL−1. At final harvesting in the field experiment set 1, maximum count was observed in soil treated with CNS (piko-2.55, alexandro-2.26, ashbyii-2.66, kings B-2.76, and PCA-2.86) which were significantly high as compared to uninoculated treated soil (piko-2.46, alexandro-2.2, ashbyii-2.56, kings B-2.48, and PCA-2.7). In a field experiment set 2, a similar pattern was observed in which maximum response was observed in soil treated with CNS (piko-2.12, alexandro-1.78, ashbyii-1.54, kings B-2.22, and PCA-2.66) and in uninoculated control (piko-2, alexandro-1.68, ashbyii-1.62, kings B-2.1, and PCA-2.54) (Supplementary Table 2). Alpha diversity of both field experiment set 1 (Figure 8A) and set 2 (Figure 8B) was measured, and S.H.E analysis was performed where S refers to species richness, H refers to hierarchal information, and E refers to evenness in a community (Figures 8C, D). S.H.E model was applied to check the goodness of fit among the population.

α-diversity indices of the microbial population with different media of field type 1 (A) and S.H.E indices of the field set 1 (B). α-diversity (C) indices and S.H.E indices (D) of the field set 2.

 3.7. Soil physico-chemical analysis

The availability of soil OC, N, P, and K was enhanced in soils treated with PSB in both field experiments. In the field experiment set 1, the highest OC was found in CNS (0.88%)-treated plots followed by L3 (0.84%)-treated plots. However, in field experiment set 2, the highest OC was observed in L3 (0.78%)-treated plots followed by CNS (0.76%)-treated plots. Available N in field experiment set 1 was highest in CNS treatment (280.44 Kgha−1) followed by L3 treatment (274.64 Kgha−1), whereas in field experiment set 2, available N was maximum in P2 treatment (242.66 Kgha−1) followed by CNS (241.65 Kgha−1). Soil available P also varies with the treatments, and in the field experiment set 1, maximum P was found in CNS-treated soil (24.2 Kgha−1) followed by L3-treated soil (22.6 Kgha−1) and in field experiment set 2 too, a similar pattern was observed. Available K content follows the pattern of available P content. The pH and EC of various treatment soils are mentioned in Supplementary Table 3.

 3.8. N, P, and K content of grains

In the field experiment set 1 (PSB + RDF), maximum %N, %P, and %K content was found in CNS-treated wheat plants in which values were 0.35, 0.3, and 1.84, respectively, which were on the higher side of compared to uninoculated control (C) in which %N, %P, and %K were found to be 0.26, 0.18, and 1.66, respectively (Figure 9A). In field experiment set 2 (PSB – RDF), a similar pattern was observed where maximum %N, %P, and %K content was found in CNS with 0.27, 0.26, and 1.46%, respectively, which is significant from the control in which %N, %P, and %K were 0.21, 0.14, and 1.22, respectively (Figure 9B).

Nutrient analysis (N, P, and K) in plant sample (A) field set 1 (PSB + RDF), (B) field set 2 (PSB – RDF).

 3.9. Chlorophyll estimation of wheat plants

Chlorophyll estimation of wheat plants directly links to nitrogen content present in a sample. In the field experiment set 1, maximum Chlorophyll A response was observed in CNS treatment (3.43), whereas in P2, T4, L3, and T3, it was found to be 2.76, 2.54, 2.33, and 2.33 mgg−1, respectively, which were on a higher side as compared to uninoculated control (2.12 mgg−1). For Chl B, the maximum response was observed in T4 (0.56 mgg−1)-treated wheat plants followed by P2 (0.35 mgg−1), T3 (0.28 mgg−1), L3 (0.18 mgg−1), and CNS (0.14 mgg−1) which are quite significant from uninoculated control (0.08 mgg−1). In the case of TC and carotenoid, the maximum response was observed in CNS-treated wheat plants which were found to be 2.31 and 0.6 mgg−1, respectively, which were significantly high as compared to uninoculated wheat plants in which TC and carotenoid content was found to be 1.72 and 0.49 mgg−1, respectively. In the field experiment set 2, chlorophyll A was maximum in CNS-treated (1.96 mgg−1) wheat plants, whereas, in L3, T4, P2, and T3, it was found to be 1.61, 1.45, 1.4, and 1.36 mgg−1, respectively. Most of the treatments, except T3, were significantly higher as compared to the uninoculated control (1.46). For Chlorophyll B, TC, and a carotenoid similar pattern was observed in which CNS-treated plants showed a maximum response and the CNS-treated values were 0.69, 2.18, and 0.59 mgg−1, respectively, which were significantly high when compared to their respective control in which the values of Chlorophyll B, TC, and carotenoid were found to be 0.19, 1.43, and 0.64 mgg−1, respectively. Chord analysis of chlorophyll and carotenoid data in both field experiment set 1 and experiment set 2 shows the maximum % increase share by CNS treatment (Figures 10A, B).

Chord analysis of chlorophyll showing response percent sharing of PSB treatments in the field set 1 (A) and field set 2 (B).

 3.10. ESI-MS analysis of PSB strains

During ESI-MS analysis, two organic acids, e.g., citric acid and maleic acid, were detected in selected L3 and P2 PSB strains, respectively (Figure 11). The production of OA has been considered one of the main mechanisms to solubilize inorganic P and these OAs have been reported earlier in studies.

ESI (Electrospray ionization)-MS results of L3, P2, and CNS for the detection of organic acid.

 3.11. Principal component analysis

Principal component analysis results revealed that CNS treatment is the most effective treatment in both the field condition (set I and set II) followed by L3 followed by P2 PSB treatments. These responses have been interpreted via plotting the performance of PSB treatments against plant vigor parameters in the component form which further results in selection of L3 and P2 PSB treatment (Figure 5).

 3.12. Different C source effects on PSB solubilization

Maximum inorganic P solubilization in two PSB strains L3 (80 μg mL−1h−1) and P2 (78 μg mL−1h−1) was observed when glucose was used as a C source. P solubilization by two strains L3 and P2 in the following sugars is as follows: fructose (L3-74 and P2-72 g mL−1h−1), lactose (L3-62 and P2-64 g mL−1h−1), and sucrose (L3-68 and P2-64 g mL−1h−1), respectively.

 3.13. Optimization of P solubilizing potential of PSB strains

Response surface methodology plot of both PSB strains, e.g., L3 and P2, clearly indicates that quantitative P solubilization correlated positively with pH, temperature, and sugar concentration (Figure 12). However, the increase in P quantification was up to some points above that further increase in variable length decrease in P release. Most optimized condition for maximum P solubilization in L3 comes out to be at temperature-18.46°C, pH-5.2 with sugar concentration of 0.8% and for P2 optimized conditions are at temperature-17°C, pH-5.0 and with sugar concentration 0.89% respectively. However, a significant amount of P solubilization has been observed in both (L3 and P2) PSB strains even at low temperature and pH conditions which suggests that L3 and P2 PSB can be used as efficient P solubilizing psychrotroph.

RSM plots for optimizing P quantification of L3 (A–C) and P2 (D–F) PSB treatment including the interaction of three variables, e.g., temperature, pH, and sugar concentration.

 4. Discussion

As phosphorous (P) is one of the major limiting factors in agriculture crop production, P solubilization via PSB provides a sustainable strategy to release the unavailable or precipitated soil P in an available form to plants (Messer and Johnson, 2000; Rawat et al., 2022). Although there is a high amount of total P content in the soil, its major portion is unavailable for plants. In this situation, the role of PSB alone or in combination (as a consortium) with inorganic or organic fertilizers as potent biofertilizers to release available P from soil to plants becomes prominent (Sharma et al., 2013). The present study established the role of potential PSB isolated from the Dalbergia forest ecosystem not only as CAM but also as PGPRs that enhance wheat growth under field conditions via influencing soil health and soil microbial communities positively when applied with RDF or without RDF. A major concern for any biofertilizer is its performance under field conditions. All four PSB were positive for in vitro P and Zn solubilization. However, maximum solubilization was in two Pseudomonas strains L3 and P2. The range of PSI and ZSI reported was similar to those reported in previous studies. Suleman et al. (2018) reported a PSI in the range of 3.2–5.8 cm for Pseudomonas sp., and Zhang et al. (2020) reported a PSI of 3.6 cm for Pseudomonas lini KM349410. The amount of P solubilized in the NBRIP broth medium ranges from 54.54 to 62.83 μgmL−1 with the strain L3 solubilizing maximum P (62.83 μgmL−1). Fahsi et al. (2021) reported that Pseudomonas lini strain KM349410 recovered from plant rhizosphere (Jujube) solubilized 69 μgmL−1 of inorganic P. The two strains L3 and P2 solubilized P at temperatures <20°C categorizing them as CAMs or psychrotrophic.

Organic acid production and protonation are two major mechanisms for inorganic P solubilization. P solubilization in L3 and P2 strains is associated with a pH drop. During P solubilization, the rate of protonation (H+) in the medium increases, and this could be the main reason for the pH drop. The pH in Pseudomonas strain P2 decreased to 4.2. A previous study reported a drop in pH (4.38) during P solubilization by a Bacillus strain (Upadhayay et al., 2022). There is an inverse relationship between soluble P and pH drop (r = −0.49) (Chen et al., 2012). Acidification via OA production is considered one most important mechanisms of P solubilization (Upadhayay et al., 2022). PSMs produce numerous organic acids that are reported to solubilize inorganic P in the medium including citric (pKa = 2.93 ± 0.28), oxalic (pKa = 1.38 ± 0.54), glycolic (pKa = 3.74 ± 0.11), and gluconic acids (pKa = 3.35 ± 0.35) (Halder et al., 1990; Akintokun et al., 2007; Patel et al., 2008; and Puente et al., 2009). In our study, two organic acids were detected during P solubilization via ESI-MS analysis (Figure 11), and these were maleic acid with an m/z value of 115.08 and citric acid with an m/z value of 191. These two organic acids might be responsible for P solubilization. Both OAs have been previously detected under P solubilizing conditions, e.g., malic acid (Patel et al., 2008) and citric acid (Puente et al., 2009). The production of carboxylic acids during p solubilization via acidification has been detected in FTIR analysis, for example, a peak corresponding to lactic acid was reported by Tang et al. (2019). Similarly, a broad peak of the carbonyl group mostly corresponding to the carboxylic group present in organic acid has been detected in four PSB strains (Figure 3B). A PSB inoculation promotes wheat development under field conditions. In both, the field experiment set 1 (PSB + RDF) and field experiment set 2 (PSB – RDF), agronomic response in PSB treatments was significantly higher as compared to the uninoculated control (C). Among the treatments, CNS treatment showed the maximum response in both field types. In field experiment set 1 (PSB + RDF), maximum SL was found in CNS treatment which is 0.33-fold higher than uninoculated control followed by L3 and P2 with 0.22- and 0.27%-fold increase, respectively. In the case of plant vigor parameters, the response of CNS treatment was maximum with a 0.13% increase in GY, 0.13 in BY, 0.14% in HI, 0.38% in ET, and 0.12% in 1,000 g wt. In addition to that, a 0.7% increase in SFW, a 0.66% increase in RL, and a 0.75 % increase in RFW were observed for control. In the field experiment set 2 (PSB – RDF), also maximum percent increase was observed in CNS treatment with 0.39% in SL, 1.69% in SFW, 0.16% in RL, and 0.21% in RFW. Plant vigor parameters in field experiment set 2 (PSB – RDF) were maximum in CNS treatment with 0.21% increase in GY, 0.06% in BY, 0.12% in HI, 0.58% in ET, and 0.32% in 1,000 g wt followed by L3 and P2 PSB treatments. The main factor that could be responsible for the maximum response of CNS is co-metabolism in which both PSB strains together (L3 and P2) metabolize or solubilize more complex compounds in the soil as compared to when used alone. This resulted in the greater mobilization of free available macro- and micronutrients and was consistent with the observation of Upadhayay et al. (2022), where treatment with consortia of strains having potential for both P and Zn solubilization showed maximum GY, BY, ET, and 1,000 g wt in rice plants. In our study performance of two Pseudomonas strains, e.g., L3 and P2 when used in combination (CNS treatment) and when used alone showed significant response as compared to uninoculated control, and these results were supported by the study where inoculation of Pseudomonas sp. MS16 resulted in increased wheat biomass and yield (Suleman et al., 2018).

The overall impact of PSB inoculation on soil health was evaluated by determining soil enzyme activities (AP, AcP, and FDA) post-inoculation. AP and AcP activities represent the conversion of organic bound P to the most available inorganic form to plants. In the field experiment set 1, maximum AP (66.64 μgmL−1h−1) and AcP (72.42 μgmL−1h−1) activity was found in CNS-treated soil which is significantly higher than in control where AP and AcP activity was 58 and 64.62 μgmL−1h−1, respectively. The observations were for field experiment set 2 (PSB – RDF). An increase in AP and AcP represents phosphatase activity and, thus, indicates the higher mobilization of free P in soil for utilization by plants. These results were also consistent with an earlier study where greater wheat biomass was accompanied by an increase in AP soil enzyme activity (Mäder et al., 2011). A similar pattern was observed for FDA activity too, FDA represents total microbial activity in the soil, and an increase in soil FDA activity could positively influence soil AP and AcP activity which in turn promotes plant growth. This explanation is supported by a previous experiment where a positive correlation between soil FDA and AP activity has been established for crop development (Dasila et al., 2018). Mobilization of P and other soil nutrients via PSB inoculation leads to better plant development.

To study soil microbial population dynamics, different media such as PCA (total bacterial count), King's B (for Pseudomonas), Pikovskaya agar (for positive phosphate-solubilizing bacteria), and Ashbyii (for N fixer) were used. The maximum α-diversity was found in PCA media and α-diversity among the treatment changes with soil status in both field experiment set 1 (PSB + RDF) and field experiment set 2 (PSB – RDF) as analyzed during S.H.E analysis (Figure 8). Microbial population dynamics change in our study has been supported by a previous study in which soil properties alter microbial population (Rfaki et al., 2018).

Soil physico-chemical properties such as OC, available N, P, and K content were maximum in CNS- and L3-treated soil in both field experiments. This observation provides an explanation that inoculation of PSB promotes soil nutrient status via solubilizing and mobilizing soil nutrients. The results of the present study were in tune with a previous study that reported an increase in soil OC content up to 1.79% upon PGPR inoculation (Shahdi Kumleh, 2020). In the present study, available soil and NPK content also increased as reported by Gupta et al. (2021), upon inoculation of Pseudomonas aeruginosa with a variable dose of fertilizers.

Plant nutrient status also changes with different treatments in both field types 1 and 2. In the field experiment set 1, % N, P, and K content in the field was maximum in CNS treatment with 0.35, 0.3, and 1.84%, respectively, which is significantly higher as compared to control in which % N, P, and K were found to be 0.26, 0.18, and 1.66, respectively. A similar pattern was observed in field experiment set 2. The main reason could be due to increased P and N availability in soil which is later utilized by the wheat plant itself for growth upon PSB inoculation. These results were supported by Elhaissoufi et al. (2020), where inoculation of PSB promotes wheat plant P content, and in another study, where PSB inoculation along with RDF in wheat plants promotes N, P, and K status in wheat plants (Sharma et al., 2013).

In both field experiment set 1 (PSB + RDF) and field experiment set 2 (PSB – RDF), there was a substantial increase in chlorophyll a, chlorophyll b, total chlorophyll (TC), and carotenoid content in the PSB-treated wheat plants, the and maximum was observed in CNS treatments (Figure 10). In the field experiment set 1 (PSB + RDF), chlorophyll a, chlorophyll b, total chlorophyll (TC), and carotenoid content in CNS treatment was 3.4, 0.14, 2.3, and 0.6 mgg−1, respectively, which is significant from control in which chlorophyll a, chlorophyll b, total chlorophyll (TC), and carotenoid content was 2.17, 0.08, 1.72, and 0.49, respectively. This may be due to the greater mobilization of soil P in wheat plants which later, in turn, promotes N content in plants (Schlichting et al., 2015; Adhikari et al., 2021), where mobilization of P causes an increase in chlorophyll and carotenoid content in wheat plants which positively correlated with N content in wheat plants. ESI-MS analysis of selected psychrotroph PSB strains L3 and P2 shows the production of two organic acids citric acid and maleic acid. P solubilization, via OA's production, is the principal mechanism adopted by PSB. Citric acid production by psychrotroph L3 PSB strain seems to have a prominent role in P solubilization and citric acid production by psychrotrophic Pseudomonas sp. which has also been previously reported, and in a similar manner, maleic acid is also reported to associate with P solubilization (Zaheer et al., 2019). Different C sources or sugars (glucose, fructose, lactose, and sucrose) were used to find out which C source supports maximum P solubilization in L3 and P2 PSB strain and that carbon source comes out to be glucose (L3-80 and P2-78). This could be because glucose is a primary carbon source and Pseudomonas sp. can use it via various metabolic pathways; these results are consistent with Suleman et al. (2018), in which Pseudomonas sp. MS16 solubilized more P in presence of glucose as a C source. Upon selecting glucose as a C source, the interaction effect of C source concentration, pH, and temperature on P solubilization potential was evaluated using RSM software. For both PSB strains, L3 and P2 conditions were optimized. Optimized conditions for both L3 and P2 PSB strains were found to be maximum at temperatures <20°C suggesting their application as biofertilizers in winter crops. RSM plots of both L3 and P2 PSB strains signify that the P solubilization rate increases with temperature and glucose concentration but decreases with an increase in pH. Upon analysis, it was found that temperature and glucose were positively correlated with P solubilization and pH is negatively correlated and these were supported by Suleman et al. (2018), a study in which similar observations were recorded.

Potential plant growth-promoting rhizobacteria including PSM play a significant role in P solubilization (Sharma et al., 2013), plant growth, and health promotion in various plants including wheat (Kalam et al., 2020; Kusale et al., 2021; Vafa et al., 2021; Ilyas et al., 2022). They have been seen as effective biostimulants (Hamid et al., 2021) and bioinoculants (Basu et al., 2021) under normal and cold stress conditions (Yadav et al., 2019; Baba et al., 2021). EPS-producing bacteria serve as potential biofloculant to mitigate cold stress (Sheikh et al., 2022).

Phosphate-solubilizing bacteria strains L3 and P2 have been submitted to the National Bureau of Agriculturally Important Microorganisms (NBAIM) culture collection with the following accession number Pseudomonas paralactis (L3)-NAIMCC-B-03170 and Pseudomonas aeruginosa (P2)-NAIMCC-B-03171, respectively.

 5. Conclusion

To the best knowledge of our case study, this is the first report of PSB isolated from an agroforestry ecosystem (Dalbergia sissoo Roxb.) and their applications in winter wheat crops for growth and development under field conditions. The present study investigated the P solubilization potential of four PSB strains (L3, P2, T3, and T4) in vivo first followed by their application in wheat crops under field conditions. The response of all the PSB treatments in terms of plant agronomic and plant vigor parameters was higher as compared to the uninoculated control. PSB inoculation also promotes soil health and mitigates soil P into wheat plants via promoting NPK content in wheat seedlings, and later, two PSB strains (L3 and P2) and consortia were selected via exhaustive statistical analysis and used for further studies. Low-temperature growth and P solubilization of potential L3 and P2 make them suitable candidates for developing psychrotroph-based biofertilizers that can be used for winter crops. RSM plots also suggest the optimized P solubilization at a temperature <20°C for both L3 and P2 PSB strains.

 Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

 Author contributions

HD: writing and original draft preparation. MS: conceptualization, methodology, and supervision. VJ, VS, and LT: editing and reviewing of manuscript. GT: inputs for framing of manuscript. KP: reviewing of manuscript. NB and AK: validation of data. TS: reviewing of manuscript. All authors contributed to the article and approved the submitted version.

 Funding

The authors would like to acknowledge the support provided by Researchers Supporting Project Number RSP2023R358, King Saud University, Riyadh, Saudi Arabia. The authors also gratefully acknowledge Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, India, for providing field and research support.

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1135693/full#supplementary-material

 References Adhikari P. Jain R. Sharma A. Pandey A. (2021). Plant growth promotion at low-temperature by phosphate-solubilizing Pseudomonas spp. isolated from high-altitude Himalayan soil. Microb. Ecol. 29, 1–11. 10.1007/s00248-021-01702-133512536 Agri U. Chaudhary P. Sharma A. Kukreti B. (2022). Physiological response of maize plants and its rhizospheric microbiome under the influence of potential bioinoculants and nanochitosan. Plant Soil 474, 451–468. 10.1007/s11104-022-05351-2 Akhtar M. S. Oki Y. Nakashima Y. Adachi T. Nishigaki M. (2016). Phosphorus stress-induced differential growth, and phosphorus acquisition and use efficiency by spring wheat cultivars. Commun. Soil Sci. Plant Anal. 47, 15–27. 10.1080/00103624.2016.1232089 Akintokun A. K. Akande G. A. Akintokun P. O. Popoola T. O. S. Babalola A. O. (2007). Solubilization of insoluble phosphate by organic acid-producing fungi isolated from Nigerian soil. Int. J. Soil Sci. 2, 301–307. 10.3923/ijss.2007.301.307 Alori E. T. Glick B. R. Babalola O. O. (2017). Microbial phosphorus solubilization and its potential for use in sustainable agriculture. Front. Microbiol. 8, 971. 10.3389/fmicb.2017.0097128626450 Araya M. A. Valenzuela T. Inostroza N. G. Maruyama F. Jorquera M. A. Acuña J. J. (2021). Isolation and characterization of cold-tolerant hyper-ACC-degrading bacteria from the rhizosphere, endosphere, and phyllosphere of antarctic vascular plants. Microorganisms 8, 1788. 10.3390/microorganisms811178833202619 Baba Z. A. Hamid B. Sheikh T. A. Alotaibi S. H. El Enshasy H. A. Ansari M. J. . (2021). Psychrotolerant Mesorhizobium sp. isolated from temperate and cold desert regions solubilizes potassium and produces multiple plant growth promoting metabolites. Molecules 26, 5758. 10.3390/molecules2619575834641302 Basu A. Prasad P. Das S. N. Kalam S. Sayyed R. Z. Reddy M. S. . (2021). Plant growth promoting rhizobacteria (PGPR) as green bioinoculants: recent developments, constraints, and prospects. Sustainability 13, 1140. 10.3390/su13031140 Bower C. A. Wilcox L. V. (1965). Soluble salts. Methods of Soil Analysis: Part 2. Chem. Microbiol. Prop. 9, 933–951. 10.2134/agronmonogr9.2.c1132503433 Chai L. Zhang F. She Y. Banat I. Hou D. (2015). Impact of a microbial-enhanced oil recovery field trial on microbial communities in a low-temperature heavy oil reservoir. Nat. Environ. Pollut. Technol. 14, 455–462. Chen Y. Liu R. Sun C. Zhang P. Feng C. Shen Z. (2012). Spatial and temporal variations in nitrogen and phosphorous nutrients in the Yangtze River Estuary. Mar. Pollut. Bull. 64, 2083–2089. 10.1016/j.marpolbul.2012.07.02022910333 Dasila H. Anjul R. Damini D. Anamika R. Manvika S. Salil T. (2018). Interaction between Dalbergia sissoo Roxb. and Pseudomonas koreensis AS15 strain is cultivar specific. Int. J. Curr. Microbiol. Appl. Sci. 7, 297–306. 10.20546/ijcmas.2018.710.031 Dasila H. Sah V. K. Jaggi V. Sahgal M. (2022). Phosphate solubilizing bacteria (PSB) a potential tool to enhance soil health and wheat vigor parameters in pot trial experiment. Pharma Innovat. 11, 1829–1835. Dolkar D. Dolkar P. Angmo S. Chaurasia O. P. Stobdan T. (2018). Stress tolerance and plant growth promotion potential of Enterobacter ludwigii PS1 isolated from Seabuckthorn rhizosphere. Biocatal. Agricul. Biotechnol. 14, 438–443. 10.1016/j.bcab.2018.04.012 Elhaissoufi W. Khourchi S. Ibnyasser A. Ghoulam C. Rchiad Z. Zeroual Y. (2020). Phosphate solubilizing rhizobacteria could have a stronger influence on wheat root traits and aboveground physiology than rhizosphere P solubilization. Front. Plant Sci. 11, 979. 10.3389/fpls.2020.0097932765544 Fahsi N. Mahdi I. Mesfioui A. Biskri L. Allaoui A. (2021). Phosphate solubilizing rhizobacteria isolated from jujube ziziphus lotus plant stimulate wheat germination rate and seedlings growth. PeerJ 9, e11583. 10.7717/peerj.1158334249493 Firdous S. Iqbal S. Anwar S. (2017). Optimization and modeling of glyphosate biodegradation by a novel Comamonas odontotermitis P2 through response surface methodology. Pedosphere 30, 618–627. 10.1016/S1002-0160(17)60381-3 Fiske C. H. Subbarow Y. (1925). The colorimetric determination of phosphorus. J. Biol. Chem. 66, 375–400. 10.1016/S0021-9258(18)84756-1 Gupta R. Noureldeen A. Darwish H. (2021). Rhizosphere mediated growth enhancement using phosphate solubilizing rhizobacteria and their tri-calcium phosphate solubilization activity under pot culture assays in Rice (Oryza sativa.). Saudi J. Biol. Sci. 28, 3692–3700. 10.1016/j.sjbs.2021.05.05234220220 Ha S. Tran L. S. (2014). Understanding plant responses to phosphorus starvation for improvement of plant tolerance to phosphorus deficiency by biotechnological approaches. Crit. Rev. Biotechnol. 34, 16–30. 10.3109/07388551.2013.78354923586682 Halder A. K. Mishra A. K. Bhattacharyya P. Chakrabartty P. K. (1990). Solubilization of rock phosphate by Rhizobium and Bradyrhizobium. J. Gen. Appl. Microbiol. 36, 81–92. 10.2323/jgam.36.81 Hamid S. Lone R. Mohamed H. I. (2021). Production of Antibiotics From PGPR and Their Role in Biocontrol of Plant Diseases. Plant Growth-Promoting Microbes for Sustainable Biotic and Abiotic Stress Management (Jammu and Kashmir: Springer), 441–461. Hanway J. J. Heidel H. (1952). Soil analysis methods as used in Iowa statecollege soil testing laboratory. Iowa State Coll. Agric. Bull. 57, 1–31. Hiscox J. D. Israelstam G. F. (1979). A method for the extraction of chlorophyll from leaf tissue without maceration. J. Bot. 57, 1332–1334. 10.1139/b79-163 Ilyas N. Akhtar N. Yasmin H. Sahreen S. Hasnain Z. Kaushik P. . (2022). Efficacy of citric acid chelate and Bacillus sp. in amelioration of cadmium and chromium toxicity in wheat. Chemosphere 290, 133342. 10.1016/j.chemosphere.2021.13334234922965 Inbar Y. Boehm M. J. Hoitink H. A. (1991). Hydrolysis of fluorescein diacetate in sphagnum peat container media for predicting suppressiveness to damping-off caused by Pythium ultimum. Soil Biol. Biochem. 23, 479–483. 10.1016/0038-0717(91)90013-A Jackson M. L. (1967). Soil Chemical Analysis. New Delhi: Prentice-Hall of India Pvt. Ltd. Jacoby R. Peukert M. Succurro A. Koprivova A. Kopriva S. (2017). The role of soil microorganisms in plant mineral nutrition—current knowledge and future directions. Front. Plant Sci. 8, 1617. 10.3389/fpls.2017.0161728974956 Jain D. Sharma J. Kaur G. Bhojiya A. A. Chauhan S. Sharma V. . (2021). Phenetic and molecular diversity of nitrogen fixating plant growth promoting Azotobacter isolated from semiarid regions of India. Biomed Res. Int. 167, 705–712. 10.1155/2021/6686283 Joshi S. Jaggi V. Tiwari S. Sah V. K. Sahgal M. (2019). Multitrate phosphate solubilizing bacteria fromDalbergia sissoo Roxb. Rhizosphere in natural forests of Indian Central Himalayas. Environ.Ecol. 37, 894–908. Kadioglu G. B. Koseoglu M. S. Ozdal M. U. Sezen A. L. Ozdal O. G. Algur O. F. (2018). Isolation of cold tolerant and ACC deaminase producing plant growth promoting rhizobacteria from high altitudes. Rom. Biotechnol. Lett. 23, 13479–13486. Kalam S. Basu A. Podile A. R. (2020). Functional and molecular characterization of plant growth promoting Bacillus isolates from tomato rhizosphere. Heliyon 6, e04734. 10.1016/j.heliyon.2020.e0473432904284 Kirk J. T. O. Allen R. L. (1965). Dependence of total carotenoids and chlorophyll “a” and “b” of leaf extracts in different solvents. Biochem. Soc. Trans. 603, 591. Kusale S. P. Attar Y. C. Sayyed R. Z. El Enshasy H. Hanapi S. Z. Ilyas N. . (2021). Inoculation of Klebsiella variicola alleviated salt stress and improved growth and nutrients in wheat and maize. Agronomy 11, 927. 10.3390/agronomy11050927 Kushwaha P. Kashyap P. L. Pandiyan K. Bhardwaj A. K. (2020). Zinc-solubilizing microbes for sustainable crop production: current understanding, opportunities, and challenges. Phytobiomes 281–298. 10.1007/978-981-15-3151-4_11 Lin H. Chu L. Wang X. Yao Z. Liu F. Ai Y. . (2016). Boron, nitrogen, and phosphorous ternary doped graphene aerogel with hierarchically porous structures as highly efficient electrocatalysts for oxygen reduction reaction. New J. Chem. 40, 6022–6029. 10.1039/C5NJ03390J Mäder P. Kaiser F. Adholeya A. Singh R. Uppal H. S. Sharma A. K. . (2011). Inoculation of root microorganisms for sustainable wheat-rice and wheat-black gram rotations in India. Soil Biol. Biochem. 43, 609–619. 10.1016/j.soilbio.2010.11.031 Messer J. W. Johnson C. H. (2000). “Total viable counts. Pour plate technique,” in Encyclopedia of Food Microbiology, Vol. 3, eds R. K. Robinson, C. A. Batt, and P. D. Patel (London: Academic Press), 2154–2158. Patel D. K. Archana G. Kumar G. N. (2008). Variation in the nature of organic acid secretion and mineral phosphate solubilization by Citrobacter sp. DHRSS in the presence of different sugars. Curr. Microbiol. 56, 168–174. 10.1007/s00284-007-9053-017965911 Patni B. Panwar A. S. Negi P. Joshi G. K. (2018). Plant growth promoting traits of psychrotolerant bacteria: a boon for agriculture in hilly terrains. Plant Sci. Today. 5, 24–28. 10.14719/pst.2018.5.1.352 Pikovskaya R. I. (1948). Mobilization of phosphorous in soil in connection with the vital activity of some microbial species. Mikrobiologiya 17, 362–370. Prasad A. A. Babu S. (2017). Compatibility of Azospirillum brasilense and Pseudomonas fluorescens in growth promotion of groundnut (Arachis hypogea L.). An. Acad. Bras. Cienc. 89, 1027–1040. 10.1590/0001-376520172016061728489199 Puente M. E. Li C. Y. Bashan Y. (2009). Rock-degrading endophytic bacteria in cacti. Environ. Exp. Bot. 66, 389–401. 10.1016/j.envexpbot.2009.04.010 Rawat P. Sharma A. Shankhdhar D. Shankhdhar S. C. (2022). Improvement of phosphorus uptake, phosphorus use efficiency, and grain yield of upland rice (Oryza sativa L.) in response to phosphate-solubilizing bacteria blended with phosphorus fertilizer. Pedosphere. 32, 752–763. Rfaki A. Zennouhi O. Nassiri L. Ibijbijen J. (2018). Soil properties related to the occurrence of rock phosphate-solubilizing bacteria in the rhizosphere soil of faba bean (Vicia faba L.) in Morocco. Soil Syst. 2, 31. 10.3390/soilsystems2020031 Roshani K. A. Singh A. V. Upadhayay V. U. Prasad B. (2020). Development of potential microbial consortia and their assessment on wheat (Triticum aestivum) seed germination. Environ. Ecol. 38, 6–16. Saravanan V. S. Madhaiyan M. Thangaraju M. (2007). Solubilization of zinc compounds by the diazotrophic, plant growth promoting bacterium Gluconacetobacter diazotrophicus. Chemosphere 66, 1794–1798. 10.1016/j.chemosphere.2006.07.06716956644 Schlichting A. F. Bonfim-Silva E. M. Silva M. D. C. Pietro-Souza W. da Silva T. J. Farias L. D. N. (2015). Efficiency of portable chlorophyll meters in assessing the nutritional status of wheat plants. Rev. Bras. Engenharia Agríc. Ambiental 19, 1148–1151. 10.1590/1807-1929/agriambi.v19n12p1148-1151 Selvakumar G. Joshi P. Suyal P. Mishra P. K. Joshi G. K. Venugopalan R. . (2013). Rock phosphate solubilization by psychrotolerant Pseudomonas spp. and their effect on lentil growth and nutrient uptake under polyhouse conditions. Ann. Microbiol. 63, 1353–1362. 10.1007/s13213-012-0594-5 Shahdi Kumleh A. (2020). Effect of plant growth promoting rhizobacteria (PGPRs) on soil chemical properties in a clover-rice cropping system. J. Soil Water Conserv. 9, 89–106. Sharma A. Shankhdhar D. Sharma A. Shankhdhar S. C. (2014). Growth promotion of the rice genotypes by pgprs isolated from rice rhizosphere. J. Soil Sci. Plant Nutr. 14, 505–517. 10.4067/S0718-9516201400500004027315006 Sharma D. Pramanik A. Agrawal P. K. (2016). Evaluation of bioactive secondary metabolites from endophytic fungus Pestalotiopsis neglecta BAB-5510 isolated from leaves of Cupressus torulosa D. Don. 3 Biotech. 6, 1–14. 10.1007/s13205-016-0518-328330281 Sharma S. B. Sayyed R. Z. Trivedi M. H. Gobi T. A. (2013). Phosphate solubilizing microbes: sustainable approach for managing phosphorus deficiency in agricultural soils. Springerplus 2, 1–14. 10.1186/2193-1801-2-58725674415 Sharma S. K. Sharma M. P. Ramesh A. Joshi O. P. (2012). Characterization of zinc-solubilizing Bacillus isolates and their potential to influence zinc assimilation in soybean seeds. J. Microbiol. Biotechnol. 352–359. 10.4014/jmb.1106.0506322450791 Sheikh T. Hamid B. Baba Z. Iqbal S. Yatoo A. Fatima S. . (2022). Extracellular polymeric substances in psychrophilic cyanobacteria: a potential bioflocculant and carbon sink to mitigate cold stress. Biocatal. Agric. Biotechnol. 42, 102375. 10.1016/j.bcab.2022.102375 Siedliska A. Baranowski P. Pastuszka-Wozniak J. Zubik M. Krzyszczak J. (2021). Identification of plant leaf phosphorus content at different growth stages based on hyperspectral reflectance. BMC Plant Biol. 21, 1–17. 10.1186/s12870-020-02807-433413120 Subbiah B. V. Asija G. L. (1956). A rapid procedure for estimation of available nitrogen in soils. Curr. Sci. 25, 259–260. Suleman M. Yasmin S. Rasul M. Yahya M. Atta B. M. Mirza M. S. (2018). Phosphate solubilizing bacteria with glucose dehydrogenase gene for phosphorus uptake and beneficial effects on wheat. PLoS ONE 13, e0204408. 10.1371/journal.pone.020440830240432 Tabatabai M. A. Bremner J. M. (1969). Use of p-nitrophenyl phosphate for assay of soil phosphatase activity. Soil Biol. Biochem. 1, 301–307. 10.1016/0038-0717(69)90012-1 Tahat M. Alananbeh K. A. Othman Y. I. Leskovar D. (2020). Soil health and sustainable agriculture. Sustainability 12, 4859. 10.3390/su12124859 Tang R. Tang T. Tang G. Liang Y. Wang W. Yang J. . (2019). Pyrimethanil ionic liquids paired with various natural organic acid anions for reducing its adverse impacts on the environment. J. Aagri. food Chem. 67, 11018–11024. 10.1021/acs.jafc.9b0364331512873 Tariq M. Hameed S. Malik K. A. Hafeez F. Y. (2007). Plant root associated bacteria for zinc mobilization in rice. Pak. J. Bot. 39, 245–253. Upadhayay V. K. Singh A. V. Khan A. Singh J. Pareek N. Raghav A. (2022). FE-SEM/EDX based zinc mobilization analysis of Burkholderia cepacia and Pantoea rodasii and their functional annotation in crop productivity, soil quality, and zinc biofortification of paddy. Front. Microbiol. 13, 852192. 10.3389/fmicb.2022.85219235602065 Vafa Z. N. Sohrabi Y. Sayyed R. Z. Suriani N. L. Datta R. (2021). Effects of combinations of Rhizobacteria, mycorrhizae, and seaweeds on growth and yields in wheat cultivars under the influence of supplementary irrigation. Plants 10, 811. 10.3390/plants1004081133924128 Vance C. P. Uhde-Stone C. Allan D. L. (2003). Phosphorus acquisition and use: critical adaptations by plants for securing a nonrenewable resource. New Phytol. 157, 423–447. 10.1046/j.1469-8137.2003.00695.x33873400 Vyas P. Joshi R. Sharma K. C. Rahi P. Gulati A. Gulati A. (2010). Cold-adapted and rhizosphere-competent strain of Rahnella sp. with broad-spectrum plant growth-promotion potential. J. Microbiol. Biotechnol. 20, 1724–1734.21193830 Walkley A. J. Black I. A. (1934). Estimation of soil organic carbon by the chromic acid titration method. Soil Sci. 37, 29–38. 10.1097/00010694-193401000-00003 Xiong Y. Liao B. Proffitt E. Guan W. Sun Y. Wang F. . (2018). Soil carbon storage in mangroves is primarily controlled by soil properties: a study at Dongzhai Bay, China. Sci. Total Environ. 619, 1226–1235. 10.1016/j.scitotenv.2017.11.18729734601 Yadav A. N. Kour D. Sharma S. Sachan S. G. Singh B. Chauhan V. S. . (2019). Psychrotrophic Microbes: Biodiversity, Mechanisms of Adaptation, and Biotechnological Implications in Alleviation of Cold Stress in Plants. Plant Growth Promoting Rhizobacteria for Sustainable Stress Management: Volume 1: Rhizobacteria in Abiotic Stress Management (Himachal Pradesh: Springer), 219–253. Zaheer M. S. Raza M. A. S. Saleem M. F. Khan I. H. Ahmad S. Iqbal R. . (2019). Investigating the effect of Azospirillum brasilense and Rhizobium pisi on agronomic traits of wheat (Triticum aestivum L.). Arch. Agron. Soil Sci. 65, 705–712.10.1080/03650340.2019.1566954 Zhang M. Yang L. Hao R. Bai X. Wang Y. Yu X. (2020). Drought-tolerant plant growth-promoting rhizobacteria isolated from jujube (Ziziphus jujuba) and their potential to enhance drought tolerance. Plant Soil. 452, 423–440. 10.1007/s11104-020-04582-5