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

ORIGINAL RESEARCH article

Front. Microbiol., 06 November 2025

Sec. Terrestrial Microbiology

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

This article is part of the Research TopicSoil microbial ecology and soil ecosystem functions for addressing climate change and maintaining ecological sustainabilityView all 6 articles

Soil core microbial taxa maintain community resistance to drive soil ecosystem multifunctionality under Alternanthera philoxeroides invasion

Xiao Lin&#x;Xiao LinXinyu Zhang&#x;Xinyu ZhangShujing LiShujing LiYongming WangYongming WangHaiyang YuHaiyang YuWenxiu Qin
Wenxiu Qin*Jiaoyang Zhang
&#x;Jiaoyang Zhang*
  • College of Resources and Environment, Anhui Agricultural University, Hefei, China

Alien plant invasion exerts profound impacts on local ecosystems, like biodiversity, stability, and overall functionality. Understanding the ecological stability of soil microbial communities is essential for elucidating how plant invasions drive changes in soil ecosystem functions. In this study, soil ecosystem multifunctionality (EMF) was examined in relation to the resistance of microbial communities across five Alternanthera philoxeroides–invaded sites (Xiaoxian, Hefei, Tongling, Anqing, and Huangshan) within different ecological functional zones in Anhui Province, China, with contrasting soil properties and environmental conditions. Soil samples were collected from invaded and uninvaded plots, microbial core taxa were identified, and structural equation modeling was applied to examine the relationships among invasion, microbial traits, and ecosystem functions. The results showed that invasion significantly increased soil EMF at Xiaoxian by 34.2% and at Hefei by 24.1%, which was primarily by increasing in nitrogen- and phosphorus-related functions. Bacterial communities exhibited consistently higher resistance than fungal communities, owing to their more complex interaction networks. Core bacterial taxa were positively associated with microbial resistance, supporting nutrient cycling and maintaining ultimately EMF. Structural equation modeling revealed that invasion-driven changes in soil properties influenced EMF indirectly through bacterial community resistance mediated by core microbial taxa. These findings demonstrated that bacterial resistance, underpinned by the stabilizing role of core microbial taxa, played a central role in maintaining EMF during plant invasion. The results highlight the importance of core microbial taxa as ecological stabilizers and provide new insights into the belowground mechanisms linking biological invasions to ecosystem resilience.

1 Introduction

Invasive alien plants have attracted increasing global attention due to their profound impacts on terrestrial biodiversity and stability (Lai et al., 2025). These impacts extend to key ecosystem processes, like biogeochemical cycling. For instance, the invasion of Ligustrum in the Andes significantly accelerated biomass accumulation in successional forests, thereby enhancing carbon sequestration (Nguyen et al., 2020). At a broader scale, a global meta-analysis revealed that invasive plants could exhibit higher rhizosphere nitrogen mineralization rates, greater soil available nitrogen content, and increased activity of nitrogen-releasing enzymes compared with native species, ultimately enhancing soil nitrogen cycling (Negesse et al., 2025). In addition to altering carbon and nitrogen dynamics, invasive species could influence phosphorus turnover. For example, Bidens pilosa showed approximately 30% higher alkaline phosphomonoesterases activity in the rhizosphere than native plants, directly promoting the decomposition of organic phosphorus (Sun et al., 2022). While such studies provide valuable insights into individual soil element cycling function, they often overlook the reality that ecosystems operate through multiple, interconnected and trade-off functions. Understanding how these functions occur simultaneously is essential for capturing the full complexity of ecosystem responses. Collectively, biogeochemical cycling serves as a fundamental basis for ecosystem functions, underpinning processes such as productivity, nutrient turnover, and carbon sequestration. EMF provides an integrative framework to evaluate these processes, enabling a more comprehensive assessment of invasion consequences by capturing both trade-offs and synergies among ecosystem functions. This system-level perspective not only links community structural changes with ecosystem processes but also offers valuable guidance for ecosystem management and restoration under invasion pressure. Yet, comprehensive evaluations of EMF that incorporate multiple, interrelated soil functions are still lacking, leaving the ecosystem-level consequences of plant invasion poorly understood. Addressing this knowledge gap, the present study examines how plant invasion affects soil EMF by integrating multiple biogeochemical and microbial properties, thereby providing a comprehensive assessment of ecosystem responses to invasion.

Soil microbial communities, with their highly complex species composition, are sensitive and adaptable to environmental change due to traits such as metabolic flexibility, large population sizes, and rapid adaptation (Allison and Martiny, 2008). These attributes underpin their resistance, the capacity to maintain structural, compositional, and functional stability under stress (Shade et al., 2012), a key aspect of ecosystem stability. However, resistance can vary with disturbance intensity and environmental context, and may be altered by global change drivers such as climate, land use, and nutrient availability (Jiao et al., 2022; Zheng et al., 2025; Zhou et al., 2020). Plant invasion represents a major biotic disturbance that can reshape microbial stability and community composition. Invasive plants influence soil microorganisms by plant secondary metabolites or changes in higher trophic interactions, referring to resource-driven influences such as enhanced nutrient acquisition and root exudation (Xiao et al., 2014; Zhang et al., 2020). These processes can accelerate bacterial community succession, reduce stochasticity in community assembly, and lead to uneven shifts in microbial composition (Zhao et al., 2024). For example, invasive plants enhance nitrogen turnover through interactions with rhizosphere microorganisms, thereby strengthening their competitive advantage (Ji et al., 2024). And the invasion of Mikania micrantha has been shown to alter microbial metabolic functions in subtropical forests, creating conditions more favorable for its own growth (Zhao et al., 2023). Yet, it remains unclear how such invasion-driven changes in microbial resistance influence soil ecosystem functions.

Most research on invasion–microbe interactions have emphasized overall diversity or community composition, while largely overlooking the contributions of specific functional taxa, which directly underpin key ecosystem processes. Microbial communities are structured into distinct taxa with differing ecological importance, among which core and keystone microbial taxa are particularly critical for maintaining function and network stability (Jacquiod et al., 2020). These taxa are consistently present across diverse habitats, often exhibit stable genetic traits, contribute disproportionately to multiple microbial functions (Jiao et al., 2021; Shade and Handelsman, 2012; Sun et al., 2024), and occupy broader ecological niches (Scheller and Parajuli, 2018; Vietorisz et al., 2025), thereby enhancing ecosystem resilience (Coban et al., 2022; Jiao et al., 2022). Alterations in these taxa may play a pivotal role in mediating the impacts on EMF. For example, shifts in arbuscular mycorrhizal fungi, a keystone microbial taxa in farmland and saline-alkali wasteland ecosystems, can alter soil nutrient supply (Zhu et al., 2025), while changes in denitrifying or cellulose degrading bacteria can influence nitrogen retention and organic matter turnover (Hu et al., 2019; Luo et al., 2020), ultimately affecting both the resistance and functional stability of soil ecosystems. Nevertheless, empirical evidence from invasive ecosystems remains limited (He et al., 2025). Introducing the perspective of core microbial taxa into invasion ecology is therefore novel, as it moves beyond documenting diversity loss or compositional shifts to reveal the stabilizing mechanisms that buffer invasion-induced disturbances. Unlike agricultural or forest ecosystems, where core taxa are often studied in relation to long-term management practices or successional stability, invaded ecosystems provide a unique context in which core taxa mediate sudden, non-equilibrium disturbances caused by exotic plants, thereby shaping soil EMF. Understanding how core and keystone microbial taxa respond to invasion and contribute to EMF is therefore essential for elucidating the mechanisms that link above- and belowground processes.

To address this gap, we focused on the perennial invasive species Alternanthera philoxeroides, which can alter soil ecosystem processes through complex interactions with soil microbial communities (Fan et al., 2016). Plant and soil samples were collected from abandoned fields (>2 years) that had been invaded by A. philoxeroides, together with adjacent non-invaded fields serving as controls. Through integrated analyses, we aimed to combine invasion, changes in EMF, and shifts in core microbial taxa. Specifically, we tested the following hypotheses: (1) A. philoxeroides invasion enhances soil EMF; (2) core and keystone microbial taxa, particularly bacteria, are critical for maintaining soil microbial stability in response to invasion; (3) the core microbial taxa and the resistance they mediate are closely related to the improvement of soil EMF. Clarifying these relationships is essential for understanding the mechanisms by which biological invasions reshape belowground ecosystems and for informing strategies to preserve or restore soil EMF under global change.

2 Materials and methods

2.1 Site description

To maximize regional representativeness, we selected the Yangtze–Huaihe River Basin, a north–south climatic transition zone with pronounced edaphoclimatic heterogeneity. The study area spans five ecological sub regions: (i) The Plains Ecological Zone north of the Huaihe River, (ii) The Hilly and Ridge Ecological Zone between the Yangtze and Huaihe Rivers, (iii) The Dabie Mountain Ecological Zone in western Anhui, (iv) The Yangtze River Plain Ecological Zone, and (v) The Southern Anhui Mountainous and Hilly Ecological Zone. Based on prior work and field reconnaissance, we confirmed that A. philoxeroides are widely invasive across these zones. The locations of the sampling sites are shown in Figure 1a, and their edaphoclimatic contrasts are summarized in supporting information.

Figure 1
Map with labeled regions and invasion indicators (a), bar chart comparing soil ecosystem multifunctionality under invasion versus non-invasion across five sites (b), and a triangular plot illustrating the multifunctionality index distribution with color-coded sites (c).

Figure 1. The changes of soil functional indicators after the invasion of Alternanthera philoxeroides. (a) The position information of different sites; (b) The histogram of the effect of A. philoxeroides invasion on soil multifunctionality; (c) The ternary phase diagram of C, N, P cycling function of A. philoxeroides invasion background soil ecosystem, all functional values were standardized using the “Z-score” method.

Xiaoxian, located in a warm-temperate semi-humid monsoon climate and dominated by sandy soil, had the lowest mean annual precipitation (MAP, 854 mm) and mean annual temperature (MAT, 14.3 °C). Hefei lies in a subtropical humid monsoon climate with yellow-brown soils and received 1,000 mm MAP, with a MAT of 15.7 °C. Tongling, situated in a mid-subtropical humid monsoon climate dominated by yellow-brown and brownish-red soils, recorded 1,412 mm MAP and a MAT of 15.9 °C. Anqing, positioned at the transition between subtropical humid and warm-temperate semi-humid zones, with yellow-brown, yellowish-red, and mountain meadow soils, received 1,427 mm MAP and had a MAT of 14.5 °C. Huangshan, located in the transition from northern to central subtropical climate and characterized by red, yellow, dark yellow-brown, and mountain meadow soils, recorded the highest MAP (1706 mm) with a MAT of 16.1 °C.

2.2 Experimental design and soil samples collection

After finalizing the sampling sites, we conducted fieldwork in September–October 2022 across five ecological zones. In each zone, we established four paired sites—A. philoxeroides–invaded versus uninvaded controls—on land that had been abandoned for >2 years (The land had been used for tillage before being abandoned, after which A. philoxeroides began to invade). Each sample site was no more than 30 m between the A. philoxeroides invasion-affected and non-invasion-affected sample plots to minimize the effects of soil background differences. Four pairs of sample plots were separated by a distance of no less than 50 m. Each sample plot was 1 × 1 m2 and was used for investigating the basic characteristics of the vegetation and soil sample collection. The aboveground portions of plants within each small quadrat were harvested at ground level and sorted by species. Each species was then stored separately in envelopes and brought back to the laboratory for the determination of vegetation aboveground biomass. Soil samples from 0 to 20 cm depth were collected at multiple points in an S-shaped pattern within each quadrat. Soil samples from the same quadrat were mixed, and stones, plant roots, litter, and visible soil fauna were removed. The fresh soil was then divided into two portions. One portion was air-dried and sieved through a 0.25 mm sieve for the determination of basic soil physicochemical properties. The other fresh sample was sieved through a 2 mm sieve and further divided into two subsamples. One of these subsamples was stored in a 4 °C refrigerator for the measurement of microbial biomass carbon, nitrogen, and phosphorus, soluble organic carbon, mineralization rate, soil phosphatase, β-glucosidase, and N-acetyl-β-glucosidase. The other subsample was stored in a −80 °C freezer for DNA extraction followed by microbial community characterization by marker gene analysis.

2.3 Determination of soil physical and chemical properties

Soil samples were analyzed for various parameters using specific methods. The potassium dichromate oxidation-external heating method was used to determine soil organic carbon (SOC) content (Nóbrega et al., 2015), the Kjeldahl method for total nitrogen (TN) (Bremner and Mulvaney, 1982), as well as the sodium hydroxide alkali fusion-molybdenum antimony colorimetric method for total phosphorus (TP) contents (Schade et al., 2003), respectively. To measure soil ammonium nitrogen (NH4+–N) and nitrate nitrogen (NO3–N) content, extraction was performed using a KCl solution (1 mol·L−1). Soil pH was determined using solutions with a ratio of 1:2.5 (soil:water), the pH electrode type is a glass electrode, and the calibration temperature is 25 °C. Soil nitrogen mineralization rate (NMR) was determined by anaerobic method (Morecroft et al., 1992). Soil available phosphorus (AP) was extracted using NaHCO3 and measured with a spectrophotometer (Olsen, 1954). Additionally, Soil microbial biomass, including soil microbial biomass carbon, nitrogen, phosphorus (MBC, MBN, MBP), was determined by chloroform fumigation method (Harris et al., 1997; Oberson et al., 1997).

2.4 Soil enzymatic activity analysis

In this study, soil enzyme activity analysis included soil phosphatase, soil β-glucosidase (S-β-G-C), and soil N-acetyl-β-D-glucosidase (NAG). All of them were detected by Beijing Sheng gong soil enzyme activity detection kit, the unit of enzyme activity is U/g (1 μmol of p-nitrophenol is produced per gram of soil every day). Soil phosphatase activity detection method: 0.1 g soil sample was taken in 2 mL centrifuge tube, 0.05 mL toluene was added, 15 min later, substrate was added, and 24 h was catalyzed in 37 °C environment. After mixing and centrifugation, the supernatant was added to 96 microporous plates, and the absorbance value was measured at 660 nm wavelength by multifunctional microplate reader, and the enzyme activity index was calculated according to the standard curve. Soil β-glucosidase and soil N-acetyl-β-D-glucosidase detection method, 30 mg of soil was taken in a 1.5 mL centrifuge tube, and the corresponding substrates were added, respectively. The reaction was catalyzed at 37 °C for 1 h, and then immediately taken out and placed in a boiling water bath for 5 min. The process must ensure that the centrifuge tube is closed. After natural cooling, the supernatant was centrifuged and added to a 96-well microplate. The absorbance was measured at a wavelength of 400 nm using a multifunctional microplate reader, and the enzyme activity index was calculated according to the standard curve. In the above method, each soil sample needs to be set up a control group and a blank group.

2.5 Microbial community analysis

To isolate microbial DNA from soil samples, 0.25 g of soil were processed using the TIANamp Soil DNA Kit (Tiangen Biotech Co., Ltd., Beijing, China). The purity and concentration of the extracted DNA were rigorously evaluated using a spectrophotometer. For the purpose of amplifying bacterial and fungal genes, the specific regions of interest were the V3–V4 region of bacterial genes and the ITS1 region of fungal genes. A polymerase chain reaction (PCR) was conducted with tailored primers: 341F/806R (Wang et al., 2019) for bacteria and 1737F/2043R (Hao et al., 2016) for fungi. The PCR reaction mixture was composed of 2 μL of sterile ultrapure water, 15 μL of Phusion Master Mix (2× concentration), 3 μL of 6 μM primers, and 10 μL of template DNA (ranging from 5 to 10 ng). The amplification protocol initiated with a denaturation step at 98 °C for 60 s, followed by 30 cycles, each comprising denaturation at 98 °C for 10 s, annealing at 50 °C for 30 s, and extension at 72 °C for 30 s. A final extension step was executed at 72 °C for 5 min. Subsequently, the PCR products were purified with a Qiagen Gel Extraction Kit (Qiagen, Germany) and subjected to sequencing on the Illumina HiSeq2500 platform (Majorbio Bio-Pharm, Shanghai).

The raw sequence data underwent rigorous processing and chimera screening utilizing QIIME (Caporaso et al., 2010). During this step, chimeric sequences were identified and excluded, while the remaining sequences were clustered into operational taxonomic units (OTUs) based on a 97% similarity threshold using UPARSE (Edgar, 2013). Taxonomic assignment of the OTUs was then performed utilizing the Green Genes Database 13_5 and the RDP classifier (DeSantis et al., 2006).

2.6 Calculation of EMF

In assessing the functional indices of soil ecosystems, we selected 14 ecosystem functional indices modulated by soil microorganisms. These were consolidated into three primary functional indices pertaining to soil carbon, nitrogen, and phosphorus cycling, respectively. Specifically, the soil C multifunctionality (C_Mfa) indices encompassed SOC, DOC, S-β-G-C, and MBC. Similarly, the soil N multifunctionality (N_Mfa) indices comprised TN, NH4+–N, NO3–N, NMR, NAG, and MBN. Lastly, the soil P multifunctionality (P_Mfa) indices included TP, AP, MBP, and SAP. Additionally, soil pH, soil moisture content (SMC), and aboveground biomass (AGB), which are recognized as factors influencing soil ecosystem function, were measured and subsequently subjected to a joint analysis alongside the multifunctional indices of the ecosystem. Prior to the calculation of emf, all functional values were standardized using the “Z-score” method, i.e., standardized using the mean and standard deviation (measured values minus the mean divided by the standard deviation) (Maestre et al., 2012). Afterward, the multifunctionality index MFa was calculated using the following formula: according to the following formula:

MFa=1Fi=1Fg(ri(fi))

Here, F refers to the amount of functions determined; fi refers to the determined value of function i; ri is a mathematical function that converts fi into a positive value; and g refers to the standardization of all determined values.

2.7 Statistical analysis

Microbial community resistance is typically defined as the extent to which microorganisms exhibit tolerance in response to disturbances. In this study, microbial community resistance was specifically quantified by assessing the similarity of microbial communities between invaded and non-invaded sites within each ecological zones, utilizing the “vegan” package in R 4.3.1 for data processing.

Microbial network diagrams can assist us in better observing the changes in microbial community structure in response to disturbances. By screening the top 500 species at the genus level of bacteria and fungi, we can select OTUs, species, genera, and phyla with significant interaction effects, as well as the interaction information among them. The data are then processed in Gephi 0.9.7 for visual analysis of the interaction relationships within the microbial community.

Utilizing R, the raw OTU files of bacteria and fungi are processed to identify eligible microbial functional groups. Core microbial taxa are selected based on OTUs that account for the top 10% of abundance across all samples and are ubiquitously present in 50% or more of soil samples, thus recognized as core microbial taxa. Key taxa are identified by first calculating the Spearman rank correlation coefficients among OTUs to form a relationship matrix, followed by selecting OTUs with absolute correlation coefficients ≥ 0.8, along with their corresponding genera and phyla, which are then recognized as key taxa.

Mantel analysis facilitates a deeper understanding of the relationships between different microbial functional groups and environmental factors. Mantel test was conducted in R using the “linkET,” “ggplot2,” and “dplyr” packages, and the results were visualized graphically. Performing regression analysis between microbial communities and ecosystem functions, facilitated by SPSS 22 for data processing and Origin 2019 for plotting regression models, can enhance our understanding of the correlations between individual communities and specific ecosystem functions. Furthermore, structural equation modeling (SEM) assesses the linkages among environmental factors, soil microorganisms, and ecosystem functions, involving data processing in SPSS 22 and subsequent model construction in SPSS Amos 26.

3 Results

3.1 Response of soil EMF and elemental cycling related function effects of A. philoxeroides invasion

Overall, the invasion of A. philoxeroides exhibited region-specific effects on soil EMF (Figure 1b). EMF increased significantly in Xiaoxian and Hefei after invasion, whereas changes in Tongling, Anqing, and Huangshan were not significant. Trade-off analysis of C-, N-, and P-related functions revealed that carbon- and nitrogen-associated multifunctionality (C_Mfa and N_Mfa) dominated overall EMF, contributing 38.23 and 37.84%, respectively, while phosphorus-associated multifunctionality (P_Mfa) accounted for 23.93% (Figure 1c). Most samples were distributed along a balanced gradient between C_Mfa and N_Mfa. Invasion did not act uniformly but instead rebalanced functional contributions along the C–N–P axes. In Xiaoxian and Anqing, C_Mfa shifted toward P_Mfa; in Hefei, it shifted toward N_Mfa; in Huangshan, N_Mfa shifted toward P_Mfa; and in Tongling, invaded and uninvaded samples remained closely clustered, reflecting a comparatively balanced functional structure.

3.2 Responses of microbial community and microbial taxa to A. philoxeroides invasion

Community resistance analysis showed that bacterial communities exhibited significantly higher resistance than fungal communities across all regions (Figure 2a). For bacteria, community resistance differed significantly between Xiaoxian, Hefei, and Tongling regions with Anqing and Huangshan regions, showing a decreasing trend with increasing latitude. For fungal communities, the highest resistance was observed in Huangshan, which was significantly greater than in Xiaoxian, Tongling, and Anqing, but not significantly different from Hefei.

Figure 2
(a) Boxplots and scatter plots showing community resistance of bacteria and fungi with respect to different conditions and their core relationships. (b) Network diagrams for bacteria before and after invasion, showing positive and negative interactions among various taxa. (c) Network diagrams for fungi before and after invasion, illustrating interaction shifts. Graphs show changes in nodes, edges, and interaction ratios. Color-coded legends indicate different taxa and interaction types.

Figure 2. Changes in bacterial and fungal communities. (a) The difference in resistance between soil bacterial and fungal communities; The regression analysis between bacterial community resistance and EMF; (b) The interaction effect of species in the bacterial community during A. philoxeroides invasion, in which the nodes are labeled with different colors at the gate level, and the information of the corresponding phylum is labeled in the lower left corner of the network diagram. The pink connection line indicates a positive interaction, and the blue connection line indicates a negative interaction. Nodes represent the number of nodes, Edges represent the number of connection lines, Positive represents the positive interaction rate, and Negative represents the negative interaction rate, the same below; (c) The interaction effect of species in the fungal community during A. philoxeroides invasion, in which the nodes are labeled with different colors at the gate level, and the information of the corresponding phylum is labeled in the lower left corner of the network diagram. And the pink connection line indicates positive interaction and the green connection line indicates negative interaction.

In the co-occurrence networks constructed from the 500 most abundant OTUs, the bacterial networks included phyla such as Actinobacteriota, Proteobacteria, Acidobacteriota, Chloroflexi, and Firmicutes, and the fungal networks included Ascomycota, Basidiomycota, Mortierellomycota, and unclassified fungi (Figures 2b,c). For bacteria, the pre-invasion network contained 462 OTUs and 6,250 edges (57.02% positive, 42.98% negative), which decreased to 431 OTUs and 5,698 edges (56.84% positive, 43.16% negative) after invasion. For fungi, the pre-invasion network contained 402 OTUs and 3,120 edges (95.38% positive, 4.62% negative), which decreased to 364 OTUs and 1,382 edges (70.26% positive, 29.74% negative) after invasion.

Mantel tests indicated that core bacterial taxa were significantly correlated with soil available phosphorus (AP), total nitrogen (TN), EMF, and N_Mfa (p < 0.001), and also related to total phosphorus (TP), phosphatase activity, nitrogen mineralization rate (NMR), and aboveground biomass (AGB) (p < 0.05) (Figure 3b). Keystone bacterial taxa were strongly correlated with ammonium nitrogen (NH₄+–N), microbial biomass nitrogen (MBN), phosphatase activity, and soil β-glucosidase activity. In the fungal community, core microbial taxa were strongly associated with soil organic carbon (SOC) and significantly correlated with microbial biomass phosphorus (MBP), NMR, and P_Mfa. Keystone fungal taxa were significantly related to TP, TN, and soil pH (Figure 3c).

Figure 3
(a) Four scatterplots show relationships between the bacterial core, fungal core, bacterial community resistance, fungal community resistance, and ecosystem multifunctionality. Trend lines with confidence intervals indicate correlations.(b) Heatmap for bacteria displays correlations between various soil and ecosystem variables, with color intensity indicating correlation strength. Arrows connect core and key variables.(c) Heatmap for fungi mirrors the bacterial heatmap, illustrating similar variable correlations. Arrows indicate connections between core and key variables.

Figure 3. Linking bacterial and fungal core, key taxa and ecosystem multifunctionality. (a) The regression analysis between bacterial and fungi core, key taxa and environmental factors; (b) The correlation between bacterial core taxa and community resistance and EMF. Line width corresponds to the partial Mantel’s r statistic. Pairwise comparisons of environmental factors are also shown, with a color gradient denoting Pearson’s correlation coefficient. Asterisks indicate the statistical significance (***p < 0.001; **p < 0.01; and *p < 0.05), the same below; (c) The correlation between fungal core taxa and community resistance and EMF.

3.3 Community resistance mediated by microbial taxa drives changes in soil EMF

Furthermore, both bacterial and fungal community resistance were positively associated with EMF (Figure 2a), indicating that higher microbial resistance was associated with enhanced EMF. Significant positive correlations were observed between selected core microbial taxa (both bacterial and fungal taxa) and community resistance (Figure 3a). Structural equation modeling (SEM) further clarified these relationships (Figure 4, Structural equation modeling overall fit indices: CFl = 0.96, TLI = 0.90, RMSEA = 0.08, SRMR = 0.02, x2/df = 1.27). Invasion of A. philoxeroides and changes in soil properties significantly affected the abundance of core bacterial taxa. Soil water content exerted a strong negative effect on these taxa (path coefficient = −0.50, p < 0.001). The core bacterial taxa were positively associated with the resistance of both bacterial and fungal communities. Notably, bacterial community resistance had a direct and significant positive effect on soil EMF (path coefficient = 1.05, p < 0.05), indicating that the stabilization of microbial communities played a substantial role in maintaining EMF.

Figure 4
Diagram (a) shows a structural path model with multiple factors like Invasion, PH, SWC, Bcore, Fcore, Brs, Frs, and EMF connected by arrows indicating relationships and standardized coefficients. Bar graph (b) displays standardized effect sizes for variables Invasion, SWC, PH, Fcore, Bcore, Brs, and Frs in three colors: green for total, yellow for direct, and orange for indirect effects.

Figure 4. (a) Structural equation modeling (SEM) assessing the effects of Invasion, PH, SWC (soil water content), B and Fcore (bacterial and fungal core taxa), B and Frs (bacterial and fungal community resistance) in EMF (ecosystem multifunctionality). Numbers adjacent to arrows show standardized path coefficients. Red and blue lines indicate positive and negative relationships. ***p < 0.001; **p < 0.01; and *p < 0.05. (b) Standardized total, direct and indirect effects based on the SEM for the concentrations of Invasion, PH, SWC, B and Fcore, B and Frs. The green column is the total effect, the yellow column is the direct effect, and the orange column is the indirect effect.

4 Discussion

4.1 Impact of A. philoxeroides invasion on soil EMF

Soil EMF is a comprehensive indicator of the soil’s capacity to perform multiple ecosystem services and functions simultaneously, supporting nutrient cycling, carbon sequestration, water regulation, biodiversity maintenance, and pollutant degradation (Delgado-Baquerizo et al., 2016; Liu et al., 2023). Our results showed that the invasion of A. philoxeroides significantly enhanced soil EMF in both Xiaoxian and Hefei sites (Figure 1b), with nitrogen and phosphorus cycling being particularly affected. Multifunctionality trade-off analyses for carbon, nitrogen, and phosphorus revealed that in Xiaoxian, the invasion increased the proportion of P_Mfa, whereas in Hefei, it enhanced N_Mfa, both contributing to elevated soil EMF (Figure 1c). However, unlike the general consensus that plant invasions often suppress soil functions (Negesse et al., 2025), A. philoxeroides invasion appears to promote soil EMF. This outcome likely reflects the combined influence of plant traits, soil properties, and environmental conditions (Hu et al., 2024). In Xiaoxian, lower MAP, sandy soil and A. philoxeroides’s complex root network create a loose soil texture (Li et al., 2009). Improved aeration accelerates organic matter decomposition (Paul, 2016), releasing phosphorus bound in organic matter (Feng et al., 2022; Yang et al., 2022), which enhances phosphorus turnover and elevates EMF. In Hefei, the strong stress tolerance, rapid proliferation, and extensive root system of A. philoxeroides facilitate root secretions that provide carbon sources for rhizosphere nitrogen-fixing bacteria (Li et al., 2023; Yang et al., 2020). This symbiotic interaction promotes nitrogen cycling, further increasing EMF. These region-specific pathways suggest that elevated soil EMF is regulated by differential nutrient turnover processes across sites. Given that soil microorganisms are central to nutrient cycling and biogeochemical regulation, invasion-induced microbial restructuring likely mediated these effects (Long et al., 2025; Xu et al., 2025). While invasive plants often reduce microbial diversity and network stability (Chen et al., 2019; Yin et al., 2020), the core microbial taxa associated with A. philoxeroides are robust to such disturbances and can functionally reconfigure soil nutrient dynamics. Taken together, our findings highlight that the enhancement of soil EMF by A. philoxeroides results from the interplay between altered soil properties and microbial community restructuring, supporting our first hypothesis.

4.2 Effects of A. philoxeroides invasion on core-taxa-driven microbial resistance

Our results showed that bacterial communities exhibited stronger resistance than fungal communities across regions (Figure 2a), consistent with previous findings (Delgado-Baquerizo et al., 2016; Jiao et al., 2022). This difference likely reflects contrasting ecological strategies: bacteria, with denser interaction networks and greater interspecies connectivity (Figure 2b), form cooperative and functionally versatile assemblages (Wang et al., 2022; Zhang et al., 2018). These relationships, coupled with high metabolic flexibility, enable rapid adjustment to environmental fluctuations (Wang et al., 2022), enhancing resistance to invasion-driven disturbances. In contrast, fungi often occupy narrower niches and rely on more specialized resource-use strategies, which slow their responses to change and reduce short-term contributions to community resistance (Chen et al., 2020; Jing et al., 2015).

The capacity of microbial communities to resist invasion is closely linked to the presence and function of core microbial taxa, which are consistently detected across habitats and disproportionately shape community structure and ecosystem processes (Jacquiod et al., 2020; Sun et al., 2024). In our study, these taxa were positively correlated with microbial resistance (Figure 3a), supporting their role as stabilizing agents within microbial networks (Wertz et al., 2007; Zhu et al., 2016). Mechanistically, their stabilizing effect likely arises from three complementary traits: stable diversity, which maintains functional groups across gradients; functional redundancy, which buffers against loss of individual taxa; and niche complementarity, which promotes efficient resource partitioning and minimizes competition. Together, these mechanisms explain why core microbial taxa are pivotal in maintaining microbial community stability under invasion pressure, supporting our second hypothesis.

4.3 Core microbial taxa mediating microbial resistance as a driver of soil EMF

Our results revealed a positive relationship between bacterial community resistance and soil EMF (Figure 2a), consistent with previous findings that microbial resistance benefited ecosystem functions. This relationship stems from the capacity of microbial communities to maintain compositional and functional stability, thereby supporting ecosystem processes. In this study, bacterial communities exhibited broad functional stability owing to their complex interaction networks, particularly within core microbial taxa such as Proteobacteria and Actinobacteria. Proteobacteria, one of the most functionally diverse bacterial phyla, contribute to carbon fixation, organic matter decomposition, and plant–microbe symbioses, thereby enhancing nutrient acquisition (Peter et al., 2011; Wongkiew et al., 2022). Actinobacteria, in turn, degrade recalcitrant compounds, suppress soil pathogens, and facilitate nutrient turnover (Zhang et al., 2020; Rousk and Bååth, 2011). In parallel, fungi also play a crucial role in long-term ecosystem functioning by decomposing lignin and other recalcitrant organic compounds, stabilizing soil organic carbon through oregano–mineral associations, and enhancing nutrient cycling via mycorrhizal symbioses. Such contributions highlight their importance in sustaining nutrient turnover and carbon storage over extended timescales. In our study, their stable functional traits ultimately drove soil EMF under invasion conditions. These findings support our third hypothesis by demonstrating that community resistance mediated by core bacterial taxa plays a central role in sustaining soil EMF under invasion pressure. Taken together, they highlight the importance of core bacterial groups as ecological stabilizers that safeguard multifunctionality in invaded ecosystems.

The structural equation modeling further clarified the cascading pathway in which A. philoxeroides invasion and associated changes in soil properties significantly influenced soil EMF through community resistance mediated by core bacterial taxa (Figure 4). Under invasion pressure, core bacterial groups may strengthen resistance via two complementary mechanisms. First, their ability to utilize diverse carbon and nitrogen sources sustains metabolic activity under fluctuating resources, buffering community functioning against nutrient imbalances. Second, their high connectivity within microbial symbiotic networks facilitates resource redistribution, stabilizing both bacterial and fungal communities. This mechanism aligns with previous findings that cooperative interactions among dominant bacterial groups mitigate disturbance impacts on microbial structure and function (Wang et al., 2021). Moreover, core bacterial taxa were more responsive to environmental change than core fungal taxa, likely because fungi often exhibit narrower niches and slower turnover rates, limiting their capacity to rapidly compensate for functional losses (Jing et al., 2015; Xue et al., 2023). Ultimately, bacterial community resistance exerted a direct and significant positive effect on soil EMF, highlighting the pivotal role of core bacterial groups as ecological stabilizers in maintaining EMF under invasion. Our results suggest that core microbial taxa can serve as effective bio indicators of soil stability under invasion, offering a practical tool for monitoring belowground resilience. Their consistent contribution to soil EMF also highlights them as potential targets for microbial interventions, such as inoculation or soil amendments, to enhance resistance and accelerate ecosystem recovery. Integrating core microbial taxa into land management frameworks may therefore improve the prediction and restoration of soil functions in invaded ecosystems.

5 Conclusion

The invasion of A. philoxeroides significantly enhanced soil EMF at specific sites, such as Xiaoxian and Hefei, mainly driven by N_Mfa and P_Mfa, respectively. Integrating microbial resistance analysis, co-occurrence networks, and structural equation modeling, we found that bacterial communities exhibited stronger resistance than fungi, and this resistance was largely mediated by core bacterial taxa, including Proteobacteria and Actinobacteria. These taxa function as ecological stabilizers buffering invasion-induced disturbances and sustaining soil EMF. Notably, bacterial resistance exerted a direct and significant positive effect on EMF, highlighting its critical role in regulating EMF under invasion pressure. These findings, through cross-regional comparisons and comparative analyses of bacteria and fungi, have provided new mechanistic insights into the microbial basis of ecosystems’ responses to biological invasions, and it provides information for the control of invasive species, microbial inoculation or soil remediation measures, and at the same time offers a basis for predicting and managing the situation of soil electromagnetic fields in the context of future global changes. This study is based on a single sampling campaign, and future temporal monitoring would help capture seasonal dynamics and further refine predictions of soil EMF under invasion.

Data availability statement

The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive (Genomics, Proteomics & Bioinformatics 2025) in National Genomics Data Center (Nucleic Acids Res 2025), China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences (GSA: CRA032627 and GSA: CRA031326) that are publicly accessible at https://ngdc.cncb.ac.cn/gsa.

Author contributions

XL: Writing – original draft, Writing – review & editing. XZ: Writing – original draft. SL: Data curation, Writing – review & editing. YW: Software, Writing – review & editing. HY: Methodology, Writing – review & editing. WQ: Formal analysis, Writing – review & editing. JZ: Funding acquisition, Project administration, Software, Writing – review & editing.

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 National Natural Science Foundation of China (41771557 and 42307389).

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 author(s) declare that no Gen AI was 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.

Supplementary material

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

References

Allison, S. D., and Martiny, J. B. H. (2008). Resistance, resilience, and redundancy in microbial communities. Proc. Natl. Acad. Sci. USA 105, 11512–11519. doi: 10.1073/pnas.0801925105

PubMed Abstract | Crossref Full Text | Google Scholar

Bremner, J. M., and Mulvaney, C. S. (1982). Nitrogen-Total. In: Methods of soil analysis. Part 2. Chemical and microbiological properties, (eds.) A. L. Page, R. H. Miller, and D. R. Keeney, (, Madison, Wisconsin: American Society of Agronomy, Soil Science Society of America), 595–624.

Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). Qiime allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, Q.-L., Ding, J., Zhu, D., Hu, H.-W., Delgado-Baquerizo, M., Ma, Y.-B., et al. (2020). Rare microbial taxa as the major drivers of ecosystem multifunctionality in long-term fertilized soils. Soil Biol. Biochem. 141:107686. doi: 10.1016/j.soilbio.2019.107686

Crossref Full Text | Google Scholar

Chen, L., Fang, K., Zhou, J., Yang, Z.-P., Dong, X.-F., Dai, G.-H., et al. (2019). Enrichment of soil rare bacteria in root by an invasive plant Ageratina adenophora. Sci. Total Environ. 683, 202–209. doi: 10.1016/j.scitotenv.2019.05.220

PubMed Abstract | Crossref Full Text | Google Scholar

Coban, O., De Deyn, G. B., and van der Ploeg, M. (2022). Soil microbiota as game-changers in restoration of degraded lands. Science 375:abe0725. doi: 10.1126/science.abe0725

PubMed Abstract | Crossref Full Text | Google Scholar

Delgado-Baquerizo, M., Maestre, F. T., Reich, P. B., Jeffries, T. C., Gaitan, J. J., Encinar, D., et al. (2016). Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat. Commun. 7:10541. doi: 10.1038/ncomms10541

PubMed Abstract | Crossref Full Text | Google Scholar

DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05

PubMed Abstract | 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

Fan, S., Yu, H., Dong, X., Wang, L., Chen, X., Yu, D., et al. (2016). Invasive plant Alternanthera philoxeroides suffers more severe herbivory pressure than native competitors in recipient communities. Sci. Rep. 6:36542. doi: 10.1038/srep36542

Crossref Full Text | Google Scholar

Feng, J., He, K., Zhang, Q., Han, M., and Zhu, B. (2022). Changes in plant inputs alter soil carbon and microbial communities in forest ecosystems. Glob. Chang. Biol. 28, 3426–3440. doi: 10.1111/gcb.16107

PubMed Abstract | Crossref Full Text | Google Scholar

Hao, D., Song, S., Mu, J., Hao, D. C., Song, S. M., Hu, W. L., et al. (2016). Unearthing microbial diversity of Taxus rhizosphere via MiSeq high-throughput amplicon sequencing and isolate characterization. Sci. Rep. 6:22006. doi: 10.1038/srep22006

PubMed Abstract | Crossref Full Text | Google Scholar

Harris, D., Voroney, R. P., and Paul, E. A. (1997). Measurement of microbial biomass N:C by chloroform fumigation-incubation. Can. J. Soil Sci. 77, 507–514. doi: 10.4141/s96-064

Crossref Full Text | Google Scholar

He, Y., Li, J., Siemann, E., Li, B., Xu, Y., and Wang, Y. (2025). Plant invasion increases soil microbial biomass carbon: meta-analysis and empirical tests. Glob. Change Biol. 31:e70109. doi: 10.1111/gcb.70109

PubMed Abstract | Crossref Full Text | Google Scholar

Hu, Z., Delgado-Baquerizo, M., Fanin, N., Chen, X., Zhou, Y., Du, G., et al. (2024). Nutrient-induced acidification modulates soil biodiversity-function relationships. Nat. Commun. 15:2858. doi: 10.1038/s41467-024-47323-3

PubMed Abstract | Crossref Full Text | Google Scholar

Hu, T., Wang, X., Zhen, L., Gu, J., Zhang, K., Wang, Q., et al. (2019). Effects of inoculation with lignocellulose-degrading microorganisms on antibiotic resistance genes and the bacterial community during co-composting of swine manure with spent mushroom substrate. Environ. Pollut. 252, 110–118. doi: 10.1016/j.envpol.2019.05.078

PubMed Abstract | Crossref Full Text | Google Scholar

Jacquiod, S., Puga-Freitas, R., Spor, A., Mounier, A., Monard, C., Mougel, C., et al. (2020). A core microbiota of the plant-earthworm interaction conserved across soils. Soil Biol. Biochem. 144:107754. doi: 10.1016/j.soilbio.2020.107754

Crossref Full Text | Google Scholar

Ji, L., Sheng, S., Shen, F., Yang, L., Wen, S., He, G., et al. (2024). Stochastic processes dominated the soil bacterial community assemblages along an altitudinal gradient in boreal forests. Catena 237:107816. doi: 10.1016/j.catena.2024.107816

Crossref Full Text | Google Scholar

Jiao, S., Peng, Z., Qi, J., Gao, J., and Wei, G. (2021). Linking bacterial-fungal relationships to microbial diversity and soil nutrient cycling. mSystems 6:e01052-20. doi: 10.1128/msystems.01052-20

PubMed Abstract | Crossref Full Text | Google Scholar

Jiao, S., Qi, J., Jin, C., Liu, Y., Wang, Y., Pan, H., et al. (2022). Core phylotypes enhance the resistance of soil microbiome to environmental changes to maintain multifunctionality in agricultural ecosystems. Glob. Chang. Biol. 28, 6653–6664. doi: 10.1111/gcb.16387

PubMed Abstract | Crossref Full Text | Google Scholar

Jing, X., Sanders, N. J., Shi, Y., Chu, H., Classen, A. T., Zhao, K., et al. (2015). The links between ecosystem multifunctionality and above- and belowground biodiversity are mediated by climate. Nat. Commun. 6:8159. doi: 10.1038/ncomms9159

PubMed Abstract | Crossref Full Text | Google Scholar

Lai, Y., Yan, Z., Xiao, F., Zhou, Y., Liu, Y., Li, Y., et al. (2025). Temporal changes in plant–soil feedbacks between the invasive Phytolacca americana and congeneric native and non-invasive alien plants. Oikos 2025:e10772. doi: 10.1111/oik.10772

Crossref Full Text | Google Scholar

Li, X., Huang, C., Pang, J., and He, Z. (2009). Dust source of the Holocene loess-soil and pedogenic environmental changes in the upper Huaihe River. J. Geogr. Sci. 19, 107–117. doi: 10.1007/s11442-009-0107-z

Crossref Full Text | Google Scholar

Li, X., Zhang, Y., Kong, F.-L., Naz, M., Zhou, J.-Y., Qi, S.-S., et al. (2023). Invasive plant Alternanthera philoxeroides benefits more competition advantage from rhizosphere bacteria regardless of the host source. Plants 12:2085. doi: 10.3390/plants12112085

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Z., Zhuang, J., Zheng, K., and Luo, C. (2023). Differential response of the soil nutrients, soil bacterial community structure and metabolic functions to different risk areas in lead-zine tailings. Front. Microbiol. 14:1131770. doi: 10.3389/fmicb.2023.1131770

PubMed Abstract | Crossref Full Text | Google Scholar

Long, X., Li, J., Liao, X., Wang, J., Zhang, W., Wang, K., et al. (2025). Stable soil biota network enhances soil multifunctionality in agroecosystems. Glob. Change Biol. 31:e70041. doi: 10.1111/gcb.70041

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, G., Xue, C., Jiang, Q., Xiao, Y., Zhang, F., Guo, S., et al. (2020). Soil Carbon, Nitrogen, and Phosphorus cycling microbial populations and their resistance to global change depend on soil C:N:P stoichiometry. mSystems 5:e00162-20. doi: 10.1128/mSystems.00162-20

PubMed Abstract | Crossref Full Text | Google Scholar

Maestre, F. T., Castillo-Monroy, A. P., Bowker, M. A., and Ochoa-Hueso, R. (2012). Species richness effects on ecosystem multifunctionality depend on evenness, composition and spatial pattern. J. Ecol. 100, 317–330. doi: 10.1111/j.1365-2745.2011.01918.x

Crossref Full Text | Google Scholar

Morecroft, M. D., Marrs, R. H., and Woodward, F. I. (1992). Altitudinal and seasonal trends in soil nitrogen mineralization rate in the Scottish highlands. J. Ecol. 80, 49–56. doi: 10.2307/2261062

Crossref Full Text | Google Scholar

Negesse, Z., Pan, K., Guadie, A., Justine, M. F., Azene, B., Pandey, B., et al. (2025). Plant invasions alter soil biota and microbial activities: a global meta-analysis. Plant Soil 513, 1031–1050. doi: 10.1007/s11104-025-07227-7

Crossref Full Text | Google Scholar

Nguyen, T. H., Jones, S. D., Soto-Berelov, M., Haywood, A., and Hislop, S. (2020). Monitoring aboveground forest biomass dynamics over three decades using Landsat time-series and single-date inventory data. Int. J. Appl. Earth Obs. Geoinf. 84:101952. doi: 10.1016/j.jag.2019.101952

Crossref Full Text | Google Scholar

Nóbrega, G. N., Ferreira, T. O., Artur, A. G., Mendonça, E. S., Leão, R. A. O., and Teixeira, A. S. (2015). Evaluation of methods for quantifying organic carbon in mangrove soils from semi-arid region. J. Soils Sediments 15, 282–291. doi: 10.1007/s11368-014-1019-9

Crossref Full Text | Google Scholar

Oberson, A., Friesen, D. K., Morel, C., and Tiessen, H. (1997). Determination of phosphorus released by chloroform fumigation from microbial biomass in high P sorbing tropical soils. Soil Biol. Biochem. 29, 1579–1583. doi: 10.1016/S0038-0717(97)00049-7

Crossref Full Text | Google Scholar

Olsen, S. R. (1954). Estimation of available phosphorus in soils by extraction with sodium bicarbonate. Washington DC: US Government Printing Office.

Google Scholar

Paul, E. A. (2016). The nature and dynamics of soil organic matter: plant inputs, microbial transformations, and organic matter stabilization. Soil Biol. Biochem. 98, 109–126. doi: 10.1016/j.soilbio.2016.04.001

Crossref Full Text | Google Scholar

Peter, H., Ylla, I., Gudasz, C., Romaní, A. M., Sabater, S., and Tranvik, L. J. (2011). Multifunctionality and diversity in bacterial biofilms. PLoS One 6:e23225. doi: 10.1371/journal.pone.0023225

PubMed Abstract | Crossref Full Text | Google Scholar

Rousk, J., and Bååth, E. (2011). Growth of saprotrophic fungi and bacteria in soil. FEMS Microbiol. Ecol. 78, 17–30. doi: 10.1111/j.1574-6941.2011.01106.x

PubMed Abstract | Crossref Full Text | Google Scholar

Schade, J. D., Kyle, M., Hobbie, S. E., Fagan, W. F., and Elser, J. J. (2003). Stoichiometric tracking of soil nutrients by a desert insect herbivore. Ecol. Lett. 6, 96–101. doi: 10.1046/j.1461-0248.2003.00409.x

Crossref Full Text | Google Scholar

Scheller, R. M., and Parajuli, R. (2018). Forest management for climate change in New England and the Klamath ecoregions: motivations, practices, and barriers. Forests 9:626. doi: 10.3390/f9100626

Crossref Full Text | Google Scholar

Shade, A., and Handelsman, J. (2012). Beyond the Venn diagram: the hunt for a core microbiome. Environ. Microbiol. 14, 4–12. doi: 10.1111/j.1462-2920.2011.02585.x

PubMed Abstract | Crossref Full Text | Google Scholar

Shade, A., Peter, H., Allison, S. D., Baho, D., Berga, M., Buergmann, H., et al. (2012). Fundamentals of microbial community resistance and resilience 3:417. doi: 10.3389/fmicb.2012.00417

Crossref Full Text | Google Scholar

Sun, L., Li, G., Zhao, J., Zhang, T., Liu, J., and Zhang, J. (2024). Core microbiota drive multi-functionality of the soil microbiome in the Cinnamomum camphora coppice planting. BMC Microbiol. 24:18. doi: 10.1186/s12866-023-03170-8

PubMed Abstract | Crossref Full Text | Google Scholar

Sun, F., Zeng, L., Cai, M., Chauvat, M., Forey, E., Tariq, A., et al. (2022). An invasive and native plant differ in their effects on the soil food-web and plant-soil phosphorus cycle. Geoderma 410:115672. doi: 10.1016/j.geoderma.2021.115672

Crossref Full Text | Google Scholar

Vietorisz, C. R., Nash, J. A., Siggers, J. A., Leander, E. J., Bock, B. M., Camuy-Vélez, L. A., et al. (2025). Pine-fungal co-invasion alters whole-ecosystem properties of a native eucalypt forest. New Phytol. 247, 2342–2356. doi: 10.1111/nph.70363

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Q., Chen, Z., Zhao, J., Ma, J., Yu, Q., Zou, P., et al. (2022). Fate of heavy metals and bacterial community composition following biogas slurry application in a single rice cropping system. J. Soils Sediments 22, 968–981. doi: 10.1007/s11368-021-03117-4

Crossref Full Text | Google Scholar

Wang, J., Wang, X., Liu, G., Zhang, C., and Wang, G. (2021). Bacterial richness is negatively related to potential soil multifunctionality in a degraded alpine meadow. Ecol. Indic. 121:106996. doi: 10.1016/j.ecolind.2020.106996

Crossref Full Text | Google Scholar

Wang, H., Wang, S., Wang, R., Zhang, Y., Wang, X., and Li, J. (2019). Direct and indirect linkages between soil aggregates and soil bacterial communities under tillage methods. Geoderma 354:113879. doi: 10.1016/j.geoderma.2019.113879

Crossref Full Text | Google Scholar

Wertz, S., Degrange, V., Prosser, J. I., Poly, F., Commeaux, C., Guillaumaud, N., et al. (2007). Decline of soil microbial diversity does not influence the resistance and resilience of key soil microbial functional groups following a model disturbance. Environ. Microbiol. 9, 2211–2219. doi: 10.1111/j.1462-2920.2007.01335.x

PubMed Abstract | Crossref Full Text | Google Scholar

Wongkiew, S., Chaikaew, P., Takrattanasaran, N., and Khamkajorn, T. (2022). Evaluation of nutrient characteristics and bacterial community in agricultural soil groups for sustainable land management. Sci. Rep. 12:7368. doi: 10.1038/s41598-022-09818-1

PubMed Abstract | Crossref Full Text | Google Scholar

Xiao, H. F., Feng, Y. L., Schaefer, D. A., and Yang, X. D. (2014). Soil fungi rather than bacteria were modified by invasive plants, and that benefited invasive plant growth. Plant Soil 378, 253–264. doi: 10.1007/s11104-014-2040-x

Crossref Full Text | Google Scholar

Xu, D., Yu, X., Zhang, Y., Liu, Y., Chen, C., Li, L., et al. (2025). Effect of compost as a soil amendment on the structure and function of fungal diversity in saline-alkali soil. Curr. Res. Microb. Sci. 9:100405. doi: 10.1016/j.crmicr.2025.100405

PubMed Abstract | Crossref Full Text | Google Scholar

Xue, H., Wang, W., Wang, X., Zhang, B., Shi, F., Kurakov, A. V., et al. (2023). The potential of mycorrhizal fungi to increase terrestrial ecosystem carbon sink: a review. Eurasian Soil Sci. 56, 1724–1738. doi: 10.1134/S1064229323601178

Crossref Full Text | Google Scholar

Yang, X., Kong, Y., Guo, E., Chen, X., and Li, L. (2022). Organic acid regulation of inorganic phosphorus release from Mollisols with different organic matter contents. Soil Use Manag. 38, 576–583. doi: 10.1111/sum.12710

Crossref Full Text | Google Scholar

Yang, R., Li, K., Fang, J., Quan, Q., Zhang, C., and Liu, J. (2020). The invasion of Alternanthera philoxeroides increased soil organic carbon in a river and a constructed wetland with different mechanisms. Front. Ecol. Evol. 8:574528. doi: 10.3389/fevo.2020.574528

Crossref Full Text | Google Scholar

Yin, L., Liu, B., Wang, H., Zhang, Y., Wang, S., Jiang, F., et al. (2020). The rhizosphere microbiome of Mikania micrantha provides insight into adaptation and invasion. Front. Microbiol. 11:1462. doi: 10.3389/fmicb.2020.01462

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, H.-Y., Goncalves, P., Copeland, E., Qi, S.-S., Dai, Z.-C., Li, G.-L., et al. (2020). Invasion by the weed Conyza canadensis alters soil nutrient supply and shifts microbiota structure. Soil Biol. Biochem. 143:107739. doi: 10.1016/j.soilbio.2020.107739

Crossref Full Text | Google Scholar

Zhang, J., Li, M., Jia, K., Zheng, G., and Long, X.-E. (2018). Seasonal variation rather than stand age determines bacterial diversity in the rhizosphere of wolfberry (Lycium barbarum L.) associated with soil degradation. J. Soils Sediments 18, 1518–1529. doi: 10.1007/s11368-017-1854-6

Crossref Full Text | Google Scholar

Zhao, P., Liu, B., Zhao, H., Lei, Z., and Zhou, T. (2023). Significant changes in soil microbial community structure and metabolic function after Mikania micrantha invasion. Sci. Rep. 13:1141. doi: 10.1038/s41598-023-27851-6

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, M., Yang, Y., Zhang, H., Li, Q., Zhao, X., Guo, X., et al. (2024). Asymmetric succession in soil microbial communities enhances the competitive advantage of invasive alien plants. Microbiome 12:265. doi: 10.1186/s40168-024-01989-5

PubMed Abstract | Crossref Full Text | Google Scholar

Zheng, B., Hui, N., Jumpponen, A., Lu, C., Pouyat, R., Szlavecz, K., et al. (2025). Urbanization leads to asynchronous homogenization of soil microbial communities across biomes. Environ. Sci. Ecotechnol. 25:100547. doi: 10.1016/j.ese.2025.100547

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, Z., Wang, C., and Luo, Y. (2020). Meta-analysis of the impacts of global change factors on soil microbial diversity and functionality. Nat. Commun. 11:3072. doi: 10.1038/s41467-020-16881-7

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, H., Fazliddin, K., Li, Q., Wang, C., Chen, P., Yang, J., et al. (2025). Contrasting adaptations of soil prokaryotes and arbuscular mycorrhizal fungi in saline wildland and non-saline farmland. Fundam. Res. doi: 10.1016/j.fmre.2025.02.009

Crossref Full Text | Google Scholar

Zhu, D. H., Wang, P., Zhang, W. Z., Yuan, Y., Li, B., and Wang, J. (2016). Correction: sampling and complementarity effects of plant diversity on resource use increases the invasion resistance of communities. PLoS One 11:e0150128. doi: 10.1371/journal.pone.0150128

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: A. philoxeroides, microorganisms, soil ecosystem multifunctionality, core microbial taxa, community resistance

Citation: Lin X, Zhang X, Li S, Wang Y, Yu H, Qin W and Zhang J (2025) Soil core microbial taxa maintain community resistance to drive soil ecosystem multifunctionality under Alternanthera philoxeroides invasion. Front. Microbiol. 16:1707273. doi: 10.3389/fmicb.2025.1707273

Received: 17 September 2025; Accepted: 20 October 2025;
Published: 06 November 2025.

Edited by:

Yanfu Bai, Sichuan Agricultural University, China

Reviewed by:

Huitong Li, Yan’an University, China
Xu Yehong, Jiangsu Academy of Agricultural Sciences (JAAS), China

Copyright © 2025 Lin, Zhang, Li, Wang, Yu, Qin and Zhang. 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: Wenxiu Qin, d3hxaW5AYWhhdS5lZHUuY24=; Jiaoyang Zhang, amlhb3lhbmc5NjY2QGFoYXUuZWR1LmNu

These authors have contributed equally to this work

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.