Next Article in Journal
Assessing the Ecological Value of an Urban Forest Park: A Case Study of Sinhua Forest Park in Taiwan
Previous Article in Journal
Development of Genomic SSR for the Subtropical Hardwood Tree Dalbergia hupeana and Assessment of Their Transferability to Other Related Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 12
Issue 6
10.3390/f12060805
forests-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 thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Soil Properties and Plant Diversity on Soil Microbial Community Composition and Diversity during Secondary Succession

by
Ruiguang Shang
1,2,
Shuaifeng Li
1,2,
Xiaobo Huang
1,2,
Wande Liu
1,2,
Xuedong Lang
1,2 and
Jianrong Su
1,2,*
1
Research Institute of Resource Insects, Chinese Academy of Forestry, Kunming 650224, China
2
Observation and Research Station of Pu’er Forest Ecosystem, National Forestry and Grassland Administration of China, Pu’er 665007, China
*
Author to whom correspondence should be addressed.
Forests 2021, 12(6), 805; https://doi.org/10.3390/f12060805
Submission received: 18 May 2021 / Revised: 11 June 2021 / Accepted: 15 June 2021 / Published: 18 June 2021
Browse Figures
Versions Notes

Abstract

:
Soil microbial communities play an important role in maintaining the ecosystem during forest secondary succession. However, the underlying mechanisms that drive change in soil microbial community structures during secondary succession remain poorly defined in species-rich subtropical coniferous forests. In this study, Illumina high-throughput sequencing was used to analyze the variations in soil microbial community structures during forest secondary succession in subtropical coniferous forests in China. The role of soil properties and plant diversity in affecting soil bacterial and fungal communities was determined using random forest and structural equation models. Highly variable soil microbial diversity was observed in different stages of secondary succession. Bacterial community diversity rose from early to middle and late successional stages, whereas fungal community diversity increased from early to middle successional stages and then declined in the late stage. The relative abundance of Acidobacteria, Gemmatimonadetes, Eremiobacterota (WPS-2), Rokubacteria, and Mortierellomycota increased during succession, whereas the relative abundance of Ascomycota and Mucoromycota decreased. The community composition and diversity of the soil microbial community were remarkably influenced by plant diversity and soil properties. Notably, tree species richness (TSR) displayed a significant and direct correlation to the composition and diversity of both bacterial and fungal communities. The carbon-to-nitrogen (C:N) ratio had a direct impact on the bacterial community composition and diversity, and pH had a marked impact on the fungal community composition and diversity. Furthermore, succession stage and plant diversity indirectly impacted the composition and diversity of soil bacterial and fungal communities via soil properties. Overall, it can be concluded that soil intrinsic properties and plant diversity might jointly drive the changes in soil microbial community composition and diversity during secondary succession of subtropical coniferous forests.

1. Introduction

Secondary succession is becoming one of the research focuses of terrestrial ecosystems with global land use changes [1]. Forest succession is driven by aboveground–belowground interactions [2,3]. Because soil microbes have essential roles in key ecosystem processes [4,5,6], such as plant productivity [7], nutrient cycling, water regulation, and litter decomposition [8], they are increasingly regarded as vital edaphic variables affecting secondary succession [9,10]. To clarify the roles of soil microbial communities in secondary succession, it is necessary to identify the factors that drive the changes in soil microbial community structures in the successional process [11]. Therefore, how the soil microbial community structure varies during long-term succession needs to be determined.
During ecological succession, soil microbial communities can vary dramatically [12,13,14]. For example, microbial communities change substantially in the transition from grassland to pioneer forest, but less so in the transition from pioneer stage to later stage [15]. Liu et al. [10] found linear increases in the diversity of bacteria and fungi along secondary succession, whereas the abundance increased in the early stages and then plateaued. There are major shifts in the composition of both bacterial and fungal communities during succession in a tropical forest [16], semiarid abandoned farmland [17], grassland [18,19], and a desert steppe [20]. However, it is difficult to identify consistent patterns for the changes in soil microbial communities during succession because of environmental heterogeneity. In addition, whether some microbial species can be regarded as biomarkers for secondary succession remains unclear [10]. Therefore, because of the large heterogeneity among subtropical forest types, gaps remain in the understanding of changes in soil microbial communities during secondary succession in those forests.
Soil microbial community composition and diversity are affected by multiple biotic and abiotic factors, such as forest succession, plant diversity, carbon:nitrogen (C:N) ratio, pH, nutrient availability, and land use type [21,22,23,24]. In one study, soil properties rather than the plant community played the decisive part in shaping the soil microbial community [25], whereas in others, the plant community was decisive [26,27]. However, there is accumulating evidence that soil abiotic properties and plant community characteristics are simultaneously important in regulating soil microbial communities [2,28,29]. The composition and diversity of tree species can directly affect the functions of soil microbial communities by altering the quantity and quality of resources [23,30,31]. It is thus necessary to find out the relations between the composition and diversity of soil bacterial and fungal communities, and plant diversity and soil properties, in order to predict the changes in soil microbial community structures during ecosystem development [8]. Accordingly, the relations between vegetation and soil biotic and abiotic properties during long-term succession need to be studied further [10].
In this study, to identify the key factors that drive changes in soil microbial communities in different stages of succession, we evaluated the effects of plant diversity and soil properties on the composition and diversity of soil bacterial and fungal communities during forest secondary succession in subtropical China. Three secondary successional stages were examined, including a Pinus kesiya forest, a mixed P. kesiya and evergreen broad-leaved forest, and an evergreen broad-leaved forest. Typically, a P. kesiya forest is the first stage in the restoration of evergreen broad-leaved forests after clear-cutting, and there is a predictable sequence of these forest types in succession [32,33]. The composition and diversity of soil bacterial and fungal communities were measured using high-throughput sequencing. The objectives were: (1) to investigate how the soil bacterial and fungal communities changed during secondary succession in a subtropical forest; (2) to identify the plant and soil factors that shaped the soil bacterial and fungal communities during secondary succession. We hypothesized that: (1) the composition and diversity of soil bacterial and fungal communities would clearly shift during succession, with greater changes occurring in fungal communities than in bacterial communities because of the changes in plant species with succession [27]; (2) because plant diversity and soil abiotic properties are the primary regulators of variations in multifunctionality [34], they would be the drivers of alterations in soil microbial communities during succession.

2. Materials and Methods

2.1. Study Sites

The study was conducted in Pu’er (22°34′ N–22°53′ N, 100°56′ E–101°9′ E) in Yunnan Province in Southwest China. The area has a typical monsoon climate, characterized by distinct dry and rainy seasons [35]. The average annual temperature in the region is 17.7 °C. The mean annual precipitation is 1490 mm, and the precipitation is mainly concentrated from May to October. The soil type is latosol [33]. The climax vegetation community in the area is monsoon evergreen broad-leaved forest [32]. Similar to other regions in subtropical China [36,37], the study area has a distinct successional gradient composed of pine forest (early stage), pine and broad-leaved mixed forest (middle stage), and monsoon evergreen broad-leaved forest (late stage).

2.2. Field Sampling

On the basis of species composition, three secondary forest stands in different successional stages were selected. By evaluating tree rings and consulting historical documents, the age of the three secondary forest communities were approximately 40, 70, and 90 years old. Detailed descriptions of the stands are presented in Table 1 and Table S1. In May 2018, seven 30 m × 30 m plots were established in each successional stage, for a total of 21 plots. Within each plot, all living woody plants with a diameter at breast height ≥ 1 cm were recorded. The plants were identified by species and then numbered. Three 1 m × 1 m subplots were established in each plot to determine the species richness and abundance in the herb layer. During the vegetation survey, five soil samples were collected from the depth of 0–10 cm in each plot using a five-point sampling method and then mixed into one composite sample. All composite soil samples were sieved through a mesh 2 mm in size and then divided into two parts. One part was transported to the lab in a portable refrigerator at −20 °C, and then stored at −80 °C until the next step for molecular analysis. The other part was air dried for physicochemical analyses.

2.3. Soil Physicochemical Properties

Standard methods were used to test soil physicochemical properties [38], including the maximum soil water holding capacity (SWHC), pH, soil organic carbon (SOC), total nitrogen (TN), total phosphorus (TP), total potassium (TK), hydrolyzable N (HN), available P (AP), and available K (AK). The C:N ratio was calculated as a mass ratio of SOC and TN.

2.4. DNA Extraction, Illumina Sequencing, and Data Processing

The TIANamp Soil DNA Kit (Tiangen BiotECH Co., Beijing, China) was used to extract the soil DNA from 0.25 g of fresh soil according to the guidelines. The composition and diversity of soil bacterial and fungal communities were measured using 16S rRNA and ITS genes high-throughput sequencing [39]. The 16S rRNA was amplified by the V3–V4 region with the universal primers 338F (5′-ACTCCTACGGGAGGCAGCA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). Targeting the fungal community, the broad-spectrum primers ITS1 (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2 (5′-GCTGCGTTCTTCATCGATGC-3′) were used. According to standard protocols, an Illumina HiSeq2500 platform (Biomarker Technologies, Beijing, China) was used to conduct amplicon sequencing. UCHIME software was used to identify and remove chimera sequences [2]. USEARCH was used to cluster tags at the 97% similarity level to obtain operational taxonomic units (OTUs) [40], and then taxonomic annotation of OTUs was carried out based on the SILVA (bacteria) and UNITE (fungi) databases [41,42]. The raw sequencing data are available at the NCBI Sequence Read Archive (SRA) under the accession number SRP315258: PRJNA722794.

2.5. Data Analyses

The alpha diversity of soil bacterial and fungal communities in the three succession stages was compared using one-way ANOVA, including the number of OTUs (richness) and the Chao1 and Shannon indexes [43,44]. Additionally, nonmetric multidimensional scaling (NMDS), based on the Bray–Curtis dissimilarity matrices of OTUs of bacterial and fungal communities (Hellinger transformed), was used to test for differences in beta diversity using the “vegan” packages [31,43]. The NMDS analysis had certain reliability when Stress < 0.2. The composition of soil bacterial and fungal communities in each plot was represented by the averaged values of the Bray–Curtis dissimilarity matrices [8]. In addition, analysis of similarities (ANOSIM) was used to determine the differences in soil microbial community composition [25,45]. The effects of plant diversity and soil properties on the composition and diversity of soil bacterial and fungal communities in the different secondary successional stages of the P. kesiya forest were analyzed using random forest and structural equation models (SEMs) in “randomForest” and “lavaan” packages, respectively [39]. First, random forest models measured the relative contributions of successional stage, plant diversity, and the different soil properties to the soil bacterial and fungal richness and community composition to determine the main response variables during the forest secondary succession. Second, an SEM measured the direct and indirect effects of successional stage, plant diversity, and soil properties on the soil bacterial and fungal richness and community composition based on the following theoretical framework: (1) successional stage, plant diversity, and soil properties directly affect the soil bacterial and fungal richness and community composition; (2) successional stage indirectly affects the soil bacterial and fungal richness and community composition by affecting plant diversity and soil properties; (3) plant diversity indirectly affects the soil bacterial and fungal richness and community composition by affecting soil properties. The Shipley test of D-separation based on the “piecewiseSEM” package was used to check the missing paths [46]. Finally, the best-fit SEM was determined using a chi-square test (p > 0.05), root mean square error of approximation (RMSEA < 0.05), and goodness-of-fit index (GFI > 0.95) [4,47]. All analyses were performed in R3.6.2 [48].

3. Results

3.1. Soil Properties and Vegetation Characteristics

Soil properties (Table 1) and the composition of plant community (Table S1) were different in the three stages of succession. The SOC, TN, TP, TK, HN, AK, and SWHC increased significantly during secondary succession (p < 0.05). However, the C:N ratio decreased remarkably by 28.18% from the early stage to the late stage (p < 0.05). The pH decreased significantly from the early to middle stage (p < 0.05), but then remained approximately the same in the late stage. Available P was not significantly different among the succession stages. According to the vegetation data, tree species richness (TSR) was markedly higher in the middle and late stages than in the early succession stage (p < 0.05). However, herb species richness (HSR) was similar among the three succession stages.

3.2. Soil Bacterial and Fungal Community Composition and Diversity

For bacteria, 1139 OTUs were obtained and, for fungi, 667 OTUs were obtained. ANOVA analysis of alpha diversity showed remarkable differences among the three succession stages, except for the Shannon index (Figure 1). The OTU and Chao1 index of the bacterial community increased significantly from the early stage to the middle and late stages of succession (p < 0.05, Figure 1a,b). The OTU and Chao1 index of the fungal community was also significantly higher in the middle and late stages than in the early stage, but the index was notably higher in the middle stage than in the late stage (p < 0.05, Figure 1c,d).
The NMDS analysis indicated that the bacterial and fungal communities were different in the three succession stages (Stress = 0.0956 and 0.1109, respectively) (Figure 2a,c). Similarly, according to the ANOSIM, the composition of both bacterial (R = 0.448, p = 0.001) and fungal (R = 0.425, p = 0.001) communities differed significantly among succession stages (Figure 2b,d).
Representative OTU sequences were aligned to the microbial reference database to obtain the classifications of the corresponding species. Twenty-two bacterial phyla and eight fungal phyla were identified. Generally, Acidobacteria, Proteobacteria, and Actinobacteria played a leading role in the bacterial communities, which accounted for more than 88% of the total bacterial sequences (Figure 3a). Notably, succession stage had no significant effect on the relative abundance of the top seven bacterial phyla (except Acidobacteria), but the relative abundance of the less-abundant phyla, including Gemmatimonadetes, WPS-2, and Rokubacteria, increased significantly from the early to late stage of succession (p < 0.05, Table S2). The dominant fungal phylum was Basidiomycota (58.33%), which was followed by Ascomycota (15.84%), Mortierellomycota (15.55%), and Mucoromycota (3.78%) (Figure 3b). The relative abundance of the dominant fungal phyla (except Basidiomycota) was distinctly influenced by successional stage (p < 0.05, Table S2). The relative abundance of Ascomycota and Mucoromycota decreased strongly and significantly from the early to late successional stages, whereas the abundance of Mortierellomycota increased significantly from the early and middle stages to the late stage (Table S2).

3.3. Effects of Plant Diversity and Soil Properties on the Composition and Diversity of Soil Microbial Communities

Random forest analyses were used for indentifying the soil properties and vegetation characteristics that were important in predicting composition and diversity of the soil microbial communities (Figure 4). In the bacterial communities, random forest models explained 29.20% of the variance in composition and 47.45% in diversity. In the fungal communities, the models explained 12.36% of the variance in composition and 32.39% in diversity. Tree species richness and the C:N ratio were the primary predictors of bacterial community diversity, whereas succession stage, TN, HN, C:N ratio, and TSR were the primary predictors of bacterial community composition. The primary predictors of fungal community diversity were succession stage, TN, HN, SWHC, and SOC, whereas SOC, TN, HN, and AP were the primary predictors of fungal community composition.
The SEM was used to gain deeper insights into the effects of succession stage, soil properties, and vegetation characteristics on composition and diversity, and of soil microbial communities. The best-fit model was selected with χ2 = 0.281; n = 1; p = 0.596; RMSEA < 0.001; and GFI = 1. Multiple abiotic and biotic pathways regulated soil microbial community composition and diversity (Figure 5). For example, TSR had a notable direct influence on the composition and diversity of both bacterial and fungal communities. The C:N ratio had a significant negative direct effect on bacterial community diversity, but a significant positive direct effect on bacterial community composition (Figure 5a,b). Soil pH had a significant positive direct effect on bacterial community diversity, and AK had a significant positive direct effect on bacterial community composition. In fungal communities, soil pH and TP had striking negative direct effects on diversity, whereas TK and AP had significant positive direct effects (Figure 5c). Soil pH and TK had significant positive direct effects on fungal community composition (Figure 5d).
Indirect effects of succession stage, plant diversity, and soil properties on the composition and diversity of bacterial and fungal communities were also identified. For example, HSR indirectly affected the bacterial community composition by regulating AK (Figure 5b), and succession stage indirectly affected soil microbial communities through changes in TSR. In addition, succession stage indirectly affected the soil microbial community composition and diversity via soil pH, TP, TK, AK, and C:N ratio (Figure 5).

4. Discussion

Our results showed forest secondary succession altered the diversity of soil bacterial and fungal communities. This finding is in agreement with those of previous studies in temperate and subtropical forests [31,49]. In this study, the diversity of the soil bacterial community increased from the early to middle and late stages, while fungal community diversity was highest in the middle successional stage. This difference might be due to the different nutritional preferences of bacteria and fungi [10]. The direction of secondary succession was from a Pinus forest to an evergreen broad-leaved forest in the present study [50], and changes were observed in plant species, forest type, and soil fertility. Such changes could provide different environmental conditions for bacteria and fungi [23].
Generally, soil bacteria exhibit faster growth and turnover rates than fungi during early successional stages [51,52]. In P. kesiya forests in the early successional stage, biomass production is relatively high and soil fertility is relatively low [50]. Because these conditions lead to greater competition for available resources and lower plant diversity, the result can be lower soil bacterial and fungal community diversity in early succession, as observed in this study. In addition, the decomposition of conifer species litter is slower than that of evergreen broad-leaved species [53]. Thus, in early succession, the environment did not support high soil microbial diversity. However, because plants and soil fungi usually form symbioses, the increase in plant diversity with succession significantly positively affected soil fungal community diversity. These results provided further evidence that secondary succession could indirectly affect the soil microbial community structure [49].
The findings of this study suggested that TSR was the key driver shaping soil fungal community diversity, with the path coefficient from plant diversity to soil microbial diversity greater for soil fungi than for bacteria (Figure 5). Previous studies have demonstrated that plant communities significantly affect soil microbial communities via litter chemical properties and root exudates [22,23,28]. The diversity of soil bacterial and fungal communities was affected by different successional stages [21,52]. Our results showed that the differences in fungal community diversity between successional stages were greater than those in soil bacterial community diversity, with the highest fungal community diversity in the middle successional stage (Figure 1). This might be attributed to the middle successional stage featuring more plant species with a particular soil fungal guild, such as the ectomycorrhizal fungi specifically associated with P. kesiya, Alnus nepalensis, and Fagaceae species [54,55]. Furthermore, the composition of soil microbial communities may be determined by the dominant plant species [27,56]. A previous study indicated that plant attributes (such as plant cover and functional traits) were also important factors affecting soil microbial communities [57]. In fact, the vegetation of middle successional stages contains more plant species because it includes both early and late successional stages [31]. With more plant species, more fungal species can coexist, which increases soil fungal community diversity. These fungi include the Basidiomycota and Ascomycota, which can easily metabolize the organic substrates in rhizodeposition [52]. Moreover, NMDS analysis showed that the distance in fungal community composition between successional stages was larger than that in bacterial community composition (Figure 2a,c), which suggested that fungi were more sensitive to environmental variations than the bacteria [17]. Taken together, these results confirm our hypothesis that the shift of fungi with succession was more obvious than bacteria.
Soil properties also significantly affected the composition and diversity of soil bacterial and fungal communities. In this study, soil pH and C:N ratio played critical roles in shaping soil microbial communities. Soil pH is a major factor that affects the soil bacterial community structure [8,56]. Extensive studies have shown that forest soil acidification on decadal scales is caused by trees [29]. The relative abundance of Acidobacteria increased during secondary succession. This phylum tends to be more abundant at lower pH levels [23,58] and, therefore, the increase in Acidobacteria could be explained by the decrease in soil pH in the middle and late successional stages. Soil C:N ratio and P concentration are usually considered limiting factors for soil microbial diversity [8,59]. In addition, secondary succession affects soil C and N cycles [56]. The soil C:N ratio significantly affected bacterial community diversity and composition, which suggests that the determinants of the soil bacterial community could be the variables related to soil carbon and nitrogen transformation. Soil fungal community structure is also influenced by soil pH [11,23,60]; similar findings were observed in the present research (Figure 5c,d). TK was found to be a strong factor influencing the fungal community structure, which may be attributed to the easy leaching of potassium from soil [44]. Furthermore, the solubility and availability of potassium could be partially affected by soil fungi [61]. Total P limited soil fungal community diversity, whereas the effect of AP was positive and significant. This may be attributed to Mortierellomycota being able to dissolve phosphorus [62]. We found a significant decrease in the abundance of Ascomycota with the accumulation of soil nutrients, which was in line with the works by Liu et al. [10]. The phenomenon might be explained by the fact that Ascomycota was able to grow in low-nitrogen conditions and had negative correlations with nitrate and nitrite [63]. Taken together, these findings demonstrate that indices of soil quality can be predictors of changes in soil microbial community diversity during pedogenesis [8].

5. Conclusions

In summary, the structure of soil bacterial and fungal communities showed significant differences at the three stages in the secondary succession of a P. kesiya forest, and the changes in plant diversity and those of most soil properties during succession could explain the differences in the microbial communities. In addition, plant diversity was a better predictor of fungal community diversity than bacterial community diversity. The relative abundance of dominant soil microbial phyla could also predict the direction of forest secondary succession, indicating that some phyla could be an important reference for future afforestation management in most subtropical zones in China. Future research should focus on the effects of specific plant species and plant attributes on functional groups of soil microbial communities to provide more precise reference for the management of forest secondary succession.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/f12060805/s1, Table S1: dominant plant species in different stages of secondary succession of a P. kesiya forest, Table S2: relative abundance of phyla in bacterial and fungal communities in different stages of secondary succession of a P. kesiya forest.

Author Contributions

Conceptualization, J.S. and W.L.; formal analysis, R.S. and S.L.; funding acquisition, J.S.; investigation, R.S., S.L., X.H. and X.L.; writing (original draft), R.S.; writing (review and editing), J.S. and S.L. All authors have read and agreed to the published version of the manuscript.

Funding

The Fundamental Research Funds of the Chinese Academy of Forestry (CAFYBB2021ZA002), the Biosafety and Genetic Resource Management Program of the National Forestry and Grassland Administration (KJZXSA2018001, KJZXSA2019048), the National Key R&D Program of China (2017YFC0505005), and the Basic Research Foundation of Yunnan Province (202001AS070005) supported this work.

Acknowledgments

We would like to thank Chonghua Xu and Pengbin Han for their assistance in the field work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, Z.; Wang, C.; Jiang, L.; Luo, Y. Trends in soil microbial communities during secondary succession. Soil Biol. Biochem. 2017, 115, 92–99. [Google Scholar] [CrossRef]
  2. Chai, Y.; Cao, Y.; Yue, M.; Tian, T.; Yin, Q.; Dang, H.; Quan, J.; Zhang, R.; Wang, M. Soil abiotic properties and plant functional traits mediate associations between soil microbial and plant communities during a secondary forest succession on the Loess Plateau. Front. Microbiol. 2019, 10, 895. [Google Scholar] [CrossRef]
  3. Li, J.-J.; Zheng, Y.-M.; Yan, J.-X.; Li, H.-J.; He, J.-Z. Succession of plant and soil microbial communities with restoration of abandoned land in the Loess Plateau, China. J. Soils Sediments 2013, 13, 760–769. [Google Scholar] [CrossRef]
  4. Delgado-Baquerizo, M.; Maestre, F.T.; Reich, P.B.; Jeffries, T.C.; Gaitan, J.J.; Encinar, D.; Berdugo, M.; Campbell, C.D.; Singh, B.K. Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat. Commun. 2016, 7, 10541. [Google Scholar] [CrossRef] [Green Version]
  5. Delgado-Baquerizo, M.; Eldridge, D.J.; Ochoa, V.; Gozalo, B.; Singh, B.K.; Maestre, F.T. Soil microbial communities drive the resistance of ecosystem multifunctionality to global change in drylands across the globe. Ecol. Lett. 2017, 20, 1295–1305. [Google Scholar] [CrossRef] [Green Version]
  6. Fierer, N. Embracing the unknown: Disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 2017, 15, 579–590. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, Q.-L.; Ding, J.; Zhu, Y.-G.; He, J.-Z.; Hu, H.-W. Soil bacterial taxonomic diversity is critical to maintaining the plant productivity. Environ. Int. 2020, 140, 105766. [Google Scholar] [CrossRef]
  8. Delgado-Baquerizo, M.; Bardgett, R.D.; Vitousek, P.M.; Maestre, F.T.; Williams, M.A.; Eldridge, D.J.; Lambers, H.; Neuhauser, S.; Gallardo, A.; García-Velázquez, L.; et al. Changes in belowground biodiversity during ecosystem development. Proc. Natl. Acad. Sci. USA 2019, 116, 6891–6896. [Google Scholar] [CrossRef] [Green Version]
  9. Banning, N.C.; Gleeson, D.B.; Grigg, A.H.; Grant, C.D.; Andersen, G.L.; Brodie, E.L.; Murphy, D.V. Soil microbial community successional patterns during forest ecosystem restoration. Appl. Environ. Microbiol. 2011, 77, 6158–6164. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, J.; Jia, X.; Yan, W.; Zhong, Y.; Shangguan, Z. Changes in soil microbial community structure during long-term secondary succession. Land Degrad. Dev. 2020, 31, 1151–1166. [Google Scholar] [CrossRef]
  11. Hendershot, J.N.; Read, Q.D.; Henning, J.A.; Sanders, N.J.; Classen, A.T. Consistently inconsistent drivers of microbial diversity and abundance at macroecological scales. Ecology 2017, 98, 1757–1763. [Google Scholar] [CrossRef] [Green Version]
  12. Cutler, N.A.; Chaput, D.L.; van der Gast, C.J. Long-term changes in soil microbial communities during primary succession. Soil Biol. Biochem. 2014, 69, 359–370. [Google Scholar] [CrossRef] [Green Version]
  13. Jangid, K.; Whitman, W.B.; Condron, L.M.; Turner, B.L.; Williams, M.A. Soil bacterial community succession during long-term ecosystem development. Mol. Ecol. 2013, 22, 3415–3424. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, L.; Zhu, K.; Krause, S.M.B.; Li, S.; Wang, X.; Zhang, Z.; Shen, M.; Yang, Q.; Lian, J.; Wang, X.; et al. Changes in assembly processes of soil microbial communities during secondary succession in two subtropical forests. Soil Biol. Biochem. 2021, 154, 108144. [Google Scholar] [CrossRef]
  15. Zhong, Y.; Yan, W.; Wang, R.; Wang, W.; Shangguan, Z. Decreased occurrence of carbon cycle functions in microbial communities along with long-term secondary succession. Soil Biol. Biochem. 2018, 123, 207–217. [Google Scholar] [CrossRef]
  16. Smith, A.P.; Marín-Spiotta, E.; Balser, T. Successional and seasonal variations in soil and litter microbial community structure and function during tropical postagricultural forest regeneration: A multiyear study. Glob. Chang. Biol. 2015, 21, 3532–3547. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, J.; Liu, G.; Zhang, C.; Wang, G.; Fang, L.; Cui, Y. Higher temporal turnover of soil fungi than bacteria during long-term secondary succession in a semiarid abandoned farmland. Soil Tillage Res. 2019, 194, 104305. [Google Scholar] [CrossRef]
  18. Kuramae, E.E.; Gamper, H.A.; Yergeau, E.; Piceno, Y.M.; Brodie, E.L.; DeSantis, T.Z.; Andersen, G.L.; van Veen, J.A.; Kowalchuk, G.A. Microbial secondary succession in a chronosequence of chalk grasslands. ISME J. 2010, 4, 711–715. [Google Scholar] [CrossRef]
  19. Lupatini, M.; Suleiman, A.K.A.; Jacques, R.J.S.; Lemos, L.N.; Pylro, V.S.; Van Veen, J.A.; Kuramae, E.E.; Roesch, L.F.W. Moisture is more important than temperature for assembly of both potentially active and whole prokaryotic communities in subtropical grassland. Microb. Ecol. 2019, 77, 460–470. [Google Scholar] [CrossRef]
  20. Fan, M.; Li, J.; Tang, Z.; Shangguan, Z. Soil bacterial community succession during desertification in a desert steppe ecosystem. Land Degrad. Dev. 2020, 31, 1662–1674. [Google Scholar] [CrossRef]
  21. Hanif, M.A.; Guo, Z.; Moniruzzaman, M.; He, D.; Yu, Q.; Rao, X.; Liu, S.; Tan, X.; Shen, W. Plant taxonomic diversity better explains soil fungal and bacterial diversity than functional diversity in restored forest ecosystems. Plants 2019, 8, 479. [Google Scholar] [CrossRef] [Green Version]
  22. Lladó, S.; López-Mondéjar, R.; Baldrian, P. Drivers of microbial community structure in forest soils. Appl. Microbiol. Biotechnol. 2018, 102, 4331–4338. [Google Scholar] [CrossRef]
  23. Nakayama, M.; Imamura, S.; Taniguchi, T.; Tateno, R. Does conversion from natural forest to plantation affect fungal and bacterial biodiversity, community structure, and co-occurrence networks in the organic horizon and mineral soil? For. Ecol. Manag. 2019, 446, 238–250. [Google Scholar] [CrossRef]
  24. Slabbert, E.; Kongor, R.Y.; Esler, K.J.; Jacobs, K. Microbial diversity and community structure in Fynbos soil. Mol. Ecol. 2010, 19, 1031–1041. [Google Scholar] [CrossRef] [PubMed]
  25. Kuramae, E.; Gamper, H.; van Veen, J.; Kowalchuk, G. Soil and plant factors driving the community of soil-borne microorganisms across chronosequences of secondary succession of chalk grasslands with a neutral pH. FEMS Microbiol. Ecol. 2011, 77, 285–294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Krüger, C.; Kohout, P.; Janoušková, M.; Püschel, D.; Frouz, J.; Rydlová, J. Plant communities rather than soil properties structure arbuscular mycorrhizal fungal communities along primary succession on a mine spoil. Front. Microbiol. 2017, 8, 719. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Urbanová, M.; Šnajdr, J.; Baldrian, P. Composition of fungal and bacterial communities in forest litter and soil is largely determined by dominant trees. Soil Biol. Biochem. 2015, 84, 53–64. [Google Scholar] [CrossRef]
  28. Bai, Z.; Wu, X.; Lin, J.-J.; Xie, H.-T.; Yuan, H.-S.; Liang, C. Litter-, soil- and C:N-stoichiometry-associated shifts in fungal communities along a subtropical forest succession. Catena 2019, 178, 350–358. [Google Scholar] [CrossRef]
  29. 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]
  30. Cline, L.C.; Zak, D.R. Soil microbial communities are shaped by plant-driven changes in resource availability during secondary succession. Ecology 2015, 96, 3374–3385. [Google Scholar] [CrossRef]
  31. Li, S.; Huang, X.; Shen, J.; Xu, F.; Su, J. Effects of plant diversity and soil properties on soil fungal community structure with secondary succession in the Pinus yunnanensis forest. Geoderma 2020, 379, 114646. [Google Scholar] [CrossRef]
  32. Liu, W.; Su, J. Successional leaf traits of monsoon evergreen broad-leaved forest, Southwest China. iForest 2017, 10, 391–396. [Google Scholar] [CrossRef] [Green Version]
  33. Liu, W.; Su, J.; Li, S.; Lang, X.; Huang, X. Non-structural carbohydrates regulated by season and species in the subtropical monsoon broad-leaved evergreen forest of Yunnan Province, China. Sci. Rep. 2018, 8, 1083. [Google Scholar] [CrossRef]
  34. Lucas-Borja, M.E.; Delgado-Baquerizo, M. Plant diversity and soil stoichiometry regulates the changes in multifunctionality during pine temperate forest secondary succession. Sci. Total Environ. 2019, 697, 134204. [Google Scholar] [CrossRef]
  35. Li, S.; Su, J.; Lang, X.; Liu, W.; Ou, G. Positive relationship between species richness and aboveground biomass across forest strata in a primary Pinus kesiya forest. Sci. Rep. 2018, 8, 2227. [Google Scholar] [CrossRef]
  36. Huang, W.; Liu, J.; Wang, Y.P.; Zhou, G.; Han, T.; Li, Y. Increasing phosphorus limitation along three successional forests in southern China. Plant Soil 2013, 364, 181–191. [Google Scholar] [CrossRef]
  37. Zhou, G.; Guan, L.; Wei, X.; Zhang, D.; Zhang, Q.; Yan, J.; Wen, D.; Liu, J.; Liu, S.; Huang, Z.; et al. Litterfall production along successional and altitudinal gradients of subtropical monsoon evergreen broadleaved forests in Guangdong, China. Plant Ecol. 2007, 188, 77–89. [Google Scholar] [CrossRef]
  38. Bao, S. Soil and Agricultural Chemistry Analysis; China Agriculture Press: Beijing, China, 2000. [Google Scholar]
  39. Li, S.; Huang, X.; Lang, X.; Shen, J.; Xu, F.; Su, J. Cumulative effects of multiple biodiversity attributes and abiotic factors on ecosystem multifunctionality in the Jinsha River valley of southwestern China. For. Ecol. Manag. 2020, 472, 118281. [Google Scholar] [CrossRef]
  40. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef]
  41. Nilsson, R.H.; Larsson, K.-H.; Taylor, A.F.S.; Bengtsson-Palme, J.; Jeppesen, T.S.; Schigel, D.; Kennedy, P.; Picard, K.; Glöckner, F.O.; Tedersoo, L.; et al. The UNITE database for molecular identification of fungi: Handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 2019, 47, D259–D264. [Google Scholar] [CrossRef] [PubMed]
  42. Zhou, H.; Zhang, D.; Jiang, Z.; Sun, P.; Xiao, H.; Yuxin, W.; Chen, J. Changes in the soil microbial communities of alpine steppe at Qinghai-Tibetan Plateau under different degradation levels. Sci. Total Environ. 2019, 651, 2281–2291. [Google Scholar] [CrossRef]
  43. Montagna, M.; Berruti, A.; Bianciotto, V.; Cremonesi, P.; Giannico, R.; Gusmeroli, F.; Lumini, E.; Pierce, S.; Pizzi, F.; Turri, F.; et al. Differential biodiversity responses between kingdoms (plants, fungi, bacteria and metazoa) along an Alpine succession gradient. Mol. Ecol. 2018, 27, 3671–3685. [Google Scholar] [CrossRef]
  44. Wang, X.; Zhang, Z.; Yu, Z.; Shen, G.; Cheng, H.; Tao, S. Composition and diversity of soil microbial communities in the alpine wetland and alpine forest ecosystems on the Tibetan Plateau. Sci. Total Environ. 2020, 747, 141358. [Google Scholar] [CrossRef]
  45. Huang, F.; Liu, Z.; Mou, H.; Zhang, P.; Jia, Z. Effects of different long-term farmland mulching practices on the loessial soil fungal community in a semiarid region of China. Appl. Soil Ecol. 2019, 137, 111–119. [Google Scholar] [CrossRef]
  46. Shipley, B. The AIC model selection method applied to path analytic models compared using a d-separation test. Ecology 2013, 94, 560–564. [Google Scholar] [CrossRef]
  47. Rosseel, Y. lavaan: An R package for structural equation modeling. J. Stat. Softw. 2012, 48, 1–36. [Google Scholar] [CrossRef] [Green Version]
  48. R Core Computing Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: https://www.R-project.org/ (accessed on 22 March 2021).
  49. Yan, B.; Sun, L.; Li, J.; Liang, C.; Wei, F.; Xue, S.; Wang, G. Change in composition and potential functional genes of soil bacterial and fungal communities with secondary succession in Quercus liaotungensis forests of the Loess Plateau, western China. Geoderma 2020, 364, 114199. [Google Scholar] [CrossRef]
  50. Li, S.; Liu, W.; Lang, X.; Huang, X.; Su, J. Species richness, not abundance, drives ecosystem multifunctionality in a subtropical coniferous forest. Ecol. Indic. 2021, 120, 106911. [Google Scholar] [CrossRef]
  51. Zechmeister-Boltenstern, S.; Keiblinger, K.M.; Mooshammer, M.; Peñuelas, J.; Richter, A.; Sardans, J.; Wanek, W. The application of ecological stoichiometry to plant–microbial–soil organic matter transformations. Ecol. Monogr. 2015, 85, 133–155. [Google Scholar] [CrossRef] [Green Version]
  52. Zhong, Z.; Zhang, X.; Wang, X.; Fu, S.; Wu, S.; Lu, X.; Ren, C.; Han, X.; Yang, G. Soil bacteria and fungi respond differently to plant diversity and plant family composition during the secondary succession of abandoned farmland on the Loess Plateau, China. Plant Soil 2020, 448, 183–200. [Google Scholar] [CrossRef]
  53. Li, C.; Shi, L.-L.; Ostermann, A.; Xu, J.; Li, Y.; Mortimer, P.E. Indigenous trees restore soil microbial biomass at faster rates than exotic species. Plant Soil 2015, 396, 151–161. [Google Scholar] [CrossRef]
  54. Frąc, M.; Hannula, S.E.; Bełka, M.; Jędryczka, M. Fungal biodiversity and their role in soil health. Front. Microbiol. 2018, 9, 707. [Google Scholar] [CrossRef] [Green Version]
  55. Goldmann, K.; Schöning, I.; Buscot, F.; Wubet, T. Forest management type influences diversity and community composition of soil fungi across temperate forest ecosystems. Front. Microbiol. 2015, 6, 1300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. 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] [Green Version]
  57. Delgado-Baquerizo, M.; Fry, E.L.; Eldridge, D.J.; de Vries, F.T.; Manning, P.; Hamonts, K.; Kattge, J.; Boenisch, G.; Singh, B.K.; Bardgett, R.D. Plant attributes explain the distribution of soil microbial communities in two contrasting regions of the globe. New Phytol. 2018, 219, 574–587. [Google Scholar] [CrossRef] [Green Version]
  58. Cui, J.; Wang, J.; Xu, J.; Xu, C.; Xu, X. Changes in soil bacterial communities in an evergreen broad-leaved forest in east China following 4 years of nitrogen addition. J. Soils Sediments 2017, 17, 2156–2164. [Google Scholar] [CrossRef]
  59. Cao, Y.; Fu, S.; Zou, X.; Cao, H.; Shao, Y.; Zhou, L. Soil microbial community composition under Eucalyptus plantations of different age in subtropical China. Eur. J. Soil Biol. 2010, 46, 128–135. [Google Scholar] [CrossRef]
  60. Chen, W.; Wang, J.; Meng, Z.; Xu, R.; Chen, J.; Zhang, Y.; Hu, T. Fertility-related interplay between fungal guilds underlies plant richness–productivity relationships in natural grasslands. New Phytol. 2020, 226, 1129–1143. [Google Scholar] [CrossRef] [PubMed]
  61. Benito, B.; González-Guerrero, M. Unravelling potassium nutrition in ectomycorrhizal associations. New Phytol. 2014, 201, 707–709. [Google Scholar] [CrossRef]
  62. Grzadziel, J.; Galazka, A. Fungal biodiversity of the most common types of polish soil in a long-term microplot experiment. Front. Microbiol. 2019, 10, 6. [Google Scholar] [CrossRef] [Green Version]
  63. Wang, Y.; Sen, K.; He, Y.; Xie, Y.; Wang, G. Impact of environmental gradients on the abundance and diversity of planktonic fungi across coastal habitats of contrasting trophic status. Sci. Total Environ. 2019, 683, 822–833. [Google Scholar] [CrossRef] [PubMed]
Forests 12 00805 g001 550
Figure 1. Alpha diversity of soil bacterial (a,b) and fungal (c,d) communities in different stages of secondary succession in a P. kesiya forest. Different letters above the box plots indicate significant differences (p < 0.05) between succession stages.
Figure 1. Alpha diversity of soil bacterial (a,b) and fungal (c,d) communities in different stages of secondary succession in a P. kesiya forest. Different letters above the box plots indicate significant differences (p < 0.05) between succession stages.
Forests 12 00805 g001
Forests 12 00805 g002 550
Figure 2. Nonmetric multidimensional scaling (NMDS) (a,c) and analysis of similarities (b,d) based on Bray–Curtis dissimilarities of soil bacterial (a,b) and fungal (c,d) communities during secondary succession in a P. kesiya forest.
Figure 2. Nonmetric multidimensional scaling (NMDS) (a,c) and analysis of similarities (b,d) based on Bray–Curtis dissimilarities of soil bacterial (a,b) and fungal (c,d) communities during secondary succession in a P. kesiya forest.
Forests 12 00805 g002
Forests 12 00805 g003 550
Figure 3. Relative abundance (%) of the phyla in soil bacterial (a) and fungal (b) communities in different stages of secondary succession in a P. kesiya forest.
Figure 3. Relative abundance (%) of the phyla in soil bacterial (a) and fungal (b) communities in different stages of secondary succession in a P. kesiya forest.
Forests 12 00805 g003
Forests 12 00805 g004 550
Figure 4. Main predictors of the composition and diversity of soil bacterial and fungal communities based on random forest models. Abbreviations: SWHC, soil water holding capacity; SOC, soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium; HN, hydrolyzable nitrogen; AP, available phosphorus; AK, available potassium; TSR, tree species richness; HSR, herb species richness.
Figure 4. Main predictors of the composition and diversity of soil bacterial and fungal communities based on random forest models. Abbreviations: SWHC, soil water holding capacity; SOC, soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium; HN, hydrolyzable nitrogen; AP, available phosphorus; AK, available potassium; TSR, tree species richness; HSR, herb species richness.
Forests 12 00805 g004
Forests 12 00805 g005 550
Figure 5. Effects of plant diversity, succession stage, and soil properties on the diversity (a,c) and composition (b,d) of soil bacterial (a,b) and fungal (c,d) communities based on a structural equation model. Solid lines represent significant paths (* p < 0.05; ** p < 0.01); dashed lines represent nonsignificant paths (p ≥ 0.05). Abbreviations: SWHC, soil water holding capacity; SOC, soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium; HN, hydrolyzable nitrogen; AP, available phosphorus; AK, available potassium; TSR, tree species richness; HSR, herb species richness.
Figure 5. Effects of plant diversity, succession stage, and soil properties on the diversity (a,c) and composition (b,d) of soil bacterial (a,b) and fungal (c,d) communities based on a structural equation model. Solid lines represent significant paths (* p < 0.05; ** p < 0.01); dashed lines represent nonsignificant paths (p ≥ 0.05). Abbreviations: SWHC, soil water holding capacity; SOC, soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium; HN, hydrolyzable nitrogen; AP, available phosphorus; AK, available potassium; TSR, tree species richness; HSR, herb species richness.
Forests 12 00805 g005
Table
Table 1. Soil properties and stand attributes in different stages of secondary succession of a P. kesiya forest.
Table 1. Soil properties and stand attributes in different stages of secondary succession of a P. kesiya forest.
PropertiesStages of Succession
EarlyMiddleLate
SOC (g·kg−1)18.99 ± 1.83 c28.80 ± 2.09 b43.81 ± 2.39 a
TN (g·kg−1)0.90 ± 0.07 c1.65 ± 0.10 b2.86 ± 0.13 a
TP (g·kg−1)0.19 ± 0.02 c0.27 ± 0.02 b0.39 ± 0.03 a
TK (g·kg−1)4.37 ± 0.79 b6.90 ± 1.01 ab8.59 ± 1.21 a
C:N ratio21.07 ± 1.02 a17.52 ± 0.80 b15.29 ± 0.43 b
HN (mg·kg−1)80.47 ± 6.08 c139.99 ± 5.85 b246.64 ± 24.35 a
AP (mg·kg−1)7.35 ± 1.41 a9.65 ± 1.09 a9.21 ± 2.74 a
AK (mg·kg−1)91.63 ± 14.51 b108.19 ± 11.36 ab140.87 ± 9.18 a
pH4.28 ± 0.10 a4.07 ± 0.04 b4.10 ± 0.04 ab
SWHC (%)32.80 ± 1.54 c41.01 ± 1.25 b52.22 ± 3.34 a
TSR26 ± 1 b34 ± 2 a34 ± 2 a
HSR13 ± 1 a14 ± 1 a14 ± 1 a
Different lowercase letters in the same row indicate significant differences at p < 0.05. Abbreviations: SOC, soil organic carbon; TN, total nitrogen; TP, total phosphorus; TK, total potassium; HN, hydrolyzable nitrogen; AP, available phosphorus; AK, available potassium; SWHC, soil water holding capacity; TSR, tree species richness; HSR, herb species richness.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shang, R.; Li, S.; Huang, X.; Liu, W.; Lang, X.; Su, J. Effects of Soil Properties and Plant Diversity on Soil Microbial Community Composition and Diversity during Secondary Succession. Forests 2021, 12, 805. https://doi.org/10.3390/f12060805

AMA Style

Shang R, Li S, Huang X, Liu W, Lang X, Su J. Effects of Soil Properties and Plant Diversity on Soil Microbial Community Composition and Diversity during Secondary Succession. Forests. 2021; 12(6):805. https://doi.org/10.3390/f12060805

Chicago/Turabian Style

Shang, Ruiguang, Shuaifeng Li, Xiaobo Huang, Wande Liu, Xuedong Lang, and Jianrong Su. 2021. "Effects of Soil Properties and Plant Diversity on Soil Microbial Community Composition and Diversity during Secondary Succession" Forests 12, no. 6: 805. https://doi.org/10.3390/f12060805

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