Skip to main content

Elevated atmospheric CO2 alters the microbial community composition and metabolic potential to mineralize organic phosphorus in the rhizosphere of wheat



Understanding how elevated atmospheric CO2 (eCO2) impacts on phosphorus (P) transformation in plant rhizosphere is critical for maintaining ecological sustainability in response to climate change, especially in agricultural systems where soil P availability is low.


This study used rhizoboxes to physically separate rhizosphere regions (plant root-soil interface) into 1.5-mm segments. Wheat plants were grown in rhizoboxes under eCO2 (800 ppm) and ambient CO2 (400 ppm) in two farming soils, Chromosol and Vertosol, supplemented with phytate (organic P). Photosynthetic carbon flow in the plant-soil continuum was traced with 13CO2 labeling. Amplicon sequencing was performed on the rhizosphere-associated microbial community in the root-growth zone, and 1.5 mm and 3 mm away from the root.


Elevated CO2 accelerated the mineralization of phytate in the rhizosphere zones, which corresponded with increases in plant-derived 13C enrichment and the relative abundances of discreet phylogenetic clades containing Bacteroidetes and Gemmatimonadetes in the bacterial community, and Funneliformis affiliated to arbuscular mycorrhizas in the fungal community. Although the amplicon sequence variants (ASVs) associated the stimulation of phytate mineralization under eCO2 differed between the two soils, these ASVs belonged to the same phyla associated with phytase and phosphatase production. The symbiotic mycorrhizas in the rhizosphere of wheat under eCO2 benefited from increased plant C supply and increased P access from soil. Further supportive evidence was the eCO2-induced increase in the genetic pool expressing the pentose phosphate pathway, which is the central pathway for biosynthesis of RNA/DNA precursors.


The results suggested that an increased belowground carbon flow under eCO2 stimulated bacterial growth, changing community composition in favor of phylotypes capable of degrading aromatic P compounds. It is proposed that energy investments by bacteria into anabolic processes increase under eCO2 to level microbial P-use efficiencies and that synergies with symbiotic mycorrhizas further enhance the competition for and mineralization of organic P.

Video Abstract


Phosphorus (P) is fundamentally important to soil biota as a major building block of life [1, 2]. Although many soils have extensive stocks of total P, orthophosphate, low in many soils, is the only form of P available to plants and microbes [3]. Non-labile P in soils requires mobilization to increase P availability and improve P nutrition of plants and other biota [4, 5]. Since organic P comprises up to 80% of total P in soils [6, 7], the mineralization of organic P by soil microorganisms could have the potential to be a prominent process in P transformation, especially the mineralization of a dominant component of organic P, such as phytate accounting for 50–60% of the organic P.

Climate change has the potential to impact organic P transformation. One of most important climate change factors, elevated atmospheric CO2 concentration (eCO2), would considerably accelerate the mineralization of soil organic P [8, 9]. The acceleration of organic P mineralization is attributed to modified plant-soil-microbe interactions due to changes in the plant carbon (C) flow belowground [10,11,12]. Elevated CO2 is known to enhance microbial growth and activity in the rhizosphere of many plant species, as a consequence of plant C efflux in the forms of root exudates and mucilages. Synchronously, plant P demand increases under eCO2 as well [10]. The competition for P between plants and soil microbes in P-limited soils, and balancing C/P stoichiometry in microorganisms with increasing P-use efficiency for microbial population growth likely intensify the mobilization of soil organic P [10,11,12,13]. In a study where crops were grown under eCO2 (550 ppm) over a 7-year period, Jin et al. [11] reported a significant decrease in organic P concentration in two soils, demonstrating the microbial contribution to organic P mineralization.

However, the microbial mechanisms on the acceleration of soil organic P mineralization under eCO2 remain largely unknown. While plant roots release extracellular phosphatases, much of this enzyme activity is restricted to the root surface [14,15,16]. Thus, active mineralization of soil organic P mainly takes place where soil microorganisms interact with plant-derived C. The phylogenetic structure of soil microbial communities has been found to be altered by eCO2 in a number of natural and agricultural systems [17,18,19], which would be associated with the alternation of metabolic potentials in terms of recalcitrant carbon degradation and nutrient mineralization [11, 20]. Recent studies revealed that Actinomycetales, Rhizobiales, Acidobacteriales, and Solibacterales were the potential drivers of P turnover in a beech (Fagus sylvatica)-dominated forest soil [5], and Advenella and Cellulosimicrobium were identified as phytase-producing bacteria to mineralize phytate in soils [21]. Furthermore, in a fungal community, numerous mycorrhizal fungi are able to hydrolyse organic P [22]. Nevertheless, knowledge about the mineralization of organic P under eCO2 is mostly restricted to what specific microbial phylotypes contribute to the mineralization of major organic P compounds in the rhizosphere of crops grown under eCO2 and whether those phylotypes that enhance P turnover are universal across different soils. Relevant studies help us to better understand the impacts of climate change on the microbial ecological service on P availability to crops grown in P-deficient soils, given that approximately 5.7 billion ha of arable lands around the world are P-deficient and the natural P reserves for fertilizers are diminishing [23, 24].

This study quantified the microbial contribution to mineralization of organic P under eCO2 by exploring the community-wide genetic profiling of soil microorganisms, including bacterial and fungal communities and their functions in relation to mineralization of a major organic P compound, phytate, in the rhizosphere of wheat. We hypothesized that (1) eCO2 would increase mineralization of organic P compounds in the rhizosphere, resulting from changes in relevant microbial phylotypes and metabolic potentials of soil microbial communities; (2) independent of soil type, eCO2 would intensify microbial P use in P-deficient soils to such an extent that microbial phylogeny and functional potentials are significantly altered in the process.


Photosynthetic C allocation and phytate mineralization under elevated CO2 (eCO2)

Elevated CO2 significantly increased total biomass of wheat and plant P contents in both Chromosol and Vertosol (Fig. S1). Furthermore, eCO2 resulted in increased C fixation by wheat plants which was indicated by increased 13C enrichment in the shoot and root by 32% and 50%, respectively, across both soils.

Elevated CO2 further increased C deposition belowground as shown by a 1.6-fold of increase in 13C enrichment in the rhizosphere soil across the Chromosol and Vertosol (Fig. S1). Furthermore, eCO2 significantly increased the amount of mineralized phytate in the rhizosphere of wheat (9% and 45% increase in the Chromosol and Vertosol, respectively) (Fig. 1A). By contrast, eCO2 decreased the concentration of Olsen P by 11% and 13% across the rhizosphere compartments for the Chromosol and Vertosol, respectively (Fig. 1B).

Fig. 1
figure 1

The mineralized phytate (A) and Olsen P (B) in the rhizosphere compartments: 0, 1.5, and 3 mm away from the root growth zone of wheat grown for 10 weeks in Chromosol and Vertosol under elevated CO2 (eCO2, 800 ppm) and ambient CO2 (aCO2, 400 ppm). Error bars are standard error (n=6). The dotted lines represent the corresponding value in the bulk soil. The statistical significance for the main effects of CO2 and rhizosphere (Rhizo) and their interaction are presented

Microbial C and P, and respiration in relation to phytate mineralization

Elevated CO2 altered the microbiological characteristics in the rhizosphere. Microbial biomass C increased from 258 μg g-1 (±12.2 SE) under aCO2 to 301 μg g-1 (±13.1 SE) under eCO2 across the rhizosphere compartments of wheat grown in the Chromosol. In the Vertosol, eCO2 increased microbial biomass C from 152 (±7.9 SE) to 172 μg g-1 (±6.6 SE) (Fig. 2A). However, eCO2 did not increase microbial P in the Chromosol decreased it in the rhizosphere in the Vertosol (Fig. 2B). The microbial respiration rate in the rhizosphere significantly (p < 0.01) increased under eCO2 in both soils (Fig. 2C). With the considerable increases in microbial biomass C, eCO2 dramatically increased the microbial-C-to-P ratio and microbial respiration rate per unit microbial P in both soils (Fig. 2D, E). The mineralized P was correlated negatively with Olsen P but positively with the microbial-C-to-P ratio and microbial respiration rate per unit microbial P in two soils (Fig. 3).

Fig. 2
figure 2

The microbial biomass C (A), microbial P (B), soil respiration rate (C), microbial-C-to-P ratio (D), and respiration rate per unit microbial P (E) in the rhizosphere compartments: 0, 1.5, and 3 mm away from the root growth zone of wheat grown in Chromosol and Vertosol for 10 weeks under elevated CO2 (eCO2, 800 ppm) and ambient CO2 (aCO2, 400 ppm). Error bars are standard error (SE, n=6). The dotted lines represent the corresponding values in the bulk soil. The statistical significance for the main effects of CO2 and rhizosphere (Rhizo) and their interaction are presented

Fig. 3
figure 3

The relationships between mineralized phytate and Olsen P (A), microbial-C-to-P ratio (B), and microbial respiration per unit microbial P (C) in the rhizosphere of wheat plants grown for 10 weeks in Chromosol and Vertosol under elevated CO2 (eCO2, 800 ppm) and ambient CO2 (aCO2, 400 ppm). *, **, and *** indicate significance at p < 0.05, p < 0.01, and p < 0.001, respectively

Bacterial and fungal diversity under eCO2 and their association with phytate mineralization

Microbial diversity was more clearly affected by eCO2 in the Chromosol than in the Vertosol. Elevated CO2 significantly increased bacterial species richness by 33% across the 0, 1.5 and 3 mm rhizosphere compartments in the Chromosol (Fig. 4). A similar trend was observed with the Shannon diversity index which increased from 6.5 to 6.8 on average across the rhizosphere compartments in response to eCO2. The bacterial evenness in the rhizosphere of wheat grown in the Chromosol decreased under eCO2 relative to ambient CO2. In the Vertosol, however, species richness and Shannon and evenness indices in the rhizosphere were not significantly affected by eCO2.

Fig. 4
figure 4

The observed bacterial and fungal species richness, Shannon and Pielou evenness indices across rhizosphere compartments of 0, 1.5, and 3 mm away from the root growth zone. Wheat plants were grown in Chromosol and Vertosol for 10 weeks under elevated CO2 (eCO2, 800 ppm) or ambient CO2 (aCO2, 400 ppm). The diversity indices were calculated from rarefied abundance of filtered amplicon sequence variants (ASVs). Bars show the maximum (top edge) and minimum (lower edge) percentiles, and boxes the 25% and 75% percentiles. The median (50%) percentile is represented by the horizontal line within the box. ns, *, **, and *** indicate significance of two-sample Wilcoxon Mann-Whitney tests at p > 0.05, p < 0.05, p < 0.01 and p < 0.001, respectively

Similar to the bacterial community, the fungal community had greater species richness and Shannon indices in the 0, 1.5, and 3-mm rhizosphere compartments when wheat plants were grown under eCO2 in the Chromosol (Figs. 4, S2). However, in the Vertosol, CO2 treatment did not significantly affect any indicators of α-diversity of fungal community.

The non-metric multidimensional scaling (NMDS) analysis further demonstrated that the bacterial and fungal communities in the rhizosphere zone differed from those in the bulk soil and had significant (p < 0.05) responses to eCO2 (Fig. 5) based on the test of permutational multivariate analysis of variance (PERMANOVA). The dissimilarities of the bacterial and fungal communities in the rhizosphere from 0 to 3 mm under eCO2 were clearly separated from those under aCO2 in the Chromosol, indicating that eCO2 changed microbial community composition (Fig. 5). The effect of eCO2 on microbial community composition in the Vertosol followed a similar trend.

Fig. 5
figure 5

Non-metric multidimensional scaling (NMDS) of soil bacterial (A and B) and fungal (C and D) community composition (Bray-Curtis dissimilarities) in wheat rhizosphere compartments of 0, 1.5, and 3 mm away from the root growth zone and the bulk soil. Wheat plants were grown in Chromosol and Vertosol for 10 weeks under elevated CO2 (eCO2, 800 ppm) and ambient CO2 (aCO2, 400 ppm). Ellipses represent 95% confidence contours of samples grouped by ambient and elevated CO2. PERMANOVA results are presented based on Bray-Curtis dissimilarities for the main effects of rhizosphere and CO2 and their interaction on the soil community composition

In order to present phylogenetic-scale-underlying patterns in the relative abundance data of ASVs in this study, we employed the phylofactor analysis to further confirm that both soils contained groups of Bacteroidetes which responded to eCO2 (Model 1) and were associated with increased phytate mineralization (Model 2) (Fig. 6). Through this analysis, the phylogenetic groups or phylofactors were identified and numbered based on the dominant axes of variation corresponding to contrasts of ASVs separated by an edge [25]. Phylofactor 8, 9 and 15 in the Chromosol, which represented the abundances of 112 ASVs in phylum Bacteroidetes, and phylofactor 1 in the Vertosol which represented the abundances of 270 ASVs mainly in phylum Bacteroidetes (and also included Gemmatimonadetes and Chloroflexi), all significantly responded to eCO2 (p < 0.05) (Fig. 6). The majority of these Bacteroidetes (96% of phylofactors 8, 9, and 15 of the Chromosol, 70% of phylofactor 1 of the Vertosol) belonged to families Sphingobacteriaceae, Chitinophagaceae, and Microscillaceae and in the case of the Vertosol also Hymenobacteraceae. In addition, phylofactor analysis showed that a group of Gemmatimonadaceae consisting of 47 ASVs in the Chromosol (Phylofactor 5) and a group of 105 ASVs in the Vertosol (Phylofactors 1 and 6) that mostly belonged to Gemmatimonadaceae, all significantly responded to eCO2 (p < 0.05). The relative abundances of Acidobacteria declined consecutively from the bulk soil to 3-, 1.5-, and 0-mm rhizosphere compartments while Bacteroidetes increased with proximity to roots (Fig. S3).

Fig. 6
figure 6

Effect of elevated CO2 on bacterial phylogenetic groups (phylofactors) based on modeling the isometric log ratios (ILR) of aggregated abundances as response and CO2 treatments (ambient or elevated) as explanatory variables (Model 1) using the package Phylofactor. Samples of all three rhizosphere compartments in Chromosol (A) and Vertosol (B) were combined (n = 18; Chromosol, n = 18; Vertosol). Edge colors in phylogenetic trees indicate their phylum affiliation (left tree). Highlighted groups (right tree) represent selected phylofactors with abundances that were associated with CO2. Boxplots (right) with ILRs of these phylofactors are presented. Taxa encased in dashed red rectangles highlight taxa which were associated with increased phytate mineralization according to a secondary phylofactor model (Model 2)

Moreover, in order to elucidate the association of microbial abundances with phytate mineralization, the regression-phylofactorization model [25] was adopted for the additive analysis of phylofactors which were regressed with mineralized phytate as an explanatory variable (Model 2). The results showed that the combined abundances of all Bacteroidetes and most Gemmatimonadetes were associated with increased phytate mineralization in the rhizosphere, a feature that was consistent across both soils (Phylofactor 1 in Chromosol and Vertosol) (Fig. 6).

Across fungal genera, eCO2 significantly increased the abundances of Basidiomycota genus Agaricus, and the genera Claroideoglomus and Funneliformis of Glomeromycota in the rhizosphere of wheat in the Chromosol, and unidentified Basidiomycota fungi in the order Auriculariales in the Vertosol (Fig. 7).

Fig. 7
figure 7

Effect of elevated CO2 on the abundance (centered-log ratios) of fungal genera in the rhizosphere. Wheat plants were grown in Chromosol (A) and Vertosol (B) for 10 weeks under elevated CO2 (800 ppm) and ambient CO2 (400 ppm). Genera responding significantly (p < 0.05, Holm corrected) to elevated CO2 with a log-fold change > 0.5 or < −0.5 are presented. Bars represent standard errors (n=9)

Microbial metabolic pathway in response to eCO2

It was estimated from bacterial marker genes that eCO2 significantly increased the genetic pool of bacteria for glycolysis (from glucose 6-phosphate) and for the pentose phosphate pathway in the Chromosol and Vertosol (Fig. S4). The pentose phosphate pathway generates pentoses and ribose 5-phosphate, a precursor for the synthesis of nucleotides, and the glycolysis pathway converts glucose to pyruvate, which produces high-energy adenosine triphosphate (ATP) [26, 27]. In addition, the metabolic potentials for salvage pathways of deoxyribonucleosides and ribonucleosides increased with eCO2 in the Chromosol and Vertosol (Fig. S4). As sparse partial least squares regression (sPLS) is an efficient method to deal with genomic selection data in terms of dimension reduction and variable selection [28], this approach has shown that in the Vertosol, the increases in the genetic potentials of glycolysis and the pentose phosphate pathways (among other metabolic pathways) were predictive of increased phytate mineralization (Fig. 8). Moreover, it was observed that in general, potentials of metabolic pathways that were predictive of increased phytate mineralization were associated with decreased available N and increased microbial C (Fig. 8).

Fig. 8
figure 8

Heatmaps of relative abundances of microbial metabolic pathways (MetaCyc database) with soil variables of mineralized phytate, microbial biomass C, and available N across the rhizosphere compartments of wheat grown in Chromosol and Vertosol for 10 weeks under elevated CO2 (eCO2, 800 ppm) in comparison to ambient CO2 (aCO2, 400 ppm). The color indicates the strength of negative (blue to purple) or positive (orange to red) correlations of explanatory variables (pathway potentials) and response variables (soil variables) for the first two components of partial least squares regression. Dendrograms indicate distances of hierarchical clustering


Increased microbial P demand favors P mineralization

Elevated CO2 raised the microbial P-use efficiency for population growth in the rhizosphere of wheat. The increased microbial P-use efficiency was attributable to the increase in bioavailable C in response to increased root exudation under eCO2. A greater amount of plant-derived C in the rhizosphere under eCO2 corresponded with increased microbial growth (Figs. 2 and 4), which is consistent with previous studies using other plant species [29,30,31]. Thus, the increase in microbial biomass within the rhizosphere resulted in more P demands (Fig. 2). However, the increased P-uptake by wheat plants due to eCO2 in this study (Fig. S1) likely enhanced the competition for P between plants and microbes in the rhizosphere, especially in P-deficient soils. The dramatic decrease in Olsen P concentration in the rhizosphere of wheat grown under eCO2 further indicated that microbial competitiveness for P in the rhizosphere would be enhanced (Fig. 1B). Consequently, the microbial C/P ratio in the rhizosphere considerably increased under eCO2 (Fig. 2), implying that the nutrient status of microorganisms in the rhizosphere turned from the C-limiting to P-limiting, and that the P-use efficiency of their population improved. Moreover, the increased respiration per unit microbial P under eCO2 further suggest that eCO2 improved the microbial P use at the metabolic level and this was most pronounced in the Vertosol.

The increase in microbial P-use efficiency under eCO2 may favor the mineralization of phytate. In this study, the mineralization of phytate was accelerated in the rhizosphere of wheat grown under eCO2 (Fig. 1), which is supported by previous FACE (Free Air CO2 Enrichment) studies reporting increased microbial access to the organic P fraction in a number of farming soils in response to eCO2 [10, 11]. Additionally, our present study showed that the increased phytate mineralization and the ensuing P-limitations to the microbial biomass (Figs. 2 and 3) favored the microbial population with a more efficient P metabolism in the rhizosphere. This was most obvious in the Vertosol where the response of phytate mineralization to the microbial C/P ratio was more significant (Fig. 3B).

Bacterial taxonomic and metabolic attributes to P mineralization and microbial P-use efficiency

The actual mechanism of the eCO2-induced increase in phytate mineralization is about how the microbial community composition and function shifted in response to the increased plant C efflux into the rhizosphere under eCO2 and the amount of phytate in soil. In our study, β-diversity rather than α-diversity (species and abundance distributions on sample-level) of the microbial community corresponded with changes in phytate mineralization under eCO2. Species richness, Shannon’s diversity and Pielou’s evenness of bacteria did not change under eCO2 in the rhizosphere of wheat grown in the Vertosol, which was the more P-limiting soil (Fig. 4).

The β-diversity, on the other hand, was more responsive to P mineralization in the rhizosphere under eCO2. The results showed that eCO2 shifted the bacterial community composition in the rhizosphere of wheat which was significant in the Chromosol and marginally significant in the Vertosol (Fig. 5). Hence, it was suspected that across both soils some bacterial taxonomic groups with common metabolic traits responded to eCO2 and the ensuing P-limitations in the rhizosphere. Phylofactor analysis subsequently showed that in both the Chromosol and Vertosol, abundances of two major bacterial phyla, Bacteroidetes and Gemmatimonadetes were relatively enriched in the rhizosphere under eCO2 and positively associated with the mineralization of phytate (Figs. 6, S3). Taken together, this group of bacteria consisted of 375 ASVs (i.e., putative species) in the Vertosol and 159 ASVs in the Chromosol. Although those ASVs were different between the two soils, they likely performed similar metabolic functions regarding the C and P transformation, which will be discussed further below.

In this study, Bacteroidetes and Gemmatimonadetes were likely involved in the processes of C and P transformation in the rhizosphere under eCO2. The most prevalent Bacteroidetes families, Chitinophagaceae and Microscillaceae, followed by Sphingobacteriaceae and Hymenobacteraceae in this group may play such roles across both soils. For example, Chitinophagaceae have been found to be able to mineralize complex organic compounds such as chitin and cellulose [32, 33] and produce β-glucosidase [34]. Moreover, members of the Gemmatimonadaceae bacteria are adapted to low-nutrient environments and are known for their intracellular accumulation of polyphosphate and are indeed utilized for phosphate immobilization in wastewater treatment systems [33]. The similar response of Gemmatimonadaceae to organic P in the water and soil in this present study implies that Gemmatimonadacea play an important role for efficient phytate mineralization of the microbial population across various ecosystems.

A number of microbial families of the Bacteroidetes and Gemmatimonadetes associated with P mineralization had a similar response to eCO2 in the rhizosphere, indicating that similar phylotypes and genes may be involved in the eCO2 response. There were significant associations of P-mineralization with a number of functional genes coding for methylgallate, protocatchuate, and gallate degradation pathways in both soils (Fig. 8). These pathways are involved in the degradation of aromatic compounds as part of the terrestrial C cycle. As phytate is structured with the similar C skeleton to aromatic compounds, the mineralization of C and P might be a co-metabolic process. Spohn and Kuzyakov [35] also demonstrated that P mineralization in several forest soils was driven by microbial need for C, supporting the view of this point.

Elevated CO2 further altered the microbial metabolisms leading to an increase in microbial P-use efficiency. The metabolic mechanisms by which bacteria maximized their growth under eCO2, involved the central pathways for biosynthesis of RNA/DNA precursors. Phosphorus is a crucial constituent in nucleotides assembling to form nucleic acids. For both soils in this study, it was predicted from the ASV abundances that the genetic pool shifted towards the pentose phosphate pathway and pyrimidine deoxyribonucleoside salvage pathways under eCO2 (Fig. S4). The pentose phosphate pathway produces equivalents for anabolic reactions (Nicotinamide adenine dinucleotide phosphate (NDPH) and protects microbial cells against toxic reactive oxygen species (ROS). Furthermore, the de novo synthesis of ribonucleotides is an energy-consuming process and hence bacteria with an advantageous cell membrane morphology and the ability to express key enzymes to facilitate the pyrimidine deoxyribonucleoside salvage pathway under eCO2 have a distinct advantage. The results therefore indicated that when energy supply was abundant (e.g., increased root exudation under eCO2), such phylotypes became more competitive in the P-limiting environment through an increased ability to synthesize and maintain (phosphate-containing) DNA and RNA building blocks. The eCO2-induced shift of microbial metabolisms towards anabolic processes in the rhizosphere is likely to be an adaptative ability of microbial community to energy-enriched but P-deficient environments, which may result in the increase of microbial P-use efficiency under eCO2.

Fungal community contribution to P mineralization

Greater fungal Shannon index was associated with increases in phytate mineralization in both the Chromosol and the Vertosol (Fig. S5). This implied a special role of fungi in the organic P mineralization under eCO2 in both soils. Under eCO2, the relative higher abundance of Agaricus in the Chromosol and Auriculariales in the Vertosol (Fig. 7), which are affiliated to Ectomycorrhiza, likely accessed diverse P pools and organic matters [36, 37]. Moreover, eCO2-induced enrichment of Claroideoglomus and Funneliformis in the rhizosphere in the Chromosol are affiliated with arbuscular mycorrhizas, which may improve the plant capability to compete for P in soil or stimulated bacteria-mediated organic P mineralization processes with increasing the expression of phosphatase genes in bacteria [38]. However, it is unclear whether the fungal community directly contributed to the phytate mineralization or accelerated the depletion of orthophosphate by plants via mycorrhizal hypha, resulting in the bacterial community change towards mining more organic P. In this study, the increase in the relative abundance of those fungal genera, i.e., root-colonizing arbuscular mycorrhizas, under eCO2 was positively associated with mineralized phytate and negatively with Olsen P, suggesting that the contribution of fungal community to the phytate mineralization might be indirectly via colonizing plant roots to deplete orthophosphate, and enhancing rhizobacteria to mineralize organic P.

Moreover, there were several more fungal genera associated with phytate mineralization which might be involved in the decomposition of organic matters for nutrients. In this study, Ascobolus and Arachnomyces were likely relevant to the phytate mineralization (Fig. 7), of which Ascobolus can access nutrients in dung [39] and Arachnomyces have been reported to be able to decompose lignin and cellulose [40]. Therefore, a number of genera in the fungal community might also directly accelerate phytate mineralization by decomposing organic compounds in the soil in the elevated CO2 environment.


Elevated CO2 accelerated the mineralization of phytate in the rhizosphere of wheat. Under eCO2 environments, the increased plant C efflux into the rhizosphere resulted in the stimulation of microbial activity, leading to stronger competition for P in the microbial community. The eCO2-induced increase in the relative abundances of the bacterial genera being able to degrade aromatic compounds likely contributed to the stimulation of phytate mineralization in the rhizosphere. Moreover, the enriched mycorrhizas might form a symbiotic relationship with plants to further deplete labile P in the rhizosphere, which enhanced the rhizobacterial capability to mineralize organic P (Fig. 9). The enhanced microbial anabolic metabolism under eCO2 may reflect an adaptative ability of the microbial community in the rhizosphere to respond to increasing plant-C flow under P-deficient conditions, resulting in the increases of microbial P-use efficiency and organic P mineralization. This microbial ability to metabolically mineralize organic P appears to be similar across the soils. The mineralization of organic P in soil is likely to be accelerated in the rhizosphere of crops grown under eCO2 environments due to the stimulation of microbial metabolisms on C and P transformation. Future study should clarify the link of bacterial and fungal integrative metabolisms with organic P mineralization in the rhizosphere under eCO2.

Fig. 9
figure 9

Diagram of microbial ability to access phytate in the rhizosphere of wheat grown under elevated CO2. The arrows (↑) and (↓) indicate increase and decrease, respectively. The red and blue arrows represent C flow and P dynamics, respectively. The dot arrow indicates unknown process of P transformation

Materials and methods

Experimental design, soils, and growth conditions

A rhizobox experiment was established in growth cabinets where wheat (Triticum aestivum. L., cultivar Yitpi) was grown in two soils, Chromosol and Vertosol [41], under two CO2 concentrations (ambient CO2 (aCO2) and elevated CO2 (eCO2)) in a randomized block design. Twelve rhizoboxes were set up for each soil and planted with wheat, with six replicates grown under aCO2 and six under eCO2. Additionally, four no-plant controls were set up for each soil. Plant materials and rhizosphere soils were destructively harvested and processed after 10 weeks of growth.

The Chromosol soil was collected from the 0−15 cm topsoil layer of a long-term experimental site at the Agriculture Victoria research center in Hamilton, Victoria. (37° 50′ S, 142° 05′ E) [42]. The pasture plots, from which the soil was taken, had received an average of 15 kg P ha-1 year-1 in the form of triple superphosphate since 1977. Prior to the experiment, the soil had a total C of 43 g kg-1, total P of 200 mg kg-1, and Olsen P of 24.5 mg kg-1. The Vertosol soil was collected from the 0−15 cm topsoil at a long-term experiment run by the International Plant Nutrition Institute (IPNI) at Dahlen Victoria, Australia (36° 38′ S and 142° 08′ E). The long-term experiment was established in 1996. The plots were planted in a canola-wheat-barley-pulse rotation, received 80 kg N ha1 year-1 applied as urea and 18 kg P ha-1 year-1 applied as triple superphosphate. The Vertosol had a total soil C of 12.2 g kg-1, total P 172 mg kg-1, and Olsen P 15.0 mg kg-1. Other basic soil physicochemical properties are shown in Table S1.

Uniform wheat seeds were selected for 1-day germination on moistened filter paper under 25°C. Each rhizobox was sown with six geminated seeds and five plants were kept after emergence. All rhizoboxes were allocated into four plant growth cabinets (Fitotron SGC120, Weiss Gallenkamp Ltd., Leicestershire, UK) with two cabinets set at aCO2 (400 ppm) and two at eCO2 (800 ppm). The eCO2 concentration selected was in the range of predictive levels by the end of this century (IPCC 2021). Temperature was set at 22°C (day, 12 h) and 18°C (night), and light intensity (over the waveband 400–700 nm) was at 550 μmol photons m−2 s−1. Soil water content was maintained at 80±5% of field capacity. Two watering tubes were installed in the bulk soil on each side of rhizospheric compartments in each rhizobox. Via these tubes, the deionized water was added into soil daily to achieve the designed soil moisture uniformly in soil by weighing.

Rhizobox setup

The rhizobox was constructed for collecting rhizosphere soil at precise distances from the root growth zone, which was developed by Youssef and Chino [43]. The dimension of the rhizobox was 120 mm in width, and 200 mm in length and height. Each rhizobox comprised two compartment types, i.e., root growth and rhizosperic compartments, which vertically sat inside of the rhizobox. Two rhizospheric compartments were attached to each side of the root-growth compartment (Fig. S6). The thicknesses of root-growth and rhizospheric compartments were 3 mm and 1.5 mm, respectively. Each compartment was separated by PVC frames covered with nylon cloth (mesh size < 26 μm), which allows the movement of water, nutrients, and microorganisms across the compartments.

Phytate (C6H6O24P6Ca6, Sigma Aldrich, USA) at a rate of 70 mg P kg-1 soil was mixed with soil thoroughly before loaded into compartments. Basal nutrients were simultaneously added as well in the following composition (mg kg-1): urea, 60; K2SO4, 147; CaCl2, 186; MgSO4.7H2O, 122; MnSO4.H2O, 6; ZnSO4.7H2O, 8; CuSO4.5H2O, 6; CoCl2.6H2O, 0.4; FeCl3, 0.6; NaB407.10H2O, 1.6 and Na2MoO4.2H2O, 0.4.

13CO2 labeling

Wheat plants were labeled with 13CO2 for 10 days in air-tight clear plastic chambers (area 500 × 540 mm2, height 680 mm) prior to harvest. 13CO2 concentrations of 500 and 900 ppm were attained for aCO2 and eCO2, respectively, after injection of 9.2 M sulfuric acid (H2SO4) into Na213CO3 inside the chamber (Fig. S6). Once the CO2 concentration dropped by 200 ppm, new injections occurred to maintain the 13CO2 levels to an average of 400 and 800 ppm over a 6-h labeling period per day. The frequency of injection depended on the reduction rate of 12CO2 inside the chambers in a preliminary test using a CO2 meter [44]. A fan was placed inside each chamber to homogenize the air. During labeling, half of the pots in each treatment were sealed in chambers while the other half were not subjected to 13CO2 labeling as unlabeled controls.

Harvest and biochemical measurements

After 10 weeks of plant growth, wheat plants were harvested by cutting shoots at the ground level and processed according to Tang et al. [45]. The rhizosphere compartments were carefully pulled out of the rhizobox and dismantled into five individual compartments: one root-growth compartment at 0 mm, and two rhizospheric compartments at 1.5 mm and the other two at 3 mm away from both sides of the root-growth compartment (Fig. S6). Soil samples were collected from the rhizospheric compartments and considered as rhizosphere soil. Soils in those compartments at the same distance from the root-growth compartment on both sides were combined and homogenized. Approximately 5 g of soil was immediately put into liquid nitrogen for 15 min, and then stored at −80°C for DNA extraction. A sample of soil (40 g) was kept at 4°C for available N, microbial biomass and respiration measurements while approximately 30 g of soil was air-dried for measurements of chemical properties.

Shoots and roots were washed and dried at 70°C for 72 h and dry weights recorded. Plant materials were finely ground using a ball mill (Retsch MM400, Germany), and C and P concentrations were measured using a dry-combustion analyzer (PerkinElmer EA2400, Shelton, CT, USA), and a spectrometer [46] after acid digestion [47].

Air-dried soil samples were used to measure Olsen P concentration according to Olsen et al. [48], and pH using a Wettler Toledo 320 pH meter (soil to water = 1:5). The stable 13C isotope ratios in soil and plant samples were analyzed with an isotope ratio mass spectrometer (Deltaplus, Finnigan MAT GmbH, Bremen, Germany). To measure phytate concentrations, 1 g soil was extracted with 0.1 mol L-1 NaOAc with 50 mmol L-1 ethylenediaminetetraacetic acid (EDTA) for 4 h of shaking before centrifuged at 5000 g for 20 min. Afterwards 1 mL of extract was added to 1 mL of 0.1 mol L-1 sodium acetate together with wheat phytase (20 g L-1) (Sigma Aldrich) and incubated at 55°C overnight. An equal volume of 1 mol L-1 hydrochloric acid was used to cease reactions. Phytate P concentrations were quantified as inorganic P after hydrolysis minus the controls without addition of phytase [49, 50]. The malachite green method was used to measure the concentration of inorganic P in solutions [51].

With fresh soil samples, NH4+ and NO3- were determined by a continuous flow analytical system (SKALAR SAN++, Skalar, the Netherlands) after extraction in 2 M KCl for 1 h [52]. The microbial biomass carbon (MBC) was calculated as the difference in total organic C (TOC) concentrations in extracts between fumigated and non-fumigated soils [53]. The TOC concentration in the non-fumigated control was considered as dissolved organic C (DOC) [54]. Microbial P was measured according to Brookes et al. [55]. Using an Infra-red Gas Analyzer (Servomex 4210 Industrial Gas Analyzer, Cowborough, UK) [56, 57], microbial respiration was measured after 10 g fresh soil was incubated in a sealed half-pint (237 ml) Mason jar at 25°C for 24 h.

Marker gene sequencing and bioinformatic processing

Soil DNA was extracted using a DNeasy Powersoil Pro extraction kit (Qiagen, Hilden, Germany) and stored at −20°C. DNA extracts were dissolved in 100 μl of TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8.0). The V4 hypervariable region of bacterial 16S rRNA genes were amplified with the primer pair 515F (GTGYCAGCMGCCGCGGTAA)/806R (GGACTACNVGGGTWTCTAAT) [58]. The Internal Transcribed Spacer (ITS) region 2 of fungi was amplified with primers FITS7 (GTGARTCATCGAATCTTTG)/ITS4 (TCCTCCGCTTATTGATATGC) [59]. Library preparation was performed for 16S and ITS amplicon sequencing using Nextera XT indices according to Illumina reference guides and then sequenced on an Illumina MiSeq platform (2 × 300) [60, 61].

The default pipeline settings of Quantitative Insights Into Microbial Ecology (QIIME2) (version 2020.2; including the plugins cutadapt [62] and dada2 [63] were deployed to assess read quality, trim primers, denoise, filter, and dereplicate sequences. The sequences were truncated to 290 bp. A total of 3,017, 587 and 942,797 reads were retained at a median frequency of 53,362 and 16,955 per sample for bacteria and fungi, respectively. Primer specific taxonomy classifiers were then used (Silva132 at 99% similarity and UNITE v8 dynamic) to assign taxonomies to amplicon sequence variants (ASVs) with the qiime2 plugin “classify- sklearn.” Amplicon sequence variants with less than 10 reads were filtered out resulting in a total of 5 565 bacterial and 1 562 fungal ASVs for diversity analyses.

Enzyme metagenomes and abundances of metabolic pathways (MetaCyc database) [64] were predicted from bacterial ASVs using PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States 2) [65] with default settings. ASVs that were not present in at least two samples and with a frequency of less than 10 reads were removed. The filtered abundance table comprised of 423 pathways.

Diversity and statistical analyses

Analysis of variance (ANOVA) with a block design model [66] in Genstat 19 (VSN International, Hemel Hemspstead, UK) was used to estimate the effects of eCO2 and rhizosphere, and their interactions on the concentrations of phytate, available N (NO3- and NH4+), DOC, and Olsen P, MBC, microbial P, microbial C/P, microbial respiration, and microbial respiration per unit microbial P.

The α-diversity indices were calculated after rarefying abundances to a depth of minimum sample size (27,038 and 6262 reads for bacteria and fungi, respectively). Common diversity indices such as species richness (N0) and Shannon index (H') were calculated with function “estimate_richness” from the Phyloseq package [67] while Pielou evenness (J') was calculated as following [68]:

$${J}^{\prime }=\frac{H^{\prime }}{\ln (N0)}$$

Differences in β-diversity of bacterial and fungal community compositions were visualized on two axes with nonmetric dimensional scaling (NMDS) with Bray-Curtis dissimilarities using the “ordinate” function of the Phyloseq package [67]. Permutational multivariate analysis of variance (PERMANOVA) was then performed using the “adonis” function of the vegan package [69] on the same dissimilarities to assess significant differences of community composition between the rhizosphere compartments and CO2 treatments.

Moreover, phylofactorization was performed to study associations of eCO2 and mineralized phytate with taxonomic groups of bacteria. Bacterial ASVs were regressed into phylogenetic factors based on significant F-statistics at the adjusted p < 0.05 level using ratios of relative abundances as response variable and either CO2 treatments (Model 1) or mineralized phytate (Model 2) as explanatory variables. First, a phylogenetic tree was created using the default options of the Qiime2 plugin q2-fragment-insertion with the silva reference tree (sepp-refs-silva-128, Relative abundances which were transformed to isometric log ratios (ILR) were then regressed along edges of the phylogenetic tree with the package phylofactor [25] as described in the package tutorial ( after pruning the data set to rhizosphere samples only (n = 18 per soil per soil type) and filtering amplicon sequence variants (ASVs) to keep those with a minimum of 50 reads (to minimize undue influence of minor abundances on phylofactors). Briefly, the package phylofactor was applied to break apart the phylogeny with a variety of contrasts and objective functions, summarize the splits, and visualize the circle tree with ASVs generated from rhizosphere soils. Selected phylogenetic groups which were associated with CO2 or mineralized phytate were visualized on a phylogenetic tree and the corresponding ILRs of their abundances (i.e., the ratio from contrast basis elements of abundances of two taxonomic groups) were plotted.

Phylogenetic methods for fungal ITS-based amplicons were not used in this study. Instead, analysis of Compositions of Microbiomes with Bias Correction (ANCOM-BC) [70] was deployed to evaluate the effect of CO2 treatment on the abundance of fungal genera using the package ANCOMBC.

Lastly, the mixOmics package was utilized for sparse partial least squares (sPLS) regression ( to predict mineralized phytate, microbial C, and available N (Y matrix) from Picrust2-predicted metabolic pathway abundances in the rhizosphere as explanatory variables (X matrix). Heatmaps of correlation clusters between important pathways and mineralized phytate, microbial C, and available N of the first two regressed components were subsequently created to visualize their associations. Prior to analysis, the pathway potentials were pruned to those pathways with a minimum relative abundance of 0.011% resulting in a total of 331 pathways. Furthermore, Wilcoxon Mann-Whitney tests were performed on relative abundances (centered log ratio transformed) of each pathway (n = 18/soil) to identify the effect of eCO2 on the bacterial metabolic potentials in the rhizosphere. Accordingly, metabolic pathways that increased in both soils as a result of eCO2 were selected for further review.

Availability of data and materials

Please contact author for data requests.







Dissolved organic C


Total organic C

aCO2 :

Ambient CO2

eCO2 :

Elevated atmospheric CO2


Free air CO2 enrichment


Amplicon sequence variants


Non-metric multidimensional scaling


Permutational multivariate analysis of variance


Analysis of variance


Isometric log ratios


Sparse partial least squares regression


Adenosine triphosphate


Nicotinamide adenine dinucleotide phosphate


Reactive oxygen species


Internal transcribed spacer


Compositions of microbiomes with bias correction


  1. Vitousek PM, Porder S, Houlton BZ, Chadwick OA. Terrestrial phosphorus limitation: mechanisms, implications, and nitrogen–phosphorus interactions. Ecol Appl. 2010;20:5–15.

    Article  PubMed  Google Scholar 

  2. Gong H, Li Y, Yu T, Zhang S, Gao J, Zhang S, et al. Soil and climate effects on leaf nitrogen and phosphorus stoichiometry along elevational gradients. Glob Ecol Conserv. 2020;23:e01138.

    Article  Google Scholar 

  3. Rodriguez H, Fraga R. Phosphate solubilizing bacteria and their role in plant growth promotion. Biotechnol Adv. 1999;17:319–39.

    Article  CAS  PubMed  Google Scholar 

  4. Richardson AE, Simpson RJ. Soil microorganisms mediating phosphorus availability update on microbial phosphorus. Plant Physiol. 2011;156:989–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Bergkemper F, Schöler A, Engel M, Lang F, Krüger J, Schloter M, et al. Phosphorus depletion in forest soils shapes bacterial communities towards phosphorus recycling systems. Environ Microbiol. 2016;18:1988–2000.

    Article  CAS  PubMed  Google Scholar 

  6. Huang L, Jia X, Zhang G, Shao M. Soil organic phosphorus transformation during ecosystem development: a review. Plant Soil. 2017;417:17–42.

    Article  CAS  Google Scholar 

  7. McLaren T, Smernik R, McLaughlin MJ, McBeath TM, McCaskill M, Robertson FA, et al. Soil phosphorus pools with addition of fertilizer phosphorus in a long-term grazing experiment. Nutr Cycl Agroecosys. 2020;116:151–64.

    Article  CAS  Google Scholar 

  8. Fleischer K, Rammig A, De Kauwe MG, Walker AP, Domingues TF. FuchsluegerL, et al. Amazon forest response to CO2 fertilization dependent on plant phosphorus acquisition. Nat Geosci. 2019;12:736–41.

    Article  CAS  Google Scholar 

  9. Jiang M, Caldararu S, Zhang H, Fleicher K, Crous KY, Yang J, et al. Low phosphorus supply constrains plant responses to elevated CO2: a meta-analysis. Glob Chang Biol. 2020;26:5856–73.

    Article  PubMed  Google Scholar 

  10. Jin J, Armstrong R, Tang C. Long-term impact of elevated CO2 on phosphorus fractions varies in three contrasting cropping soils. Plant Soil. 2017;419:257–67.

    Article  CAS  Google Scholar 

  11. Jin J, Wood J, Franks A, Armstrong R, Tang C. Long-term CO2 enrichment alters the diversity and function of the microbial community in soils with high organic carbon. Soil Biol Biochem. 2020;144:107780.

    Article  CAS  Google Scholar 

  12. Samal SK, Dwivedi SK, Rao KK, Choubey AK, Prakash V, Kumar S, et al. Five years’ exposure of elevated atmospheric CO2 and temperature enriched recalcitrant carbon in soil of subtropical humid climate. Soil Till Res. 2020;203:104707.

    Article  Google Scholar 

  13. Osanai Y, Janes JK, Newton PCD, Hovenden MJ. Warming and elevated CO2 combine to increase microbial mineralization of soil organic matter. Soil Biol Biochem. 2015;85:110–8.

    Article  CAS  Google Scholar 

  14. Dodd JC, Burton CC, Burns RG, Jeffries P. Phosphatase activity associated with the roots and rhizosphere of plants infected with vesicular-arbuscular mycorrhizal fungi. New Phytol. 1987;107:163–72.

    Article  CAS  Google Scholar 

  15. Dinkelaker B, Marschner H. In vivo demonstration of acid phosphatase activity in the rhizosphere of soil-grown plants. Plant Soil. 1992;144:199–205.

    Article  CAS  Google Scholar 

  16. Giles CD, Dupuy L, Boitt G, Condron LM BLK, Darch T, et al. Root development impacts on the distribution of phosphatase activity: improvements in quantification using soil zymography. Soil Biol Biochem. 2018;116:158–66.

    Article  CAS  Google Scholar 

  17. Deng Y, Jiang YH, Yang Y, He Z, Luo F, Zhou J. Molecular ecological network analyses. BMC Bioinformatics. 2012;13:113.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Zhou JZ, Deng Y, Luo F, He ZL, Yang YF. Phylogenetic molecular ecological network of soil microbial communities in response to elevated CO2. MBio. 2011;2:e00122–11.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Yu Z, Li Y, Hu X, Jin J, Wang G, Tang C, et al. Elevated CO2 increases the abundance but simplifies networks of soybean rhizosphere fungal community in Mollisol soils. Agr Ecosyst Environ. 2018;264:94–8.

    Article  CAS  Google Scholar 

  20. Butterly CR, Phillips LA, Wiltshire JL, Franks AE, Armstrong RD, Chen D, et al. Long-term effects of elevated CO2 on carbon and nitrogen functional capacity of microbial communities in three contrasting soils. Soil Biol Biochem. 2016;97:157–67.

    Article  CAS  Google Scholar 

  21. Singh P, Kumar V, Agrawal S. Evaluation of phytase producing bacteria for their plant growth promoting activities. Int J Microbiol. 2014;426483.

  22. Koide RT, Kabir Z. Extraradical hyphae of the mycorrhizal fungus Glomus intraradices can hydrolyse organic phosphate. New phytol. 2020;148:511–7.

    Article  Google Scholar 

  23. Hinsinger P. Bioavailability of soil inorganic P in the rhizosphere as affected by root-induced chemical changes: A review. Plant Soil. 2001;237:173–95.

    Article  CAS  Google Scholar 

  24. Dhillon J, Torres G, Driver E, Figueiredo B, Raun W. World phosphorus use efficiency in cereal crops. Agron J. 2017;109:1670–7.

    Article  CAS  Google Scholar 

  25. Washburne AD, Silverman JD, Morton JT, Becker DJ, Crowley D, Mukherjee S, et al. Phylofactorization: a graph partitioning algorithm to identify phylogenetic scales of ecological data. Ecol Monogr. 2019;89:e01353.

    Article  Google Scholar 

  26. Romano AH, Conway T. Evolution of carbohydrate metabolic pathways. Res Microbiol. 1996;147:448–55.

    Article  CAS  PubMed  Google Scholar 

  27. Kruger NJ, von Schaewen A. The oxidative pentose phosphate pathway: structure and organisation. Curr Opin Plant Biol. 2003;6:236–46.

    Article  CAS  PubMed  Google Scholar 

  28. Chun H, Keleş S. Sparse partial least squares regression for simultaneous dimension reduction and variable selection. J R Statist Soc B. 2010;72:3–25.

    Article  Google Scholar 

  29. Drigo B, Kowalchuk GA, Knapp BA, Pijl AS, Boschker HTS, van Veen JA. Impacts of 3 years of elevated atmospheric CO2 on rhizosphere carbon flow and microbial community dynamics. Glob Chang Biol. 2013;19:621–36.

    Article  PubMed  Google Scholar 

  30. Da Costa DP, Dias ACF, Cotta SR, Vilela D, de Andrade PAM, Pellizari VH, et al. Changes of bacterial communities in the rhizosphere of sugarcane under elevated concentration of atmospheric CO2. Glob Chang Biol. 2018;10:137–45.

    Article  Google Scholar 

  31. Wang Y, Yu Z, Li Y, Wang G, Tang C, Liu X, et al. Elevated CO2 alters the structure of the bacterial community assimilating plant-derived carbon in the rhizosphere of soya bean. Eur J Soil Sci. 2019;70:1212–20.

    CAS  Google Scholar 

  32. Wieczorek AS, Schmidt O, Chatzinotas A, von Bergen M, Gorissen A, Kolb S. Ecological functions of agricultural soil bacteria and microeukaryotes in chitin degradation: a case study. Front Microbiol. 2019;10:1293.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Hanada S, Sekiguchi Y. The phylum gemmatimonadetes. In: Rosenberg E, Delong E, Lory S, Stackebrandt E, Thompson F. (eds). The Prokaryotes: Other Major Lineages of Bacteria and the Archaea. (Heidelberg Springer, Berlin, 2014).

  34. Bailey VL, Fansler SJ, Stegen JC, McCue LA. Linking microbial community structure to β-glucosidic function in soil aggregates. ISME J. 2013;7:2044–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Spohn M, Kuzyakov Y. Distribution of microbial- and root-derived phosphatase activities in the rhizosphere depending on P availability and C allocation – coupling soil zymography with 14C imaging. Soil Biol Biochem. 2013;67:106–13.

    Article  CAS  Google Scholar 

  36. Huang L, Chen D, Zhang H, Song Y, Chen H, Tang M. Funneliformis mosseae enhances root development and Pb phytostabilization in Robinia pseudoacacia in Pb-contaminated soil. Front Microbiol. 2019;12:2591.

    Article  Google Scholar 

  37. Köhler J, Yang N, Pena R, Raghavan V, Polle A, Meier IC. Ecomycorrhizal fungao diversity increases phosphorus uptake efficiency of European beech. New Phytolog. 2018;220:1200–10.

    Article  Google Scholar 

  38. Zhang L, Geng G, Declerck S. Signal beyond nutrient, fructose, exuded by an arbuscular mycorrhizal fungus triggers phytate mineralization by an phosphate solubilizing bacterium. ISME J. 2018;12:2339–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Bell A (eds). Dung fungi: an illustrated guide to coprophilous fungi in New Zealand. (Victoria University Press, Wellington, New Zealand, 1983).

  40. Ali A, Ghani MI, Ding H, Iqbal M, Cheng Z, Cai Z. Garlic substrate induces cucumber growth development and decreases Fusarium wilt through regulation of soil microbial community structure and diversity in replanted disturbed soil. Int J Mol Sci. 2020;21:6008.

    Article  CAS  PubMed Central  Google Scholar 

  41. Isbell RF. The Australian soil classification. (CSIRO Publishing, Melbourne, Victoria); 1996.

    Google Scholar 

  42. McCaskill MR, Mitchell ML, Zollinger R, Armstrong RD, Partington D. Dry matter and nutritive value responses of native, naturalised and sown pasture species to soil Olsen P. Crop Pasture Sci. 2019;70:1097–109.

    Article  CAS  Google Scholar 

  43. Youssef RA, Chino M. Development of a new rhizobox system to study the nutrient status in the rhizosphere. Soil Sci Plant Nutr. 1988;34:461–5.

    Article  Google Scholar 

  44. Subedi KD, Ma BL, Liang BC. New method to estimate root biomass in soil through root-derived carbon. Soil Biol Biochem. 2006;38:2212–8.

    Article  CAS  Google Scholar 

  45. Tang C, Robson AD, Dilworth MJ. A split-root experiment shows that iron is required for nodule initiation in Lupinus angustifolius L. New Phytol. 1990;115:61–7.

    Article  CAS  Google Scholar 

  46. Motomizu S, Wakimoto T, Toei K. Spectrophotometric determination of phosphate in river waters with molybdate and malachite green. Analyst. 1983;108:361–7.

    Article  CAS  Google Scholar 

  47. Yuen SH, Pollard AG. Determination of nitrogen in agricultural materials by the Nessler Reagent. II. Micro-determinations in plant tissue and in soil extracts. J Sci Food Agr. 1954;5:364–9.

    Article  CAS  Google Scholar 

  48. Olsen SR, Dean LA. Phosphorus. In: Black CA, editor. Methods of soil analysis. (Soil Science Society of America, Wisconsin); 1965. p. 1044–7.

    Google Scholar 

  49. Zhang L, Fan J, Ding X, He X, Zhang F, Feng G. Hyphosphere interactions between an arbuscular mycorrhizal fungus and a phosphate solubilizing bacterium promote phytate mineralization in soil. Soil Biol Biochem. 2014;74:177–83.

    Article  CAS  Google Scholar 

  50. Stout L, Nguyen TT, Jaisi DP. Relationship of phytate, phytate-mineralizing bacteria, and beta-propeller phytase genes along a coastal tributary to the chesapeake bay. Soil Sci Soc Am J. 2016;80:84–96.

    Article  CAS  Google Scholar 

  51. Motomizu S, Wakimoto T, Toei K. Spectrophotometric determination of phosphate in river waters with molybdite and malachite green. Analyst. 1980;108:361–7.

    Article  Google Scholar 

  52. Miranda KM, Espey MG, Wink DA. A rapid, simple spectrophotometric method for simultaneous detection of nitrate and nitrite. Nitric Oxide. 2001;5:62–71.

    Article  CAS  PubMed  Google Scholar 

  53. Vance ED, Brookes PC, Jenkinson DS. An extraction method for measuring soil microbial biomass C. Soil Biol Biochem. 1987;19:703–7.

    Article  CAS  Google Scholar 

  54. Domanski G, Kuzyakov Y, Siniakina SV, Stahr K. Carbon flows in the rhizosphere of ryegrass (Lolium perenne). J Plant Nutr Soil Sc. 2001;164:381–7.

    Article  CAS  Google Scholar 

  55. Brookes PC, Landman A, Pruden G, Jenkinson D. Chloroform fumigation and the release of soil nitrogen: a rapid direct extraction method to measure microbial biomass nitrogen in soil. Soil Biol Biochem. 1985;17:837–42.

    Article  CAS  Google Scholar 

  56. Zibilske LM. Carbon mineralization. In: RW Waever, S Angle, P Bottomley, D Bezdicek, S Smith, A Tabatabai, A Wollum (eds). Methods of soil analysis part 2: microbiological and chemical properties. (Soil Science Society of America, Madison, WI, 1994.) pp 835– 863.

  57. Rukshana F, Butterly CR, Baldock JA, Xu JM, Tang C. Model organic compounds differ in priming effects on alkalinity release in soils through carbon and nitrogen mineralisation. Soil Biol Biochem. 2012;51:35–43.

    Article  CAS  Google Scholar 

  58. Apprill A, McNally S, Parsons R, Weber L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat Microb Ecol. 2015;75:129–37.

    Article  Google Scholar 

  59. Blaalid R, Kumar S, Nilsson RH, Abarenkov K, Kirk PM, Kauserud H. ITS1 versus ITS2 as DNA metabarcodes for fungi. Mol Ecol Resour. 2013;13:218–24.

    Article  CAS  PubMed  Google Scholar 

  60. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. PNAS. 2011;108:4516–22.

    Article  CAS  PubMed  Google Scholar 

  61. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011;17:10–2.

    Article  Google Scholar 

  63. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Karp PD, Riley M, Paley SM, Pellegrini-Toole A. The MetaCyc database. Nucleic Acids Res. 2002;30:59–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Douglas GM, Maffei VJ, Zaneveld J, Yurgel SN, Brown JR, Taylor CM, et al. PICRUSt2: an improved and extensible approach for metagenome inference. BioRxiv. 2019.

  66. Steel RG, Torrie JH. (eds) Principles and procedures of statistics: a biometrical approach. 2nd ed. McGraw-Hill. New York. .

  67. McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of micdrobiome census data. Plos One. 2013;8:e61217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Borcard D, Gillet F, Legendre P. (eds). Numerical ecology with R. (Springer, New York, 2011).

  69. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community ecology package. R package version 24−1 Available at:

  70. Lin H, Peddada SD. Analysis of compositions of microbiomes with bias correction. Nat Commun. 2020;11:3514.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


We appreciate that Dr. Rob Norton at the International Plant Nutrition Institute (IPNI) kindly provided the soil from a long-term experiment at Dahlen, Victoria, Australia


This research was supported by La Trobe University Research Focus Area project (RFA) (2000002982) and Australian Research Council Project (DP210100775).

Author information

Authors and Affiliations



J.J., C.T., and A.F. conceived the idea and designed the research. J.J. and M.M. collected the samples. J.J., C.K., S.B., and X.W. conducted the experiments. J.J., C.K., J.W., S.P., S.B., and Z.X. analyzed the data. J.J. and C.K. wrote the first draft of the manuscript, and the authors further reviewed and edited the manuscript. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jian Jin or Caixian Tang.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1

The biomass of (A), and P contents (B), and 13C atom‰ excess in shoot and root (C), and soil 13C atom‰ excess in the rhizosphere compartments of 1.5 and 3 mm away from the root growth zone (D). Wheat plants were grown in Chromosol and Vertosol for 10 weeks under elevated CO2 (eCO2, 800 ppm) and ambient CO2 (aCO2, 400 ppm). Error bars are standard errors (n=6). *, ** and *** indicate significant differences between aCO2 and eCO2 treatments at p < 0.05, p < 0.01 and p < 0.001, respectively. The significance levels of the main effects of CO2 and rhizosphere (Rhizo) and their interactions on 13C atom‰ excess in Chromosol and Vertosol are presented. Figure S2 Boxplots of species richness, Shannon indices and Pielou evenness in the rhizosphere compartments at 0 mm, 1.5 mm and 3 mm, and in the bulk soil of Chromosol and Vertosol. Wheat plants were grown in Chromosol and Vertosol for 10 weeks under elevated CO2 (eCO2, 800 ppm) and ambient CO2 (aCO2, 400 ppm). Figure S3 The effects of elevated CO2 (A, C) and rhizosphere (B, D) on the relative abundances of bacterial phyla in Chromosol (A, B) and Vertosol (C, D). Wheat plants were grown for 10 weeks in a rhizobox comprising rhizosphere compartments of 0, 1.5 and 3 mm away from the root growth zone and the bulk soil under elevated CO2 (eCO2, 800 ppm) and ambient CO2 (aCO2, 400 ppm). Figure S4 Relative abundances (centred log-ratios) of selected microbial metabolic pathways in the rhizosphere of wheat. Plants were grown in Chromosol and Vertosol under elevated CO2 (eCO2, 800 ppm) and ambient CO2 (aCO2, 400 ppm) for 10 weeks. ns, *, ** and *** indicate significance of two-sample Wilcoxon Mann-Whitney tests at p > 0.05, p < 0.05, p < 0.01 and p < 0.001, respectively. Figure S5 Scatterplots of mineralized phytate (μg g-1) to species richness, Shannon indices and Pielou evenness in the rhizosphere compartments of 0 mm, 1.5 mm and 3 mm, and the bulk soil of Chromosol and Vertosol. Figure S6 The photos showing (A) plants grown for 10 weeks under ambient CO2 (aCO2, 400 ppm) and elevated CO2 (eCO2, 800 ppm), and (B) a rhizo-compartment, and (C) a schematic structural diagram of the rhizobox in the 13C-labelling device used in the experiment. Table S1 The concentrations of total organic C, total P and Olsen P, soil pH and soil particle composition (texture) of the Chromosol and Vertosol used in the study

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jin, J., Krohn, C., Franks, A.E. et al. Elevated atmospheric CO2 alters the microbial community composition and metabolic potential to mineralize organic phosphorus in the rhizosphere of wheat. Microbiome 10, 12 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: