Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 08 January 2026

Sec. Terrestrial Microbiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1698121

Effects of tree mycorrhizal dominance on soil microbial community structure and microbial nutrient limitation

Yajie XuYajie Xu1Longfei Hao
Longfei Hao1*Yongjie Yue
Yongjie Yue1*Runhong GaoRunhong Gao1Lingze ZhangLingze Zhang1Kai ZhaoKai Zhao2Zhenghui ZhaoZhenghui Zhao1
  • 1College of Forestry, Inner Mongolia Agricultural University, Hohhot, China
  • 2College of Desert Control Science and Engineering, Inner Mongolia Agricultural University, Hohhot, China

Introduction: Mycorrhizal fungi play a central role in nutrient cycling in forest ecosystems. The functional differences between arbuscular mycorrhizal (AM) and ectomycorrhizal (EcM) tree species significantly affect soil microbial community structure and patterns of microbial nutrient limitation, with substantial implications for ecosystem stability and biogeochemical cycling under global changes. However, the regulatory mechanisms of different dominant mycorrhizal tree species and their mixed mycorrhizal configurations on microbial nutrient limitation remain unclear.

Methods: This study investigated a typical AM tree species (Ulmus pumila) and an EcM tree species (Pinus sylvestris var. mongolica) in the southern Horqin Sandy Land. We compared rhizosphere soil nutrient status, extracellular enzyme activities, and microbial community structure among pure U. pumila stands, pure P. sylvestris var. mongolica stands, and U. pumilaP. sylvestris mixed stands.

Results: Results showed that AM pure forests exhibited extremely high C-acquiring enzyme activities but reduced soil nutrient content and microbial biomass, maintaining higher bacterial diversity. By contrast, the activities of N- and P-acquiring enzymes, soil nutrient contents, and fungal diversity in the EcM pure stands were significantly higher than those in the AM pure stands (p < 0.05). Mixed forests improved soil nutrient status through complementary mixed-mycorrhizal strategies, promoted microbial biomass accumulation, and modulated extracellular enzyme activities, thereby significantly alleviating microbial nutrient limitations. Moreover, their microbial networks exhibited greater complexity and stability than pure stands. Structural equation modeling further revealed that tree mycorrhizal dominance and microbial biomass were the primary factors alleviating microbial nutrient limitation: microbial biomass and mycorrhizal dominance showed significant negative effects on vector length (path coefficients −0.71 and 0.33, p < 0.01), whereas mycorrhizal dominance exerted a highly significant positive effect on vector angle (path coefficient 0.70, p < 0.001).

Discussion: In conclusion, mixed mycorrhizal strategies alleviate microbial nutrient limitations by enhancing soil nutrient status, microbial community structure, and extracellular enzyme activities, providing theoretical support for ecosystem restoration and sustainable development in arid regions.

1 Introduction

Almost all tree species are associated with either arbuscular mycorrhizal (AM) fungi or ectomycorrhizal (EcM) fungi (Phillips et al., 2013). These fungi form symbiotic relationships with trees, exchanging carbon (C) as a nutrient source and protecting trees from environmental stresses (Brundrett and Tedersoo, 2018). Tree species associated with AM and EcM fungi differ in litter quality, root and hyphal traits, and nutrient uptake efficiency, which strongly influence biogeochemical processes (C and nutrient cycling) in ecosystems (Zheng et al., 2023; Bahram et al., 2020; Tedersoo and Bahram, 2019). In general, EcM-associated trees are characterized by lower-quality litter, slower decomposition, higher forest floor C storage, and elevated soil C/N ratios compared with AM-associated trees (Craig et al., 2018; Jo et al., 2019). According to the mycorrhiza-associated nutrient economy framework, AM trees are linked to an inorganic nutrient economy with rapid nutrient cycling, whereas EcM trees are linked to an organic nutrient economy with slow turnover of plant-derived C and strong root–microbe interactions (Yu et al., 2023). Differences in mycorrhizal type are also consistent with pronounced shifts in soil microbial composition, activity, and function, which can alter soil organic carbon (SOC) and nutrient turnover (Heděnec et al., 2020; Tang et al., 2024a). Currently, most studies remain focused on the ecological effects of single dominant mycorrhizal types, whereas systematic understanding of the role of AM–EcM mixed forests in regulating soil biogeochemical processes and microbial nutrient limitations is still lacking. In recent years, a gradual shift in forest management from monoculture plantations to mixed-species stands has emerged as an important strategy to accelerate the restoration of degraded forests and enhance ecosystem services (Deng et al., 2020). Studies suggest that litter from AM-associated trees may accelerate the decomposition of needle litter by adding labile nutrients and stimulating microbial metabolism, thereby improving soil ecosystem stability (Polyakova and Billor, 2007; Zhang et al., 2020). Tang et al. (2024b) further demonstrated that introducing native broadleaf trees into coniferous plantations promotes biogeochemical cycling in Masson pine forests and strengthens ecosystem functions. Thus, combining different mycorrhiza-associated tree species may facilitate efficient resource use, enhancing forest ecosystem stability.

Soil microbes regulate C and nutrient transformations and cycling through the secretion of extracellular enzymes, thereby obtaining energy and nutrients to meet their metabolic requirements. Moreover, shifts in microbial extracellular enzyme activities significantly affect microbial metabolic functions, nutrient demands, and stoichiometric balance (Zhou et al., 2020; Cui et al., 2021). Previous studies have shown that β-1,4-glucosidase (BG), β-1,4-N-acetylglucosaminidase (NAG), leucine aminopeptidase (LAP), and acid or alkaline phosphatase (AP) serve as indicators of microbial C, N, and P demands, respectively (Sinsabaugh et al., 2009; Liu et al., 2023; Feng et al., 2019). Ecoenzymatic stoichiometry connects microbial nutrient acquisition strategies with soil nutrient availability, reflecting the supply–demand balance of nutrients for soils and microbes. Over the past decades, enzyme activity and ecoenzymatic stoichiometry, including vector analysis, have been widely applied to assess microbial metabolic activities and resource limitations in ecosystems (Wang et al., 2023). For example, Moorhead et al. (2013) proposed calculating vector length and angle based on (BG/[BG + ALP]) and (BG/[BG + NAG + LAP]) to quantify relative C versus nutrient acquisition (vector length) and P versus N acquisition (vector angle). Mycorrhizal dominance of tree species (i.e., pure stands dominated by one mycorrhizal type vs. mixed stands composed of AM- and EcM-associated trees) can lead to distinct soil nutrient regimes and microbial metabolic strategies, triggering specific resource acquisition strategies (Luo et al., 2023). Therefore, exploring the relationships between mycorrhizal dominance and microbially mediated soil processes is critical for optimizing stand structure, enhancing forest ecosystem functions, and responding to global environmental changes.

Semi-arid regions account for 37% of China’s land area and play a critical role in national ecological security and ecological civilization. Aohan County, Chifeng City, Inner Mongolia, is a typical semi-arid region known as China’s “first county of afforestation,” with a forest cover of 44.07%. However, local forests face widespread soil fertility decline and reduced ecosystem functions due to prolonged drought, extreme climate, nutrient-poor soils, and long-term monoculture plantations. In light of the ongoing escalation of global climate change and regional ecological pressure, the stability and adaptability of artificial forests must be urgently improved through scientific means. Thus, a thorough understanding of mycorrhizal dominance and microbially mediated nutrient cycling is essential for restoring soil fertility in degraded plantations, optimizing stand composition, and improving ecological functions and sustainable forest management. The present study focused on the typical AM-associated tree species (Ulmus pumila) and EcM-associated tree species (Pinus sylvestris var. mongolica), as well as their mixed stands, at the southern edge of the Horqin Sandy Land. We analyzed the rhizosphere soil microbial community structure and nutrient limitation patterns. We proposed the following hypotheses: (1) compared with EcM-associated trees, soils under AM-associated trees exhibit stronger microbial C limitation and greater bacterial dominance; and (2) resource complementarity strategies are present in AM + EcM mixed forests, which can improve the overall resource-use efficiency of the stands and mitigate microbial nutrient limitations.

2 Materials and methods

2.1 Overview of the study area

The study area is located in Aohan Banner, Chifeng City, Inner Mongolia Autonomous Region, at the southern margin of the Horqin Sandy Land (E 119°30′–120°53′, N 41°42′–43°02′). The region extends 176 km from north to south and 122 km from east to west, with a total land area of 829 km2. The climate is mainly continental monsoon type, with a large temperature difference between day and night. Spring is dry and windy with rapid warming; summer is hot with concentrated rainfall; autumn is marked by a sharp temperature drop; and winter is cold. The mean annual precipitation is approximately 351.8 mm, with a mean annual temperature of 8 °C. The extreme minimum and maximum temperatures are −31 °C and 39.7 °C, respectively, with a frost-free period of 130–150 days. The accumulated annual temperature ranges from 2,600 °C to 3,200 °C. The dominant soil types are brown soil, chestnut soil, and aeolian sandy soil, among which chestnut soil and sandy soil account for 52.6% of the total land area. The region is mainly covered by plantations, with P. sylvestris var. mongolica, U. pumila, and Populus spp. as the dominant species. Prunus sibirica and Pinus tabuliformis are also present (Figure 1).

Figure 1
Topographic map of Chifeng City, Inner Mongolia, showing elevation with a color gradient from green (256 meters) to red (2054 meters). The inset map highlights Chifeng City's location. A black dot marks the study site. The map includes administrative boundaries and scale bars for reference.

Figure 1. Location of the study area.

2.2 Site selection and sample collection

We selected three forest stand types planted in 2008 at the Shuangjing Forest Farm in Aohan Banner, Inner Mongolia: P. sylvestris var. mongolica pure stands (EcM-associated tree species), U. pumila pure stands (AM-associated tree species), and mixed stands in which the two species were planted at a 1:1 ratio. To ensure comparability, we established plots under similar site conditions (elevation, slope, and aspect). For each stand type, three replicate plots (50 m × 50 m) were established, resulting in nine plots in total, with an inter-plot distance of approximately 1,000 m. Soil samples were collected in late September 2024, ensuring that no rainfall had occurred for at least seven consecutive days prior to sampling. Using a tree-by-tree inventory, we selected 15 standard trees from each pure stand plot. In the mixed stand, 15 P. sylvestris var. mongolica and 15 U. pumila individuals were selected as standard trees in each plot. We took the standard trees as the center and collected rhizosphere soil samples within a horizontal radius of 1 m and a vertical soil layer of 0–20 cm by the shaking method. This soil layer represents the zone with the highest fine-root density and strongest rhizosphere effects (Cleveland and Liptzin, 2007), and it also exhibits the greatest extracellular enzyme activity and microbial activity, making it the most sensitive to microbial resource limitation (Saiya-Cork et al., 2002). Subsequently, the soil samples obtained from every 5 standard trees within the same standard plot were evenly mixed to obtain 1 soil sample. For each standard plot, 15 standard trees of the same tree species were mixed to form 3 soil samples. A total of 9 soil samples were obtained from each treatment. Soil samples from P. sylvestris var. mongolica pure stands, U. pumila pure stands, and from P. sylvestris var. mongolica and U. pumila in the mixed stands were labeled Ps.P, Up.P, Ps.M, and Up.M, respectively. Samples were sieved (2 mm) to remove roots and stones and then divided into three subsamples: one stored at 4 °C for microbial biomass, enzyme activity, and available nutrient analyses; one frozen at −80 °C for microbial community analysis; and one air-dried for soil total nutrients and organic matter determination.

2.3 Soil nutrient determination

The desktop electrode pH meter (Sartorius Basic pH Meter PB-10) was used to measure the pH of the soil at a ratio of 1:25 (w/v) with 0.01 M CaCl2 solution. Total carbon (TC) and total nitrogen (TN) were determined with a Vario EL elemental analyzer (Elementar, Hanau, Germany) (Zeng et al., 2010), whereas total phosphorus (TP) was measured using the molybdenum–antimony anti-spectrophotometric method (Tighe et al., 2018). SOC, available nitrogen (AN), and available phosphorus (AP) were measured by the potassium dichromate oxidation method, alkaline hydrolysis diffusion, and Olsen method, respectively (Mebius, 1960; Cornfield, 1960; Page et al., 1982).

2.4 Measurement of microbial biomass

Microbial biomass carbon (MBC), microbial biomass nitrogen (MBN), and microbial biomass phosphorus (MBP) were determined using the chloroform fumigation–extraction method (Shu et al., 2023). The C and N concentrations in K2SO4extracts from fumigated and non-fumigated samples were analyzed with a total organic carbon analyzer equipped with an N unit (TOC-V and TNM-1, Shimadzu, Japan). The P concentration in NaHCO3 extracts was determined with a continuous flow analyzer (BRAN+LUEBBE AA3, Germany). Conversion factors of kEC = 0.45, kEN = 0.54, and kEP = 0.40 were applied for MBC, MBN, and MBP, respectively (Wu et al., 1990; Brookes et al., 1985).

2.5 Determination of soil enzyme activity and quantification of microbial nutrient limitation

We selected four extracellular enzymes that represent the cycling of C, N, and P substances: BG, NAG, LAP, and ALP. Among them, BG is a C cycle-related enzyme, NAG and LAP are N cycle-related enzymes, and ALP is a P cycle-related enzyme. Enzyme assays followed the method of Saiya-Cork et al. (2002). A vector analysis was used to quantify the microbial metabolic limitation based on untransformed enzymatic activities (Xu et al., 2022). The length of the vector represents the microbial carbon limitation, and the microbial C limitation increases with the increase in vector length. An angle of the vector lower than 45° indicates microbial N limitation, whereas an angle of the vector higher than 45° indicates microbial P limitation. The microbial N limitation decreases with the increase in vector angle, and the microbial phosphorus limitation increases with the increase in vector angle. The vector length and angle can be calculated as follows (Moorhead et al., 2013):

Vector length = ( ln ( BG ) / ln ( NAG + LAP ) ) 2 + ( ln ( BG ) / ln ( AP ) ) 2
Vector length = DEGREES { ATAN 2 [ ( ln ( BG ) / ln ( AP ) ) , ( ln ( BG ) / ln ( NAG + LAP ) ) ] }

2.6 Determination of microbial community composition

According to the manufacturer’s instructions, the microbial community gene DNA was extracted from the soil samples using HiPure Soil DNA kits. The ITS2 region of ITS and the V3 + V4 region of 16S rDNA were amplified using specific primers with barcodes. The primer sequences were as follows: ITS3_KYO2: GATGAAGAACGYAGYRAA; ITS4: TCCTCCGCTTATTGATATGC and 341F: CCTACGGGNGGCWGCAG; and 806R: GGACTACHVGGGTATCTAAT. The amplified products (i.e., amplicons) were ligated with sequencing adapters, and a sequencing library was constructed. Sequencing was performed on the Illumina platform. The operational taxonomic units (OTUs) with 97% similarity were clustered using the UPARSE method (Edgar, 2013), and chimeric sequences were identified and removed. The OTUs obtained were checked using the MaarjAM database1 (Öpik et al., 2010), and unclassified OTUs were manually checked and filtered in the NCBI database (Jiang et al., 2018). The sequencing depth of the fungal ITS data in this study ranged from 67,329 to 78,875 valid sequences per sample, with an average of 73,428 sequences per sample; the sequencing depth of the bacterial 16S data ranged from 114,221 to 130,978 valid sequences per sample, with an average of 121,412 sequences per sample. These sequence data also served as the basis for the subsequent construction of microbial co-occurrence networks and the calculation of network topological properties. To eliminate the influence of sequencing depth differences on the α diversity analysis, the ITS and 16S data sets were, respectively, subjected to sparsification to achieve the minimum sequencing depth that could retain all samples in each dataset. The sparsified OTU table was used to calculate the Observed richness, Chao1, Shannon–Wiener, and Simpson indices.

2.7 Statistical analysis

Data analyses were performed using DPS (version 9.01), and data visualization was conducted in Origin (version 9.4). One-way analysis of variance (ANOVA) followed by Fisher’s least significant difference (LSD) test at a 5% significance level was used to evaluate rhizosphere soil physicochemical properties, extracellular enzyme activities, microbial biomass, and fungal and bacterial diversity. Principal component analysis (PCA) of soil extracellular enzyme activities was conducted using Origin (version 9.4). OTUs with a mean relative abundance <0.01% and occurring in less than 50% of the total samples were filtered out in R (version 4.2.1) Spearman’s rank correlations among OTUs were calculated using the “Hmisc package, and only strong (|r| > 0.9) and significant (p < 0.05) associations were retained for constructing the co-occurrence network. Undirected weighted microbial co-occurrence networks were then constructed using the “igraph” and “ggClusterNet” packages, and network topological properties (including node number, edge number, average degree, average path length, clustering coefficient, and modularity) were calculated. The ecological interpretations of these topological metrics are as follows: the numbers of nodes and edges reflect the overall network size and the richness of potential interactions; the average degree represents the overall connectivity among taxa; the average path length indicates the compactness of the network and the potential efficiency of interactions; the clustering coefficient describes local connectivity and the aggregation of taxa; and modularity reflects the extent to which the network is organized into ecological niches or functional guilds (Ma et al., 2020). Network visualization was performed using Gephi v0.9.2 (Ding et al., 2025). The “mgcv package in R (version 4.2.1) was used to fit generalized additive models (GAMs) to examine the relationships between microbial C and N limitations and soil pH, soil nutrients, and microbial biomass. Partial least squares path modeling (PLS-PM) was performed using the “plspm” package in R (version 4.2.1) to assess pathways influencing microbial resource limitation, including mycorrhizal dominance (i.e., rhizosphere soils of Pinus sylvestris var. mongolica monoculture, Ulmus pumila monoculture, and both species in the mixed stand), pH, soil nutrients (i.e., total carbon, total nitrogen, total phosphorus, organic carbon, available nitrogen, and available phosphorus), fungal diversity (i.e., fungal observed richness and Chao1 index), bacterial diversity (i.e., bacterial observed richness and Chao1 index), and microbial biomass (i.e., microbial biomass C, N, and P). To ensure model robustness, variables with loading values below 0.6—including the fungal Shannon–Wiener index, fungal Simpson index, bacterial Shannon–Wiener index, and bacterial Simpson index—were excluded (Liu et al., 2024).

3 Results

3.1 Influence of tree mycorrhizal dominance on soil chemical properties

The dominance of tree mycorrhizae had a significant impact on the chemical properties of the rhizosphere soil (p < 0.05, Table 1). TC, TN, TP, SOC, AN, and AP all reached their lowest values at Up.P, suggesting significant differences compared with other treatments (p < 0.05). By contrast, Up.M showed the highest TC (9.85 g kg−1), TN (0.58 g kg−1), SOC (4.05 g kg−1), and AN (114.91 mg g−1), indicating that the mixed forest significantly increased the organic carbon input and nitrogen accumulation levels in the AM-related rhizosphere soil, thereby promoting soil fertility. AP reached the highest value at Ps.P (2.73 mg g−1), with significant differences observed across each treatment (p < 0.05), suggesting that the rhizosphere soil of exogenous mycorrhizal tree species have certain advantages in obtaining and maintaining P. The mixed forest effectively improved the carbon and nitrogen nutrient conditions of the rhizosphere soil through the combination of AM and EcM tree species, and it exerted a certain regulatory effect on soil pH, facilitating the availability of soil nutrients and microbial functions.

Table 1
www.frontiersin.org

Table 1. Soil physicochemical properties under different dominant mycorrhizal forests.

3.2 Impacts of tree mycorrhizal dominance on soil extracellular enzyme activity and microbial resource limitation

The activity of C-acquiring enzymes was highest in Up.P (29.60 μmol h−1), which was 1.94 and 1.92 times that of Ps.M and Up.M, respectively (p < 0.05, Figure 2A). N-acquiring enzymes reached the highest activity in Ps.P (71.49 μmol h−1) (p < 0.05, Figure 2B). The activity of P-acquiring enzymes was elevated in Ps.M and Up.M soils at 43.37 and 43.57 μmol h−1, respectively (p < 0.05, Figure 2C). Mixed AM and EcM forests enhance the utilization efficiency of C resources through complementary C demands but also intensify microbial competition for P. Principal component analysis of extracellular enzyme activities revealed the highest microbial C limitation in Up.P (Figure 2D). The distribution of data points further indicated that Up.P experienced stronger microbial C and N limitations than Ps.P, whereas microbial resource limitations were alleviated in AM–EcM mixed stands. Vector length and vector angle, calculated from the relative proportions of extracellular enzyme activities, quantified microbial C and N limitations. All treatments had vector angles below 45°, suggesting that microbes were primarily limited by N rather than P. Compared with Up.P, vector length and angle in Up.M decreased from 1.34 to 0.98 and increased from 39.66° to 42.02°, respectively; compared with Ps.P, vector length and angle in Ps.M decreased from 1.08 to 0.97 and increased from 40.00° to 42.91°, respectively, with vector angles shifting toward 45° (Figures 2EG). Compared with the pure stands, the vector angles in the mixed stands shifted toward 45°, accompanied by reduced vector lengths. This pattern indicates that microbial N limitation was alleviated, the N: P supply–demand balance became more equilibrated, and microbial C limitation was also weakened under mixed stands. It is important to note that all vector angles remained <45°, demonstrating that the mixed stands did not shift to P limitation. Instead, they transitioned from a strongly N-limited state toward a more balanced N–P stoichiometric condition.

Figure 2
Box and scatter plots show various enzyme activities and vector analyses for Ps.P, Up.P, Ps.M, and Up.M groups. Panels A to C present box plots of BG, NAG+LAP, and ALP activities, respectively, with statistical details. Panel D displays PCA biplot with enzymes BG, NAG, LAP, and ALP. Panel E shows a scatter plot of BG/(BG+ALP) against BG/(BG+NAG+LAP) with vector illustrations. Panels F and G depict box plots of vector length and angle, respectively, with significant differences among groups.

Figure 2. Effects of different dominant mycorrhizal types on soil enzyme activities and microbial resource limitations. (A) C-acquiring enzyme: β-1,4-glucosidase (BG); (B) N-acquiring enzymes: β-1,4-N-acetylglucosaminidase (NAG), leucine aminopeptidase (LAP); (C) P-acquiring enzyme: alkaline phosphatase (ALP); (D) Principal component analysis of soil extracellular enzyme activities; (E) Enzymatic stoichiometry of the relative proportions of C to N acquisition versus C to P acquisition; (F) Vector length; and (G) Vector angle. Ps.P, Pinus sylvestris pure stand; Up.P, Ulmus pumila pure stand; Ps.M, Pinus sylvestris in mixed stand; Up.M, Ulmus pumila in mixed stand. Significant differences (p < 0.05) between treatments are indicated by different lowercase letters.

3.3 Effects of tree mycorrhizal dominance on microbial biomass and community composition

MBC, MBN, and MBP reached the highest values in Up.M, with 158.91, 7.52, and 59.53 mg g−1, respectively, which were significantly higher than in other stand types (p < 0.05, Figure 3). By contrast, the lowest values of MBC, MBN, and MBP were observed in Up.P (89.13, 4.81, and 40.25 mg g−1, respectively). Mixed AM–EcM stands significantly enhanced microbial biomass in the rhizosphere soils compared with AM or EcM monocultures (p < 0.05). Fungal community diversity was higher in Ps.P than in Up.P but lower than in Ps.M and Up.M (p < 0.05, Figures 4AD). Conversely, bacterial community diversity was highest in Up.P compared with Ps.P, Ps.M, and Up.M (p < 0.05, Figures 4EH). Network analysis revealed that the rhizosphere microbial network of Up.P had the lowest values in terms of node number, edge number, average degree, average path length, and average clustering coefficient (528, 1,002, 3.79, 8.82, and 0.38, respectively), compared with Ps.P (530, 1,012, 3.82, 8.94, and 0.42, respectively) (Figure 5, Table 2). Among the four stand types, Up.M showed the most complex rhizosphere microbial network, with the highest node number, edge number, average degree, and clustering coefficient (542, 1,164, 4.30, and 0.43, respectively), followed by Ps.M (530, 1,025, 3.86, and 0.43), and it exhibited the highest modularity (0.73). These results suggest that the combination of AM and EcM tree species enhances the complexity and stability of co-occurrence networks within rhizosphere microbial communities.

Figure 3
Three box plots labeled A, B, and C compare MBC, MBN, and MBP levels, respectively, across four treatments: Ps.P, Up.P, Ps.M, and Up.M. Each plot shows significant differences in values with the highest in Up.M and the lowest in Up.P. The F-values and p-values indicate significant differences in all plots (P < 0.0001).

Figure 3. Soil microbial biomass in different dominant mycorrhizal forest types. (A) Microbial biomass C (MBC). (B) Microbial biomass N (MBN). (C) Microbial biomass P (MBP). Significant differences (p < 0.05) between treatments are indicated by different lowercase letters.

Figure 4
Bar charts comparing microbial diversity indices under different conditions for fungi and bacteria. Panels E and H show significant differences in observed richness and Chao1 index for bacteria, with respective F and P-values. Other panels show less significant results. Conditions include Ps.P, Up.P, Ps.M, and Up.M.

Figure 4. Fungal and bacterial community diversity indices in soils of different dominant mycorrhizal forest types. (A) Fungus observed richness; (B) fungus Shannon–Wiener index; (C) fungus Simpson index; (D) fungus Chao1 index; (E) bacteria observed richness; (F) bacteria Shannon–Wiener index; (G) bacteria Simpson index; and (H) bacteria Chao1 index. Significant differences (p < 0.05) between treatments are indicated by different lowercase letters.

Figure 5
Four circular diagrams labeled A, B, C, and D display network graphs with pink and green nodes representing bacteria and fungi, respectively. Pink and green nodes are interconnected, illustrating their relationships. Each diagram depicts a different arrangement and density of nodes, with green nodes mostly clustered centrally. A legend at the bottom identifies the colors: pink for bacteria and green for fungi.

Figure 5. Microbial co-occurrence networks in the rhizosphere soils of forests dominated by different mycorrhizal types. The size of the dot represents the size of the microbial group. (A) Pinus sylvestris pure stand (Ps.P); (B) Ulmus pumila pure stand (Up.P); (C) Pinus sylvestris in mixed stand (Ps.M); and (D) Ulmus pumila in mixed stand (Up.M).

Table 2
www.frontiersin.org

Table 2. Topological parameters of microbial co-occurrence networks in forests dominated by different mycorrhizal types.

3.4 Relationships among soil chemical properties, microbial biomass, and microbial resource limitation

The GAM results showed that pH exerted a relatively weak effect on vector length and displayed an approximately linear increasing trend (Figure 6A). In contrast, soil nutrient content PC1 and microbial biomass had the highest explanatory power for vector length (R2 = 0.939 and 0.945, respectively) and exhibited strongly nonlinear responses (edf = 5.671 and 4.322, Figures 6B,C), indicating that they are the primary drivers of microbial C limitation. Changes in vector angle were mainly controlled by pH (R2 = 0.631, edf = 1D), whose influence was stronger than that of soil nutrient content PC1 (R2 = 0.589, Figure 6E) and microbial biomass (R2 = 0.544, Figure 6F).

Figure 6
Six scatter plots labeled A to F show relationships between variables. Data points are colored by category: green (Ps.P), pink (Up.P), orange (Ps.M), and red (Up.M). Lines with shaded confidence intervals indicate trends. Key metrics include R² values, p-values, and edf. Plot A relates vector length to pH, B to principal component one of nutrient content, C to total microbial biomass. Plot D relates vector angle to pH, E to principal component one of nutrient content, F to total microbial biomass. All plots highlight significant correlations (P < 0.001).

Figure 6. Microbial nutrient limitation in relation to soil resources under different mycorrhizal dominance: (A) Vector length: soil pH; (B) vector length: PC1 of soil nutrient content; (C) vector length: total microbial biomass; (D) vector angle: soil pH; (E) vector angle: PC1 of soil nutrient content; and (F) vector angle: total microbial biomass.

We further applied PLS-PM to evaluate the effects of dominant mycorrhizal type, soil nutrients, pH, fungal diversity, bacterial diversity, and microbial biomass on vector length (microbial C limitation) and vector angle (microbial N limitation), with the model achieving a goodness of fit index of 0.81. The model demonstrated that stand configuration significantly influenced soil nutrients, pH, and microbial processes through direct and indirect pathways, explaining 81% of microbial nutrient limitation (Figure 7A). The direct effects indicated that tree mycorrhizal dominance was a crucial factor of microbial nutrient limitation. Specifically, mycorrhizal dominance significantly decreased soil pH and vector length (path coefficients = −0.41 and −0.33, p < 0.01) but strongly increased soil nutrient content and vector angle (path coefficients = 0.61 and 0.70, p < 0.001). The total effects of PLS-PM further revealed that pH exerted an overall positive effect on vector length (path coefficient = 0.31) and a negative effect on vector angle (path coefficient = −0.02) by influencing nutrient content (path coefficient = −0.51, p < 0.001) and microbial biomass (path coefficient = −0.11). Soil nutrients and microbial biomass showed direct negative effects on vector length but positive effects on vector angle, with microbial biomass significantly reducing vector length (path coefficient = −0.77, p < 0.01, Figure 7B, C). These findings suggest that mixing AM- and EcM-associated tree species can effectively improve soil nutrient status and alleviate microbial carbon and nitrogen limitations. For clarity, Figure 7 displays only the direct path coefficients, whereas the total effects reported in the text represent the overall effects (i.e., the sum of direct and indirect influences).

Figure 7
Diagram (A) is a structural equation model illustrating relationships between mycorrhizal dominance, pH, soil nutrients, fungi, bacteria, soil microbial biomass, vector length, and vector angle, with paths demonstrating direct and indirect effects. Diagrams (B) and (C) are bar graphs showing vector length and vector angle for each variable, respectively, with positive and negative values compared across different factors.

Figure 7. Partial least squares path model (PLS-PM) illustrating the direct and indirect effects of tree mycorrhizal dominance, soil pH, soil nutrients, fungal diversity, bacterial diversity, and soil microbial biomass on vector length (microbial C limitation) and vector angle (microbial N limitation). Solid and dashed arrows represent positive and negative effects, respectively. The significance of standardized path coefficients is indicated by numbers on the arrows. (*, **, and *** denote significant differences at the p < 0.05, p < 0.01, and p < 0.001 levels, respectively). R2 represents the explained variance of individual variables. Mycorrhizal dominance is treated as categorical, representing stand types (i.e., Pinus sylvestris var. mongolica monoculture, Ulmus pumila monoculture, and Pinus × Ulmus mixed stands); soil nutrient (i.e., total carbon, total nitrogen, total phosphorus, organic carbon, available nitrogen, and available phosphorus); fungus (i.e., fungus observed richness and fungus Chao1 index); bacteria (i.e., bacteria observed richness and bacteria Chao1 index); and soil microbial biomass (i.e., soil microbial biomass carbon, soil microbial biomass nitrogen, and soil microbial biomass phosphorus). Figure 7 presents direct path coefficients only; total effects reported in the main text represent the sum of direct and indirect effects.

4 Discussion

4.1 Effects of AM and EcM monocultures on rhizosphere soil

The competition for nutrients between different mycorrhizal types regulates soil nutrient cycling (Seidel et al., 2024). EcM ecosystems typically rely on the mineralization of organic nutrients and fungus-mediated decomposition processes. This nutrient acquisition strategy satisfies plant nitrogen and phosphorus demands and inhibits the mineralization of organic carbon, thereby promoting SOC accumulation and maintaining elevated nutrient availability (Averill et al., 2016; Tunlid et al., 2022; Guan et al., 2025). By contrast, AM ecosystems primarily depend on the direct absorption of inorganic nutrients, exhibiting an increased demand for labile N and P. Consequently, soil nutrient dynamics in AM systems are more susceptible to environmental fluctuations than those in EcM systems (Faghihinia et al., 2024). Consistent with these mechanisms, our results showed that soils in AM monocultures had lower nutrient contents but significantly higher pH than those in EcM monocultures (p < 0.05). In EcM-dominated ecosystems, the competition for organic N and suppression of saprotrophic decomposers (the “Gadgil effect”) impede organic matter decomposition, enhancing the retention of soil C and organic nutrient pools (Gadgil and Gadgil, 1971; Orwin et al., 2011). Similarly, Averill et al. (2014) demonstrated that EcM-dominated ecosystems store more SOC than AM-dominated ones. The study area is located in a typical arid-to-semiarid ecological zone. The soil, which is alkaline and has low humidity, tends to maintain AM fungal communities that are resistant to high pH and have rapid metabolism. However, it is also more susceptible to microbial limitations, thereby restricting the effective supply capacity of soil nutrients (Barea et al., 2011). Enzyme activities and microbial community structures further revealed these contrasting mechanisms. Brzostek et al. (2015) showed that EcM fungi preferentially produce key enzymes involved in the degradation of complex organic matter, such as NAG, LAP, and ALP, to acquire organic nutrients. Similarly, our study found that NAG, LAP, and ALP activities; microbial biomass C, N, and P; and fungal diversity were significantly higher in EcM monocultures than in AM monocultures (p < 0.05). Conversely, AM monocultures exhibited higher BG activity, stronger microbial C limitation (vector length), N limitation (vector angle) (p < 0.05), and higher bacterial diversity than EcM monocultures. This pattern may stem from the higher pH and lower soil moisture observed in AM pure forests. In such alkaline and low-moisture environments, bacterial taxa with greater osmotic-stress tolerance and rapid metabolic capacity are more likely to be selected and maintain higher diversity. Consequently, bacterial diversity in AM pure forests is higher than that in EcM pure forests. Zheng et al. (2022) reported that AM tree species in European temperate forests are associated with increased microbial activity and bacterial diversity. The discrepancy in our results may be due to differences in climate, soil pH, and site conditions. The study area is located in a typical arid–semiarid ecosystem, characterized by generally low soil moisture, limited total inorganic nitrogen, and a relatively high soil C: N stoichiometry. These abiotic constraints inherently restrict microbial decomposition and nutrient regeneration rates, thereby exacerbating microbial demands for C and N and resulting in stronger microbial nutrient limitations. In contrast, the study site investigated by Zheng et al. is situated in a temperate European forest with higher soil fertility and a more neutral pH, where AM-associated tree species can sustain greater microbial biomass under more nutrient-rich conditions. Therefore, differences in microbial activity and nutrient limitation observed across regions—or even between forest types dominated by AM versus EcM trees—may partly arise from environmental filtering effects imposed by local abiotic stressors on extracellular enzyme investment patterns and microbial metabolic efficiency. This mechanism indicates that environmental constraints and mycorrhizal type jointly shape microbial nutrient acquisition strategies and patterns of resource limitation. Microbial network analysis further revealed that EcM monocultures had higher node number, edge number, average degree, average path length, clustering coefficient, and modularity than AM monocultures, indicating complex and highly connected microbial communities. This may be attributed to the extensive EcM hyphal networks that promote the decomposition of complex organic matter and nutrient release, thereby enhancing microbial biomass and network stability. By contrast, the simplified microbial networks in AM systems are likely driven by high soil pH and nutrient scarcity. Overall, the contrasting patterns observed in AM and EcM rhizospheres not only reflect the fundamental differences in their nutrient acquisition strategies but also highlight their far-reaching impacts on soil C, N, and P cycling, as well as on the stability of microbial interaction networks.

4.2 Effects of AM–EcM mixed stands on rhizosphere soil properties

Compared with monocultures, mixed forests usually exhibit enhanced environmental heterogeneity, which influences soil nutrient cycling and microbial processes at both the rhizosphere and stand scales (Kang et al., 2022; Jiang et al., 2017). The coexistence of multiple tree species in mixed stands creates diverse microhabitats and resource niches, thereby supporting a broad range of soil microorganisms, particularly fungi (Gould et al., 2016). This increased diversity often results in higher structural and functional complexity, yielding more stable and complementary ecological niches. At the stand scale, mixed forests exhibit greater litter diversity, with pronounced differences in stoichiometric traits, lignin content, and decomposition rates. Such heterogeneity in litter inputs can alter the composition of soil organic substrates, enhance microbial decomposition potential, and ultimately increase nutrient regeneration efficiency at the ecosystem level (Hättenschwiler et al., 2005; Prescott and Grayston, 2013). At the rhizosphere scale, AM and EcM trees display markedly complementary nutrient acquisition strategies: EcM fungi degrade organic substrates and promote the mineralization of nitrogen and phosphorus, whereas AM fungi efficiently utilize inorganic N and P. This complementary mechanism diversifies the nutrient sources available in mixed forests and enhances soil organic matter accumulation. Previous studies have shown that the coexistence of AM and EcM tree species can enhance soil resource-use efficiency (Deng et al., 2023; Singavarapu et al., 2022). Our study found that, except for SOC, the nutrient contents in the rhizosphere soil of P. sylvestris var. mongolica in mixed stands were lower than those in its pure stands (P < 0.05), while soil pH was higher than in the pure stands. This pattern may be attributed to enhanced nutrient transfer into plant and microbial biomass, as well as the release of alkaline cations such as Ca and Mg during litter decomposition, which can increase soil pH (Vesterdal et al., 2013). In mixed forests, the nutrient contents in the rhizosphere soil of U. pumila were all significantly higher than those in its pure stands, whereas soil pH was significantly reduced. This may result from AM tree species utilizing EcM-mediated organic nutrient mineralization, thereby enhancing P and N uptake efficiency. In addition, AM trees can lower soil pH through increased cation uptake, proton release, and organic acid secretion (Clarholm et al., 2015). The complementary stoichiometric traits of mixed litter and differential decomposition processes may accelerate SOC accumulation and nutrient turnover, particularly in AM-dominated ecosystems where nutrient limitation is pronounced. Luo et al. (2023) also demonstrated that, although fungal communities dominated by extracellular enzyme-producing groups have been shown to decompose organic matter slowly in any soil environment, the presence of herbaceous AM plants in the understory of such enzyme-dominated forests can accelerate the decomposition of soil organic matter. If AM tree species similarly promote decomposition, nutrient cycling in mixed forests may exceed that in stands dominated by a single mycorrhizal type. Furthermore, functional differences between AM and ECM may lead to complementary effects, enhancing productivity and alleviating microbial resource limitations in forests that adopt a mixed mycorrhizal strategy. Our study showed that vector length in mixed stands significantly decreased but vector angle significantly increased than in monocultures (p < 0.05). This indicates that mixed mycorrhizal strategies not only accelerate nutrient cycling and optimize the stoichiometric structure of organic matter decomposition and nutrient release but also move microbial nutrient acquisition close to the theoretical balance of C: N: P (1:1:1). As a result, microbial nutrient limitations are alleviated, promoting coordinated and balanced microbial resource use (Cheeke et al., 2017; Wu et al., 2023). This pattern may stem from the complementary functions of AM and EcM fungi jointly regulating organic matter decomposition and nutrient regeneration. EcM fungi enhance the mineralization of organic nitrogen and phosphorus, whereas AM fungi efficiently absorb inorganic P, thereby increasing the overall availability of soil N and P. At the same time, mixed forests increase SOC and microbial biomass, improving the structure of substrate supply and allowing microbes to acquire C, N, and P in a more balanced manner rather than investing disproportionately in a single class of extracellular enzymes. Previous studies have also noted that mixed-mycorrhizal systems can enhance enzymatic diversity and complementarity, enabling microbes to maintain stoichiometric balance, minimize metabolic costs, and improve resource-use efficiency (Camenzind et al., 2018; Moorhead and Sinsabaugh, 2006). Several studies (Banerjee et al., 2018; Wagg et al., 2019) have reported that mixed mycorrhizal strategies can increase fungal diversity and microbial biomass and enhance interactions and metabolic complementarity among functional groups, thereby forming highly interconnected microbial networks. Consistent with this, we found that MBC, MBN, MBP, and fungal diversity were significantly higher in mixed stands than in monocultures (p < 0.05), and microbial networks exhibited high complexity. By contrast, bacterial diversity was slightly reduced (p > 0.05), possibly due to shifts in organic matter inputs and strengthened fungal functions, which may partly replace bacterial ecological niches (Fierer et al., 2009).

4.3 Relationships among microbial nutrient limitation, mycorrhizal dominance, soil chemistry properties, and microbial biomass

Mixed forests, through a mixed-mycorrhizal strategy, promote complementary nutrient use between AM and EcM trees, thereby improving soil nutrient supply and microbial resource acquisition (Guo and Gong, 2024; Chen et al., 2020). By integrating the efficient inorganic nutrient uptake of AM fungi with the organic matter decomposition capacity of EcM fungi, mixed forests optimize soil nutrient availability and restructure microbial communities. This restructuring enhances fungal diversity and microbial network complexity, which in turn alleviates microbial C and N limitations. For example, Camenzind et al. (2018) found that mixed forests in the tropics increase microbial biomass and alter extracellular enzyme allocation patterns, thereby reducing microbial nutrient limitation. Our PLS-PM results further revealed that mycorrhizal dominance exerts a significant direct regulatory effect on microbial nutrient limitation. Mycorrhizal dominance had a highly significant positive direct effect on vector angle (path coefficient = 0.70, p < 0.001), suggesting that the dominant mycorrhizal type can directly alter the balance of microbial investment in N-acquiring extracellular enzymes. Microbial biomass and mycorrhizal dominance exerted strong negative effects on vector length (path coefficients = −0.71 and −0.33, respectively, p < 0.01), indicating that the dominant mycorrhizal type can directly modify the balance of microbial investment in N-acquiring extracellular enzymes. These results suggest that changes in microbial biomass and mycorrhizal dominance alter the distribution of extracellular enzymes among C, N, and P, serving as the principal factors in alleviating microbial resource limitations. Mycorrhizal dominance not only exerts direct effects but also indirectly regulates microbial resource balance by altering soil nutrient status and microbial biomass. The dominance of mycorrhizae exerted a significant negative effect on soil pH (path coefficient = −0.41, p < 0.01). The total effects of PLS-PM revealed that pH exerted an overall positive effect on vector length (path coefficient = 0.31) and a negative effect on vector angle (path coefficient = −0.02) by influencing nutrient content (path coefficient = −0.51, p < 0.001) and microbial biomass (path coefficient = −0.11). This result was aligned with the findings reported by Yuan et al. (2023), who found that mixed forests enhance nutrient availability by regulating litter decomposition and rhizosphere processes that modify soil acidity. In addition, mycorrhizal dominance had a strong positive effect on soil nutrients (0.61, p < 0.001), and nutrient enrichment indirectly promoted fungal diversity (0.83, p < 0.001) and microbial biomass (0.51, p < 0.001). These findings indicate that nutrient availability directly enhances fungal involvement in organic matter decomposition and nutrient uptake, as well as overall microbial metabolic potential. Notably, fungal diversity had a significant negative effect on vector angle (−0.68, p < 0.05), suggesting that increased fungal proportions shift extracellular enzyme allocation toward P acquisition, thereby increasing N availability; this response was likely associated with the ability of EcM fungi to access organic P (Sannakajsa et al., 2014). Furthermore, microbial biomass exhibited a negative effect on vector length but a positive effect on vector angle. The GAM analysis further revealed the environmental response patterns underlying microbial resource limitation. The results showed that pH exhibited linear effects on both vector length and vector angle (edf = 1), indicating that soil pH acts as a stable environmental gradient that drives microbial extracellular enzyme investment strategies. Nutrient PC1 achieved the highest explanatory power for vector length and vector angle (R2 > 0.90), demonstrating that the supply structure of soil C, N, and P fundamentally determines the direction of microbial resource-limitation responses. Meanwhile, the negative effect of microbial biomass on vector length suggests that higher metabolic capacity can alleviate C limitation, although such alleviation relies on external resource supply rather than increased microbial demand. Notably, GAM captures the overall response along environmental gradients, emphasizing a “supply-driven” mechanism—meaning that nutrient supply structure sets the background pattern of microbial limitation. In contrast, PLS-PM, after controlling for environmental factors, reveals a “demand-driven” mechanism—where microbial biomass and mycorrhizal function actively mitigate resource limitation through adjustments in extracellular enzyme allocation. Therefore, the combined evidence from this study indicates that microbial nutrient limitation is shaped not by supply-side or demand-side factors alone, but by a dual regulatory mechanism involving both nutrient supply conditions and microbial metabolic demand. These findings highlight that mixed forests, compared with monocultures, alleviate microbial resource limitations through multiple regulatory pathways, emphasizing the critical role of mixed-mycorrhizal strategies in maintaining ecosystem function and nutrient cycling balance. It should be noted that this study is based on rhizosphere soil, whose microenvironment is strongly influenced by root exudates, mycorrhizal symbiosis, and the interaction between rhizosphere microorganisms. Therefore, the microbial resource limitation patterns and enzyme activity regulation mechanisms revealed in this study mainly reflect the characteristics of the rhizosphere. In large soil masses or deeper soil layers, due to the significantly weakened influence of roots, the microbial community structure, enzyme investment strategies, and their responses to C, N, and P limitations may present different patterns. Therefore, the explanations of the related mechanisms in this study should be limited to the rhizosphere scale.

5 Conclusion

Our study highlights the regulatory mechanisms governing soil microbial community structure and nutrient limitations under different tree mycorrhizal dominances. We found that AM pure forests, characterized by high BG activity, maintained higher bacterial diversity but were accompanied with intensified microbial C and N limitations. By contrast, EcM pure forests demonstrated elevated NAG + LAP and ALP activities. Mixed forests of U. pumila and P. sylvestris var. mongolica demonstrated a synergistic effect of mixed mycorrhizal associations, which significantly enhanced soil SOC and microbial biomass, regulated extracellular enzyme activities, alleviated microbial nutrient limitations, and promoted the formation of complex and stable microbial interaction networks. PLS-PM further revealed that mycorrhizal dominance and microbial biomass were the primary factors alleviating microbial nutrient limitations. Overall, tree mycorrhizal dominance influences microbial nutrient limitations through variations in extracellular enzyme activities and microbial community assembly. In particular, AM–EcM mixed forests can achieve resource complementarity via mixed mycorrhizal strategies, thereby alleviating microbial nutrient limitations and enhancing the stability of forest soil ecosystems. These findings provide important theoretical support for the restoration of degraded forests, the optimization of mixed-forest configurations, and ecological rehabilitation in sandy land regions.

Data availability statement

Raw data have been deposited to National Center for Biotechnology Information (NCBI) under the BioProject number PRJNA1379231.

Author contributions

YX: Software, Visualization, Writing – original draft, Formal analysis, Conceptualization, Methodology, Data curation, Investigation. LH: Data curation, Methodology, Supervision, Validation, Writing – review & editing. YY: Funding acquisition, Writing – review & editing, Project administration, Supervision. RG: Resources, Writing – review & editing, Supervision, Methodology, Validation. LZ: Writing – review & editing, Investigation, Data curation, Supervision. KZ: Investigation, Visualization, Methodology, Writing – review & editing, Software. ZZ: Investigation, Validation, Supervision, Software, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work is supported by National Natural Science Foundation of China (NSFC) Project (32560307), Central Guidance Local Scientific and Technological Development Fund Project (2025ZY0038), Inner Mongolia Autonomous Region Science and Technology Plan Project (2022YFDZ0030), the “Unveiling and Leading” Project of Inner Mongolia (2024JBGS0021-3-1), Basic Research Operating Funds Project of Universities under the Direct Administration of the Inner Mongolia Autonomous Region (BR221013).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

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

Footnotes

References

Averill, C., Hawkes, C. V., and Bardgett, R. (2016). Ectomycorrhizal fungi slow soil carbon cycling. Ecol. Lett. 19, 937–947. doi: 10.1111/ele.12631,

PubMed Abstract | Crossref Full Text | Google Scholar

Averill, C., Turner, B. L., and Finzi, A. C. (2014). Mycorrhiza-mediated competition between plants and decomposers drives soil carbon storage. Nature 505, 543–545. doi: 10.1038/nature12901,

PubMed Abstract | Crossref Full Text | Google Scholar

Bahram, M., Netherway, T., Hildebrand, F., Pritsch, K., Drenkhan, R., Loit, K., et al. (2020). Plant nutrient acquisition strategies drive topsoil microbiome structure and function. New Phytol. 227, 1189–1199. doi: 10.1111/nph.16598,

PubMed Abstract | Crossref Full Text | Google Scholar

Banerjee, S., Schlaeppi, K., and van der Heijden, M. G. A. (2018). Keystone taxa as drivers of microbiome structure and functioning. Nat. Rev. Microbiol. 16, 567–576. doi: 10.1038/s41579-018-0024-1,

PubMed Abstract | Crossref Full Text | Google Scholar

Barea, J. M., Palenzuela, J., Cornejo, P., Sánchez-Castro, I., Navarro-Fernández, C., López-García, A., et al. (2011). Ecological and functional roles of mycorrhizas in semi-arid ecosystems of Southeast Spain. J. Arid Environ. 75, 1292–1301. doi: 10.1016/j.jaridenv.2011.06.001

Crossref Full Text | Google Scholar

Brookes, P. C., Landman, A., Pruden, G., and Jenkinson, D. S. (1985). Chloroform fumigation and the release of soil nitrogen: a rapid direct extraction method to measure microbial biomass nitrogen in soil. Soil Biol. Biochem. 17, 837–842. doi: 10.1016/0038-0717(85)90144-0

Crossref Full Text | Google Scholar

Brundrett, M. C., and Tedersoo, L. (2018). Evolutionary history of mycorrhizal symbioses and global host plant diversity. New Phytol. 220, 1108–1115. doi: 10.1111/nph.14976,

PubMed Abstract | Crossref Full Text | Google Scholar

Brzostek, E. R., Dragoni, D., Brown, Z. A., and Phillips, R. P. (2015). Mycorrhizal type determines the magnitude and direction of root-induced changes in decomposition in a temperate forest. New Phytol. 206, 1274–1282. doi: 10.1111/nph.13303,

PubMed Abstract | Crossref Full Text | Google Scholar

Camenzind, T., Hättenschwiler, S., Treseder, K. K., Lehmann, A., and Rillig, M. C. (2018). Nutrient limitation of soil microbial processes in tropical forests. Ecol. Monogr. 88, 4–21. doi: 10.1002/ecm.1279

Crossref Full Text | Google Scholar

Cheeke, T. E., Phillips, R. P., Brzostek, E. R., Rosling, A., Bever, J. D., and Fransson, P. (2017). Dominant mycorrhizal association of trees alters carbon and nutrient cycling by selecting for microbial groups with distinct enzyme function. New Phytol. 214, 432–442. doi: 10.1111/nph.14343,

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, W., Koide, R. T., Adams, T. S., DeForest, J. L., Cheng, L., and Eissenstat, D. M. (2020). Root morphology and mycorrhizal symbioses together shape nutrient foraging strategies of temperate trees. Proc. Natl. Acad. Sci. 113, 8741–8746. doi: 10.1073/pnas.1601006113,

PubMed Abstract | Crossref Full Text | Google Scholar

Clarholm, M., Skyllberg, U., and Rosling, A. (2015). Organic acid induced release of nutrients from metal-stabilized soil organic matter – the unbutton model. Soil Biol. Biochem. 84, 168–176. doi: 10.1016/j.soilbio.2015.02.019

Crossref Full Text | Google Scholar

Cleveland, C. C., and Liptzin, D. (2007). C:N:P stoichiometry in soil: is there a “Redfield ratio” for the microbial biomass? Biogeochemistry 85, 235–252. doi: 10.1007/s10533-007-9132-0

Crossref Full Text | Google Scholar

Cornfield, A. H. (1960). Ammonia released on treating soils with N sodium hydroxide as a possible means of predicting the nitrogen-supplying power of soils. Nature 187, 260–261. doi: 10.1038/187260a0

Crossref Full Text | Google Scholar

Craig, M. E., Turner, B. L., Liang, C., Clay, K., Johnson, D. J., and Phillips, R. P. (2018). Tree mycorrhizal type predicts within-site variability in the storage and distribution of soil organic matter. Glob. Chang. Biol. 24, 3317–3330. doi: 10.1111/gcb.14132,

PubMed Abstract | Crossref Full Text | Google Scholar

Cui, Y., Bing, H., Fang, L., Jiang, M., Shen, G., Yu, J., et al. (2021). Extracellular enzyme stoichiometry reveals the carbon and phosphorus limitations of microbial metabolisms in the rhizosphere and bulk soils in alpine ecosystems. Plant Soil 458, 7–20. doi: 10.1007/s11104-019-04159-x

Crossref Full Text | Google Scholar

Deng, M., Hu, S., Guo, L., Jiang, L., Huang, Y., Schmid, B., et al. (2023). Tree mycorrhizal association types control biodiversity-productivity relationship in a subtropical forest. Sci. Adv. 9:eadd4468. doi: 10.1126/sciadv.add4468,

PubMed Abstract | Crossref Full Text | Google Scholar

Deng, C., Zhang, S., Lu, Y., Froese, R. E., Xu, X., Zeng, J., et al. (2020). Thinning effects on forest evolution in Masson pine (Pinus massoniana lamb.) conversion from pure plantations into mixed forests. For. Ecol. Manag. 477:118503. doi: 10.1016/j.foreco.2020.118503

Crossref Full Text | Google Scholar

Ding, Y., Gao, X., Shu, D., et al. (2025). Enhancing soil health and nutrient cycling through soil amendments: improving the synergy of bacteria and fungi. Sci. Total Environ. 923:171332. doi: 10.1016/j.scitotenv.2024.171332

Crossref Full Text | Google Scholar

Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604,

PubMed Abstract | Crossref Full Text | Google Scholar

Faghihinia, M., Halverson, L. J., Hršelová, H., Bukovská, P., Rozmoš, M., Kotianová, M., et al. (2024). Nutrient-dependent cross-kingdom interactions in the hyphosphere of an arbuscular mycorrhizal fungus. Front. Microbiol. 14:1284648. doi: 10.3389/fmicb.2023.1284648,

PubMed Abstract | Crossref Full Text | Google Scholar

Feng, J., Wei, K., Chen, Z., Lu, X., Tian, J., Wang, C., et al. (2019). Coupling and decoupling of soil carbon and nutrient cycles across an aridity gradient in the drylands of northern China: evidence from ecoenzymatic stoichiometry. Glob. Biogeochem. Cycles 33, 559–569. doi: 10.1029/2018GB006112

Crossref Full Text | Google Scholar

Fierer, N., Strickland, M. S., Liptzin, D., Bradford, M. A., and Cleveland, C. C. (2009). Global patterns in belowground communities. Ecol. Lett. 12, 1238–1249. doi: 10.1111/j.1461-0248.2009.01360.x,

PubMed Abstract | Crossref Full Text | Google Scholar

Gadgil, R. L., and Gadgil, P. D. (1971). Mycorrhiza and litter decomposition. Nature 233:133. doi: 10.1038/233133a0,

PubMed Abstract | Crossref Full Text | Google Scholar

Gould, I. J., Quinton, J. N., Weigelt, A., De Deyn, G. B., and Bardgett, R. D. (2016). Plant diversity and root traits benefit physical properties key to soil function in grasslands. Ecol. Lett. 19, 1140–1149. doi: 10.1111/ele.12652,

PubMed Abstract | Crossref Full Text | Google Scholar

Guan, X., Jiang, J., Classen, A. T., Ullah, S., and Wang, G. (2025). Disentangling the contribution of mycorrhizal fungi to soil organic carbon storage. Soil Biol. Biochem. 209:109900. doi: 10.1016/j.soilbio.2025.109900

Crossref Full Text | Google Scholar

Guo, Q., and Gong, L. (2024). Compared with pure forest, mixed forest alters microbial diversity and increases the complexity of interdomain networks in arid areas. Microbiol. Spectr. 12:e0264223. doi: 10.1128/spectrum.02642-23,

PubMed Abstract | Crossref Full Text | Google Scholar

Hättenschwiler, S., Tiunov, A. V., and Scheu, S. (2005). Biodiversity and litter decomposition in terrestrial ecosystems. Annu. Rev. Ecol. Evol. Syst. 36, 191–218. doi: 10.1146/annurev.ecolsys.36.112904.151932

Crossref Full Text | Google Scholar

Heděnec, P., Nilsson, L. O., Zheng, H., Sihi, D., Vásquez-Gassibe, P., Gómez-Brandón, M., et al. (2020). Mycorrhizal association of common European tree species shapes biomass and metabolic activity of bacterial and fungal communities in soil. Soil Biol. Biochem. doi: 10.1016/j.soilbio.2020.107933

Crossref Full Text | Google Scholar

Jiang, S., Liu, Y., Luo, J., Qin, M., Johnson, N. C., Öpik, M., et al. (2018). Dynamics of arbuscular mycorrhizal fungal community structure and functioning along a nitrogen enrichment gradient in an alpine meadow ecosystem. New Phytol. 220, 1222–1235. doi: 10.1111/nph.15112,

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, J., Moore, J. A. M., Priyadarshi, A., and Classen, A. T. (2017). Plant-mycorrhizal interactions mediate plant community coexistence by altering resource demand. Ecology 98, 187–197. doi: 10.1002/ecy.1630,

PubMed Abstract | Crossref Full Text | Google Scholar

Jo, I., Fei, S., Oswalt, C. M., Domke, G. M., and Phillips, R. P. (2019). Shifts in dominant tree mycorrhizal associations in response to anthropogenic impacts. Sci. Adv. 5:eaav6358. doi: 10.1126/sciadv.aav6358,

PubMed Abstract | Crossref Full Text | Google Scholar

Kang, P., Pan, Y., Yang, P., Hu, J., Zhao, T., Zhang, Y., et al. (2022). A comparison of microbial composition under three tree ecosystems using the stochastic process and network complexity approaches. Front. Microbiol. 13:1018077. doi: 10.3389/fmicb.2022.1018077,

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, T., Hao, L., and Bai, S. (2024). Ecoenzymatic stoichiometry and microbial nutrient limitation in response to shrub planting patterns and arbuscular mycorrhizal fungal inoculation. Appl. Soil Ecol. 200:105449. doi: 10.1016/j.apsoil.2024.105449

Crossref Full Text | Google Scholar

Liu, S., Xu, G., Chen, H., Zhang, M., Cao, X., Chen, M., et al. (2023). Contrasting responses of soil microbial biomass and extracellular enzyme activity along an elevation gradient on the eastern Qinghai-Tibetan plateau. Front. Microbiol. 14:974316. doi: 10.3389/fmicb.2023.974316,

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, S., Phillips, R. P., Jo, I., Fei, S., Liang, J., Schmid, B., et al. (2023). Higher productivity in forests with mixed mycorrhizal strategies. Nat. Commun. 14:1377. doi: 10.1038/s41467-023-36888-0,

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, B., Wang, Y., Ye, S., Liu, S., Stirling, E., Gilbert, J. A., et al. (2020). Earth microbial co-occurrence network reveals interconnection pattern across microbiomes. Microbiome 8:82. doi: 10.1186/s40168-020-00857-2,

PubMed Abstract | Crossref Full Text | Google Scholar

Mebius, L. J. (1960). A rapid method for the determination of organic carbon in soil. Anal. Chim. Acta 22, 120–124. doi: 10.1016/S0003-2670(00)88254-9

Crossref Full Text | Google Scholar

Moorhead, D. L., Rinkes, Z. L., Sinsabaugh, R. L., and Weintraub, M. N. (2013). Dynamic relationships between microbial biomass, respiration, inorganic nutrients and enzyme activities: informing enzyme-based decomposition models. Front. Microbiol. 4:223. doi: 10.3389/fmicb.2013.00223,

PubMed Abstract | Crossref Full Text | Google Scholar

Moorhead, D. L., and Sinsabaugh, R. L. (2006). A theoretical model of litter decay and microbial interaction. Ecol. Monogr. 76, 151–174. doi: 10.1890/0012-9615(2006)076[0151:atmold]2.0.co;2

Crossref Full Text | Google Scholar

Öpik, M., Vanatoa, A., Vanatoa, E., Moora, M., Davison, J., Kalwij, J. M., et al. (2010). The online database MaarjAM reveals global and ecosystemic distribution patterns in arbuscular mycorrhizal fungi (Glomeromycota). New Phytol. 188, 223–241. doi: 10.2307/40856395

Crossref Full Text | Google Scholar

Orwin, K. H., Kirschbaum, M. U. F., St John, M. G., and Dickie, I. A. (2011). Organic nutrient uptake by mycorrhizal fungi enhances ecosystem carbon storage: a model-based assessment. Ecol. Lett. 14, 493–502. doi: 10.1111/j.1461-0248.2011.01611.x,

PubMed Abstract | Crossref Full Text | Google Scholar

Page, A. L., Miller, R. H., and Keeney, D. R. (1982). Methods of soil analysis. Part 2. Chemical and microbiological properties. Madison, WI: American Society of Agronomy, 403–430.

Google Scholar

Phillips, R. P., Brzostek, E., and Midgley, M. G. (2013). The mycorrhizal-associated nutrient economy: a new framework for predicting carbon–nutrient couplings in temperate forests. New Phytol. 199, 41–51. doi: 10.1111/nph.12221,

PubMed Abstract | Crossref Full Text | Google Scholar

Polyakova, O., and Billor, N. (2007). Impact of deciduous tree species on litterfall quality, decomposition rates and nutrient circulation in pine stands. For. Ecol. Manag. 253, 11–18. doi: 10.1016/j.foreco.2007.06.049

Crossref Full Text | Google Scholar

Prescott, C. E., and Grayston, S. J. (2013). Tree species influence on microbial communities in litter and soil: current knowledge and research needs. For. Ecol. Manag. 309, 19–27. doi: 10.1016/j.foreco.2013.02.034

Crossref Full Text | Google Scholar

Saiya-Cork, K. R., Sinsabaugh, R. L., and Zak, D. R. (2002). The effects of long term nitrogen deposition on extracellular enzyme activity in an Acer saccharum forest soil. Soil Biol. Biochem. 34, 1309–1315. doi: 10.1016/S0038-0717(02)00074-3

Crossref Full Text | Google Scholar

Sannakajsa, M., Velmala, T., Heinonsalo, J., Taylor, A. F. S., and Pennanen, T. (2014). Profiling functions of ectomycorrhizal diversity and root structuring in seedlings of Norway spruce (Picea abies) with fast- and slow-growing phenotypes. New Phytol. 201, 610–622. doi: 10.2307/newphytologist.201.2.610

Crossref Full Text | Google Scholar

Seidel, F., Castaño, C., Alday, J. G., Lopez, M. L., and Bonet, J. A. (2024). Contrasting fungal functional groups influence nutrient cycling across four Japanese cool-temperate forest soils. Appl. Soil Ecol. 198:105360. doi: 10.1016/j.apsoil.2024.105360

Crossref Full Text | Google Scholar

Shu, X., Hu, Y., Liu, W., Wang, H., Chen, L., and Zhou, H. (2023). Linking between soil properties, bacterial communities, enzyme activities, and soil organic carbon mineralization under ecological restoration in an alpine degraded grassland. Front. Microbiol. 14:1131836. doi: 10.3389/fmicb.2023.1131836,

PubMed Abstract | Crossref Full Text | Google Scholar

Singavarapu, B., Beugnon, R., Bruelheide, H., Cesarz, S., Du, J., Eisenhauer, N., et al. (2022). Tree mycorrhizal type and tree diversity shape the forest soil microbiota. Environ. Microbiol. 24, 4236–4255. doi: 10.1111/1462-2920.15690,

PubMed Abstract | Crossref Full Text | Google Scholar

Sinsabaugh, R. L., Hill, B. H., and Follstad Shah, J. J. (2009). Ecoenzymatic stoichiometry of microbial organic nutrient acquisition in soil and sediment. Nature 462, 795–798. doi: 10.1038/nature08632,

PubMed Abstract | Crossref Full Text | Google Scholar

Tang, H., Li, Q., Han, Y., Dong, Z., Zhou, Y., Jiang, Z., et al. (2024a). Interplay of soil characteristics and arbuscular mycorrhizal fungi diversity in alpine wetland restoration and carbon stabilization. Front. Microbiol. 15:1376418. doi: 10.3389/fmicb.2024.1376418,

PubMed Abstract | Crossref Full Text | Google Scholar

Tang, X., Yang, J., Lin, D., Lin, H., Xiao, X., Chen, S., et al. (2024b). Community assembly of ectomycorrhizal fungal communities in pure and mixed Pinus massoniana forests. J. Environ. Manag. 362:121312. doi: 10.1016/j.jenvman.2024.121312,

PubMed Abstract | Crossref Full Text | Google Scholar

Tedersoo, L., and Bahram, M. (2019). Mycorrhizal types differ in ecophysiology and alter plant nutrition and soil processes. Biol. Rev. 94, 1857–1880. doi: 10.1111/brv.12538,

PubMed Abstract | Crossref Full Text | Google Scholar

Tighe, M., Edwards, M. M., and Cluley, G. (2018). Colorimetrically determining total antimony in contaminated waters and screening for antimony speciation. J. Hydrol. 563, 84–91. doi: 10.1016/j.jhydrol.2018.05.056

Crossref Full Text | Google Scholar

Tunlid, A., Floudas, D., de Op Beeck, M., Wang, T., and Persson, P. (2022). Decomposition of soil organic matter by ectomycorrhizal fungi: mechanisms and consequences for organic nitrogen uptake and soil carbon stabilization. Front. For. Glob. Change 5:934409. doi: 10.3389/ffgc.2022.934409

Crossref Full Text | Google Scholar

Vesterdal, L., Clarke, N., Sigurdsson, B. D., and Gundersen, P. (2013). Do tree species influence soil carbon stocks in temperate and boreal forests? For. Ecol. Manag. 309, 4–18. doi: 10.1016/j.foreco.2013.01.017

Crossref Full Text | Google Scholar

Wagg, C., Schlaeppi, K., Banerjee, S., Kuramae, E. E., and van der Heijden, M. G. A. (2019). Fungal–bacterial diversity and microbiome complexity predict ecosystem functioning. Nat. Commun. 10:4841. doi: 10.1038/s41467-019-12798-y,

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, X., Li, Y., Wang, L., Duan, Y., Yao, B., Chen, Y., et al. (2023). Soil extracellular enzyme stoichiometry reflects microbial metabolic limitations in different desert types of northwestern China. Sci. Total Environ. 874:162504. doi: 10.1016/j.scitotenv.2023.162504,

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, J., Joergensen, R. G., Pommerening, B., Chaussod, R., and Brookes, P. C. (1990). Measurement of soil microbial biomass C by fumigation-extraction—an automated procedure. Soil Biol. Biochem. 22, 1167–1169. doi: 10.1016/0038-0717(90)90046-3

Crossref Full Text | Google Scholar

Wu, S.-W., Shi, Z.-Y., Huang, M., Yang, S., Yang, W.-Y., and Li, Y.-J. (2023). Influence of mycorrhiza on C:N:P stoichiometry in senesced leaves. J. Fungi 9:588. doi: 10.3390/jof9050588,

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, Z., Wang, Y., Sun, D., Li, C., Liu, L., Zhang, Q., et al. (2022). Soil nutrients and nutrient ratios influence the ratios of soil microbial biomass and metabolic nutrient limitations in mountain peatlands. Catena 219:106528. doi: 10.1016/j.catena.2022.106528.

Crossref Full Text | Google Scholar

Yu, Q., Ma, S., Ni, X., Jiang, L., Zhou, Z., Zhu, J., et al. (2023). A test of the mycorrhizal-associated nutrient economy framework in two types of tropical rainforests under nutrient enrichments. For. Ecosyst. 10:100083. doi: 10.1016/j.fecs.2022.100083.

Crossref Full Text | Google Scholar

Yuan, C., Wu, F., Peng, Y., Wu, Q., Zhu, G., Zhao, Z., et al. (2023). Pain or gain: the dual role of afforestation effects on soil pH at the global scale. Sci. Total Environ. 493:163609. doi: 10.1007/s11104-023-06254-6

Crossref Full Text | Google Scholar

Zeng, H. D., Du, Z. X., Yang, Y. S., Li, X. B., Zhang, Y. C., and Yang, Z. F. (2010). Effects of land cover change on soil organic carbon and light fraction organic carbon at river banks of Fuzhou urban area. Yingyong Shengtai Xuebao 21, 701–706. doi: 10.13287/j.1001-9332.2010.0092

Crossref Full Text | Google Scholar

Zhang, Y., Li, X., Zhang, D., Qin, Y., Zhou, Y., Song, S., et al. (2020). Characteristics of fungal community structure during the decomposition of mixed foliage litter from Pinus massoniana and broadleaved tree species in southwestern China. J. Plant Ecol. 13, 574–588. doi: 10.1093/jpe/rtz057

Crossref Full Text | Google Scholar

Zheng, H., Heděnec, P., Rousk, J., Schmidt, I. K., Peng, Y., and Vesterdal, L. (2022). Effects of common European tree species on soil microbial resource limitation, microbial communities and soil carbon. Soil Biol. Biochem. 172:108754. doi: 10.1016/j.soilbio.2022.108754

Crossref Full Text | Google Scholar

Zheng, H., Phillips, R. P., Rousk, J., Yue, K., Schmidt, I. K., Peng, Y., et al. (2023). Imprint of tree species mycorrhizal association on microbial-mediated enzyme activity and stoichiometry. Funct. Ecol. 37, 1366–1376. doi: 10.1111/1365-2435.14311,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, L., Liu, S., Shen, H., Zhao, M., Xu, L., Xing, A., et al. (2020). Soil extracellular enzyme activity and stoichiometry in China's forests. Funct. Ecol. 34, 1461–1471. doi: 10.1111/1365-2435.13555

Crossref Full Text | Google Scholar

Keywords: arbuscular mycorrhizal fungi, ectomycorrhizal fungi, mixed-mycorrhizal strategy, microbial nutrient limitation, soil extracellular enzyme activity, microbial community structure

Citation: Xu Y, Hao L, Yue Y, Gao R, Zhang L, Zhao K and Zhao Z (2026) Effects of tree mycorrhizal dominance on soil microbial community structure and microbial nutrient limitation. Front. Microbiol. 16:1698121. doi: 10.3389/fmicb.2025.1698121

Received: 04 September 2025; Revised: 01 December 2025; Accepted: 08 December 2025;
Published: 08 January 2026.

Edited by:

Juntao Wang, Western Sydney University, Australia

Reviewed by:

Tong Li, Zhengzhou University, China
Dipti Rai, North Carolina Agricultural and Technical State University, United States

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

*Correspondence: Longfei Hao, aGFvbG9uZ2ZlaV8wMEAxMjYuY29t; Yongjie Yue, d29sb25neXVlMTIzQDE2My5jb20=

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