Phosphate-Solubilizing Fungi and Alkaline Phosphatase Trigger the P Solubilization During the Co-composting of Sorghum Straw Residues With Burkina Faso Phosphate Rock

Phosphate rocks (PR), the primary source of phosphorus (P), are often co-composted with organic materials to enhance P availability. However, the mechanisms of P solubilization in PR-enriched composts are not well elucidated. This study investigated such mechanisms by monitoring the changes in P fractions during composting and by determining the relationships between the physicochemical and biological parameters. Sorghum straw residues were composted alone (Comp), or with 10% PR (P-Comp), or with 10% PR and 10% rhizosphere soil (P-Comp-Soil), and samples were collected at 45, 60, and 180 days for analysis. The labile-P composed of H2O- and NaHCO3-extractable inorganic P (Pi) and organic P (Po), the moderately labile-P extracted by NaOH (Pi + Po), and the unavailable P formed of the HCl-P and residual-P, increased with the progress of the composting. At 180 days, P-Comp-Soil contained the highest amount of labile-P. There were strong and positive correlations between labile-P and the abundance of total fungi, phosphate-solubilizing fungi (PSF), alkaline phosphatase phoD, phosphonatase phnX, acid phosphatase aphA. Although total fungi were much fewer than total bacteria, the PSF mainly triggered the mineral P solubilization. The alkaline phosphatase phoD was the main enzyme leading the organic P mineralization, while the contribution of phosphonatase phnX, acid phosphatase aphA, and siderophore entA to the organic P solubilization was minor. Besides, the bacterial specific-transporter (pstS) gene increased with the increase of labile-P, allowing for immobilization of little fractions of P in microbial cells. This study highlighted the significant role of PSF and alkaline phosphatase in the P solubilization of PR-enriched composts. Furthermore, it showed the benefit of supplementing the PR-enriched composts with rhizosphere soil, a niche of ecologically important source of beneficial microbes.


INTRODUCTION
The rapid increase in the world's human population expected to reach 9.1 billion by 2050 requires to increase food production by 70% (Krishnaraj and Dahale, 2014). Among the efforts to increase agricultural production, sustainable management of soil fertility will play a significant role. Phosphorus (P), the second most essential nutrient after nitrogen for plant growth, is, however, present at deficient low levels of available fractions in many ecosystems (Barber et al., 1963). This P deficiency has become a threat to soil fertility and crop production throughout the world (Ochwoh et al., 2005), and particularly in the humid tropics and subtropics, and sandy soils of the semi-arid tropics (Sanchez et al., 1997). P availability in soils depends on the P fractions, which influence the primary productivity of agricultural systems (Williams et al., 2013). Organic P (Po) fractions are found in plant remains, in composts, and microbial tissue. In contrast, inorganic P (Pi), which exists in various forms such as H 2 PO 4 can form complexes with iron, aluminum, and calcium in soils and becomes unavailable for plant uptake (Theodorou and Plaxton, 1993;Khan and Joergensen, 2009). On the other hand, Sharma et al. (2013), supported by Pii et al. (2015), reported that organic P and inorganic P are present in enough concentrations. Still, their solubility and thus their availability is deficient. The low P availability in soils is corrected by applying phosphatic fertilizers with high levels of available P. These phosphatic fertilizers use phosphate rock as the primary source of P 2 O 5 . However, the supply of phosphate rock (PR) is limited in the whole world. Moreover, the high price of chemical fertilizers and phosphate rock (Sharma et al., 2013) added to the severe environmental threat due to the exploitation of phosphate rock reserves require to consider alternative strategies to sustain soil fertility and crop production. In this perspective, animal manures, straw, and other solid wastes are used as substrates for PR-enriched composts.
There have been suggestions to increase the nutrient bioavailability of PR-enriched composts by inoculations with phosphate-solubilizing microorganisms (Biswas and Narayanasamy, 2006) or phosphate-mineralizing microorganisms (Sunita, 2014). Their higher nutrient content makes these PR-enriched composts more efficient for plant growth, and they are more environmentally safe (Wei et al., 2015). However, the use of microorganisms as bio-fertilizers remains a hard-to-adopt technology, especially for smallholder farmers in developing countries. It requires access to inocula, appropriate storage, and proper handling ability during inoculation, which is not easily applicable by many farmers. Fortunately, it is well known that rhizosphere soils are reservoirs of numerous beneficial microbes (Barea et al., 2005;Vassilev et al., 2006), including P-solubilizing microorganisms. The addition of rhizosphere soil instead of pure microbial inocula may be an alternative way to supply PR-enriched composts with substantial P-solubilizing microbes and help improve P solubilization. However, this remains a hypothesis and should be further studied to clearly understand the magnitude of P solubilization in such PR-enriched composts with rhizosphere soil and the mechanisms that sustain it. Elucidating the P solubilization mechanisms is an essential step forward in the fabrication of better-quality composts with higher available P concentrations. In soil, microbial mineral P solubilization has long been associated with the production of low molecular weight organic acids (Bolan et al., 1997). Gluconic acid has been identified as the most frequently produced organic acid by phosphate-solubilizing bacteria (PSB) (Goldstein, 1995;Richardson and Simpson, 2011). It chelates divalent cations (e.g., Ca 2+ ) in poorly soluble mineral phosphate forms, such as hydroxyapatite or tricalcium phosphate, thus releasing phosphate forms such as HPO 4 2− and H 2 PO − available for plant uptake (Khan et al., 2010;Castagno et al., 2011). Gluconic acid produced from glucose oxidation requires the glucose dehydrogenase (GDH) enzyme (Deppenmeier et al., 2002) released by the activation of quinoprotein glucose dehydrogenase genes (gcd) (Chen et al., 2015). The activity of GDH, however, requires the redox cofactor pyrroloquinoline quinone (PQQ), whose biosynthesis involves a PQQ operon consisting of about 5-11 genes (Choi et al., 2008). Therefore, we targeted the pqq genes to study the abundance of PSB in this study. Since the mineral P solubilization is also carried out by fungi, we included the characterization of phosphatesolubilizing fungi (PSF). In soil, the solubilization of organic P depends much on the release of several enzymes such as acid and alkaline phosphatases, phosphonatases, phytases, and siderophores (Bergkemper et al., 2016). The determination of these enzymes by qPCR-based molecular techniques will give insights about their relative contribution to P solubilization of PR-enriched composts.
We designed this work to determine the changes in phosphorus fractions during the co-composting of sorghum straw residues with phosphate rock in the presence or absence of rhizosphere soil addition. We aimed at understanding the mechanisms that trigger such changes, focusing on the relationships between the physicochemical and biological properties.

Preparation and Monitoring of Composts
In this study, we used sorghum straw as an organic material substrate for the fabrication of three compost types. One of them was a sole compost (Comp) mainly composed of sorghum straw residues.  (Nakamura et al., 2015). One BPR-enriched compost was without rhizosphere soil (P-Comp) while the other (P-Comp-Soil) received rhizosphere soil collected from farmers' sorghum fields as a microbial source supplement. Just after harvesting the grains, several sorghum stands were uprooted. The loose soil was removed after shaking the root systems, and the soil adhering to the roots, defined as rhizosphere soil, was collected with care and used as inoculum during the composting. In detail, Comp comprised five layers of 20 kg oven-dried sorghum straw each and 460 g urea each. P-Comp was the same as Comp except that we added 10% BPR above each layer. P-Comp-Soil had the same composition as P-Comp, with the only difference that there was an addition of 10% of rhizosphere above each layer. Before making the compost piles, we chopped the sorghum straw to approximately 10 cm segments. An appropriate amount of urea was calculated and added to drop down the initial C/N ratio of sorghum straw from 55/1 to 25/1, which is crucial to initiate organic matter decomposition. There were three replications per compost type, giving a total of nine piles of 1.5 m × 1.5 m × 0.6 m size. The compost piles set on plastic sheets had their moisture levels adjusted to 65% before they were covered with the same plastic sheets. Compost piles were made on 22 March 2018 in INERA (National Institute for Environmental and Agricultural Research) Kamboinse (Burkina Faso) and monitored until 22 September 2018. Every 2 weeks, the piles were uncovered, and temperatures at the 20 cm deep layer were recorded in triplicates, using a Fluorescent portable dissolved oxygen meter DOP-01 (Automatic System Research Co., Ltd., Saitama, Japan). The moisture level was readjusted to 65% after the measurement. When returning the piles, we took necessary care not to contaminate a pile with the content of the others.

Compost Sampling
On days 45, 60, and 180 of composting, we collected a composite compost sample from each pile by combining aliquots of 3 subsamples taken after turning the pile. Immediately after collection, each fresh sample (˜500 g) was immediately put in tight plastic bags and stored at −20 • C until chemical and molecular analysis. Just after the last sampling at 180 days, the piles were uncovered and left to air-dry entirely. The air-dried composts were put in nylon bags and stored at room temperature until further use.

Extraction of DNA From Compost, Soil, and Sorghum Straw Samples
Before DNA extraction, compost samples stored at −20 • C were left at room temperature for 1 h to melt the attached ice. Following that, 200 mg of fresh compost per sample, 250 mg of the rhizosphere soil, and 100 mg sorghum straw were used for DNA extraction using DNeasy Soil Kit (Qiagen). We confirmed the presence of DNA bands by performing a 2% TAE agarose gel electrophoresis and assessed the DNA concentration by Qubit HS DNA Assay (Invitrogen, Carlsbad, CA). The moisture content of compost, soil, and sorghum straw at the time of DNA extraction was calculated after drying in a heat block an aliquot of known water content at 100 • C for 3 h. The% moisture intervenes in subsequent calculations, especially during the processing of qPCR data. DNA was extracted at INERA (Burkina Faso) and transported to Japan in cold conditions together with air-dried composts, sorghum straw, and cold-stored soil used for chemical analysis. The transportation of samples from Burkina Faso to Japan was possible after obtaining transportation permits issued by the indicated authorities of both countries, accompanied by a Material Transfer Agreement.

Chemical Analysis
Nitrate (NO 3 − ), Ammonia (NH 4 + ), pH, Total Carbon, and Total Nitrogen Air-dried compost and sorghum straw materials were ground into fine particles in a Vibrating Sample Mill TI-100 model (CMT Co., Ltd., Tokyo, Japan) for 1 min. Following that, 1 g ground air-dried material of each compost type or sorghum straw of known moisture content (determined after oven-drying at 70 • C for 48 h) was added to 10 ml of 1M KCl in a 50 ml falcon tube, and incubated at 150 rpm, at 25 • C, for 30 min to determine NH 4 -N and NO 3 -N according to Sarr et al. (2019). The soil's NH 4 -N and NO 3 -N were assessed in the same way, except that the soil stored at −20 • C was used. The pH and EC from air-dried compost samples were measured at a 10:1 distilled water-to-compost slurry by the glass electrode method (Thompson et al., 2001). However, soil EC was measured from a saturated paste of the soil (1:1 distilled water:soil). Total nitrogen (N) and total carbon (C) of all materials were determined by the dry combution method using an NC analyzer (Sumigraph NC-220: Sumika Chemical Analysis Service, Osaka, Japan).

P Fractionation of Composts
P in composts was sequentially extracted using a modified procedure of Hedley et al. (1982) and Tiessen and Moir (2008). A 0.5 g dried sample in duplicate was extracted sequentially using 30 ml deionized water (H 2 O), 0.5 M NaHCO 3 (pH 8.5), 0.1 M NaOH, and 1 M HCl. A portion of NaHCO 3 and NaOH extracts was acidified to precipitate the suspended organic matter, and then the pH was adjusted. These supernatants were used to analyze the inorganic P (Pi). Total P content of NaHCO 3 , and NaOH extracts were determined after digestion by an acidified ammonium persulfate in an autoclave (103.4 kPa, 121 • C). The difference between the total P and Pi in extracts was the organic P (Po). Residual P was extracted by HNO 3 -H 2 SO 4 acid using high-performance microwave digestion unit (Microwave Laboratory System; Milestone mls 1,200 mega). As a check, we also determined separately the total in each compost type the ignition of the sample at 550 • C for 4 h and extracted by 0.5 M HCl. The Pi in all extracts and digested solutions was measured colorimetrically according to the molybdenum-blue method of Murphy and Riley (1962), using UV-1800 (UV-VIS spectrophotometer, Hitachi). Since P-Com and P-Comp-Soil received BPR or BPR plus soil during composting, the amount of organic matter in a 1 g material would differ. Therefore, we determined the decomposition rate of organic matter by estimating the crude ash content based on TMECC protocol (Thompson et al., 2001).

Determination of Microbial Abundance and Genes Involved in P Solubilization by qPCR
To perform qPCR, we first prepared plasmid standards. The target genes corresponded to 16S rRNA bacteria, ITS fungi, phosphate-solubilizing bacteria (Pyrroloquinoline Quinones pqqB and pqqE), PSF, alkaline phosphatase (phoD), acid phosphatase (aphA), phosphonatase (phnX), siderophore-mediated enterobactin (entA), phosphate specific transporter (pstS). At first, a soil sample was collected from the experimental field of the Japan International Research Center for Agricultural Sciences (JIRCAS) in Tsukuba, Japan (36.0535 • N and 140.0792 • E). The soil DNA was extracted by the Fast DNA Spin Kit for Soil (Mo-Bio, Carlsbad, CA) and used to amplify the target genes. Amplifications were carried out in triplicates for each gene in 20 µl reaction mixtures composed of 10 µl of 2 × Ampdirect Plus, 0.1 µl Biotaq HS, 10 µM of each primer, 7.7 µl sterile water and 2 µl DNA (˜30 ng) purified using Nucleo Spin gDNA Clean-up (Macherey-Nagel). The PCR conditions consisted of an initial denaturation of 95 • C for 10 min, followed by 35 cycles of 94 • C for 30s, annealing temperature for 1 min, 72 • C for 1 min, and a final elongation of 72 • C for 7 min. The best annealing temperature ( Table 1) for each primer set was determined by first performing a gradient PCR. PCR products were purified using Wizard SV Gel and PCR Clean-Up System (Promega, CA, United States) and quantified by NanoDrop. A 2% gel agarose electrophoresis was performed, and PCR products showing bands with the same size as the size of the target gene were ligated to pGEM-T-Easy Vector system 2 following the indicated instructions and cloned to JM109 competent cells. About 6-10 positive clones per gene were incubated for 16 h in LB Amp + liquid medium, and their plasmids were extracted by the QIAprep Spin Miniprep Kit (Qiagen). The integrity of plasmids was checked on 2% agarose gel, and concentrations were determined by NanoDrop method. Plasmids were sequenced by Fasmac Ltd., sequencing company (Kanagawa, Japan). Raw sequences were processed using MEGA version X  to generate the exact nucleotide sequences of each gene after deleting the included sequences of the vector. The obtained edited sequences were subjected to a homology search via BLAST (Altschul et al., 1997) at DDBJ (the DNA Data Bank of Japan) 1 . For each gene, the plasmids showing sequences closer to that of the expected reference were used as standards for qPCR. The initial copy number of plasmid standards is necessary for downstream calculations. It was calculated as follow: The length of the used vector (pGEM-T Easy system 2) is 3,015 bp. The insert's length was determined after processing the sequences data in MEGA X as indicated above. The Avogadro's number unit is in "number/mole." The Dalton representing the value for dsDNA is in "g/mole of bp." 10E + 9 is the 1 http://www.ddbj.nig.ac.jp

pstS-R TGCCGACGGGCCAITYNWC
Primers used to amplify bacterial 16S rRNA ( † ) for plasmid fabrication were different from those ( ‡ ) for qPCR. For the other genes, primers for normal PCR and qPCR were the same. Ta, annealing temperature: the first Ta for a gene is the Ta during normal PCR and the second is the Ta during qPCR.

Statistical Analyses
All physicochemical and biological properties obtained in the three sampling periods were subjected to a two-way ANOVA analysis to estimate the impact of treatments and composting days on these parameters, using CropStat ver. 7.2 (IRRI, Philippines). When no significant interaction was found between treatments and sampling days, the average value of the treatment at the three sampling dates was considered (n = 9). In such a condition, the average of the three treatments at each sampling period was included. However, in the presence of a significant interaction between treatments and the sampling days, a one-way ANOVA was applied to each sampling day separately to estimate the difference level between treatments. When a significant effect was detected (p < 0.05) in both two-way and one-way ANOVA, Fisher's LSD was used to test the differences between treatments. Furthermore, all physicochemical and biological parameters were subjected to a principal component analysis (PCA) method using Kyplot version 5.0 (Kyence Co., Ltd., Japan). The results gathered in the significant principal components (PC1 and PC2) were discussed. The key parameters relating to P fractions and genes' abundances were extracted and analyzed for Spearman's rank correlations, PAST v.2.17 (Hammer et al., 2001). For the Spearman's rank correlation analysis, we considered all the retrieved parameter data collected at 40, 60, and 180 days.

Changes of P Fractions During Composting
The two-way ANOVA results showed non-significant interactions between treatments and sampling days for the three organic P fractions H 2 O-Po, NaHCO 3 -Po, and NaOH-Po. Therefore, each treatment's combined data over the three sampling periods were considered, and they showed no significant differences among treatments, nor sampling days (Supplementary Table S1). For the remaining P fractions, significant interactions were found between treatments and sampling days, following two-way ANOVA. Thus, we applied a one-way ANOVA to treatments at each sampling day, and the results are reported in Table 2. At each sampling period (45, 60, and 180 days), significant p-values (p < 0.05) were obtained for all reported parameters, where the two BPR-enriched composts (P-Comp, P-Comp-Soil) showed significantly higher concentrations than the compost (Comp). Although not significant, the magnitude of the differences between P-Comp and P-Comp-Soil for these parameters increased by the composting age, to be relatively higher at 180 days than at 60 days. Along the composting period, we observed an increase of each P fraction's concentration from 45 to 180 days ( Table 2), which was more pronounced in BPR-enriched composts. For instance, total P at 180 days in P-Com and P-Comp-Soil was 2-and 2.3-fold, respectively, than their corresponding values at 45 days. These two BPR-enriched composts showed similar amounts of total P at 45 days. Their total P levels increased by 18.5 and 29% in P-Comp and P-Comp-Soil, respectively, at 60 days. At 180 days, the increase further reached 104 and 128% that at 45 days for P-Comp and P-Comp-Soil, respectively. As illustrated in Figure 1A, 86-88% of total P In BPRenriched composts corresponded to the HCl-P fraction regardless of the composting period. The compost treatment (Comp) contained only 17-23% of HCl-P but had higher proportions of the other P fractions. Figure 1B shows the cumulative concentrations of labile-P fractions. While at 60 days, P-Comp and P-Comp-Soil showed similar total labile-P levels, the labile-P at 180 days was 13.05% higher in P-Comp-Soil than in P-Comp. Total labile-P in BPRenriched composts at 180 days was twofold that at 45 days. Moreover, 60-67% of the labile-P in BPR-enriched composts was H 2 O_Pi.

Total and Inorganic Nitrogen, Total Carbon, C/N, pH, and EC During Composting
The results of the two-way ANOVA showed significant interactions between treatments and sampling days. The average data for each treatment over the three sampling periods are thus reported in Table 3. It revealed significant differences among treatments for total nitrogen (p < 0.05), total carbon (p < 0.001), and pH (p < 0.05), and among sampling periods for C/N ratio (p < 0.05). P-Comp and P-Comp-Soil contained significantly higher total nitrogen levels than Comp, and significantly lower pH than Comp. However, P-Comp-Soil had significantly less organic carbon contents, followed by Comp, then P-Comp-Soil. Considering the averages of parameters in the three compost treatments at each sampling date, we found an overall significantly less C/N ratio at 180 days than at 45 days. Still the C/N ratio between 60 days was like that at 45 days and that at 180 days.

Abundance of Microbes and Genes Involved in P Solubilization
The two-way ANOVA revealed non-significant interactions between treatments and sampling days for bacteria (16S rRNA) and PSF. The average data over the three sampling periods (n = 9) are thus reported in Supplementary Table S2. While there was no significant difference among treatments for the abundance of total bacteria, P-Comp-Soil contained significantly higher PSF abundance than Comp. However, the differences in PSF abundance between P-Comp and Comp, and P-Comp and P-Comp-Soil, was not significant, respectively. Table 5 shows that the overall abundance of bacteria and PSF was the same among sampling days.
Other than bacteria and PSF, the two-way ANOVA applied to all other biological parameters showed significant interactions between treatments and sampling days. Therefore, a one-way ANOVA was performed for each sampling period's data, and the results are reported in Table 4. At 45 days, the significant LSD differences were obtained for entA, pstS, phnX, pqqB, and aphA. The siderophore entA gene's abundance was substantially higher in Comp than in BRR-enriched composts. The phosphate specific transporter pstS gene was more abundant in P-Comp than in the two other treatments. The abundance of the phosphonatase phnX gene was significantly higher in Comp than in P-Comp, which contained more phnX abundance than the P-Comp-Soil. The PSB pqqB gene was more abundant in P-Comp, while pqqE abundance was higher in Comp. The acid phosphatase aphA gene was significantly higher in BPR-enriched composts (P-Comp and P-Comp-Soil).
At 60 days, the pattern of the significance levels differed from that at 45 days. The significant probability levels (p < 0.05) were found in total fungi, phoD, and entA genes. The total fungal abundance in P-Comp was the same as in Comp and P-Comp-Soil, but P-Comp-Soil contained a significantly higher abundance of fungi than Comp. The phoD gene's abundance was markedly higher in P-Comp-Soil than in P-Comp and Comp. In contrast, P-Comp and Comp held significantly more copy numbers of the entA gene than P-Comp-Soil. At 180 days, except for entA and pqqB genes, all statistical analyses were significant (p < 0.05). The abundance of total fungi was significantly higher in the two BPR-enriched composts than in Comp. The phoD copy number in P-Comp-Soil was considerably higher than that in P-Comp and Comp. The pstS gene was substantially more abundant in P-Comp-Soil than in P-Comp, which had more pstS than Comp. Moreover, P-Comp-Soil contained significantly higher abundances of phnX, pqqE, and aphA genes than P-Comp and Comp. Figure 2A showed the predominance of bacteria over fungi in all compost treatments regardless of the composting stage. At 45 days, the bacteria themselves accounted for 97-98.2% of the total bacteria and total fungi combined. Total bacteria and total fungi were in the same range at 60 and 45 days. However, the proportion of fungi increased in the three composts at 180 days. While at 45 days, the fungi accounted for only 1.8-3% of total bacteria total fungi combined, they reached 11-17.5% in the three compost treatments at 180 days. Although total bacteria dominated total fungi in composts, PSF significantly outnumbered the PSB genes pqqB and pqqE ( Figure 2B). The Figure 2B revealed that 98.25-99.75% of P-solubilizing microbes in composts, regardless of the sampling day, were the PSF. The PSF's proportion increased toward the compost maturation phase, where the ratio of pqqE and pqqB genes decreased. Figure 3 illustrates the distribution of phoD, phnX, and aphA genes, involved in the mineralization of organic P. It shows that the abundance of the alkaline phosphatase phoD was significantly higher than that of the phosphonatase phnX and acid phosphatase aphA. The phoD gene represented at least 93% of the total of these three genes. The proportion of phnX in composts was significantly higher than that of aphA which showed a very little abundances during the composting period. The Figure 3 shows that the percentage of phoD gene in Comp and P-Comp Two-way ANOVA was applied to test the significance of interactions between sampling periods and treatments CropStat ver. 7.2. Since no significant interaction was found for these variables, the average data of each compost treatment during the three sampling periods are reported (n = 9) in the upper panel. The lower panel contain averages of the three compost treatments per sampling date (n = 9). Values (followed by standard errors) with different letters in a column (above each Sign. Trait) are significantly different at p < 0.001 (***), p < 0.05 (*), ns, not significant. Sign. Trait, significant trait. Comp, compost with sorghum straw as organic material; P-Comp, Burkina phosphate rock (BPR)-enriched compost; P-Comp-Soil, BPR-enriched compost with addition of rhizosphere soil; TN, total nitrogen; TC, total carbon; EC, electric conductivity.
slightly increases over the composting period, to be higher at 180 days. The proportion of phnX followed an opposite trend since it was less at 180 days in these two compost types. In contrast, the percentage of phnX gene increased in P-Comp-Soil at 180 days compared to 45 and 60 days. In sorghum straw, the phnX (41%) and aphA (59%) were the only amplified enzymes, while the soil contained mostly the phoD gene (data not shown).

Principal Component Analysis and Correlation Analysis of Significant Parameters
To determine which factor influences P fractions in composts, we performed a PCA using a total of 29 parameters, comprising the sampling dates, seventeen physicochemical parameters (Supplementary Table S1 and Tables 2, 3), and ten biological parameters (Supplementary Table S2 and Table 4). After the PCA, the principal components PC1 and PC2, which showed significant p levels (p < 0.001), were selected. Supplementary  Table S3 shows the factor-loading matrix of the two chosen PCs. PC1 explained a positive very high factor loading for the 16S bacteria (0.997), and moderate factor loading for NaOH-Pi, HCl-Pi, Residual-P, total P, and pstS. This PC1 can be defined as the "bacterial-mediated changes in total P and recalcitrant P forms" even if the bacterial influence in these P forms was weak (Supplementary Figures S1G,H). Furthermore, Supplementary Figures S1A-F shows that PC1 explained by bacterial changes is weakly correlated to labile-P, or not at all related to ITS fungi, PSF, phoD, and aphA. The PC2 had higher factor loadings for most physicochemical properties, and biological factors such as fungi with a factor loading as high as 0.955, PSF, phosphatases (phoD, aphA), and phosphonatase (phnX). This PC2 can be defined as "compost properties change as influenced by fungi and enzymatic activities." The major difference between the two PCs is that while bacteria lead the changes in PC1, fungi and the analyzed enzymes explain the properties changes in PC2. In the PC2, all P fractions, total P, total nitrogen, NH 4 + , EC, and sampling date, provided significant positive values, like the fungi, PSF, and phosphatases (phoD, aphA), and phosphonatase (phnX). In contrast, the C/N ratio, NO 3 − , and pH values were significantly negative. Although PC1 showed significant loads for NaOH-Pi, HCl-Pi, Residual-P, and total P, the loads of these parameters in the PC2 were much higher. Moreover, the PC2 had a more substantial load for the labile-P than PC1. Since PC2 explains powerfully the P fractions changes, including the labile-P, we plotted the PC2 scores obtained during PCA to some selected significant parameters such as labile-P, total fungi (ITS), PSF, phoD, aphA and phnX gene's abundances (Figure 4) to assess how strong these properties influence individually the P fraction changes. Similarly to the results in Supplementary Table S3, Figure 4A confirms that labile-P positively, significantly correlates with PC2 scores. The same significant positive relationships existed between PC2 scores and total fungi (Figure 4B), PC2 scores and PSF (Figure 4C), PC2 scores and phoD (Figure 4D), PC2 scores and aphA (Figures 4E,F), and PC2 scores and phnX (Figures 4G,H). However, the significance of the relationships between PC2 scores and phoD, aphA, and phnX is higher when only data obtained at 180 days are considered (red dots of Figures 4D,F,H).
The parameters for which we did not find significant differences among compost treatments after the two-way ANOVA, and those with non-significant loads in the PCA, all corresponding to the properties that do not directly influence the P fractions, were excluded from the Spearman's correlation analysis. Spearman's correlation test was performed to assess the details relationships between the significant parameters. The correlation's results of these theoretically selected parameters are reported in Table 5. We found meaningful positive relationships between the sampling date and total nitrogen, labile-P, total P, microbial abundances, phoD, and phnX. The labile-P, which is one of the essential fractions in this study, was positively, strongly correlated to total nitrogen, and moderately negatively related to the pH. Furthermore, it was positively, very strongly interrelated to the concentration of total P. Moreover, Spearman's correlations between labile-P and total fungi, PSF, alkaline phosphatase phoD, and acid phosphatase aphA were strongly positive. However, its positive relationship with total bacteria was moderate. Table 5 shows that the pH, other than with the labile-P, was also negatively correlated with total P (strongly) and fungi (moderately). The total nitrogen was positively, strongly to very strongly related to P, fungi, phoD, and moderately to total bacteria. There were positive and strong relationships between total P and total bacteria, total fungi, and phoD, while total P was moderately and weakly related with phnX, and PSF, respectively. We found that bacteria were positively strongly correlated with phnX, and weekly with fungi. Finally, the correlation was positive and strong between phoD and fungi, and moderate between phoD and PSF. These PSFs were positively strongly related to the acid phosphatase aphA. Such positive relationships were missing between total bacteria and phoD and aphA.

DISCUSSION
In this study, we established three compost types using sorghum straw as organic material, supplemented with Burkina Faso phosphate rock (BPR) and a rhizosphere soil collected from farmer's sorghum fields. Based on qPCR analysis, we compared the change evolution during composting of several microbial genes including 16S rRNA total bacteria, ITS fungi, PSF, phosphate solubilizing bacteria (pqqB and pqqE genes), alkaline phosphatase (phoD gene), acid phosphatase (aphA), phosphonatase (phnX gene), siderophore (enterobactin entA gene), and phosphate specific transporter (pstS gene). Also, we determined the different P fractions in composts, the total and available nitrogen levels, the pH, and the electric conductivity from samples collected at 45, 60, and 180 days incubation. Multivariate analyses were then performed to assess the relationships between compost P availability and functional genes.

Evolution of Compost's Physicochemical Parameters
The solid negative Spearman's correlation between sampling dates and C/N ratio (Table 5) indicates a progress of the compost's decomposition through composting time. At 180 days, the lower C/N ratio in the compost treatments enriched with Burkina Faso rock phosphate (BPR) revealed a better compost decomposition in these treatments, which led to higher N concentrations, inorganic P fractions, residual P, HCl-P, and total P as shown in Table 2. Wei et al. (2015) attributed the increase in total P to the so-called "concentration effect" occurring when C, hydrogen, and N are lost with the exit gas as CO 2 , H 2 O, and NH 3 , respectively, but P is retained in the sample. In this study, the addition of rhizosphere soil to the compost enriched with BPR led to some increases in total P and available P at 180 days. Still the increases were not significant compared to the BPRenriched compost without soil. At this late composting phase, the labile-P in P-Comp-Soil was 13% higher than that in P-Comp, although these two treatments had the same concentration of available P at 60 days. Although the difference in available P was not significant between the two treatments, this result further supports that the release of available P in BPR-enriched composts is sustained better by the addition of rhizosphere soil during composting. The rhizosphere which is the narrow zone in the soil near roots (Kumar and Dubey, 2020) provides various ecosystem services, including the cycling of nutrients (Adle, 2016), and it is an enormous reservoir of microbial community and small soil-inhabiting animals (Malla et al., 2018;Dubey et al., 2019). However, the functionality of this rhizosphere microbiome depends on several factors such as the intrinsic natural abundance of microbes in the environment which also relates to the soil physicochemical properties, the type of plant since the nature of the released secondary molecules triggering microbial activities differ from plant to plant, the period of the season owing that microbes may function more actively during wet seasons than dry seasons (Soni et al., 2017;Vyas et al., 2018). In the present study, the sorghum rhizosphere soil used as inoculum was collected toward the end of the growing phase. It is plausible that at this period, the microbiome was less active than if it was collected during the vegetative stage, justifying why there was an increase of labile-P in BPR-enriched compost with soil, but not significant compared to BPR-enriched compost alone. Therefore, it may be necessary to optimize the soil addition rates to obtain a more substantial increase of available P in BPRenriched composts. The labile inorganic P (Pi) corresponding to the sum of H 2 O-Pi and NaHCO 3 -Pi, and which led the variations of total labile-P in treatments compared to the labile organic P (Po) (Supplementary Table S1 and Table 2) amounted half of the available P in the sole compost and over 80% of the labile-P in BPR-enriched composts. We found an increase in the levels of labile Pi in treatments with the progress of the composting, supporting Sunita (2014), who also reported a rise of Pi during composting. NaOH-extracted P, considered as the moderately soluble phosphorus pool (Ahmed et al., 2019), was not included in the readily available P in our study. In compost, NaOH extracted mostly the organic P fraction, while in BPR-enriched composts, it extracted 4-5 times more Pi than Po. The NaOH-extractable P fractions, including Fe-, Al-bound P, humic P, Ca-, Mg-bound P, and moderately labile-P (Lu et al., 2013), increased with the composting age ( Table 2). The last two P fractions HCl-P and residual P are mainly composed of insoluble and stable forms of phosphorus, such as Ca-, Fe-, Al-bound P (Aulakh et al., 2003;Zicker et al., 2018). Our data showed that this non-labile-P pool followed the same increasing trend as the labile-P and the moderately labile-P (Table 2 and Figure 1A) and was higher in the BPR-enriched composts. The non-labile-P accounted for 30-35% of the total P in compost and 92-93% of that in BPR-enriched composts. Therefore, the addition of BPR during composting also increases the chance of supplying higher amounts of insoluble phosphorus. However, Ahmed et al. (2019) stated that HCl-P and residual P fractions of compost might be involved in long-term phosphorus cycling in soil. The higher NaOH-P and HCl-P levels at the end of the composting could be related to P immobilization by microorganisms, which can cause a slower buildup of residual P (Hedley et al., 1982). The mineral solubilization of the highly present CaO (34.39%) in BPR may release free Ca 2+ , which bounds to free Pi and increase the level of residual P when the composting proceeds.

Relationships Between Compost's P Fractions and Microbiological Properties
The primary goal of this work was to determine the changes in P fractions during composting and to relate such changes with microbial functional genes in composts. Identifying the key factors contributing to increase the availability of P in natural ecosystems is crucial, as P is thought to override all other elements, including C and N (Elser, 2012). Adequate P availability positively influences the net productivity in agricultural systems (Williams et al., 2013), and it lays down the primordia of plant reproductive parts during the early phase of plant development (Satyaprakash et al., 2017). Although many soils can contain extensive total P stocks, it is widely accepted that the bioavailability of soluble orthophosphate, which can be used by most biota, is low (Rodriguez and Fraga, 1999).  For each parameter, compost data at 45, 60, and 180 days, and in sorghum straw and soil are considered. Significant Spearman's (rs) correlations are in bold and preceded by probability levels (***p < 0.001, **p < 0.01, *p < 0.5). The strength of the correlations is considered "very weak" if rs values range between "0.00 and 0.19," "weak" between "0.20 and 0.39," "moderate" between "0.40 and 0.59," "strong" between "0.60 and 0.79," and "very strong" between "0.80 and 1.00." TN, Total nitrogen; PSF, phosphate-solubilizing fungi; 16S, 16S rRNA (total bacteria); ITS, total fungi.
Fortunately, microorganisms have proven their ability to increase P availability by mineralization of organic-P and solubilization of inorganic-P (Richardson and Simpson, 2011). Recently, the study of microbes and genomes encoding enzymes responsible for P-cycling has been facilitated by the advance of high-throughput sequencing methods (Yao et al., 2018), eliminating the bias encountered in previous microbial isolation techniques. The qPCR-based molecular analyses applied in the present study provided substantial information on the evolution of P-cyclingrelated genes during the composting period.
We observed that as early as 45 days incubation, the total bacterial abundance in the three composts was already 4.85to 9.85-fold than that of the used sorghum straw and soil (data not shown). This result indicates that bacteria may have increased tremendously during the thermophilic (55-60 • C in average) phase in the first 2 weeks (microbial abundance not quantified during this period) to reach a pick as observed at 45 days. After that, they maintained a more a less stable number up to the end of the composting. Partanen et al. (2010) reported high concentrations of Actinobacteria and Bacillus spp. during the thermophilic phase of composting, which decompose hemicellulose, cellulose, and lignin, as also reported by Yu et al. (2007). The total fungi, which did not show much change during the first 2 months, increased significantly at 180 days in each compost treatment. This result reveals that fungi express well in composts even at lower temperatures. While the average daily temperature was 32 • C, the temperature in composts dropped to 45 • C at 45 days, and 35 • C at 110 days (data not shown). The higher nutrient content in the BPRenriched composts at 180 days could also partly explain the higher abundance of fungi in these treatments. This relationship between compost's nutrients richness and fungal abundance was further supported by the highly significant correlation (r 2 = 0.91) between PC2 scores, which explained most P fractions changes, including the labile-P, and the abundance of ITS fungi (Figure 4B). This observation, in addition to the positive strong Spearman's correlations between the abundance of fungi and total N, labile-P, and total P ( Table 5), justifies the positive interdependence between nutrient availability and fungal abundance. Furthermore, the study showed the significant contribution of PSF in the increase of labile-P (Table 5), particularly at the late stage of the composting ( Figure 4C). The ability of PSFs to transform insoluble organic and inorganic phosphorus into soluble forms has been reported in many studies (Sharma et al., 2013;Kaur and Reddy, 2017;Wang et al., 2018).
In this study, although we found moderate but positive relationships between the abundance of total bacteria and the total nitrogen and labile-P (Table 5), the PCA (Supplementary Table S3) failed to relate the changes in P fractions to any of the analyzed phosphate solubilizing bacterial (PSB) genes pqqB and pqqE. This absence of correlation may arise from non-adequate PSB gene targets or the types of primers used, which might justify why we barely obtained noticeable differences in the abundance of these genes among the compost treatments. Zheng et al. (2017) previously indicated the difficulty of finding homologous genes among all mineral PSB (mPSB) regulating inorganic phosphate solubilization. However, since pqqA, pqqB, pqqC, pqqD, and pqqE genes are conserved and arranged in that order in the pqq operon (pqqABCDE) (Shen et al., 2012), we selected the pqqB and pqqE genes to study the mPSB. Given the substantial wide distribution of PSB in natural ecosystems including soils, freshwater, seawater, and sediments Liang et al., 2020;Wan et al., 2020), we cannot completely rule out a possible contribution of phosphate solubilizing bacteria in the increase of P fractions in our compost treatments. The moderate but positive correlation of 16S bacteria with the labile-P supports the possibility of bacteria to intervene in the P solubilization in the composts. As stated earlier, a lot of PSBs are responsible for the release of organic acids such as gluconic acid (Richardson and Simpson, 2011), which highly contribute to releasing available P for plant uptake from recalcitrant phosphorus forms (Castagno et al., 2011). Since the production of gluconic acid from glucose requires the oxidation of the glucose dehydrogenase (GDH) enzyme released by the activation of the gcd gene (Chen et al., 2015), and that the pyrroloquinoline quinone (PQQ) synthesized by the activation of the PQQ operon with 5-11 genes (Choi et al., 2008) acts as a cofactor, the obtained little concentrations of pqqB and pqqE genes are likely sufficient to trigger the activation of gcd genes and the production of gluconic acid by PSBs. Unfortunately, our attempt to estimate the abundance of gcd genes in our compost treatments was unsuccessful, probably due to non-adequate selected gcd primers or PCR conditions (data not shown). Moreover, we did not determine the levels and identity of organic acids, which might have indicated the presence or absence of PSB in the composts. After all, it is also possible that mPSB via the organic acids production pathway may be minor in this study. In this respect, Schinner (1992, 1995), reported a solubilization of calcium phosphate by Pseudomonas even in the absence of organic acid production. Parks et al. (1990) also postulated that H + excretion originating from NH 4 + assimilation and or respiratory H 2 CO 3 production could be the alternative mechanism of mineral phosphate solubilization. Furthermore, Kim et al. (1997) suggested that the release of organic acids may play an essential role in phosphate solubilization but is not the sole reason for the increase in Pi concentration. On the other hand, some reports (Kucey, 1983;Gyaneshwar et al., 1998;Sharma et al., 2013) pointed out that PSFs possess higher abilities to release P from recalcitrant inorganic P when compared to bacteria. These reports support our finding on the substantial role PSFs play in the mineralization of P in the compost treatments.
In general, the ability of phosphate solubilizing microbes to enhance the availability of P during P-cycling depends not only on the microbial abundance, but also on the nutritional richness of the substrate, the physiological and growth status of the organism (Alori et al., 2017), and the environmental conditions. Seshachala and Tallapragada (2012) revealed that among other factors, the interactions with other microorganisms in the substrate, and the physicochemical properties such as organic matter and pH, also influence microbial phosphate solubilization. In the present study, the clustering of the moderately negative pH load to the PC2 (Supplementary Table S3), indicated that a lower pH may also enhance the dissolution of P in BPR-enriched composts. Indeed, the two-way ANOVA had shown significantly lower pH values in the two BPR-enriched composts (Table 3). However, they were close to neutrality (6.86 and 6.91 in P-Comp-Soil and P-Comp, respectively). Even though lower pH values (acidic-oriented pH status) seemed to trigger P mineralization in composts, the abundance of the acid phosphatase aphA gene was very low compared to that of the alkaline phosphatase phoD gene. Alori et al. (2017) mentioned that in most soils, the pH where phosphate activities were reported ranges from neutral to acidic. Our investigation, which focused on composts and not on soil, also finds a similar range of pH from neutral to acidic, suitable to initiate acid phosphatase activities, even if the acid phosphatase gene abundance was low in this study. In contrast, the high abundance of the phoD gene in composts, and more markedly in the BPR-enriched compost with rhizosphere soil at 180 days, highlights the significant contribution of the alkaline phosphatase gene in P mineralization. Since, the primers of phoD used to encode the alkaline phosphatase covers a high diversity of bacteria and fungi, with a stronger incidence on bacteria (Bergkemper et al., 2016), our results on phoD abundance and its relationships with the compost's physicochemical properties (Supplementary Table S3 and Table 5), indicate that not only the PSF but also phosphate-solubilizing bacteria may take parts in the mineralization of P, as we already stated above. However, Kutu et al. (2018) also reported that the release of alkaline phosphatase was tightly linked to the activity of fungi.
The phosphonatase gene, phnX in the present study, encodes the phosphonatase enzyme, which cleaves the relatively stable carbon-phosphorus bonds occurring in natural and synthetic organophosphonates via C-P lyases and phosphonoacetaldehyde hydrolases (Hsieh and Wanner, 2010). Our PCA results (Supplementary Table S3) showed that apart from phosphatases, the gene encoding the phosphonatase enzyme also plays a significant role in the P mineralization in composts. The effect of phnX on P fractions may be more active toward compost maturity ( Figure 4H). The siderophore entA gene and the phosphate specific transporter pstS were not clustered to the PC2. Thus, their involvement in the P fractions changes during the composting period was nil (entA) or minor (pstS). Instead, the pstS load in the PC1 was significant same as that of total bacteria. Because phosphate specific transporters are harbored by bacteria (Hudek et al., 2016), the significance of the 16S bacteria load in the PC1 may have pulled up that of pstS. Moreover, the labile-P load on the PC1 that also explains pstS, was significant but weak (0.382), indicating the existence relationship between labile-P and pstS. In this respect, we found that at 180 days, the abundance of pstS was significantly higher in P-Comp-Soil, followed by P-Comp, then Comp. It is known that available P is a useful nutrient for microbes. Therefore, the labile-P in composts can activate the expression of the bacterial phosphate-specific transporter (pstS) gene to import the nutrient in microbial cells. In this regard, Vuppada et al. (2018) reported that bacterial cells must maintain intracellular Pi pools for optimal growth and have developed strategies to sense Pi and control the expression of genes to fit best their environmental circumstances. Pst (phosphate specific transporter) and Pit (phosphate inorganic transporter) are the two forms of Pi import systems in bacteria. The Pit import system is activated when the external phosphate concentration is lower than 20 µM (Hudek et al., 2016), while the Pst import system is induced at this concentration (Rao and Torriani, 1990). Since the labile-Pi levels of our compost treatments ( Table 2) were all above this concentration, the periplasmic substrate-binding protein (pstS) can be well activated for the transport and immobilization of Pi in microbes. The immobilized P in bacterial cells can later by solubilized and made available in soil upon compost application to agricultural fields.

CONCLUSION
This study evaluated the mechanisms of P solubilization of Burkina Faso phosphate rock-enriched composts, using sorghum straw as organic material. It demonstrated that the addition of rhizosphere as a microbial source to composts helps improve the solubility of the used low-grade phosphate rock. Quantitative PCR data revealed that although bacteria are substantially more abundant than fungi during the composting process, the PSF significantly increased during the late composting phase and were the leading players of the mineral P solubilization. Furthermore, the alkaline phosphatase (phoD), which dominated the genes encoding for the enzymes involved in organic P mineralization increased its abundance with the progress of the composting and played a substantial role in the organic P mineralization. The increase of inorganic P in compost solution favored the transport of Pi to microbial cells. However, such immobilized P can be later decomposed and released in soil upon application of composts in fields. Based on the observed increase of labile-P in phosphate rock-enriched compost following soil addition, but which was not very significant compared to the enriched compost without soil, we hypothesize that a soil collected just at crop harvest or collected during crop cultivation and stored at the appropriate temperature and light conditions, may provide a superior solubilization potential due to the existence of a more active microbial community.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.