Pengyuan Sun1
Weihua Guo
94% of researchers rate our articles as excellent or good
Learn more about the work of our research integrity team to safeguard the quality of each article we publish.
Find out more
- 1Qingdao Key Laboratory of Ecological Protection and Restoration, Ministry of Natural Resources Key Laboratory of Ecological Prewarning, Protection and Restoration of Bohai Sea, School of Life Sciences, Shandong University, Qingdao, China
- 2State Key Laboratory of Water Environment Simulation, School of Environment, Beijing Normal University, Beijing, China
Soil microbial communities are integral to almost all terrestrial biogeochemical cycles, which are essential to coastal wetland functioning. However, how soil bacterial community assembly, composition, and structure respond to native and non-native plant invasions in coastal wetlands remains unclear. In this study of the coastal wetlands of the Yellow River Delta in China, the assembly, community composition, and diversity of soil bacterial communities associated with four wetland plant species (Phragmites australis, Spartina alterniflora, Suaeda salsa, and Tamarix chinensis) and four soil depths (0–10 cm, 10–20 cm, 20–30 cm, and 30–40 cm) were characterized using high-throughput sequencing. Plant species identity, as well as environmental factors, rather than soil depth, was found to play predominant roles in shaping the diversity and structure of wetland soil bacterial communities. S. alterniflora invasion altered bacterial community structure and increased bacterial diversity. Phragmites australis-associated bacterial communities were enriched with sulfate-reducing bacteria such as Desulfurivibrio and Desulfuromonas. In comparison, S. alterniflora-associated bacterial communities were enriched with both sulfate-reducing bacteria (SEEP-SRB1) and sulfate-oxidizing bacteria (Sulfurimonas), which maintained a dynamic balance in the local sulfur-cycle, and thereby enhanced S. alterniflora growth. In addition, stochastic processes dominated the assembly of soil bacterial communities associated with all four plant species, but were most important for the S. alterniflora community. The S. alterniflora-associated bacterial community also showed stronger interactions and more extensive connections among bacterial taxa; a co-occurrence network for this community had the greatest average clustering coefficient, average degree, modularity, and number of links and nodes, but the lowest average path length. Altogether, individual plant species had distinct effects on soil bacterial community assembly and structure, with the invasive species having the strongest impact. These results provide insights into microbial ecology and inform management strategies for coastal wetland restoration.
1 Introduction
Wetlands are among the most sensitive ecosystems to global climate change. Despite covering only 5–8% of the global land area, wetlands store 20–30% of all soil carbon (C) in terrestrial ecosystems (Tan et al., 2020; Wang et al., 2019). As dynamic zones bridging terrestrial and marine environments, coastal wetlands are critical for biodiversity maintenance, biogeochemical cycling, and climate regulation (Li et al., 2024; Osland et al., 2022). However, land-use changes and sea-level rise over the past few decades have resulted in the loss of large areas of coastal wetlands globally, severely weakening their capacity to provide ecosystem services (Osland et al., 2022; Wang et al., 2019). Soil microbial communities are key drivers of biogeochemical cycles in terrestrial ecosystems, and understanding their responses to environmental changes is crucial for maintaining the functioning and stability of coastal wetlands (Zhou et al., 2020).
Soil bacterial communities are highly sensitive to environmental change (Hermans et al., 2020) and are often shaped by the local plant community, with shifts in the abundance of specific taxa explaining the “home-field advantage” effect on litter decomposition (Philippot et al., 2013). Plants directly affect soil microbes via the provision of rhizodeposits and leaf litter (Cline and Zak, 2015; Philippot et al., 2013), as well as indirectly by altering the rhizosphere environment (Zhalnina et al., 2018). In coastal wetlands, native plant communities are threatened by invasive plant species and significant losses of native plant diversity (Zhang G. et al., 2021). These changes in wetland plant community composition may greatly alter the soil microenvironment, for example by affecting nutrient availability and the soil pH (Sardans et al., 2017; Zhang et al., 2019), which then influences soil bacterial community composition and structure. In addition, plant community effects on soil bacterial communities may vary with soil depth due to variability in the distribution of plant roots, thereby affecting the amount of plant litter and rhizodeposits across soil layers (Leewis et al., 2022). Variation in oxygen availability across soil layers further shapes bacterial community composition (Tian et al., 2022). Therefore, investigating how bacterial communities respond to changes in the vegetation type and soil depth is crucial for preserving the ecological functions of wetlands.
Community-assembly studies help explain species coexistence and the maintenance of diversity (Zhou and Ning, 2017). Microbial community assembly may be influenced by both deterministic and stochastic ecological processes (Zhou et al., 2014; Chen et al., 2019). Deterministic processes include environmental selection (i.e., habitat filtering) and species interactions such as antagonism and competition (Webb et al., 2002). Stochastic processes include random birth-death events, probabilistic dispersal, and unpredictable disturbances (Chase and Myers, 2011; Zhou and Ning, 2017). Soil bacterial community assembly has been widely studied, but study findings are not currently conclusive, with both deterministic (Gao et al., 2021; Huang et al., 2022; Zhang et al., 2022) and stochastic assembly processes (Hussain et al., 2023; Yue et al., 2022; Zhang et al., 2023) identified in coastal wetlands. Other factors, including diverse soil physical and chemical properties, may also impact bacterial community assembly (Huo et al., 2023; Ni et al., 2021). Additionally, changes in plant species composition can affect soil processes and thereby regulate the composition and diversity of soil microbial community composition and diversity directly or indirectly (Ward et al., 2015; Shigyo et al., 2019). The relative contribution of deterministic versus stochastic processes may differ among bacterial communities varying in diversity, with deterministic processes tending to dominate in less diverse communities (Xun et al., 2019); therefore, changes in plant species composition may alter bacterial community assembly. Further investigation is needed to assess how soil bacterial community assembly responds to changes in plant species composition.
The Yellow River Delta (YRD) is a biodiversity hotspot in China that contains extensive wetland ecosystems (Jiang et al., 2022). The YRD encompasses diverse vegetation types, with a natural successional gradient from low-to-high tide zones; wetland plant species include Phragmites australis, Suaeda salsa, and Tamarix chinensis. However, due to climate change and human activities, the YRD is under threat from agricultural land reclamation, oil pollution, and non-native plant invasions, all of which put local biodiversity at risk (Yu et al., 2015; Gao et al., 2016). Since its introduction to the YRD in 1990, the rapid expansion of Spartina alterniflora has reduced the habitat available for native plant species, thereby threatening local biodiversity (Ren et al., 2019). Therefore, the S. alterniflora invasion, in addition to the natural successional gradient for native plants in the YRD, provides a good opportunity to examine differences in the assembly and structure of soil bacterial communities.
In this study, four sites containing different plant species (P. australis, S. salsa, T. chinensis, and the invasive S. alterniflora) in the coastal wetlands of YRD were selected to investigate how plant species composition and soil depth (0–10 cm, 10–20 cm, 20–30 cm, and 30–40 cm) affect soil bacterial community assembly and structure. For each site, the soil physicochemical properties were analyzed, high-throughput sequencing was performed for soil bacterial DNA, and a soil bacteria co-occurrence network was constructed. The study asked: (1) Which has a stronger influence on soil bacterial community structure, plant species identity or soil depth? and (2) Does the soil bacterial community associated with S. alterniflora differ from that of native plant species?
2 Materials and methods
2.1 Site description and soil collection
This study was conducted in the Yellow River Delta (YRD: 118° 33′ E–119° 20′ E, 37° 35′ N–38° 12′ N), located in Shandong Province, China. The YRD experiences a warm temperate continental monsoon climate with four seasons; the mean annual temperature ranges from 11.5–12.4°C, and the average annual precipitation from 530 to 630 mm (Lu et al., 2022). The study area represented a typical salt marsh. Four sampling sites were established in the salt marsh wetland (Figure 1), varying in their successional stage and consisting of either S. salsa (early successional stage, hereinafter referred to as SS), T. chinensis (mid-successional stage, TC), P. australis (late successional stage, PA), or S. alterniflora (invasive species, SA) (Zhang et al., 2007). These vegetation types are distributed along a sea-to-land gradient. S. alterniflora is closest to the sea and experiences tidal flooding during every high tide, whereas S. salsa and T. chinensis are in transitional zones with less frequent flooding. P. australis, located farthest from the sea near the Yellow River bank, is rarely affected by tidal flooding. Soil salinity, inferred from electrical conductivity (EC) (Liu et al., 2023), decreases along this gradient. Spartina alterniflora shows lower salinity due to frequent seawater flushing, while P. australis exhibits the lowest overall salinity near the freshwater source. This spatial and environmental gradient provides a unique natural laboratory for studying the ecological and biogeochemical processes in coastal wetlands.
Figure 1
Figure 1. Maps of the study location and sampling sites. “PA,” “TC,” “SS,” and “SA” represent the sampling sites for Phragmites australis, Tamarix chinensis, Suaeda salsa, and Spartina alterniflora, respectively.
Soil samples were collected from six plots (1 m × 1 m) within each site in September 2020; these were sectioned into four layers (representing soil depths of 0–10 cm, 10–20 cm, 20–30 cm, and 30–40 cm). Using a five-point sampling method, five subsamples were obtained per plot at each soil depth, and these were then thoroughly mixed to create a single representative sample. To avoid contamination with plant fragments, we first carefully removed plant litter from the soil surface before sampling. After the samples were brought back to the laboratory, any visible plant fragments were manually removed before further analysis. In total, 96 soil samples (four sampling sites × four soil depths × six replicates [plots]) were obtained and transported to the laboratory on ice. Each soil sample was divided into two parts: one was air-dried for later soil physicochemical analysis, and the second was frozen (−80°C) prior to microbial DNA extraction and high-throughput sequencing.
2.2 Soil physicochemical property analysis
The air-dried soil samples were sieved through 2 mm mesh to remove visible plant roots and pebbles. Soil samples were then oven-dried at 105°C to constant weight, after which they were weighed to determine the soil water content (SWC). Distilled water soil extracts (5:1 volume: weight) were used to measure the soil pH (FiveEasy Plus™, Mettler Toledo, Switzerland) and electrical conductivity (EC, FiveEasy Plus™, Mettler Toledo, Switzerland). Soil organic carbon (SOC) was determined via potassium dichromate oxidation with external heating (Nelson and Sommers, 1982). Soil dissolved organic carbon (DOC) was determined using a TOC analyzer (TOC-L CPH, Shimadzu, China) after extracting fresh soil samples with 2 mol L−1 KCl. Soil total nitrogen (TN) and total phosphorus (TP) were determined by the Kjeldahl method and molybdenum antimony colorimetric method, respectively. The soil nitrate nitrogen (NO3−) and ammonium nitrogen (NH4+) were quantified via a spectrophotometry-based colorimetric method after leaching samples in a 2 mol L−1 KCl solution.
2.3 DNA extraction and high-throughput sequencing
For microbial analysis, the soil was first passed through a 40-mesh (0.5 mm) sieve to remove finer plant debris, and then the DNA of soil microorganisms was extracted using a commercial DNA extraction kit (MagPure Soil DNA LQ Kit, Guangzhou Magen Biotechnology Co., Ltd.). The quality of the extracted DNA was assessed using a NanoDrop® ND-1000 spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA, United States) and 0.8% agarose gel electrophoresis. The primers 343F-5’-TACGGRAGGCAGCAG-3′ and 798R-5′-AGGGTA TCTAA TCCT-3′ were used to amplify the non-conserved region (16S V3–V4) of the 16S rRNA gene (Nossa et al., 2010). To ensure the quality of the downstream analysis, raw data obtained from Illumina MiSeq sequencing were processed to remove ambiguous bases in paired-end reads. Low-quality sequences with an average quality score of less than 20 and sequences of <50 bp were also removed. Next, quality-filtered reads were merged to obtain complete paired-end sequences with a maximum overlap of 200 bp. QIIME (version 1.8.0) and UCHIME (version 2.4.2) were used to remove sequences containing N-bases, reads with single base repetitions of more than eight, reads shorter than 200 bp, and reads containing chimeric sequences. Clean reads were clustered into OTUs using Vsearch (version 2.4.2) (Rognes et al., 2016) with a 97% similarity cutoff. The most abundant sequence for each OTU was selected as the representative sequence for that OTU. The Ribosomal Database Project Classifier, a naive Bayesian classification algorithm (Wang et al., 2007), was used to compare the representative sequences and annotate these via the Silva database (version 123) (Quast et al., 2012).
2.4 Network analysis and visualization
Microbial networks were constructed for each study site using the relative abundance of each bacterial genus (n = 24). Genera that were detected in more than 18 out of 24 samples and that had a relative abundance >0.01% were included in the network analysis. Connections between genera were evaluated using Spearman rank correlation tests as implemented in the R (version 4.2.1) package “Hmisc”. Significantly-related genera (r > 0.7, p < 0.05) were retained for network construction. Network topological properties, including the average path length, clustering coefficient, density, diameter, and modularity, were calculated in Gephi (version 0.9.3), which was also used for visualization of the co-occurrence networks. A random network analysis was performed using the R package “igraph”.
2.5 Statistical analysis
All statistical tests were considered significant at p < 0.05. Differences in bacterial community α-diversity were assessed using one-way analyses of variance (ANOVAs) in SPSS 26.0 Statistics (SPSS, Chicago, IL, United States). A principal coordinate analysis (PCoA) based on Bray-Curtis distances was performed using the “vegdist” function in the “vegan” package in R (version 4.2.1) to visualize differences in soil bacterial community composition among study sites. To determine which bacterial taxa varied significantly in abundance among sites, a linear discriminant analysis (LDA) effect size (LEfSe) analysis was run online.1 Pearson correlation coefficients were used to characterize how bacterial α-diversity was related to soil physicochemical properties, while Mantel tests were implemented in the “linkET” package in R to assess how soil physicochemical properties shaped bacterial community composition. The neutral model proposed by Sloan et al. (2006) was used to assess how stochastic processes affected bacterial community assembly. Spearman’s rank correlation was used to assess the relationships between the top 30 genera in all soil bacterial communities and various soil physicochemical properties across the four study sites, as shown in Supplementary Figure S2.
3 Results
3.1 Soil physicochemical properties
At a given soil depth, soil DOC, EC, NO3−, pH, SOC, and TN varied significantly among sites (p < 0.05, Supplementary Table S1). The soil pH, DOC, and TN were higher in PA than the other three sites, while the soil C/N and EC were lowest in PA. In addition, the soil DOC content was higher in PA and SA versus the other study sites. Compared to other sites, the soil NO3− content was the highest and the NH4+ content the lowest in SA. Within given site, differences in soil physicochemical properties were relatively small across soil depths (Supplementary Table S1). Several factors decreased with soil depth, including soil DOC in SS; SOC in TC; NH4+, SOC, and TN in PA; and NH4+ in SA.
3.2 Bacterial community diversity
A total of 33,158–66,617 high quality sequences were obtained after quality control and chimera removal, with the number of OTUs ranging from 3,608 to 8,134 per sample. Study site had a significant effect on the Chao1 index of the soil microbial community (p < 0.05), but no effect of soil depth was observed (p > 0.05). Specifically, SA had the highest Chao1 index at 8,505.75 ± 1,169.08, which was significantly higher than that of the other three sites (p < 0.05) (Figure 2A). However, the differences among PA, SS, and TC in the Chao1 index were not significant (p > 0.05). In contrast to study site, soil depth did not affect the Chao1 index (p > 0.05) (Figure 2A).
Figure 2
Figure 2. Chao 1 index (α-diversity) for bacterial communities from four sampling sites and soil depths (A). Lines between pairs of boxes indicate significant differences (p < 0.05) between sites (left) or soil depths (right). A principal coordinate analysis (PCoA) was performed for bacterial community composition (based on Bray-Curtis distances) to illustrate variation among sites and soil depths (B). SS, Suaeda salsa; TC, Tamarix chinensis; PA, Phragmites australis; SA, Spartina alterniflora. *p < 0.05.
To assess how study site and soil depth affected bacterial community β–diversity, a principal coordinate analysis (PCoA) was performed based on Bray-Curtis distances. Study site had a stronger effect on bacterial community composition than soil depth (ANOVA, Site: R2 = 0.254, p < 0.001; Soil depth: R2 = 0.046, p = 0.024). The bacterial communities from SA differed most significantly from those at the other three sites. In addition, the bacterial communities from PA were moderately different from those at SS and TC (Figure 2B).
3.3 Bacterial community composition
Taxonomic assignment of representative sequences for each OTU revealed a total of 59 phyla, of which the Proteobacteria (13.4–71.1%), Bacteroidetes (4.7–66.6%), Acidobacteria (1.2–13.7%), and Actinobacteria (1.8–13.8%) were the dominant phyla, together accounting for more than 80% of all reads across samples (Figure 3A). As expected, the effect of study site was stronger than that of soil depth for the Actinobacteria, Gemmatimonadetes, Nitrospirae, and Patescibacteria (p < 0.05; Supplementary Table S2). Notably, TC contained more Acidobacteria, Gemmatimonadetes, and Patescibacteria than the other sites (p < 0.05), while SS had a higher proportion of Actinobacteria (p < 0.05) and PA a greater relative abundance of Nitrospirae (p < 0.05). Finally, Acidobacteria and Calditrichaeota were significantly more abundant in the 30–40 cm soil layer (p < 0.05).
Figure 3
Figure 3. Bacterial community composition across study sites and soil depths. Relative abundance of the top ten most common bacterial phyla across sampling sites and soil depths (A). Most abundant bacterial genera at each sampling site (B); Heatmap of bacterial genera across sampling sites for the top 30 most abundant genera (C). SS, Suaeda salsa; TC, Tamarix chinensis; PA, Phragmites australis; SA, Spartina alterniflora.
A total of 1,694 genera were identified from the soil samples, and a heat map was created for the 30 most abundant genera (Figure 3C). The five most common genera across all samples were Desulfuromonas, Woeseia, Subgroup_10, Marinobacterium, and Nitrospira. However, bacterial genera showed differences in abundance among the four study sites (Figure 3C). The LEfSe analysis indicated that seven genera were biomarkers for SS (Figure 3B), with Marinobacterium, Pelobacter, and Pelagibius the most important of these. Subgroup_10, Kordiimonas, Marinobacter, Halomoans, and jdFR_76 were identified as biomarkers for TC, whereas six genera emerged as biomarkers for SA, including Subgroup_23, Gillisia, SEEP-SRB1, and Sulfurimonas. In addition, Desulfurivibrio, Nitrospira, Desulfuromonas, and Desulfurivibrio were the dominant genera at PA, while the relative abundance of Calorithrix was relatively low.
3.4 Assembly of the bacterial community
The soil bacterial communities at each study site were fit to a neutral community model. The model goodness-of-fit was better for SA (explaining 75.70% of the variation in composition) than the other three sites (Figure 4). The R2 of the null models at all four sites was greater than 0.50, indicating that stochastic processes strongly shaped bacterial community assembly in the coastal wetlands of the YRD. Moreover, the migration rate was higher for bacterial communities in SA (Nm = 11,149) than other sites (Figure 4D), indicating that bacterial genera in SA were less limited by dispersal.
Figure 4
Figure 4. The predicted frequency of bacterial OTUs based on a neutral community model for SS (A), TC (B), PA (C), and SA (D). Solid blue lines indicate the best fit to the neutral model, and dashed blue lines represent 95% confidence intervals around model predictions. R2 denotes the coefficient of determination (a measure of model fit), and Nm represents the product of metacommunity size and the migration rate. SS, Suaeda salsa; TC, Tamarix chinensis; PA, Phragmites australis; SA, Spartina alterniflora.
3.5 Co-occurrence networks for the bacterial communities
Four co-occurrence networks were constructed with genera plotted as nodes and the correlations between genera plotted as links (Figure 5). Overall, the number of positive links (74.27–88.94%) exceeded that of negative links (11.06–25.73%), indicating a preference for co-existence rather than co-exclusion among bacterial genera (Table 1). The co-occurrence network for SA had the greatest number of nodes (242) and links (1,287), and highest average degree (10.562), modularity (0.606), and average clustering coefficient (0.541), but the lowest average path length (3.793), suggesting more bacterial interactions and closer relationships among genera at this site. Meanwhile, the co-occurrence network for SS had the lowest modularity (0.379), reflecting significant differences in co-occurrence network topology for SS compared to the other three sites. Finally, the proportion of positive links (74.27%) was lowest in the PA co-occurrence network.
Figure 5
Figure 5. Co-occurrence networks for bacterial communities from SS (A), TC (B), PA (C), and SA (D) (n = 24). Different modules are shown in different colors. Lines connecting two nodes (edges) represent strong (r ≥0.7) and significant correlations (p < 0.05). Further details on network topology are provided in Table 1. SS, Suaeda salsa; TC, Tamarix chinensis; PA, Phragmites australis; SA, Spartina alterniflora.
Table 1
Table 1. Topological properties of the co-occurrence networks.
3.6 Relationships between soil physicochemical properties and bacterial community structure
The relationships between soil bacterial α-diversity or community composition and the soil physicochemical properties varied among sites (Figure 6). In TC, the bacterial community composition and Chao1 index were both significantly correlated with the soil EC and pH, while community composition was also correlated with TP. In SA, both the Chao1 index and community composition were correlated with NO3− and NH4+, while community composition was also significantly correlated with the soil depth and EC. In PA, the Chao1 index and community composition were significantly correlated with aspects of C stoichiometry (C/P for the Chao1 index and C/N for community composition). By contrast, the Chao1 index of soil bacterial communities was not affected by any soil physicochemical properties in SS, whereas the composition of the SS communities was correlated with soil depth, DOC, and the DOC/IN ratio.
Figure 6
Figure 6. Associations between the α-diversity and composition of soil bacterial communities and various soil physicochemical properties across all study sites. Mantel tests examining connections between the α-diversity (A) and composition (B) of bacterial communities and various soil physicochemical properties at four study sites containing different plant species. DOC, dissolved organic carbon; EC, electrical conductivity; NH4+, ammonium; NO3−, nitrate; SOC, soil organic carbon; SWC, soil water content; TN, total nitrogen; TP, total phosphorus; TC, Tamarix chinensis; SA, Spartina alterniflora; PA, Phragmites australis; SS, Suaeda salsa.
4 Discussion
4.1 Microbial community composition and diversity across study sites
The α-diversity of the soil bacterial community varied among study sites (Figure 2), with SA having higher α-diversity than the other three sites. Spartina alterniflora, an invasive plant with enhanced root exudates and litter inputs, greatly improved nutrient availability in the soil and thus supported greater bacterial α-diversity, in addition to affecting the composition and structure of the soil bacterial community (Song et al., 2020; Lin et al., 2022). Consistent with these results, Acidobacteria and Actinobacteria, widely recognized as r-strategists, were significantly more abundant than K-strategists (e.g., Bacteroidetes and Proteobacteria) in SA, supporting the finding that S. alterniflora affected the soil bacterial community by improving the nutrient content (Supplementary Table S2; Supplementary Figure S1). In addition, Mantel tests revealed that bacterial α-diversity in SA was correlated with the inorganic nitrogen content, which may be due to the fact that S. alterniflora significantly increased carbon resources in the soil (Figure 6A; Supplementary Table S1), leading to N-limitations for soil microorganisms.
In addition, substantial differences in bacterial community composition were observed among sites. Key bacterial genera involved in C-, N-, and S-cycling in coastal wetlands, including Desulfuromonas, Marinobacterium, Nitrospira, Subgroup_10, and Woeseia, were abundant at all four sites. However, the bacterial community composition varied across sites, especially for functions related to the N- and S-cycles. Using Mantel tests, differences in environmental drivers of community composition were found among sites. This suggests that variation in soil bacterial community composition was associated with site differences, including in local soil nutrient pools.
Among the bacterial taxa associated with N-cycling, the relative abundance of Calorithrix was relatively low in PA, while that of Nitrospira was elevated compared to other sites (Figures 3B,C). Calorithrix is associated with denitrification, while Nitrospira primarily participates in nitrification (Zhang X. et al., 2021; Daims and Wagner, 2018). Nitrification is regarded as a pivotal process in the production of soil N2O, potentially representing the primary mechanism (Lan et al., 2014). Here, the C/N ratio was the lowest among sites at PA, measuring 6.79: 1 (Supplementary Table S1), which is much lower than the average soil C/N ratio (10.1: 1) for coastal wetlands in China (Xia et al., 2023). Mantel test findings have shown that the soil bacterial composition at PA was correlated with the C/N ratio (Figure 6B). Sequently, Excess N may be removed by enhancing Nitrospira-based nitrification to reduce the adverse effects of a low C/N ratio (e.g., such as those caused by shifts in Calorithrix and Nitrospira abundance at PA). These results were confirmed by correlation test, which showed that higher TN and N/P and lower C/N at PA increased relative abundance of Nitrospira (Supplementary Figure S2).
Considering S-cycling, Desulfurivibrio and Desulfuromonas were most abundant in PA, whereas SEEP-SRB1 was most common in SA (Figure 3B). Desulfurivibrio, Desulfuromonas, and SEEP-SRB1 belong to the sulfate-reducing bacteria (SRB) (Petro et al., 2019; Yu et al., 2023) and are widely distributed in coastal wetlands, playing important roles in the S-cycle and participating in the release of wetland biosulfur gases and thereby mitigating sulfate stress in plants (Cui et al., 2017). In addition, sulfate can be reduced to sulfide via assimilatory sulfate reduction (ASR), an energy-consuming process, or it can be used as an electron acceptor in dissimilatory sulfate reduction (DSR) to generate energy (Yu et al., 2023). Previous studies have shown that Desulfovibrio and SEEP-SRB1 participate in DSR and ASR, respectively (Yu et al., 2023; Murali et al., 2023). In our study, the enrichment of these three genera was primarily driven by soil nitrogen content. For instance, the high abundance of Desulfurivibrio and Desulfuromonas in PA sites was associated with elevated TN and NO₃− levels, as well as higher N/P and lower C/N. Similarly, the high abundance of SEEP-SRB1 was linked to elevated NO₃− and reduced NH₄+ levels (Supplementary Figure S2). Interestingly, Sulfurimonas were also enriched in SA. Sulfurimonas are deep-sea chemoautotrophic bacteria capable of oxidizing sulfides (Wang et al., 2020). The co-occurrence of both sulfur-reducing and sulfur-oxidizing bacteria in SA is likely linked to the plant’s sulfur demands. Frequent tidal inundation in the low tidal zone, where S. alterniflora is located, promotes the accumulation of sulfur compounds, including acid-volatile sulfides, thereby supporting a higher abundance of sulfur-oxidizing bacteria such as Sulfurimonas (Li M. et al., 2021; Zheng et al., 2017). In contrast, other vegetation types like S. salsa and T. chinensis, located farther from the sea, experience less frequent tidal flooding and lower sulfur accumulation, which results in a lower abundance of sulfur-oxidizing bacteria. This dynamic balance between sulfur-reducing and sulfur-oxidizing bacteria in SA likely plays a key role in maintaining the sulfur cycle in its rhizosphere. Following S. alterniflora invasion and native plant succession, both plant species diversity and soil physicochemical properties were altered (Supplementary Table S1), resulting in changes in bacterial community composition among sites that were related to those in plant species composition.
4.2 Assembly of the bacterial community at different study sites
The null model results showed that deterministic and stochastic processes jointly drove bacterial community assembly at the four study sites, but stochastic processes dominated (Figure 4). The significance of stochastic processes for bacterial community assembly has been previously illustrated (Huo et al., 2023; Huang et al., 2022). Smaller organisms are more plastic in their metabolic abilities and have broader environmental tolerances; thus, bacterial communities are more affected by stochastic processes, according to the “size-plasticity” hypothesis (Wu et al., 2018; Zhu L. et al., 2024). However, it is worth noting that some studies have come to the opposite conclusion, showing that bacteria tend to have faster generation times and higher dispersal rates, which meaning they can respond to environmental changes more quickly and may be more susceptible to deterministic processes such as environmental filtering and/or biological interactions (Liu et al., 2023; Huo et al., 2023). The importance of deterministic versus stochastic processes may also depend on the specific microbial community and ecosystem being considered.
A conceptual model has been proposed to link the alternative hypotheses of ecosystem development to ecological assembly processes (Dini-Andreote et al., 2015). In this model, bacterial community composition is initially controlled by stochastic processes, while the role of deterministic selection increases with succession. As a relatively young wetland, the YRD continues to expand, forming new land (Qin et al., 2010). Relatively-infertile soils, repeated tidal action, and the lack of strong environmental filters allow stochastic processes to dominate microbial community assembly under these conditions. Environmental factors can then explain only a small fraction of the variation in bacterial community composition due to ongoing drift and dispersal (Fillinger et al., 2019; Danczak et al., 2018).
The strength of microbial responses to plant invasion is closely related to plant species identity, as well as to the duration and intensity of the invasion (Dostál et al., 2013; Zhang et al., 2020). Here, the amount of variation explained by stochastic processes was higher for SA versus the other sites, consistent with the findings of a previous study (Lin et al., 2022). The dominance of stochastic processes suggests that, when bacterial taxa have similar abilities to compete for common resources, taxa may coexist in highly overlapping ecological niches; this supports higher bacterial diversity as individual species are not eliminated by competition (Chase and Myers, 2011; Jiao et al., 2020). The assembly of high-diversity bacterial communities is mainly mediated by stochastic processes, while deterministic processes are more important when bacterial diversity is low, likely due to the reduction in specialized functions correlated with specific bacterial taxa (Xun et al., 2019). At SA, greater bacterial abundance and diversity enhanced the environmental tolerance and metabolic plasticity of the bacterial community, thus reducing the influence of environmental factors and increasing the importance of stochastic processes (Chen and Wen, 2021). The importance of stochasticity in bacterial community assembly may also increase with resource availability (Chase, 2010; Zheng et al., 2021). Thus, the higher nutrient availability and bacterial diversity associated with S. alterniflora invasion might favor stochastic processes in bacterial community assembly in coastal wetland soils.
4.3 Microbial co-occurrence networks for different study sites
Network analysis can be applied to identify species co-occurrence patterns and potential interactions in complex biotic communities (Faust and Raes, 2012; Zhu L. et al., 2024). Here, the bacterial co-occurrence network for SA had more nodes and edges, higher modularity values, and shorter average path lengths than the networks created for the other sites (Table 1), indicating greater complexity. This is because high modularity values indicate that while the nodes within a module are highly connected, effects on the rest of the network are reduced; this helps to maintain a community of complex interactions (Qiao et al., 2024). In addition, species in a module may be ecologically similar with overlapping niches (Huo et al., 2023). Niche overlap induces complex interactions in the community, explaining why the co-occurrence network of SA had more nodes and edges. The average path length is defined as the average number of steps along the shortest path between each pair of nodes, and represents a measure of network efficiency (Zhou et al., 2010). The low average path length in the SA network indicates faster information transmission within the network, allowing the microbial community to respond more quickly to environmental change (Zhou et al., 2010). S. alterniflora, located in the low tidal zone, is subject to frequent tidal inundation, which results in rapid transitions between aerobic and anaerobic conditions in the soil (Gao et al., 2018; Wei et al., 2020). Additionally, the high levels of sulfur accumulation associated with S. alterniflora invasion may cause stress to the soil bacterial community (Li M. et al., 2021). Consequently, faster information transmission within bacterial networks could enhance soil bacterial survival, enabling the bacterial community to better cope with stressful environmental conditions.
Positive interactions between microorganisms include cooperation, commensalism, and mutualism, while negative interactions include competition, parasitism, and predation (Zhu P. et al., 2024). The larger proportion of positive versus negative correlations in all four co-occurrence networks indicates the presence of cooperative or syntrophic interactions, suggesting that mutualistic interactions were common among soil bacterial taxa (Gao et al., 2021; Li C. et al., 2021). The coastal wetlands of YRD are characterized by relatively infertile soils and frequent tidal action, creating stressful conditions that promote mutualistic interactions within soil bacterial communities (Hernandez et al., 2021). Nevertheless, the PA network exhibited a greater number of negative interactions than the other three sites, suggesting a more competitive microbial environment. The supratidal zone and late successional plant species in PA contributed to improvements in soil nutrient availability, as evidenced by the relatively high SOC, TN, and TP, as well as the low C/N ratio. These factors likely contributed to the development of more competitive associations within the bacterial community.
5 Conclusion
In the YRD coastal wetlands, plant species identity was more influential than soil depth for shaping the soil bacterial community. However, it should be noted that other environmental factors, such as tidal inundation, also significantly contribute to microbial community composition. Spartina alterniflora invasion strongly impacted the soil microbial community, altering its structure and increasing bacterial α-diversity. Additionally stochastic processes largely drove bacterial community assembly, with S. alterniflora invasion further enhancing their dominance. The co-occurrence networks at the study sites primarily featured positive interactions, with S. alterniflora exhibiting higher complexity and modularity. These findings underscore the importance of plant communities in shaping soil bacterial diversity while also considering the influence of environmental factors like tidal flooding.
Data availability statement
The 16S rRNA sequencing data presented in this study are deposited in the NCBI SRA repository, BioProject accession number PRJNA1217841.
Author contributions
PS: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing. YW: Conceptualization, Investigation, Methodology, Visualization, Writing – original draft. PZ: Methodology, Resources, Visualization, Writing – review & editing. JW: Investigation, Visualization, Writing – review & editing. XY: Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing. WG: Funding acquisition, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (No. 32401428, No. U22A20558) and Key Laboratory of Ecological Prewarning, Protection and Restoration of Bohai Sea, Ministry of Natural Resources (No. 2022102).
Acknowledgments
We thank Dayou Zhou from Shandong University for assistance with samples collection.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1525632/full#supplementary-material
Footnotes
References
Chase, J. M. (2010). Stochastic community assembly causes higher biodiversity in more productive environments. Science 328, 1388–1391. doi: 10.1126/science.1187820
Chase, J. M., and Myers, J. A. (2011). Disentangling the importance of ecological niches from stochastic processes across scales. Philos. Trans. R. Soc. Lond. B Biol. Sci. 366, 2351–2363. doi: 10.1098/rstb.2011.0063
Chen, W., Ren, K., Isabwe, A., Chen, H., Liu, M., and Yang, J. (2019). Stochastic processes shape microeukaryotic community assembly in a subtropical river across wet and dry seasons. Microbiome 7, 138–116. doi: 10.1186/s40168-019-0749-8
Chen, W., and Wen, D. (2021). Archaeal and bacterial communities assembly and co-occurrence networks in subtropical mangrove sediments under Spartina alterniflora invasion. Environ. Microbiome 16:10. doi: 10.1186/s40793-021-00377-y
Cline, L. C., and Zak, D. R. (2015). Soil microbial communities are shaped by plant-driven changes in resource availability during secondary succession. Ecology 96, 3374–3385. doi: 10.1890/15-0184.1
Cui, J., Chen, X., Nie, M., Fang, S., Tang, B., Quan, Z., et al. (2017). Effects of Spartina alterniflora invasion on the abundance, diversity, and community structure of sulfate reducing bacteria along a successional gradient of coastal salt marshes in China. Wetlands 37, 221–232. doi: 10.1007/s13157-016-0860-6
Daims, H., and Wagner, M. (2018). Nitrospira. Trends Microbiol. 26, 462–463. doi: 10.1016/j.tim.2018.02.001
Danczak, R. E., Johnston, M. D., Kenah, C., Slattery, M., and Wilkins, M. J. (2018). Microbial community cohesion mediates community turnover in unperturbed aquifers. mSystems 3, e00066–e00018. doi: 10.1128/msystems.00066-18
Dini-Andreote, F., Stegen, J. C., Van Elsas, J. D., and Salles, J. F. (2015). Disentangling mechanisms that mediate the balance between stochastic and deterministic processes in microbial succession. Proc. Natl. Acad. Sci. USA 112, E1326–E1332. doi: 10.1073/pnas.1414261112
Dostál, P., Müllerová, J., Pyšek, P., Pergl, J., and Klinerová, T. (2013). The impact of an invasive plant changes over time. Ecol. Lett. 16, 1277–1284. doi: 10.1111/ele.12166
Faust, K., and Raes, J. (2012). Microbial interactions: from networks to models. Nat. Rev. Microbiol. 10, 538–550. doi: 10.1038/nrmicro2832
Fillinger, L., Hug, K., and Griebler, C. (2019). Selection imposed by local environmental conditions drives differences in microbial community composition across geographically distinct groundwater aquifers. FEMS Microbiol. Ecol. 95:fiz160. doi: 10.1093/femsec/fiz160
Gao, H., Bai, J., Deng, X., Lu, Q., and Ye, X. (2018). Short-term effects of tidal flooding on soil nitrogen mineralization in a Chinese tidal salt marsh. Phys. Chem. Earth 103, 3–10. doi: 10.1016/j.pce.2017.04.002
Gao, J. H., Feng, Z. X., Chen, L., Wang, Y. P., Bai, F., and Li, J. (2016). The effect of biomass variations of Spartina alterniflora on the organic carbon content and composition of a salt marsh in northern Jiangsu Province, China. Ecol. Eng. 95, 160–170. doi: 10.1016/j.ecoleng.2016.06.088
Gao, G. F., Peng, D., Zhang, Y., Li, Y., Fan, K., Tripathi, B. M., et al. (2021). Dramatic change of bacterial assembly process and co-occurrence pattern in Spartina alterniflora salt marsh along an inundation frequency gradient. Sci. Total Environ. 755:142546. doi: 10.1016/j.scitotenv.2020.142546
Hermans, S. M., Buckley, H. L., Case, B. S., Curran-Cournane, F., Taylor, M., and Lear, G. (2020). Using soil bacterial communities to predict physico-chemical variables and soil quality. Microbiome 8:79. doi: 10.1186/s40168-020-00858-1
Hernandez, D. J., David, A. S., Menges, E. S., Searcy, C. A., and Afkhami, M. E. (2021). Environmental stress destabilizes microbial networks. ISME J. 15, 1722–1734. doi: 10.1038/s41396-020-00882-x
Huang, L., Bai, J., Wang, J., Zhang, G., Wang, W., Wang, X., et al. (2022). Different stochastic processes regulate bacterial and fungal community assembly in estuarine wetland soils. Soil Biol. Biochem. 167:108586. doi: 10.1016/j.soilbio.2022.108586
Huo, X., Ren, C., Wang, D., Wu, R., Wang, Y., Li, Z., et al. (2023). Microbial community assembly and its influencing factors of secondary forests in Qinling Mountains. Soil Biol. Biochem. 184:109075. doi: 10.1016/j.soilbio.2023.109075
Hussain, S., Chen, M., Liu, Y., Mustafa, G., Wang, X., Liu, J., et al. (2023). Composition and assembly mechanisms of prokaryotic communities in wetlands, and their relationships with different vegetation and reclamation methods. Sci. Total Environ. 897:166190. doi: 10.1016/j.scitotenv.2023.166190
Jiang, S., Zhang, C., Chen, L., Liu, C., Yan, L., and Li, B. (2022). Effects of smooth cordgrass Spartina alterniflora invasion on macrobenthic fauna in the Yellow River Delta. Wetlands 42:13. doi: 10.1007/s13157-021-01527-1
Jiao, S., Yang, Y., Xu, Y., Zhang, J., and Lu, Y. (2020). Balance between community assembly processes mediates species coexistence in agricultural soil microbiomes across eastern China. ISME J. 14, 202–216. doi: 10.1038/s41396-019-0522-9
Lan, T., Han, Y., Roelcke, M., Nieder, R., and Car, Z. (2014). Sources of nitrous and nitric oxides in paddy soils: nitrification and denitrification. J. Environ. Sci. 26, 581–592. doi: 10.1016/s1001-0742(13)60453-2
Leewis, M. C., Lawrence, C. R., Schulz, M. S., Tfaily, M. M., Ayala-Ortiz, C. O., Flores, G. E., et al. (2022). The influence of soil development on the depth distribution and structure of soil microbial communities. Soil Biol. Biochem. 174:108808. doi: 10.1016/j.soilbio.2022.108808
Li, M., Fang, A., Yu, X., Zhang, K., He, Z., Wang, C., et al. (2021). Microbially-driven sulfur cycling microbial communities in different mangrove sediments. Chemosphere 273:128597. doi: 10.1016/j.chemosphere.2020.128597
Li, C., Wang, L., Ji, S., Chang, M., Wang, L., Gan, Y., et al. (2021). The ecology of the plastisphere: microbial composition, function, assembly, and network in the freshwater and seawater ecosystems. Water Res. 202:117428. doi: 10.1016/j.watres.2021.117428
Li, M., Zhou, W., Sun, M., Shi, W., Lun, J., Zhou, B., et al. (2024). Decoupling soil community structure, functional composition, and nitrogen metabolic activity driven by salinity in coastal wetlands. Soil Biol. Biochem. 198:109547. doi: 10.1016/j.soilbio.2024.109547
Lin, Y., Hu, H. W., Yang, P., and Ye, G. (2022). Spartina alterniflora invasion has a greater impact than non-native species, Phragmites australis and Kandelia obovata, on the bacterial community assemblages in an estuarine wetland. Sci. Total Environ. 822:153517. doi: 10.1016/j.scitotenv.2022.153517
Liu, L., Wu, Y., Yin, M., Ma, X., Yu, X., Guo, X., et al. (2023). Soil salinity, not plant genotype or geographical distance, shapes soil microbial community of a reed wetland at a fine scale in the Yellow River Delta. Sci. Total Environ. 856:159136. doi: 10.1016/j.scitotenv.2022.159136
Lu, K., Han, G., and Wu, H. (2022). Effects of Spartina alterniflora invasion on the benthic invertebrate community in intertidal wetlands. Ecosphere 13:e 3963. doi: 10.1002/ecs2.3963
Murali, R., Yu, H., Speth, D. R., Wu, F., Metcalfe, K. S., Crémière, A., et al. (2023). Physiological potential and evolutionary trajectories of syntrophic sulfate-reducing bacterial partners of anaerobic methanotrophic archaea. PLoS Biol. 21:e3002292. doi: 10.1371/journal.pbio.3002292
Nelson, D. W., and Sommers, L. E. (1982). Total carbon, organic carbon, and organic matter. Methods Soil Anal. 9, 539–579. doi: 10.2134/agronmonogr9.2.2ed.c29
Ni, Y., Yang, T., Ma, Y., Zhang, K., Soltis, P. S., Soltis, D. E., et al. (2021). Soil pH determines bacterial distribution and assembly processes in natural mountain forests of eastern China. Glob. Ecol. Biogeogr. 30, 2164–2177. doi: 10.1111/geb.13373
Nossa, C. W., Oberdorf, W. E., Yang, L., Aas, J. A., Paster, B. J., DeSantis, T. Z., et al. (2010). Design of 16S rRNA gene primers for 454 pyrosequencing of the human foregut microbiome. World J. Gastroenterol. 16, 4135–4144. doi: 10.3748/wjg.v16.i33.4135
Osland, M. J., Chivoiu, B., Enwright, N. M., Thorne, K. M., Guntenspergen, G. R., Grace, J. B., et al. (2022). Migration and transformation of coastal wetlands in response to rising seas. Sci. Adv. 8:eabo5174. doi: 10.1126/sciadv.abo5174
Petro, C., Jochum, L. M., Schreiber, L., Marshall, I. P., Schramm, A., and Kjeldsen, K. U. (2019). Single-cell amplified genomes of two uncultivated members of the deltaproteobacterial SEEP-SRB1 clade, isolated from marine sediment. Mar. Genomics 46, 66–69. doi: 10.1016/j.margen.2019.01.004
Philippot, L., Raaijmakers, J. M., Lemanceau, P., and Van Der Putten, W. H. (2013). Going back to the roots: the microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799. doi: 10.1038/nrmicro3109
Qiao, Y., Wang, T., Huang, Q., Guo, H., Zhang, H., Xu, Q., et al. (2024). Core species impact plant health by enhancing soil microbial cooperation and network complexity during community coalescence. Soil Biol. Biochem. 188:109231. doi: 10.1016/j.soilbio.2023.109231
Qin, Y., Yang, Z., and Yang, W. (2010). A novel index system for assessing ecological risk under water stress in the Yellow River delta wetland. Procedia Environ. Sci. 2, 535–541. doi: 10.1016/j.proenv.2010.10.058
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Ren, G. B., Wang, J. J., Wang, A. D., Wang, J. B., Zhu, Y. L., Wu, P. Q., et al. (2019). Monitoring the invasion of smooth cordgrass Spartina alterniflora within the modern Yellow River Delta using remote sensing. J. Coast. Res. 90, 135–145. doi: 10.2112/SI90-017.1
Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahé, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584
Sardans, J., Bartrons, M., Margalef, O., Gargallo-Garriga, A., Janssens, I. A., Ciais, P., et al. (2017). Plant invasion is associated with higher plant–soil nutrient concentrations in nutrient-poor environments. Glob. Chang. Biol. 23, 1282–1291. doi: 10.1111/gcb.13384
Shigyo, N., Umeki, K., and Hirao, T. (2019). Plant functional diversity and soil properties control elevational diversity gradients of soil bacteria. FEMS Microbiol. Ecol. 95:fiz025. doi: 10.1093/femsec/fiz025
Sloan, W. T., Lunn, M., Woodcock, S., Head, I. M., Nee, S., and Curtis, T. P. (2006). Quantifying the roles of immigration and chance in shaping prokaryote community structure. Environ. Microbiol. 8, 732–740. doi: 10.1111/j.1462-2920.2005.00956.x
Song, S., Zhang, C., Gao, Y., Zhu, X., Wang, R., Wang, M., et al. (2020). Responses of wetland soil bacterial community and edaphic factors to two-year experimental warming and Spartina alterniflora invasion in Chongming Island. J. Clean. Prod. 250:119502. doi: 10.1016/j.jclepro.2019.119502
Tan, L., Ge, Z., Zhou, X., Li, S., Li, X., and Tang, J. (2020). Conversion of coastal wetlands, riparian wetlands, and peatlands increases greenhouse gas emissions: a global meta-analysis. Glob. Chang. Biol. 26, 1638–1653. doi: 10.1111/gcb.14933
Tian, J., Liu, L., Chen, H., Zhong, L., Zhou, X., Jiang, L., et al. (2022). Aerobic environments in combination with substrate additions to soil significantly reshape depth-dependent microbial distribution patterns in Zoige peatlands, China. Appl. Soil Ecol. 170:104252. doi: 10.1016/j.apsoil.2021.104252
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, S., Jiang, L., Liu, X., Yang, S., and Shao, Z. (2020). Sulfurimonas xiamenensis sp. nov. and Sulfurimonas lithotrophica sp. nov., hydrogen-and sulfur-oxidizing chemolithoautotrophs within the Epsilonproteobacteria isolated from coastal sediments, and an emended description of the genus Sulfurimonas. Int. J. Syst. Evol. Microbiol. 70, 2657–2663. doi: 10.1099/ijsem.0.004087
Wang, W., Sardans, J., Wang, C., Zeng, C., Tong, C., Chen, G., et al. (2019). The response of stocks of C, N, and P to plant invasion in the coastal wetlands of China. Glob. Chang. Biol. 25, 733–743. doi: 10.1111/gcb.14491
Ward, S. E., Orwin, K. H., Ostle, N. J., Briones, M. J., Thomson, B. C., Griffiths, R. I., et al. (2015). Vegetation exerts a greater control on litter decomposition than climate warming in peatlands. Ecology 96, 113–123. doi: 10.1890/14-0292.1
Webb, C. O., Ackerly, D. D., McPeek, M. A., and Donoghue, M. J. (2002). Phylogenies and community ecology. Annu. Rev. Ecol. Syst. 33, 475–505. doi: 10.1146/annurev.ecolsys.33.010802.150448
Wei, S., Han, G., Chu, X., Song, W., He, W., Xia, J., et al. (2020). Effect of tidal flooding on ecosystem CO2 and CH4 fluxes in a salt marsh in the Yellow River Delta. Estuar. Coast. Shelf Sci. 232:106512. doi: 10.1016/j.ecss.2019.106512
Wu, W., Lu, H. P., Sastri, A., Yeh, Y. C., Gong, G. C., Chou, W. C., et al. (2018). Contrasting the relative importance of species sorting and dispersal limitation in shaping marine bacterial versus protist communities. ISME J. 12, 485–494. doi: 10.1038/ismej.2017.183
Xia, S., Song, Z., Singh, B. P., Guo, L., Bolan, N., Wang, W., et al. (2023). Contrasting patterns and controls of soil carbon and nitrogen isotope compositions in coastal wetlands of China. Plant Soil 489, 483–505. doi: 10.1007/s11104-023-06034-2
Xun, W., Li, W., Xiong, W., Ren, Y., Liu, Y., Miao, Y., et al. (2019). Diversity-triggered deterministic bacterial assembly constrains community functions. Nat. Commun. 10:3833. doi: 10.1038/s41467-019-11787-5
Yu, X., Hu, R., Tao, M., Qian, L., Wang, F., Wang, S., et al. (2023). Microbially driven sulfur cycling in the river–wetland–ocean continuum. Ocean Land Atmos. Res. 2:0027. doi: 10.34133/olar.0027
Yu, X., Yang, J., Liu, L., Tian, Y., and Yu, Z. (2015). Effects of Spartina alterniflora invasion on biogenic elements in a subtropical coastal mangrove wetland. Environ. Sci. Pollut. Res. 22, 3107–3115. doi: 10.1007/s11356-014-3568-2
Yue, Y., Tang, Y., Cai, L., Yang, Z., Chen, X., Ouyang, Y., et al. (2022). Co-occurrence relationship and stochastic processes affect sedimentary archaeal and bacterial community assembly in estuarine–coastal margins. Microorganisms 10:1339. doi: 10.3390/microorganisms10071339
Zhalnina, K., Louie, K. B., Hao, Z., Mansoori, N., Da Rocha, U. N., Shi, S., et al. (2018). Dynamic root exudate chemistry and microbial substrate preferences drive patterns in rhizosphere microbial community assembly. Nat. Microbiol. 3, 470–480. doi: 10.1038/s41564-018-0129-3
Zhang, G., Bai, J., Jia, J., Wang, W., Wang, X., Zhao, Q., et al. (2019). Shifts of soil microbial community composition along a short-term invasion chronosequence of Spartina alterniflora in a Chinese estuary. Sci. Total Environ. 657, 222–233. doi: 10.1016/j.scitotenv.2018.12.061
Zhang, G., Bai, J., Tebbe, C. C., Huang, L., Jia, J., Wang, W., et al. (2021). Spartina alterniflora invasions reduce soil fungal diversity and simplify co-occurrence networks in a salt marsh ecosystem. Sci. Total Environ. 758:143667. doi: 10.1016/j.scitotenv.2020.143667
Zhang, G., Bai, J., Tebbe, C. C., Huang, L., Jia, J., Wang, W., et al. (2022). Plant invasion reconstructs soil microbial assembly and functionality in coastal salt marshes. Mol. Ecol. 31, 4478–4494. doi: 10.1111/mec.16600
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
Zhang, G., Jia, J., Zhao, Q., Wang, W., Wang, D., and Bai, J. (2023). Seasonality and assembly of soil microbial communities in coastal salt marshes invaded by a perennial grass. J. Environ. Manag. 331:117247. doi: 10.1016/j.jenvman.2023.117247
Zhang, G. S., Wang, R. Q., and Song, B. M. (2007). Plant community succession in modern Yellow River Delta, China. J Zhejiang Univ Sci B 8, 540–548. doi: 10.1631/jzus.2007.B0540
Zhang, X., Zhang, Z., Li, Z., Li, M., Wu, H., and Jiang, M. (2021). Impacts of Spartina alterniflora invasion on soil carbon contents and stability in the Yellow River Delta, China. Sci. Total Environ. 775:145188. doi: 10.1016/j.scitotenv.2021.145188
Zheng, Y., Bu, N. S., Long, X. E., Sun, J., He, C. Q., Liu, X. Y., et al. (2017). Sulfate reducer and sulfur oxidizer respond differentially to the invasion of Spartina alterniflora in estuarine salt marsh of China. Ecol. Eng. 99, 182–190. doi: 10.1016/j.ecoleng.2016.11.031
Zheng, W., Zhao, Z., Lv, F., Wang, R., Wang, Z., Zhao, Z., et al. (2021). Assembly of abundant and rare bacterial and fungal sub-communities in different soil aggregate sizes in an apple orchard treated with cover crop and fertilizer. Soil Biol. Biochem. 156:108222. doi: 10.1016/j.soilbio.2021.108222
Zhou, J., Deng, Y., Luo, F., He, Z., Tu, Q., and Zhi, X. (2010). Functional molecular ecological networks. mBio 1, 10–1128. doi: 10.1128/mbio.00169-10
Zhou, J., Deng, Y. E., Zhang, P., Xue, K., Liang, Y., Van Nostrand, J. D., et al. (2014). Stochasticity, succession, and environmental perturbations in a fluidic ecosystem. Proc. Natl. Acad. Sci. 111, E836–E845. doi: 10.1073/pnas.1324044111
Zhou, J., and Ning, D. (2017). Stochastic community assembly: does it matter in microbial ecology? Microbiol. Mol. Biol. Rev. 81, 10–1128. doi: 10.1128/mmbr.00002-17
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
Zhu, P., Liu, C., Wei, W., Wu, Y., Sardar, M. F., Yu, X., et al. (2024). Stochastic processes limit the effect of organic fertilizer application on soil bacterial community composition in salt marsh Suaeda salsa. J. Clean. Prod. 441:141034. doi: 10.1016/j.jclepro.2024.141034
Keywords: assembly process, bacterial community, coastal wetlands, co-occurrence network, Spartina alterniflora invasion
Citation: Sun P, Wu Y, Zhu P, Wang J, Yu X and Guo W (2025) Spartina alterniflora invasion significantly alters the assembly and structure of soil bacterial communities in the Yellow River Delta. Front. Microbiol. 16:1525632. doi: 10.3389/fmicb.2025.1525632
Edited by:
Francesca Mapelli, University of Milan, ItalyReviewed by:
Abhinav Aeron, Chonbuk National University, Republic of KoreaNagaraju Yalavarthi, Central Silk Board, India
Milan Farkas, Hungarian University of Agricultural and Life Sciences, Hungary
Copyright © 2025 Sun, Wu, Zhu, Wang, Yu and Guo. 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: Weihua Guo, Z3Vvd2hAZW1haWwuc2R1LmVkdS5jbg==; Xiaona Yu, ZWNvX3NkdUAxMjYuY29t