Next Article in Journal
Approaches to Forest Site Classification as an Indicator of Teak Volume Production
Next Article in Special Issue
Seasonal Variation and Soil Texture-Related Thinning Effects on Soil Microbial and Enzymatic Properties in a Semi-Arid Pine Forest
Previous Article in Journal
Opportunities and Prospects for the Implementation of Reforestation Climate Projects in the Forest Steppe: An Economic Assessment
Previous Article in Special Issue
Arthropod Recolonization of Soil Surface Habitat in Post-Fire Mulch Treatments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Early Effect of Plant Density on Soil Physicochemical Attributes and Bacterial and Understory Plant Diversity in Phoebe zhennan Plantations

1
Forestry Ecological Engineering in the Upper Reaches of the Yangtze River Key Laboratory of Sichuan Province & National Forestry and Grassland Administration Key Laboratory of Forest Resources Conservation and Ecological Safety on the Upper Reaches of the Yangtze River & Rainy Area of West China Plantation Ecosystem Permanent Scientific Research Base, College of Forestry, Sichuan Agricultural University, Chengdu 611130, China
2
Sichuan Key Laboratory of Ecological Restoration and Conservation for Forest and Wetland, Sichuan Academy of Forestry, Chengdu 610081, China
3
Sichuan Academy of Grassland Sciences, Chengdu 611731, China
*
Authors to whom correspondence should be addressed.
Forests 2023, 14(8), 1612; https://doi.org/10.3390/f14081612
Submission received: 8 May 2023 / Revised: 21 July 2023 / Accepted: 6 August 2023 / Published: 10 August 2023

Abstract

:
The effect of stand density on the soil bacterial community and diversity remains unclear. Spectrophotometry and full-length 16S rRNA sequences were used to determine the effects of planting density on soil physicochemical attributes and the diversity of soil bacterial and understory vegetation in a young Phoebe zhennan plantation at five densities. The findings showed that stand density had significant effects on the total nitrogen, ammonium nitrogen ( NH 4 + - N ), nitrate-nitrogen ( NO 3 - N ), organic carbon, and the dominance and evenness of shrubs. Candidatus Udaeobacter and Candidatus Soilbacter were the two most common genera across the five stand densities. The density D5 (850 stems/hm2) demarcated from the others with a lower diversity of soil bacteria. Overall, the relatively low- and middle-density plantations were more conducive to complex and stable understory vegetation, bacterial communities, and soil nutrient cycles. The functional categories of the bacterial communities revealed a high proportion associated with chemoheterotrophy, aerobic chemoheterotrophy, and nitrogen fixation. Bacterial diversity and function were significantly influenced by soil pH, NH 4 + - N , NO 3 - N , total phosphorus, and available phosphorus. However, there were no significant correlations between soil physicochemical attributes, understory vegetation, and bacterial diversity. Therefore, we speculated that the key drivers of the soil bacterial community were the soil physicochemical attributes and that stand density affected the soil bacterial community diversity by changing the soil physicochemical attributes. Overall, P. zhennan plantations with densities below 600 stems/hm2 were conducive to complex and stable soil bacterial communities and nutrient cycles.

1. Introduction

In plantation ecosystems, stand density is an essential forest management strategy that influences soil microbial diversity, understory vegetation, and soil characteristics [1,2,3]. The diversity and composition of soil microbial communities vary considerably across temporal and geographical scales in terrestrial ecosystems, and they drive most biogeochemical processes and regulate soil nitrogen cycling [4,5,6]. Previous studies have supposed that soil physicochemical attributes, such as pH, soil organic, organic carbon, and soil nitrogen, are closely related to the communities of soil microbes [4,7,8,9,10,11,12,13]. Research on the driving factors and functions of soil microbes has been the focus of forest management strategies [3]. Some studies have shown that the nutrient cycle, understory vegetation diversity, and microbial community are also influenced by stand density [1,14,15,16,17]. In Chinese fir plantations, high-density plantations decreased the composition and diversity of microbes and the contents of soil pH, N, and P [2,3]. Stand density also affects the metabolic function of soil microbes in plantation systems [15]. Similar results were also reported in the forests of Catalonia, which showed that an increase in density decreased soil bacterial diversity by altering the understory vegetation and litter input [18]. In summary, stand density can influence the community and diversity of soil microbes. Nevertheless, the response of soil physicochemical attributes, understory vegetation, and soil microbes of plantations to density regulation is different with different tree species and tree ages, and the difference in optimal planting density is significant. For instance, Chinese fir plantations (37 years old) with a stand density below 3333 stems/ha are essential for maintaining soil fertility and improving forest ecosystem productivity [3]. Reducing the density (1212 stems/ha) in Pinus massoniana plantations (62 years old) can promote the diversity of understory vegetation and soil fertility [1]. However, 4300 stems/ha is the suitable stand density for 8-year-old Dendrocalamus minor var. amoenus plantations, which is more conducive to soil quality [14]. Therefore, Chen et al. [19] noted a fundamental change in the biomass of understory vegetation after 17 years of density-reduction treatment in Picea spp. forests. As such, a deeper understanding is required to clarify how stand density affects soil physicochemical attributes, understory vegetation, and the composition and diversity of soil microbes. The results could provide fundamental information for plantation management.
Phoebe zhennan S. Lee et F. N. Wei is a precious tree species in China, and its wood is very valuable and famous with a golden color, which is called the Golden-thread nanmu [20,21,22]. The price of P. zhennan wood in China is higher than 5000 RMB/m3, and the price can reach 10,000 RMB/m3 if the diameter at breast height (DBH) of the tree exceeds 30 cm [23]. Consequently, the area of P. zhennan plantations has increased rapidly due to the high value of the wood. However, unscientific density management has induced low production and poor value in the process of P. zhennan plantation management. Previous studies on plantation ecosystems regarding stand density have mainly focused on fast-growing and high-yielding coniferous plantations [2,3,15,24,25]. Differences in litter between coniferous and broadleaf forests result in different responses to changes in stand density [26]. Given the importance of stand density in soil microbial composition, a better understanding of the soil nutrient cycle would be extremely valuable in P. zhennan plantations. However, no comprehensive analysis of the effects of stand density on soil physicochemical attributes, understory vegetation, and soil microbes or any detailed observations of the correlations between biological and physicochemical properties have been performed in P. zhennan plantations. Therefore, we aimed to provide a scientific reference for density management in P. zhennan plantations to enhance soil-fertility maintenance techniques. Consequently, our work had two objectives: (1) to determine how soil bacterial compositions, soil physicochemical properties, and understory vegetation diversity vary with stand density and (2) to identify the main factors influenced by plant density that drive the soil bacterial community composition and diversity.

2. Materials and Methods

2.1. Study Site

The experimental site was in Pengshan County, Sichuan Province, China (103.969021° E, 30.304569° N), with a mean elevation of 453 m and a mean slope of 10°. The mean annual temperature, mean annual precipitation, and mean annual relative humidity of the experimental site were 18.3 °C, 1110.3 mm, and 76%, respectively.

2.2. Experimental Plot Design and Soil Sampling

P. zhennan plantations were planted in 2015 using 2-year-old seedlings and a randomized block design. Five planting densities were designed, 400 (D1), 500 (D2), 600 (D3), 700 (D4), and 850 stems/hm2 (D5), with similar site conditions of topography, altitude, soil type, and slope position. A plantation of Eucalyptus grandis was built before the construction of the P. zhennan plantation in this study site, and the P. zhennan plantation was constructed after felling. Routine tending was carried out three years before afforestation, and there was no artificial interference in the later stage. Three 20 m × 20 m plots per density (three repetitions per density) were selected for understory vegetation investigation and soil sampling (Table 1). The height (H, m), diameter at 1.3 m height (DBH, cm), and tree volume (V, m3) were measured in November 2021. The tree volume was calculated using the equation V = 0.000074954DBH1.884652825H0.881513308 [27]. In April 2021, five soil samples (0–20 cm) per plot were collected using the “five points sampling” method with a 5 cm diameter earth screw, and then the five soil samples were mixed together to make one composite sample [2,3,28]. The soil samples were divided into two parts: one part was used for 16S rRNA sequencing after passing through a 2 mm sieve to remove plant tissues and rocks, and the other part was used to analyze soil chemical properties after air-drying, grinding, and passing through a 0.15 mm sieve.

2.3. Understory Plant Community Inventory

Fifteen subplots (1 m × 1 m) per density were designed to investigate the herbs. For shrubs, the whole plot was investigated. The shrubs and herbs were identified, and their height, quantities, and coverage were measured. All individual species and species coverage in the shrub and herb plots were identified by three observers working together. The importance value (IV) and Shannon–Wiener (H′), Pielou (Jsw), Simpson (D), and Margalef (DMG) indices, which reflect species richness and diversity, were estimated through the following formulas [29,30]:
IV = (RH + RC + RF)/3
H′ = −∑(Pi)ln(Pi)
D = 1 − ∑(Pi)2
Jsw = −∑PilnPi/lnS
DMG = (S − 1)/lnN
where RH, RC, RF, N, S, and ln indicate the relative height, relative coverage, relative frequency, number of individuals, number of herbs/shrubs, and log base-e, respectively. Pi is the proportion of individuals or the abundance of the ith herb/shrub expressed as the proportion of the total in the community.

2.4. Soil Physicochemical Analyses

Soil pH, soil moisture (SM), soil bulk density (SBD), organic carbon (SOC), total nitrogen (TN), ammonium nitrogen ( NH 4 + - N ), nitrate-nitrogen ( NO 3 - N ), total phosphorus (TP), and available phosphorus (AP) were tested through standard soil physical and chemical analyses [31]. A pH meter was used to measure the soil pH. The ring-knife soil samples were dried in a constant temperature oven to determine the SM and calculate the SBD. The Walkley–Black K2Cr2-H2SO4 wet oxidation method was used to measure SOC. TN was measured using the Kjeldahl method. The phenol sulfonic acid colorimetric process and the KCl extraction–indophenol blue colorimetric technique were utilized to quantify NO 3 - N and NH 4 + - N . TP was measured using HClO4-H2SO4 colorimetry, and AP was determined using acid solution–molybdenum antimony resistance.

2.5. Amplification of Full-Length 16S rRNA Genes

Despite previous studies that used specific regions of the 16S rRNA gene to assign taxonomic assignments showing that taxonomic assignments are highly sensitive to specific regions, high-throughput analysis revealed unintended missing classifications with low accuracy, especially at the genus and species levels [32]. Therefore, more precise identification methods are required to determine the complexities of soil microbial communities. Recently, full-length 16S rRNA sequencing has greatly contributed to obtaining better resolution in soil microbial taxonomy [33]. Hence, the soil bacterial communities and diversity in P. zhennan plantations were analyzed by full-length 16S rRNA sequencing. The soil genomic DNA was extracted using a Soil Genomic DNA Kit (TIANGEN, Beijing, China), and then the full-length 16S rRNA gene was amplified. The PCR for each sample was run in triplicate and performed in a total volume of 30 μL containing 1.5 μL DNA, 10.5 μL H2O, 15.0 μL KOD OneTM PCR Master Mix (Beijing Biolink Biotechnology Co., Ltd., Beijing, China), and 3.0 μL primers (27F: 5′-AGRGTTTGATYNTGGCTCAG-3′ and 1492R: 5′-TASGGHTACCTTGTTASGACTT-3′). The PCR protocol was as follows: 95 °C for 2 min, 25 cycles at 98 °C for 10 s, 55 °C for 30 s, 72 °C for 90 s, and 72 °C for 2 min. The contaminants and undesired PCR products were removed through gel-based size selection with 1.8% gel electrophoresis. DNA damage repair, end repair, and ligation of hairpin adapters were performed according to the manufacturer’s instructions. The binding concentration and the on-plate loading concentration of DNA template libraries were optimized using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Then, Sequel Binding Kit 2.0 (Pacific Biosciences, Menlo Park, CA, USA) was used to bind the libraries to sequel polymerase 2.0. The data collection per sample was performed in a single Sequel SMRT Cell 1 M v2 with a 30 h movie time using Sequel Sequencing Kit 2.0 (Pacific Biosciences, Menlo Park, CA, USA).

2.6. Sequencing Analysis

Original subreads were corrected to generate circular consensus sequencing (CCS) reads using SMRT Link (version 8.0). Then, the CCS reads were screened using Mothur v.1.39.5 [34], and sequences outside the range <1200 bp to >1650 bp were removed using cutadapt v2.7 [35]. Chimeras were removed using UCHIME v4.2 [34,36]. Operational taxonomic unit (OTU) classification was performed with a 97% identity cutoff using USEARCH version 10.0, and the conservative threshold for OTU filtration was 0.005% [37]. Taxonomic alignment was performed using the SILVA 132 reference with a 90% cutoff value [38]. Alpha diversity analysis included the calculation of various diversity indices (Simpson and Shannon diversity index), richness abundance estimators (Chao1 and ACE), and Good’s coverage, calculated using QIIME [39]. Principal coordinate analysis (PCoA) was used to evaluate beta diversity with weighted and unweighted UniFrac distance metrics [33]. The functional profiles of the microbial communities were predicted by FAPROTAX [40].

2.7. Statistical Analysis

A general linear model (GLM) test was used to assess the effects of stand density on soil physicochemical attributes, understory vegetation diversity, and the diversity and relative abundance of soil bacteria. Then, an LSD test was used to calculate significant differences (α = 0.05). The abundant taxa used to characterize the microbiomes associated with different stand densities were identified by linear discriminant analysis effect size (LEfSe) [41]. Pearson’s correlation coefficient was calculated to assess the correlation of the relative abundance of the bacteria at the genus level (p = 0.05), and the soil physicochemical attributes and understory vegetation diversity and microbial diversity indices were also determined (p = 0.05). Spearman rank correlation analysis was performed on the difference in all species abundances among samples to identify data with a correlation coefficient larger than 0.1 and p value smaller than 0.05, which were used to construct the network. Redundancy analysis (RDA) was conducted using a Bray–Curtis distance matrix plot calculated by the vegan package in R to understand the degree to which the soil physicochemical attributes elucidated bacterial and understory vegetation diversity [42].

3. Results

3.1. Understory Community Structure

The understory vegetation of the P. zhennan plantation was composed of 63 genera of herbs and shrubs (Tables S1 and S2). There were no significant differences in herb community diversity among the five stand densities (Table 2). Within the herb layer, the Simpson value (D) in the D5 stand was obviously higher than that in the D2, D3, and D4 stands. The H′, Jsw, and DMG values first increased and then decreased with increasing stand density. There were significant differences in the Jsw and D values at the shrub level. The Jsw value in the D1, D4, and D5 stands was significantly higher than that in the D3 stand. In contrast, the D value in the D5 stand was significantly higher than that in the D1, D4, and D5 stands. The DMG and H′ values of shrubs tended to decrease first and then increase with increasing stand density.

3.2. Soil Physicochemical Attributes

There were no significant differences in soil pH, soil bulk density (SBD), soil moisture (SM), total phosphorus (TP), or available phosphorus (AP) among the five stand densities (Table 3). The soil of the P. zhennan plantation is acidic soil, and obvious differences were found among the five stand densities (pH = 5.08–5.85). Significant differences were found in total nitrogen (TN), NH 4 + - N , NO 3 - N , and organic carbon (SOC) among the five density plantations. The concentrations of SOC in the D1, D2, and D3 stands were significantly higher than those in the D4 stands. The concentrations of TN in the D1 and D3 stands were significantly higher than that in the other stands. However, the concentration of NH 4 + - N in the D5 stand was significantly higher than those in the other stands, and NO 3 - N in the D5 and D3 stands was significantly higher than that in the other stands.

3.3. Overview of Full-Length 16S rRNA Amplification

A total of 174,192 effective CCS reads were separated into 15,849 OTUs across all samples, and the mean number of reads per sample was 9941 sequences (range 9108–10,439 sequences) (Table S3). Good’s coverage values were higher than 97%, revealing that most phylogenetic groups in all samples were captured (Table 4). No significant differences were estimated in the richness estimator among the five stand densities. However, the ACE, Chao1, and Shannon values first increased gradually (from D1 to D4) and then decreased (from D4 to D5). The ACE, Chao1, and Shannon values in the D5 stand were obviously lower than those in the D2 and D4 stands (Table 4). Additionally, the PCoA results also indicated that the distribution of the soil bacterial communities was random (Figure 1).

3.4. Soil Bacterial Community Composition

In total, 33 phyla, 154 orders, 359 genera, and 434 species were identified in the P. zhennan plantation. At the phylum level, the dominant bacterial communities were Proteobacteria, Acidobacteria, Verrucomicrobia, and Bacteroidetes, which comprised more than 80% of the sequences under the five stand densities (Figure 2A). The relative abundances of Acidobacteria in D5 and Verrucomicrobia and Bacteroidetes in D2 were higher than those in the other stands (Figure 2B). Candidatus Udaeobacter, Candidatus Soilbacter, and ADurb. Bin063-1 were the predominant genera in all samples (Figure 2C). The relative abundances of Occallatibacter and Bryobacter in D3 and Candidatus Udaeobacter in D5, D3, and D2 were higher than those in the other stands (Figure 2D). The results of the correlation analysis showed that the relative abundance of Bryobacter was significantly positively correlated with Occallatibacter and Bradyrhizobium (p ˂ 0.05); nevertheless, there was no significant correlation between Occallatibacter and Bradyrhizobium (r = 0.26, p ˃ 0.05) (Table S4, Figure 3). The dominant species in the P. zhennan plantation were Bradyrhizobium jicamae, Sphingomonas parvus, and Paraburkholderia terrae (Figure S1). Overall, there were no significant differences in the relative abundance of the bacterial community between the five density stands at every classification level (Figure 3). At the genus level, the relative abundance of the top 20 bacterial communities also showed no significant differences among the five density stands (Table 5).

3.5. Predictive Functional Profiles of Bacterial Communities

Among the putative functions predicted by FAPROTAX, the highest relative abundance of the putative function was related to chemoheterotrophy (28.40%–31.08%) (Figure 4A). The relative abundances of bacteria with aerobic chemoheterotrophic and nitrogen fixation functions were estimated to be 26.93% and 6.87%, respectively. The stand density could affect the putative functions of bacterial communities with putative functions, although the difference in bacteria with putative functions was not significant among the five densities (Figure S2). Overall, the D4 stand had the highest relative abundance of bacteria with nitrification and aerobic ammonia oxidation functions. The relative abundance of bacteria that functioned as animal parasites or symbionts in the D5 stand was higher than that in the other density plantations. The results of the correlation analysis showed that the relative abundance of bacteria with nitrification functions and predatory or exoparasitic bacteria was significantly negatively correlated with NO 3 - N and AP but significantly positively correlated with the pH value (Figure 4B).

3.6. Relationship between Soil Physicochemical Attributes and Bacterial and Understory Diversity

The redundancy analysis (RDA) results revealed that there were close relationships between the soil physicochemical attributes and the compositions of soil bacteria, herbs, and shrubs under density control (Figure 5). The RDA explained 35.34%, 33.83%, and 35.44% of the variation in the bacterial composition, herbs, and shrubs, respectively. From the perspective of the relationship between the understory vegetation and soil physicochemical attributes, the distribution frequency of most herbs and shrubs was positively related to the concentrations of NO 3 - N and NH 4 + - N (Figure 5B,C). The bacteria were concentrated in the third and fourth quadrants at the genus level, and the NH 4 + - N , AP, SOC, and TP contents had positive correlations with the genera Bryobacter, Candidatus Solibacter, and Bradyrhizobium. The alpha diversity of bacteria was significantly influenced by pH, NH 4 + - N , NO 3 - N , TP, and SOC (Figure 6A). However, there were no significant correlations between the soil physicochemical attributes and understory vegetation (Figure 6B,C). The correlation between the understory vegetation and bacterial diversity was also not significant (Figure 6D,E).

4. Discussion

In multiple studies, stand density was one of the key factors that affected soil physicochemical attributes in plantations [1,16]. In Phoebe zhennan plantations, the influence of stand density on total nitrogen (TN), ammonium nitrogen ( NH 4 + - N ), nitrate-nitrogen ( NO 3 - N ), and organic carbon (SOC) contents was significant. Plants at higher densities need to use and uptake greater amounts of nutrients [43]. Therefore, the contents of TN and SOC in the relatively high-density plantations were lower than those in the relatively low-density plantations of P. zhennan. We also found that the soil pH decreased with increasing density, indicating that soil acidification increased with increasing density, which was caused by organic acids produced by roots in relatively high-density plantations [44]. In acidic soils, phosphate easily precipitates with iron, aluminum, and manganese in the soil, reducing the available phosphorus [45]. Consistent with the results of a previous study, our current study found that the available phosphorus (AP) content reached a maximum with the lowest pH at the highest density.
Stand density and soil physicochemical attributes explained most of the variation in the diversity and composition of the bacterial communities in plantations [4]. Every plant tends to be associated with a specific soil microbial community. In P. zhennan plantations, the dominance of soil bacteria at every taxonomic level is different from that in other tree species plantations [2,3,4]. This may be caused by species specificity and the planting environment. For instance, Verrucomicrobia and Bacteroidetes were the dominant bacteria at the phylum level, but these two bacteria were not dominant in Chinese fir plantations [3]. At the genus level, Subgroup-2-norank, JG37-AG-4-norank, and Acidobacteriaceae-Subgroup-1-uncultured were the dominant bacteria in Hevea brasiliensis plantations [46], but these were not determined in our study. Plantations with higher stand density changed plant detrital inputs and subsequently enhanced the activity of soil microbial communities [4]. There was no significant difference in the relative abundance of bacterial communities among the five stand densities, indicating that the soil bacterial communities are resilient and not sensitive to density control in young P. zhennan plantations. Multiple studies have suggested that density is one of the most important driving factors of soil microbial composition [2,3,4]. Consistent with those of previous studies, our results showed that the diversity of the bacterial community in the relatively low- and middle-density plantations (400–700 stems/ha) was higher than that in the high-density plantations, indicating that P. zhennan plantations with low densities are conducive to more diverse bacterial communities. This may be partly explained by changes in canopy, litter species, and understory vegetation [47,48].
A previous study considered that density control significantly affected the distribution and relative abundance of soil microbes [3,4]. In this study, the diversity of the bacterial community was significantly related to soil pH, NH 4 + - N , NO 3 - N , TP, and SOC, suggesting that various soil physicochemical attributes can affect the diversity of bacterial communities in P. zhennan plantations. Previous reports have suggested that soil pH is a predictor of soil microbial communities [4,49]. We also observed a significant positive correlation between soil pH and the diversity of the bacterial community. The relative abundance of the bacterial community with putative functions was also significantly positively correlated with the soil pH, such as nitrification and aerobic ammonia oxidation. In terms of nutrients, higher soil N and P contents indicate a higher efficiency of soil microbial decomposition. In contrast, TP, NO 3 - N , and NH 4 + - N showed a negative correlation with the diversity of bacteria, indicating that N and P were limited during the bacterial decomposition process in this study. Soil bacteria and understory vegetation compete for available N and P [16]. However, there was no significant correlation between the diversity of bacteria and understory vegetation, suggesting that the effects of the diversity of understory vegetation on soil bacterial community diversity in the early stage of P. zhennan plantations were not obvious. Compared with understory vegetation, litter decomposition and nutrient return played a key role in topsoil physiochemical properties and bacterial diversity, at least in the early stage of P. zhennan plantations. The results of the correlation analysis also showed that there were no significant correlations between the soil physicochemical attributes and the diversity of the understory vegetation. Consequently, bold speculation suggests that soil physicochemical attributes are the main drivers of soil bacterial community composition and diversity in P. zhennan plantation ecosystems.
Microorganisms are important factors that influence nutrient availability in plant–soil feedback systems [50,51]. We found that density management significantly influenced the composition and diversity of soil microbial communities along with soil physicochemical attributes, which is consistent with the results from previous studies [2,3,4,15]. We also found that a handful of key bacterial species were affected by density-dependent P. zhennan plantation soil feedback, although soil bacterial communities had no significant response to changes in stand density. For instance, the abundances of Acidobacteria and Acidobacteriales in the high-density stands were higher than those in the low-density stands, and they were also the dominant bacteria at the phylum and order levels. The activity of Acidobacteria affected soil acidity and nutrient availability and had a negative response to increases in carbon and pH [52,53,54], suggesting the presence of higher soil acidity and lower carbon availability in relatively high-density stands than in relatively low-density stands in P. zhennan plantations. In this study, changes in stand density with changes in soil pH and SOC may have resulted in differences in the relative abundance of Acidobacteria in plantations with different stand densities. Wang et al. [2] suggested that low-density plantations may be more beneficial to soil nutrient cycles. In our study, the relative abundance of bacteria with nitrification and nitrogen fixation functions at low and moderate densities was higher than that at high densities, which indicates that plantations with an intermediate density had a preference for the soil N cycle. The relative abundance of bacteria with putative functions was influenced differently by soil chemical properties [55]. For example, the relative abundance of bacteria with a nitrogen fixation function was significantly positively correlated with TN, but there was no difference among the stand densities. The relative abundance of bacteria involved in nitrification was significantly positively correlated with soil pH, indicating that the management of soil pH is important for soil fertility and the N cycle. Consistent with the literature [15,16,55], these results implied that stand density can alter soil processes via functional bacteria and that soil physicochemical attributes are important for the functional engineering of soil bacteria, at least in early-stage P. zhennan plantations.

5. Conclusions

The main goal of the current study was to determine the effects of shifts in stand density on soil physicochemical attributes, understory vegetation, and the bacterial community across young P. zhennan plantations under different stand densities. Low- and middle-density plantations showed higher contents of soil TN, TP, and SOC and a higher diversity in herb and bacterial communities than high-density plantations, suggesting that the soil ecosystem functions of low- and middle-density plantations are more stable than those of high-density plantations. Therefore, to maintain good soil ecosystem function, we suggest that the planting density of P. zhennan should be controlled below 600 stems/ha. We also found that the effects of pH, TN, NH 4 + - N , NO 3 - N , and AP on herb diversity and the effects of NH 4 + - N on shrub diversity were significant. The results of this study provide essential insights that soil physicochemical attributes are indeed factors that drive the formation of understory vegetation and bacterial communities in P. zhennan plantations. This is also an important factor for how stand density influences the diversity of understory vegetation and the bacterial community.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/f14081612/s1, Figure S1: Bacterial community patterns at the species level at different stand densities; Figure S2: Average relative abundance of bacteria contributing to the top 10 relative abundances in different putative functions; Table S1: Herb species importance values (IVs) in understory woody species of Phoebe zhennan plantations with different densities; Table S2: Shrub species importance values (IVs) in understory woody species of Phoebe zhennan plantations with different densities; Table S3: Statistics of bacterial high-quality sequences and OTUs; Table S4: Correlation of dominant genera in bacteria.

Author Contributions

Conceptualization, H.Y. and Y.G.; methodology, Y.C., H.Y. and T.J.; formal analysis, Y.C., J.P. and H.G.; writing—original draft preparation, H.Y., J.P. and Y.C.; writing—review and editing, H.Y.; funding acquisition, H.Y. and Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Research Funds of Key Research and Development Project of Sichuan Province (2021YFYZ0032), Natural Science Foundation of Sichuan Province (2022NSFSC1062), Forest Ecosystem Improvement in the Upper Reaches for Forest and Wetland (2022SKT4-02), and Forest Ecosystem Improvement in the Upper Reaches of Yangtze River Basin Program (510201202038467).

Data Availability Statement

All data analyzed during this study are included in this article and its additional files.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ali, A.; Dai, D.; Akhtar, K.; Teng, M.; Yan, Z.; Urbina-Cardona, N.; Mullerova, J.; Zhou, Z. Response of understory vegetation, tree regeneration, and soil quality to manipulated stand density in a Pinus massoniana plantation. Glob. Ecol. Conserv. 2019, 20, e00775. [Google Scholar] [CrossRef]
  2. Wang, C.; Xue, L.; Dong, Y.; Hou, L.; Wei, Y.; Jiao, R. Responses of soil microbial community structure to stand densities of Chinese fir plantations. J. For. Res. 2019, 24, 162–167. [Google Scholar] [CrossRef]
  3. Wang, C.; Xue, L.; Dong, Y.; Jiao, R. Effects of stand density on soil microbial community composition and enzyme activities in subtropical Cunninghamia lanceolate (Lamb.) Hook plantations. For. Ecol. Manag. 2021, 479, 118559. [Google Scholar] [CrossRef]
  4. Li, S.; Huang, X.; Lang, X.; Xu, F.; Li, H.; Zheng, M.; Su, J. Effect of selective logging on soil microbial communities in a Pinus yunnanensis forest. Land Degrad. Dev. 2020, 31, 2268–2280. [Google Scholar] [CrossRef]
  5. Cui, Y.; Wang, X.; Zhang, X.; Ju, W.; Duan, C.; Guo, X.; Wang, Y.; Fang, L. Soil moisture mediates microbial carbon and phosphorus metabolism during vegetation succession in a semiarid region. Soil Biol. Biochem. 2020, 147, 107814. [Google Scholar] [CrossRef]
  6. Heydari, M.; Eslaminejad, P.; Kakhki, F.V.; Mirab-balou, M.; Omidipour, R.; Prévosto, B.; Kooch, Y.; Lucas-Borja, M.E. Soil quality and mesofauna diversity relationship are modulated by woody species and seasonality in semiarid oak forest. For. Ecol. Manag. 2020, 473, 118332. [Google Scholar] [CrossRef]
  7. Yang, F.; Niu, K.; Collins, C.G.; Yan, X.; Ji, Y.; Ling, N.; Zhou, X.; Du, G.; Guo, H.; Hu, S. Grazing practices affect the soil microbial community composition in a Tibetan alpine meadow. Land Degrad. Dev. 2019, 30, 49–59. [Google Scholar] [CrossRef] [Green Version]
  8. 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]
  9. Courtney, R.; Harris, J.A.; Pawlett, M. Microbial Community Composition in a Rehabilitated Bauxite Residue Disposal Area: A Case Study for Improving Microbial Community Composition. Restor. Ecol. 2014, 22, 798–805. [Google Scholar] [CrossRef]
  10. 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]
  11. Delgado-Baquerizo, M.; Oliverio Angela, M.; Brewer Tess, E.; Benavent-González, A.; Eldridge David, J.; Bardgett Richard, D.; Maestre Fernando, T.; Singh Brajesh, K.; Fierer, N. A global atlas of the dominant bacteria found in soil. Science 2018, 359, 320–325. [Google Scholar] [CrossRef] [Green Version]
  12. Li, D.; Liu, J.; Chen, H.; Zheng, L.; Wang, K. Soil microbial community responses to forage grass cultivation in degraded karst soils, Southwest China. Land Degrad. Dev. 2018, 29, 4262–4270. [Google Scholar] [CrossRef]
  13. Liu, C.; Jin, Y.; Hu, Y.; Tang, J.; Xiong, Q.; Xu, M.; Bibi, F.; Beng, K.C. Drivers of soil bacterial community structure and diversity in tropical agroforestry systems. Agric. Ecosyst. Environ. 2019, 278, 24–34. [Google Scholar] [CrossRef]
  14. Zheng, J.M.; Chen, X.Y.; Chen, L.G.; He, T.Y.; Rong, J.D.; Lin, Y.; Zheng, Y.S. Comprehensive evaluation of soil quality at different stand densities of Dendrocalamus minor var. amoenus plantations. Appl. Ecol. Environ. Res. 2020, 18, 5985–5996. [Google Scholar] [CrossRef]
  15. Wang, C.; Xue, L.; Dong, Y.; Wei, Y.; Jiao, R. Unravelling the functional diversity of the soil microbial community of Chinese fir plantations of different densities. Forests 2018, 9, 532. [Google Scholar] [CrossRef] [Green Version]
  16. Zhang, Y.; Liu, T.; Guo, J.; Tan, Z.; Dong, W.; Wang, H. Changes in the understory diversity of secondary Pinus tabulaeformis forests are the result of stand density and soil properties. Glob. Ecol. Conserv. 2021, 28, e01628. [Google Scholar] [CrossRef]
  17. Tripathi, B.M.; Edwards, D.P.; Mendes, L.W.; Kim, M.; Dong, K.; Kim, H.; Adams, J.M. The impact of tropical forest logging and oil palm agriculture on the soil microbiome. Mol. Ecol. 2016, 25, 2244–2257. [Google Scholar] [CrossRef]
  18. Settineri, G.; Mallamaci, C.; Mitrović, M.; Sidari, M.; Muscolo, A. Effects of different thinning intensities on soil carbon storage in Pinus laricio forest of Apennine South Italy. Eur. J. For. Res. 2018, 137, 131–141. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, Y.M.; Cao, Y. Response of tree regeneration and understory plant species diversity to stand density in mature Pinus tabulaeformis plantations in the hilly area of the Loess Plateau, China. Ecol. Eng. 2014, 73, 238–245. [Google Scholar] [CrossRef]
  20. Ding, X.; Xiao, J.H.; Li, L.; Conran, J.G.; Li, J. Congruent species delimitation of two controversial gold-thread nanmu tree species based on morphological and restriction site-associated DNA sequencing data. J. Syst. Evol. 2019, 57, 234–246. [Google Scholar] [CrossRef]
  21. Xiao, J.-H.; Ding, X.; Li, L.; Ma, H.; Ci, X.-Q.; van der Merwe, M.; Conran, J.G.; Li, J. Miocene diversification of a golden-thread nanmu tree species (Phoebe zhennan, Lauraceae) around the Sichuan Basin shaped by the East Asian monsoon. Ecol. Evol. 2020, 10, 10543–10557. [Google Scholar] [CrossRef]
  22. Zhu, Y.; An, W.; Peng, J.; Li, J.; Gu, Y.; Jiang, B.; Chen, L.; Zhu, P.; Yang, H. Genetic diversity of nanmu (Phoebe zhennan S. Lee. et F. N. Wei) breeding population and extraction of core collection using nSSR, cpSSR and phenotypic markers. Forests 2022, 13, 1320. [Google Scholar] [CrossRef]
  23. Yang, H.; An, W.; Gu, Y.; Peng, J.; Jiang, Y.; Li, J.; Chen, L.; Zhu, P.; He, F.; Zhang, F.; et al. Integrative metabolomic and transcriptomic analysis reveals the mechanism of specific color formation in Phoebe zhennan heartwood. Int. J. Mol. Sci. 2022, 23, 13569. [Google Scholar] [CrossRef] [PubMed]
  24. Fang, X.-M.; Christenson, L.M.; Wang, F.-C.; Zeng, J.-P.; Chen, F.-S. Pine caterpillar outbreak and stand density impacts on nitrogen and phosphorus dynamics and their stoichiometry in Masson pine (Pinus massoniana) plantations in subtropical China. Can. J. For. Res. 2016, 46, 601–609. [Google Scholar] [CrossRef]
  25. Wilson, D.S.; Puettmann, K.J. Density management and biodiversity in young Douglas-fir forests: Challenges of managing across scales. For. Ecol. Manag. 2007, 246, 123–134. [Google Scholar] [CrossRef]
  26. Zhou, T.; Wang, C.; Zhou, Z. Impacts of forest thinning on soil microbial community structure and extracellular enzyme activities: A global meta-analysis. Soil Biol. Biochem. 2020, 149, 107915. [Google Scholar] [CrossRef]
  27. Lu, C.; Feng, Y.; Huang, Y. Sample stem analysis and its applycation in construction of volume table. For. Resour. Manag. 2004, 6, 31–35. [Google Scholar]
  28. Zhang, Y.; Li, Z.; Hou, L.; Song, L.; Yang, H.; Sun, Q. Effects of stand density on understory species diversity and soil nutrients in Chinese fir plantation. Acta Pedol. Sin. 2020, 57, 239–250. [Google Scholar]
  29. Fang, J.; Wang, X.; Shen, Z.; Tang, Z.; He, J.; Yu, D.; Jiang, Y.; Wang, Z.; Zheng, C.; Zhu, J.; et al. Methods and protocols for plant community inventory. Biodivers. Sci. 2009, 17, 533–548. [Google Scholar]
  30. Magurran, A. Ecological Diversity and Its Measurement; Princeton University Press: Princeton, NJ, USA, 1988. [Google Scholar]
  31. Lu, R. Analtical Methods of Soil Agricultural Chemistry; China Agricultural Science and Technology Press: Beijing, China, 1999. [Google Scholar]
  32. Soergel, D.A.W.; Dey, N.; Knight, R.; Brenner, S.E. Selection of primers for optimal taxonomic classification of environmental 16S rRNA gene sequences. ISME J. 2012, 6, 1440–1444. [Google Scholar] [CrossRef] [Green Version]
  33. Lam, T.Y.C.; Mei, R.; Wu, Z.; Lee, P.K.H.; Liu, W.-T.; Lee, P.-H. Superior resolution characterisation of microbial diversity in anaerobic digesters using full-length 16S rRNA gene amplicon sequencing. Water Res. 2020, 178, 115815. [Google Scholar] [CrossRef]
  34. Schloss Patrick, D.; Westcott Sarah, L.; Ryabin, T.; Hall Justine, R.; Hartmann, M.; Hollister Emily, B.; Lesniewski Ryan, A.; Oakley Brian, B.; Parks Donovan, H.; Robinson Courtney, J.; et al. Introducing mothur: Open-Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial Communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [Green Version]
  35. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  36. Edgar, R.C.; Haas, B.J.; Clemente, J.C.; Quince, C.; Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 2011, 27, 2194–2200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef]
  38. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glöckner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  39. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  40. Louca, S.; Parfrey Laura, W.; Doebeli, M. Decoupling function and taxonomy in the global ocean microbiome. Science 2016, 353, 1272–1277. [Google Scholar] [CrossRef]
  41. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [Green Version]
  42. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D. vegan: Community Ecology Package. R Package Version 2.4-1. 2016. Available online: https://CARAN.R-project.org/package=vegan (accessed on 12 October 2022).
  43. Will, R.E.; Narahari, N.V.; Shiver, B.D.; Teskey, R.O. Effects of planting density on canopy dynamics and stem growth for intensively managed loblolly pine stands. For. Ecol. Manag. 2005, 205, 29–41. [Google Scholar] [CrossRef]
  44. Rodríguez-Loinaz, G.; Onaindia, M.; Amezaga, I.; Mijangos, I.; Garbisu, C. Relationship between vegetation diversity and soil functional diversity in native mixed oak forests. Soil Biol. Biochem. 2008, 40, 49–60. [Google Scholar] [CrossRef]
  45. Lei, J.; Du, H.; Duan, A.; Zhang, J. Effect of stand density and soil layer on soil nutrients of a 37-year-old Cunninghamia lanceolata plantation in naxi, Sichuan Province, China. Sustainability 2019, 11, 5410. [Google Scholar] [CrossRef] [Green Version]
  46. Zhou, Y.; Li, J.H.; Ross, F.C.; Wang, H.F. Variation of soil bacterial communities in a chronosequence of Rubber tree (Hevea brasiliensis) plantations. Front. Plant Sci. 2017, 8, 849. [Google Scholar] [CrossRef] [Green Version]
  47. Bagchi, R.; Swinfield, T.; Gallery, R.E.; Lewis, O.T.; Gripenberg, S.; Narayan, L.; Freckleton, R.P. Testing the Janzen-Connell mechanism: Pathogens cause overcompensating density dependence in a tropical tree. Ecol. Lett. 2010, 13, 1262–1269. [Google Scholar] [CrossRef]
  48. Otsing, E.; Barantal, S.; Anslan, S.; Koricheva, J.; Tedersoo, L. Litter species richness and composition effects on fungal richness and community structure in decomposing foliar and root litter. Soil Biol. Biochem. 2018, 125, 328–339. [Google Scholar] [CrossRef]
  49. 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]
  50. Miki, T. Microbe-mediated plant–soil feedback and its roles in a changing world. Ecol. Res. 2012, 27, 509–520. [Google Scholar] [CrossRef]
  51. Miki, T.; Ushio, M.; Fukui, S.; Kondoh, M. Functional diversity of microbial decomposers facilitates plant coexistence in a plant–microbe–soil feedback model. Proc. Natl. Acad. Sci. USA 2010, 107, 14251–14256. [Google Scholar] [CrossRef]
  52. Navarrete, A.A.; Venturini, A.M.; Meyer, K.M.; Klein, A.M.; Tiedje, J.M.; Bohannan, B.J.M.; Nusslein, K.; Tsai, S.M.; Rodrigues, J.L.M. Differential response of Acidobacteria subgroups to forest-to-pasture conversion and their biogeographic patterns in the Western Brazilian Amazon. Front. Mcrobiol. 2015, 6, 1443. [Google Scholar] [CrossRef] [Green Version]
  53. Fierer, N.; Bradford, M.A.; Jackson, R.B. Toward an ecological classification of soil bacteria. Ecology 2007, 88, 1354–1364. [Google Scholar] [CrossRef]
  54. Ward Naomi, L.; Challacombe Jean, F.; Janssen Peter, H.; Henrissat, B.; Coutinho Pedro, M.; Wu, M.; Xie, G.; Haft Daniel, H.; Sait, M.; Badger, J.; et al. Three genomes from the phylum Acidobacteria provide insight into the lifestyles of these microorganisms in soils. Appl. Environ. Microbiol. 2009, 75, 2046–2056. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Eo, J.; Kim, M.-H.; Kim, M.-K.; Choi, S.-K. Shift of dominant species in plant community and soil chemical properties shape soil bacterial community characteristics and putative functions: A case study on topographic variation in a mountain pasture. Microorganisms 2021, 9, 961. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Principal coordinate analysis using the UniFrac metric (weighted and unweighted) to measure the phylogenetic community distance of bacterial communities.
Figure 1. Principal coordinate analysis using the UniFrac metric (weighted and unweighted) to measure the phylogenetic community distance of bacterial communities.
Forests 14 01612 g001
Figure 2. Bacterial community patterns at the phylum (A) and genus (C) levels in different stand densities. Note: * indicates a significant difference at the level of 0.05. (B,D) represent the comparison of differences in specific phylum and genus at different stand densities.
Figure 2. Bacterial community patterns at the phylum (A) and genus (C) levels in different stand densities. Note: * indicates a significant difference at the level of 0.05. (B,D) represent the comparison of differences in specific phylum and genus at different stand densities.
Forests 14 01612 g002
Figure 3. The taxonomic differences among different stand densities and specific networks at the genus level. Note: The top 50 correlated genera are presented in the figure. The size of the circle represents the abundance, the edges represent the correlation between the two species, the thickness of the edge represents the strength of the correlation, and with the color of the line, orange represents a positive correlation, while green represents a negative correlation.
Figure 3. The taxonomic differences among different stand densities and specific networks at the genus level. Note: The top 50 correlated genera are presented in the figure. The size of the circle represents the abundance, the edges represent the correlation between the two species, the thickness of the edge represents the strength of the correlation, and with the color of the line, orange represents a positive correlation, while green represents a negative correlation.
Forests 14 01612 g003
Figure 4. Putative functions of soil bacteria at the genus level showed the top 10 relative abundances according to FAPROTAX and their correlation with soil chemical properties. Note: (A) the relative abundance of bacteria; (B) Spearman’s correlation coefficient, * and ** indicate significant correlations at p ˂ 0.05 and p ˂ 0.01, respectively. SM: soil moisture, SBD: soil bulk density, SOC: soil oxidizable organic carbon, TN: total nitrogen, NH 4 + - N : ammonium nitrogen, NO 3 - N : nitrate-nitrogen, TP: total phosphorus, AP: available phosphorus, the same as below.
Figure 4. Putative functions of soil bacteria at the genus level showed the top 10 relative abundances according to FAPROTAX and their correlation with soil chemical properties. Note: (A) the relative abundance of bacteria; (B) Spearman’s correlation coefficient, * and ** indicate significant correlations at p ˂ 0.05 and p ˂ 0.01, respectively. SM: soil moisture, SBD: soil bulk density, SOC: soil oxidizable organic carbon, TN: total nitrogen, NH 4 + - N : ammonium nitrogen, NO 3 - N : nitrate-nitrogen, TP: total phosphorus, AP: available phosphorus, the same as below.
Forests 14 01612 g004
Figure 5. Redundancy analysis (RDA) results of bacterial and understory species communities. Notes: The scales on the horizontal and vertical coordinates are the values generated by the regression analysis calculation of soil physicochemical attributes for each understory species or bacterium at the genus level. The dot represents the plot, the blue dotted line represents the genus of bacteria with the top 10 abundance (A) and the herb (B) or shrub species (C), and the black arrows represent different soil physicochemical attributes.
Figure 5. Redundancy analysis (RDA) results of bacterial and understory species communities. Notes: The scales on the horizontal and vertical coordinates are the values generated by the regression analysis calculation of soil physicochemical attributes for each understory species or bacterium at the genus level. The dot represents the plot, the blue dotted line represents the genus of bacteria with the top 10 abundance (A) and the herb (B) or shrub species (C), and the black arrows represent different soil physicochemical attributes.
Forests 14 01612 g005
Figure 6. Pearson’s correlation between soil physicochemical attributes, understory, and microbial diversity index. Note: (A) represents the correlation between soil physicochemical attributes and microbial diversity index, (B) represents the correlation between soil physicochemical attributes and diversity index of herb, (C) represents the correlation between soil physicochemical attributes and diversity index of shrub, (D) represents the correlation between microbial diversity index and diversity index of herb, (E) represents the correlation between microbial diversity index and diversity index of shrub. H_herb: Shannon index of herb, Jsw_herb: Pielou index of herb, C_herb: Simpson index of herb, DMG_herb: Margalef richness index of herb, H_shrub: Shannon index of shrub, Jsw_shrbu: Pielou index of shrub, C_shrub: Simpson index of shrub, DMG_shrub: Margalef richness index of shrub, SM: soil moisture, SBD: soil bulk density, SOC: soil oxidizable organic carbon, TN: total nitrogen, NH 4 + - N : ammonium nitrogen, NO 3 - N : nitrate-nitrogen, TP: total phosphorus, AP: available phosphorus, the same as below. * and ** indicate significant correlations at p ˂ 0.05 and p ˂ 0.01, respectively.
Figure 6. Pearson’s correlation between soil physicochemical attributes, understory, and microbial diversity index. Note: (A) represents the correlation between soil physicochemical attributes and microbial diversity index, (B) represents the correlation between soil physicochemical attributes and diversity index of herb, (C) represents the correlation between soil physicochemical attributes and diversity index of shrub, (D) represents the correlation between microbial diversity index and diversity index of herb, (E) represents the correlation between microbial diversity index and diversity index of shrub. H_herb: Shannon index of herb, Jsw_herb: Pielou index of herb, C_herb: Simpson index of herb, DMG_herb: Margalef richness index of herb, H_shrub: Shannon index of shrub, Jsw_shrbu: Pielou index of shrub, C_shrub: Simpson index of shrub, DMG_shrub: Margalef richness index of shrub, SM: soil moisture, SBD: soil bulk density, SOC: soil oxidizable organic carbon, TN: total nitrogen, NH 4 + - N : ammonium nitrogen, NO 3 - N : nitrate-nitrogen, TP: total phosphorus, AP: available phosphorus, the same as below. * and ** indicate significant correlations at p ˂ 0.05 and p ˂ 0.01, respectively.
Forests 14 01612 g006
Table 1. Characteristics of Phoebe zhennan plantations of different stand densities.
Table 1. Characteristics of Phoebe zhennan plantations of different stand densities.
NumberStand Density (stems/hm2)H (m)DBH (cm)P (m)V (m3)Altitude (m)Aspect and Slope
D14007.4 ± 0.212.4 ± 0.63.4 ± 0.40.050 ± 0.005490S 8°
D25007.4 ± 0.212.4 ± 0.03.4 ± 0.10.051 ± 0.002492S 9°
D36007.8 ± 0.312.8 ± 0.93.5 ± 0.30.056 ± 0.009493S 8°
D47007.4 ± 0.112.5 ± 0.63.3 ± 0.20.051 ± 0.004505S 8°
D58507.5 ± 0.112.0 ± 0.73.1 ± 0.20.048 ± 0.006510S 9°
Note: H—height, DBH—diameter at breast height, P—crown diameter, V—volume, S—the southern slope.
Table 2. Species diversity indices of the shrub and herb layers under five P. zhennan plantation stand densities.
Table 2. Species diversity indices of the shrub and herb layers under five P. zhennan plantation stand densities.
Stand DensitiesHerb LayerShrub Layer
H′JswDDMGH′JswDDMG
D11.86 ± 0.140.69 ± 0.050.25 ± 0.042.32 ± 0.261.00 ± 0.470.87 ± 0.08 a0.44 ± 0.18 b1.12 ± 0.69
D22.02 ± 0.260.74 ± 0.090.21 ± 0.092.58 ± 0.360.59 ± 0.610.47 ± 0.44 ab0.67 ± 0.34 ab0.52 ± 0.48
D31.93 ± 0.170.75 ± 0.080.22 ± 0.062.20 ± 0.580.16 ± 0.280.15 ± 0.25 b0.92 ± 0.14 a0.25 ± 0.43
D41.76 ± 0.250.67 ± 0.120.24 ± 0.062.24 ± 0.381.22 ± 0.260.89 ± 0.06 a0.34 ± 0.10 b1.30 ± 0.52
D51.38 ± 0.340.51 ± 0.110.40 ± 0.152.39 ± 0.181.29 ± 0.180.78 ± 0.05 a0.37 ± 0.05 b1.41 ± 0.39
F value3.1463.3872.5350.4813.4005.5174.1032.808
p value0.0640.0540.1060.7490.0530.0130.0320.085
Note: Data are reported as the mean ± standard error (n = 3). H′—Shannon index, Jsw—Pielou index, D—Simpson index, DMG—Margalef richness index. Within a column, values followed by the same lowercase letter indicate that they did not differ significantly (p ˂ 0.05).
Table 3. Soil physicochemical characteristics in Phoebe zhennan plantations of different stand densities.
Table 3. Soil physicochemical characteristics in Phoebe zhennan plantations of different stand densities.
Stand DensitiespHSBD (g/cm3)SM (%)TN (g/kg) NH 4 + - N NO 3 - N TP (mg/kg)AP (mg/kg)SOC (g/kg)
D15.49 ± 0.481.76 ± 0.150.17 ± 0.0011.63 ± 1.65 a6.38 ± 1.03 b8.45 ± 0.61 b0.37 ± 0.022.13 ± 0.5923.73 ± 3.86 a
D25.31 ± 0.131.77 ± 0.070.17 ± 0.039.31 ± 0.27 b6.70 ± 0.80 b8.36 ± 0.48 b0.30 ± 0.022.01 ± 0.2521.26 ± 1.00 a
D35.26 ± 0.511.72 ± 0.180.15 ± 0.0411.44 ± 1.28 a3.63 ± 0.77 c11.31 ± 0.84 a0.30 ± 0.102.36 ± 0.8923.53 ± 3.22 a
D45.85 ± 0.121.83 ± 0.080.17 ± 0.029.43 ± 0.71 b6.50 ± 0.76 b7.16 ± 0.21 b0.25 ± 0.051.79 ± 0.1716.25 ± 1.89 b
D55.08 ± 0.291.76 ± 0.140.18 ± 0.026.23 ± 0.36 c9.67 ± 0.10 a12.84 ± 1.67 a0.30 ± 0.112.76 ± 0.2120.27 ± 2.43 ab
F value2.1490.2640.53514.06023.70020.1201.0151.5843.969
p value0.1490.8940.7130.0000.0000.0000.4450.2530.035
Note: Significant differences between the five density soils were determined using a one-way ANOVA at p ˂ 0.05. The data are shown as the mean ± SD (n = 3). SM: soil moisture, SBD: soil bulk density, SOC: soil oxidizable organic carbon, TN: total nitrogen, NH 4 + - N : ammonium nitrogen, NO 3 - N : nitrate nitrogen, TP: total phosphorus, AP: available phosphorus. Within a column, values followed by the same lowercase letter indicate that they did not differ significantly (p ˂ 0.05).
Table 4. Richness and diversity indices of bacterial communities in different stand densities.
Table 4. Richness and diversity indices of bacterial communities in different stand densities.
Stand DensitiesACEChao1SimpsonShannonCoverage
D11234.57 ± 85.041244.27 ± 86.370.99 ± 0.008.44 ± 0.120.97 ± 0.00
D21357.93 ± 86.651345.38 ± 92.840.99 ± 0.008.52 ± 0.100.97 ± 0.00
D31234.59 ± 177.421238.40 ± 172.100.99 ± 0.008.52 ± 0.270.97 ± 0.01
D41365.42 ± 54.391355.03 ± 56.330.99 ± 0.008.58 ± 0.140.97 ± 0.00
D51141.21 ± 55.341135.91 ± 47.420.99 ± 0.008.27 ± 0.090.98 ± 0.00
F value2.5762.3740.9161.7951.469
p value0.1030.1220.4920.2060.283
The data are shown as the mean ± SD (n = 3).
Table 5. Relative abundance of the top 20 bacterial genera in different stand densities.
Table 5. Relative abundance of the top 20 bacterial genera in different stand densities.
PhylumClassOrderFamilyGenusD1D2D3D4D5F Valuep Value
VerrucomicrobiaVerrucomicrobiaeChthoniobacteralesChthoniobacteraceaeCandidatus_Udaeobacter23.8727.4819.6825.6620.421.0740.419
PedosphaeralesPedosphaeraceaeADurb.Bin063_17.006.496.155.835.070.8680.516
ChthoniobacteralesXiphinematobacteraceaeCandidatus_Xiphinematobacter1.042.031.582.076.283.9300.036
ProteobacteriaAlphaproteobacteriaRhizobialesXanthobacteraceaeBradyrhizobium4.994.596.115.625.590.3290.852
Pseudolabrys2.812.282.012.511.681.1010.408
SphingomonadalesSphingomonadaceaeSphingomonas3.522.873.273.471.782.8250.083
ReyranellalesReyranellaceaeReyranella1.221.080.850.851.271.0940.411
GammaproteobacteriaBurkholderialesBurkholderiaceaeParaburkholderia1.411.471.591.683.622.4600.113
BetaproteobacterialesNitrosomonadaceaeEllin60671.762.251.181.881.252.8960.079
XanthomonadalesRhodanobacteraceaeRhodanobacter1.071.181.451.321.550.4920.742
DeltaproteobacteriaMyxococcalesHaliangiaceaeHaliangium1.461.471.531.560.792.0470.163
AcidobacteriaBlastocatellia_Subgroup_4PyrinomonadalesPyrinomonadaceaeRB412.962.843.646.960.602.4460.115
AcidobacteriiaAcidobacterialesAcidobacteriaceae_Subgroup_1Occallatibacter2.372.363.211.184.862.1390.150
Edaphobacter1.811.541.761.073.051.9580.177
KoribacteraceaeCandidatus_Koribacter2.161.542.131.222.671.1080.405
SolibacteralesSolibacteraceae_Subgroup_3Bryobacter2.422.353.152.472.890.5350.713
Candidatus_Solibacter8.877.449.557.559.540.2830.882
NitrospirotaNitrospiriaNitrospiralesNitrospiraceaeNitrospira1.311.432.122.200.991.4460.289
BacteroidetesBacteroidiaFlavobacterialesFlavobacteriaceaeFlavobacterium1.022.221.711.270.942.1520.149
GemmatimonadetesGemmatimonadetesGemmatimonadalesGemmatimonadaceaeGemmatimonas1.311.011.520.931.240.5020.735
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cheng, Y.; Peng, J.; Gu, Y.; Guo, H.; Jiang, T.; Yang, H. The Early Effect of Plant Density on Soil Physicochemical Attributes and Bacterial and Understory Plant Diversity in Phoebe zhennan Plantations. Forests 2023, 14, 1612. https://doi.org/10.3390/f14081612

AMA Style

Cheng Y, Peng J, Gu Y, Guo H, Jiang T, Yang H. The Early Effect of Plant Density on Soil Physicochemical Attributes and Bacterial and Understory Plant Diversity in Phoebe zhennan Plantations. Forests. 2023; 14(8):1612. https://doi.org/10.3390/f14081612

Chicago/Turabian Style

Cheng, Yilun, Jian Peng, Yunjie Gu, Hongying Guo, Tianyi Jiang, and Hanbo Yang. 2023. "The Early Effect of Plant Density on Soil Physicochemical Attributes and Bacterial and Understory Plant Diversity in Phoebe zhennan Plantations" Forests 14, no. 8: 1612. https://doi.org/10.3390/f14081612

APA Style

Cheng, Y., Peng, J., Gu, Y., Guo, H., Jiang, T., & Yang, H. (2023). The Early Effect of Plant Density on Soil Physicochemical Attributes and Bacterial and Understory Plant Diversity in Phoebe zhennan Plantations. Forests, 14(8), 1612. https://doi.org/10.3390/f14081612

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