Basic Composition and Structure of the Phyllosphere Microbia Community of Rubber Tree Powdery Mildew
Four samples from each sampling site were sequenced by MiSeq high-throughput sequencing. A total of 16 samples from four areas in Hainan, namely Baisha (BS), Danzhou (DZ), Wuzhishan (WZS) and Wanning (WN), were sequenced for phyllosphere microbial DNA of powdery mildew. A total of 567,903 high quality sequence fragments were read after processing the raw data, with the number of sequences in a single sample ranging from 30, 261 to 44, 775, with an average of 35,494, and a total of 529 OTUs detected. The average sequence length was found to be 429.50 bp, with 89.14% of the sequences in the range of 421–440 bp.
The dilution curves showed differences in the diversity of the different sample species. The flattening of the dilution curves for each sample indicates that the amount of sequenced data is reasonable and the depth of sequencing has covered all species in the samples (Fig. 2).
The community composition of each sample was counted at different taxonomic levels and all samples were identified to a total of 24 phyla, 49 class, 100 orders, 172 families and 256 genus of bacteria. With the exception of WZS_3,15, other 15 of the samples from the four regions were dominated at the phylum level by Cyanobacteria, Poteobacteria and Actinobacteria. The dominant phylum in individual regions also includes Firmicutes and Bacteroidetes, which have different relative abundances in different regions (Figs. 3 and 4).
Except for WZS_3, in the 15 sampls the relative abundance of Cyanobacteria ranged from 71.24–98.29%; Proteobacteria from 1.55–23.13%; and Actinobacteria from 0.82–1.23% in the remaining 15 samples. The relative abundance of WZS_3 was 0.04%, 28.92% and 10.59%, respectively.
Except for WZS_3, at the genus level, the dominant genus in Baisha (BS) were Cyanobacteria, Mitochondria and Rhodococcus; in Danzhou (DZ), Cyanobacteria, Mitochondria and Curtobacterium; in Wanning (WN), Cyanobacteria, Pseudomonas and Pantoea; and in Wuzhishan (WZS), Cyanobacteria, Pseudomonas and Pantoea. The dominant genus in Wanning (WN) are Cyanobacteria, Pseudomonas, Pantoea; the dominant genus in Wuzhishan (WZS) sample WZS_3 are Acidobacteria and other unknown bacteria, the dominant genus in WZS_1 and WZS_2 are Cyanobacteria, Pantoea, and the dominant genus in WZS_CK are Cyanobacteria and Rhodobacteria. Cyanobacteria and Rhodococcus (Fig. 5, Fig. 6)
At the genus level, the most dominant communities at the healthy leave (level 0) of BS, DZ, WN and WZS were Cyanobacteria, with relative abundances of 96.93%, 93.78%, 93.93% and 83.04%, respectively, while the most dominant community at the infected leave (level 3) was also Cyanobacteria, with relative abundances (mean) of 97.34%, 91.2%, 77.64% and 49.88%, respectively (Table 4).
Table 4
Haplotypes information of rubber tree powdery mildew samples
Haplotype | No. of Samples | Herbarium accession No. | Country/Countries |
H1 | 2 | YN-306, HN-509 | China |
H2 | 1 | GD-102 | China |
H3 | 1 | HB1 | Viet Nam |
H4 | 11 | HB2, HO-73, HN-101, HN-205, HN-404, HN-407, MUMH2545, MUMH2602, HKAS 94620, HKAS 94622, HN-501 | Viet Nam, China, Malaysia,Thailand |
H5 | 1 | HKAS 95625 | Sri Lanka |
H6 | 7 | HKAS 94629, HKAS 94633, YN-305, HKAS 94638, HN-102, YN-206, YN-302 | Sri Lanka, China |
H7 | 1 | HN-103 | China |
H8 | 1 | HN-104 | China |
H9 | 1 | HN-105 | China |
H10 | 1 | HN-201 | China |
H11 | 1 | HN-202 | China |
H12 | 1 | HN-203 | China |
H13 | 1 | HN-204 | China |
H14 | 1 | HN-206 | China |
H15 | 1 | HN-207 | China |
H16 | 1 | HN-301 | China |
H17 | 1 | HN-302 | China |
H18 | 1 | HN-303 | China |
H19 | 1 | HN-304 | China |
H20 | 1 | HN-305 | China |
H21 | 1 | HN-306 | China |
H22 | 1 | HN-307 | China |
H23 | 1 | HN-308 | China |
H24 | 1 | HN-402 | China |
H25 | 1 | HN-405 | China |
H26 | 1 | HN-406 | China |
H27 | 1 | HN-408 | China |
H28 | 1 | HN-502 | China |
H29 | 1 | HN-507 | China |
H30 | 1 | HO-373 | China |
H31 | 1 | HO-473 | China |
H32 | 2 | MUMH2418, MUMH2419 | Brazil |
H33 | 1 | YN-101 | China |
H34 | 1 | YN-102 | China |
H35 | 1 | YN-103 | China |
H36 | 1 | YN-105 | China |
H37 | 1 | YN-201 | China |
H38 | 1 | YN-203 | China |
H39 | 1 | YN-204 | China |
H40 | 1 | YN-303 | China |
For the subdominant genus, Mitochondria and Rhodococcus were reported to be found on both uninfected and Level 3 infected leaves of BS. The relative abundance on uninfected leaves was 2.51% and 0.2%, while it was 0.58% and 0.72% on Level 3 infected leaves. For samples of DZ, the relative abundance of uninfected leaves was 2.51% and 0.2%, while the relative abundance of Level 3 affected leaves (at average) was 0.58% and 0.72%. For samples collected from DZ, the subdominant genus on healthy leaves was Mitochondria with the relative abundance at 6.03%, while the relative abundance of Level 3 infected leaves was 1.41% and 2.18% with Mitochondria and Curtobacterium as subdominant genus. For samples from WN, the subdominant genus were Mitochondria and Bacillus (Bacillus spp.) on healthy leaves with relative abundances of 3.59% and 1.87% respectively, while Mitochondria, Pseudomonas, Pantoea and Curtobacterium were found in Level 3 affected leaves with relative abundances (average) of 3.36%, 9.39%, 6.72% and 1.45%. Mitochondria and Rhodococcus were found as subdominant on healthy leaves from WZS with relative abundances of 1.05% and 6.63%, while Mitochondria, Pantoea, Acidobacteria and Enterobacteriaceae were found in Level 3. The WZS healthy leaves had Mitochondria and Rhodococcus, with relative abundances of 1.05% and 6.63%, while the Level 3 affected leaves had Mitochondria, Pantoea, Acidobacteria and Enterobacteriaceae as the subdominant genus, with relative abundances (at average) of 2.59%, 6.25%, 4.47% and 3.51% respectively (Table 5).
Table 5
The relative abundance of bacterial community in rubber phyllosphere at genus level
Genus | Community relative abundance of BS (%) | Community relative abundance of DZ (%) | Community relative abundance of WN (%) | Community relative abundance of WZS (%) |
Level 3(Average) | Level 0 | Level 3(Average) | Level 0 | Level 3(Average) | Level 0 | Level 3(Average) | Level 0 |
Cyanobacteria_norank | 97.34 ± 1.32 | 96.93 | 91.2 ± 3.12 | 93.78 | 77.64 ± 9.10 | 93.93 | 49.88 ± 43.18 | 83.04 |
Pantoea | 0.05 ± 0.04 | 0.00 | 0.84 ± 1.27 | 0.02 | 6.72 ± 7.51 | 0.02 | 6.25 ± 5.44 | 0.36 |
Mitochondria_norank | 0.58 ± 0.36 | 2.51 | 1.41 ± 0.03 | 6.03 | 3.36 ± 2.14 | 3.59 | 2.59 ± 2.59 | 1.05 |
Pseudomonas | / | / | 0.26 ± 0.03 | 0.02 | 9.39 ± 7.87 | 0.01 | 1.38 ± 2.0 | 1.16 |
Acidobacteria_norank | / | / | / | / | / | / | 4.47 ± 7.7 | 0.13 |
Curtobacterium | / | / | 2.18 ± 1.29 | 0.00 | 1.45 ± 0.58 | 0.16 | 0.64 ± 0.75 | 0.21 |
Rhodococcus | 0.72 ± 0.48 | 0.20 | 0.44 ± 0.08 | 0.01 | 0.23 ± 0.1 | 0.10 | 0.29 ± 0.31 | 6.63 |
Enterobacteriaceae_unclassified | / | / | / | / | / | / | 3.51 ± 5.22 | 0.01 |
Ralstonia | 0.34 ± 0.21 | 0.13 | 0.92 ± 0.05 | 0.03 | 0.17 + 0.06 | 0.14 | 0.58 ± 0.59 | 0.43 |
Sphingomonas | / | / | 0.85 ± 0.72 | 0.04 | / | / | 0.93 ± 1.1 | 0.15 |
Bacillus | / | / | / | / | 0.08 ± 0.1 | 1.87 | 0.65 ± 0.97 | 0.00 |
Others | 0.89 ± 0.56 | 0.20 | 1.87 ± 0.57 | 0.07 | 0.76 ± 0.09 | 0.19 | 9.41 ± 12.92 | 6.69 |
Note: Bacterial communities of < 0.02% were not counted. |
In summary, except for the WZS samples, the dominant genus of healthy leaves in all three regions were Cyanobacteria and Mitochondria, while the composition of the dominant genus on in Level 3 infected leaves did not vary much, as they all included Cyanobacteria and Mitochondria. Other dominant genus were different in different regions, but their relative abundance was not high. Therefore, the composition of the genus-level communities of the three levels of infected and healthy leaves in the same region was similar, except for differences in relative abundance, while the composition of the genus-level communities differed significantly between regions.
A total of 520 OTUs were detected in the four regions after quality control, 102 OTUs in BS, 113 OTUs in DZ, 81 OTUs in WN and 479 OTUs in WZS (Fig. 7). The proportion of OTUs unique to each of the four regions was 3.27% for BS, 1.54% for DZ, 0.77% for WN and 67.69% for WZS, showing that WZS > BS > DZ > WN. There were 36 overlapping OTUs in the four regions, accounting for 6.92% of the total number of samples. This indicates that although the samples were all collected from rubber powdery mildew phyllosphere environment, there were significant differences in OTU levels between the four regions.
Diversity Index Analysis Of Phyllosphere Microbia Of Rubber Tree Powdery Mildew
The abundance and diversity of the microbial community can be studied by analysing the diversity of individual samples (Alpha diversity) using a range of statistical indices such as Chao1, Ace, Shannon Index and Simpson Index diversity (Table 6).
Table 6
Diversity indexes for bacterial community in rubber phyllosphere
Sample ID | Reads | 97% identity |
OTU | ACE | Chao 1 | Coverage | Shannon | Simpson |
BS_1 | 12214 | 36 | 49 | 54 | 99.89% | 0.14 | 0.96 |
BS_2 | 12214 | 73 | 163 | 126 | 99.73% | 0.31 | 0.91 |
BS_3 | 12214 | 33 | 144 | 76 | 99.84% | 0.16 | 0.96 |
BS_CK | 12214 | 24 | 86 | 47 | 99.89% | 0.17 | 0.94 |
DZ_1 | 12214 | 65 | 110 | 138 | 99.75% | 0.59 | 0.81 |
DZ_2 | 12214 | 60 | 82 | 95 | 99.83% | 0.68 | 0.78 |
DZ_3 | 12214 | 60 | 83 | 76 | 99.83% | 0.4 | 0.88 |
DZ_CK | 12214 | 18 | 34 | 63 | 99.92% | 0.26 | 0.88 |
WN_1 | 12214 | 41 | 85 | 61 | 99.87% | 0.9 | 0.58 |
WN_2 | 12214 | 42 | 133 | 105 | 99.84% | 0.62 | 0.78 |
WN_3 | 12214 | 50 | 170 | 89 | 99.80% | 1.04 | 0.53 |
WN_CK | 12214 | 25 | 48 | 29 | 99.93% | 0.32 | 0.88 |
WZS_1 | 12214 | 104 | 136 | 130 | 99.74% | 1.28 | 0.55 |
WZS_2 | 12214 | 78 | 282 | 146 | 99.66% | 1.15 | 0.57 |
WZS_3 | 12214 | 324 | 329 | 329 | 99.89% | 5.09 | 0.01 |
WZS_CK | 12214 | 128 | 132 | 135 | 99.90% | 1.02 | 0.69 |
Note: BS is Baisha, DZ is Danzhou, WN is Wanning, WZS is Wuzhishan; fungal communities which are less than 0.02 % are not included in the statistics.
The OTU statistics were analysed at a 97% similarity level and the OTU numbers for the 16 samples ranged from 24–324, with the Ace and Chao indices used to estimate the total number of species at 34–329 and 29–329 respectively. Shannon Index and Simpson Index were used to estimate microbial diversity in the samples, ranging from 0.14–5.09 and 0.01–0.96 respectively.
The coverage of the samples was greater than 99%, proving that the sequencing results represent the true picture of the microorganisms in the samples. Despite distinct algorithms, the Ace and Chao Index were generally consistent and did not differ significantly from the measured OTU numbers.
Both the Shannon and Simpson Index indicated that WZS > WN > DZ > BS in terms of the alpha diversity of species in the four regions.
Pca Analysis And Nmds Analysis Of Phyllosphere Microbia In Samples From Different Regions
The study compares differences in the composition of phyllosphere bacterial communities of rubber trees from four regions in Hainan. The composition of the DZ and BS samples was similar, and these two areas were close together in the PCA diagram due to their geographical proximity to each other. The samples from WN and WZS were clearly distinguishable from each other, except for the WZS_3 sample, which was more dispersed but still mostly clustered together (Fig. 8). The results indicate that bacterial communities differed significantly between the four regions and that different geographical locations are influential in accounting for the variation of community composition and structure of phyllosphere bacteria in rubber trees.
Non-metric multidimensional scaling showed that, except for the WZS_3 samples, samples that had close geographical proximity or similar climates shared more similarities in composition, as analysis results were clustered together (Fig. 4). For example, DZ and BS has a short distance between one another with a dry climate, whereas WN and WZS are close with a wet climate. These results suggest that the different geographical locations or climate types can lead to the separation of different microbial communities, and that samples with geographical proximity or similar climate types illustrate large clustering results.
Analysis Of Phyllosphere Microbia Taxa Lefse In Samples From Different Regions
A Kruskal-Wallis (KW) sum-rank test was performed using the LEfSe ((Linear discriminant analysis effect size) algorithm to compare the distribution of OTUs between the leaves of rubber trees in different regions. The results showed that five OTUs were distributed in BS, three in DZ, five in WN and 24 in WZS. Cyanobacteria were the most abundant taxa in the BS, Oxalobacteraceae in the DZ, Pseudomonadales in the WN, and Proteobacteria and Acidobacteria in the WZS (Fig. 9).
The results of the Plot cladogram analysis also reveal the communities or species that have a significant differential impact on division of sample by region (Fig. 10).
The communities with significant differential impact are Cyanobacteria in BS, Methylobacterium in DZ, are Firmicutes in WN and Proteobacteria and Acidobacteria in WZS.
Correlation Analysis Between Clustering Of Phyllosphere Microbia Otus And Environmental Factors
The correlations between the environmental factors (temperature, humidity, altitude) sampled at each location and the OTUs were calculated. The achieved R and P values were obtained to cluster the OTUs and environmental factors separately (Fig. 11). The results confirmed that OTU280, OTU257, OTU146, OTU27, OTU437, OTU196, OTU385, OTU30 and OTU510 were significantly correlated with temperature and clustered into one large branch, OTU430 was significantly correlated with temperature and humidity and was in this large branch; OTU144 was significantly correlated with temperature and humidity but was a separate branch. Altitude has no effect on OTUs. OTU144 is highly correlated with temperature and humidity, OTU196 and OTU30 are highly correlated with temperature and OTU 430 is moderately correlated with humidity.
Basic Composition And Structure Of The Phyllosphere Fungal Community Of Rubber Tree Powdery Mildew
A total of 593,237 high quality sequence fragments were obtained from 16 samples, with an average of 37, 077 sequences in the range of 30, 129–44, 619 and 1,406 OTUs detected in a single sample. 293.01 bp was the average sequence length, with 57.38% of the sequences concentrated between 301–320 bp and 26.99% between 241–260 bp.
The dilution curves show the differences in the diversity of the different sample species. The flattening of the dilution curves for each sample indicates that the amount of sequenced data is reasonable and the depth of sequencing has largely covered all species in the samples (Fig. 13).
The community composition of each sample was calculated at different taxonomic levels and a total of 6 phyla, 26 orders, 69 families, 157 families and 351 genus of fungi were identified in all samples. At the phylum level, the community composition of the samples from all four regions was similar, with Ascomycota, Fungi_unclassified and Basidiomycota dominating, with different relative abundances in the samples from different regions (Fig. 14). The relative abundance of Ascomycota ranged from 43.76–99.80%, Fungi_unclassified from 0.02–53.32% and Basidiomycota from 0.17–10.91% in the 16 samples.
Except for WZS_3, the community composition of the samples was genuslly consistent at the generic level, but varied in abundance, while the community composition and relative abundance differed significantly between regions (Fig. 15). The dominant genus in all four regions include Erysiphe, Cladosporium and Fungi_unclassified. In addition, the different regions have varied dominant genus. The dominant genus in the Baisha region (BS) are Alternaria, in the Danzhou region (DZ) are Sporidiobolales, in the Wanning region (WN) are Dothioraceae, and in the Wuzhishan region (WZS) are Uwebraunia, as shown in Fig. 16.
For infected leaves at Level 3, the dominant genus of BS were Erysiphe, Cladosporium, Fungi_unclassified and Alternaria with relative abundances (mean) of 71.81%, 17.75%, 1.19% and 1.8%, respectively. The dominant genus were Erysiphe, Cladosporium and Sporidiobolales in DZ.
In WN, the dominant genus were Erysiphe and Cladosporium, with relative abundances (mean) of 93.24% and 1.16%, respectively. In WZS, the dominant genus were Erysiphe and Cladosporium, Fungi_unclassified and Ascomycota, and the relative abundances (mean) were 63.97%, 1.66%, 14.99% and 5.66%, respectively (Table 7).
In addition, healthy leaves (Level 0) and infected leaves at Level 3 in the four regions also included a number of less abundant communities such as Epicoccum, Mycosphaerellaceae, Pseudocercospora, Pleosporales and Pseudozym, whose relative abundances ranged from 63.97%, 1.66%, 14.99% and 5.66%, respectively (Table 7). The relative abundances of Pseudocercospora, Pleosporales and Pseudozym range from 0.12–0.98%, 0.02–1.88%, 0.02–1.44%, 0.04–1.75% and 0.09–2.75%, respectively. The relative abundance of fungi (others) whose taxonomic status is unclear or not yet identified ranged from 1.04–10.58% (Table 7).
In summary, except for samples from WZS, the dominant genus in the remaining three regions for both Level 3 infected leaves and healthy leaves included Erysiphe, Cladosporium and Fungi_unclassified, with 1–3 different subdominant genus in different regions, but all with low relative abundance. Therefore, the genus-level community composition of the three levels of infected and healthy leaves in the same area did not vary much, but only in relative abundance, whereas the genus-level community composition varied between areas. Powdery mildew had little effect on the structure of the phyllosphere fungal community, but had an effect on the relative abundance of each community.
A total of 1398 OTUs were detected in the four regions after quality control, 420 OTUs in BS, 295 OTUs in DZ, 500 OTUs in WN and 875 OTUs in WZS (Fig. 17). the proportion of unique OTUs in the four regions was 11.09% in BS, 1.93% in DZ, 13.23% in WN and 44.49%, showing WZS > WN > BS > DZ. There were 90 overlapping OTUs in the four regions, accounting for 6.44% of the total number of OTUs in the samples. This indicates that although the samples were all sourced from rubber powdery mildew phyllosphere, there were significant differences in OTU levels between the four regions.
Diversity Index Analysis Of The Phyllosphere Fungal Flora Of Rubber Tree Powdery Mildew
The abundance and diversity of the microbial community can be studied by analysing the diversity of individual samples (Alpha diversity) using a range of statistical indices such as Chao1, Ace, Shannon index and Simpson index diversity (Table 8). The OTU statistics were analysed at 97% similarity level and the OTU numbers for the 16 samples ranged from 71–692, with Ace and Chao Index used to estimate the total number of species, 121–723 and 120–713, respectively. Both the Shannon Index and Simpson Index indicated that the alpha diversity of species in the four regions was WZS > BS > DZ > WN.
Table 8
Diversity indexes for fungal community in rubber phyllosphere
Sample ID | Reads | 97% identity |
OTU | ACE | Chao 1 | Coverage | Shannon | Simpson |
BS_1 | 29913 | 105 | 121 | 120 | 99.92% | 0.53 | 0.82 |
BS_2 | 29913 | 193 | 204 | 204 | 99.92% | 1.84 | 0.31 |
BS_3 | 29913 | 171 | 236 | 217 | 99.80% | 0.79 | 0.73 |
BS_CK | 29913 | 240 | 251 | 251 | 99.92% | 2.65 | 0.15 |
DZ_1 | 29913 | 132 | 191 | 183 | 99.84% | 1.53 | 0.34 |
DZ_2 | 29913 | 124 | 179 | 163 | 99.86% | 1.61 | 0.28 |
DZ_3 | 29913 | 159 | 336 | 269 | 99.77% | 0.93 | 0.65 |
DZ_CK | 29913 | 147 | 215 | 204 | 99.81% | 1.09 | 0.44 |
WN_1 | 29913 | 196 | 296 | 321 | 99.73% | 0.47 | 0.88 |
WN_2 | 29913 | 243 | 316 | 291 | 99.76% | 0.67 | 0.83 |
WN_3 | 29913 | 186 | 266 | 267 | 99.77% | 0.52 | 0.86 |
WN_CK | 29913 | 381 | 421 | 415 | 99.79% | 2.82 | 0.17 |
WZS_1 | 29913 | 72 | 169 | 115 | 99.89% | 0.33 | 0.89 |
WZS_2 | 29913 | 71 | 230 | 142 | 99.87% | 0.24 | 0.92 |
WZS_3 | 29913 | 692 | 723 | 713 | 99.78% | 4.64 | 0.02 |
WZS_CK | 29913 | 182 | 225 | 215 | 99.84% | 1.79 | 0.30 |
Lefse Analysis Of Phyllosphere Fungal Taxa In Samples From Different Regions
A Kruskal-Wallis (KW) sum-rank test was performed using the LEfSe ((Linear discriminant analysis effect size) algorithm to compare the distribution of OTUs between the leaves of rubber trees in different regions. The results showed that 11 OTUs were distributed in BS, 3 OTUs in DZ, 22 OTUs in WN and 3 OTUs in WZS. Trichosphaeriales is the most abundant taxon in the BS, Sporidiobolaceae is the most abundant taxon in the DZ, Dothideales and Dothioraceae are the most abundant taxa in the WN, and Chaetomiaceae is the most abundant taxon in the WZS (Fig. 18). The results of the Plot cladogram analysis also revealed the communities or species that significantly influenced the delimitation of the samples in different regions (Fig. 19).
The communities significantly differentially affected by BS are Trichosphaeriales, Botryosphaeriaceae, Capnodiales, Montagnulaceae, Pleosporales, Sclerotiniaceae and Cystofilobasidiaceae; the communities significantly differentially affected by DZ are Botryosphaeriales, Dothideomycetes and Cystofilobasidiaceae. DZ The significantly different communities are Botryosphaeriales, Dothideomycetes and Sporidiobolaceae; WZS The significantly different communities are Chaetomiaceae, Sordariales and Psathyrellaceae; The communities affected by significant differences in WN were Teratosphaeriaceae, Dothioraceae, Dothideomycetes, Diaporthales, Cyphellaceae, Schizophyllaceae, Schizophyllaceae, and WN. Schizophyllaceae, Auriculariales, Phanerochaetaceae, Polyporales, Peniophoraceae, Russulales, Basidiomycota, and Ustilaginomycetes.
Hierarchical Clustering Tree Of Phyllosphere Fungal Samples From Different Regions
Sample Hierarchical clustering tree was constructed by Hierarchical cluatering analysis based on beta diversity distance matrix, which can compare similarity and difference among multiple samples (Fig. 20). The results showed that the community composition of WN_1, WN_2 and WN_3 was similar to each other and clustered into a single branch. The similarity between WZS_1 and WZS_2 is high, and DZ_1 and DZ_2 are clustered together. The similarity of WZS_CK, WN_CK, BS_CK and DZ_CK community composition was not high, and the similarity between Wzs_CK and the other three samples in this area was low, and the difference was obvious.
Correlation Analysis Of Otus Clustering Of Interleaf Fungi With Environmental Factors
The correlations between the environmental factors (temperature, humidity, altitude) sampled at each location and the OTUs were calculated. The achieved R and P values were obtained to cluster the OTUs and environmental factors separately (Fig. 21). The results showed that OTU1391, OTU202, OTU953, OTU985, OTU1073 and OTU769 were significantly correlated with temperature and humidity and clustered into one group; OTU420, OTU419, OTU173 and OTU856 were significantly correlated with temperature and clustered into one branch; OTU652, OTU758, OTU1289 and OTU12 It was significantly correlated with altitude and clustered into one group, OTU223, OTU76, OTU1066 and OTU96 were significantly correlated with altitude and clustered into one cluster; OTU919 was significantly correlated with temperature. Among them, OTU76 and OTU1066 are highly correlated with altitude factor, OTU1391 is highly correlated with humidity, and OTU420, OTU856 and OTU1073 are highly correlated with temperature.