Next Article in Journal
Inhibitory Effect of Lactococcus and Enterococcus faecalis on Citrobacter Colitis in Mice
Next Article in Special Issue
Real-Time PCR (qtPCR) to Discover the Fate of Plant Growth-Promoting Rhizobacteria (PGPR) in Agricultural Soils
Previous Article in Journal
Repeated Occurrence of Mobile Colistin Resistance Gene-Carrying Plasmids in Pathogenic Escherichia coli from German Pig Farms
Previous Article in Special Issue
The Bacterial and Fungal Compositions in the Rhizosphere of Asarum heterotropoides Fr. Schmidt var. mandshuricum (Maxim.) Kitag. in a Typical Planting Region
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 12
Issue 4
10.3390/microorganisms12040728
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association between Soil Physicochemical Properties and Bacterial Community Structure in Diverse Forest Ecosystems

by
Bing Yang
1,*,
Wanju Feng
1,
Wenjia Zhou
1,
Ke He
2 and
Zhisong Yang
1
1
Sichuan Academy of Giant Panda, Chengdu 610041, China
2
Key Laboratory of Southwest China Wildlife Resources Conservation (Ministry of Education), China West Normal University, Nanchong 637002, China
*
Author to whom correspondence should be addressed.
Microorganisms 2024, 12(4), 728; https://doi.org/10.3390/microorganisms12040728
Submission received: 18 March 2024 / Revised: 29 March 2024 / Accepted: 30 March 2024 / Published: 3 April 2024
(This article belongs to the Special Issue Microbial Ecology of Soil Microorganisms, Networks with Plants and Ecosystem Services)
Browse Figures
Review Reports Versions Notes

Abstract

:
Although the importance of the soil bacterial community for ecosystem functions has long been recognized, there is still a limited understanding of the associations between its community composition, structure, co-occurrence patterns, and soil physicochemical properties. The objectives of the present study were to explore the association between soil physicochemical properties and the composition, diversity, co-occurrence network topological features, and assembly mechanisms of the soil bacterial community. Four typical forest types from Liziping Nature Reserve, representing evergreen coniferous forest, deciduous coniferous forest, mixed conifer-broadleaf forest, and its secondary forest, were selected for this study. The soil bacterial community was analyzed using Illumina MiSeq sequencing of 16S rRNA genes. Nonmetric multidimensional scaling was used to illustrate the clustering of different samples based on Bray–Curtis distances. The associations between soil physicochemical properties and bacterial community structure were analyzed using the Mantel test. The interactions among bacterial taxa were visualized with a co-occurrence network, and the community assembly processes were quantified using the Beta Nearest Taxon Index (Beta-NTI). The dominant bacterial phyla across all forest soils were Proteobacteria (45.17%), Acidobacteria (21.73%), Actinobacteria (8.75%), and Chloroflexi (5.06%). Chao1 estimator of richness, observed ASVs, faith-phylogenetic diversity (faith-PD) index, and community composition were distinguishing features of the examined four forest types. The first two principal components of redundancy analysis explained 41.33% of the variation in the soil bacterial community, with total soil organic carbon, soil moisture, pH, total nitrogen, carbon/nitrogen (C/N), carbon/phosphorous (C/P), and nitrogen/phosphorous (N/P) being the main soil physicochemical properties shaping soil bacterial communities. The co-occurrence network structure in the mixed forest was more complex compared to that in pure forests. The Beta-NTI indicated that the bacterial community assembly of the four examined forest types was collaboratively influenced by deterministic and stochastic ecological processes.

1. Introduction

Forest ecosystems play a vital role in supporting global biodiversity and contributing to various ecosystem services that are essential for human well-being and the health of our planet. Bacteria, the most abundant and diverse microorganisms in the soil [1,2,3,4], participate in ecological processes and provide ecosystem services [1,5,6]. Bacterial communities in forest soils, driven by various ecological interactions and feedback loops, contribute significantly to nutrient cycling, organic matter decomposition, and plant health. Additionally, soil bacterial communities are sensitive to environmental changes. Numerous studies have shown that soil bacterial communities on a local scale change with vegetation types and succession stages [7,8] and soil prosperities, such as soil water content (SWC) [9,10], pH [9,11,12,13], soil organic carbon (SOC) [9,14], total phosphorous (TP), total nitrogen (TN), and available phosphorous (AP) concentrations [15], and C/N ratio [16]. However, due to high forest soil heterogeneity and complicated plant-soil interactions, the environmental factors driving the abundance, diversity, and composition of the soil bacterial community remain unclear. This greatly hampers our predictions about the responses of the soil bacterial community to global climate change. Under the Paris agreement on climate change reached in 2015, China promised to halt the rise of its carbon emissions by 2030 and pledge “carbon neutrality” by 2060. Since the soil bacterial community plays a crucial role in organic matter mineralization, nutrient cycling, plant production, and greenhouse gas emissions [17,18], a deep understanding of the associations among forest type, soil bacterial community, and related functions is of utmost importance.
Over the past few decades, advances in high-throughput sequencing technologies and computational methods have revolutionized the study of soil microbial communities, allowing researchers to delve deeper into the complex interactions and dynamics that govern their assembly. High-throughput sequencing of 16S rDNA in environmental samples enables assessment of bacterial beta diversity and determinants of its biogeography [12] and has been applied to assess the diversity and functions of bacterial communities in forests [4,19,20,21,22,23]. The development of statistical methods makes it possible to gain a comprehensive understanding of the responses of the soil bacterial community to environmental change. For instance, the co-occurrence network analysis can offer new insights into the putative interactions between taxa in complex microbial communities [9,20,24,25]. Besides, it has proven powerful in identifying keystone taxa that have the largest influence in a community [20,25,26,27]. To date, co-occurrence analysis has been widely applied to exploring the responses of the microbial community to environmental changes in forest soils [20,21,23,28]. Additionally, community assembly offers a framework to explore the mechanisms driving microbial community dynamics. Community assembly theories, such as the niche theory [29] and the neutral theory [30], have provided conceptual frameworks to elucidate the underlying mechanisms driving the composition and structure of soil microbial communities [31,32]. Various studies have suggested that soil bacterial community assembly is combined and driven by deterministic and stochastic processes. However, there is no consensus yet on the relative importance of deterministic and stochastic processes in shaping soil bacterial community assembly, particularly in forest ecosystems.
The original intention of establishing Giant Panda National Park was to preserve the giant panda (Ailuropoda melanoleuca) and its habitat and to restore and enhance ecosystem multifunctionality. Studies have demonstrated that giant pandas’ habit are found in primary forests [33] and in secondary forests [34], whereas giant pandas are more abundant in traces in primary forests in comparison to secondary forests [35]. The soil bacterial community is a critical participant in the ecological process and thus an important resource to utilize in ecological restoration. However, its diversity and driving force in giant panda habitat remain poorly understood, greatly restricting our manipulation of the microbiome to improve soil health and enhance ecosystem multifunctionality.
The objectives of the present study are to explore the association between forest stand-type soil physicochemical properties and the composition, diversity, co-occurrence network topological features, and assembly mechanisms of the soil bacterial community. We hypothesize that (1) forest type significantly influences composition, diversity, co-occurrence topological features, and community assembly; (2) soil properties directly shape soil bacterial community features and co-occurrence pattern; and (3) ecological community assembly processes of soil bacterial communities are dominated by stochastic processes.

2. Material and Method

2.1. Site Description

The Liziping National Nature Reserve is a core distribution area and critical habitat for the giant panda. This reserve encompasses 47,885 hm2 of rugged ridges and narrow valleys at elevations ranging from 1330 to 4550 m. As reported, the mean annual temperature and precipitation are 11.7–14.4 °C and 800–1250 mm, respectively. The soil type of this area is braunerde, or brown forest soil.

2.2. Soil Sampling

Soil samples were acquired in October 2020 from the four selected forest types, namely coniferous (Con) forest, evergreen broadleaf (Bro) forest, secondary (Sec) forest, natural deciduous broadleaved and evergreen broadleaved mixed (Mix) forest. Within each forest type, nine sampling plots were established; dominant tree species, canopy cover, and dominant understory species were considered among the three plots to avoid unrepresentative areas. Within each plot, three 5 m × 5 m sized quadrates were established, and soil samples (0–20 cm) were collected using a 5 cm diameter soil auger from 10 locations along an S-shaped path after removing superficial debris (e.g., leaves, dry vegetation, and litter), mixed uniformly, and pooled to ensure that the soil samples were representative. The soil samples were then transferred to sterilized press seal bags and stored in an icebox. Next, the samples were taken to the laboratory and immediately sieved to remove stones and plant materials using a 2 mm mesh. One portion of the samples was air-dried and used to determine soil physicochemical properties, whereas the remainder was kept at −80 °C before soil bacterial community analysis.

2.3. Soil DNA Isolation and PCR Amplification

DNA extraction was carried out on 0.5 g frozen bulk soil samples using the E.Z.N.A.® soil DNA Kit (OMEGA Bio-Tek, Norcross, GA, USA) following the provided instructions. Subsequently, the concentration and purity of the extracted DNA were determined through the utilization of a NanoDrop 2000 UV-Vis spectrophotometer (Thermo Scientific, Wilmington, NC, USA) and 1% agarose gel. The target sequence for analysis was the 16s rRNA V3-V4 region of bacteria, for which PCR amplification was conducted using specific primer pairs 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) on an ABI GeneAmp® 9700 PCR thermocycler (ABI, Los Angeles, CA, USA). The PCR mixtures were meticulously prepared, comprising various components such as high-fidelity DNA polymerase, reaction buffer, dNTPs, DNA templates, and specific primers. The amplification process involved a series of cycles with specified temperature and duration parameters, ultimately leading to the desired amplification of the target sequence in the extracted DNA. The PCR mixtures were prepared using 0.25 μL of Q5 high-fidelity DNA polymerase, 5 μL of reaction buffer (5×), 5 μL of high GC buffer (5×), 2 μL of dNTPs (10 mM), 2 μL of DNA templates, 1 μL each of forward (10 μm) and reverse (10 μm) primers, and 8.75 μL of ddH2O PCR. The specified amplification conditions, including denaturation, annealing, and extension steps, were diligently followed to facilitate the amplification of the target sequence. The amplification conditions were an initial denaturation at 95 °C for 3 min, followed by 27 cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s, and extension at 72 °C for 45 s, with a final extension at 72 °C for 10 min.

2.4. Illumina Miseq Sequencing

The PCR amplification product was extracted using a 2% agarose gel and purified using an AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA). Purified amplicons were pooled in equimolar amounts and paired-end sequenced on an Illumina Miseq PE 300 platform (Illumina, San Diego, CA, USA) according to the standard protocols of Shanghai Personalbio Technology Co., Ltd. (Shanghai, China). The raw reads for all samples were sent to the NCBI Sequence Read Archive (SRA) database (accession number: PRJNA1005035).

2.5. Bioinformatic Analyses

The analysis of raw sequences was carried out using the QIIME2 2019.4 platform [36], where the sequences were trimmed, quality filtered, denoised, merged, and chimeras removed with the DADA2 pipeline [37]. Subsequently, sequences with a 100% match were grouped into amplicon sequence variants (ASVs) [37,38]. ASVs that were not singleton were aligned using MAFFT [39], and a phylogenetic tree was constructed using fasttree2 [40]. Following this, diversity metrics such as observed species, Chao1, Shannon, Simpson, and Faith’s PD were computed using QIIME2 [41].

2.6. Soil Physicochemical Parameters Measurements

Soil physicochemical properties were assayed with a standardized method. Specifically, SWC was determined with ten grams of field-moist soil samples at 105 °C for 24 h. Soil pH was measured using a glass electrode pH meter (the soil-to-water ratio was 1:2.5). The SOC content was measured using the high-temperature external heat dichromate oxidation capacity method. TN and AN were determined using the Kjeldahl digestion method and the alkaline hydrolysis method, respectively. TP was measured using the Mo-Sb anti-spectrophotometric method. Available phosphorus was extracted using an HCl-NH4F solution and determined using the molybdenum-antimony resistance colorimetric method. Soil C/N, C/P, and N/P were calculated based on SOC, TN, and TP.

2.7. Statistical Analyses

The effects of forest type on soil physicochemical parameters and bacterial community diversity were examined using one-way ANOVAs or non-parametric tests, depending on the normality and homogeneity of the variance. If a significant effect was detected, the differences were further assessed using the TukeyHSD test or pairwise. t test with the “Benjaminiand-Hochberg” method to adjust the p value. Spearman’s correlation analysis was utilized to assess the relationship between bacterial alpha diversity and soil physicochemical properties in R studio.
The effects of forest type on soil bacterial community composition and structure were examined using the following protocol: initially, Venn diagrams were generated to visualize the unique and shared amplicon sequence variants (ASVs) among the samples in R v4.2.1 using the “Venn Diagram” package. Additionally, heatmaps of the top 50 genera per sample were created in R using the “ggplot” and “pheatmap” packages. Furthermore, linear discriminant analysis (LDA) effect size (LEfSe) analysis was conducted to identify bacterial indicator taxa based on a normalized relative abundance matrix across groups using default parameters. Nonmetric multidimensional scaling (NMDS) based on Bray–Curtis matrices was performed to characterize the composition of the soil bacterial community using the “vegan” package. Analysis of similarities (ANOSIM) and permutation multivariate analysis of variance (PERMANOVA) were carried out to assess the statistically significant differences among bacterial communities using Bray–Curtis distances and 999 permutations. The Bray–Curtis distance was chosen due to its popularity as a similarity index for abundance data [16]. Redundancy analysis (RDA) based on Hellinger distance, a commonly used method for quantifying the similarity between two probability distributions [16], and the Mantel test with 999 permutations were employed to determine the effects of soil properties on the bacterial communities using the “vegan” and “ade4” packages in R studio.
The co-occurrence patterns of soil bacterial communities among the different forest types were determined using network analysis based on correlations of the 100 most abundant ASVs in the soil bacterial community. Co-occurring networks based on Spearman correlation analysis in this study were conducted using the “psych” package in R studio. The co-occurrence patterns of soil bacterial communities were studied based on strong correlations (r > 0.6) and significant correlations (p < 0.01). Network analysis and visualization were performed using igraph “package” [42].
The assemblage mechanisms of soil bacterial communities were analyzed through combined quantifying niche width, Sloan neutral community model [43] and null model [44]. Niche width was measured according to Levins’ coefficient:
B i = 1 ∑ j = 1 r P i j 2
where Bi is the habitat niche width of ASVi, and Pij is the proportion of ASVi in the total ASVs within a given forest type j. The average number of all ASVs’ B was calculated to represent the niche width of the soil bacterial community. The neutral community model was fitted by the nonlinear least-square fitting method, and the 95% confidence interval was predicted by the “Hmisc” package [45]. In terms of the null model, Beta Taxon Index (βNTI) and Raup–Crick (RC-Bray) were calculated to represent phylogenetic and taxonomic diversity [45,46,47]. |βNTI| > 2 indicates the dominance of deterministic processes, while |βNTI| < 2 indicates the dominance of stochastic processes. |βNTI| < 2 and RCBray < 0.95 represent homogenizing dispersal. |βNTI| < 2 and RCBray > 0.95 represent dispersal limitations. |βNTI| < 2 and |RCBray| < 0.95 represent “undominated” assembly (mainly consists of weak selection, weak dispersal, diversification, and/or drift). βNTI < 2 represents homogeneous selection. βNTI > 2 represents variable selection.

3. Results

3.1. Soil Physicochemical Properties

The soil moisture, SOC, TN, AN, TP, C:P, and N:P ratios varied greatly depending on forest type (Figure 1). The Bro Forest had the highest SWC, TN, AN, C/P, and N/P. The Con Forest had the lowest SWC, TN, AN, C/P, and N/P, but the highest TP.

3.2. Composition and Diversity

The full dataset contained 3,924,845 sequences from 36 successfully amplified soil samples. After processing within Qiime 2, a total of 3,288,127 sequences and 31,340 ASVs were retained for further analysis. According to the Venn diagram of the soil bacterial community, the total number of ASVs of soil bacteria shared among the four forest types was 2222 (Figure 2), accounting for approximately 7.09% of all ASVs. The Sec Forest exhibited the largest number of unique ASVs, with 6428 specific ASVs accounting for 20.51% of the total ASVs. The largest number of shared soil bacterial ASVs was found in the pair of Sec and Con forests. There were 1395 shared soil bacterial ASVs between these two forests, accounting for approximately 4.45% of the total ASVs. The dominant bacterial phyla across all forest soil samples were Proteobacteria (relative abundance 45.17%), Acidobacteria (21.73%), Actinobacteria (8.75%), and Chloroflexi (5.06%), whereas Bacteroidetes, Planctomycetes, Verrucomicrobia, Thaumarchaeota, Gemmatimonadetes, Rokubacteria, and other bacterial species occupied only a minor fraction of the bacterial community composition in bulk soils at the topsoil layer (Figure 3). Further analysis showed that the family with the greatest relative abundance of sequences was the Gammaproteobacteria (Phylum: Proteobacteria), and different taxa showed distinct associations with the same soil physiochemical property (Table A1 in Appendix A). Forest type significantly affected the Chao1 estimator of richness, observed ASVs, and faith-PD index (Figure 4). Specifically, the Chao1 estimator of richness, observed ASVs, and faith-PD index of the soil bacterial community in Sec Forest were higher than those in the Mix, Bro, and Con forests. Further analysis showed they were positively related to TN, AN, and AP contents, whereas they were negatively correlated with soil pH (Figure 4; Table A2 in Appendix A).
The NMDS ordination (Figure 5) and similarity analysis (Table 1) showed that bacterial community structures in bulk soils at the topsoil layer differed significantly across forest types.
RDA indicated that the soil bacterial community structure was significantly correlated with soil physicochemical properties (Figure 6, Table A3 in Appendix A). At the phylum level, the first two principal components of RDA explained 41.33% of the variations in the soil bacterial community (Figure 6). Envfit analysis showed that SWC, pH, SOC, TN, C/N, N/P, and C/P were significantly associated with RDA1 and RDA2 (Table A2 in Appendix A).

3.3. Co-Occurrence Network Topological Index of the Soil Bacterial Community

In general, all co-occurrence networks of soil bacterial communities were dominated by positive edges. Additionally, the network topological features of the four examined forest types varied greatly (Table 2; Figure 7; Figure A4; and Table A3 in Appendix A). Soil bacterial community co-occurrence networks in Mix and Sec forests are more complex and stable than those in Bro and Con forests. As observed, the number of nodes and edges are higher in Mix forests compared to pure forests, regardless of vegetation function type. Additionally, positive links in soil bacterial networks were greater than negative links (Table 2). In the networks, a total of 2, 4, 17, and 18 ASVs were keystones in Bro, Con, Sec, and Mix forests, respectively (Figure A4; Table A3 in Appendix A).

3.4. Assembly of the Soil Bacterial Community

As presented, the niche width (2.92–8.85) of the soil bacterial community varied greatly across forest types. Specifically, the niche width of the soil bacterial community in Con was significantly discrete than that in the Bro, Mix, and Sec forests (Figure 8a). The βNTI of soil bacterial communities were between 2 and −2 (Figure 8b), and dispersal-limited ecological process predominated in community assembly (Figure 8c). The neutral model fitted well for all soil bacterial communities in Bro, Con, Sec, and Mix forests. In addition, it successfully allowed the estimation of the relationship between the frequency of ASVs and their relative abundance, yielding R2 values ranging from 0.53 to 0.698 (Figure 8d).

4. Discussion

4.1. Dominant Taxa

In the present study, Proteobacteria, Acidobacteriota, Actinobacteria, and Chloroflexi were the dominant phyla across forest types (Figure 3). The dominance of these phyla has also been observed in forest soils in different regions. In Changbai Nature Reserve of China, Acidobacteria was found to be dominant in mixed conifer-broadleaf forest, while Proteobacteria showed a higher proportion in the coniferous forest ecosystem [2]. In the Barren Hills of North China, Proteobacteria, Acidobacteria, and Actinobacteria were the dominant phyla [48]. In the Qinling Mountains, the dominant bacterial phyla in secondary forest soils were Proteobacteria, Acidobacteria, Firmicutes, and Verrucomicrobia [8]. The prevailing dominance of these phyla implies their strong adaptability and critical ecological functioning in the forest ecosystem. For instance, Proteobacteria can thrive in various soil conditions, exhibit morphological, physiological, and metabolic diversity, and be involved in the degradation of soil organic matter and the biogeochemical cycling of carbon, nitrogen, and sulfur [49]. Acidobacteriota play a crucial role in maintaining soil health and nutrient cycling and are responsible for the degradation of plant- and microorganism-based polysaccharides. Their abundance is strongly correlated to soil pH [50,51,52,53] and nitrogen content [51,54]. Actinobacteria play an important role as symbionts and as pathogens in plant-associated microbial communities [55], and their abundance is positively correlated with soil organic matter and pH [56,57,58]. Their presence helps maintain soil fertility and health. The phylum Chloroflexi participates in hydrolyzing polysaccharides such as cellulose, xylan, and chitin [58]. The phylum Chloroflexi participates in hydrolyzing polysaccharides such as cellulose, xylan, and chitin [58]. Their presence in forest soils plays a critical role in nutrient availability and ecosystem productivity.

4.2. Soil Bacterial Community Diversity and Similarity

Our results, combined with previous studies, demonstrate that forest vegetation types select specific bacterial communities [2,59,60,61]. We observed a higher Chao1 index, observed ASVs, and faith-PD index of the soil bacterial community in Sec Forest compared to other forest vegetation types (Figure 2). Moreover, they were positively correlated with TN, AN, and AP contents, whereas they were negatively correlated with soil pH (Figure 2; Table A1 in Appendix A). The difference in soil bacterial community diversity and similarity between the Con and Bro forests is not surprising and can be attributed to several factors, including differences in carbon and nitrogen availability and root exudates. Con forest soils are also low in nutrient availability, which favors slow-growing bacterial species that are adapted to nutrient-poor environments. Previous studies have shown divergent correlations between soil bacterial richness and soil physicochemical properties. For instance, soil bacterial richness is correlated with soil C and N contents as well as the C/N ratio [62]. Additionally, bacterial alpha diversity is largely influenced by SOC and TN contents [63,64]. Our findings, combined with those of other studies, suggest that SOC and TN contents were not significantly correlated with bacterial diversity indices [13,65]. Contrary to the reported negative correlations between soil bacterial diversity and soil C/N [46,66], no evident correlation was found between diversity index and soil C/N in the present study.
In agreement with previous studies [60,61], our study identified significant differences in the bacterial community composition across Bro forests, Mix forests, and Con forests. The most likely explanation would be the differences in amount and quality of litter and root exudates. However, it is also worth noting that the observed dissimilarity in soil bacterial community across forest types can be explained by differences in soil moisture, pH, SOC, TN, C/N, C/P, and N/P (Figure 6, Table A3 in Appendix A), implying intricate associations between the soil microbial community and soil physicochemical properties. Similar findings have been reported elsewhere. Soil pH is a critical factor that affects the structure of the soil bacterial community [67,68,69,70,71]. Soil pH, TC, TN, AP, and AK were found to be the main abiotic factors structuring the bacterial communities in the forest soils of the Baishilazi Nature Reserve, China [65]. Similarly, soil pH, SOC, and TN were found to be the most important factors affecting the bacterial community in Purple Mountain Park forest soils in China [60]. Soil pH, SOC, AP, and ammonia nitrogen have been shown to be the best indicators for explaining the changes in the soil bacterial community in the Barren Hills of North China. Furthermore, soil pH and AP were found to be the main determinants of the soil bacterial community [72]. Additionally, aluminum, nitrogen, calcium, nutrient availability, and pH also contribute to the dissimilarity in the soil bacterial community in beech forests [73]. Finally, soil resource elemental stoichiometry was also found to play an essential role in bacterial diversity and community composition [9,18,74,75].

4.3. Co-Occurrence Network

The most noticeable finding of the present study was that the co-occurrence networks of soil bacterial communities in Mix and Sec forests are more complex and stable than those in Bro and Con forests (Figure 7). Firstly, the number of nodes and edges is higher in Mix forests compared to pure forests, regardless of vegetation function type. As is known to us, in a microbial co-occurrence network, nodes and edges indicate microbes and statistically significant associations between nodes, respectively; modules may represent different niches and ecological processes driving community structure [76,77]. Nodes with high betweenness centrality scores are crucial for maintaining network stability and structure [77]. A high betweenness centrality indicates that more microbes bridge connections between modules in the sub-networks; edge overlap among sub-networks indicates a similarity in the microbial co-occurrence patterns among these environments [77]. Negative edges suggest varying intensities of competition or niche differentiation in different environments. Low proportions of negative edges in sub-networks suggest a prevalence of collaboration or niche sharing [77]. Higher clustering coefficients and modularity values indicate that ecosystems have fragmented niches for the environmental adaptation of soil microorganisms [66,78]. Additionally, positive interactions reflect cooperation and niche differentiation among species in the networks, whereas negative interactions indicate competition and niche overlap among species [77,79]. An increase in positive associations implies that the co-occurrence network between synergistic groups was enhanced [66,80,81]. We observed that there were more positive links than negative links in the soil bacterial networks. This finding is consistent with previous studies [82,83]. The likely reasons why positive linkages might outweigh negative linkages among soil bacterial taxa are as follows: (1) the presence of certain bacteria promotes the growth or activity of others, leading to the formation of cooperative relationships; (2) forest soils are rich and diverse in organic matter, providing a wide range of substrates for bacterial growth, and different bacterial taxa specialize in utilizing specific organic compounds, leading to a complementary utilization of resources and enhanced overall ecosystem functioning; (3) some bacterial taxa may create favorable microenvironments for other species to thrive. For instance, one species might produce compounds that enhance the growth of neighboring species, leading to a positive feedback loop; (4) groups of bacteria work together to perform complex metabolic processes; and (5) the sampling bias or experimental design may result in certain taxa being more easily detected or studied, leading to an apparent dominance of positive interactions. Further studies are needed to elucidate the underlying mechanism of the prevalence of positive linkages among soil bacterial taxa in the co-occurrence network of the soil bacterial community in forest soils.

4.4. Bacterial Community Assembly Process

Identifying the relative contribution of deterministic and stochastic processes in soil microbial community assembly is a focal point in soil ecology. Deterministic processes, including environmental selection, species interactions, and niche differentiation, underpin the roles of biotic and abiotic filtering, leading to significant variation in community composition under different environmental conditions. Stochastic processes, such as dispersal limitation, undominated dispersal, and homogenizing dispersal, highlight the contribution of probabilistic dispersal and ecological drift to community composition patterns [46,71]. Recent evidence also supports that life-history traits are important determinants governing soil microbial community assembly in terrestrial ecosystems [53,84]. The intricate network of soil bacterial communities is a fundamental component of forest ecosystems, playing a crucial role in nutrient cycling, organic matter decomposition, and overall ecosystem functioning. Understanding the assembly processes that govern the structure and composition of soil bacterial communities in forest ecosystems is crucial for devising effective strategies for sustainable forest management, conservation, and the preservation of essential ecosystem services.
Previous studies suggest that soil bacterial community assembly is driven by a combination of deterministic and stochastic processes [8,46], but the relative importance of these two processes varies greatly depending on forest stand age [8], environmental factors [53,85,86,87], plant diversity [8], SOC [85,88], soil pH [86,89], soil moisture [87,89], and soil fungi [90] as critical mediators. Numerous pieces of evidence support that the soil bacterial community is primarily dominated by stochastic processes [53,71,84,91]. In contrast, the findings of the present study highlight the dominance of deterministic processes. The likely reasons are as follows: Firstly, we did not account for the potential effects of vegetation composition and stand age on the soil bacterial community. In previous studies, vegetation, associated soil physicochemical properties, and stand age have been shown to influence bacterial community composition [70,72,81]. Additionally, we did not consider the effects of altitude on the soil bacterial community, despite the pronounced altitudinal patterns observed in soil bacterial communities [70,89,92,93]. Finally, the observed patterns are based on a single sampling event. However, the effects of forest vegetation type on the soil bacterial community were seasonal-dependent [94], and the assembly processes of soil bacterial communities are dynamic and context-dependent [8,31,87]. Despite these limitations, the present study demonstrates that forest type is an important factor shaping the soil bacterial community, and its community assembly processes were shaped by a combination of deterministic and stochastic ecological processes.

5. Conclusions and Implications

In conclusion, the findings provide solid empirical evidence of close interlinkages among forest vegetation type, soil physicochemical properties, and the soil bacterial community, supporting our first and second hypotheses. Furthermore, Mix and Sec forests foster complex and robust co-occurrence networks compared with Bro and Con forests, implying increased biodiversity can enhance ecosystem resilience and functioning. However, the soil bacterial community of the four examined forests is consistently primarily controlled by deterministic processes, leading to the rejection of our third hypothesis. These findings enhance the holistic understanding of the interlinkages among forest vegetation type, soil physicochemical properties, and the soil bacterial community, which is essential for effective forest management, biodiversity conservation, ecological restoration, and sustainable ecosystem functioning. Studying the relationships between forest vegetation type, soil properties, and soil bacteria allows forest managers to assess the overall health and functioning of the ecosystem. This information can guide conservation efforts and sustainable management practices. Managing these interlinkages effectively enables forest managers to promote biodiversity conservation, protect endangered species, and implement strategies to enhance soil health and productivity. Knowledge of the interlinkages can also guide ecological restoration efforts in degraded forest ecosystems. By restoring the balance between vegetation, soil properties, and soil bacteria; restoration projects can be more effective in rebuilding healthy and resilient ecosystems. Finally, managing these interlinkages can help mitigate the effects of climate change by promoting carbon storage in soils and reducing emissions, as soil bacteria contribute to carbon sequestration and greenhouse gas emissions.

Author Contributions

B.Y. designed the research and acquired funding grants. B.Y. and Z.Y. supervised the project. B.Y., W.F. and Z.Y. refined the experimental method. B.Y., W.F., W.Z., K.H. and Z.Y. performed the research and conducted the data analyses. B.Y. drafted the manuscript. W.F., W.Z., K.H. and Z.Y. reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was combined supported by Sichuan Science and Technology Program (2021YFH0158), the Innovation Team Program of Sichuan Forestry and Grassland Bureau (CXTD202300601) and National Natural Science Foundation of China (U21A20193).

Data Availability Statement

The datasets presented in this study can be found in online repositories.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Summary of the Spearman’s rank correlation between soil bacterial taxa’s relative abundance and soil physiochemical properties.
Table A1. Summary of the Spearman’s rank correlation between soil bacterial taxa’s relative abundance and soil physiochemical properties.
TaxonSWCpHAPSOCTPTNANC/NN/PC/P
Proteobacteria−0.110.52 ***0.17−0.01−0.14−0.20−0.190.14−0.100.02
Acidobacteria0.25−0.69 ***0.100.130.150.290.26−0.110.220.12
Actinobacteria−0.48 **0.41 **0.03−0.31−0.23−0.40 *−0.40 **−0.03−0.24−0.23
Chloroflexi−0.070.11−0.32−0.030.17−0.13−0.090.27−0.29−0.19
Bacteroidetes−0.150.250.17−0.050.02−0.13−0.170.02−0.17−0.04
Planctomycetes0.000.100.06−0.10−0.51 **−0.06−0.05−0.160.340.17
Verrucomicrobia0.27−0.06−0.110.240.160.36 *0.40 *−0.090.180.10
Thaumarchaeota0.30−0.57 ***−0.130.150.130.320.31 **−0.090.210.10
Gemmatimonadetes−0.46 **0.33 *0.00−0.39 *0.18−0.46 **−0.46 **0.05−0.69 ***−0.49 ***
Rokubacteria−0.34 *0.61 ***−0.20−0.19−0.02−0.38 *−0.35 *0.28−0.38 **−0.15
Note: *, p ≤ 0.05; **, p < 0.01; and ***, p < 0.001. AN is available nitrogen, AP is available phosphorus, TN is soil total nitrogen, TP is total phosphorus, SWC is soil water content, SOC is soil total organic carbon, C/N is the ratio of soil total organic carbon to total nitrogen, C/P is the ratio of soil total organic carbon to total phosphorus, and N/P is the ratio of total nitrogen to total phosphorus.
Table A2. Summary of the Spearman’s rank correlation between soil bacterial diversity indices and soil physiochemical properties.
Table A2. Summary of the Spearman’s rank correlation between soil bacterial diversity indices and soil physiochemical properties.
VariableSWCpHAPSOCTPTNANC/NN/PC/P
Chao1−0.33 *0.53 ***−0.05−0.25−0.41 *−0.37 *−0.40 *0.08−0.170.00
Faith_pd−0.240.69 ***−0.18−0.18−0.43 **−0.38 *−0.41 *0.29−0.110.10
Goods_coverage−0.01−0.150.02−0.11−0.12−0.04−0.02−0.080.07−0.11
Observed_species−0.33 *0.53 ***−0.05−0.25−0.41 *−0.37 *−0.40 *0.08−0.170.00
Pielou_e−0.04−0.260.05−0.28−0.21−0.07−0.10−0.26−0.01−0.16
Shannon−0.14−0.040.05−0.26−0.23−0.17−0.21−0.15−0.11−0.10
Simpson−0.06−0.310.05−0.28−0.13−0.05−0.08−0.29−0.04−0.19
Note: *, p ≤ 0.05; **, p < 0.01; and ***, p < 0.001. AN is available nitrogen, AP is available phosphorus, TN is soil total nitrogen, TP is total phosphorus, SWC is soil water content, SOC is soil total organic carbon, C/N is the ratio of soil total organic carbon to total nitrogen, C/P is the ratio of soil total organic carbon to total phosphorus, and N/P is the ratio of total nitrogen to total phosphorus.
Table A3. Summary of the envfit analysis (envfit function used with 999 permutations) shows soil properties variables that correlate to bacterial community structure.
Table A3. Summary of the envfit analysis (envfit function used with 999 permutations) shows soil properties variables that correlate to bacterial community structure.
VariableRDA1RDA2r2Pr (>r)
SWC0.37300.92780.530.001
pH−0.9950−0.10020.810.001
SOC0.12310.99240.220.013
TN0.50180.86500.470.001
C/N−0.9625−0.27120.160.048
N/P0.40590.91390.720.001
C/P0.06240.99800.360.001
Note: TN is soil total nitrogen, SWC is soil water content, SOC is soil total organic carbon, C/N is the ratio of soil total organic carbon to total nitrogen, C/P is the ratio of soil total organic carbon to total phosphorus, and N/P is the ratio of total nitrogen to total phosphorus.
Microorganisms 12 00728 g0a1
Figure A1. The circular packing represents the hierarchical taxonomic tree to organize discovered microbial taxa (red circles) at five different taxonomic levels: phylum, class, order, family, and genus.
Figure A1. The circular packing represents the hierarchical taxonomic tree to organize discovered microbial taxa (red circles) at five different taxonomic levels: phylum, class, order, family, and genus.
Microorganisms 12 00728 g0a1
Microorganisms 12 00728 g0a2
Figure A2. Cladogram of the soil bacterial communities in the five vegetation types of forest soils with liner discriminant analysis (LDA) > 2.65. The circles represent bacterial taxa from phyla to genera, starting from the center. Bro is evergreen broadleaved forest, Con is coniferous forest, Mix is deciduous broadleaved and evergreen broadleaved mixed forest, and Sec is secondary deciduous broadleaved and coniferous mixed forest.
Figure A2. Cladogram of the soil bacterial communities in the five vegetation types of forest soils with liner discriminant analysis (LDA) > 2.65. The circles represent bacterial taxa from phyla to genera, starting from the center. Bro is evergreen broadleaved forest, Con is coniferous forest, Mix is deciduous broadleaved and evergreen broadleaved mixed forest, and Sec is secondary deciduous broadleaved and coniferous mixed forest.
Microorganisms 12 00728 g0a2
Microorganisms 12 00728 g0a3
Figure A3. Phylogenetic tree showing the top 50 amplicon sequence variants (ASVs) with the highest cumulative relative abundance. The taxonomic affiliation at the family level of different Phyla of each ASV is identified by the color range in the below panel and within the tree.
Figure A3. Phylogenetic tree showing the top 50 amplicon sequence variants (ASVs) with the highest cumulative relative abundance. The taxonomic affiliation at the family level of different Phyla of each ASV is identified by the color range in the below panel and within the tree.
Microorganisms 12 00728 g0a3
Microorganisms 12 00728 g0a4
Figure A4. Scatter plot of Zi (within-module connectivity) and Pi (among-module connectivity).
Figure A4. Scatter plot of Zi (within-module connectivity) and Pi (among-module connectivity).
Microorganisms 12 00728 g0a4

References

  1. Bardgett, R.D.; van der Putten, W.H. Belowground biodiversity and ecosystem functioning. Nature 2014, 515, 505–511. [Google Scholar] [CrossRef] [PubMed]
  2. Li, H.; Ye, D.; Wang, X.; Settles, M.L.; Wang, J.; Hao, Z.; Zhou, L.; Dong, P.; Jiang, Y.; Ma, Z. Soil bacterial communities of different natural forest types in Northeast China. Plant Soil 2014, 383, 203–216. [Google Scholar] [CrossRef]
  3. Delgado-Baquerizo, M.; Oliverio, A.M.; Brewer, T.E.; Benavent-González, A.; Eldridge, D.J.; Bardgett, R.D.; Maestre, F.T.; Singh, B.K.; Fierer, N. A global atlas of the dominant bacteria found in soil. Science 2018, 359, 320–325. [Google Scholar] [CrossRef] [PubMed]
  4. Wei, H.; Peng, C.; Yang, B.; Song, H.; Li, Q.; Jiang, L.; Wei, G.; Wang, K.; Wang, H.; Liu, S.; et al. Contrasting Soil Bacterial Community, Diversity, and Function in Two Forests in China. Front. Microbiol. 2018, 9, 1693. [Google Scholar] [CrossRef]
  5. Fierer, N.; Jackson, R.B. The diversity and biogeography of soil bacterial communities. Proc. Natl. Acad. Sci. USA 2006, 10, 626–631. [Google Scholar] [CrossRef] [PubMed]
  6. Hermans, S.M.; Buckley, H.L.; Case, B.S.; Curran-Cournane, F.; Taylor, M.; Lear, G.; Loeffler, F.E. Bacteria as Emerging Indicators of Soil Condition. Appl. Environ. Microbiol. 2017, 83, e02826-16. [Google Scholar] [CrossRef] [PubMed]
  7. Sawada, K.; Inagaki, Y.; Sugihara, S.; Funakawa, S.; Ritz, K.; Toyota, K. Impacts of conversion from natural forest to cedar plantation on the structure and diversity of root-associated and soil microbial communities. Appl. Soil Ecol. 2021, 167, 104027. [Google Scholar] [CrossRef]
  8. Huo, X.; Ren, C.; Wang, D.; Wu, R.; Wang, Y.; Li, Z.; Huang, D.; Qi, H. Microbial community assembly and its influencing factors of secondary forests in Qinling Mountains. Soil Biol. Biochem. 2023, 184, 109075. [Google Scholar] [CrossRef]
  9. Xue, L.; Ren, H.; Li, S.; Leng, X.; Yao, X. Soil Bacterial Community Structure and Co-occurrence Pattern during Vegetation Restoration in Karst Rocky Desertification Area. Front. Microbiol. 2017, 8, 2377. [Google Scholar] [CrossRef]
  10. Banerjee, S.; Helgason, B.; Wang, L.; Winsley, T.; Ferrari, B.C.; Siciliano, S.D. Legacy effects of soil moisture on microbial community structure and N2O emissions. Soil Biol. Biochem. 2016, 95, 40–50. [Google Scholar] [CrossRef]
  11. Rousk, J.; Bååth, E.; Brookes, P.C.; Lauber, C.L.; Lozupone, C.; Caporaso, J.G.; Knight, R.; Fierer, N. Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 2010, 4, 1340–1351. [Google Scholar] [CrossRef]
  12. Landesman, W.J.; Nelson, D.M.; Fitzpatrick, M.C. Soil properties and tree species drive ß-diversity of soil bacterial communities. Soil Biol. Biochem. 2014, 76, 201–209. [Google Scholar] [CrossRef]
  13. Lin, Y.T.; Hu, H.W.; Whitman, W.B.; Coleman, D.C.; Chiu, C.Y. Comparison of soil bacterial communities in a natural hardwood forest and coniferous plantations in perhumid subtropical low mountains. Bot. Stud. 2014, 55, 50. [Google Scholar] [CrossRef] [PubMed]
  14. Jiang, Y.; Lei, Y.; Yang, Y.; Korpelainen, H.; Niinemets, Ü.; Li, C. Divergent assemblage patterns and driving forces for bacterial and fungal communities along a glacier forefield chronosequence. Soil Biol. Biochem. 2018, 118, 207–216. [Google Scholar] [CrossRef]
  15. Liu, G.y.; Chen, L.l.; Shi, X.r.; Yuan, Z.y.; Yuan, L.Y.; Lock, T.R.; Kallenbach, R.L. Changes in rhizosphere bacterial and fungal community composition with vegetation restoration in planted forests. Land Degrad. Dev. 2019, 30, 1147–1157. [Google Scholar] [CrossRef]
  16. Kuramae, E.E.; Yergeau, E.; Wong, L.C.; Pijl, A.S.; Veen, J.A.; Kowalchuk, G.A. Soil characteristics more strongly influence soil bacterial communities than land-use type. FEMS Microbiol. Ecol. 2012, 79, 12–24. [Google Scholar] [CrossRef]
  17. Razanamalala, K.; Razafimbelo, T.; Maron, P.-A.; Ranjard, L.; Chemidlin, N.; Lelièvre, M.; Dequiedt, S.; Ramaroson, V.H.; Marsden, C.; Becquer, T.; et al. Soil microbial diversity drives the priming effect along climate gradients: A case study in Madagascar. ISME J. 2018, 12, 451–462. [Google Scholar] [CrossRef] [PubMed]
  18. Delgado-Baquerizo, M.; Reich, P.B.; Khachane, A.N.; Campbell, C.D.; Thomas, N.; Freitag, T.E.; Abu Al-Soud, W.; Sørensen, S.; Bardgett, R.D.; Singh, B.K. It is elemental: Soil nutrient stoichiometry drives bacterial diversity. Environ. Microbiol. 2017, 19, 1176–1188. [Google Scholar] [CrossRef]
  19. Sun, H.; Terhonen, E.; Koskinen, K.; Paulin, L.; Kasanen, R.; Asiegbu, F.O. Bacterial diversity and community structure along different peat soils in boreal forest. Appl. Soil Ecol. 2014, 74, 37–45. [Google Scholar] [CrossRef]
  20. Li, X.; He, Z.; Yao, M.; Kou, Y.; Li, C.; Li, J. Distinct mechanisms shape soil bacterial and fungal co-occurrence networks in a mountain ecosystem. FEMS Microbiol. Ecol. 2020, 96, fiaa030. [Google Scholar] [CrossRef]
  21. Tu, Q.; Yan, Q.; Deng, Y.; Michaletz, S.T.; Buzzard, V.; Weiser, M.D.; Waide, R.; Ning, D.; Wu, L.; He, Z.; et al. Biogeographic patterns of microbial co-occurrence ecological networks in six American forests. Soil Biol. Biochem. 2020, 148, 107897. [Google Scholar] [CrossRef]
  22. Xu, M.; Lu, X.; Xu, Y.; Zhong, Z.; Zhang, W.; Ren, C.; Han, X.; Yang, G.; Feng, Y. Dynamics of bacterial community in litter and soil along a chronosequence of Robinia pseudoacacia plantations. Sci. Total Environ. 2020, 703, 135613. [Google Scholar] [CrossRef] [PubMed]
  23. Qiang, W.; He, L.; Zhang, Y.; Liu, B.; Liu, Y.; Liu, Q.; Pang, X. Aboveground vegetation and soil physicochemical properties jointly drive the shift of soil microbial community during subalpine secondary succession in southwest China. Catena 2021, 202, 105251. [Google Scholar] [CrossRef]
  24. Barberán, A.; Bates, S.T.; Casamayor, E.O.; Fierer, N. Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J. 2012, 6, 343–351. [Google Scholar] [CrossRef] [PubMed]
  25. Ma, B.; Wang, H.; Dsouza, M.; Lou, J.; He, Y.; Dai, Z.; Brookes, P.C.; Xu, J.; Gilbert, J.A. Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China. ISME J. 2016, 10, 1891–1901. [Google Scholar] [CrossRef] [PubMed]
  26. Banerjee, S.; Kirkby, C.A.; Schmutter, D.; Bissett, A.; Kirkegaard, J.A.; Richardson, A.E. Network analysis reveals functional redundancy and keystone taxa amongst bacterial and fungal communities during organic matter decomposition in an arable soil. Soil Biol. Biochem. 2016, 97, 188–198. [Google Scholar] [CrossRef]
  27. Banerjee, S.; Schlaeppi, K.; van der Heijden, M.G. Keystone taxa as drivers of microbiome structure and functioning. Nat. Rev. Microbiol. 2018, 16, 567–576. [Google Scholar] [CrossRef] [PubMed]
  28. He, D.; Shen, W.; Eberwein, J.; Zhao, Q.; Ren, L.; Wu, Q.L. Diversity and co-occurrence network of soil fungi are more responsive than those of bacteria to shifts in precipitation seasonality in a subtropical forest. Soil Biol. Biochem. 2017, 115, 499–510. [Google Scholar] [CrossRef]
  29. Vandermeer, J.H. Niche Theory. Annu. Rev. Ecol. Syst. 1972, 3, 107–132. [Google Scholar] [CrossRef]
  30. Hubbell, S.P. Neutral Theory in Community Ecology and the Hypothesis of Functional Equivalence. Funct. Ecol. 2005, 19, 166–172. [Google Scholar] [CrossRef]
  31. Langenheder, S.; Lindström, E.S. Factors influencing aquatic and terrestrial bacterial community assembly. Environ. Microbiol. Rep. 2019, 11, 306–315. [Google Scholar] [CrossRef] [PubMed]
  32. Zhao, P.; Liu, J.; Jia, T.; Wang, Y.; Chai, B. Environmental filtering drives bacterial community structure and function in a subalpine area of northern China. J. Basic Microbiol. 2018, 59, 337–347. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, Z.; Swaisgood, R.R.; Zhang, S.; Nordstrom, L.A.; Wang, H.; Gu, X.; Hu, J.; Wei, F. Old-growth forest is what giant pandas really need. Biol. Lett. 2011, 7, 403–406. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, Y.; Wang, X.; Zheng, X.; Gong, Y.; Chen, M.; Qiu, L.; Zhou, H.; Wei, W.; Han, H. Space use and microhabitat selection of wild giant pandas in Meigu Dafengding National Nature Reserve, China. Front. Ecol. Evol. 2022, 10, 1000841. [Google Scholar] [CrossRef]
  35. Hong, M.; Yuan, S.; Yang, Z.; Yang, X.; Gu, X.; Huang, F.; Zhang, Z. Comparison of microhabitat selection and trace abundance of giant pandas between primary and secondary forests in Liziping Nature Reserve, China: Effects of selective logging. Mamm. Biol. 2015, 80, 373–379. [Google Scholar] [CrossRef]
  36. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. QIIME 2: Reproducible, interactive, scalable, and extensible microbiome data science. PeerJ 2018, preprints. [Google Scholar] [CrossRef]
  37. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, Y.-F.; Chen, P.; Wang, F.-H.; Han, W.-X.; Qiao, M.; Dong, W.-X.; Hu, C.-S.; Zhu, D.; Chu, H.-Y.; Zhu, Y.-G. The ecological clusters of soil organisms drive the ecosystem multifunctionality under long-term fertilization. Environ. Int. 2022, 161, 107133. [Google Scholar] [CrossRef] [PubMed]
  39. Katoh, K.; Misawa, K.; Kuma, K.I.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef]
  40. Poon, A.F.Y.; Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree 2—Approximately Maximum-Likelihood Trees for Large Alignments. PLoS ONE 2010, 5, e9490. [Google Scholar] [CrossRef]
  41. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  42. Csardi, M.G. igraph_ Network Analysis and Visualization. Available online: https://igraph.org/ (accessed on 26 January 2024).
  43. Sloan, W.T.; Lunn, M.; Woodcock, S.; Head, I.M.; Nee, S.; Curtis, T.P. Quantifying the roles of immigration and chance in shaping prokaryote community structure. Environ. Microbiol. 2006, 8, 732–740. [Google Scholar] [CrossRef] [PubMed]
  44. Stegen, J.C.; Lin, X.; Konopka, A.E.; Fredrickson, J.K. Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J. 2012, 6, 1653–1664. [Google Scholar] [CrossRef]
  45. Dong, Y.; Wu, S.; Deng, Y.; Wang, S.; Fan, H.; Li, X.; Bai, Z.; Zhuang, X. Distinct Functions and Assembly Mechanisms of Soil Abundant and Rare Bacterial Taxa Under Increasing Pyrene Stresses. Front. Microbiol. 2021, 12, 689762. [Google Scholar] [CrossRef] [PubMed]
  46. Stegen, J.C.; Lin, X.; Fredrickson, J.K.; Chen, X.; Kennedy, D.W.; Murray, C.J.; Rockhold, M.L.; Konopka, A. Quantifying community assembly processes and identifying features that impose them. ISME J. 2013, 7, 2069–2079. [Google Scholar] [CrossRef] [PubMed]
  47. Zhou, J.Z.; Ning, D.L. Stochastic Community Assembly: Does It Matter in Microbial Ecology? Microbiol. Mol. Biol. Rev. 2017, 81, e00002-17. [Google Scholar] [CrossRef]
  48. Qiu, Z.; Shi, C.; Zhao, M.; Wang, K.; Zhang, M.; Wang, T.; Shi, F. Improving Effects of Afforestation with Different Forest Types on Soil Nutrients and Bacterial Community in Barren Hills of North China. Sustainability 2022, 14, 1202. [Google Scholar] [CrossRef]
  49. Spain, A.M.; Krumholz, L.R.; Elshahed, M.S. Abundance, composition, diversity and novelty of soil Proteobacteria. ISME J. 2009, 3, 992–1000. [Google Scholar] [CrossRef]
  50. Eichorst, S.A.; Breznak, J.A.; Schmidt, T.M. Isolation and Characterization of Soil Bacteria That Define Terriglobus gen. nov., in the Phylum Acidobacteria. Appl. Environ. Microbiol. 2007, 73, 2708–2717. [Google Scholar] [CrossRef]
  51. Jones, R.T.; Robeson, M.S.; Lauber, C.L.; Hamady, M.; Knight, R.; Fierer, N. A comprehensive survey of soil acidobacterial diversity using pyrosequencing and clone library analyses. ISME J. 2009, 3, 442–453. [Google Scholar] [CrossRef]
  52. Rodrigues, J.L.M.; Pellizari, V.H.; Mueller, R.; Baek, K.; Jesus, E.d.C.; Paula, F.S.; Mirza, B.; Hamaoui, G.S.; Tsai, S.M.; Feigl, B.; et al. Conversion of the Amazon rainforest to agriculture results in biotic homogenization of soil bacterial communities. Proc. Natl. Acad. Sci. USA 2012, 110, 988–993. [Google Scholar] [CrossRef] [PubMed]
  53. Chen, Y.; Yang, X.; Fu, W.; Chen, B.; Hu, H.; Feng, K.; Geisen, S. Conversion of natural grassland to cropland alters microbial community assembly across northern China. Environ. Microbiol. 2022, 24, 5630–5642. [Google Scholar] [CrossRef] [PubMed]
  54. Fierer, N.; Bradford, M.A.; Jackson, R.B. Toward an Ecological Classification of Soil Bacteria. Ecology 2007, 88, 1354–1364. [Google Scholar] [CrossRef] [PubMed]
  55. Barka, E.A.; Vatsa, P.; Sanchez, L.; Gaveau-Vaillant, N.; Jacquard, C.; Klenk, H.-P.; Clément, C.; Ouhdouch, Y.; van Wezel, G.P. Taxonomy, Physiology, and Natural Products of Actinobacteria. Microbiol. Mol. Biol. Rev. 2016, 80, 1–43. [Google Scholar] [CrossRef] [PubMed]
  56. Lauber, C.L.; Hamady, M.; Knight, R.; Fierer, N. Pyrosequencing-Based Assessment of Soil pH as a Predictor of Soil Bacterial Community Structure at the Continental Scale. Appl. Environ. Microbiol. 2009, 75, 5111–5120. [Google Scholar] [CrossRef]
  57. Li, X.; Sun, M.; Zhang, H.; Xu, N.; Sun, G. Use of mulberry–soybean intercropping in salt–alkali soil impacts the diversity of the soil bacterial community. Microb. Biotechnol. 2016, 9, 293–304. [Google Scholar] [CrossRef] [PubMed]
  58. Lin, Y.-T.; Whitman, W.B.; Coleman, D.C.; Chiu, C.-Y. Effects of Reforestation on the Structure and Diversity of Bacterial Communities in Subtropical Low Mountain Forest Soils. Front. Microbiol. 2018, 9, 1968. [Google Scholar] [CrossRef]
  59. Hackl, E.; Pfeffer, M.; Donat, C.; Bachmann, G.; Zechmeisterboltenstern, S. Composition of the microbial communities in the mineral soil under different types of natural forest. Soil Biol. Biochem. 2005, 37, 661–671. [Google Scholar] [CrossRef]
  60. Qu, Z.L.; Liu, B.; Ma, Y.; Xu, J.; Sun, H.; Gallery, R. The response of the soil bacterial community and function to forest succession caused by forest disease. Funct. Ecol. 2020, 34, 2548–2559. [Google Scholar] [CrossRef]
  61. Sui, X.; Li, M.; Frey, B.; Wang, M.; Weng, X.; Wang, X.; Chen, F.; Li, X.; Du, Z.; Yang, L.; et al. Climax forest has a higher soil bacterial diversity but lower soil nutrient contents than degraded forests in temperate northern China. Ecol. Evol. 2022, 12, e9535. [Google Scholar] [CrossRef]
  62. Wang, Y.; Chen, L.; Xiang, W.; Ouyang, S.; Zhang, T.; Zhang, X.; Zeng, Y.; Hu, Y.; Luo, G.; Kuzyakov, Y. Forest conversion to plantations: A meta-analysis of consequences for soil and microbial properties and functions. Glob. Chang. Biol. 2021, 27, 5643–5656. [Google Scholar] [CrossRef] [PubMed]
  63. Siles, J.A.; Margesin, R. Abundance and Diversity of Bacterial, Archaeal, and Fungal Communities Along an Altitudinal Gradient in Alpine Forest Soils: What Are the Driving Factors? Microb. Ecol. 2016, 72, 207–220. [Google Scholar] [CrossRef] [PubMed]
  64. Ren, C.; Zhang, W.; Zhong, Z.; Han, X.; Yang, G.; Feng, Y.; Ren, G. Differential responses of soil microbial biomass, diversity, and compositions to altitudinal gradients depend on plant and soil characteristics. Sci. Total Environ. 2018, 610–611, 750–758. [Google Scholar] [CrossRef] [PubMed]
  65. Deng, J.; Yin, Y.; Zhu, W.; Zhou, Y. Variations in Soil Bacterial Community Diversity and Structures Among Different Revegetation Types in the Baishilazi Nature Reserve. Front. Microbiol. 2018, 9, 2874. [Google Scholar] [CrossRef] [PubMed]
  66. Wang, R.; Wang, Y.; Zheng, W.; Hou, F.; Hu, Y.; Guo, S. Converting croplands to orchards changes soil microbial community composition and co-occurrence patterns. Land Degrad. Dev. 2021, 32, 2509–2519. [Google Scholar] [CrossRef]
  67. Ferrenberg, S.; O’Neill, S.P.; Knelman, J.E.; Todd, B.; Duggan, S.; Bradley, D.; Robinson, T.; Schmidt, S.K.; Townsend, A.R.; Williams, M.W.; et al. Changes in assembly processes in soil bacterial communities following a wildfire disturbance. ISME J. 2013, 7, 1102–1111. [Google Scholar] [CrossRef]
  68. Shen, C.; Xiong, J.; Zhang, H.; Feng, Y.; Lin, X.; Li, X.; Liang, W.; Chu, H. Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biol. Biochem. 2013, 57, 204–211. [Google Scholar] [CrossRef]
  69. Guo, S.; Han, X.; Li, H.; Wang, T.; Tong, X.; Ren, G.; Feng, Y.; Yang, G. Evaluation of soil quality along two revegetation chronosequences on the Loess Hilly Region of China. Sci. Total Environ. 2018, 633, 808–815. [Google Scholar] [CrossRef]
  70. Li, J.; Shen, Z.; Li, C.; Kou, Y.; Wang, Y.; Tu, B.; Zhang, S.; Li, X. Stair-Step Pattern of Soil Bacterial Diversity Mainly Driven by pH and Vegetation Types Along the Elevational Gradients of Gongga Mountain, China. Front. Microbiol. 2018, 9, 569. [Google Scholar] [CrossRef]
  71. He, G.; Peng, T.; Guo, Y.; Wen, S.; Ji, L.; Luo, Z. Forest succession improves the complexity of soil microbial interaction and ecological stochasticity of community assembly: Evidence from Phoebe bournei-dominated forests in subtropical regions. Front. Microbiol. 2022, 13, 1021258. [Google Scholar] [CrossRef]
  72. Bi, B.; Zhang, H.; Yuan, Y.; Wu, Z.; Wang, Y.; Han, F. Dynamic changes of soil microbial community in Pinus sylvestris var. mongolica plantations in the Mu Us Sandy Land. J. Environ. Manag. 2021, 287, 112306. [Google Scholar] [CrossRef] [PubMed]
  73. Jeanbille, M.; Buée, M.; Bach, C.; Cébron, A.; Frey-Klett, P.; Turpault, M.P.; Uroz, S. Soil Parameters Drive the Structure, Diversity and Metabolic Potentials of the Bacterial Communities Across Temperate Beech Forest Soil Sequences. Microb. Ecol. 2015, 71, 482–493. [Google Scholar] [CrossRef] [PubMed]
  74. Högberg, M.N.; Högberg, P.; Myrold, D.D. Is microbial community composition in boreal forest soils determined by pH, C-to-N ratio, the trees, or all three? Oecologia 2006, 150, 590–601. [Google Scholar] [CrossRef] [PubMed]
  75. Wan, X.; Huang, Z.; He, Z.; Yu, Z.; Wang, M.; Davis, M.R.; Yang, Y. Soil C:N ratio is the major determinant of soil microbial community structure in subtropical coniferous and broadleaf forest plantations. Plant Soil 2014, 387, 103–116. [Google Scholar] [CrossRef]
  76. Röttjers, L.; Faust, K. From hairballs to hypotheses–biological insights from microbial networks. FEMS Microbiol. Rev. 2018, 42, 761–780. [Google Scholar] [CrossRef] [PubMed]
  77. Ma, B.; Wang, Y.; Ye, S.; Liu, S.; Stirling, E.; Gilbert, J.A.; Faust, K.; Knight, R.; Jansson, J.K.; Cardona, C.; et al. Earth microbial co-occurrence network reveals interconnection pattern across microbiomes. Microbiome 2020, 8, 82. [Google Scholar] [CrossRef] [PubMed]
  78. Faust, K.; Raes, J. Microbial interactions: From networks to models. Nat. Rev. Microbiol. 2012, 10, 538–550. [Google Scholar] [CrossRef] [PubMed]
  79. Ghoul, M.; Mitri, S. The Ecology and Evolution of Microbial Competition. Trends Microbiol. 2016, 24, 833–845. [Google Scholar] [CrossRef] [PubMed]
  80. Zhang, B.; Zhang, J.; Liu, Y.; Shi, P.; Wei, G. Co-occurrence patterns of soybean rhizosphere microbiome at a continental scale. Soil Biol. Biochem. 2018, 118, 178–186. [Google Scholar] [CrossRef]
  81. Yu, Z.; Liang, K.; Huang, G.; Wang, X.; Lin, M.; Chen, Y.; Zhou, Z. Soil Bacterial Community Shifts Are Driven by Soil Nutrient Availability along a Teak Plantation Chronosequence in Tropical Forests in China. Biology 2021, 10, 1329. [Google Scholar] [CrossRef]
  82. Wei, Y.; Chen, L.; Feng, Q.; Xi, H.; Zhang, C.; Gan, K.; Yong, T. Structure and assembly mechanism of soil bacterial community under different soil salt intensities in arid and semiarid regions. Ecol. Indic. 2024, 158, 111631. [Google Scholar] [CrossRef]
  83. Kang, H.; Chai, A.; Lin, Z.; Shi, Y.; Xie, X.; Li, L.; Fan, T.; Xiang, S.; Xie, J.; Li, B. Deciphering Differences in Microbial Community Diversity between Clubroot-Diseased and Healthy Soils. Microorganisms 2024, 12, 251. [Google Scholar] [CrossRef]
  84. Luan, L.; Jiang, Y.; Cheng, M.; Dini-Andreote, F.; Sui, Y.; Xu, Q.; Geisen, S.; Sun, B. Organism body size structures the soil microbial and nematode community assembly at a continental and global scale. Nat. Commun. 2020, 11, 6406. [Google Scholar] [CrossRef] [PubMed]
  85. Dini-Andreote, F.; Stegen, J.C.; van Elsas, J.D.; Salles, J.F. Disentangling mechanisms that mediate the balance between stochastic and deterministic processes in microbial succession. Proc. Natl. Acad. Sci. USA 2015, 112, E1326–E1332. [Google Scholar] [CrossRef]
  86. Tripathi, B.M.; Stegen, J.C.; Kim, M.; Dong, K.; Adams, J.M.; Lee, Y.K. Soil pH mediates the balance between stochastic and deterministic assembly of bacteria. ISME J. 2018, 12, 1072–1083. [Google Scholar] [CrossRef] [PubMed]
  87. Gao, G.-F.; Peng, D.; Zhang, Y.; Li, Y.; Fan, K.; Tripathi, B.M.; Adams, J.M.; Chu, H. Dramatic change of bacterial assembly process and co-occurrence pattern in Spartina alterniflora salt marsh along an inundation frequency gradient. Sci. Total Environ. 2021, 755, 142546. [Google Scholar] [CrossRef]
  88. Zhang, Q.; Wang, X.; Zhang, Z.; Liu, H.; Liu, Y.; Feng, Y.; Yang, G.; Ren, C.; Han, X. Linking soil bacterial community assembly with the composition of organic carbon during forest succession. Soil Biol. Biochem. 2022, 173, 108790. [Google Scholar] [CrossRef]
  89. Zou, Z.; Yuan, K.; Ming, L.; Li, Z.; Yang, Y.; Yang, R.; Cheng, W.; Liu, H.; Jiang, J.; Luan, T.; et al. Changes in Alpine Soil Bacterial Communities With Altitude and Slopes at Mount Shergyla, Tibetan Plateau: Diversity, Structure, and Influencing Factors. Front. Microbiol. 2022, 13, 839499. [Google Scholar] [CrossRef]
  90. Jiao, S.; Chu, H.; Zhang, B.; Wei, X.; Chen, W.; Wei, G. Linking soil fungi to bacterial community assembly in arid ecosystems. iMeta 2022, 1, e2. [Google Scholar] [CrossRef]
  91. Kang, L.; Chen, L.; Zhang, D.; Peng, Y.; Song, Y.; Kou, D.; Deng, Y.; Yang, Y. Stochastic processes regulate belowground community assembly in alpine grasslands on the Tibetan Plateau. Environ. Microbiol. 2021, 24, 179–194. [Google Scholar] [CrossRef]
  92. Sun, H.; Wu, Y.; Zhou, J.; Bing, H.; Zhu, H. Climate influences the alpine soil bacterial communities by regulating the vegetation and the soil properties along an altitudinal gradient in SW China. Catena 2020, 195, 104727. [Google Scholar] [CrossRef]
  93. Yang, Y.; Wu, P. Soil bacterial community varies but fungal community stabilizes along five vertical climate zones. Catena 2020, 195, 104841. [Google Scholar] [CrossRef]
  94. Xu, T.; Shen, Y.; Ding, Z.; Zhu, B. Seasonal dynamics of microbial communities in rhizosphere and bulk soils of two temperate forests. Rhizosphere 2023, 25, 100673. [Google Scholar] [CrossRef]
Microorganisms 12 00728 g001
Figure 1. Soil physicochemical characteristics of the topsoil layer sampled from four selected forest types. (A) soil moisture, (B) pH, (C) soil organic carbon content, (D) total nitrogen content, (E) total phosphorus, (F) available nitrogen content, (G) available phosphorus content, (H) soil organic carbon/total nitrogen, (I) total nitrogen/total phosphorus, and (J) soil organic carbon/total phosphorus. Points represent the individual measures of the sequenced microbial communities, whereas the boxes span the interquartile range (IQR) of all measures, with the middle line being the median of the calculated indices. Different lowercase letters indicate a significant difference between forest types at p = 0.05. If a nonparametric test was used, the differences between groups were further examined using the Benjamini-Hochberg (BH) method to adjust the p-value. Otherwise, the differences between groups were further examined using the Tukey’s Honestly Significant Difference (HSD) test. The organ, blue, cyan and green represented natural deciduous broadleaved and evergreen broadleaved mixed forest, secondary (Sec) forest, evergreen broadleaf (Bro) forest and coniferous (Con) forest, respectively.
Figure 1. Soil physicochemical characteristics of the topsoil layer sampled from four selected forest types. (A) soil moisture, (B) pH, (C) soil organic carbon content, (D) total nitrogen content, (E) total phosphorus, (F) available nitrogen content, (G) available phosphorus content, (H) soil organic carbon/total nitrogen, (I) total nitrogen/total phosphorus, and (J) soil organic carbon/total phosphorus. Points represent the individual measures of the sequenced microbial communities, whereas the boxes span the interquartile range (IQR) of all measures, with the middle line being the median of the calculated indices. Different lowercase letters indicate a significant difference between forest types at p = 0.05. If a nonparametric test was used, the differences between groups were further examined using the Benjamini-Hochberg (BH) method to adjust the p-value. Otherwise, the differences between groups were further examined using the Tukey’s Honestly Significant Difference (HSD) test. The organ, blue, cyan and green represented natural deciduous broadleaved and evergreen broadleaved mixed forest, secondary (Sec) forest, evergreen broadleaf (Bro) forest and coniferous (Con) forest, respectively.
Microorganisms 12 00728 g001
Microorganisms 12 00728 g002
Figure 2. Venn diagram of the soil bacterial community of the evergreen broadleaved forest (Bro), coniferous forest (Con), natural deciduous broadleaved and evergreen broadleaved mixed forest (Mix), and secondary forest (Sec). The numbers within circles represent the specific bacterial amplicon sequence variants (ASVs) in those habitats, and the numbers in the overlaps represent the common bacterial ASVs between habitats. The core number represents the common bacterial ASVs present in all habitats.
Figure 2. Venn diagram of the soil bacterial community of the evergreen broadleaved forest (Bro), coniferous forest (Con), natural deciduous broadleaved and evergreen broadleaved mixed forest (Mix), and secondary forest (Sec). The numbers within circles represent the specific bacterial amplicon sequence variants (ASVs) in those habitats, and the numbers in the overlaps represent the common bacterial ASVs between habitats. The core number represents the common bacterial ASVs present in all habitats.
Microorganisms 12 00728 g002
Microorganisms 12 00728 g003
Figure 3. Relative abundances (%) of the main soil bacterial phyla. Low abundance phyla with less than 0.5% of the total sequences across all samples are grouped into “Others”. The colors of the bars represent the bacterial phyla, and the length represents the relative abundance of the bacterial phyla in the corresponding sample of evergreen broadleaved forest (Bro), coniferous forest (Con), natural deciduous broadleaved and evergreen broadleaved mixed forest (Mix), and secondary deciduous broadleaved and coniferous mixed forest (Sec).
Figure 3. Relative abundances (%) of the main soil bacterial phyla. Low abundance phyla with less than 0.5% of the total sequences across all samples are grouped into “Others”. The colors of the bars represent the bacterial phyla, and the length represents the relative abundance of the bacterial phyla in the corresponding sample of evergreen broadleaved forest (Bro), coniferous forest (Con), natural deciduous broadleaved and evergreen broadleaved mixed forest (Mix), and secondary deciduous broadleaved and coniferous mixed forest (Sec).
Microorganisms 12 00728 g003
Microorganisms 12 00728 g004
Figure 4. Bacterial richness, alpha diversity indices, and their correlations with the soil physicochemical properties of the topsoil layer sampled from four different forest types. Richness is estimated by the abundance-based coverage estimator (ACE), which is a nonparametric richness estimator based on the distribution of abundant (>10) and rare (<10) ASVs, and the Chao1 estimator of richness, which is a nonparametric richness estimator based on the distribution of singletons and doubletons. (A) Chao 1 index, (B) Faith-pd index, (C) observed-ASVs, (D) Shannon index, (E) Simpson index, (F) ACE index, (G) Pielou-e index, and (H) spearman correlations between soil bacterial community diversity and soil physicochemical properties. Different lowercase letters indicate a significant difference between forest types at p = 0.05. The organ, blue, cyan and green represented natural deciduous broadleaved and evergreen broadleaved mixed forest, secondary (Sec) forest, evergreen broadleaf (Bro) forest and coniferous (Con) forest, respectively. *, p ≤ 0.05; **, and p < 0.01.
Figure 4. Bacterial richness, alpha diversity indices, and their correlations with the soil physicochemical properties of the topsoil layer sampled from four different forest types. Richness is estimated by the abundance-based coverage estimator (ACE), which is a nonparametric richness estimator based on the distribution of abundant (>10) and rare (<10) ASVs, and the Chao1 estimator of richness, which is a nonparametric richness estimator based on the distribution of singletons and doubletons. (A) Chao 1 index, (B) Faith-pd index, (C) observed-ASVs, (D) Shannon index, (E) Simpson index, (F) ACE index, (G) Pielou-e index, and (H) spearman correlations between soil bacterial community diversity and soil physicochemical properties. Different lowercase letters indicate a significant difference between forest types at p = 0.05. The organ, blue, cyan and green represented natural deciduous broadleaved and evergreen broadleaved mixed forest, secondary (Sec) forest, evergreen broadleaf (Bro) forest and coniferous (Con) forest, respectively. *, p ≤ 0.05; **, and p < 0.01.
Microorganisms 12 00728 g004
Microorganisms 12 00728 g005
Figure 5. Comparisons of soil bacterial communities by non-metric multidimensional scaling (NMDS) plots based on Bray–Curtis similarity in evergreen broadleaved forest (Bro), coniferous forest (Con), natural deciduous broadleaved and evergreen broadleaved mixed forest (Mix), and secondary forest (Sec).
Figure 5. Comparisons of soil bacterial communities by non-metric multidimensional scaling (NMDS) plots based on Bray–Curtis similarity in evergreen broadleaved forest (Bro), coniferous forest (Con), natural deciduous broadleaved and evergreen broadleaved mixed forest (Mix), and secondary forest (Sec).
Microorganisms 12 00728 g005
Microorganisms 12 00728 g006
Figure 6. Biplot of site scores for the first two redundancy analysis (RDA) axes showing the associations between soil physicochemical properties and the Bray–Curtis distance of soil bacterial communities. The values of RDA axes 1 and 2 are the percentages explained by the corresponding axis. TN is total nitrogen, SWC is soil water content, SOC is total soil organic carbon, C/N is the ratio of soil total organic carbon to total nitrogen, C/P is the ratio of soil total organic carbon to total phosphorus, and N/P is the ratio of total nitrogen to total phosphorus. Bro is evergreen broadleaved forest, Con is coniferous forest, Mix is natural deciduous broadleaved and evergreen broadleaved mixed forest, and Sec is secondary forest.
Figure 6. Biplot of site scores for the first two redundancy analysis (RDA) axes showing the associations between soil physicochemical properties and the Bray–Curtis distance of soil bacterial communities. The values of RDA axes 1 and 2 are the percentages explained by the corresponding axis. TN is total nitrogen, SWC is soil water content, SOC is total soil organic carbon, C/N is the ratio of soil total organic carbon to total nitrogen, C/P is the ratio of soil total organic carbon to total phosphorus, and N/P is the ratio of total nitrogen to total phosphorus. Bro is evergreen broadleaved forest, Con is coniferous forest, Mix is natural deciduous broadleaved and evergreen broadleaved mixed forest, and Sec is secondary forest.
Microorganisms 12 00728 g006
Microorganisms 12 00728 g007
Figure 7. Co-occurrence network structure (a) and network robustness (b) of soil bacterial communities from evergreen broadleaved forest (Bro), coniferous forest (Con), natural deciduous broadleaved and evergreen broadleaved mixed forest (Mix), and secondary forest (Sec). Nodes are colored according to microbial phylum, and the nodes with a larger size show the potential keystone genera (correlation coefficient ≥ 0.6, p < 0.01). Circle and square node shapes represent bacteria and genera, respectively. Edges indicate correlations among nodes; the red and blue edges represent positive and negative correlations, respectively.
Figure 7. Co-occurrence network structure (a) and network robustness (b) of soil bacterial communities from evergreen broadleaved forest (Bro), coniferous forest (Con), natural deciduous broadleaved and evergreen broadleaved mixed forest (Mix), and secondary forest (Sec). Nodes are colored according to microbial phylum, and the nodes with a larger size show the potential keystone genera (correlation coefficient ≥ 0.6, p < 0.01). Circle and square node shapes represent bacteria and genera, respectively. Edges indicate correlations among nodes; the red and blue edges represent positive and negative correlations, respectively.
Microorganisms 12 00728 g007
Microorganisms 12 00728 g008
Figure 8. The niche width (a), the relative importance of different ecological processes based on the null model (b,c), and the relative abundance-frequency relationships (d) of the soil bacterial community under different forest types. Values assigned with the same letter were not significantly different by the post hoc Tukey honestly significant difference (HSD) test (p ≤ 0.05). The dashed blue line represents the 95% confidence interval above and below the prediction (the solid blue line). R2 indicates the coefficient of the neutral fit, and Nm indicates the meta-community size times immigration.
Figure 8. The niche width (a), the relative importance of different ecological processes based on the null model (b,c), and the relative abundance-frequency relationships (d) of the soil bacterial community under different forest types. Values assigned with the same letter were not significantly different by the post hoc Tukey honestly significant difference (HSD) test (p ≤ 0.05). The dashed blue line represents the 95% confidence interval above and below the prediction (the solid blue line). R2 indicates the coefficient of the neutral fit, and Nm indicates the meta-community size times immigration.
Microorganisms 12 00728 g008
Table 1. The overall bacterial community structure comparison using three nonparametric statistical methods.
Table 1. The overall bacterial community structure comparison using three nonparametric statistical methods.
ObjectAdonisANOSIMMRPP
FR2Adjusted-pRAdjusted-pAδAdjusted-p
Among Group2.95590.21700.010.33640.010.07960.56960.001
Mix vs. Sec2.56660.13820.0120.26850.0140.04890.63450.009
Mix vs. Con1.27140.11950.0230.20030.0260.03870.59570.020
Mix vs. Bro2.46310.13340.0060.28570.0080.04390.62950.008
Sec vs. Con2.79650.14880.0240.28500.0190.04750.56410.022
Sec vs. Bro4.83520.23210.0060.66020.0060.09770.62790.003
Con vs. Bro3.05450.16030.0060.38580.0060.05810.56880.003
Note: Adonis, nonparametric multivariate analysis of variance (MANOVA) with the adonis function; ANOSIM, analysis of similarity; MRPP, multi-response permutation procedure. Bro, evergreen broadleaved forest; Con, coniferous forest; Mix, natural deciduous broadleaved and evergreen broadleaved mixed forest; Sec, secondary forest.
Table 2. Topological characteristics of the co-occurrence network of soil bacterial communities in soils with different dominant species.
Table 2. Topological characteristics of the co-occurrence network of soil bacterial communities in soils with different dominant species.
TreatmentsNodesLinksPositive EdgesNegative EdgesAPL ACCAverage DegreeDiameterDensityModularity
Mix869546447327325.02370.471712.575412.96970.01452.4978
Sec8268420596420955.17280.467914.440715.78530.01751.0741
Bro5517375541836.82420.420813.435413.43540.00490.6006
Con5507075121955.24400.442117.223917.22390.00470.4483
Note: Bro, evergreen broadleaved forest; Con, coniferous forest; Mix, natural deciduous broadleaved and evergreen broadleaved mixed forest; Sec, secondary deciduous broadleaved and coniferous mixed forest. ACC, average clustering coefficient; and APL, average path length.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, B.; Feng, W.; Zhou, W.; He, K.; Yang, Z. Association between Soil Physicochemical Properties and Bacterial Community Structure in Diverse Forest Ecosystems. Microorganisms 2024, 12, 728. https://doi.org/10.3390/microorganisms12040728

AMA Style

Yang B, Feng W, Zhou W, He K, Yang Z. Association between Soil Physicochemical Properties and Bacterial Community Structure in Diverse Forest Ecosystems. Microorganisms. 2024; 12(4):728. https://doi.org/10.3390/microorganisms12040728

Chicago/Turabian Style

Yang, Bing, Wanju Feng, Wenjia Zhou, Ke He, and Zhisong Yang. 2024. "Association between Soil Physicochemical Properties and Bacterial Community Structure in Diverse Forest Ecosystems" Microorganisms 12, no. 4: 728. https://doi.org/10.3390/microorganisms12040728

APA Style

Yang, B., Feng, W., Zhou, W., He, K., & Yang, Z. (2024). Association between Soil Physicochemical Properties and Bacterial Community Structure in Diverse Forest Ecosystems. Microorganisms, 12(4), 728. https://doi.org/10.3390/microorganisms12040728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop