Abstract
The ecological significance of rhizosphere microbiomes in forest ecosystems is increasingly recognized. This study provides comparative analysis of microbial communities in Pinus massoniana-Quercus acutissima mixed forests versus monoculture systems. Mixed stands exhibited superior rhizosphere nutrient conditions and supported more diverse microbial populations, particularly with respect to Proteobacteria, Actinobacteria, and Basidiomycota. Principal component analysis revealed clear separation between rhizosphere and non-rhizosphere soil communities, as well as distinct clustering patterns between mixed and pure forest types. The functional analysis revealed conserved metabolic pathways across forest stands, with bacterial metabolic processes and fungal saprophytic functions representing dominant community roles. Network topology analysis demonstrated enhanced connectivity in mixed forest rhizosphere systems, featuring Proteobacteria, Acidobacteria, and Actinobacteria as bacterial network hubs, while Basidiomycota emerged as central fungal network components. Environmental drivers exhibited differential influences, with bacterial assemblages responding primarily to soil pH, organic carbon content, and phosphorus availability, whereas fungal communities showed stronger associations with organic carbon and potassium levels. These findings collectively demonstrate that mixed-species plantations foster robust microbial networks through microenvironmental regulation, offering valuable insights for sustainable forest management practices.
1 Introduction
The rhizosphere, as a critical plant–soil interface, represents the soil microenvironment predominantly shaped by root exudates and microbial activity (Smith, 1976; Han et al., 2023). Characterized by intensive biogeochemical cycling and energy transformation, this dynamic zone serves as a key regulator of nutrient fluxes in forest ecosystems (Kuzyakov and Razavi, 2019). Recent studies have demonstrated that variations in forest vegetation composition—particularly between pure (monoculture) and mixed-species stands—play a critical role in shaping soil microbial community structure and diversity (Gilliam et al., 2023; Gilliam, 2025). Although numerous investigations confirm that mixed-species plantations modify soil microbial resource distribution, existing research predominantly treat soil as a homogeneous entity while neglecting the distinct responses of microbial metabolic activity in rhizospheric versus non-rhizospheric soils to tree species mixing. Substantial evidence indicates that rhizospheric soils exhibit distinct characteristics compared to non-rhizospheric soils, particularly in C/N/P stoichiometry, microbial assemblage structure, and enzymatic activities governing nutrient transformation processes (Gan et al., 2021; Ling et al., 2022; Ma et al., 2023). Understanding the mechanisms underlying mixed-species effects on soil stoichiometric traits (rhizosphere vs. non-rhizosphere) and microbial functional metabolism is essential for deciphering tree diversity-driven biogeochemical processes, particularly in microbial niche differentiation and carbon-nutrient coupling dynamics.
The introduction of multi-species tree mixtures generally enhances soil microbial biomass and diversity while modifying community composition (Pereira et al., 2019; Li et al., 2021). Notably, mixed plantations of Pinus massoniana and Schima superba demonstrate elevated soil nitrogen mineralization rates, coupled with increased rhizosphere microbial biomass and a higher fungal-to-bacterial ratio (Wu et al., 2022). Comparative studies reveal that mixed poplar (Populus spp.) and black locust (Robinia pseudoacacia) stands exhibit greater rhizosphere microbial abundance and biochemical activity than monocultures (Hu et al., 2002). Similarly, soil bacterial diversity in monoculture Cunninghamia lanceolata plantations is significantly lower than in mixed conifer-broadleaf forests, though fungal diversity displays an inverse trend (Liu et al., 2010). Furthermore, soil fungal communities in pure Pinus stands are dominated by saprophytic taxa, whereas symbiotic fungi become predominant when pines are interplanted with broadleaf species (Wu et al., 2019; Li et al., 2021). These findings suggest that mixed forest management fosters cross-kingdom microbial interactions (e.g., bacteria-fungi and saprotrophic-symbiotic fungal guilds), which drive the reorganization of soil microbial communities into novel structural and functional configurations compared to monocultures (Garau et al., 2019; Guo et al., 2023). Nevertheless, mechanistic understanding of microbial network dynamics in mixed forest soils remains limited. Recent network analyses reveal strengthened positive edge correlations in mixed forest microbial co-occurrence patterns, implying tree species diversity preferentially enhances microbial cooperation over competitive exclusion (Ding et al., 2022).
P. massoniana, the predominant timber species and pioneer tree in subtropical China’s silviculture systems, has become extensively naturalized across this ecoregion due to its exceptional growth vigor, superior wood properties, and broad ecological plasticity. Its adaptive traits—including drought tolerance, low soil fertility tolerance, and phenotypic resilience—underpin both large-scale commercial plantations and critical ecosystem services like erosion control and degraded land rehabilitation (Tang et al., 2024; He et al., 2025; Tang et al., 2025). However, silvicultural practices emphasizing monoculture plantations of P. massoniana—particularly successive rotation regimes—have precipitated ecological degradation mirroring global monoculture syndromes: biodiversity depletion, constrained nutrient cycling, pathogen susceptibility, chronic decline in stand productivity, and compromised carbon sequestration capacity. These systemic dysfunctions fundamentally constrain multi-objective sustainable management frameworks for masson pine ecosystems (Wen et al., 2019; Li et al., 2014). Ecological restoration trials demonstrate that admixing native broadleaved species (e.g., Quercus acutissima) through near-natural silvicultural conversion enhances stand structural complexity, accelerates nutrient turnover, and reinstates carbon sink resilience (Zhao et al., 2016; Zhang et al., 2023).
The introduction of Q. acutissima into monoculture P. massoniana plantations induces significant alterations in litterfall dynamics, root biomass allocation, and rhizodeposition profiles (Zhang et al., 2023; Zhao et al., 2016). These shifts are hypothesized to drive stoichiometric reorganization in soil matrices and restructure microbial functional guilds across rhizospheric versus bulk soil compartments. This study addresses two critical knowledge gaps: (1) While mixed-species systems are known to enhance soil nutrient availability, their specific mechanisms in modulating microbial functional redundancy and keystone taxa selection remain unresolved; (2) Microbial community assembly drivers between root-influenced versus bulk soil zones under mixed-species regimes requires to be identified. This study investigates the ecological effects of P. massoniana-Q. acutissima mixed plantations by comparatively analyzing rhizosphere and bulk soils from three distinct forest management systems. Through comprehensive characterization of soil physicochemical properties (pH, nutrient availability) and microbial community profiles, we hypothesize that: (i) interspecific facilitation in mixed stands will enhance soil N/P availability and subsequently shift microbial community structure; (ii) Tree species mixing strengthens microbial network connectivity.
2 Materials and methods
2.1 Overview of the study area
The study was conducted at the Nanwan Experimental Forest Farm in Xinyang City, Henan Province, China (31°49′–32°14′N, 113°45′–114°10′E) (Figure 1). The forest farm spans 8,487 ha, with a forest coverage rate of 98.6% and a total standing timber volume of 876,000 m3. Dominant tree species include Cunninghamia lanceolata, Pinus massoniana, Quercus acutissima, and Liquidambar formosana. This region has a subtropical continental monsoon climate, featuring four distinct seasons: spring (dry, windy), summer (hot, rainy), autumn (cool), and winter (cold, snowy). Climatic data indicate an average annual temperature of 15.1°C and mean yearly precipitation of 1109.1 mm. The dominant soil type is yellow-brown soil.
Figure 1
2.2 Investigation object
The study focused on three forest types: pure P. massoniana stands, pure Q. acutissima stands, and mixed stands of P. massoniana and Q. acutissima (hereafter referred to as mixed forests). The P. massoniana pure stands are 47 years old, with a stand density of 625 trees·ha−1, mean tree height of 19.5 m, and mean diameter at breast height (DBH) of 26.2 cm. The Q. acutissima pure stands are 42 years old, with a stand density of 833 trees·ha−1, mean height of 13.3 m, and mean DBH of 18.1 cm. In the mixed stands (42 years old, 833 trees·ha−1), P. massoniana exhibits a mean height of 19.1 m and mean DBH of 24.9 cm, while Q. acutissima averages 12.8 m in height and 13.2 cm in DBH.
2.3 Soil sample collection
In late July 2024, five randomly selected plots were established in each of the three forest stands, maintaining a minimum spacing of 20 m between plots. Each plot measured 900 m2 (30 × 30 m), resulting in a total of 15 sampling plots. From each pure forest plot, five healthy, vigorously growing trees were randomly chosen, and both rhizosphere soil (0–30 cm depth) and non-rhizosphere soil samples were collected for analysis. Rhizosphere soil was collected via the shaking method with maintained sterility (Bai et al., 2024). Two soil types from the same plot were individually mixed, sieved (0.85 mm-sieve), split into two portions, and labeled. Part of the soil samples were maintained at ambient temperature for physical and chemical characterization, whereas the rest were immediately transferred to the lab and frozen at −80°C for microbial analysis. In the mixed forest plot, rhizosphere and non-rhizosphere soil samples were collected from the intermediate zone between P. massoniana and Q. acutissima (Supplementary Figure S1). Sampling procedures and analytical measurements were consistent with those used in pure stands. Three random sampling points per plot were selected for soil moisture determination via the ring knife method. Post-drying measurements revealed moisture contents of 15.37% (pure P. massoniana), 15.34% (pure Q. acutissima), and 16.71% (mixed forest), demonstrating slightly higher water retention in mixed stands.
2.4 Determination of soil chemical properties
Soil chemical properties were analyzed following the methodology outlined by the literature (Lu, 2002). The measured parameters included soil pH, organic carbon (SOC), total nitrogen (TN), total phosphorus (TP), total potassium (TK), alkali-hydrolyzable nitrogen (AN), available phosphorus (AP), and available potassium (AK). Soil pH was measured potentiometrically. Soil organic matter (SOM) was quantified via the potassium dichromate oxidation-colorimetric method, and SOC content was derived by dividing SOM by the conventional factor of 1.724. TN was determined using the Kjeldahl method, while TP and TK were measured via the molybdenum antimony colorimetric method and hydrofluoric-perchloric acid digestion, respectively. AN was analyzed using the alkaline hydrolysis diffusion method. AP was extracted with hydrochloric acid-ammonium fluoride, and AK was determined via ammonium acetate extraction followed by flame photometry.
2.5 Soil microbial determination
2.5.1 DNA extraction, library construction, and metagenomic sequencing
Total genomic DNA was extracted from 30 samples using the E.Z.N.A.® Soil DNA Kit (Omega Bio-tek, Norcross, GA, U.S.) according to manufacturer’s instructions. Concentration and purity of extracted DNA were determined with TBS-380 and NanoDrop2000, respectively. DNA extract quality was checked on 1% agarose gel.
DNA extract was fragmented to an average size of about 400 bp using Covaris M220 (Gene Company Limited, China) for paired-end library construction. Paired-end library was constructed using NEXTflex™ Rapid DNA-Seq (Bioo Scientific, Austin, TX, USA). Adapters containing the full complement of sequencing primer hybridization sites were ligated to the blunt-end of fragments. Paired-end sequencing was performed on Illumina Dnbseq-T7 (Illumina Inc., San Diego, CA, USA) at IGENEBOOK Biotechnology Co., Ltd. (Wuhan, China) using Dnbseq Reagent Kits according to the manufacturer’s instructions.1
2.5.2 Sequence quality control and genome assembly
The raw reads from metagenome sequencing were used to generate clean reads by removing adaptor sequences, trimming and removing low-quality reads (reads with N bases, a minimum length threshold of 50 bp and a minimum quality threshold of 20) using the fastp (Chen et al., 2018) (https://github.com/OpenGene/fastp, version 0.20.0). These high-quality reads were then assembled to contigs using MEGAHIT (Li et al., 2015) (parameters: kmer_min = 47, kmer_max = 97, step = 10) (https://github.com/voutcn/megahit, version 1.1.2), which made use of succinct de Bruijn graphs. Contigs with the length being or over 300 bp were selected as the final assembling result.
2.5.3 Gene prediction, taxonomy, and functional annotation
Open reading frames (ORFs) in contigs were identified using MetaGene (Noguchi et al., 2006).2 The predicted ORFs with length being or over 100 bp were retrieved and translated into amino acid sequences using the NCBI translation table.3
A non-redundant gene catalog was constructed using CD-HIT (Fu et al., 2012) (http://www.bioinformatics.org/cd-hit/, version 4.6.1) with 90% sequence identity and 90% coverage. Reads after quality control were mapped to the non-redundant gene catalog with 95% identity using SOAPaligner (Li et al., 2008) (http://soap.genomics.org.cn/, version 2.21), and gene abundance in each sample were evaluated.
Representative sequences of non-redundant gene catalog were annotated based on the NCBI NR database using blastp as implemented in DIAMOND v0.9.19 with e-value cutoff of 1e−5 using Diamond (http://www.diamondsearch.org/index.php, version 0.8.35) for taxonomic annotations (Buchfink et al., 2015). Cluster of orthologous groups of proteins (COG) annotation for the representative sequences were performed using Diamond (http://www.diamondsearch.org/index.php, version 0.8.35) against eggNOG database (version 4.5.1) with an e-value cutoff of 1e−5 (Buchfink et al., 2015). The KEGG annotation was conducted using Diamond (http://www.diamondsearch.org/index.php, version 0.8.35) against the Kyoto Encyclopedia of Genes and Genomes database (http://www.genome.jp/keeg/, version 94.2) with an e-value cutoff of 1e−5 (Buchfink et al., 2015).
2.6 Statistical analysis
Multiple comparisons of soil physicochemical properties, microbial relative abundance, and alpha diversity indices were conducted using two-way analysis of variance (ANOVA) and Tukey’s HSD test (p < 0.05). Variations in soil microbial community composition among treatments were evaluated through principal component analysis (PCA) (Vázquez-Baeza et al., 2013). The associations between microbial communities and soil environmental factors were determined via redundancy analysis (RDA) (Bokulich et al., 2018). Microbial co-occurrence patterns were investigated using network analysis. Functional predictions for soil bacteria and fungi were performed using PICRUSt2 and FUNGuild, respectively (Nguyen et al., 2016; Langille et al., 2013).
3 Results
3.1 Soil indicators
Among all measured soil parameters (excluding pH), the rhizosphere soil exhibited higher values compared to the non-rhizosphere soil (Figure 2). Within the same soil type, mixed forests demonstrate significantly greater TP, AP, and AN contents compared to pure stands—with pure PM forests ranking intermediate and pure QA forests displaying the lowest values (Figures 2D,F,G). Mixed forests demonstrated greater SOC, TN, TK, and AK concentrations than both single-species forests (Figures 2B,C,E,H). Compared to pure QA and PM stands, the mixed forest showed significant increases in rhizosphere soil nutrients—32.6 and 5.8% for AN, and 51.4 and 24.0% for AP, respectively (Figures 2F,G).
Figure 2
3.2 Composition of soil microbial community
The predominant bacterial phyla in the soil were Proteobacteria (29.1–39.1%) and Acidobacteria (28.2–36.5%), whereas the fungal community was dominated by Basidiomycota (66.1–84.2%), with Ascomycota (6.6–25.0%) as a secondary major group (Figure 3). Microbial populations in rhizosphere soil showed a greater relative abundance of Proteobacteria, Acidobacteria, Actinobacteria, Basidiomycota, and Ascomycota than those in non-rhizosphere soil (Figures 4A–C,E,F). Conversely, Verrucomicrobia displayed an opposite trend in their distribution between the two soil types (Figure 4D). The rhizosphere soil of mixed forests showed a greater relative abundance of Proteobacteria, Acidobacteria, Actinobacteria, and Basidiomycota than that of the two pure forest types (Figures 4A–C,E). Notably, Proteobacteria abundance increased by 20.5% (vs. QA) and 12.6% (vs. PM), while Actinobacteria increased by 15.6% (vs. QA) and 14.6% (vs. PM) in mixed forests (Figures 4A,C).
Figure 3
Figure 4
3.3 Community structure of soil microbes
In the bacterial community structure, the first two principal components (PC1 and PC2) explained 29.3 and 13.0% of the total variation, respectively (Figure 5A). For fungal communities, PC1 captured a substantially higher proportion (62.3%), followed by PC2 (16.2%) (Figure 5B). The bacterial and fungal community structures exhibited significant differences between rhizosphere and non-rhizosphere soils across all tree species. Additionally, distinct variations in microbial community composition were observed between mixed forests and pure forests, regardless of soil type (Figures 5A,B).
Figure 5
3.4 Diversity of soil microbes
The Shannon and Simpson index of rhizosphere soil microbes were significantly higher than those in non-rhizosphere soil (Figures 6A,B). The Shannon index of rhizosphere soil microbes exhibited increases of 11.6% (QA), 7.9% (PM), and 3.3% (Mix) compared to non-rhizosphere soil (Figure 6A). Similarly, the Simpson index showed respective increases of 4.9% (QA), 2.8% (PM), and 2.6% (Mix) (Figure 6B). However, the Shannon and Simpson index of soil microbial communities in mixed forests showed no significant differences compared to those in the two pure forests (p > 0.05) (Figures 6A,B).
Figure 6
3.5 Linkages between soil microbes and soil environmental factors
Redundancy analysis (RDA) revealed that the first two axes (RDA1 and RDA2) explained 40.15 and 18.53% of the total variance in the bacterial community, respectively (Figure 7A). For the fungal community, RDA1 accounted for 69.60%, while RDA2 contributed 18.62% of the explained variation (Figure 7B). The abundance of bacterial communities was significantly influenced by pH, SOC, TP, AN, and AP, whereas fungal community abundance was primarily affected by SOC, TP, AK, and AP (Figures 7, 8). There was a significant negative correlation between pH and SOC. In contrast, SOC was significantly positively correlated with both AN and AP. Moreover, TK and AK displayed a significant positive correlation (Figure 8).
Figure 7
Figure 8
3.6 Network analysis of soil microbes
The rhizosphere soil exhibited a higher number of bacterial connections compared to non-rhizosphere soil, with increases of 38.5% (QA), 28.6% (PM), and 71.2% (Mix), respectively (Figure 9). Across both soil types, mixed forests demonstrated greater bacterial network connectivity than pure forests, while PM forests maintained a higher number of connections than QA forests (Figure 9). Among all soil types, the predominant bacterial phyla—Proteobacteria, Acidobacteria, and Actinobacteria—serve as keystone taxa (Supplementary Table S1). Fungal network connectivity was significantly higher in rhizosphere soil compared to non-rhizosphere soil (Figure 10). However, the difference in fungal connectivity between Mix and pure QA stands was marginal. The dominant fungal phyla across all soil types included Basidiomycota, Ascomycota, and Mucoromycota (Supplementary Table S2).
Figure 9
Figure 10
3.7 Functional prediction of soil microbes
The core functional roles of soil bacteria encompassed cellular processes, environmental information processing, genetic information processing, human diseases, metabolism, and organismal systems. Among these, metabolism dominated, representing 50.5–52.1% of the total functional abundance (Figure 11). Notably, minimal functional variation was observed across different tree species and soil types (Figure 11). The secondary functions of soil bacteria primarily encompassed xenobiotics biodegradation, lipid metabolism, translation, membrane transport, signal transduction, prokaryotic cellular communities, cofactor and vitamin metabolism, energy metabolism, aminoacid metabolism, and carbohydrate metabolism (Figure 12). Among these, aminoacid metabolism and carbohydrate metabolism represented the dominant secondary functions. Notably, variations in secondary functional profiles across different tree species and soil types were minimal (Figure 12).
Figure 11
Figure 12
Soil fungi were primarily classified into three functional groups: pathotrophs, saprotrophs, and symbiotrophs. Saprotrophs were the most abundant, comprising 39.1–42.8% of the fungal community, followed by symbiotrophs (19.4–27.8%). Notably, variations in fungal functional composition between different tree species and soil types were minimal (Figure 13).
Figure 13
4 Discussion
4.1 Changes in soil nutrients
Soil carbon (C), nitrogen (N), and phosphorus (P) constitute fundamental nutrient elements that govern plant growth and mediate biogeochemical cycling at the plant–soil-microbe interface (Lu et al., 2021). The rhizosphere, as a critical zone of plant-environment interaction, exhibits distinct physicochemical properties from bulk soil through pronounced rhizosphere effects (Thacker and Quideau, 2021). Our findings demonstrated a reduction in rhizosphere pH relative to non-rhizosphere soil in mixed forests, aligning with observations by the literature of Xu et al. (2021). This acidification phenomenon primarily stemed from two mechanisms: (1) proton excretion by rhizosphere microbiota during metabolic activity, and (2) release of low-molecular-weight organic acids (e.g., citric, malic acids) through root exudation processes (Yuan et al., 2020). Our findings demonstrated that the mixed plantation of P. massoniana (coniferous) and Q. acutissima (broad-leaved) significantly enhanced soil nutrient pools (e.g., SOC, TN, and TP) compared to monocultures. These results aligned with the findings of Bai (2023) but contrasted with the report of Huang et al. (2025) on C. lanceolata-Castanopsis fissa mixed forests. The discrepancy may stem from stand developmental stages: Huang’s study examined 10-year-old Chinese fir during peak growth phases, where high nutrient demand coupled with reduced litter input likely created temporary nutrient depletion (Li et al., 2023). In our system, the substantial litter accumulation and accelerated decomposition rates in P. massoniana-Q. acutissima stands appeared to drive continuous nutrient replenishment through microbial-mediated mineralization processes.
4.2 Microbial community and composition
Recent studies demonstrated that Acidobacteria, a ubiquitous phylum in soil bacterial communities, exhibited substantial variations in relative abundance across different ecosystems (Huber et al., 2016). Members of this phylum possessed remarkable metabolic versatility, particularly in degrading recalcitrant carbohydrates through flexible carbon utilization pathways, thereby playing a pivotal role in soil organic matter cycling (Kalam et al., 2020). Notably, Acidobacteria demonstrated exceptional ecological competence in oligotrophic environments, as evidenced by their predominance in nutrient-poor soils through efficient resource acquisition strategies (Zelenev et al., 2005; Dennis et al., 2010). Our results demonstrated significantly higher relative abundances of Acidobacteria and Proteobacteria in mixed forest rhizospheres compared to pure stands. This enrichment likely reflected two synergistic mechanisms: (1) Enhanced acidobacterial populations promoting organic matter mineralization, thereby elevating nutrient availability (Bai, 2023); and (2) Plant-mediated selection through root exudation, which shaped distinct rhizosphere microbiomes (Garau et al., 2019). Notably, the copiotrophic Proteobacteria showed particularly strong rhizosphere enrichment, consistent with their known capacity for rapid utilization of labile root-derived carbon (Philippot et al., 2013). These findings suggested mixed-species planting created a positive feedback loop where diversified exudate profiles stimulate nutrient-cycling bacteria that in turn supported plant growth.
Our results demonstrated a significant enrichment of Basidiomycota in mixed forest soils compared to pure stands. This fungal predominance correlated strongly with elevated lignin content and acid-unhydrolyzable residue concentrations. As key decomposers of recalcitrant organic matter, Basidiomycota secreted a diverse array of carbohydrate-active enzymes (CAZymes), particularly lignin peroxidases and manganese peroxidases, which facilitated the breakdown of aromatic polymers (Rytioja et al., 2014). This enzymatic machinery not only drove litter decomposition but also promoted the formation of stable soil organic carbon (SOC) through humification processes, potentially explaining the higher SOC stocks observed in mixed stands. Empirical evidence indicated that introducing broad-leaved tree species into coniferous stands significantly enhanced Basidiomycota colonization in the rhizosphere zone, while concurrently altering fungal community structure (Li et al., 2021). These basidiomycetous fungi demonstrated remarkable phosphorus-mobilizing capacity through two synergistic mechanisms: (1) secretion of acid phosphatases and (2) production of low-molecular-weight organic acids (e.g., oxalic and citric acids) that chelated metal-phosphates (Taylor and Sinsabaugh, 2015). Consequently, the observed elevation in bioavailable phosphorus within mixed forest soils likely resulted from these microbially-mediated solubilization processes.
4.3 Microbial structure and function
Principal component analysis (PCA) revealed distinct clustering patterns among microbial communities, demonstrating significant structural differentiation. Bacterial communities exhibited stronger spatial separation than fungal communities across microhabitats, consistent with their higher sensitivity to microenvironmental variations (Li et al., 2023). Our data showed pronounced divergence in both bacterial and fungal community structures between rhizosphere and bulk soils. Furthermore, mixed-forest soils harbored unique microbial assemblages that significantly differed from pure stands, highlighting the dual regulatory effects of soil compartment and forest type on microbial community organization. Functional annotation of soil bacterial communities revealed six predominant metabolic guilds, though their relative abundances showed minimal variation across soil and forest types. Fungal functional groups were dominated by saprotrophic and symbiotic lifestyles across all stands. The saprotrophic fungi, as primary decomposers, drive ecosystem processes through: (1) lignocellulolytic enzyme production (e.g., laccases, lignin peroxidases) facilitating organic matter turnover, and (2) hyphal networks enhancing soil aggregation and xenobiotic degradation (Zhou et al., 2015; Sherman et al., 2014). Symbiotic fungi formed mutualistic associations characterized by: (i) extended nutrient acquisition networks, (ii) induced systemic resistance, and (iii) carbon-nutrient exchange (Xu et al., 2024; Ren et al., 2011; Wall et al., 2025). Notably, neither functional group exhibited significant abundance differences among forest or soil types, suggesting potential functional redundancy or unmeasured stabilizing factors requiring further study.
4.4 Network analysis
Network topological analysis revealed that microbial co-occurrence patterns reflected fundamental ecological relationships, with positive correlations suggesting potential mutualistic interactions or niche overlap, while negative correlations may indicate competitive exclusion (Faust and Raes, 2012; Ding et al., 2022). Our results demonstrated significantly higher network complexity in rhizosphere soils compared to bulk soils. Notably, mixed forest bacterial networks exhibited enhanced connectivity relative to pure conifer and broadleaf stands, supporting the hypothesis that tree species diversity promotes microbial interaction intensity. Network analysis revealed that increased topological complexity and intensified competitive interactions collectively enhanced the stability of rhizosphere microbial communities in mixed forests, as evidenced by higher robustness indices compared to pure stands (Yuan et al., 2021). Interestingly, the keystone taxa maintaining the interactions remained consistent across forest and soil types, suggesting potential functional redundancy or uncharacterized stabilizing mechanisms that warrant further investigation through targeted metatranscriptomics.
5 Limitations and prospects
This study had several limitations: (1) temporal variations in microbial functional dynamics were not assessed, (2) specific microbe-microbe interaction mechanisms remained uncharacterized, and (3) root exudate-mediated microbiome modulation was not quantitatively evaluated. Future research should integrate multi-omics approaches (metatranscriptomics and metabolomics) to elucidate functional gene expression patterns, implement in situ enzymatic activity profiling across seasons, and develop synthetic microbial communities for controlled gnotobiotic experiments to validate key drivers of observed community patterns.
6 Conclusion
This study confirmed that the mixed forest of P. massoniana and Q. acutissima significantly reshaped the microbial community structure by enhancing rhizosphere soil nutrients (TP, AP, AN, SOC, etc.). The key soil factors in the rhizosphere (SOC, TP, AP, AN) drived the enrichment of dominant bacterial communities (Proteobacteria, Acidobacteria) and fungal communities (Basidiomycota), and enhanced the strength of microbial network interactions. Although the changes in α-diversity were not significant, the unique β-diversity characteristics and more complex community structure of mixed forests indicated the formation of functionally specialized microbial communities. This change, by promoting nutrient cycling and enhancing ecological stability, confirmed that mixed forest transformation was an effective strategy for improving the soil ecological function of artificial forests.
Statements
Data availability statement
The data generated for this study are available on request to the corresponding author. Raw amplicon sequences were deposited in the Sequence Read Archive (SRA) and assigned the following BioProject accession number: PRJNA1270166.
Author contributions
LQ: Methodology, Funding acquisition, Writing – original draft, Conceptualization. ZG: Formal analysis, Writing – original draft, Investigation, Writing – review & editing. FR: Investigation, Formal analysis, Visualization, Writing – review & editing. TM: Funding acquisition, Writing – review & editing, Supervision.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Central Finance Forestry Science and Technology Promotion Project (GTH [2024]15), and Xinyang Agriculture and Forestry University Science and Technology Innovation Team Project (XNKJTD004).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1646535/full#supplementary-material
Footnotes
2.^http://metagene.cb.k.u-tokyo.ac.jp/
3.^http://www.ncbi.nlm.nih.gov/Taxonomy/taxonomyhome.html/index.cgi?chapter=tgencodes#SG1
References
1
BaiY. (2023). Research of soil microbial communities in mixed Pinus massoniana-Castanopsis hystrix plantations and pure plantations. Nanning, China: Guangxi University.
2
BaiY.LiuR.SuJ.ShenW. (2024). Effects of tree species mixing on bulk and rhizosphere soil microbial resource limitation in stands of Pinus massoniana and Castanopsis hystrix. Acta Ecol. Sin.44, 10770–10781. doi: 10.20103/j.stxb.202310312
3
BokulichN. A.KaehlerB. D.RideoutJ. R.DillonM.BolyenE.KnightR.et al. (2018). Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2's q2-feature-classifier plugin. Microbiome6:90. doi: 10.1186/s40168-018-0470-z
4
BuchfinkB.XieC.HusonD. H. (2015). Fast and sensitive protein alignment using diamond. Nat. Methods12, 59–60. doi: 10.1038/nmeth.3176
5
ChenS.ZhouY.ChenY.GuJ. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics34, i884–i890. doi: 10.1093/bioinformatics/bty560
6
DennisP. G.MillerA. J.HirschP. R. (2010). Are root exudates more important than other sources of rhizodeposits in structuring rhizosphere bacterial communities?FEMS Microbiol. Ecol.72, 313–327. doi: 10.1111/j.1574-6941.2010.00860.x
7
DingK.ZhangY.YrjalaK.TongZ.ZhangJ. (2022). The introduction of Phoebe bournei into Cunninghamia lanceolata monoculture plantations increased microbial network complexity and shifted keystone taxa. For. Ecol. Manag.509:120072. doi: 10.1016/j.foreco.2022.120072
8
FaustK.RaesJ. (2012). Microbial interactions: from networks to models. Nat. Rev. Microbiol.10, 538–550. doi: 10.1038/nrmicro2832
9
FuL.NiuB.ZhuZ.WuS.LiW. (2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics28, 3150–3152. doi: 10.1093/bioinformatics/bts565
10
GanD.FengJ.HanM.ZengH.ZhuB. (2021). Rhizosphere effects of woody plants on soil biogeochemical processes: a meta-analysis. Soil Biol. Biochem.160:108310. doi: 10.1016/j.soilbio.2021.108310
11
GarauG.MorillasL.RoalesJ.CastaldiP.MangiaN. P.SpanoD.et al. (2019). Effect of monospecific and mixed Mediterranean tree plantations on soil microbial community and biochemical functioning. Appl. Soil Ecol.140, 78–88. doi: 10.1016/j.apsoil.2019.04.005
12
GilliamF. S. (2025). Linkage of strata of forest vegetation with forest soil microbiomes: a review. Front. Microbiol.16:1575691. doi: 10.3389/fmicb.2025.1575691
13
GilliamF. S.HargisE. A.RabinowitzS. K.DavisB. C.SweetL. L.MossJ. A. (2023). Soil microbiomes of hardwood- versus pine-dominated stands: linkage with overstory species. Ecosphere14:e4537. doi: 10.1002/ecs2.4537
14
GuoJ.WeiB.LiuJ.EissenstatD. M.YuS.GongX.et al. (2023). Linkages between plant community composition and soil microbial diversity in masson pine forests. Plan. Theory12:1750. doi: 10.3390/plants12091750
15
HanM.ChenY.SunL.YuM.LiR.LiS.et al. (2023). Linking rhizosphere soil microbial activity and plant resource acquisition strategy. J. Ecol.111, 875–888. doi: 10.1111/1365-2745.14067
16
HeQ.HuangH.LiJ.LiangY.MaH.HuangX.et al. (2025). Effects of uneven-aged mixed transformation of Pinus massoniana plantation on understory plant diversity and functional groups. Chin. J. Ecol. 44, 1–16.
17
HuY.ZhaiM.WuJ.JiaL. (2002). Dynamies of amount and bio-chemical activity of microorganism in the soil rhizosphere of the pure and mixed stands of poplar and black locust. Chin. J. Soil Sci.33, 219–222. doi: 10.19336/j.cnki.trtb.2002.03.015
18
HuangA.JiL.HeG.NiY.YangL.WangM. (2025). Ecological stoichiometric characteristics and their rhizosphere effects of soil in mixed forests of Cunninghamia lanceolata and Phoebe bournei. Chin. J. Ecol. 44, 1–11.
19
HuberK. J.GeppertA. M.WannerG.FselB. U.WstP. K.OvermannJ. (2016). The first representative of the globally widespread subdivision 6 acidobacteria, Vicinamibacter silvestris gen. Nov., sp. nov., isolated from subtropical savannah soil. Int. J. Syst. Evol. Microbiol.66, 2971–2979. doi: 10.1099/ijsem.0.001131
20
KalamS.BasuA.AhmadI.SayyedR. Z.El-EnshasyH. A.DailinD. J.et al. (2020). Recent understanding of soil acidobacteria and their ecological significance: a critical review. Front. Microbiol.11:580024. doi: 10.3389/fmicb.2020.58
21
KuzyakovY.RazaviB. S. (2019). Rhizosphere size and shape: temporal dynamics and spatial stationarity. Soil Biol. Biochem.135, 343–360. doi: 10.1016/j.soilbio.2019.05.011
22
LangilleM. G. I.ZaneveldJ.CaporasoJ. G.McDonaldD.KnightsD.ReyesJ. A.et al. (2013). Predictive functional profiling of microbial communities using 16s rRNA marker gene sequences. Nat. Biotechnol.31, 814–821. doi: 10.1038/nbt.2676
23
LiM.DuM.NieC. (2014). Research advances in soil improvement of Pinus massoniana plantation. World For. Res.27, 31–36. doi: 10.13348/j.cnki.sjlyyj.2014.05.006
24
LiW.HuangY.ChenF.LiuY.LinX.ZongY.et al. (2021). Mixing with broad-leaved trees shapes the rhizosphere soil fungal communities of coniferous tree species in subtropical forests. For. Ecol. Manag.480:118664. doi: 10.1016/j.foreco.2020.118664
25
LiR.LiY.KristiansenK.WangJ. (2008). Soap: short oligonucleotide alignment program. Bioinformatics24, 713–714. doi: 10.1093/bioinformatics/btn025
26
LiD.LiuC.LuoR.SadakaneK.LamT. (2015). Megahit: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de bruijn graph. Bioinformatics31, 1674–1676. doi: 10.1093/bioinformatics/btv033
27
LiW.SunH.CaoM.WangL.FangX.JiangJ. (2023). Diversity and structure of soil microbial communities in Chinese fir plantations and Cunninghamia lanceolata–Phoebe bournei mixed forests at different successional stages. Forests14:1977. doi: 10.3390/f14101977
28
LingN.WangT.KuzyakovY. (2022). Rhizosphere bacteriome structure and functions. Nat. Commun.13:836. doi: 10.1038/s41467-022-28448-9
29
LiuL.DuanZ.XuM.HuJ.WangS.HuZ.et al. (2010). Effect of monospecific and mixed Cunninghamia lanceolata plantations on microbial community and two functional genes involved in nitrogen cycling. Plant Soil327, 413–428. doi: 10.1007/s11104-009-0067-1
30
LuR. (2002). Agricultural and chemistry analysis of soil. Beijing: China Agricultural Science and Technology Press.
31
LuX.VitousekP. M.MaoQ.GilliamF. S.LuoY.TurnerB. L.et al. (2021). Nitrogen deposition accelerates soil carbon sequestration in tropical forests. Proc. Natl. Acad. Sci. USA118:e2020790118. doi: 10.1073/pnas.2020790118
32
MaY.YueK.HedenecP.LiC.LiY.WuQ. (2023). Global patterns of rhizosphere effects on soil carbon and nitrogen biogeochemical processes. Catena220, 1843–1852. doi: 10.1016/j.catena.2022.106661
33
NguyenN. H.SongZ.BatesS. T.BrancoS.TedersooL.MenkeJ.et al. (2016). FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol.20, 241–248. doi: 10.1016/j.funeco.2015.06.006
34
NoguchiH.ParkJ.TakagiT. (2006). MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. Nucleic Acids Res.34, 5623–5630. doi: 10.1093/nar/gkl723
35
PereiraA. P. A.DurrerA.GumiereT.GoncalvesJ. L. M.RobinA.BouilletJ.et al. (2019). Mixed eucalyptus plantations induce changes in microbial communities and increase biological functions in the soil and litter layers. For. Ecol. Manag.433, 332–342. doi: 10.1016/j.foreco.2018.11.018
36
PhilippotL.RaaijmakersJ. M.LemanceauP.Van-der PuttenW. H. (2013). Going back to the roots: the microbial ecology of the rhizosphere. Nat. Rev. Microbiol.11, 789–799. doi: 10.1038/nrmicro3109
37
RenA. Z.LiX.HanR.YinL. J.WeiM. Y.GaoY. B. (2011). Benefits of a symbiotic association with endophytic fungi are subject to water and nutrient availability in Achnatherum sibiricum. Plant Soil346, 363–373. doi: 10.1007/s11104-011-0824-9
38
RytiojaJ.HildenK.YuzonJ.HatakkaA.de-VriesR. P.MakelaM. R. (2014). Plant-polysaccharide-degrading enzymes from basidiomycetes. Microbiol. Mol. Biol. Rev.78, 614–649. doi: 10.1128/MMBR.00035-14
39
ShermanC.GrishkanI.BarnessG.SteinbergerY. (2014). Fungal community-plant litter decomposition relationships along a climate gradient. Pedosphere24, 437–449. doi: 10.1016/S1002-0160(14)60030-8
40
SmithW. H. (1976). Character and significance of forest tree root exudates. Ecology57, 324–331. doi: 10.2307/1934820
41
TangB.XuZ.YangZ.TangS.LuX.LuK.et al. (2025). Effect of pruning intensity on growths of Pinus massoniana at different stand ages. Guangxi For. Sci.54, 167–172. doi: 10.19692/j.issn.1006-1126.20250206
42
TangS.YinM.MengC.QinY.PanZ.YangZ. (2024). Effects of different pruning intensities on growths of Pinus massoniana young plantations. Guangxi For. Sci.53, 763–768. doi: 10.19692/j.issn.1006-1126.20240610
43
TaylorD. L.SinsabaughR. (2015). The soil fungi, occurrence, phylogeny and ecology. Soil Microbiol. Ecol. Biochem.4, 77–110.
44
ThackerS. J.QuideauS. A. (2021). Rhizosphere response to predicted vegetation shifts in boreal forest floors. Soil Biol. Biochem.154:108141. doi: 10.1016/j.soilbio.2021.108141
45
Vázquez-BaezaY.PirrungM.GonzalezA.KnightR. (2013). EMPeror: a tool for visualizing high-throughput microbial community data. Gigascience2:16. doi: 10.1186/2047-217X-2-16
46
WallC. B.KajiharaK.RodriguezF. E.VilonenL.YogiD.SwiftS. O. I.et al. (2025). Symbiotic fungi alter plant resource allocation independent of water availability. Am. J. Bot.112:e16459. doi: 10.1002/ajb2.16459
47
WenY.LiH.ZhouX.ZhuH.LiY.CaiD.et al. (2019). Effects of uneven-aged Pinus massoniana x Castanopsis hystrix mixed plantations on structural and functions of soil microbial community. Guangxi Sci.26, 188–198. doi: 10.13656/j.cnki.gxkx.20190419.001
48
WuZ.DuanX.LiW.ChenF.LiuY.FangX. (2022). Effects of tree mixture on rhizosphere soil nitrogen mineralization and microbial characteristics of coniferous trees in subtropical plantations. Acta Ecol. Sin.42, 8414–8424.
49
WuD.ZhangM.PengM.SuiX.SunG. (2019). Variations in soil functional fungal community structure associated with pure and mixed plantations in typical temperate forests of China. Front. Microbiol.10:1636. doi: 10.3389/fmicb.2019.01636
50
XuH.DuH.ZengF.SongT.PengW. (2021). Diminished rhizosphere and bulk soil microbial abundance and diversity across succession stages in karst area, Southwest China. Appl. Soil Ecol.158:103799. doi: 10.1016/j.apsoil.2020.103799
51
XuL.HeJ.MengY.ZhengY.LuB.ZhangJ.et al. (2024). Enhancing drought resistance in Pinus tabuliformis seedlings through root symbiotic fungi inoculation. Front. Plant Sci.15:1446437. doi: 10.3389/fpls.2024.1446437
52
YuanY.DaiX.FuX.KouL.LuoY.JiangL.et al. (2020). Differences in the rhizosphere effects among trees, shrubs and herbs in three subtropical plantations and their seasonal variations. Eur. J. Soil Biol.100:1. doi: 10.1016/j.ejsobi.2020.103218
53
YuanM. M.GuoX.WuL.ZhangY.ZhouJ. (2021). Climate warming enhances microbial network complexity and stability. Nat. Clim. Chang.11, 343–348. doi: 10.1038/s41558-021-00989-9
54
ZelenevV. V.Van-BruggenA.SemenovA. M. (2005). Modeling wave-like dynamics of oligotrophic and copiotrophic bacteria along wheat roots in response to nutrient input from a growing root tip. Ecol. Model.188, 404–417. doi: 10.1016/j.ecolmodel.2005.01.046
55
ZhangS.YeL.ZhouY.WangX.XuY.JiangJ.et al. (2023). Water use sources and its influencing factors of Pinus massoniana and Quercus acutissima community in hilly region of southern China. Chin. J. Appl. Ecol.34, 1729–1736. doi: 10.13287/j.1001-9332.202307.010
56
ZhaoJ.YangX.XuX.ChenJ.KongY. (2016). Spatial distribution characteristic of soil dissolved organic carbon and nitrogen under the Pinus massoniana-Quercus acutissima mixed forest. J. Soil Water Conserv.30, 213–219. doi: 10.13870/j.cnki.stbcxb.2016.03.037
57
ZhouY.ClarkM.SuJ.XiaoC. (2015). Litter decomposition and soil microbial community composition in three Korean pine (Pinus koraiensis) forests along an altitudinal gradient. Plant Soil386, 171–183. doi: 10.1007/s11104-014-2254-y
Summary
Keywords
rhizosphere, microbial communities, mixed forests, microbial function, Pinus massoniana-Quercus acutissima
Citation
Qiao L, Guan Z, Ren F and Ma T (2025) Comparative analysis of rhizosphere microbial communities in monoculture and mixed oak–pine forests: structural and functional insights. Front. Microbiol. 16:1646535. doi: 10.3389/fmicb.2025.1646535
Received
18 June 2025
Accepted
14 July 2025
Published
25 July 2025
Volume
16 - 2025
Edited by
Fred O. Asiegbu, University of Helsinki, Finland
Reviewed by
Xiangshi Kong, Jishou University, China
Frank S. Gilliam, University of West Florida, United States
Updates
Copyright
© 2025 Qiao, Guan, Ren and Ma.
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: Zhuizhui Guan, 20201159005@stu.scau.edu.cn; Tianxiao Ma, 646451368@qq.com
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.