- Research
- Open access
- Published:
Asymmetric succession in soil microbial communities enhances the competitive advantage of invasive alien plants
Microbiome volume 12, Article number: 265 (2024)
Abstract
Background
Biological invasions pose an escalating threat to native ecosystems. The accumulation of invasive alien plants worldwide is not saturated yet, underscoring the persistent and growing impact of invasions. Soil microorganisms play a key role in the process of alien plant invasion. However, the temporal dynamics of microbial communities has rarely been determined during the invasion owing to the dearth of long-term, in situ experimental systems.
Results
Here, we examined the temporal succession of soil microbial communities 8 years after experiment setup in a common garden. Bacterial communities displayed divergent temporal succession, with invasive plants exhibiting higher turnover rates. Invasive alien plants reduced stochasticity in bacterial communities, likely acting as an environmental filter on community assembly. Plant growth-promoting microbes underwent higher succession rates in invasive alien plants compared to native plants, suggesting that invasive alien plants may possess a distinct advantage in fostering a favorable microbiota for their own growth and establishment. In sharp contrast, native plants selectively increased succession rates of specific plant pathogens. Furthermore, the microbial co-occurrence network was more complex in invasive plants, suggesting that invasive plants foster intricate relationships among microbial communities.
Conclusions
Therefore, the asymmetric succession in soil microbial communities enables invasive plants recruit beneficial microbiota from the surrounding soil. These results deepen our understanding of the mechanism underlying plant invasion and provide novel insights into predicting the ecological consequences resulting from widespread plant invasion. This knowledge can be incorporated into management strategies to address the evolving challenges posed by invasive plants.
Video Abstract
Introduction
The incidence of biological invasion intensifies complex global environmental challenges, including biodiversity loss, biota homogenization, and habitat fragmentation [1], which has become a pervasive aspect of global change. Invasive alien plants, characterized by their rapid growth and substantial biomass accumulation, often colonize and expand swiftly, leading to reduced plant diversity and the formation of monocultural communities in new habitats [2]. Invasive alien plants can disrupt ecosystem dynamics by altering soil biogeochemical pools and energy fluxes. These disruptions can have long-term consequences, potentially leading to the degradation of ecosystems and the loss of valuable ecosystem services [3]. Previous studies have primarily focused on deciphering the mechanisms that enable invasive alien plants to establish and spread more successfully than their native counterparts [4,5,6,7].
In recent years, there has been growing recognition of the pivotal role that soil microorganisms associated with plants play in enhancing plant reproduction, fitness, nutrient acquisition, and resistance to environmental stresses [8]. Invasive plants like Mikania micrantha, Ageratina adenophora, and Flaveria bidentis can increase soil carbon and nitrogen availability by enriching microbes involved in these processes [9, 10, 2, 11]. Increased access to limiting nutrients (e.g., C, N, P) can thus stimulate the rapid growth and reproduction of exotic invasive plants [12]. Plant growth-promoting microbes, such as arbuscular mycorrhizal fungi (AMF), have been found to colonize the roots of exotic F. bidentis more extensively than those of co-occurring native plants [13], and AMF inoculation has facilitated the competition of F. bidentis [14]. Furthermore, many plant-soil feedback (PSF) experiments have demonstrated positive effects on invasive plant species in soils conditioned by the invasive plants themselves [15]. For example, Lespedeza cuneata and Prosopis juliflora benefit from a history of invasion, with these positive effects attributed to modifications in soil microbial communities17, 18. In contrast, many native plant species exhibit less growth on soil conditioned by invasive plants, which is largely due to allelochemical effects by invasive plants or accumulation of soil-borne pathogens. However, most of these studies are based on single-time observations due to the dearth of long-term, in situ experimental systems, so we cannot determine whether these phenomena are temporary or vary over time. In addition, without longitudinal data, predicting the future ecological consequences of invasion remains challenging. Therefore, studies on chronosequence-based microbial successions after plant invasion are necessary to address this uncertainty.
Alien plant invasions are temporally dynamic; however, the role of soil microorganisms in driving invasions is often overlooked, especially in terms of temporal components [15]. Owing to the uncertainty in succession trajectories, which can exhibit patterns of convergent, divergent, cyclical or linear, etc. [18], one of the major challenges in understanding biological invasion lies in the difficulty of predicting its ecological consequences. The litter input from invasive alien plants typically supports larger microbial populations [19], triggering a cascade of microbial succession. Additionally, specific root exudates, such as allelochemicals, can attract specific microbes and potentially accelerate the temporal succession of microbial communities. These dynamics introduce a level of unpredictability in forecasting whether microbial communities will exhibit divergent or convergent behaviors, making it challenging to know the ecological outcomes of plant invasions.
Previous studies have shown that alien plant invasion can trigger an increase in relative abundances of soil microbes involved in nutrient cycling [2, 9] or a decrease in pathogenic microbes21. However, it remains unclear whether all functional groups within soil microbial communities respond uniformly to plant invasion. In this context, we hypothesize that (1) alien plant invasion may result in nonuniform succession dynamics among various microbial groups, specifically activating those that promote plant growth over time, thereby enhancing the temporal successions of these microorganisms, and (2) as increased nutrient availability triggers microbial activities in soils alien plants, invasive plants may increase the complexity of microbial networks, resulting in more nodes and keystone taxa.
Materials and methods
Site description and sample collection
In 2012, we initiated a common garden experiment at Langfang Experimental Station of the Chinese Academy of Agricultural Sciences in Hebei province, China (39° 30′ 42″ N, 116° 30′ 42″ N). The experiment included four types of annual plants: invasive Ambrosia artemisiifolia L. and Bidens pilosa L and native Chenopodium serotinum L. and Setaria viridis L. The two invasive plants were selected due to their severe ecological impact and widespread distribution in China, including the experiment site [21, 22]. Two native plants frequently co-occur in areas invaded by the two invasive plants in China and have been used as control natives [13, 23] (Fig. 1a). This setup resulted in two treatments: the native treatment involved monoculture of C. serotinum or S. viridis, and the invasive treatment involved intercropping of A. artemisiifolia and C. serotinum, A. artemisiifolia and S. viridis, B. pilosa and C. serotinum, and B. pilosa and S. viridis. Additionally, bare fallow soil was used as a control. Each type of plantation had three replicate plots, each measuring 3 m in length, 2 m in width, and spaced 1 m apart. This design resulted in a total of 21 plots [(2 monocultural plantations + 4 intercropping plantations + 1 bare fallow soils) × 3 replicates], arranged in a randomized complete block design. Our aim was to investigate the impact on soil microbial communities in native species when invasive species are introduced, a scenario that frequently occurs in nature. Therefore, we focused on comparing mixtures and monocultures of native species, excluding monocultures of invasive plants from our data analysis. A total of 100 seeds (50 seeds of each plant in intercropped plots) were randomly broadcasted by hand in each plot at the beginning of the experiment, with no additional seeding or fertilizer applied in the subsequent years. Plants were removed as needed to maintain the original experimental design, and aboveground biomass was manually removed before the growth season to reduce the shading effect on plant growth.
Effects of alien plant invasion on soil microbial communities. a Illustration of the experiment treatments. The native treatment includes all samples from monocultures of C. serotinum and S. viridis; the invasive treatment includes all intercropping samples from the intercropping of A. artemisiifolia and C. serotinum, A. artemisiifolia and S. viridis, B. pilosa and C. serotinum, and B. pilosa and S. viridis. “Baseline” refers to soil samples collected at the beginning of the experiment in 2012. Then soil samples were taken yearly from 2016 to 2020. b NMDS ordination of bacterial communities based on Bray–Curtis dissimilarity, with the Adonis analysis indicating significant among treatments. c NMDS ordination of fungal communities. d Shannon diversities of bacterial and fungal communities. *P < 0.05, **P < 0.01, ***P < 0.001
Soil samples were collected from three bare fallow plots as baseline samples in October 2012. After allowing 4 years for system stabilization, soil samples were collected from each plot every October from 2016 to 2020, resulting in a total of 108 samples. The 15 samples from bare fallow soils were not involved in downstream analyses. Before collecting soil in each plot, we randomly selected five cores of 20-cm depth and 30-cm width and removed 1 cm of litter from the surface soil. Then we collected soils in cores, avoiding roots. The collected soils were homogenized, sieved through a 20-mesh sieve to remove impurities, and stored in polyethylene (PE) bags at − 80 °C for 16S rRNA and ITS2 gene sequencing. Fresh soils were used to analyze soil properties.
Measurements of environmental variables
For each sampling point, we determined the plant biomass by placing plant tissues per 0.25 m2 into paper bags and drying them at 80 °C for 72 h until a constant weight was achieved. We determined soil geochemical properties by traditional methods [24]. Briefly, we measured soil organic matter content (SOM) by a potassium dichromate oxidation external heating method26. We measured soil pH in a soil suspension (soil:water of 1:2.5) using a pH meter (Mettler-Toledo Instruments, Shanghai, China). We determined total nitrogen (TN) by the Kjeldahl method with sulfuric acid as accelerator27, total phosphorus (TP) by a molybdenum antimony colorimetry method after extraction by sodium carbonate28, and total potassium (TK) by flame photometry after extraction using sodium hydroxide29.
Experiments of microbial communities
We extracted soil genome DNA using the MoBio PowerSoil DNA Isolation Kit (MoBio Laboratories, Carlsbad, CA, USA) and purified the DNA using the MagaBio Soil DNA Purification Kit (Bioer Technology, Hangzhou, Zhejiang province, China). We then assessed DNA quality and concentration using Thermo NanoDrop One (Thermo Fisher Scientific, Waltham, MA, USA) by the absorbance of DNA at 230 nm, 260 nm, and 280 nm.
We amplified the V4 hypervariable region of the 16S rRNA gene for bacteria and archaea using the primer set 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) and the fungal ITS2 gene using the primer set ITS2F (5′-GCATCGATGAAGAACGCAGC-3′) and ITS2R (5′-TCCTCCGCTTATTGATATGC-3′). The 50-µl PCR reaction system included 2 × Premix Taq (Takara Biotechnology, Dalian Co. Ltd., Liaoning Province, China), 1 μl of each primer (10 μM), and 3 μl of DNA template (20 ng/μl) with a thermal cycling condition of an initial denaturation at 94 °C for 5 min, 30 cycles at 94 °C for 30 s, 52 °C for 30 s, 72 °C for 30 s, and a final extension at 72 °C for 10 min. We then purified PCR products using an E.Z.N.A. Gel Extraction Kit (Omega, Bio-tek, Norcross, GA, USA) and sequenced the PCR products on the Illumina MiSeq platform (Illumina Inc., San Diego, CA, USA).
Data analyses
We processed the sequencing data on a Galaxy pipeline as previously described (http://zhoulab5.rccc.ou.edu) [28]. Briefly, we trimmed sequences to 245 ~ 260 bp for the 16S rRNA gene and 250 ~ 350 bp for the ITS gene and conducted a chimera check with the UCHIME method [29]. The remaining sequences were classified into operational taxonomic units (OTUs) with 97% similarity using the UPARSE algorithm [30]. After singleton removal, sequences were rarefied to 67,259 reads for the 16S rRNA V4 region and to 28,850 reads per sample for the fungal ITS gene. We assigned taxonomic annotations of OTUs to representative sequences by a 16S rRNA training set for the 16S rRNA V4 region amplicon [31] and a UNITE database 8.2 version for the fungal ITS34 using RDP Classifier. Then we determined the relative abundance (RA%) of sequences, which were used in downstream analyses as the following:\(\textit{RAij}\%=\text{}\frac{\textit{Sij}}{\sum_\text{j=1}^\text{N}\textit{Sij}}\,\times\,100\), where Sij was the sequence number of the jth OTU in the ith sample.
Microbial time-decay relationships (TDRs) were determined by similar linear regressions between logarithmic β similarities and logarithmic temporal (year) distance as previously described [33] using the ieggr (V1.6) package in R software (R Development Core Team, R Foundation for Statistical Computing, Vienna, Austria) [34, 35]. The equation to evaluate TDRs was as follow:
ln(Ss) = c − ωln(T) + ε
where Ss was the pairwise similarity between pairwise microbial communities within the same treatment, T was the time interval (year), the slope ω measuring the temporal succession rate of microbial communities along time was the TDR value, c was the intercept, and ε was the residual. A bootstrapping of 999 times was used to determine the significance of the TDR value. The significance of differences between the native treatment and the invasive treatment was calculated by comparing observed TDR differences and those in permuted datasets. These analyses were conducted for selected function groups or genera with relative abundance > 0.05%. Functional groups for bacterial communities were predicted using the FAPROTAX method [36], and those for fungi communities were predicted using the FUNGuild method39.
To construct the co-occurrence network, Spearman’s correlation coefficients between OTUs with relative abundance > 0.10% were first calculated using the “microeco” package (version 1.3.0) in R. The Spearman correlation threshold was set using random matrix theory (RMT) with a significance level of P < 0.01. We used the igraph package (version 1.4.2) to perform various network analyses. We then performed a one-way analysis of variance (ANOVA) to determine the significance for differences in subnetwork properties between treatments. Networks were then visualized using Gephi software (version 0.10, https://gephi.org/). The value zi describes how well the node i is connected to other nodes in the same module, and Pi reflects how node i connects to different modules. Nodes with zi > 2.5 and Pi ≤ 0.62 were considered as module hubs that were highly connected to many nodes in their own modules and zi ≤ 2.5 and Pi > 0.62 as connectors that were highly linked to several modules [38].
Other statistical analyses were carried out using different packages in R software with details as follows. The β diversities were determined using the Bray–Curtis index [39]. We used a null model analysis to evaluate the relative importance of stochastic and deterministic processes in assembling microbial communities as described by Chase and Myers [40]. Microbial stochastic ratios were determined using the iCAMP (V1.5.3), the ape (V5.3), and the NST (V3.1.8) packages [41]. Simple Mantel tests were conducted with the vegan (V2.3–5) package [42], and the multiple regression on distance matrices (MRM) was conducted with the ecodist (V2.0.9) package. Random forest analysis was performed to identify the individual effects of each variable on soil microbial communities using the randomForest (V4.7–1.1) package45. The importance of each factor was estimated by the percentage increase in the mean squared error of variables using the rfPermute (V2.5.2) package by performing 500 permutations on the response variable for evaluation.
Results
Impacts of alien plant invasion on the structure and diversity of microbial communities
Alien plant invasion significantly (Adonis R2 = 0.04–0.08, P < 0.001) altered both bacterial (Fig. 1b) and fungal (Fig. 1c) community structures, as revealed by the nonmetric multidimensional scaling (NMDS) ordination. Microbial communities in both invasive and native plants were also significantly (P < 0.05) different from the baseline treatment in 2012 (Table S1). To test the species effect on microbial community structure, we conducted DMDS for A. artemisiifolia and B. pilosa separately and found different community structures between invasive and native species for both bacterial and fungal communities (Adonis R2 = 0.02–0.07, P < 0.02, Fig. S1).
Bacterial Shannon diversity was 7.67 ± 0.01 in the invasive treatment, which was significantly (P < 0.05) higher than those in the native treatment (7.62 ± 0.02, Fig. 1d, Table S2). Fungal Shannon diversity was also significantly higher in the invasive treatment (5.03 ± 0.05) than in the native treatment (4.84 ± 0.08). Similarly, other diversity indices, including richness and Simpson indices, all increased due to plant invasion (Table S2).
Impacts of alien invasion on microbial temporal succession
To understand the impacts of alien plant invasion on the temporal succession of soil microbial communities, we examined the slopes (ω) of the bacterial or fungal time-decay relationships (TDRs). Significant bacterial TDRs (P < 0.05) were observed 8 years after experiment setup, with a ω value of 0.05 for the invasion treatment, which was twice that of the native treatment (Fig. 2a). For fungal communities, alien plant invasion also resulted in a significant successive ω value of 0.07, whereas such succession was insignificant in the native treatment (Fig. 2b). To further dissect microbial temporal succession in different time periods, we estimated the annual TDRs from the 5th to the 8th year. Microbial succession in both the native and the invasive treatments was significant (P < 0.05) from the 5th to the 8th year (Fig. S2). Remarkably, bacterial ω values were smaller in the invasive treatment than those in the native treatment until the 7th year but were larger in the invasive treatment in the 8th year, underscoring the importance of temporal dynamics.
Time-decay relationships (TDRs) at various levels of microbial communities, including (a) the entire bacterial community, b the entire fungal community, c functional groups, and (d) selected genera. The significance (asterisk) of the ω-value is determined by bootstrapping 999 times. The significance of differences between treatments was then calculated by comparing the observed slope differences with those in permuted datasets. The letter “A” in diagram (c) indicates a significant difference (P < 0.05) between the native treatment and the invasion treatment. In diagram d, the rates are scaled between − 2 and 2 in the heatmap; the unscaled mean rates are represented by green bars for native plants, red bars for invasive A. artemisiifolia, and blue bars for invasive B. Pilosa; asterisks over red or blue bars indicate significantly different rates in the invasion treatment compared to the native treatment; the scatter plot refers to the mean rates in each plant. Abbreviations: A. + C. — intercropping of A. artemisiifolia and C. serotinum; A. + S. — intercropping of A. artemisiifolia and S. viridis; B. + C. — intercropping of B. pilosa and C. serotinum; B. + S. — intercropping of B. pilosa and S. viridis. *P < 0.05, **P < 0.01, ***P < 0.001
We further investigated TDRs in 12 major functional groups, including carbon and nitrogen cycling groups, arbuscular mycorrhizal fungi (AMF), and pathogenic fungi. Half of these functional groups exhibited significant (P < 0.05) TDRs in the native treatment, such as nitrogen fixation, nitrification, nitrate respiration, and plant pathogen (Fig. 2c). Alien plant invasion induced significant succession in carbon cycling groups and AMF and amplified 2–3 times (P < 0.05) higher TRDs in nitrogen fixation and denitrification groups.
We then examined microbial TDRs at the genus level in response to the invasion of A. artemisiifolia and B. pilosa, respectively. Only seven genera exhibited significant (P < 0.05) temporal succession across all treatments, including Azotobacter performing nitrogen fixation, Virgisporangium and Gp7 within Acidobacteria, and two pathogenic fungi, Aspergillus and Metarhizium (Table S3). Almost all genera in cluster 1, revealed by a heatmap of TDRs’ slopes, showed significantly higher TDRs in A. artemisiifolia or B. pilosa than in the native treatment (Fig. 2d). By contrast, most genera in cluster 3 had lower turnover rates in the invasive treatment compared to the native treatment, including Cyphellophora, Penicillium, Nocardia, and Paraphoma. Some microbial genera had significantly higher turnover rates only in the A. artemisiifolia invasive treatment, such as Ilumatobacter, Modestobacter, Actinoplanes, and Caulobacter, while some other genera, such as Anaeromyxobacter and Rhizobium, exhibited significantly higher turnover rates only in the B. pilosa invasive treatment. These results indicated that the two invasive plants differed in their abilities to modify certain microbial genera.
Community assembly processes underlying microbial successional trajectories
To determine whether microbial communities converge or diverge during succession, we calculated the distances between the native and invasive microbial communities yearly. The variations in bacterial communities increased linearly over time (r = 0.16, P = 0.004, Fig. 3a), showing a pattern of divergence. In contrast, this divergent pattern was not significant for fungal communities. For bacterial community, stochastic processes accounted for 42.53–49.78% of the community variations, with alien plant invasion significantly decreasing the stochasticity by about 5.00% (Fig. 3b). Furthermore, the relative contribution of the stochastic process in shaping bacteria communities increased over time in the native treatment but remained consistent in the invasion treatment (Fig. S3a). In contrast, stochastic ratios were similar between the native and the invasion treatments in shaping fungal communities, with stochasticity ranging from 70.12 to 79.28% (Fig. 3b, Fig. S3b).
a Community distance between the native and invasive microbial communities based on Bray–Curtis and b community stochasticity. Statistical significance in pairwise distance was determined by permutation tests. The Student t-test was used to calculate the statistical difference in stochasticity between the native treatment and the invasion treatment. *P < 0.05; ns, not significant
Impacts of alien plant invasion on microbial co-occurrence networks
Recognizing that bacteria and fungi communicate with each other in soil, we constructed co-occurrence networks that encompassed both bacterial and fungal communities. The networks differed profoundly, with the invasive network being more complex (Fig. 4). Using the same threshold for network construction, we obtained similar nodes of 135–160 in these two networks, in which 38.52–48.75% belonged to Ascomycota phyla and 28.75–36.30% belonged to Acidobacteria, Actinobacteria, and Proteobacteria. However, edges in the invasive treatment were 1256, nearly 3 times the number observed in the native treatment (Table S4). The average degree and network centralization were also about 2–3 times higher than those in the native treatment. Furthermore, the network diameter in the invasive treatment was one-third of that in the native treatment.
Co-occurrence networks of microbial communities in the native treatment (left panel) and the invasion treatment (right panel). Connections represent significant (Spearman’s r2 > 0.24, P < 0.01) correlations. The label of each node is colored according to taxonomic affiliation, and its size is proportional to the number of connections (i.e., degree). Microbial communities include both bacterial and fungal communities. Labels with ITS indicate that these OTUs are derived from fungal communities
According to the zi-Pi plot, OTU36 belonging to Ramlibacter emerged as a module hub in both native and invasive treatments (Fig. 4, Table S5). In the native treatment, there were 10 nodes identified as module hubs or connectors, including members of Gp4 and Blastocatella belonging to Acidobacteria and plant-pathogenic fungi Periconia, Cladosporium, and Penicillium. In contrast, the invasive network had 21 nodes functioning as connectors, almost all belonging to the phyla Ascomycota, Proteobacteria, and Actinobacteria.
Environmental drivers of microbial communities
Mantel tests were used to determine the relationship between each environmental variable (Table S6) and microbial communities. Soil characteristics, including soil organic matter, soil pH, total nitrogen, total phosphorus, and total potassium, were correlated with bacterial and fungal communities (P < 0.05, Table S7). The variation of bacterial community was also significantly explained by plant biomass (r = 0.15, P < 0.05). The random forest model further confirmed these observations (Fig. S4). These variables collectively contributed to 42.40% of bacterial community variation and 23.11% of fungal community variation (P < 0.05, Table S7). However, a large portion of fungal community variation could not be explained by measured variables, signifying the importance of stochastic processes (70.12–79.28%, Fig. 3b) in assembling fungal communities.
Discussion
Although numerous studies have demonstrated how soil microbes drive plant invasion through both direct routes (such as pathogens and mutualists) and indirect routes (such as mediation of nutrient availability) [7, 9, 15], it remains unclear whether these soil microbes respond dynamically to invasion over time or if microbial communities respond uniformly. Our results indicate that invasive alien plants significantly impact the composition of soil microbial communities over time. Specifically, alien plant invasion enhances microbial diversity, accelerates microbial succession, and intensifies microbial associations. Furthermore, invasive alien plants often activate beneficial microbiota that increase resource availability and suppress soil-borne pathogens.
Alien plant invasion leads to nonuniform succession dynamics in microbial communities
Alien invasive plants can recruit mutualistic microbes, such as nitrogen-fixation microbes and AMF, from native mutualist communities to support their establishment [44]. Consistently, we found that functional groups and genera with significantly higher succession rates in invasive plants compared to native plants typically play crucial roles in maintaining plant growth. For example, microbial groups involved in carbon and nitrogen cycling, genera in maintaining plant growth such as Variovorax [45] and Pseudolabrys [46], carbon-decomposing microbes like Actinomadura [47] and Preussia [48], phototrophic bacteria such as Cyanobacteria [49], and other beneficial taxa like Gemmatimonas, Pseudonocardia, and Simplicillium [50,51,52] had higher succession rates in the invasive treatment. This evidence suggests that soil microbial communities in alien invasive plants underwent succession over time to a unique edge in cultivating beneficial microbiota that increases resource availability.
Conversely, plants typically defend against pathogen attacks by utilizing specific receptors to identify pathogens [53]. It has been demonstrated that invasive alien plants are infected by 84% fewer fungi in invaded habitats than in their native ranges [54]. This observation supports the “enemy release hypothesis,” [5] which proposes that invasive plants are liberated from natural enemies, including pathogens in invasive areas. Similarly, many genera in our study with lower turnover rates in the invasive treatment were pathogenic microbes. Such examples included Cyphellophora causing sooty blotch and flyspeck disease in plants [55], soil-borne pathogen Plectosphaerella [56], Acrocalymma causing crown and root disease [57], Penicillium causing a postharvest plant disease named green mold [58], the pathogenic bacteria Nocardia [59], and the soil-borne pathogen Paraphoma [60].
The asymmetric succession dynamics between beneficial microbiota and pathogenic microbes support our hypothesis that alien plant invasion leads to a nonuniform succession dynamic in various microbial groups, specifically activating those that promote plant growth over time. One possible reason for this nonuniform succession is that plants recruit beneficial microbiota through root exudates and molecules recognized by these microbiota [53]. Certain exudates act as selective recruiters for specific microbial taxa or enrich the rhizosphere with organic carbon, attracting beneficial microbiota. Many invasive alien plants, including A. artemisiifolia and B. pilosa in our study, secrete allelopathic compounds from their roots [6, 61]. These compounds not only pose a threat to native flora due to their strong phytotoxic activities but also serve as carbon substrates in the soil, attracting specific microbiota. For instance, root residues in A. artemisiifolia are rich in sesquiterpenes [61], which increase the abundance of bacterial functional groups involved in carbon and nitrogen metabolism [62]. Additionally, the higher total carbon content in the litter of invasive alien plants compared to native plants may provide more opportunities for decomposer proliferation. For example, the carbon content in the litter of B. pilosa is 822.94 g/kg, substantially higher than in S. viridis (737.11 g/kg)65. To further investigate the reasons behind the asymmetric succession in soil microbial communities, metabonomic analyses and laboratory experiments are needed to identify chemicals that attract beneficial microbiota or suppress pathogenic microbes.
Microbial communities undergo divergent succession trajectories
Long-term experiment systems provide insights into the directionality and convergence of community assembly [33]. We expected a deterministic process to dominant the structuring microbial succession in the invasive treatment due to selective pressure from root exudates [64]. Consistently, alien plant invasion significantly increased the deterministic ratio in assembling bacterial communities. While the stochasticity is higher for fungi than for bacteria in our study, this could be explained by high asexual reproduction, rapid dispersal rates, and resistance to extinction in fungi [65]; furthermore, the much larger size of fungi compared to bacteria may result in lower cell counts and thus less diverse communities [66]. Under deterministic perspective, microbial succession is directional and dissimilar between treatments [67]. Consequently, we observed divergent successional trajectories between native and invasive treatments. Furthermore, our results demonstrate that microbial communities in the invasive treatment evolved towards a particular state that potentially benefits plant growth. This phenomenon highlights the intricate interplay between microbial dynamics and alien plant invasions, contributing to the ecological success of invasive alien species. Network analysis is a valuable tool for identifying and understanding species associations and interactions in complex systems [34, 68]. Higher average degree and average density indicate a more complex network [38], while larger network centralization reflects a more centralized topology. All these topology indices were higher in the invasive treatment than in the native treatment, suggesting that alien plant invasion compacted microbial communities, resulting in nodes being more closely connected. This increased complexity may enhance the adaptability and resilience of invasive alien plants in various environments, underscoring the potential role of microbial interactions in the ecological success of invasive species.
In ecological network analyses, module hubs and connectors are critical elements. Module hubs, with numerous connections within their modules, and connectors, with high links to multiple modules, are identified as microbial keystone taxa [69]. These keystone taxa typically play important roles in communities. In our study, Ramlibacter was a keystone taxon in both native and invasion treatments, which plays a pivotal role in soil N cycling [70]. Additionally, Periconia, Cladosporium, and Penicillium, which are plant-pathogenic fungi, were among the 10 keystone taxa of the native treatment. Therefore, alien plant invasion not only causes divergent microbial compositions but also leads to divergence in soil microbial community interactions.
Conclusion and implications
Our findings provide robust evidence supporting the dynamic aspect of plant microbiome assembly (Fig. 5). Alien plants enhance microbial diversity, accelerate microbial succession, form beneficial associations, and mitigate soil pathogen limitations. Notably, turnover rates vary among microbiota, reflecting their distinct ecological roles. This comprehensive perspective suggests that the spread of invasive alien plants is facilitated by a strategically assembled microbial support system, offering a comprehensive and multifaceted view of the invasive mechanisms. Understanding the microbial interactions that contribute to alien plant invasion can inform the development of biological control methods, such as introducing natural enemies of specialized pathogens or herbivores. While Asteraceae species are among the most predominant invasive alien plants, the asymmetric succession in soil microbial communities requires verification through additional research involving a broader range of plant categories. With sufficient information, we may be able to predict the ecological consequences of microbial community succession, particularly in scenarios where invasive plants are leading to a global homogenization of plant species.
The asymmetric succession in soil microbial communities due to alien plant invasion. The invasion of alien plants diversifies soil microorganisms and enhances their associations through litter and root effects. It also accelerates the divergent succession of microbial communities by activating beneficial microbiota related to nutrient cycling and mutualistic symbiosis
Data availability
No datasets were generated or analysed during the current study.
References
McKinney ML, Lockwood JL. Biotic homogenization: a few winners replacing many losers in the next mass extinction. Trends Ecol Evol. 1999;14(11):450–3.
Zhao M, Lu X, Zhao H, Yang Y, Hale L, Gao Q, et al. Ageratina adenophora invasions are associated with microbially mediated differences in biogeochemical cycles. Sci Total Environ. 2019;677:47–56.
Ehrenfeld JG. Ecosystem consequences of biological invasions. Ann Rev Ecology Evolution and Systematics. 2010;41(1):59–80.
Colautti R, Grigorovich I, MacIsaac H. Propagule pressure: a null model for biological invasions. Biol Invasions. 2006;8(5):1023–37.
Keane RM, Crawley MJ. Exotic plant invasions and the enemy release hypothesis. Trends Ecol Evol. 2002;17(4):164–70.
Bais HP, Vepachedu R, Gilroy S, Callaway RM, Vivanco JM. Allelopathy and exotic plant invasion: from molecules and genes to species interactions. Science. 2003;5638:301.
Callaway RM, Aschehoug ET. Invasive plants versus their new and old neighbors: a mechanism for exotic invasion. Science. 2000;290(5491):521–3.
Martin FM, Uroz S, Barker DG. Ancestral alliances: plant mutualistic symbioses with fungi and bacteria. Science. 2017;356(6340):eaad4501.
Liu B, Yan J, Li W, Yin L, Li P, Yu H, et al. Mikania micrantha genome provides insights into the molecular mechanism of rapid growth. Nat Commun. 2020;11(1):340.
Zhao P, Liu B, Zhao H, Lei Z, Zhou T. Significant changes in soil microbial community structure and metabolic function after Mikania micrantha invasion. Sci Rep. 2023;13(1):1141.
Li H, Wei Z, Huangfu C, Chen X, Yang D. Litter mixture dominated by leaf litter of the invasive species, Flaveria bidentis, accelerates decomposition and favors nitrogen release. J Plant Res. 2017;130(1):167–80.
Sardans J, Bartrons M, Margalef O, Gargallo-Garriga A, Janssens IA, Ciais P, et al. Plant invasion is associated with higher plant-soil nutrient concentrations in nutrient-poor environments. Glob Change Biol. 2017;23(3):1282–91.
Yan J, Zhang XY, Chen X, Wang Y, Zhang FJ, Wan FH. Effects of rhizosphere soil microorganisms and soil nutrients on competitiveness of Bidens pilosa with different native plants. Biodiv Sci. 2016;24(12):1381–9.
Zhang FJ, Li Q, Chen FX, Xu HY, Wan FH. Arbuscular mycorrhizal fungi facilitate growth and competitive ability of an exotic species Flaveria bidentis. Soil Biol Biochem. 2017;115:275–84.
Dawson W, Schrama M, Austin A. Identifying the role of soil microbes in plant invasions. J Ecol. 2016;104(5):1211–8.
Crawford KM, Knight TM. Competition overwhelms the positive plant–soil feedback generated by an invasive plant. Oecologia. 2016;183(1):211–20.
Ali HE, Al-Wahaibi AM, Shahid MS. Plant–soil feedback and plant invasion: effect of soil conditioning on native and invasive Prosopis species using the plant functional trait approach. Front Plant Sci. 2024;15:1321950.
Walker LR, Wardle DA, Bardgett RD, Clarkson BD. The use of chronosequences in studies of ecological succession and soil development. J Ecol. 2010;98(4):725–36.
Zhang P, Li B, Wu J, Hu S. Invasive plants differentially affect soil biota through litter and rhizosphere pathways: a meta-analysis. Ecol Lett. 2019;22(1):200–10.
Blumenthal D, Mitchell CE, Pyšek P, Jarošík V. Synergy between pathogen release and resource availability in plant invasion. Proc Natl Acad Sci. 2009;106(19):7899–904.
Lu S, Luo X, Wang H, Gentili R, Citterio S, Yang J, et al. China-US grain trade shapes the spatial genetic pattern of common ragweed in east China cities. Commun Biol. 2023;6(1):1072.
Zhang K, Ebihara A, Tong S, White JC, Shen Y. Bidens pilosa root exudates modulate Pteris multifida gametophyte development: a proteomic investigation. Ind Crop Prod. 2023;205:117499.
Zhang F, Sun J, Wang C, Li C, Chen F, Xu H, et al. Bacillus benefits the competitive growth of Ambrosia artemisiifolia by increasing available nutrient levels. Front Plant Sci. 2022;13:1069016.
Lu R. Soil agricultural chemical analysis methods. Beijing: Agricultural Sci-Tech Press; 1999.
Ciavatta C, Govi M, Antisari LV, Sequi P. Determination of organic carbon in aqueous extracts of soils and fertilizers. Commun Soil Sci Plant Anal. 1991;22(9–10):795–807.
Khan J, Wei JS, Ringnér M, Saal LH, Ladanyi M, Westermann F, et al. Classification and diagnostic prediction of cancers using gene expression profiling and artificial neural networks. Nat Med. 2001;7(6):673–9.
Olsen SR, Cole C, Watanabe FS, Dean L. Estimation of available phosphorus in soils by extraction with sodium bicarbonate, vol. 939. Washington, DC: US Department of Agriculture; 1954.
Zhao M, Sun B, Wu L, Gao Q, Wang F, Wen C, et al. Zonal soil type determines soil microbial responses to maize cropping and fertilization. mSystems. 2016;1(4):e00075-00016.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.
Kõljalg U, Nilsson RH, Abarenkov K, Tedersoo L, Taylor AF, Bahram M, et al. Towards a unified paradigm for sequence-based identification of fungi. Mol Ecol. 2013;22(21):5271–7.
Guo X, Feng J, Shi Z, Zhou X, Yuan M, Tao X, et al. Climate warming leads to divergent succession of grassland microbial communities. Nat Clim Chang. 2018;8(9):813–8.
Gao Q, Yang Y, Feng J, Tian R, Guo X, Ning D, et al. The spatial scale dependence of diazotrophic and bacterial community assembly in paddy soil. Glob Ecol Biogeogr. 2019;28(8):1093–105.
Barberan A, Bates ST, Casamayor EO, Fierer N. Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 2012;6(2):343–51.
Sansupa C, Wahdan SFM, Hossen S, Disayathanoowat T, Wubet T, Purahong W. Can we use functional annotation of prokaryotic taxa (FAPROTAX) to assign the ecological functions of soil bacteria? Appl Sci. 2021;11(2): 688.
Nguyen NH, Song Z, Bates ST, Branco S, Tedersoo L, Menke J, et al. FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 2016;20:241–8.
Deng Y, Jiang Y-H, Yang Y, He Z, Luo F, Zhou J. Molecular ecological network analyses. BMC bioinformatics. 2012;13(1):113.
Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27(4):325–49.
Chase JM, Myers JA. Disentangling the importance of ecological niches from stochastic processes across scales. Philos Trans R Soc B: Biol Sci. 2011;366(1576):2351–63.
Ning D, Deng Y, Tiedje JM, Zhou J. A general framework for quantitatively assessing ecological stochasticity. Proc Natl Acad Sci. 2019;116(34):16892–8.
Smouse PE, Long JC, Sokal RR. Multiple regression and correlation extensions of the Mantel test of matrix correspondence. Syst Zool. 1986;35:627–32.
Zhang Z, Lu Y, Wei G, Jiao S. Rare species-driven diversity–ecosystem multifunctionality relationships are promoted by stochastic community assembly. mBio. 2022;13(3):e00449–e00422.
Richardson DM, Pysek P, Rejmanek M, Barbour MG, Panetta FD, West CJ. Naturalization and invasion of alien plants: concepts and definitions. Divers Distrib. 2000;6(2):93–107.
Finkel OM, Salas-González I, Castrillo G, Conway JM, Law TF, Teixeira PJPL, et al. A single bacterial genus maintains root growth in a complex microbiome. Nature. 2020;587(7832):103–8.
Kämpfer P, Young CC, Arun A, Shen FT, Jäckel U, Rosselló-Mora R, et al. Pseudolabrys taiwanensis gen. nov., sp. nov., an alphaproteobacterium isolated from soil. Int J Syst Evol Microbiol. 2006;56(10):2469–72.
Kroppenstedt RM, Goodfellow M. The family thermomonosporaceae: actinocorallia, actinomadura, spirillospora and thermomonospora. Handbook Biol Bacteria. 2006;3:682–724.
Gonzalez-Menendez V, Martin J, Siles JA, Gonzalez-Tejero MR, Reyes F, Platas G, et al. Biodiversity and chemotaxonomy of Preussia isolates from the Iberian Peninsula. Mycol Prog. 2017;16(7):713–28.
Whitton BA, Potts M. Introduction to the cyanobacteria. In: Ecology of cyanobacteria II: their diversity in space and time. Dordrecht: Springer Netherlands; 2012. p. 1–13.
Chee-Sanford J, Tian D, Sanford R. Consumption of N2O and other N-cycle intermediates by Gemmatimonas aurantiaca strain T-27. Microbiology. 2019;165(12):1345–54.
Currie CR, Scott JA, Summerbell RC, Malloch D. Fungus-growing ants use antibiotic-producing bacteria to control garden parasites. Nature. 1999;398(6729):701–4.
Agrawal S, Saha S. The genus Simplicillium and Emericellopsis: a review of phytochemistry and pharmacology. Biotechnol Appl Biochem. 2022;69(5):2229–39.
Santoyo G. How plants recruit their microbiome? New insights into beneficial interactions. J Adv Res. 2022;40:45–58.
Mitchell CE, Power AG. Release of invasive plants from fungal and viral pathogens. Nature. 2003;421(6923):625–7.
Gao L, Ma Y, Zhao W, Wei Z, Gleason ML, Chen H, et al. Three new species of Cyphellophora (Chaetothyriales) associated with sooty blotch and flyspeck. PLoS One. 2015;10(9):e0136857.
Carlucci A, Raimondo M, Santos J, Phillips A. Plectosphaerella species associated with root and collar rots of horticultural crops in Southern Italy. Persoonia-Mol Phyl Evol Fungi. 2012;28(1):34–48.
Alcorn J, Irwin J. Acrocalymma medicaginis gen. et sp. nov. causing root and crown rot of Medicago sativa in Australia. Trans Br Mycol Soc. 1987;88(2):163–7.
Luciano-Rosario D, Keller NP, Jurick WM. Penicillium expansum: biology, omics, and management tools for a global postharvest pathogen causing blue mould of pome fruit. Mol Plant Pathol. 2020;21(11):1391–404.
Mehta HH, Shamoo Y. Pathogenic Nocardia: a diverse genus of emerging pathogens or just poorly recognized? PLoS Pathog. 2020;16(3): e1008280.
Cao S, Liang Q, Nzabanita C, Li Y. Paraphoma root rot of alfalfa (Medicago sativa) in Inner Mongolia. China Plant Pathol. 2020;69(2):231–9.
Liu Z, Zhang N, Ma X, Zhang T, Li X, Tian G, et al. Sesquiterpenes from Ambrosia artemisiifolia and their allelopathy. Front Plant Sci. 2022;13:996498.
Yu H, Le Roux JJ, Zhao M, Li W. Mikania sesquiterpene lactones enhance soil bacterial diversity and fungal and bacterial activities. Biol Invasions. 2022;25(1):237–50.
Li Q, Guo JY, Zhang H, Zhao MX. The competition between Bidens pilosa and Setaria viridis alters soil microbial composition and soil ecological function. J Integr Agr. 2024;23(1):267–82.
Steinauer K, Thakur MP, Emilia Hannula S, Weinhold A, Uthe H, van Dam NM, et al. Root exudates and rhizosphere microbiomes jointly determine temporal shifts in plant-soil feedbacks. Plant Cell Environ. 2023;46(6):1885–99.
Wang J, Liu G, Zhang C, Wang G, Fang L, Cui Y. Higher temporal turnover of soil fungi than bacteria during long-term secondary succession in a semiarid abandoned farmland. Soil Tillage Res. 2019;194:104305.
Rousk J, Bååth E. Fungal biomass production and turnover in soil estimated using the acetate-in-ergosterol technique. Soil Biol Biochem. 2007;39(8):2173–7.
Li S-p, Cadotte MW, Meiners SJ, Pu Z, Fukami T, Jiang L. Convergence and divergence in a long-term old-field succession: the importance of spatial scale and species abundance. Ecol Lett. 2016;19(9):1101–9.
Xiao N, Zhou A, Kempher ML, Zhou BY, Shi ZJ, Yuan M, et al. Disentangling direct from indirect relationships in association networks. Proceed Natl Acad Sci. 2022;119(2):e2109995119.
Eiler A, Heinrich F, Bertilsson S. Coherent dynamics and association networks among lake bacterioplankton taxa. ISME J. 2012;6(2):330–42.
Hu X, Gu H, Sun X, Wang Y, Liu J, Yu Z, et al. Distinct influence of conventional and biodegradable microplastics on microbe-driving nitrogen cycling processes in soils and plastispheres as evaluated by metagenomic analysis. J Hazard Mater. 2023;451:131097.
Acknowledgements
We thank Daliang Ning from the University of Oklahoma for assistance in the R code.
Funding
This research was supported by grants from the National Key R&D Program of China (2021YFC2600400, 2022YFC2601100) and the National Science Foundation of China (42207162).
Author information
Authors and Affiliations
Contributions
FW and MZ developed the study concept and design; QL collected the samples and completed the DNA preparation and environmental genomics experiments; MZ performed data analyses and drafted the manuscript, HZ and XZ performed the data illustration. WL, YY and XG revised the manuscript. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Zhao, M., Yang, Y., Zhang, H. et al. Asymmetric succession in soil microbial communities enhances the competitive advantage of invasive alien plants. Microbiome 12, 265 (2024). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40168-024-01989-5
Received:
Accepted:
Published:
DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40168-024-01989-5