Kai-Li Liu
Bo-Yao Chen2
Chun-Sheng Wang- 1College of Forestry, Central South University of Forestry & Technology, Changsha, China
- 2Research Institute of Tropical Forestry, Chinese Academy of Forestry, Guangzhou, China
- 3National Long-term Scientific Research Base of Central and Subtropical Forestry, Changsha, China
Introduction: Soil microorganisms are the key factors in elucidating the effects of thinning on tree growth performance, but the effects of vegetation and soil on the species composition and function of soil microorganisms after thinning are still not well elaborated.
Methods: The effects of thinning on understory vegetation diversity, soil physicochemical properties and soil microbial community composition were investigated in a thinning trial plantation of Cryptomeria japonica var. sinensis, including four thinning intensities (control: 0%, LIT: 20%, MIT: 30% and HIT: 40%), and the relationships of the microbial community structure with the understory vegetation diversity and soil properties were assessed.
Results: The results showed that thinning had a greater effect on the diversity of the shrub layer than the herb layer. The soil bulk density and the contents of soil organic matter, total potassium and nitrogen increased with increasing thinning intensities. The Shannon and Chao indices of soil bacteria and fungi were significantly lower in the LIT, MIT and HIT treatments than in the control. Thinning can significantly increase the abundance of Proteobacteria and Actinobacteria, and higher thinning intensities led to a higher relative abundance of Ascomycota and a lower relative abundance of Basidiomycota, Rozellomycota, and Mortierellomycota. Redundancy analysis indicated that soil physicochemical properties rather than understory vegetation diversity were the main drivers of microbial communities, and fungi were more sensitive to soil properties than bacteria. Functional prediction showed that thinning significantly reduced the potential risk of human diseases and plant pathogens, and the nitrogen fixation capacity of bacteria was the highest in the HIT treatment. Thinning significantly increased the relative abundance of cellulolysis and soil saprotrophs in bacteria and fungi.
Conclusion: The findings provide important insights into the effects of thinning on C. japonica var. sinensis plantation ecosystems, which is essential for developing thinning strategies to promote their ecological and economic benefits.
1. Introduction
Soil microorganisms play a crucial role in nutrient decomposition, absorption, and transformation in terrestrial ecosystems (Soong et al., 2019; Tian et al., 2019; Averill et al., 2021). The diversity and abundance of soil bacteria and fungi not only affect nutrient mineralization and cycling but are also sensitive to ecological environmental conditions, such as soil properties (Stark et al., 2012; Shen et al., 2013; Urbanová et al., 2015; Lladó et al., 2018), vegetation diversity (Liu et al., 2021; Shang et al., 2021), and climatic conditions (Bona et al., 2021). Soil properties are considered crucial factors influencing soil microbial activity and community structure, which have become a research hotspot in the field of the soil microbiome. For example, pH is frequently considered the main driver of microbial community structure in soil (Lammel et al., 2018; Ni et al., 2021; Queiroz et al., 2021); soil moisture strongly affects the fate of soil nutrient diffusion (Lei et al., 2022), and is also an important factor influencing the soil microbial community (Goransson et al., 2013; Lei et al., 2021); and the composition and diversity index of bacteria are significantly correlated with the contents of total nitrogen, available nitrogen and available phosphorus in soil (Zhang et al., 2016). However, soil nutrients are one of the factors influencing vegetation growth. In the process of nutrient recycling, nutrients are taken up and utilized, partly retained in plants, and partly returned to the soil as litter; and then, nutrients are released by the decomposition of litter materials (Osman, 2013). Soil organic matters mainly come from litter, and their decomposition can be an important energy source for soil microorganisms (Yang et al., 2017; Phillips et al., 2019; Xu et al., 2020). Therefore, vegetation may also play an important role in changes in the soil microbial community.
As a traditional technique in forest management, thinning not only has direct effects on light transmittance and microclimate in forests (Son et al., 2021) but also indirectly affects understory vegetation diversity (Ares et al., 2009) and soil physicochemical properties (Wic Baena et al., 2013). Thinning practices influence the understory vegetation diversity, and affect the functions of the soil microorganisms (Dang et al., 2018). It has been proven that the diversity and community structure of aboveground plants and belowground microbial communities in forest ecosystems are closely correlated (Chen et al., 2015; Zhou et al., 2020). To date, the relationship between soil physicochemical characteristics and vegetation has been widely studied. However, research on the relationships of soil microorganisms with understory vegetation diversity and soil properties under thinning is still scarce.
Cryptomeria japonica var. sinensis is an excellent fast-growing timber species and is one of the main plantation species in the high-altitude subtropical areas of China (Zhang et al., 2021). In particular, the planting area of C. japonica var. sinensis has been nearly 200,000 hectares in the western parts of Sichuan Province, which is located in the upper reaches of the Yangtze River. These C. japonica var. sinensis plantations are important ecological shields of the Yangtze River, and play an important role in economic development and environmental protection in this region. Due to the high density of monoculture in the initial planting, low stand quality and ecological benefits occur in the C. japonica var. sinensis forest.
To improve stand quality and promote the ecological benefits of these plantations, we arranged a thinning trial in 2014 and investigated the understory vegetation diversity, soil physicochemical properties and microbial community compositions of the trial plantations 5 years after thinning. We hypothesized that thinning of C. japonica var. sinensis plantations would affect the soil microbial composition and diversity, which would be driven by the changes in understory vegetation diversity and soil properties. This study aimed to specifically address the following: (1) to understand the changes in understory vegetation diversity and soil properties after thinning at different intensities; (2) to explore the differences in soil microorganism composition and function under thinning at different intensities; and (3) to elucidate the relationships of soil bacteria and fungi with understory vegetation and soil physicochemical properties.
2. Materials and methods
2.1. Experimental site
The thinning trial Cryptomeria japonica var. sinensis plantation is located at Yangziling Forest Farm, Yaan City, Sichuan Province, China (29°47′37″N, 102°56′18″E). It belongs to the subtropical monsoon mountain climate with abundant rainfall. The annual mean air temperature, precipitation and humidity are 13.1°C, 1,800 mm and 79%, respectively. The plantation was established with a spacing of 2.0 m × 1.5 m in April 2006 with an area of approximately seven hectares. The soil is yellow loam with a pH value of 4.3–4.7. The mean altitude of the site is 1,539 m, and the slope is approximately 15°.
2.2. Experimental design
The thinning trial was carried out in October 2014, and arranged in a randomized complete block design with three replicates and four thinning intensity treatments, including no thinning (control), light-intensity thinning (LIT: 20% of the trees removed), moderate-intensity thinning (MIT: 30% of the trees removed) and strong-intensity thinning (HIT: 40% of the trees removed). Each plot was 600 m2 in size and was surrounded by buffer zones (5 m) to reduce potential edge effects. The stand density and height (m) and diameter at breast height (DBH) (cm) of each tree were measured for each plot just after thinning (October 2014) and October 2019. The tree growth performance of each treatment is presented in Table 1.
Table 1. Information on the thinning trial plantation.
2.3. Understory vegetation investigation and diversity analysis
Four 2 m × 2 m subplots in each plot and one 1 m × 1 m quadrat in each subplot were established to assess the diversity of the shrub and herb layers, respectively. In each subplot or quadrat, the number was counted, and coverage and frequency were measured for each understory plant species in October 2019. The Shannon index and Pielou index were calculated (Wang G. et al., 2021), and the relative species abundance was calculated to indicate the species richness for the shrub and herb layers.
2.4. Soil sampling and physicochemical analyses
In each plot, soil profiles were excavated routinely at five points, and soils were sampled with a cutting ring (100 cm3) at topsoil layers (0–15 cm) to measure soil bulk density (SBD). Two topsoil samples were then collected at each point. One was placed on ice in the field and was then transported promptly to the laboratory and stored at −80°C for DNA extraction, and the other was dried at room temperature for physicochemical analysis. For each purpose, soils from five points were mixed equally as one sample for each plot before further treatment or analysis.
SBD was measured by the cutting ring method (Grossman and Reinsch, 2002), and soil moisture (SM) was determined by oven drying soils at 105°C for 24 h. Soil pH was measured using a pH meter (soil to water ratio was 1:2.5) (Kenworthy et al., 1976). The total nitrogen content (TN) was measured by the Kjeldahl method using a 2,300 Kjeltec Analyzer Unit (FOSS, Sweden) (Bremner and Mulvaney, 1982). The contents of total phosphorus (TP) and total potassium (TK) were determined via the ascorbic acid colorimetric method and atomic absorption method, respectively. The soil organic matter content (SOM) was assessed using the dichromate wet combustion method and a visible spectrophotometer (Lefroy et al., 1993). The alkaline hydrolysis method, molybdenum blue colorimetric method (Tan et al., 2014), and a flame photometer (Lu, 2000) were used to measure the contents of available nitrogen (AN), available phosphorus (AP) and available potassium (AK), respectively.
2.5. DNA extraction, PCR amplification and sequencing
Soil DNA was extracted in triplicate from 0.25 g of each sample using an E.Z.N.A.® soil DNA Kit (Omega Bio-Tek, Norcross, GA, U.S.) following the manufacturer’s instructions. The extraction quality of DNA was detected by 1% agarose gel electrophoresis, and the concentration and purity of DNA were determined by a NanoDrop2000.
The V3-V4 region of the bacterial 16S RNA gene and partial fungal ITS regions were amplified by PCR using the primers 338F (5′- ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′- GGA CTACHVGGGTWTCTAAT-3′), and ITS1F (5’-CTTGGTCATTTA GAGGAAGTAA-3′) and ITS2R (5’-GCTGCGTTCTTCATCG ATGC-3′), respectively (Manter and Vivanco, 2007). All amplifications were performed in 20 μl mixtures containing 4 μl of 5 × FastPfu Buffer, 2 μl of 2.5 mM dNTPs, 0.8 μl of each primer (5 μM), 0.4 μl of FastPfu Polymerase and 10 ng of template DNA. The amplification program included initial denaturation at 95°C for 3 min, followed by 27 cycles at 95°C for 30 s, 55°C for 30 s, and 72°C for 45 s, with a final extension at 72°C for 10 min. However, the PCR cycle number was 35.
The PCR products were purified using an AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA), and were quantified using QuantiFluor™-ST (Promega, USA) (Zhou et al., 2019). Subsequently, all the PCR products were pooled with equal molarity. Finally, the amplicons were pair-ended sequenced on an Illumina MiSeq platform (Illumina, San Diego, USA) by the Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China).
2.6. Sequence splicing and annotation
The raw sequences were quality-filtered and merged using Fastp1 (version 0.19.6) and FLASH2 (version 1.2.11), respectively. Operational taxonomic unit (OTU) clustering and chimera removal were conducted using Upraise software3 (version 11) based on the 97% similarity threshold (Edgar, 2013). The taxonomic assignments of 16S rRNA and ITS sequence reads were determined using the bacterial SILVA reference database (Release1384) and the Unite reference database (Release 8.05), respectively. For both databases, the Ribosomal Database Project (RDP) Classifier6 (version 2.11) was used to perform taxonomic annotation for OTU representative sequences, and the confidence threshold was set to 0.7 to obtain taxonomic annotation results. The sample sequence was flattened according to the minimum number to obtain standardized data for calculation of the Shannon index and Chao index according to Cao et al. (2022), and compositions of the soil microbial community were analyzed at the phylum level.
2.7. Determination of the metabolic and functional prediction
The functional annotation and prediction of metabolic or other putative ecological functions were assessed based on the Tax4Fun-KEGG (Wang et al., 2020) and FAPROTAX databases (Lu et al., 2022) for bacteria. Tax4Fun converts the SILVA-based OTU counts into functional or metabolic profiles and computed metabolic reference profiles based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. The FUNGuild database (Nguyen et al., 2016) was adopted to predict the ecological and biological functions of ITS genes for fungi, which is a flat database hosted by GitHub.7
2.8. Statistical analysis
One-way analysis of variance (ANOVA) and a least significant difference multiple range test (p < 0.05) were performed to assess the significance that thinning affected the assayed soil physicochemical properties (SM, SBD, pH, SOM, TK, TP, TN, AN, AK and AP), understory vegetation diversity (Shannon index, Pielou index, species richness), soil microbial community diversity indices (Shannon index and Chao index), soil microbial community compositions and the relative abundance of the genes associated with different functional categories. Partial least squares discriminant analysis (PLS-DA), a supervised method, was performed on OTU data to discriminate the microbial community profiles of the four thinning treatments (Pérez-Enciso and Tenenhaus, 2003; Ma et al., 2017). Biomarker analysis was performed by linear discriminant analysis (LDA) for effect size (LEfSe) using the Kruskal–Wallis test and All-Against-All (more strict) to determine the significance of differences in soil microbial species among the four thinning treatments. LDA was performed to evaluate the difference in each microbial taxon with a threshold value of 3.5. Redundancy analysis (RDA) was used to test the relationships of microbial groups with understory vegetation and soil properties using Monte Carlo permutations (999 repetitions). Rarefaction curves, PLS-DA, correlation heatmap, and RDA were conducted with R 4.2.0 software.
3. Result
3.1. Understory vegetation and soil properties
The Shannon index, Pielou index and species richness of the shrub layer under the light-(LIT), moderate-(MIT) and high-(HIT) intensity thinning treatments were remarkably higher than those under the control (p < 0.05). While significant difference was absent in the herb layer among the four treatments (p ≥ 0.05) (Table 2).
Table 2. Understory vegetation and soil properties in Cryptomeria japonica var. sinensis plantations under four thinning treatments.
Thinning significantly influenced the soil moisture (SM), soil bulk density (SBD), pH value, and contents of total nitrogen (TN), total potassium (TK), total phosphorus (TP), available potassium (AK) and available phosphorus (AP) (Table 2). The pH value, SBD and TP were significantly higher in the control than in the LIT, MIT and HIT treatments. The SM and TK in the LIT treatment were significantly higher than those in the control, MIT and HIT treatments. The TN was significantly higher in the MIT treatment than that in other treatments. The AP of LIT and HIT was higher than that of the control and MIT treatments. Significant differences were absent in the contents of soil organic matter (SOM) and available nitrogen (AN) among the four thinning treatments.
3.2. Soil microbial community diversity and composition
A total of 318,732 bacterial and 819,288 fungal sequence reads were obtained from the complete dataset of 12 samples and were clustered into 3,994 and 2,807 OTUs, respectively. The rarefaction curve for thinning treatments of Shannon index on OUT level had been well captured at the amount of randomly drawn sequencing data (Figure 1). The Shannon and Chao indices of the bacterial community under the four thinning treatments ranged from 5.677 to 6.101 and from 2,485 to 3,110, respectively, and both indices of the fungal community ranged from 3.287 to 4.571 and from 1,028 to 1,402, respectively (Table 3). All of them were strongly influenced by thinning.
Figure 1. Rarefaction curve of Shannon index of bacteria (A) and fungus (B) communities at OTU level among the four thinning treatments.
Table 3. Diversity indices of the bacterial and fungal communities among the four thinning treatments in C. japonica var. sinensis plantations.
The top three bacterial phyla in the bacterial community were Proteobacteria (relative abundance: 34.93%), Acidobacteria (21.82%), and Actinobacteria (19.57%) (Figure 2A). Rhizobiales (17.21%) and Gammaproteobacteria_Incertae_Sedis (2.1%) were the most abundant orders observed within the Alphaproteobacteria and Gammaproteobacteria classes, respectively (Supplementary Table S1). Significant differences were mostly absent in the relative abundances of Proteobacteria, Acidobacteria and Verrucomicrobiota among the four treatments, while the relative abundance of Actinobacteriota in the LIT, MIT and HIT treatments was significantly higher than that in the control (14.3 ± 3.25%) (Figure 2A).
Figure 2. Relative abundance of the dominant groups of bacterial (A) and fungal (B) communities at the phylum level among the four thinning treatments.
Ascomycota (55.69%), Basidiomycota (14.72%), and Rozellomycota (11.85%) represented most of the observed fungal phyla (Figure 2B). The relative abundances of Ascomycota in the LIT (75.05 ± 3.89%) and MIT (86.46 ± 7.40%) treatments were significantly higher than those in the HIT (59.28 ± 2.75%) and control (20.06 ± 2.62%) treatments (Figure 2B). The relative abundances of Basidiomycota, Rozellomycota and Mortierellomycota in the control were 23.33 ± 2.05%, 24.10 ± 1.51% and 13.76 ± 2.89%, respectively, and were all significantly higher than those in the LIT, MIT and HIT treatments (Figure 2B).
The PLS-DA analysis illustrated that the bacterial communities of the four thinning treatments were clearly separated from each other (Figure 3A), and the same were the fungal communities under the control and HIT treatments, while those under the LIT and MIT treatments were clustered together (Figure 3B).
Figure 3. Partial least squares discriminant analysis (PLS-DA) of bacterial (A) and fungal (B) community composition among four thinning intensity treatment sites.
The LEfSe algorithm was used to determine the taxa in the soil microbial communities (Supplementary Figure S1). In total, 22 bacterial clades presented statistically significant differences (LDA > 3.5, p < 0.05) in all soil samples. There were six differentially abundant taxa in the control (Supplementary Figure S1a), and the most important contribution to the control was the phylum Firmicutes, which accounted for 83%. The LIT treatment had the fewest biomarkers, with only two abundant bacterial clades. For the MIT treatment, the greatest contribution was made by four biomarkers, all of which came from the phylum Firmicutes, including the class Actinobacteria, with LDA scores of 4.50. The HIT treatment contained 17 biomarkers, and the categories with the highest contribution were mainly Acidimicrobia, Planococcaceae, Chitinophagales, and Rhodospirillales.
For fungi, 73 clades were determined by the LEfSe algorithm, which exhibited significant differences under four thinning treatments with an LDA threshold of 3.5 (Supplementary Figure S2). The control accounted for 72.6% of the fungal clades, while the LIT, MIT and HIT treatments had one, 15 and four clades, respectively. The proportion of the phylum Rozellomycota was highest in the control, followed by Basidiomycota (phylum) and Tremellomycetes (class). Specifically, the LIT treatment was only rich in Meripilaceae (family). The clades of the MIT and HIT treatments were mainly concentrated in Ascomycota.
Overall, the soil microorganisms composition was significantly different among the four thinning treatments, which to some extent reflects the impact of thinning on the soil microbial community.
3.3. Relationships of microbial communities with understory vegetation and soil properties
The Pielou index in the herb layer was positively correlated with microbes, especially fungi (Figure 4). The Shannon index and species richness in the shrub layer were significantly negatively correlated with the relative abundances of Zoopagomycota, Rozellomycota, Chytridiomycota, Basidiomycota and Mortierellomycota (Figure 4B). According to the redundancy analysis (RDA), the total interpretation degree of the understory vegetation index to the bacterial genus level was 54.11% (Figure 5A), and three indicators in the shrub layer showed a significant correlation with fungal community composition (Figure 5B).
Figure 4. Correlation heatmaps of the abundance of bacterial (A) and fungal (B) communities at the phylum level with understory vegetation characteristics and soil properties. HH, HP, HS, SH, SP and SS denote the Shannon index, Pielou index and species richness in the herb and shrub layers, respectively. SM, soil moisture; SBD, soil bulk density; SOM, soil organic matter content; TN, TK, TP, AN, AK and AP, contents of total nitrogen, potassium, phosphorus, available nitrogen, potassium and phosphorus, respectively.
Figure 5. RDA of abundant bacterial (A) and fungal (B) communities at the genus level and understory vegetation characteristics for the soil samples from four thinning treatments in Cryptomeria japonica var. sinensis plantations. HH, HP, HS, SH, SP and SS denote the Shannon index, Pielou index and species richness in the herb and shrub layers, respectively. *p < 0.05, **p < 0.01, ***p < 0.001.
SBD and pH value were significantly and positively correlated with the relative abundance of Chloroflexi, while SM and AP were significantly and negatively correlated with it. The relative abundance of Firmicutes was negatively correlated with TN and SOM (Figure 4A). SOM and TN were significantly and positively correlated with the relative abundance of Ascomycota. They were significantly and negatively correlated with the relative abundances of Zoopagomycota, Rozellomycota, Chytridiomycota, Basidiomycota and Mortierellomycota, respectively (Figure 4B).
Figure 6A shows that the first and second ordination axes of RDA explained 61.32 and 16.19% of the variance in the soil bacterial community, respectively, and TP and TN explained the largest proportion of the variance. The pH value, SBD, TN, TK, TP and AP appeared to be the most important environmental parameters for the fungal community, and the variance explained by the first and second ordination axes was 93.73 and 2.68%, respectively (Figure 6B).
Figure 6. RDA of abundant bacterial (A) and fungal (B) communities at the genus level and soil properties for the soil samples from four different thinning treatments in C. japonica var. sinensis plantations. SM, soil moisture; SBD, soil bulk density; SOM, soil organic matter; TN, TK, TP, AN, AK, and AP, contents of total nitrogen, total potassium, total phosphorus available nitrogen, available potassium and available phosphorus, respectively. *p < 0.05, **p < 0.01, ***p < 0.001.
3.4. Microbial function
A total of 40 groups at level 2 of KEGG orthologs (KO) were detected (Supplementary Table S2). There were 6,053 KOs across all samples, and the most abundant functional pathways were presented in Supplementary Figure S3. The abundances of metabolism, environmental information processing and genetic information processing accounted for a relatively high proportion of the four thinning treatments, with mean values of 60.07%, 22.04%, and 10.03%, respectively. Significant differences were found in bacterial functional pathways among the thinning treatments. Analysis of level 2 metabolic pathways showed that the abundances of bacterial functions such as membrane transport, signal transduction, energy metabolism, and cell motility increased significantly with increasing thinning intensity.
The abundances of nucleotide metabolism, translation, replication, and repair were opposite to the abovementioned changes (Supplementary Table S2). The proportion of cellulolysis was significantly increased in the LIT, MIT and HIT treatments. The relative abundances of nitrogen fixation and phototrophy in the LIT and MIT treatments were lower than those in the control (Figure 7A).
Figure 7. Variations in the composition of bacterial (A) and fungal (B) functional groups under different thinning treatments. The bars represent standard deviation of the means (n = 3); different lowercase letters indicate significant differences between different thinning intensities in the same index.
Most fungi could not be classified (classification = ‘Unknown’), and 56.27%, 81.46%, 82.37%, and 70.91% of the relative abundance of fungal functional groups were unclassified in the control, LIT, MIT and HIT treatments, respectively (Supplementary Figure S4). The most abundant fungal parasite, undefined saprotrophs, had a mean abundance of 14.26% in the control, almost three times as much as it did in the LIT or MIT treatment (Supplementary Figure S4). Thinning significantly enhanced the abundance of soil saprotrophs and lichenized-undefined saprotrophs, and the abundance of soil saprotrophs in HIT was 5 times higher than that in the control (Figure 7B). The abundance of endophytes, plant pathogens and ectomycorrhizal fungi decreased with increasing thinning intensity.
4. Discussion
4.1. Response of understory vegetation diversity and soil properties to thinning intensities
In the present study, the Shannon index, the Pielou index and the species richness of the shrub layer in the HIT, MIT and LIT treatments were notably higher than those in the control. The main reason is that thinning results in more or larger canopy gaps and improves understory light conditions, which increases plant diversity (Tamura and Yamane, 2017). However, the effect of thinning on the richness of the herb layer was not as obvious as those in previous reports (Son et al., 2004; Muscolo et al., 2021), and the species richness was greater in the shrub layer than in the herb layer. These results were presumably related to the fact that shrubs are taller than herbaceous plants and can embrace light first, as the amount of light reaching the understory is one of the most important limiting factors affecting the species of understory plants (Yu et al., 2022). Additionally, the lack of seed resources of understory vegetation species in and around the plot due to the relatively high tree density before thinning (Ali et al., 2019b) and the altitude of the site might be another reason for the lower number of vegetation species.
Understory vegetation cover and richness increased with increasing thinning intensities in the present, which led to an increase in the amount of litter (He and Barclay, 2000; Teste et al., 2012), followed by changes in soil physical attributes and nutrient cycling (Elliott et al., 2015; Xiao et al., 2018; Yang et al., 2019). The roots and dead parts of the newly occurred herbs and shrubs are rich in lignin and other substances that enhance the looseness of the topsoil layer (Osman, 2013), which was the main reason why the soil bulk density decreased as the thinning intensity increased in the present study. Then, the humification of litter increased the contents of soil organic acids and organic matter (Sayer, 2006; Zou et al., 2017). The soil pH value difference in the present study was related to the changes in soil organic acid generation activity after thinning, and soil nutrients, especially the contents of soil organic matter (SOM) and total nitrogen (TN) were increased after thinning. Compared with the control, the SOM and TN under the high-intensity thinning treatment (HIT) increased by 19.87% and 35.29%, respectively (Table 2). Thus, it can be seen that LIT and MIT treatments were conducive to improving the understory vegetation diversity and soil physicochemical properties of young C. japonica var. sinensis plantation.
4.2. Responses of soil microbial community diversity and composition to thinning intensities
With increasing thinning intensities, the Shannon and Chao indices of soil bacteria and fungi first decreased and then increased in the C. japonica var. sinensis forest (Table 3). This was different from the study of Dang et al. (2018), in which the Shannon and Chao indices in soil bacteria and fungi did not vary significantly with thinning intensities. Cai et al.’s (2020) study on the thinning of Larix gmelinii var. principis-rupprechtii plantation showed that the Shannon index of bacterial communities were higher in the medium-intensity than in other treatment. Lin et al. (2016) demonstrated that soil fungal communities did not change significantly after 21 months of thinning in a Japanese cedar (Cryptomeria japonica) plantation in Taiwan. The reason for inconsistent results among these studies may be related to the difference in species composition and distribution pattern of the understory vegetation since microbial community composition varies with vegetation species (Overby et al., 2015). Thinning promotes the growth of planted trees and understory vegetation, which is the main reason for changing the composition and decomposition rate of litter and therefore affects the activity of soil microorganisms. This was also demonstrated in a study of mature Pinus contorta forests, which found that the heterogeneity of understory plant species and their rhizosphere resources (e.g., root exudates, nutrients) influenced the patterns of variation in belowground microbial communities (McIntosh et al., 2013). In addition, the changes in the soil bacterial and fungal communities may be related to their adaptability to the changed microenvironment in the stands after thinning.
Thinning significantly changed the soil bacterial community structure in the C. japonica var. sinensis plantation at the phylum level, and the relative abundances of Proteobacteria and Actinobacteriota increased with increasing thinning intensity (Figure 2A). This could be explained by the enhanced soil nitrogen content under thinning treatments (Table 2), as in Wang et al.’s (2018) study on Chinese fir (Cunninghamia lanceolata) plantations, in which nitrogen addition could increase the relative abundance of Proteobacteria and Actinobacteria. The study of Cai et al. (2020) showed that these two phyla had higher abundance in the no thinning treatment of Larix plantations. The most reasonable explanations are that the bacterial composition is closely related to the soil conditions, climate and sampling time (Green et al., 2008), and the adaptability of soil microorganisms in the community to soil environmental changes is different (Navarrete et al., 2015; Ren et al., 2016). Accordingly, the phylum level of bacteria may respond differently to changes in forest density. This study found that soil bacterial communities at the phylum level were mainly composed of Proteobacteria (relative abundance 34.93%) in four treatments (Figure 2A), and the second dominant phylum was Acidobacteriota, followed by Actinobacteriota and Chloroflexi. Coincidentally, the LEfSe algorithm showed that Bacilli (Firmicutes) may be a potential biomarker of bacteria in the control (Supplementary Figure S1). With increasing thinning intensities, Gammaproteobacteria (Proteobacteria) and Acidimicrobia (Actinobacteriota) may be potential biomarkers for light-intensity (LIT) and high-intensity (HIT) thinning treatments, respectively. For moderate-(MIT) intensity thinning treatments, Actinobacteriota and Chloroflexi may be potential biomarkers. Proteobacteria show a significant positive correlation with the soil organic carbon concentration, indicating that increases in root secretions or plant litter decomposition products stimulate their growth (Huhe et al., 2017). Actinobacteria species are recognized as degraders of cellulose, chitin and other complex carbon compounds, which gives them a central role in the carbon cycle and in the turnover of organic matter (Yergeau et al., 2010; Zhou et al., 2017). Acidobacteria and Chloroflexi also play an important role in the decomposition of organic matter and nutrient cycling (Bryant and Frigaard, 2006; Eichorst et al., 2018). Firmicutes can reduce the abundance of plant pathogens (Ali et al., 2019a), and they also embrace sulfate- and iron- reduction abilities (Gupta et al., 2018). These results demonstrate that changes in soil bacterial community composition are directly related to soil characteristics, especially changes in carbon and nitrogen.
The abundance of Ascomycota increased significantly, while Basidiomycota, Rozellomycota, and Mortierellomycota decreased at the soil fungal phylum level with increasing thinning intensities (Figure 2B). Ascomycota is a potential biomarker of the LIT, MIT and HIT treatments, and as the most abundant fungus, Ascomycota accounted for up to 75% (MIT treatment) of fungi in the C. japonica var. sinensis plantation, in which Archaeorhizomycetes had the highest abundance (Supplementary Table S1), and increased with increasing thinning intensities. Archaeorhizomycetes are a type of saprophytic fungi (Rosling et al., 2011) that play an important role in the cycling of carbon derived from living or dead roots in soil (Rosling et al., 2013). The LEfSe algorithm showed that fungi had the most taxonomic groups in the control, which were 53, 3.53 and 13.25 times that of the LIT, MIT and HIT treatments, respectively. Rozellomycota, Basidiomycota, and Tremellomycetes (Basidiomycota) may be potential biomarkers for the control (Supplementary Figure S2). Basidiomycota are found in almost all terrestrial environments, and many of them obtain nutrients by decaying wood, leaf litter and other dead organic matter. Therefore, they play an indispensable role in the carbon cycle. C. japonica var. sinensis had dense branches with a dark and moist environment in the control plots, and it was easy to cause the bottom branches to wither away and fall down. Branches of conifer species are more difficult and need a longer time to be decomposed than shrubs and herbs, and the present study was conducted 5 years after thinning, which may be the reason for the higher abundance of Basidiomycota in the control.
4.3. Relationships of microbial communities with understory vegetation and soil properties
The relationship of microbial community structure with understory vegetation and soil is one of the vital goals of microbial ecology (van der Heijden et al., 2008; Fu et al., 2018). In this study, the Pielou index of the herb layer was positively correlated with the abundance of Chloroflexi, Planctomycetota and Verrucomicrobiota, and was negatively correlated with the abundance of Actinobacteria (Figure 4A). The Shannon index and species richness in the herb layer and the Shannon index, Pielou index and species richness in the shrub layers were almost negatively correlated with the abundance of the top 10 bacterial phyla, and the bacterial community was more likely to be significantly correlated with the soil properties than vegetation. This demonstrated that soil physicochemical properties had more important effects on bacterial community structure than did understory vegetation diversity. We also found that the pH value and soil bulk density (SBD) were significantly positively correlated with the relative abundance of Chloroflexi. One study on the depth profiles of microbial communities in high-elevation soils demonstrated that soil pH was an important driver of forming soil bacterial communities in a given region (Chu et al., 2016). However, Dang et al. (2018) revealed that pH was not correlated with dominant members of bacterial communities in Pinus tabuliformis Carriere plantations. The SM and available phosphorus content (AP) were the main driving factors of the soil bacterial community structure in the C. japonica var. sinensis plantation; they were positively correlated with the abundance of Proteobacteria and negatively correlated with the abundance of Chloroflexi. Proteobacteria can improve soil fertility and sustainability (Rousk et al., 2009; Niu et al., 2020). Chloroflexi is slow-growing and was once classified as an oligotrophic group, and its growth and development may be limited by soil nutrient accumulation (Xu et al., 2021; Zhang et al., 2022), which is why TN, AP and the content of total potassium (TK) were negatively correlated with the abundance of Chloroflexi.
The correlation heatmap of understory vegetation diversity and soil properties with the main phylum of fungi showed that the relative abundances of the main phylum fungi, except Ascomycota and Glomeromycota, were negatively correlated with the indices of the shrub layer and positively correlated with SBD, pH value and the content of total phosphorus (TP). RDA further confirmed that SBD, pH value, TK, TP, TN, AP and shrub diversity affected the abundance of fungal communities. Overall, soil parameters were the most important factors affecting the soil fungal community structure. This was in accordance with Adamo et al.’s (2021) study in 42 pure and mixed pine forests, which showed that soil chemistry significantly affected the variability of soil fungal communities.
Thinning had significant effects on the relative abundances of the dominant fungal communities rather than the dominant bacterial communities, which might indicate that the bacterial community has stronger resistance to changes in stand density than fungal communities (Wang C. Q. et al., 2021). Some studies have also reported that fungal communities have a more obvious response to plant–soil feedback than bacterial communities (Dang et al., 2018; Hou et al., 2021). This is mainly because the bacterial community has a smaller ecological niche in soil and a weaker symbiotic relationship with plants than the fungal community (Shan et al., 2017). Therefore, changes in soil properties, directly or indirectly caused by the intensity of thinning, were less responsive to bacterial than to fungal community diversity.
4.4. Potential metabolic pathways in soils
Understory vegetation diversity and contents of soil nutrients in the C. japonica var. sinensis plantation increased with increasing thinning intensities. Different bacterial and fungal species had different strategies for adaptation to microenvironmental changes in the thinned plantations. Analysis of microbial function prediction showed that 60% of bacteria were involved in metabolic pathways, some of which are known to cause human disease, and the remaining bacteria were involved in genetic information processing, environmental information processing, organismal systems and cellular processes. The abundances of membrane transport, nucleotide metabolism and signal transduction in the LIT, MIT and HIT plots were higher than those in the control. Bacteria often have a strong relationship with human diseases (Geng et al., 2020; Wang et al., 2020). In this study, bacterial infectious diseases accounted for the highest proportion of human diseases, up to 10%. Fortunately, the proportion of human diseases decreased with increasing thinning intensities. This inferred that the change in the bacterial community caused by thinning in the C. japonica var. sinensis plantations might reduce the possibility of harm to humans.
Thinning significantly increased the functional groups ‘soil saprotrophs’ while decreasing the functional groups ‘endophytes’ and ‘plant saprotrophs’ in fungi. This proved that soil nutrients under the C. japonica var. sinensis plantation could be improved by thinning since the abundance of soil saprotrophs was positively correlated with soil fertility (Kyaschenko et al., 2017). This was further confirmed by RDA, from which a significant positive correlation was found between soil properties and the fungal community. Based on the composition and community function of soil bacteria and fungi, it has been predicted that higher thinning intensities are beneficial for soil properties and soil microbial relationships in Cryptomeria japonica var. sinensis plantations.
5. Conclusion
This study showed that thinning enhanced the understory vegetation diversity and improved the soil physicochemical properties. In particular, thinning had a greater effect on the diversity of the shrub layer than the herb layer, and the soil bulk density (SBD) and the contents of soil organic matter (SOM) and total nitrogen (TN) increased with increasing thinning intensities of Cryptomeria japonica var. sinensis plantations. Although the Shannon and Chao indices of soil bacteria and fungi were significantly lower in the HIT, MIT and LIT treatments, the abundance of soil bacterial and fungal species varied significantly with different thinning intensities, the abundance of Actinobacteriota and Ascomycota significantly increased, and the abundance of Basidiomycota, Rozellomycota and Mortierellomycota decreased with increasing thinning intensities. The effects of thinning on microorganisms were mainly driven by soil properties such as pH value and the contents of total nitrogen and total phosphorus, especially for fungi. Fungi are more sensitive to understory vegetation than bacteria. Changes in the distribution of microbial function are a response to changes in the microbial community composition. Thinning improves membrane transport, signal transduction, and cellulolysis in bacteria and soil saprotrophs in fungi. According to the changes of understory vegetation diversity, soil physicochemical properties and microbial composition and function after thinning, LIT and MIT treatments should be adopted in young C. japonica var. sinensis plantations. This study illustrated the relationship of soil microorganisms with the understory vegetation and soil properties in plantations with different thinning densities; however, the effect of thinning on plant–soil-microorganism interactions in plantations needs to be further studied to explore the interaction mechanisms among them.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI—PRJNA928327.
Author contributions
K-LL: investigation, data curation, software, formal analysis, writing—original draft, and writing—review and editing. B-YC: investigation, data curation, and visualization. BZ: conceptualization, methodology, and investigation. R-HW: conceptualization, methodology, supervision, funding acquisition. C-SW: methodology, software, formal analysis, writing—review and editing, and supervision. All authors contributed to the article and approved the submitted version.
Acknowledgments
This work was financially supported by Hunan Provincial Innovation Foundation For Postgraduate (CX20200706) and Scientific Innovation Fund for Postgraduates of Central South University of Forestry and Technology (CX20201005).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1117384/full#supplementary-material
Footnotes
1. ^https://github.com/OpenGene/fastp
2. ^https://ccb.jhu.edu/software/FLASH/index.shtml
3. ^http://www.drive5.com/uparse/
5. ^http://unite.ut.ee/index.php
References
Adamo, I., Castaño, C., Bonet, J. A., Colinas, C., Martínez de Aragón, J., and Alday, J. G. (2021). Soil physico-chemical properties have a greater effect on soil fungi than host species in Mediterranean pure and mixed pine forests. Soil Biol. Biochem. 160:108320. doi: 10.1016/j.soilbio.2021.108320
Ali, A., Dai, D., Akhtar, K., Teng, M., Yan, Z., Urbina-Cardona, N., et al. (2019a). Response of understory vegetation, tree regeneration, and soil quality to manipulated stand density in a Pinus massoniana plantation. Glob. Ecol. Conserv. 20:e00775. doi: 10.1016/j.gecco.2019.e00775
Ali, A., Ghani, M. I., Li, Y. H., Ding, H. Y., Meng, H. W., and Cheng, Z. H. (2019b). Hiseq base molecular characterization of soil microbial community, diversity structure, and predictive functional profiling in continuous cucumber planted soil affected by diverse cropping systems in an intensive greenhouse region of northern China. Int. J. Mol. Sci. 20:2619. doi: 10.3390/ijms20112619
Ares, A., Berryman, S., and Puettmann, K. (2009). Understory vegetation response to thinning disturbance of varying complexity in coniferous stands. Appl. Veg. Sci. 12, 472–487. doi: 10.1111/j.1654-109X.2009.01042.x
Averill, C., Werbin, Z., Atherton, K., Bhatnagar, J., and Dietze, M. (2021). Soil microbiome predictability increases with spatial and taxonomic scale. Nat. Ecol. Evol. 5, 747–756. doi: 10.1038/s41559-021-01445-9
Bona, E., Nadia, M., Omrane, T., Giorgia, N., Patrizia, C., Valeria, T., et al. (2021). Climatic zone and soil properties determine the biodiversity of the soil bacterial communities associated to native plants from desert areas of north-Central Algeria. Microorganisms 9:1359. doi: 10.3390/microorganisms9071359
Bremner, J. M., and Mulvaney, C. S. (1982). “Nitrogen-Total” in Methods of soil analysis: part 2 chemical and microbiological properties. ed. A. L. Page (Madison, WI: American Society of Agronomy, Soil Science Society of America), 595–622.
Cai, M. K., Peng, X. H., Cheng, X. Q., Liu, L., Xing, S. P., Shang, T. X., et al. (2020). Soil element stoichiometry drives bacterial community composition following thinning in a Larix plantation in the subalpine regions of northern China. Forests 11:261. doi: 10.3390/f11030261
Cao, J. J., Jiao, Y. M., Che, R. X., Holden, N. M., Zhang, X. F., Biswas, A., et al. (2022). The effects of grazer exclosure duration on soil microbial communities on the Qinghai-Tibetan plateau. Sci. Total Environ. 839:156238. doi: 10.1016/j.scitotenv.2022.156238
Chen, X. L., Wang, D., Chen, X., Wang, J., Diao, J. J., Zhang, J. Y., et al. (2015). Soil microbial functional diversity and biomass as affected by different thinning intensities in a Chinese fir plantation. Appl. Soil Ecol. 92, 35–44. doi: 10.1016/j.apsoil.2015.01.018
Chu, H. Y., Sun, H. B., Tripathi, B. M., Adams, J. M., Rong, H., Jian, Z. Y., et al. (2016). Bacterial community dissimilarity between the surface and subsurface soils equals horizontal differences over several kilometers in the western Tibetan plateau. Environ. Microbiol. 18, 1523–1533. doi: 10.1111/1462-2920.13236
Dang, P., Gao, Y., Liu, J. L., Yu, S. C., and Zhao, Z. (2018). Effects of thinning intensity on understory vegetation and soil microbial communities of a mature Chinese pine plantation in the loess plateau. Sci. Total Environ. 630, 171–180. doi: 10.1016/j.scitotenv.2018.02.197
Bryant, D. A., and, Frigaard, N. U. (2006). Prokaryotic photosynthesis and phototrophy illuminated. Trends Microbiol. ; 14: 488–496. doi: 10.1016/j.tim.2006.09.001
Edgar, R. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Eichorst, S. A., Trojan, D., Roux, S., Herbold, C., Rattei, T., and Woebken, D. (2018). Genomic insights into the Acidobacteria reveal strategies for their success in terrestrial environments. Environ. Microbiol. 20, 1041–1063. doi: 10.1111/1462-2920.14043
Elliott, K. J., Vose, J. M., Knoepp, J. D., Clinton, B. D., and Kloeppel, B. D. (2015). Functional role of the herbaceous layer in eastern deciduous forest ecosystems. Ecosystems 18, 221–236. doi: 10.1007/s10021-014-9825-x
Fu, G. Y., Yan, B., Ju, X. Q., Mei, Y. X., Ke, Z., Ping, Z. X., et al. (2018). Degradation shaped bacterial and archaeal communities with predictable taxa and their association patterns in Zoige wetland at Tibet plateau. Sci. Rep. 8:3884. doi: 10.1038/s41598-018-21874-0
Geng, S. Y., Cao, W., Yuan, J., Wang, Y. Y., Guo, Y. Q., Ding, A. Z., et al. (2020). Microbial diversity and co-occurrence patterns in deep soils contaminated by polycyclic aromatic hydrocarbons (PAHs). Ecotoxicol. Environ. Saf. 203:110931. doi: 10.1016/j.ecoenv.2020.110931
Goransson, H., Godbold, D. L., Jones, D. L., and Rousk, J. (2013). Bacterial growth and respiration responses upon rewetting dry forest soils: impact of drought-legacy. Soil Biol. Biochem. 57, 477–486. doi: 10.1016/j.soilbio.2012.08.031
Green, J. L., Bohannan, B. J., and Whitaker, R. J. (2008). Microbial biogeography: from taxonomy to traits. Science 320, 1039–1043. doi: 10.1126/science.1153475
Grossman, R. B., and Reinsch, T. G. (2002). 2.1 bulk density and linear extensibility. Methods of soil. Analysis, 201–228. doi: 10.2136/sssabookser5.4.c9
Gupta, A., Dutta, A., Sarkar, J., Panigrahi, M. K., and Sar, P. (2018). Low-abundance members of the Firmicutes facilitate bioremediation of soil impacted by highly acidic mine drainage from the Malakand copper project, India. Front. Microbiol. 9:9. doi: 10.3389/fmicb.2018.02882
He, F., and Barclay, H. (2000). Long-term response of understory plant species to thinning and fertilization in a Douglas-fir plantation on southern Vancouver Island, British Columbia. Can. J. For. Res. 30, 566–572. doi: 10.1139/x99-237
Hou, L. Y., Zhang, Y. Q., Li, Z. C., Shao, G. D., Song, L. G., and Sun, Q. W. (2021). Comparison of soil properties, understory vegetation species diversities and soil microbial diversities between Chinese fir plantation and close-to-natural forest. Forests 12:632. doi: 10.3390/f12050632
Huhe, J. C. X., Jiang, H. F., Fei, W. P., and Xiang, C. Y. (2017). Bacterial and fungal community structures in loess plateau grasslands with different grazing intensities. Front. Microbiol. 8:606. doi: 10.3389/fmicb.2017.00606
Kenworthy, J. B., Allen, S. E., Grimshaw, H. M., Parkinson, J. A., and Quarmby, C. (1976). Chemical analysis of ecological materials. J. Ecol. 64:775. doi: 10.2307/2258784
Kyaschenko, J., Clemmensen, K. E., Karltun, E., and Lindahl, B. D. (2017). Below-ground organic matter accumulation along a boreal forest fertility gradient relates to guild interaction within fungal communities. Ecol. Lett. 20, 1546–1555. doi: 10.1111/ele.12862
Lammel, D. R., Barth, G., Ovaskainen, O., Cruz, L. M., Zanatta, J. A., Ryo, M., et al. (2018). Direct and indirect effects of a pH gradient bring insights into the mechanisms driving prokaryotic community structures. Microbiome 6:106. doi: 10.1186/s40168-018-0482-8
Lefroy, R. D. B., Blair, G. J., and Strong, W. M. (1993). “Changes in soil organic matter with cropping as measured by organic carbon fractions and 13C natural isotope abundance” in Plant nutrition — from genetic engineering to field practice: Proceedings of the twelfth international plant nutrition colloquium, 21–26 september 1993, Perth. ed. N. J. Barrow (Dordrecht: Springer Netherlands), 551–554.
Lei, L., MarcEstiarte, P. B., Jian, L., Asensio, D., Wallander, H., and Peñuelas, J. (2022). Drought legacies on soil respiration and microbial community in a Mediterranean forest soil under different soil moisture and carbon inputs. Geoderma 405:115425. doi: 10.1016/j.geoderma.2021.115425
Lei, L., Xiao, W. F., Zeng, L. X., Frey, B., Huang, Z. L., Zhu, J. H., et al. (2021). Effects of thinning intensity and understory removal on soil microbial community in Pinus massoniana plantations of subtropical China. Appl. Soil Ecol. 167:104055. doi: 10.1016/j.apsoil.2021.104055
Lin, W. R., Wang, P. H., Chen, W. C., Lai, C. M., and Scott, W. R. (2016). Responses of soil fungal populations and communities to the thinning of Cryptomeria Japonica forests. Microbes Environ. 31, 19–26. doi: 10.1264/jsme2.ME15127
Liu, S. Z., Yin, H. F., Li, X. J., Li, X. W., Fan, C., Chen, G., et al. (2021). Short-term thinning influences the rhizosphere fungal community assembly of Pinus massoniana by altering the understory vegetation diversity. Front. Microbiol. 12:620309. doi: 10.3389/fmicb.2021.620309
Lladó, S., López-Mondéjar, R., and Baldrian, P. (2018). Drivers of microbial community structure in forest soils. Appl. Microbiol. Biotechnol. 102, 4331–4338. doi: 10.1007/s00253-018-8950-4
Lu, R.K. (2000). Soil Science Society of China. Analytical methods for soil and agricultural chemistry. China Agricultural Science and Technology Press, Beijing, China, (in Chinese).
Lu, Z. X., Wang, P., Ou, H. B., Wei, S. X., Wu, L. C., Jiang, Y., et al. (2022). Effects of different vegetation restoration on soil nutrients, enzyme activities, and microbial communities in degraded karst landscapes in Southwest China. For. Ecol. Manag. 508:120002. doi: 10.1016/j.foreco.2021.120002
Ma, J., Tang, J. Y., Wang, S., Chen, Z. L., Li, X. D., and Li, Y. H. (2017). Illumina sequencing of bacterial 16S rDNA and 16S rRNA reveals seasonal and species-specific variation in bacterial communities in four moss species. Appl. Microbiol. Biotechnol. 101, 6739–6753. doi: 10.1007/s00253-017-8391-5
Manter, D. K., and Vivanco, J. M. (2007). Use of the ITS primers, ITS1F and ITS4, to characterize fungal abundance and diversity in mixed-template samples by qPCR and length heterogeneity analysis. J. Microbiol Methods 71, 7–14. doi: 10.1016/j.mimet.2007.06.016
McIntosh, A. C. S., Macdonald, S. E., and Quideau, S. A. (2013). Linkages between the forest floor microbial community and resource heterogeneity within mature lodgepole pine forests. Soil Biol. Biochem. 63, 61–72. doi: 10.1016/j.soilbio.2013.03.028
Muscolo, A., Settineri, G., Romeo, F., and Mallamaci, C. (2021). Soil biodiversity as affected by different thinning intensities in a Pinus laricio stand of Calabrian Apennine, South Italy. Forests 12:108. doi: 10.3390/f12010108
Navarrete, A. A., Tsai, S. M., Mendes, L. W., Faust, K., Hollander, M., Cassman, N. A., et al. (2015). Soil microbiome responses to the short-term effects of Amazonian deforestation. Mol. Ecol. 24, 2433–2448. doi: 10.1111/mec.13172
Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006
Ni, Y. Y., Yang, T., Ma, Y. Y., Zhang, K. P., Soltis, P. S., and Soltis, D. E. (2021). Soil pH determines bacterial distribution and assembly processes in natural mountain forests of eastern China. Glob. Ecol. Biogeogr. 30, 2164–2177. doi: 10.1111/geb.13373
Niu, Y., Zhang, M. Y., Hosseini Bai, S., Xu, Z. H., Liu, Y. Q., Chen, F. S., et al. (2020). Successive mineral nitrogen or phosphorus fertilization alone significantly altered bacterial community rather than bacterial biomass in plantation soil. Appl. Microbiol. Biotechnol. 104, 7213–7224. doi: 10.1007/s00253-020-10761-2
Osman, K. T. (2013). “Forest Soils” in Soils: Principles, properties and management. ed. K. T. Osman (Netherlands, Dordrecht: Springer), 229–251.
Overby, S. T., Owen, S. M., Hart, S. C., Neary, D. G., and Johnson, N. C. (2015). Soil microbial community resilience with tree thinning in a 40-year-old experimental ponderosa pine forest. Appl. Soil Ecol. 93, 1–10. doi: 10.1016/j.apsoil.2015.03.012
Pérez-Enciso, M., and Tenenhaus, M. (2003). Prediction of clinical outcome with microarray data: a partial least squares discriminant analysis (PLS-DA) approach. Hum. Genet. 112, 581–592. doi: 10.1007/s00439-003-0921-9
Phillips, M. L., Weber, S. E., Andrews, L. V., Aronson, E. L., Allen, M. F., and Allen, E. B. (2019). Fungal community assembly in soils and roots under plant invasion and nitrogen deposition. Fungal Ecol. 40, 107–117. doi: 10.1016/j.funeco.2019.01.002
Queiroz, M. E. F., Monteiro, J. S., Viana-Junior, A. B., Praxedes, C. L. B., Lavelle, P., and Vasconcelos, S. S. (2021). Litter thickness and soil pH influence the diversity of saprotrophic fungi in primary forest fragments in the Amazon. Pedobiologia 89:150771. doi: 10.1016/j.pedobi.2021.150771
Ren, C. J., Zhao, F. Z., Kang, D., Yang, G. H., Han, X. H., Tong, X. G., et al. (2016). Linkages of C:N:P stoichiometry and bacterial community in soil following afforestation of former farmland. For. Ecol. Manag. 376, 59–66. doi: 10.1016/j.foreco.2016.06.004
Rosling, A., Cox, F., Cruz-Martínez, K., Ihrmark, K., Grelet, G.-A., Lindahl, B., et al. (2011). Archaeorhizomycetes: unearthing an ancient class of ubiquitous soil fungi. Science 333, 876–879. doi: 10.1126/science.1206958
Rosling, A., Timling, I., and Taylor, D. L. (2013). “Archaeorhizomycetes: patterns of distribution and abundance in soil” in Genomics of soil- and plant-associated fungi. 36. eds. B. A. Horwitz, P. K. DMukherjee, M. Mukherjee, and C. P. Kubicek (Berlin Heidelberg, Berlin, Heidelberg: Springer), 333–349.
Rousk, J., Brookes, P., and Bååth, E. (2009). Contrasting soil pH effects on fungal and bacterial growth suggest functional redundancy in carbon mineralization. Appl. Environ. Microbiol. 75, 1589–1596. doi: 10.1128/AEM.02775-08
Sayer, E. (2006). Using experimental manipulation to assess the roles of leaf litter in the functioning of forest ecosystems. Biol. Rev. Camb. Philos. Soc. 81, 1–31. doi: 10.1017/S1464793105006846
Shan, S., Song, L., Avera, B. N., Strahm, B. D., and Badgley, B. D. (2017). Soil bacterial and fungal communities show distinct recovery patterns during forest ecosystem restoration. Appl. Environ. Microbiol. 83:e00966-17. doi: 10.1128/AEM.00966-17
Shang, R. G., Li, S. F., Huang, X. B., Liu, W. D., Lang, X. D., and Su, J. R. (2021). Effects of soil properties and plant diversity on soil microbial community composition and diversity during secondary succession. Forests 12:12. doi: 10.3390/f12060805
Shen, C. C., Xiong, J. B., Zhang, H. Y., Feng, Y. Z., Lin, X. G., Li, X. Y., et al. (2013). Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biol. Biochem. 57, 204–211. doi: 10.1016/j.soilbio.2012.07.013
Son, S. H., Lee, J. K., Hwang, H. S., Eom, T., and Rhim, S. J. (2021). Short-term effects of tree thinning on microhabitat variables and rodents in Japanese larch Larix kaempferi forest. For. Sci. Technol. 17, 39–45. doi: 10.1080/21580103.2021.1880976
Son, Y., Lee, Y. Y., Jun, Y. C., and Kim, Z.-S. (2004). Light availability and understory vegetation four years after thinning in a Larix leptolepis plantation of Central Korea. J. For. Res. 9, 133–139. doi: 10.1007/s10310-003-0071-x
Soong, J., Fuchslueger, L., Marañón-Jiménez, S., Torn, M., Janssens, I., Penuelas, J., et al. (2019). Microbial carbon limitation: the need for integrating microorganisms into our understanding of ecosystem carbon cycling. Glob. Chang. Biol. 26, 1953–1961. doi: 10.1111/gcb.14962
Stark, S., Eskelinen, A., and Mnnist, M. K. (2012). Regulation of microbial community composition and activity by soil nutrient availability, soil pH, and herbivory in the tundra. Ecosystems 15, 18–33. doi: 10.1007/s10021-011-9491-1
Tamura, A., and Yamane, M. (2017). Response of understory vegetation over 10 years after thinning in an old-growth cedar and cypress plantation overgrazed by sika deer in eastern Japan. For. Ecosyst. 4:1. doi: 10.1186/s40663-016-0088-1
Tan, X., Xie, B., Wang, J., He, W., Wang, X., and Wei, G. (2014). County-scale spatial distribution of soil enzyme activities and enzyme activity indices in agricultural land: implications for soil quality assessment. Sci. World J. 2014:535768, 1–11. doi: 10.1155/2014/535768
Teste, F. P., Lieffers, V. J., and Strelkov, S. E. (2012). Ectomycorrhizal community responses to intensive forest management: thinning alters impacts of fertilization. Plant Soil 360, 333–347. doi: 10.1007/s11104-012-1231-6
Tian, J., Dungait, J., Xiankai, L., Yang, Y., Hartley, I., Zhang, W., et al. (2019). Long-term nitrogen addition modifies microbial composition and functions for slow carbon cycling and increased sequestration in tropical forest soil. Glob. Chang. Biol. 25, 3267–3281. doi: 10.1111/gcb.14750
Urbanová, M., Šnajdr, J., and Baldrian, P. (2015). Composition of fungal and bacterial communities in forest litter and soil is largely determined by dominant trees. Soil Biol. Biochem. 84, 53–64. doi: 10.1016/j.soilbio.2015.02.011
van der Heijden, M. G., Bardgett, R. D., and van Straalen, N. M. (2008). The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x
Wang, G., Sun, Y., Zhou, M., Guan, N., Wang, Y., Jiang, R., et al. (2021). Effect of thinning intensity on understory herbaceous diversity and biomass in mixed coniferous and broad-leaved forests of Changbai Mountain. For. Ecosyst. 8:53. doi: 10.1186/s40663-021-00331-x
Wang, Q., Wang, C., Yu, W., Turak, A., Chen, D., Huang, Y., et al. (2018). Effects of nitrogen and phosphorus inputs on soil bacterial abundance, diversity, and community composition in Chinese fir plantations. Front. Microbiol. 9:9. doi: 10.3389/fmicb.2018.01543
Wang, C. Q., Xue, L., Dong, Y. H., and Jiao, R. Z. (2021). Effects of stand density on soil microbial community composition and enzyme activities in subtropical Cunninghamia lanceolate (lamb.) hook plantations. For. Ecol. Manag. 479:118559. doi: 10.1016/j.foreco.2020.118559
Wang, X. J., Zhang, Z. C., Yu, Z. Q., Shen, G. F., Cheng, H. F., and Tao, S. (2020). Composition and diversity of soil microbial communities in the alpine wetland and alpine forest ecosystems on the Tibetan plateau. Sci. Total Environ. 747:141358. doi: 10.1016/j.scitotenv.2020.141358
Wic Baena, C., Andrés-Abellán, M., Lucas-Borja, M. E., Martínez-García, E., García-Morote, F. A., Rubio, E., et al. (2013). Thinning and recovery effects on soil properties in two sites of a Mediterranean forest, in Cuenca Mountain (south-eastern of Spain). For. Ecol. Manag. 308, 223–230. doi: 10.1016/j.foreco.2013.06.065
Xiao, W. Y., Fei, F., Diao, J. J., Chen, B. J. W., and Guan, Q. W. (2018). Thinning intensity affects microbial functional diversity and enzymatic activities associated with litter decomposition in a Chinese fir plantation. J. For. Res. 29, 1337–1350. doi: 10.1007/s11676-017-0536-y
Xu, Y. X., Du, A. P., Wang, Z. C., Zhu, W. K., Li, C., and Wu, L. C. (2020). Effects of different rotation periods of eucalyptus plantations on soil physiochemical properties, enzyme activities, microbial biomass and microbial community structure and diversity. For. Ecol. Manag. 456:117683. doi: 10.1016/j.foreco.2019.117683
Xu, Y. X., Ren, S. Q., Liang, Y. F., Du, A. P., Li, C., Wang, Z. C., et al. (2021). Soil nutrient supply and tree species drive changes in soil microbial communities during the transformation of a multi-generation eucalyptus plantation. Appl. Soil Ecol. 166:103991. doi: 10.1016/j.apsoil.2021.103991
Yang, Y. G., Geng, Y. Q., Zhou, H. J., Zhao, G. L., and Wang, L. (2017). Effects of gaps in the forest canopy on soil microbial communities and enzyme activity in a Chinese pine forest. Pedobiologia 61, 51–60. doi: 10.1016/j.pedobi.2017.03.001
Yang, H., Miao, N., Li, S. C., Ma, R., Liao, Z. Y., Wang, W. P., et al. (2019). Relationship between stand characteristics and soil properties of two typical forest plantations in the mountainous area of Western Sichuan, China. J. Mountain Sci. 16, 1816–1832. doi: 10.1007/s11629-018-5265-y
Yergeau, E., Hogues, H., Whyte, L. G., and Greer, C. W. (2010). The functional potential of high Arctic permafrost revealed by metagenomic sequencing, qPCR and microarray analyses. ISME J. 4, 1206–1214. doi: 10.1038/ismej.2010.41
Yu, J., Zhang, X., Xu, C., Hao, M., Choe, C., and He, H. (2022). Thinning can increase shrub diversity and decrease herb diversity by regulating light and soil environments. Frontiers. Plant Sci. 13:13. doi: 10.3389/fpls.2022.948648
Zhang, C., Bin, L. G., Sha, X., and Liang, W. G. (2016). Soil bacterial community dynamics reflect changes in plant community and soil properties during the secondary succession of abandoned farmland in the loess plateau. Soil Biol. Biochem. 97, 40–49. doi: 10.1016/j.soilbio.2016.02.013
Zhang, Y. T., Yang, J. J., Zhu, L. J., Xue, J. Y., Hu, H. L., Cui, J. B., et al. (2021). Identification of microRNAs and their target genes related to needle discoloration of evergreen tree Chinese cedar (Cryptomeria fortunei) in cold winters. Planta 254:31. doi: 10.1007/s00425-021-03685-2
Zhang, Q. Q., Zhao, W. W., Zhou, Z. Z., Huang, G. H., Wang, X. B., Han, Q., et al. (2022). The application of mixed organic and inorganic fertilizers drives soil nutrient and bacterial community changes in teak plantations. Microorganisms 10:958. doi: 10.3390/microorganisms10050958
Zhou, X. G., Liu, J., and Wu, F. Z. (2017). Soil microbial communities in cucumber monoculture and rotation systems and their feedback effects on cucumber seedling growth. Plant Soil 415, 507–520. doi: 10.1007/s11104-017-3181-5
Zhou, T., Wang, C., and Zhou, Z. (2020). Impacts of forest thinning on soil microbial community structure and extracellular enzyme activities: a global meta-analysis. Soil Biol. Biochem. 149:107915. doi: 10.1016/j.soilbio.2020.107915
Zhou, H., Zhang, D., Jiang, Z., Sun, P., Xiao, H., Yuxin, W., et al. (2019). Changes in the soil microbial communities of alpine steppe at Qinghai-Tibetan plateau under different degradation levels. Sci. Total Environ. 651, 2281–2291. doi: 10.1016/j.scitotenv.2018.09.336
Keywords: soil microbial community composition, microbial function, thinning, soil properties, understory vegetation
Citation: Liu K-L, Chen B-Y, Zhang B, Wang R-H and Wang C-S (2023) Understory vegetation diversity, soil properties and microbial community response to different thinning intensities in Cryptomeria japonica var. sinensis plantations. Front. Microbiol. 14:1117384. doi: 10.3389/fmicb.2023.1117384
Edited by:
Luis Raul Comolli, Independent Researcher, Basel, SwitzerlandReviewed by:
Yang Xiao-yan, Dali University, ChinaEverlon Cid Rigobelo, São Paulo State University, Brazil
Copyright © 2023 Liu, Chen, Zhang, Wang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rui-Hui Wang, wang626389@163.com; Chun-Sheng Wang, wangcs87@caf.ac.cn