Distribution characteristics of bacterial and fungal communities
Bacterial OUTs (Fig. 2, a): A total of 8760 OTUs were generated by OTUs clustering at 97% similarity level. Among them, 3787 OTUs were common in roots, root epidermis, rhizosheath, and rhizosphere soil. There were 188, 347, 576 and 711 unique OTUs in roots, root epidermis, rhizosheath, and rhizosphere soil, respectively, accounting for 3.68%, 5.98%, 8.48% and 10.06% of their respective total number of OTUs. Fungal OUTs (Fig. 2, b): A total of 1423 OTUs were generated by OTUs clustering at 97% similarity level. Among them, 318 OTUs were common in roots, root epidermis, rhizosheath soil and rhizosphere soil. The number of unique OTUs in roots, root epidermis, rhizosheath soil and rhizosphere soil were 11, 19, 164 and 262, respectively, accounting for 2.28%, 3.84%, 15.96% and 22.13% of their respective total numbers of OTUs. The dilution curves of OTUs tended to be gradual, indicating that the measured data could accurately reflect the information on the plant fungal community (Fig. 2, c, d).
Table 1 shows that the number of Observed_species, Chao1 index and ACE number index were highest in rhizosphere soil, followed by rhizosheath soil, and were smallest in the root system, and there were significant differences between them. There was no significant difference in the Shannon index between rhizosphere soil and rhizosheath soil, which were both significantly higher than the Shannon index in root epidermis and roots, but the Shannon index of bacteria in root epidermis was significantly higher than in roots. The Simpson index of bacteria in rhizosphere soil, rhizosheath soil and root epidermis had no significant differences, but were significantly higher than those in the root system. The goods_coverage index of the samples was higher than 98%.
The observed_species, Shannon, Simpon, Chao1, ACE, PD_whole_tree index of fungi in rhizosphere soil were the highest, followed by rhizosheath soil, but there was no significant difference between the two groups. The observed_species, Shannon, Simpon, Chao1, ACE, index of rhizosphere soil and rhizosheath soil fungi were significantly larger than the root epidermis and root. There was no significant difference in observed_species, Shannon, Simpon, Chao1, ACE index between roots and root epidermis. The good_coverage indexes of the samples were all higher than 99%.
Table 1
Alpha diversity of bacteria and fungi in roots, rhizosheath, rhizosphere and non-rhizosphere soil
Microorganism
|
Sample name
|
Observed_
species
|
Shannon
|
Simpson
|
Chao1
|
ACE
|
Goods_
coverage
|
Bacteria
|
SG.EC
|
2769.0 ± 143.9d
|
5.41 ± 0.31c
|
0.82 ± 0.03b
|
3210.4 ± 159.4d
|
3312.3 ± 174.8d
|
0.989 ± 0.001
|
SG.RS
|
3287.4 ± 213.2c
|
8.87 ± 0.46b
|
0.99 ± 0.01a
|
3714.6 ± 248.7c
|
3798.1 ± 266.4c
|
0.989 ± 0.001
|
SG.RH
|
3919.2 ± 187.9b
|
9.77 ± 0.28a
|
1.00 ± 0.00a
|
4392.5 ± 145.6b
|
4444.4 ± 156.3b
|
0.988 ± 0.000
|
SG.SO
|
4278.6 ± 256.5a
|
10.10 ± 0.10a
|
1.00 ± 0.00a
|
4858.4 ± 512.0a
|
4938.6 ± 504.8a
|
0.986 ± 0.003
|
Fungi
|
SG.EC
|
230.2 ± 27.5b
|
2.38 ± 0.59b
|
0.62 ± 0.14b
|
264.3 ± 34.5b
|
278.1 ± 37.7b
|
0.999 ± 0.000
|
SG.RS
|
223.8 ± 29.2b
|
2.61 ± 0.75b
|
0.68 ± 0.15b
|
260.8 ± 33.4b
|
272.6 ± 31.3b
|
0.999 ± 0.000
|
SG.RH
|
532.2 ± 65.5a
|
5.11 ± 0.62a
|
0.89 ± 0.06a
|
597.7 ± 65.4a
|
613.4 ± 68.8a
|
0.998 ± 0.000
|
SG.SO
|
563.0 ± 70.9a
|
5.64 ± 0.44a
|
0.93 ± 0.04a
|
693.9 ± 163.6a
|
656.6 ± 83.6a
|
0.998 ± 0.001
|
Note: Different lowercase letters indicate significant differences between bacterial or fungal groups (P<0.05). Observed_ species and Goods_ coverage represents sequencing depth index; Shannon and Simpson represent the diversity index of the bacterial community, and Chao1 and ACE represent the index of bacterial community richness; SG.EC, SG.RS, SG.RH, SG.SO refer to roots, root epidermis, rhizosheath, and rhizosphere soil of Stipa grandis, respectively.
The top 10 phylum level classification of bacterial communities is shown in Fig. 3 (a). The relative abundance of Cyanobacteria in the root system was 46.4%, while the relative abundance of Cyanobacteria in rhizosphere, rhizosphere and non-rhizosphere soil was only 5.0%, 2.1% and 1.1%, respectively. Actinobacteria was the dominant population in root epidermis, and its relative abundance was 43.0%, which was significantly higher than that in the root system, rhizosheath soil and rhizosphere soil. The relative abundance of Proteobacteria in root epidermis, rhizosheath soil and rhizosphere soil was 30.9%, 31.3% and 30.2%, respectively, while that in roots was only 11.7%. The relative abundances of Acidobacteria, Gemmatimonadetes, Bacteroidetes and Verrucomicrobia in rhizosheath soil and rhizosphere soil were similar and significantly higher than those in root systems and root epidermis. The relative abundances in rhizosheath soil were 22.1%, 9.2%, 8.2% and 3.2%, respectively, and those in rhizosphere soil were 21.8%, 9.6%, 7.6% and 3.7%, respectively.
The top 10 phylum level classification of the sample fungal communities is shown in Fig. 3 (b). Different from bacterial communities, the analysis shows that Basidiomycota and Ascomycota have higher relative abundance in roots, root epidermis, rhizosheath and rhizosphere soil. The relative abundance of Basidiomycota in roots and root epidermis are 72.3% and 70.3%, respectively, which is significantly higher than in rhizosheath soil and rhizosphere soil. The relative abundance of Ascomycota in roots, rhizosheath and rhizosphere soil was 20.4%, 22.4% and 27.6%, respectively, while that in root epidermis was only 8.7%. The relative abundance of Others was in the order rhizosphere soil > rhizosheath soil > root epidermis > roots.
Results of non-metric multidimensional scaling (NMDS) analysis of bacteria in roots, root epidermis, rhizosheath and rhizosphere soil of Stipa grandis are shown in Fig. 3 (c). The microbial community structure of root epidermis, rhizosheath soil and rhizosphere soil are distinct and distant in the plot, while the distance between rhizosheath soil and rhizosphere soil is relatively close. The UPGMA cluster analysis of all samples also showed that the similarity of bacterial composition and relative abundance between rhizosheath soil and rhizosphere soil was high, and that there was a big difference in bacterial composition between soil, roots and root epidermis.
The results of NMDS analysis of fungi are shown in Fig. 3 (d). The fungal community structure of roots and root epidermis, and rhizosheath and rhizosphere soil were relatively close, while that of roots and root epidermis were relatively distinct from those of rhizosheath and rhizosphere soil. Through UPGMA cluster analysis of all the samples, in contrast to bacterial community clustering, rhizosheath soil and rhizosphere soil fungal communities can be clustered into one group, while root and root epidermis can be clustered into another group.
As shown in Table 2, the bacterial community composition in rhizosheath soil and rhizosphere soil of Stipa grandis was not significantly different (P > 0.05). The bacterial community composition of roots and root epidermis of Stipa grandis was significantly different at P < 0.05, and the differences in bacterial community among other groups was extremely significantly different (P < 0.01).
There was no significant difference in fungal communities between the roots and root epidermis of Stipa grandis or between rhizosheath soil and rhizosphere soil (P > 0.05). The difference of fungal community among other groups reached an extremely significant level (P < 0.01).
Table 2
Significance test table of difference in bacterial and fungal community structure between groups
Microorganism
|
vs_group
|
SS
|
df
|
MS
|
Fs
|
p-value
|
Bacteria
|
SG.EC-SG.RH-SG.RS-SG.SO
|
2.18148(0.234448)
|
3(16)
|
0.727159(0.014653)
|
49.6253
|
< 0.001**
|
SG.EC-SG.RS
|
1.12673(0.0835403)
|
1(8)
|
1.12673(0.0104425)
|
107.898
|
0.013*
|
SG.RS-SG.SO
|
0.305587(0.11343)
|
1(8)
|
0.305587(0.0141787)
|
21.5525
|
0.003**
|
SG.RH-SG.RS
|
0.286011(0.154591)
|
1(8)
|
0.286011(0.0193239)
|
14.8009
|
0.005**
|
SG.EC-SG.RH
|
1.28529(0.121018)
|
1(8)
|
1.28529(0.0151273)
|
84.9651
|
0.007**
|
SG.RH-SG.SO
|
0.0135499(0.150908)
|
1(8)
|
0.0135499(0.0188635)
|
0.718314
|
0.672
|
SG.EC-SG.SO
|
1.34579(0.0798571)
|
1(8)
|
1.34579(0.00998214)
|
134.82
|
0.006**
|
Fungi
|
SG.EC-SG.RH-SG.RS-SG.SO
|
6.63178(6.53519)
|
3(16)
|
2.21059(0.408449)
|
5.41216
|
< 0.001**
|
SG.EC-SG.RS
|
0.407402(2.75986)
|
1(8)
|
0.407402(0.344982)
|
1.18094
|
0.309
|
SG.RS-SG.SO
|
3.31283(3.25717)
|
1(8)
|
3.31283(0.407146)
|
8.1367
|
0.006**
|
SG.RH-SG.RS
|
2.17638(3.47909)
|
1(8)
|
2.17638(0.434886)
|
5.00449
|
0.002**
|
SG.EC-SG.RH
|
2.86233(3.27802)
|
1(8)
|
2.86233(0.409753)
|
6.98551
|
0.006**
|
SG.RH-SG.SO
|
0.429087(3.77533)
|
1(8)
|
0.429087(0.471917)
|
0.909244
|
0.541
|
SG.EC-SG.SO
|
4.07552(3.0561)
|
1(8)
|
4.07552(0.382013)
|
10.6685
|
0.008**
|
Note: * Indicates significant difference between groups (P<0.05); * * indicates extremely significant difference between groups (P<0.01). SS stands for total variance, also known as the sum of squares of deviation; df is degrees of freedom; MS is the mean square (difference), i.e. SS / DF; FS is the F test value; P-value is the p value, with a value less than 0.05 indicating significant difference between groups. The values corresponding to the residual items are in brackets. SG.EC, SG.RS, SG.RH, SG.SO refer to roots, root epidermis, rhizosheath, and rhizosphere soil of Stipa grandis, respectively.