Structure of Bacterial Communities in Phosphorus-Enriched Rhizosphere Soils

: Although phytoremediation is the main method for P-removal and maintaining ecosystem balance in geological phosphorus-enriched soils (PES), little is known about the structure and function of microbial communities in PES. Interactions between plants and soil microorganisms mainly occur in the rhizosphere. The aim of this work was to investigate the composition and diversity of bacterial communities found in rhizosphere soils associated with the following three dominant plant species: Erianthus ruﬁpilus , Coriaria nepalensis , and Pinus yunnanensis . In addition, we compared these rhizosphere bacterial communities with those derived from bulk soils and grassland plots in PES from the Dianchi Lake basin of southwestern China. The Illumina MiSeq platform for high-throughput sequencing of 16S rRNA was used for the taxonomy and the analysis of soil bacterial communities. The results showed higher bacterial diversity and nutrient content in rhizosphere soils as compared with bulk soils. Rhizosphere bacteria were predominantly comprised of Proteobacteria (24.43%) and Acidobacteria (21.09%), followed by Verrucomicrobia (19.48%) and Planctomycetes (9.20%). A comparison of rhizosphere soils of the selected plant species in our study and the grassland plots showed that Acidobacteria were the most abundant in the rhizosphere soil of E. ruﬁpilus ; Bradyrhizobiaceae and Rhizobiaceae in the order Rhizobiales from C. nepalensis were found to have the greatest abundance; and Verrucomicrobia and Planctomycetes were in higher abundance in P. yunnanensis rhizosphere soils and in grassland plots. A redundancy analysis revealed that bacterial abundance and diversity were mainly inﬂuenced by soil water content, soil organic matter, and total nitrogen. and diversity in PES are mainly inﬂuenced by soil water content, soil organic matter, and total nitrogen. Our results indicate that rhizosphere microbial communities could serve as


Introduction
Soil microbial communities are responsible for many ecosystem functions, including decomposition, nutrient cycling, nitrogen fixation, and soil formation [1][2][3]. Microorganisms that inhabit different soil types exhibit astonishing diversity and wide genetic variability even within species, particularly with respect to their metabolic pathways and host-interactive capabilities [4]. The majority of microbial populations in soils are concentrated in nutrient-rich niches such as rhizospheres that offer a constant supply of easily utilizable nutrients [5,6]. Recent advances in soil ecosystem microbial ecology have highlighted the interactions that occur among distinct microbial groups, now known as the "microbiome" [7].
Soil microbial communities are influenced by a series of biotic and abiotic factors [8], and specific microbial groups show high levels of diversity in different soil ecosystems [9,10]. Exogenous nutritional inputs can change structures of soil bacterial communities, promoting plant phosphorus uptake in rhizospheres [11,12]. Moreover, plant species' richness and diversity affect phosphorus uptake in rhizospheres through interactions between plant and soil bacterial communities [13]. Rhizospheric effects influence bacterial communities in rhizospheres, and therefore the diversity of rhizosphere microbial communities differs from that of the bulk soil [14].
Several studies have reported that phosphorus is a limiting factor responsible for shaping structures of soil bacterial communities across different soil habitats [15,16]. Phosphorus-rich mountains are also vulnerable and face severe degradation due to mining and other human disturbances [17]. Phosphorus removal from phosphorus-enriched soils (PES), such as in agricultural lands, grasslands, and phosphate mining regions is a major cause of eutrophication in aquatic environments and an increasing environmental problem worldwide [18]. The Dianchi Lake basin is an important industrial and agricultural region in southwest China, home to numerous PES [19,20]. Dianchi Lake is a typical plateau lake in the basin currently undergoing severe eutrophication [21]. Previous studies regarding the restoration of PES have mainly focused on the development of dominant plant communities which is a cost-effective method [17,22].
Rhizosphere microbes play important roles for sustaining ecosystems. For example, they accelerate phytoremediation processes and enhance plant ecosystem productivity [23]. Microbes inhabiting rhizospheres help plants grow and function more effectively through bolstering water retention, disease resistance, and nutrient acquisition [24]. Many studies have also shown that rhizosphere microbes promote soil remediation and restoration. Crops have been observed to reduce salt content and increase bacterial rhizosphere diversity in coastal saline soils, improving soil quality [25]. In addition, the inoculation of plant growth-promoting bacteria can change rhizosphere bacterial communities and accelerate restoration processes in eroded desert soils [26]. Rhizosphere bacteria foster arsenic bioaccumulation by plants and then improve phytoremediation of arsenic-contaminated soils [27]. However, little is known about microbes in PES. Previous studies have investigated functional bacteria taxa in PES and have identified 377 isolates that exhibited phosphate-solubilizing potential [28], although biodiversity and distribution of bacteria still remain unclear. Understanding bacterial diversity and distribution is the first step towards harnessing functional taxa to improve phytoremediation and recalibrate restoration efforts. This work highlights the urgency of studying microbial community structures in rhizospheres to support the restoration of PES.
The main objectives of this study were to investigate divergences in structures of soil bacterial communities between rhizosphere and bulk soils of plants in PES as well as understand the environmental factors that affect bacterial communities. In previous studies on phosphorus-rich mountains, efforts have centered around understanding phytoremediation. This work reports for the first time the structure of PES microbial communities influenced by dominant plants of PES.

Study Area
The sampling was conducted on phosphorus-rich mountains near Dianchi Lake in central Yunnan, China ( Figure 1). In this study area, the average annual temperature is 14.7 • C, and the annual temperature difference is 10-15 • C. The average annual rainfall is 985.5 mm, and the rainy season lasts from May to October, accounting for 84-90% of the annual rainfall. In this study area, the variation range of the total phosphorus concentration in the soil and the nitrogen/phosphorus ratio are 1.15-80.2 g/kg and 0.006-0.98, respectively, [29] both of which are typical in a phosphorus-rich region.

Study Area
The sampling was conducted on phosphorus-rich mountains near Dianchi Lake in central Yunnan, China ( Figure 1). In this study area, the average annual temperature is 14.7 °C, and the annual temperature difference is 10-15 °C. The average annual rainfall is 985.5 mm, and the rainy season lasts from May to October, accounting for 84-90% of the annual rainfall. In this study area, the variation range of the total phosphorus concentration in the soil and the nitrogen/phosphorus ratio are 1.15-80.2 g/kg and 0.006-0.98, respectively, [29] both of which are typical in a phosphorusrich region. Most vegetation in the research area comprises natural secondary forests, and the original subtropical semi-moist evergreen broad-leaved forest has been almost completely destroyed [30]  Most vegetation in the research area comprises natural secondary forests, and the original subtropical semi-moist evergreen broad-leaved forest has been almost completely destroyed [30]. The dominant vegetation species identified for the study were Pinus yunnanensis, Coriaria nepalensis, Eupatorium adenophora, Erianthus rufipilus, Rumex hastatus, Keteleeria evelyniana, Pinus armandii, and Eucalyptus urophylla. This study selected three dominant plants for research, i.e., E. rufipilus, C. nepalensis, and P. yunnanensis, based on the investigation of previous research. E. rufipilus belongs to the family Gramineae, is a drought-tolerant herbaceous plant growing between elevation 1300-2400 m, has a plant density of 2610 plant/ha and a P removal potential of 22.78 kg/ha in PES. C. nepalensis belongs to the family Coriariaceae, is an actinorhizal N-fixing species growing between elevation 400-3200 m, has a plant density of 670 plant/ha and a P removal potential of 16.91 kg/ha in PES. P. yunnanensis belongs to the family Pinaceae growing between elevation 400-3500 m, has a plant density Appl. Sci. 2020, 10, 6387 4 of 18 of 497 plant/ha, a P removal potential of 27.69 kg/ha in PES, and features secondary communities that suffer from disturbance and invasion [17,22].

Soil Sampling
Fieldwork was conducted in March 2016. This study was conducted during the dry season, when rhizosphere soil acquisition is feasible, while also limiting the occurrence of rain erosion on bacterial communities. A total of thirteen sites were selected for the study area. Sampling sites were selected according to target plant communities. The soil samples in thirteen sites (Table 1) included the following: (1) rhizosphere and bulk soil of E. rufipilus from three sites (Er-R and Er-B), (2) rhizosphere and bulk soil of C. nepalensis from three sites (Cn-R and Cn-B), (3) rhizosphere and bulk soil of P. yunnanensis from three sites (Py-R and Py-B), and (4) soils from grassland plots from four sites (CK) with removal of surface grass debris at a random 0-20 cm soil depth in the absence of vegetation coverage. We selected plots of 10,000 m 2 (100 × 100 m) for each site. Five subplots (5 × 5 m, with at least one target plant in each plot) were chosen for soil sampling using an S-shaped sampling method [31] within sampling sites. Rhizosphere soils were collected from the fine roots of E. rufipilus, C. nepalensis, and P. yunnanensis. Bulk soil samples were taken at a depth of 0-20 cm in the same plots as rhizosphere soils, and then mixed thoroughly. Rhizosphere soils were collected by gently shaking off the bulk soil that adhered to the roots at a depth of 0-20 cm, and then mixed well. Each soil sample was sieved through a 2 mm mesh to remove plant roots and other plant materials. Sampling equipment including mesh, scalpels, and spoons was cleaned with 95% alcohol, followed by heating using an alcohol lamp for 30 s. Soil collections among samples were conducted after all sampling equipment was sterilized and cooled down. All soil samples were tagged and sealed in valve bags from fine roots at 0-20 cm depth. At the same time, the soil temperature (Tem) of each plot was measured using soil temperature instruments for five replicates. Next, all 22 soil samples were manually homogenized and divided into three parts. The first portion was preserved within sterile centrifuge tubes (2 mL) and put in liquid nitrogen for DNA extraction and molecular analysis; the second portion was kept in plastic bags for water content measurement; and the third portion was stored at 4 • C refrigerator for 24 h, and then dried at 25 • C for chemical analysis.

Soil Physicochemical Analysis
Chemical properties for all the 22 soil samples from the study area were analyzed based on the following method: Soil water content (SWC) was determined using the oven-drying method [32]. The pH value (pH) was determined in a 1:2 (soil to water ratio) using a pH meter with a standard combination of electrodes [33]. Soil organic matter content (SOM) was determined through the potassium dichromate sulfuric acid oxidation method [34]. Total nitrogen (TN) was measured using Kjeldahl's method as modified by Bremner and Mulvaney [35]. Alkali solution nitrogen (AN) was determined using the alkali soluble diffusion method [36]. Total phosphorus (TP) content was measured using the ammonium molybdate spectrophotometric method [37]. Concentrations of available phosphorus (AP) were measured using the ammonium molybdate spectrophotometric method, following hydrochloric and sulfuric acid leaching. The value of nitrogen/phosphorus (N/P) was calculated by the ratio of total nitrogen to total phosphorus. All 22 soil samples from the four parts were used in chemical tests. Lab analyses were repeated three times, and final data were averaged.

Illumina MiSeq for High-Throughput Sequencing of 16S rRNA
Soil samples stored in liquid nitrogen were used for DNA extraction. An Agilent High Sensitivity DNA Kit (Agilent Technologies Inc., Lexington, MA, USA) was used with 0.5 g of soil, according to the manufacturer's instructions. Extraction protocol followed the manufacturer's instructions with recommended modifications to enhance the efficiency of cell lysis. To assess the quality and purity of DNA, crude DNA extracts (2 µL) of each sample were run on 2% agarose gel and analyzed using ultraviolet spectrophotometer (Eppendorf Corporation, Hamburg, Germany). Detection results of qualified DNA samples are shown in the Appendix A Table A1. Electrophoresis detection bands were single, indicating no degradation, and had no protein or RNA contamination.
Bacterial 16S rRNA genes were amplified using the following two primers: 520F (5-AYTGG GYDTAAAGNG-3) and 802R (5-TACNVGGGTATCTAATCC-3). PCR reactions were carried out in triplicate with 25 µL of reaction mixture, comprising 8.75 µL of sterilized ultrapure water, 5 µL of Q5 reaction buffer, 5 µL of Q5 GC high enhancer, 2 µL of 2.5 mM dNTPs, 2 µL of DNA template, 1 µL of each primer, and 0.25 µL Q5 polymerase. For the subsequent PCR cycling reaction, the following parameters were used: 98 • C and 2 min for initial denaturation, followed by 25 cycles of denaturation at 98 • C for 30 s, annealing at 50 • C for 30 s, and elongation at 72 • C for 30 s, followed by a final extension at 72 • C for 5 min. All the samples were amplified in triplicate and detected and quantified using a Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen Corporation, Waltham, MA, USA). This step ensured that all detected DNA samples had obvious and single bands before the next experiments. After quantification, paired-end sequencing of amplicons from each reaction mixture was performed using the Illumina MiSeq platform [38]. Sequencing services were provided by Personal Biotechnology Co., Ltd. Shanghai, China. The raw sequencing data of bacteria were deposited in NCBI Sequence Read Archive (SRA) under accession number PRJNA662457.

Data Analysis
At least three groups of repeats and 22 soil samples in total from the geological phosphorus-rich mountains were collected. There were seven sets of mixed samples (Er-R, Er-B, Cn-R, Cn-B, Py-R, Py-B and CK) (as seen in Table 1) for the comparative analysis of bacterial communities. The resulting data from each set of repeated samples were averaged. For the FASTQ format of paired-end sequences, we ensured the mean quality of the base was ≥Q20 and sequence length ≥150 bp. FLASH was used to pair and connect the screened paired-end sequences in order to obtain high quality sequences [39]. UCLUST [40] was used to align sequences and call OTUs (operational taxonomic units) using 97% as the cut-off sequence similarity, showing the rarefaction curves, then removing OTUs in which the abundance value was lower than 0.001% of the total sequenced quantity [41]. Subsequent microbial classification and analyses were conducted based on classifying the information of each OTU. For the richness and diversity index of bacterial communities, the Chao1 index (1) and ACE (2) index were used to describe the richness [42]: Appl. Sci. 2020, 10, 6387 where S obs is the total number of species observed in a sample and f k is the number of species each represented by k individuals in a single soil sample. The Shannon index (3) [43] and Simpson index (4) [44] were used to describe the evenness: where P i = N i /N, N i is the individual number of species (i), N is the total number of organisms of a particular species, and P i is the proportion of i species.
Raw physicochemical data properties of soil samples were processed with Excel 2007, and the results were calculated for three (rhizosphere and bulk soil of three dominant plants) or four replicate (CK) samples and expressed by the mean ± standard deviation. Statistical significance was determined by one-way analysis of variance (ANOVA), followed by Dunnett's T3 post hoc test in SPSS ver. 19.0. To compare the structure of bacterial communities across all soil samples based on the relative abundance of bacterial phyla class, principal components analysis of soil properties using redundancy analysis (RDA) by R software package was conducted. Bacterial community composition data, for each classification level, were clustered according to their number of OTUs and taxonomic data using MEGAN [45], and a phylogenetic tree based on the sequenced result of each sample and NCBI database was constructed. Heatmap and Venn diagram were drawn by hierarchal clustering performed based on Ihaka and Gentleman [46]. GraPhlAn [47] was used to construct circular representations of the taxonomic phylogenetic tree, showing the dominant taxa of rhizosphere soils and soils from grassland plots.

Physicochemical Properties
There were significant differences among the rhizosphere soils of E. rufipilus (Er-R), C. nepalensis (Cn-R), P. yunnanensis (Py-R), and grassland plots (CK) ( Table 1) with respect to most physicochemical properties ( Table 2). The amounts of SOM, AN, TP, and AP were significantly higher in Py-R as compared with Er-R, Cn-R, and CK. Gathered data demonstrate that rhizosphere soils of P. yunnanensis possess a rich nutrient content

Sequencing Results
Analyses revealed an average of 40,667 valid sequences per sample, and subsequent screening analyzed 38,072 high quality sequences (accounting for 89.3-98.05%), with a read length of 200-250 bp.

Differences in Bacterial Diversity
The rarefaction curves ( Figure 2) show that Py-R had the lowest bacterial diversity, and CK were more diverse as compared with the rhizosphere soils of Er-R and Cn-R. For richness and diversity estimations based on OTUs picker data, four indices of each sample group were analyzed using QIIME ( Table 3). As shown in Table 3, scores for indices were as follows: Chao1 (1595) and ACE (2036) for Py-R; Chao1 (2012) and ACE (2511) for CK; the Simpson index was 0.981 and the Shannon index 8.11 for Py-R; and the Simpson index was 0.992 and the Shannon index 8.93 for Er-R. The data show that CK had the highest richness and Er-R had the highest diversity, while Py-R had the lowest richness and diversity.

Sequencing Results
Analyses revealed an average of 40,667 valid sequences per sample, and subsequent screening analyzed 38,072 high quality sequences (accounting for 89.3-98.05%), with a read length of 200-250 bp.

Differences in Bacterial Diversity
The rarefaction curves ( Figure 2) show that Py-R had the lowest bacterial diversity, and CK were more diverse as compared with the rhizosphere soils of Er-R and Cn-R. For richness and diversity estimations based on OTUs picker data, four indices of each sample group were analyzed using QIIME ( Table 3). As shown in Table 3, scores for indices were as follows: Chao1 (1595) and ACE (2036) for Py-R; Chao1 (2012) and ACE (2511) for CK; the Simpson index was 0.981 and the Shannon index 8.11 for Py-R; and the Simpson index was 0.992 and the Shannon index 8.93 for Er-R. The data show that CK had the highest richness and Er-R had the highest diversity, while Py-R had the lowest richness and diversity.

Bacterial Communities and Diversity
Relative bacterial OTU richness was present in four different soil samples across phylum, class, order, family, and genus ( Figure 3). The dataset for phylogenetic analyses comprised 1 kingdom, 3 phyla, 7 classes, 7 orders, 16 families and 40 genera. The heatmap plot depicted the relative abundance of each dominant genera (variables clustering on the Y-axis) within soil samples (X-axis clustering) ( Figure 4). As shown in Figure 4, four samples from CK grouped together; Cn-R2, Cn-R1, Er-R3, Er-R2 grouped together; and Py-R2 and Py-R3 grouped together. The 50 most abundant genera were depicted. The results show that the majority of sequences belonging to Asteroleplasma, Enterococcus, Opitutus, and Alkaliphilus were presented in soil from the CK group. Dokdonella, Pseudomonas, Snaeromyxobacter, and Hyphomicrobium were presented in soil from the Er-R and Cn-R groups. Candidatus Xiphinematobacter, Mycobacterium, Bradyrhizobium, and Nostocoida were presented in soil from the Py-R groups. We observed that the class Proteobacteria had the highest abundance (24.43%) in the rhizosphere soils of the three plants and soil samples from grassland plots, followed by Acidobacteria (21.09%), Verrucomicrobia (19.48%), Planctomycetes (9.20%), and Anctinobacteria (5.53%) ( Figure 5).
A Venn diagram ( Figure 6) illustrates that the bacterial communities of E. rufipilus (Er) shared 4095 OTUs between bulk (Er-B) and rhizosphere soils (Er-R) (Figure 6a); the bacterial communities of C. nepalensis (Cn) shared 3824 OTUs between bulk (Cn-B) and rhizosphere soils (Cn-R) (Figure 6b); the bacterial communities of P. yunnanensis (Py) shared 2817 OTUs between bulk (Py-B) and rhizosphere soils (Py-R) (Figure 6c); and the bacterial communities among the four groups (Er-R, Cn-R, Py-R, and CK) shared 2516 OTUs (Figure 6d). The OTUs derived from the soil associated with the three plant species showed that rhizosphere soils had more OTUs than bulk soils (Figure 6a-c).

Impact of Soil Properties on the Relative Abundances of Bacterial Taxa
We found a significant association among bacterial community structure and physicochemical properties in the soils (Figure 7). Bacterial communities from Py-R, Py-B, and Cn-R tend to group together, whereas Er-R, Er-B, Cn-B, and CK are scattered. Soil water content (SWC), soil organic matter (SOM), and total nitrogen (TN) were the main limiting factors, and the first two axes in RDA analysis explained 47.24% of variance for the relationship between soil bacterial community composition and physicochemical factors. Each point represents a sample, the closer the distance between two points, the higher the bacterial communities' similarity. The arrow shows the influence factor, the longer the arrow, the greater the effect. The angle between the factors for acute angle shows that the two factors are positively corrected; an angle that is an obtuse angle is a negative correlation. Notes: Er-R, rhizosphere soil of Erianthus rufipilus; Er-B, bulk soil of Erianthus rufipilus; Cn-R, Figure 7. Redundancy analysis (RDA) of bacterial communities as affected by the environmental and chemical properties factors, based on the relative abundance and similarity of bacterial phyla. RDA1 explained 30.78% and RDA2 16.46% variability. Each point represents a sample, the closer the distance between two points, the higher the bacterial communities' similarity. The arrow shows the influence factor, the longer the arrow, the greater the effect. The angle between the factors for acute angle shows that the two factors are positively corrected; an angle that is an obtuse angle is a negative correlation. Notes: Er-R, rhizosphere soil of Erianthus rufipilus; Er-B, bulk soil of Erianthus rufipilus; Cn-R, rhizosphere soil of Coriaria nepalensis; Cn-B, bulk soil of Coriaria nepalensis; Py-R, rhizosphere soil of Pinus yunnanensis; Py-B, bulk soil of Pinus yunnanensis; CK, soil from grassland plots; Tem, temperature; SWC, soil water content; SOM, soil organic matter; TN, total nitrogen; AN, alkali solution nitrogen; TP, total phosphorus; AP, available phosphorus; N/P, total nitrogen/total phosphorus.

Discussion
Results of this study demonstrate the relationships among soil physicochemical properties and the bacterial communities associated with three selected plants (E. rufipilus, C. nepalensis, and P. yunnanensis), as well as grassland plots in phosphorus-rich mountains in southwestern China. The associations among soil nutrients and bacterial community structure of rhizosphere and bulk soils has been demonstrated in several previous studies [48][49][50][51][52]. Plant species play an important role in shaping rhizosphere microbial communities through their root exudates [24,[53][54][55]. In this study, we found that the content of SOM, TN, and AN was present in higher concentrations in rhizosphere soils as compared with bulk soils among three plants in PES. The contents of SOM, TN, AN, TP, and AP in rhizosphere soils were elevated 55.31%, 34.11%, 39.02%, −2.11%, and 3.50% relative to bulk soil, consistent with soils from other studies and regions [56,57]. The rhizosphere soils of P. yunnanensis have higher TP content (elevated 121.44% and 189.06% relative to E. rufipilus and C. nepalensis, respectively), AP content (elevated 115.70% and 294.24% relative to E. rufipilus and C. nepalensis, respectively), and lower nitrogen/phosphorus (N/P) (declined 50.91% and 52.63% relative to E. rufipilus and C. nepalensis, respectively) than the other two dominant plants. Differences in soil phosphorus content among land-use types in PES could possibly be due to geological reasons; for example, this study site is home to geological phosphorus-rich mountains, and phosphate deposited in Yunnan Province accounts for 80% of high-grade phosphate ore in China [17]. According to previous studies on PES, phosphorus is found in high concentrations in soils with low soil nitrogen content [21,58]. Therefore, the results of this study are congruent with previous studies.
The 16S rRNA gene sequencing survey provides exhaustive information about the relative abundance, diversity, and composition of bacterial communities. This study was the first implementation of Illumina MiSeq high-throughput sequencing technology on soil microbial communities in phosphorus-rich mountains. The relative abundance of bacteria among the soil samples showed that P. yunnanensis had the lowest bacterial OTUs and CK had the highest. In contrast, soil nutrient properties revealed that P. yunnanensis had the highest concentration of SOM, AN, TP, and AP, whereas CK had the lowest SOM. This suggests that plant species shapes soil nutrient properties and bacterial community characteristics in rhizospheres and bulk soils.
The most abundant phylum of bacterial communities in the rhizosphere soil among plants was Proteobacteria. This result corroborates the results of several previous bacterial community studies done on the rhizosphere soils of crops [15,59], mining soils [60,61], and PES bacterial community research [28]. Furthermore, the class Alphaproteobacteria occupied the majority of the phylum Proteobacteria. Rhizobiales are the most abundant order in Proteobacteria across all soil samples. The families Bradyrhizobiaceae and Rhizobiaceae in the order Rhizobiales from C. nepalensis have higher abundance than E. rufipilus, P. yunnanensis, and grassland plots; this is because C. nepalensis forms symbiotic relationships with nitrogen-fixing bacteria [62,63].
The phylum Acidobacteria is the second most abundant community among the rhizosphere of the plants and grassland plots in this study. Acidobacteria is the predominant phylum of bacteria in semi-arid and other mature soil environments [64]. Acidobacteria exists most abundantly in the rhizosphere soil of E. rufipilus; this is likely because E. rufipilus is highly drought resistant and suitable to grow in arid regions [65,66]. As a dominant species, E. rufipilus is widely distributed in phosphorus-rich regions and has a high phosphorus tolerance [63]. The next phyla in highest abundance were Verrucomicrobia and Planctomycetes. These two phyla were abundant in the rhizosphere soil of P. yunnanensis and grassland plots. Planctomycetes represented a small minority of aquatic bacteria found in seawater and freshwater environments [67,68]. However, the dry season occurs in central Yunnan, China from November to April [69], and P. yunnanensis exhibits drought tolerance [70], which may explain why Planctomycetes was present in soils linked to P. yunnanensis in study areas. According to the rhizosphere bacterial communities of Y. yunnanensis reported in other regions of Yunnan Province where phosphorus is deficient, the dominant bacterial taxa are Acitinobacteria, Alphaproteobacteria, and Acidobacteira, whereas Verrucomicrobia and Planctomycetes were found in low abundance [71]. In addition to the plants reported in this study, phosphorus also alters bacterial community composition in soybean rhizosphere soils, and phosphorus fertilization decreased the relative abundance of Bacillales and Pseudomonadales [72].
Plant root exudates also shape rhizosphere bacterial structures [73][74][75]. Previous studies about phosphorus-enriched soils have concluded that phytoremediation is an effective measure for nutrient (nitrogen and phosphorus) removal [17,76]. Fewer studies, however, has been focused on the biodiversity of soil microbial communities based on vegetation regeneration. This study showed that bacterial communities in PES are strongly influenced by the rhizospheres of different plant species and physicochemical soil parameters. The distribution of soil bacteria was affected by multiple environmental factors. Soil macro-environmental conditions, such as soil water content, has the most notable impacts on the distribution of bacteria, followed by soil organic matter and total nitrogen. Although phosphate content was high in PES, total phosphorus was found to have the lowest impact as compared with other nutrient factors on bacterial communities. The findings of this study verified the results of a previous study on low nitrogen levels in the region [19].
On a larger scale, in the phosphorus-rich and nitrogen-deficient environments of PES, most vegetation is secondary, and habitats are relatively fragile. E. rufipilus, C. nepalensis, and P. yunnanensis are widely grown in phosphorus-enriched soils of China because of their high tolerance to drought and arid conditions. The P. yunnanensis plant communities have higher plant diversity in PES, whereas bacterial communities tend to feature lower diversity under the forest soil of P. yunnanensis. Higher plant diversity and lower litter decomposition efficiency [77], as well as litter inhibitory effect on bacterial communities such as production of tannin [78] occurring under P. yunnanensis as compared with E. rufipilus and C. nepalensis likely explains this lower level of diversity. In addition, P. yunnanensis is a common native secondary coniferous species growing in low temperatures. Because of its rapid growth and higher plant community diversity beneath trees, it is a preferred species in Yunnan [37]. On the basis of the above findings, follow-up research in plant rhizosphere microbial compositions, and diversity could be fruitful. Nutrient properties and microbial communities can provide guidance for understanding biochemical processes and ecosystem functionalities in the soils of phosphorus-rich mountains.

Conclusions
The richness and diversity of the bacterial communities reported in our study reveal prospective uncultured bacteria derived from rhizosphere soils of dominant plants in PES for the first time. Differences in the physicochemistry and structures of soil bacterial communities among the dominant plants form useful guidelines for examining how tree species and soil physicochemical parameters influence local soil environments in phosphorus-rich regions. The results from our study show that both bacterial communities and soil nutrients show Rhizospheric effects; rhizosphere bacteria are mainly composed of Proteobacteria and Acidobacteria, followed by Verrucomicrobia and Planctomycetes; bacteria abundance and diversity in PES are mainly influenced by soil water content, soil organic matter, and total nitrogen. Our results indicate that rhizosphere microbial communities could serve as an important index for phytoremediation. This is of great significance for the restoration of disturbed ecosystems in PES.