3.1. Spatial distribution of soil factors
Table 2 shows that the soil pH of the sample plots under the five forest stands in Jinyun Mountain are all acidic, and the soil pH of GA-EBF, CLF, PPF, PM&GA, PM&SS are 4.1, 4.0, 4.5, 4.1, 4.0, respectively. Moreover, the differences in soil pH of GA-EBF, CLF, PM&GA, PM&SS were not significant (P < 0.05), while the pH value of PPF (4.5) was higher than that of the other four forest stands with a statistically significant difference (P > 0.05). Furthermore, the soils under the five forests have no significant differences in the contents of soils organic carbon (SOC), total nitrogen (TN), Fed and Ald, and Feo and Alo. As far as the content of available phosphorus (AP) is concerned, the performance of different forest stands is PPF > CLF > PM&SS > PM&GA > GA-EBF, and there are significant differences between GA-EBF and PPF, CLF and PM&GA, PPF and PM&SS. In terms of cation exchange capacity (CEC), GA-EBF > CLF > M&GA > PPF > PM&SS. Among them, the CEC of PPF is significantly different from that of GA-EBF and PM&SS.
The soil texture analysis results showed that the five forest stands had significant differences in the proportion of clay, silt and sand. In addition, the proportions of sand in the soil are all the largest for the five forest stands, especially the proportion of sand in GA is as high as 70%.
3.3. Structural characteristics of SMC
Variation of soil PLFA biomarkers. 53, 62, 69, 61 and 64 phospholipid fatty acids were detected in the soil samples of the sample plots in JMNR, and 38 PLFA biomarkers greater than 0.01 nmol·g-1 were selected for further analysis (Table 4). The investigation of the spatial distribution of microorganisms by Zhou et al. (2015) showed significant differences in soil microorganisms in different locations. For this study, some PLFA biomarkers are completely distributed in five sample plots, such as 12:0 anteiso, 12:0, 13:0 anteiso, 14:0 anteiso, 15:0 anteiso, 16:1 ω5c, 17: 1 ω7c 10-methyl, 18:1 ω9c, 19:0 cyclo ω7c; some biomarkers are only distributed in specific single sample plot and are incomplete distribution, such as 11:0 anteiso, 14:0 iso, 15:1 ω5c, 16:1 ω9c, 17:1 iso ω9c, 18:1 ω7c 10-methyl, 20:0.
Table 4
Content of PLFA biomarkers in different forest stands of Jinyun Mountain1 (nmol·g-1)
No | PLFA biomarkers | Microorganisms | GA | CLF | PPF | PM&GA | PM&SS |
1 | 11:0 anteiso | G+ | — | — | 0.096 ± 0.043ab | 0.177 ± 0.076bc | 0.175 ± 0.079c |
2 | G+ | G+ | 0.144 ± 0.078a | 0.191 ± 0.061a | 0.178 ± 0.044a | 0.243 ± 0.080a | 0.201 ± 0.091a |
3 | 12:0 | Bacteria | 0.122 ± 0.048a | 0.130 ± 0.024a | 0.109 ± 0.039a | 0.166 ± 0.071a | 0.211 ± 0.049a |
4 | 13:0 anteiso | G+ | 0.154 ± 0.028a | 0.158 ± 0.042a | 0.200 ± 0.040a | 0.207 ± 0.050a | 0.182 ± 0.045a |
5 | 14:0 iso | G+ | 0.047 ± 0.006ab | 0.076 ± 0.027b | — | — | 0.086 ± 0.046b |
6 | 14:0 anteiso | G+ | 0.178 ± 0.017a | 0.155 ± 0.058a | 0.237 ± 0.054a | 0.200 ± 0.021a | 0.196 ± 0.023a |
7 | 14:0 | Bacteria | 0.230 ± 0.067a | 0.343 ± 0.156a | 0.252 ± 0.090a | 0.317 ± 0.170a | 0.472 ± 0.291a |
8 | 15:0 iso | G+ | 0.973 ± 0.293a | 1.177 ± 0.720a | 1.022 ± 1.005a | 1.322 ± 0.709a | 1.563 ± 0.776a |
9 | 15:0 anteiso | G+ | 0.497 ± 0.111a | 0.531 ± 0.199a | 0.469 ± 0.338a | 0.568 ± 0.224a | 0.646 ± 0.205a |
10 | 15:1 ω7c | G- | 0.119 ± 0.018a | 0.118 ± 0.021a | 0.119 ± 0.021a | 0.107 ± 0.024a | 0.121 ± 0.023a |
11 | 15:1 ω5c | G- | 0.649 ± 0.022b | 0.336 ± 0.300a | — | 0.666 ± 0.041b | 0.699 ± 0.075b |
12 | 15:0 | Bacteria | 0.130 ± 0.033a | 0.156 ± 0.093a | 0.131 ± 0.064a | 0.182 ± 0.109a | 0.213 ± 0.094a |
13 | 15:0 DMA | Anaerobic bacteria | 0.117 ± 0.042a | 0.146 ± 0.048a | 0.172 ± 0.084a | 0.192 ± 0.063a | 0.221 ± 0.089a |
14 | 16:3 ω6c | Eukaryote | 0.139 ± 0.044a | 0.128 ± 0.042a | 0.087 ± 0.033a | 0.165 ± 0.036a | 0.141 ± 0.034a |
15 | 16:0 iso | G+ | 0.583 ± 0.199a | 0.652 ± 0.420a | 0.717 ± 0.469a | 0.945 ± 0.676a | 1.136 ± 0.607b |
16 | 16:0 anteiso | G+ | 0.163 ± 0.015a | 0.173 ± 0.019a | 0.210 ± 0.048a | 0.175 ± 0.012a | 0.174 ± 0.010a |
17 | 16:1 ω9c | G- | 0.079 ± 0.022a | 0.117 ± 0.029a | — | 0.084 ± .036a | 0.081 ± 0.034a |
18 | 16:1 ω7c | G- | 0.263 ± 0.098a | 0.397 ± 0.238a | 0.463 ± 0.306b | 0.358 ± 0.198a | 0.321 ± 0.114a |
19 | 16:1 ω5c | AM fungi | 0.261 ± 0.103a | 0.295 ± 0.155b | 0.363 ± 0.205b | 0.293 ± 0.174a | 0.351 ± 0.109b |
20 | 16:0 | G- | 3.131 ± 0.691a | 3.711 ± 1.463b | 2.286 ± 1.562a | 4.024 ± 1.313c | 4.571 ± 0.940c |
21 | 16:0 10-methyl | Actinomycetes | 0.913 ± 0.311a | 0.870 ± 0.510a | 0.855 ± 0.901a | 1.245 ± 0.680a | 1.601 ± 0.589a |
22 | 17:1 iso ω9c | G+ | — | 1.102 ± 0.120b | 1.023 ± 0.345b | — | — |
23 | 17:0 iso | G+ | 0.284 ± 0.104a | 0.379 ± 0.227a | 0.291 ± 0.254a | 0.357 ± 0.167a | 0.502 ± 0.211a |
24 | 17:0 anteiso | G+ | 0.346 ± 0.053a | 0.447 ± 0.115a | 0.369 ± 0.208a | 0.392 ± 0.082a | 0.574 ± 0.275a |
25 | 17:0 cyclo ω7c | G- | 0.183 ± 0.059a | 0.269 ± 0.161a | 0.254 ± 0.155a | 0.287 ± 0.129a | 0.353 ± 0.118a |
26 | 17:0 | Bacteria | 0.094 ± 0.038b | 0.111 ± 0.064ab | — | 0.121 ± 0.052ab | 0.145 ± 0.042b |
27 | 17:1 ω7c 10-methyl | Actinomycetes | 0.224 ± 0.064a | 0.309 ± 0.112b | 0.267 ± 0.164a | 0.275 ± 0.108a | 0.373 ± 0.166b |
28 | 17:0 10-methyl | Actinomycetes | 0.181 ± 0.083a | 0.200 ± 0.084a | 0.172 ± 0.118a | 0.295 ± 0.128a | 0.243 ± 0.103a |
29 | 18:3 ω6c | G- | 0.198 ± 0.036a | 0.274 ± 0.077a | 0.192 ± 0.055a | 0.265 ± 0.083a | 0.214 ± 0.038a |
30 | 18:2 ω6c | Fungi | 0.392 ± 0.074a | 0.477 ± 0.193a | 0.279 ± 0.119a | 0.449 ± 0.217a | 0.513 ± 0.097a |
31 | 18:1 ω9c | Fungi | 0.826 ± 0.343a | 0.897 ± 0.560a | 0.830 ± 0.490a | 1.017 ± 0.557a | 1.052 ± 0.281a |
32 | 18:1 ω7c | G- | 0.469 ± 0.204a | 0.757 ± 0.442b | 0.512 ± 0.283a | 0.624 ± 0.356a | 0.505 ± 0.250a |
33 | 18:0 | Bacteria | 2.113 ± 0.411a | 2.317 ± 0.453a | 1.727 ± 0.616a | 2.397 ± 0.311a | 2.776 ± 0.298a |
34 | 18:1 ω7c 10-methyl | Actinomycetes | — | 0.123 ± 0.058b | — | 0.138 ± 0.043b | 0.118 ± 0.063b |
35 | 18:0 10-methyl | Actinomycetes | 0.177 ± 0.083ab | 0.280 ± 0.144c | — | 0.171 ± 0.042ab | 0.259 ± 0.118c |
36 | 19:0 cyclo ω7c | G- | 1.853 ± 0.744a | 1.815 ± 1.016a | 1.844 ± 1.189a | 2.169 ± 1.099a | 3.108 ± 1.589a |
37 | 20:0 | Bacteria | 0.059 ± 0.014b | — | — | 0.080 ± 0.022b | 0.136 ± 0.029c |
38 | 20:0 10-methyl | Actinomycetes | 0.174 ± 0.014a | 0.295 ± 0.155a | 0.236 ± 0.126a | 0.339 ± 0.092a | 0.287 ± 0.135a |
Total content of different PLFA biomarkers | 16.464 | 20.010 | 15.963 | 21.013 | 24.518 |
1 Different lowercase letters indicate the significance of the content differences of specific PLFA biomarker for different forest sample plots (P < 0.05); "—" means that the PLFA biomarker is not detected. |
For the soils of the five forest sample plots, the most abundant PLFA biomarkers are 16:0 (indicator of G-bacteria), 18:0 (indicator of bacteria), 19:0 cyclo ω7c (indicator of G-bacteria), 15 :0 iso (indicator of G + bacteria) and 16:0 10-methyl (indicator of actinomycetes), indicating that their major roles in different soils of the five forest sample plots. Moreover, the contents of the four PLFA biomarkers 19:0 cyclo ω7c, 18:0, 15:0 iso, and 16:0 10-methyl are not significantly different in the sample plots, whose content of PM&SS is the highest, followed by PM&GA. The same situation applies to the PLFA biomarker of 16:0. Regarding the changes in phospholipid fatty acids, the contents of the five PLFA biomarkers (16:0, 18:0, 19:0 cyclo ω7c, 15:0 iso and 16:0 10-methyl) in GA are 3.131, 2.113, 1.853, 0.973 and 0.913 nmol·g-1, accounting for 54.56% of the total fatty acid content of GA. The contents of 16:0, 18:0, 19:0 cyclo ω7c, 15:0 iso and 16:0 10-methyl in CLF are 3.711, 2.317, 1.815, 1.177 and 0.870 nmol·g-1, accounting for 49.42% of the total fatty acids in CLF, and their contents in PPF are 2.286, 1.727, 1.844, 1.022 and 0.855 nmol·g-1, accounting for 48.45% of the total fatty acids in PPF. The contents of 16:0, 18:0, 19:0 cyclo ω7c, 15:0 iso and 16:0 10-methyl in PM&GA and PM&SS are relatively high: their contents in PM&GA are 4.024, 2.397, 2.169, 1.322 and 1.245 nmol·g-1, accounting for 55.54% of its total fatty acids; while their contents in PM&SS are 4.571, 2.776, 3.108, 1.563 and 1.601 nmol·g-1, accounting for 52.45% of its total fatty acids. Moreover, the order of the total PLFA biomarkers of the five forest stands is as follows: PM&SS > PM&GA > CLF > PM&GA > PPF. In short, different forest stands have similar types of fatty acid with high content, but the distribution of PLFA biomarker content in each plot shows that the soil microbial PLFAs have certain spatial distribution characteristics for the forest stands in JMNR.
PCA of SMC structure. According to the PCA of SMC in different forest stand sample plots, the cumulative contribution rate of the two principal components related to the diversity of soil microbial PLFA community reached 68.62% (Fig. 2). Specifically, the first principal component (PC1) and the second principal component (PC2) explain 46.95% and 21.67% of the variance of the variables, respectively. Furthermore, 21 microbial PLFAs that play a significant role in PC1 have been identified, where 17:1 iso ω9c and 20:0 10-methyl are positively correlated with PC1, while 20:0, 15:1 ω5c, 11:0 anteiso, 17:0, 18:0 10-methyl, 16:1 ω9c, 16:3 ω6c, 14:0 iso, 12:0 and 18:0 10-methyl are negatively correlated with PC1. The microbial PLFAs that play a significant role in PC2 are 16:0 anteiso, 15:0 anteiso, 15:1 ω7c, 18:1 ω7c, 16:1 ω5c, 13:0 anteiso and 18:3 ω6c, all of which are positively correlated with PC2. The plots of the five forest stands are clearly distinguished. GA is mainly located at the negative end of PC1 and PC2, CLF is basically located at the positive end of PC1 and PC2, PPF is located at the positive end of PC1 and the positive and negative ends of PC2, and PM&GA and PM&SS are both located at the negative end of PC1 and the positive and negative ends of PC2. It is suggested that the main PLFAs affecting soil microorganisms in the five forest stands of JMNR are 17:1 iso ω9c, 20:0 10-methyl, 16:0 anteiso, 15:0 anteiso, 15:1 ω7c, 18:1 ω7c, 16: 1 ω5c, 13:0 anteiso and 18:3 ω6c.
a.11:0 anteiso,b.12:0 anteiso,c.12:0,d.13:0 anteiso,e.14:0 iso,f.14:0 anteiso,g.14:0,h.15:0 iso,i.15:0 anteiso,j.15:1 ω7c,k.15:1 ω5c,l.15:0,m.15:0 DMA,n.16:3 ω6c,o.16:0 iso,p.16:0 anteiso,q.16:1 ω9c,r.16:1 ω7c,s.16:1 ω5c,t.16:0,u.16:0 10-methyl,v.17:1 iso ω9c,w.17:0 iso,x.17:0 anteiso,y. 17:0 cyclo ω7c,z.17:0,A.17:1 ω7c 10-methyl,B.17:0 10-methyl,C.18:3 ω6c,D.18:2 ω6c,E.18:1 ω9c,F.18:1 ω7c,G.18:0,H.18:1 ω7c 10-methyl,I.18:0 10-methyl,J.19:0 cyclo ω7c,K.20:0,L.20:0 10-methyl; Numbers 1 ~ 4 represent the four sampling points of each of the five forest stands.
3.5. Correlation between SMC diversity and soil factors
RDA is used to analyze the correlation between SMC diversity, soil physical and chemical properties, and root system characteristics (Fig. 3). The two ranking axes explain a total of 95.82% of the SMC structure differences, of which the first axis and the second axis explain 91.5% and 4.32% of the SMC structure differences, respectively. Figure 3 shows that G + and AM fungi are significantly positively correlated with SOC, Ald, Alo, and Fed, and actinomycetes and fungi are significantly positively correlated with SOC, CEC, TN, and Fed. The above four microbial groups were significantly negatively correlated with AP, Sand, and the proportion of fine roots (0.2-1). Moreover, G- is significantly positively correlated with SOC, CEC, and TN, and significantly negatively correlated with pH, RLD, the proportion of very fine roots (0-0.2), and AP. The G+/G- ratio is significantly positively correlated with RLD, pH, AP and the proportion of very fine roots, and significantly negatively correlated with AD, TN, CEC, SOC and the proportion of fine roots. The ratio of F/B was significantly positively correlated with TN and AD, and significantly negatively correlated with AP, pH, and the proportion of very fine roots.