Comparisons of Soil Properties , Enzyme Activities and Microbial Communities in Heavy Metal Contaminated Bulk and Rhizosphere Soils of Robinia pseudoacacia L . in the Northern Foot of Qinling Mountain

The toxic effects of heavy metal (HM) contamination on plant metabolism and soil microorganisms have been emphasized recently; however, little is known about the differences in soil physical, chemical, and biological properties between bulk and rhizosphere soils contaminated with HMs in forest ecosystem. The present study was conducted to evaluate the rhizosphere effect on soil properties, enzyme activities and bacterial communities associated with Robinia pseudoacacia L. along a HM contamination gradient. Soil organic matter (SOM), available nitrogen (AN) and phosphorus (AP) contents were significantly higher in rhizosphere soil than those in bulk soil at HM contaminated sites (p < 0.05). Compared to bulk soil, activities of four soil enzymes indicative of C cycle (β-glucosidase), N cycle (protease, urease) and P cycle (alkaline phosphatase) in rhizosphere soil across all study sites increased by 47.5%, 64.1%, 52.9% and 103.8%, respectively. Quantitative PCR (qPCR) and restriction fragment length polymorphism (RFLP) were used to determine the relative abundance, composition and diversity of bacteria in both bulk and rhizosphere soils, respectively. The copy number of bacterial 16S rRNA gene in bulk soil was significantly lower than that in rhizosphere soil (p < 0.05), and it had significantly negative correlations with total/DTPA-extractable Pb concentrations (p < 0.01). Alphaproteobacteria, Gammaproteobacteria and Firmicutes were the most dominant groups of bacteria at different study sites. The bacterial diversity index of Species richness (S) and Margalef (dMa) were significantly higher in rhizosphere soil compared with those in bulk soil, although no difference could be found in Simpson index (D) between bulk and rhizosphere soils (p > 0.05). Redundancy analysis (RDA) results showed that soil pH, EC, SOM and total/DTPA-extractable Pb concentrations were the most important variables affecting relative abundance, composition and diversity of bacteria (p < 0.05). Our study highlights the importance of rhizosphere effect on soil nutrient content, enzyme activity, bacterial abundance and community in HM contaminated forest soils. Further study is still required to understand the specific processes in the rhizosphere to achieve a suitable rhizosphere biotechnology for restoration of degraded forest ecosystem.


Introduction
Soil heavy metal (HM) contamination has become one of the most severe environmental issues in many parts of the world because of their toxic effects on the environment and human health [1][2][3].HMs in soil are derived from both the natural (parent material) and various anthropogenic sources [4].The major anthropogenic sources of this pollution are mining, smelting, industrial activities and the application of fertilizers or pesticides [5].In recent years, with the rapid development of society and economy, soil polluted with HMs has become a serious and widespread environmental issue in China [6,7].Unlike organic pollutants, HMs are more toxic, non-biodegradable in the environment and can easily accumulated and magnified in organisms [8].This kind of pollutant not only degrades ecosystem services but also threatens human health through the food chain [9].
Soil enzymes are mainly of microbial origin, being derived from intracellular, cell-associated or free enzymes [10].They are important in catalyzing several vital reactions necessary for the life processes of microorganisms in soils and the stabilization of soil structure, the decomposition of organic wastes, organic matter formation, and nutrient cycling, hence playing an important role in maintaining soil ecology, physical and chemical properties, fertility, and soil health [11].Soil enzyme activities, reacting faster than physical variables and/or after any chemical change in the soil, have been generally considered as one of a good bioindicators of soil fertility quality and can quantify changes as a result of the natural and anthropogenic disturbance in the soil ecosystem [12].Therefore, a better understanding of the effects of HM on soil enzyme activities will potentially provide an opportunity for an integrated assessment of soil biology.Vig et al. [13] indicated that HMs could inhibit soil enzyme activities through interacting with the enzyme active sites and substrate complexes, and denaturing the enzyme protein.
The results of Brookes [14] reported that dehydrogenase activity could be considered as an indicator of soil contaminated with HM.Liao and Huang [15] have succeeded in using urease, protease and acid phosphatase as indices of soil contamination by HMs.Although there are increasing studies on soil enzymes, little information has been reported on the link between soil enzyme activities and bacterial community structure along the gradient of HM concentration in black-locust (Robinia pseudoacacia L.) forest on the Loess Plateau.
However, the influence of HM on soil enzyme activity depends on pH, nutrient form and amount, HM concentration and availability, enzyme type, etc. [16].Acosta-Martinez and Tabatabai [17] reported that β-glucosidase was sensitive to pH changes, which can be used as a good biochemical indicator for evaluating ecological changes resulting from soil acidification.Apart from being good indicators of soil fertility, phosphatase secretion from both plant roots and soil microorganisms can increase to improve the solubilization and remobilization of phosphate and affect the ability of the plant to cope with P-stressed conditions [18].Studies have shown that urease was very sensitive to toxic concentrations of HMs, and its activity generally increased with increasing temperature [19].
On the other hand, the rhizosphere is an interface or zone of soil surrounding to plant roots and has important effects on plant health and soil quality.The living plant roots are able to secrete large quantity of organic compounds that promote the growth of bacteria, resulting in a higher bacterial community, diversity and enzyme activities in rhizosphere compared with bulk soil [20].In turn, these beneficial microorganisms are capable of improving plant growth through enhancing nutrient uptake, producing auxins-like substances and protecting them from abiotic stresses [21].The net carbon and nitrogen mineralization and nitrification rates, activities of extracellular enzymes for nutrient cycling in rhizosphere are generally higher than these in bulk soil [22,23].Therefore, to better understand the mechanism of phytoremediation and highlight the importance of rhizosphere effects on ecological restoration, it is essential to compare the differences in microbial composition, soil properties and extracellular enzyme activities between rhizosphere and bulk soils in HM contaminated soils.
In our previous investigation, R. pseudoacacia had been identified as lead/zinc tolerant dominant tree species widely distributed in HM contaminated soils in Feng County [24].This tree was considered as a suitable candidate for phytoremediation owing to its developed root system, high biomass productivity, environmental stress tolerance and ability of atmospheric nitrogen fixation [24][25][26].
However, the information about the rhizosphere effects of R. pseudoacacia on bacterial community, soil properties and activities extracellular enzymes is still unknown, which limits our understanding of phytoremediation processes in HM polluted soils.Therefore, the objectives of this study were to (1) determine the differences of soil properties, enzyme activities, bacterial abundance and diversity between bulk and rhizosphere soils; (2) evaluate the impacts of rhizosphere and HM pollution on soil properties, enzyme activities, bacterial abundance and diversity; and (3) reveal the influences of soil properties on bacterial composition in both HM contaminated rhizosphere and bulk soils in Feng County, Northwest China.The findings of this study attempted to provide a better understanding of and new insights into the influences of rhizosphere effects on soil biology in Pb-Zn mine area.

Study Area and Sampling
This study was conducted in Qiandongshan Pb-Zn mine region, located in the Northern foot of Qinling Mountain, Shaanxi Province, Northwest part of China (106 • 38 E, 33 • 49 N).This region is in the warm temperate semiarid zone with an annual average temperature of 11.4 • C. The annual average rainfall is 613.2 mm with a frost-free period of 188 days.The soil texture was classified as cinnamon and brunisolic soil, according to the traditional soil genesis classification standard in China [27].Qiandongshan Pb-Zn region is the largest and the most typical five Pb and Zn producers in China, and its total reserve is approximately 1.12 × 10 6 tons [28].
The HM polluted area had the same soil type and was divided into polluted levels of light (S1, 33  S1), according to our previous investigation and the Chinese environmental quality standard (Grade II) for soils (GB 15618-1995) [27,[29][30][31].In August 2012, four study plots (15 × 15 m) dominated by target tree species (R. pseudoacacia) were identified at each site for soil sampling.The pits at three subplots within each plot were excavated to collect sufficient amount of rhizosphere soil and bulk soil.Samples collected from subplots were finally combined into one composite sample, and the parameters were measured in quadruplicate for each study site.Individual soil samples were kept in separate plastic bags and placed in a cooler with ice until they arrived at the laboratory for processing.In the laboratory, the volume of soil adhering to closely to plant fine roots defined as rhizosphere soil, while the remaining soil classified as bulk soil.Portions of rhizosphere and bulk soil samples were used for DNA extraction and analysis of enzyme activity, and the remaining soil samples were used for soil properties determination.

Soil Properties
After air-drying, the soil samples were ground to a fine powder and sieved through an 80-mesh sieve using either a mortar and pestle or freezer mill.Soil pH is determined on a 1:2.5 soil:deionized water suspension using a digital pH meter (Leici PHS-3C, Shanghai, China), and soil electronic conductivity (EC) was determined by a conductivity meter (DDSJ-308A, Hangzhou, Zhejiang, China).Soil organic matter (SOM) was measured colorimetrically after wet digestion by K 2 Cr 2 O 7 [32].The determination of soil organic carbon (SOC) was carried out using total organic carbon analyzer (TOC-VCPH, Shimadzu, Japan) according to dry combustion method [33].Total soil phosphorus (TP) was digested by HF-HClO 4 and then measured by the flame photometric method [34].Available phosphorus (AP) and available nitrogen (AN) were determined according to the method of Page [35].Total concentration of HMs was determined by aqua regia (HNO 3 /HCl = 1:3) and HClO 4 extraction method and then determined by flame atomic absorption spectrometry (FAAS, Hitachi Z-2000, Tokyo, Japan).The DTPA-extractable Pb and Zn concentrations were extracted from soil with DTPA solution (0.005 mol L −1 diethylene triamine penta-acetic acid (DTPA), 0.01 mol L −1 CaCl 2 , 0.1 mol L −1 triethanolamine) with pH adjusted to 7.3, and the extracts were analyzed for Pb and Zn using FAAS method [36].The quality assurance was achieved using blank reagent and standard reference soil, and the results were determined in triplicate for each soil sample.

Soil Enzyme Activities
Activities of four soil enzymes indicative of C cycle (β-glucosidase), N cycle (protease, urease) and P cycle (alkaline phosphatase) were determined.The activity of β-glucosidase (µg PNP g −1 h −1 ) was measured by release and detection of p-nitrophenol (PNP) using p-nitrophenyl-β-glucopyranoside as substrate at 400 nm in an UV spectrophotometer [37].Protease activity (µg TYR g −1 h −1 ) was determined as describe by Ladd and Butler [38].After the addition of casein (2% w/v) to 1 g of substrate, tyrosine (TYR) formed and reacted with Folin-phenol reagent to generate a blue complex that was determined spectrophotometrically at 700 nm.The urease activity (µg NH 4 was assayed according to the method of Kandeler and Gerber [39].After the addition of aqueous (controls) or a buffered urea solution (samples) to 5 g of soil samples, incubation took place for 2 h at 37 • C. The measurement was based on the reaction of sodium salicylate with NH 3 in the presence of sodium dichloroisocyanurate, which forms a green-colored complex under alkaline pH conditions.The concentration of NH 4 + -N extracted with 2 M KCl solution was determined spectrophotometrically at 690 nm.The activity of alkaline phosphatase (µg PNP g −1 h −1 ) was measured according to Tabatabai and Bremner [40] after 1 h soil incubation at 37 • C with 0.025 M p-nitrophenol phosphate disodium in borate buffer at pH 10.The p-nitrophenol (PNP) in the filtrate was determined colorimetrically at 410 nm using a spectrophotometer.

DNA Extraction
Rhizosphere and bulk soil samples were milled to powder using a ceramic pestle and mortar in liquid nitrogen, respectively.DNA was extracted from approximately 0.5 g of soil using E.Z.N.A. soil DNA kit (Omega Bio-tek Inc., Norcross, GA, USA) according to manufacturer s instructions.The DNA extract quantity and quality were determined by Smartspec Plus Spectrophotometer (Bio-Rad, Hercules, CA, USA) and agarose gel electrophoresis (Invitrogen, Carlsbad, CA, USA).The DNA was then diluted 10-fold with distilled ddH 2 O and used for PCR template.

Quantitative PCR
Bacterial small subunit rRNA gene amplified from total soil DNA extracts with the forward primer Eub338 and reverse primer Eub518 using the method described by Fierer et al. [41].To estimate bacterial small-subunit rRNA gene abundances, standard curve was generated using a 10-fold serial dilution of a plasmid containing a full-length copy of Escherichia coli 16S rRNA gene.The 25 µL qPCR reactions contained 12.5 µL 2 × Taq SYBR Premix Ex Taq™ II (Takara Biotechnology Co., Ltd., Dalian, China), 0.5 µL of each 10 µM forward and reverse primers (Invitrogen), and 9.5 µL sterile RNase-free water.Standard and environmental DNA samples were added at 2.0 µL per reaction.The reaction was carried out on a CFX Connect Real-Time PCR Detection System (Bio-Rad, USA) using a program of 95 • C for 10 min, followed by 40 cycles of 95 • C for 30 s, 55 • C for 30 s and 72 • C for 30 s according to our previous study [42].Specific amplification of the target gene was confirmed on a 1% (w/v) agarose gel showing the band of the expected size and melting curve analysis resulting in a single peak.Bacterial gene copy numbers were generated using a regression equation for each assay relating the cycle threshold (C t ) value to the known number of copies in the standards.The final copy numbers of bacterial 16S rRNA gene were obtained by calibrating against the total DNA concentration extracted and the soil water content (dry weight soil, DWS).All qPCR reactions were run in quadruplicate with the DNA extracted from each soil sample.

RFLP Analysis
Partial bacterial 16S rRNA genes (about 1500 bp) were amplified from total soil DNA extracts with the universal primer 1492R and primer 27F [43].The amplifications were performed in a reaction volume of 25 µL including 12.5 µL 2 × Taq PCR Master Mix (CoWin Biotech Co., Ltd., Beijing, China), 11 µL sterile distilled, 0.5 µL each primer and 0.5 µL DNA template.The thermal cycling conditions consisted of an initial denaturation at 94 • C for 4 min, and 30 cycles of 94 • C for 60 s, 55 • C for 60 s, 72 • C for 90 s with a final extension at 72 • C for 10 min and was run in BioRad PCR cycler (USA).The resulting PCR products were analyzed by 1.5% agarose gel electrophoresis and purified using a universal DNA purification kit (Tiangen Biotech Co., Ltd., Beijing, China).The purified PCR products were ligated into pGEM-T cloning vector according to manufacturer's instructions (Tiangen Biotech Co., Ltd., Beijing, China).A total of 24 clone libraries (12 for rhizosphere soil and 12 for bulk soil) were constructed in this study.For each library, approximately 50 white positive clones were picked up randomly and then screened for the presence of a complete insert (1500 bp) by reamplification and gel electrophoresis.Only clones with complete inserts were digested (3 h, 37 • C) with 5 U of the restriction enzymes MspI and HaeIII and visualized on 4% agarose gels.Respective fragments were visualized under UV light and classified into two categories depending on the absence and present of bands for statistical analyses (0, no band; 1, has band).Representative clones of each RFLP patterns were selected and sequenced unidirectionally by GenScript Corporation (Nanjing, China).

Data Analysis
The restriction fragments with different lengths were compared among samples.Although the use of multiple enzyme digests increases the specificity of data generated by RFLP analysis, different organisms can still produce fragments of similar lengths.Therefore, each fragment is assumed to represent different unique operational taxonomic units (OTU), but not necessarily a distinct species.Bacterial 16S rRNA sequence with > 97% match were classified as the species, the sequences with 95-97% similarity were identified to genus, and the sequences with < 95% match were considered as unknown bacteria based on sequences matched in the NCBI GenBank database.Phylogenetic tree was generated with MEGA version 5.0, using the neighbor-joining (NJ) and bootstrap values from 1000 replicates.The 31 gene sequences obtained in the current study had been deposited into the GenBank database with the following accession numbers: KP162186-KP162216.
OTU richness was determined for each sample as the total number of fragments identified in that sample.OTU diversity and evenness for each sample were calculated with the Margalef Index (d Ma ) and Simpson Index (D), which were calculated using the following Equations: where S is the total number of phylotypes found at each site, N is the total number of clones selected from each clone library, n i is the number of the ith representative sequence in each clone library.A t-test was performed to compare the average values of soil properties, enzyme activities, bacterial abundances between bulk and rhizosphere soils, while pearson correlation coefficients were used to determine their correlations, which were performed by SPSS 16.0 in Windows 7 (SPSS Inc., Chicago, IL, USA).Data expressed in the figures and tables are presented as the means with a standard deviation (SD).A multivariate redundancy analysis (RDA) was performed by Canoco software (version 4.5, Centre for Biometry, Wageningen, The Netherlands) as a direct gradient analysis to explore the relationships among bacterial community compositions and environmental factors.The effect and significance of each environmental factor on the bacterial abundance and diversity were calculated through forward selection and the Monte Carlo permutation test.

Soil Properties
The main chemical properties of soil samples collected from different study sites are shown in Table 1.The soils were all slightly alkaline and soil pH in rhizosphere soil was significantly lower than bulk soil, ranging from 7.99 (S4) to 8.62 (S2) and from 7.33 (S4) to 7.58 (S2) in bulk and rhizosphere soils, respectively.The EC values at S3 (0.80) and S4 (0.82) in rhizosphere soil were significant higher compared with these in bulk soil.The SOM, SOC, TN, AN and AP contents in both bulk and rhizosphere soils at S5 were significant lower than these at S1, and the contents of AN and AP at all study sites tended to be higher in rhizosphere than in bulk soil except for AN at S1 (14.9 mg/kg for bulk soil and 16.9 mg/kg for rhizosphere soil).However, no significant differences could be detected in EC, SOM, SOC, TN, AN and TP between rhizosphere and bulk soils at S1 (p > 0.05).The soil HM concentrations had a large variability among different HM pollution levels (Table 1).S4 had the highest HM contamination level, and the total concentrations of Pb and Zn were 9.5 and 1.4 times Grade II of the national environmental quality standard for soils (GB 15618-1995), while S1 was only slightly polluted by Pb (223 mg/kg) and Zn (88.4 mg/kg).The DTPA-extractable Pb and Zn showed the same pattern as total Pb and Zn.There was an increasing trend in total and DTPA-extractable Pb followed the order S1 < S2 < S3 < S4, while the total and DTPA-extractable Zn decreased in the different order of S3 > S4 > S2 > S1 in both rhizosphere and bulk soil.The concentrations of DTPA-extractable in relation to the total concentrations were 5.53% for Pb and 14.3% for Zn when all soil samples were taken into account.
The significantly negative correlations among SOM, SOC, TN, AN, AP and total/ DTPA-extractable Pb could be found in both bulk and rhizosphere soils (p < 0.05) except for TN and total Pb in bulk soil (p > 0.05), however, no correlation could be detected among SOM, AN and total/DTPA-extractable Zn in rhizosphere soil (p > 0.05).Soil EC was positively correlated with total/DTPA-extractable Pb and Zn in both bulk and rhizosphere soils, while pH only exhibited significantly negative correlation with total and DTPA-extractable Pb in bulk soil (p < 0.05) (Table S1).

Soil Enzyme Activities
Variations of the enzyme activities in bulk and rhizosphere soils samples at different study sites are shown in Figure 1.For bulk soil, the activities of β-glucosidase, protease and urease decreased by 60.8%, 48.3% and 59.5%, respectively, from light pollution level to high pollution level.On the contrary, the activities of β-glucosidase, protease and urease reduced by 34.0%, 36.9% and 37.9%, respectively.However, the activities of alkaline phosphatase at bulk and rhizosphere soils were only decreased by 24.0% and 17.1%, respectively.Both urease and alkaline phosphatase activities in rhizosphere soil were significantly larger than those in bulk soil at different study sites except for activity of urease at medium pollution level.However, no difference could be found in β-glucosidase activity between light and low pollution levels.
Forests 2017, 8, 430 7 of 18 and total Pb in bulk soil (p > 0.05), however, no correlation could be detected among SOM, AN and total/DTPA-extractable Zn in rhizosphere soil (p > 0.05).Soil EC was positively correlated with total/DTPA-extractable Pb and Zn in both bulk and rhizosphere soils, while pH only exhibited significantly negative correlation with total and DTPA-extractable Pb in bulk soil (p < 0.05) (Table S1).

Soil Enzyme Activities
Variations of the enzyme activities in bulk and rhizosphere soils samples at different study sites are shown in Figure 1.For bulk soil, the activities of β-glucosidase, protease and urease decreased by 60.8%, 48.3% and 59.5%, respectively, from light pollution level to high pollution level.On the contrary, the activities of β-glucosidase, protease and urease reduced by 34.0%, 36.9% and 37.9%, respectively.However, the activities of alkaline phosphatase at bulk and rhizosphere soils were only decreased by 24.0% and 17.1%, respectively.Both urease and alkaline phosphatase activities in rhizosphere soil were significantly larger than those in bulk soil at different study sites except for activity of urease at medium pollution level.However, no difference could be found in β-glucosidase activity between light and low pollution levels.The enzyme activities of β-glucosidase, protease, and urease were positively correlated to pH, SOM, SOC, AN and AP (p < 0.05), but had a negatively correlation with EC (p < 0.01) in bulk soil.The significantly positive correlations among the enzyme activities of β-glucosidase, protease, urease and AN, TP and AP in rhizosphere soil (p < 0.05) were identified, but no relationship could be detected among them and pH (p > 0.05).Total/DTPA-extractable Pb and Zn greatly inhibited the enzyme activities of β-glucosidase, protease and urease in both bulk and rhizosphere soils (p < 0.05), The enzyme activities of β-glucosidase, protease, and urease were positively correlated to pH, SOM, SOC, AN and AP (p < 0.05), but had a negatively correlation with EC (p < 0.01) in bulk soil.The significantly positive correlations among the enzyme activities of β-glucosidase, protease, urease and AN, TP and AP in rhizosphere soil (p < 0.05) were identified, but no relationship could be detected among them and pH (p > 0.05).Total/DTPA-extractable Pb and Zn greatly inhibited the enzyme activities of β-glucosidase, protease and urease in both bulk and rhizosphere soils (p < 0.05), however, no effects could be found on alkaline phosphatase activity (p > 0.05) (Table S2).

Soil Bacterial Abundance
The abundance of the soil bacterial community and changes at different HM pollution levels were evaluated using qPCR (Figure 2).The soil bacterial copy number decreased by 73.8% at bulk soil and 41.6% at rhizosphere soil, respectively, from light to high pollution levels.The abundance of the soil bacterial community in bulk soil was significantly lower compared with that in rhizosphere soil at low, medium and high pollution levels (p < 0.01).However, no difference could be found in bacterial copy number between bulk and rhizosphere soils at light pollution levels (p > 0.05).

Soil Bacterial Abundance
The abundance of the soil bacterial community and changes at different HM pollution levels were evaluated using qPCR (Figure 2).The soil bacterial copy number decreased by 73.8% at bulk soil and 41.6% at rhizosphere soil, respectively, from light to high pollution levels.The abundance of the soil bacterial community in bulk soil was significantly lower compared with that in rhizosphere soil at low, medium and high pollution levels (p < 0.01).However, no difference could be found in bacterial copy number between bulk and rhizosphere soils at light pollution levels (p > 0.05).Soil pH, SOM, SOC were significantly positively correlated with bacterial abundance in bulk soil (p < 0.05); however, no correlation could be detected among them and bacterial abundance in rhizosphere soil (p > 0.05).The copy numbers of bacterial 16S rRNA gene in both bulk and rhizosphere soils were significantly inhibited by total/DTPA-extractable Pb concentration (p < 0.05), but there was no correlation between bacterial abundance and total/DTPA-extractable Zn in rhizosphere soil (p > 0.05) (Table S3).

Phylogenetic Analysis
The bacterial occurrence was determined in 12 bulk and 12 rhizosphere soil samples.Template DNA was successfully amplified with the universal PCR primers 27F and 1492R, yielding PCR products of expected length.In total, 1200 clones were screened and 1171 clones containing insert of correct size of 16S rRNA gene were digested by MspI and HaeIII.One to five RFLP representatives were sequenced, yielding 31 bacterial sequences (Table S4).All representatives were deposited to the NCBI database and used for phylogenetic analysis (Figure 3).
The phylogenetic tree based on 31 representatives of the predominant RFLP types indicated that the bacterial sequences were closely related to 11 known classes: Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Deltaproteobacteria, Bacteroidetes, Bacillacea, Firmicutes, Actinobacteria, Acidobacteria, Chloroflexi and Planctomycetes.Three major classes, Soil pH, SOM, SOC were significantly positively correlated with bacterial abundance in bulk soil (p < 0.05); however, no correlation could be detected among them and bacterial abundance in rhizosphere soil (p > 0.05).The copy numbers of bacterial 16S rRNA gene in both bulk and rhizosphere soils were significantly inhibited by total/DTPA-extractable Pb concentration (p < 0.05), but there was no correlation between bacterial abundance and total/DTPA-extractable Zn in rhizosphere soil (p > 0.05) (Table S3).

Phylogenetic Analysis
The bacterial occurrence was determined in 12 bulk and 12 rhizosphere soil samples.Template DNA was successfully amplified with the universal PCR primers 27F and 1492R, yielding PCR products of expected length.In total, 1200 clones were screened and 1171 clones containing insert of correct size of 16S rRNA gene were digested by MspI and HaeIII.One to five RFLP representatives were sequenced, yielding 31 bacterial sequences (Table S4).All representatives were deposited to the NCBI database and used for phylogenetic analysis (Figure 3).represented, accounting for 54.0-84.0% of the sequenced flora.Whatever the pollution level, these two classes were found in both bulk and rhizosphere soils samples with high percentage at four study sites (Figure 3).Another class, namely Bacillacea was found in rhizosphere soil samples, accounting for 28.7-37.3% of the sequenced flora.The phylogenetic tree based on 31 representatives of the predominant RFLP types indicated that the bacterial sequences were closely related to 11 known classes: Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Deltaproteobacteria, Bacteroidetes, Bacillacea, Firmicutes, Actinobacteria, Acidobacteria, Chloroflexi and Planctomycetes.Three major classes, namely Alphaproteobacteria, Gammaproteobacteria and Firmicutes, were particularly well represented, accounting for 54.0-84.0% of the sequenced flora.Whatever the pollution level, these two classes were found in both bulk and rhizosphere soils samples with high percentage at four study sites (Figure 3).Another class, namely Bacillacea was found in rhizosphere soil samples, accounting for 28.7-37.3% of the sequenced flora.

Bacterial Communities in Rhizosphere Soil
The bacterial composition was dominated by Alphaproteobacteria, Gammaproteobacteria, and Firmicutes across all 12 soil samples analyzed, while other taxa presented at relatively lower abundance (Figure 4).Of the bacterial classes detected by the RFLP analysis, eight were found in rhizosphere soil associated with R. pseudoacacia grown at light pollution level, with an average of 2.7 classes per sample, while only five classes were detected at high pollution level, with an average 1.7 classes per sample.The class of Alphaproteobacteria and Bacillacea followed a clear decreasing gradient with increasing HM pollution levels, ranging from 7.3-26.0%and 14.7-19.3% of all clones, respectively.However, the percentage of Gammaproteobacteria and Firmicutes increased with increasing of soil HM concentrations, ranging from 18.0-34.7%and 10.7-22.0% of all clones, respectively.Actinobacteria and Chloroflexi could not be found in the rhizosphere soil at different study sites.It is interesting to notice that Betaproteobacteria and Planctomycetes were only found at light and low pollution levels, however, with increasing of HM concentrations, they disappeared at medium and high pollution levels (Figure 4).pollution levels.The access number for each division is shown on the tree.Bootstrap values (1000 replicates) higher than 50% are shown.Scale bar represents the 5% substitution percentage.

Bacterial Communities in Rhizosphere Soil
The bacterial composition was dominated by Alphaproteobacteria, Gammaproteobacteria, and Firmicutes across all 12 soil samples analyzed, while other taxa presented at relatively lower abundance (Figure 4).Of the bacterial classes detected by the RFLP analysis, eight were found in rhizosphere soil associated with R. pseudoacacia grown at light pollution level, with an average of 2.7 classes per sample, while only five classes were detected at high pollution level, with an average 1.7 classes per sample.The class of Alphaproteobacteria and Bacillacea followed a clear decreasing gradient with increasing HM pollution levels, ranging from 7.3-26.0%and 14.7-19.3% of all clones, respectively.However, the percentage of Gammaproteobacteria and Firmicutes increased with increasing of soil HM concentrations, ranging from 18.0-34.7%and 10.7-22.0% of all clones, respectively.Actinobacteria and Chloroflexi could not be found in the rhizosphere soil at different study sites.It is interesting to notice that Betaproteobacteria and Planctomycetes were only found at light and low pollution levels, however, with increasing of HM concentrations, they disappeared at medium and high pollution levels (Figure 4).

Bacterial Communities in Bulk Soil
Bulk soil presented a characteristic bacterial community that was different from that of the rhizosphere soil.Significant differences were found in the relative abundance of four bacterial classes.Class Gammaproteobacteria is of particular interest as it was found more frequently in bulk soils compared with that in rhizosphere soil samples, accounting for 46.7% of all clones at high pollution level.However, Bacillacea was seldom found in bulk soil samples, and only 2.0-10.0% of all clones had this class at different HM pollution levels (S1-S4).Two classes, namely Bacteroidetes and Planctomycetes could not be found in the bulk soil samples collected from different study sites.

Bacterial Communities in Bulk Soil
Bulk soil presented a characteristic bacterial community that was different from that of the rhizosphere soil.Significant differences were found in the relative abundance of four bacterial classes.Class Gammaproteobacteria is of particular interest as it was found more frequently in bulk soils compared with that in rhizosphere soil samples, accounting for 46.7% of all clones at high pollution level.However, Bacillacea was seldom found in bulk soil samples, and only 2.0-10.0% of all clones had this class at different HM pollution levels (S1-S4).Two classes, namely Bacteroidetes and Planctomycetes could not be found in the bulk soil samples collected from different study sites.The three most abundant bacterial classes in rhizosphere soil samples collected from four study sites were Alphaproteobacteria, Gammaproteobacteria and Firmicutes.The class of Alphaproteobacteria followed a clear decreasing gradient with increasing of HM pollution levels, ranging from 40.0-9.3% of all clones.However, the percentage of Gammaproteobacteria and Firmicutes increased with increasing of soil HM concentrations, ranging from 11.7-46.3%and 7.3-28.0% of all clones, respectively.Two classes, namely Deltaproteobacteria and Acidobacteria, could be found in the rhizosphere soil samples collected from light, low and medium pollution levels, however, with increasing HM concentrations, these two classes disappeared at high pollution level (Figure 4).

Bacterial Diversity Index
The bacterial diversity index of Species richness (S), Margalef (d Ma ) and Simpson index (D) varied in either bulk or rhizosphere soil associated with R. pseudoacacia grown at different HM pollution levels (Figure 5).The S and d Ma in bulk and rhizosphere soils decreased by 73.8% and 41.6%, respectively, from light and high pollution levels.There was no difference in S, d Ma and D between bulk and rhizosphere soils at light pollution level (p > 0.05), although the significant difference in D could be detected between bulk and rhizosphere soils at medium pollution level (p < 0.05).
Forests 2017, 8, 430 11 of 18 increasing of soil HM concentrations, ranging from 11.7-46.3%and 7.3-28.0% of all clones, respectively.Two classes, namely Deltaproteobacteria and Acidobacteria, could be found in the rhizosphere soil samples collected from light, low and medium pollution levels, however, with increasing HM concentrations, these two classes disappeared at high pollution level (Figure 4).

Bacterial Diversity Index
The bacterial diversity index of Species richness (S), Margalef (dMa) and Simpson index (D) varied in either bulk or rhizosphere soil associated with R. pseudoacacia grown at different HM pollution levels (Figure 5).The S and dMa in bulk and rhizosphere soils decreased by 73.8% and 41.6%, respectively, from light and high pollution levels.There was no difference in S, dMa and D between bulk and rhizosphere soils at light pollution level (p > 0.05), although the significant difference in D could be detected between bulk and rhizosphere soils at medium pollution level (p < 0.05).Potential differences in parameters between bulk and rhizosphere soils were analyzed using t-test (** p < 0.01; * p < 0.05; ns, not significant).

Discussion
The maximum total concentrations of Pb and Zn in both bulk and rhizosphere soils at S1 did not exceed the corresponding limits of the China Environmental Quality Standard for Soils (GB15618-1995, Grade II), suggesting that the HM concentrations at S1 were normal (Table 1).However, total Pb concentration at S2, S3 and S4 were 2.63, 6.87 and 9.70 times in bulk soil or 2.14, 4.51 and 7.68 times in rhizosphere soil above grade II of GB15618-1995 respectively, indicating that these sites were heavily polluted by HM according to either national grade II standard.Soil Zn concentration satisfied grade II quality at S1 and S2, but exceed the limits of GB15618-1995 at S3 and S4 in both bulk and rhizosphere soils.Our results reflected that significant accumulations of HM in in bulk and rhizosphere soils associated with R. pseudoacacia grown at different HM pollution levels.Potential differences in parameters between bulk and rhizosphere soils were analyzed using t-test (* p < 0.05; ns, not significant).

Discussion
The maximum total concentrations of Pb and Zn in both bulk and rhizosphere soils at S1 did not exceed the corresponding limits of the China Environmental Quality Standard for Soils (GB15618-1995, Grade II), suggesting that the HM concentrations at S1 were normal (Table 1).However, total Pb concentration at S2, S3 and S4 were 2.63, 6.87 and 9.70 times in bulk soil or 2.14, 4.51 and 7.68 times in rhizosphere soil above grade II of GB15618-1995 respectively, indicating that these sites were heavily polluted by HM according to either national grade II standard.Soil Zn concentration satisfied grade II quality at S1 and S2, but exceed the limits of GB15618-1995 at S3 and S4 in both bulk and rhizosphere soils.Our results reflected that significant accumulations of HM in soils could be attributed to smelting and mining activities, which probably pose potential risks to human health and the environment [6].
It has been widely accepted that availability rather than total HM concentration determines the toxic effects of a metal on biological systems [44].In the current study, HM availability was considered as the ratio of the DTPA-extractable concentration to the corresponding total concentration to characterize their toxicity.Pb availability in either bulk soil (12.7%) or rhizosphere soil (16.0%) was significantly lower compared with that of Zn (Figure S2), suggesting the higher solubility of Zn than that of Pb in soils, which was consistent with the finding of Fritsch et al. [45] in HM contaminated soils around mining area.However, it is interesting to note that the availability of Pb was significantly higher in bulk soil than that in rhizosphere soil at medium and high pollution levels, while a contrary tendency for Zn could be detected at medium pollution level (Figure S2).The results indicated that different types of HMs had various patterns of behavior and mobility between bulk and rhizosphere soils and along HM contamination gradient.Zn is an essential metal for plant growth and plays an important role in photosynthesis and enzyme composition for protein synthesis [46]; however, Pb is a non-essential element and has toxic effects on plant growth (membrane damage and oxidative stress) [47].Therefore, plants may develop their own strategy to survive in HM polluted soil: reducing availability of toxic HM (Pb) in rhizosphere and preventing the translocation of Pb from soil to the aerial parts of plants, which was consistent with our previous study [24,31].Various studies have indicated that plant root exudates could react with HM ions and affected metal solubility, mobility and phytoavailability [48].Luo et al. [49] reported that plant had ability to release low molecular weight organic acids such as oxalate into the soil, forming complexes with Pb and then immobilizing its availability in rhizosphere soil.
Soil enzymes are the main regulators of the biochemical processes within the soil environment and have been widely used as indicators for environmental quality, such as productivity, fertility and nutrient cycling potential [50].Our study showed that the activities of β-glucosidase, protease and urease in both bulk and rhizosphere soils were significantly inhibited by total/DTPA-extractable Pb and Zn concentrations (Table S2), suggesting that the active center of these enzymes might be inactivated by metal onions, and thereby preventing them from carrying out their hydrolytic action [51].This was consistent with the results found by other studies for slightly alkaline soil, which reported that β-glucosidase, protease and urease activities are sensitive bioindicators of soil degradation caused by HM contamination [52].However, alkaline phosphatase activity was less affected by heavy mental pollution level (Table S2), suggesting that in tolerant bacteria this enzyme was well protected and kept its functional activity in HM contaminated soils [53].On the other hand, the activities of β-glucosidase, protease and alkaline phosphatase in rhizosphere soil were significantly higher than those in bulk soil (Figure 1).The rhizosphere effect was likely due to physiological activities of plant roots under HM stress condition, which could cause the roots to release large amounts of enzymes into the rhizosphere soil [54].Meanwhile, microorganisms adhering plant roots were able to secrete various enzymes to the surrounding soil unceasingly.Therefore, higher enzyme activity may be attributed to larger bacterial copy number in rhizosphere soil (Figure 2), resulting in different enzyme activity profiles between the bulk and rhizosphere soils [55].
HM had toxic effects on the growth, morphology and metabolism of microorganism in soil, through functional disturbance, protein denaturation and destruction of the integrity of cell membranes.However, some metal-tolerant bacteria can still survive in soils heavily polluted with HMs, showing great potential for biotechnological and bioremediation applications [56].The relative abundance of the clones related to each of the major groups was calculated based on the number of clone sequences, and the results showed that the bacterial community structure was highly diverse and heterogeneous in bulk and rhizosphere soils at different HM pollution levels (Figure 4).The dominant bacterial classes found in both bulk and rhizosphere soils were Alphaproteobacteria, Gammaproteobacteria and Firmicutes.The results were consistent with Mendez et al. [57], who found that the soil bacterial communities were dominated by Proteobacteria (Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria) and Firmicutes in an abandoned semiarid lead-zinc mine tailing.The class of Alphaproteobacteria was also widely detected in soils highly polluted by toxic metals, such as Cu, Zn and Cd [58,59].Zhu et al. [60] investigated the bacterial composition of sediments contaminated with HMs in Xiangjiang River (Zhuzhou, Hunan Province, China), and found that the dominant families belonged to Alphaproteobacteria, Betaproteobacteria and Firmicutes.The current study confirmed that they were also the most dominant classes in semiarid lead-zinc mine region in the Northwest of China (Figure 4).Gammaproteobacteria also presented to be high percentage (26.3%for bulk soil, 26.7% for rhizosphere soil) in HM polluted soils.Gomez-Balderas et al. [61] reported a similar proportion of Gammaproteobacteria (30%) in soils polluted with Zn and Cd.Our results were also consistent with that a previous study which reported that Gammaproteobacteria isolated from the rhizosphere soil of willows had high resistance to Zn contamination [62].
In the present study, we confirmed that rhizosphere soil harbored much larger abundant and diversity of bacterial community compared with bulk soil (Figures 2, 4 and 5).The bacterial 16S rRNA copy number, Species richness (S), Margalef (d Ma ) were significantly higher in rhizosphere soil than these in bulk soil.However, no difference could be found in Simpson index (D) between bulk and rhizosphere soils.The phenomenon could be partly explained by that plant roots could release different types of organics (exudates, secretions, sloughed off cells) to rhizosphere soil, and provide a relatively stable environment for the growth of microorganisms [63].Besides, HM could be immobilized in soil when it formed complexes with root exudates, particularly mucilage components, producing much less or non-toxic HM [49].This explanation was supported by our results that the availability of Pb in rhizosphere soil was significantly lower compared with that in bulk soil, as shown in Figure S2.However, the detailed information about the components of root exudates from R. pseudoacacia required to be further studied to reveal the mechanisms of root exudates for immobilizing HMs in rhizosphere soil.
Although the large differences in bacterial abundance, composition and diversity between bulk and rhizosphere soils and among different study sites could be detected, redundancy analysis (RDA) indicated that HM concentration was not the only environmental parameters influencing bacterial abundance and composition in the present study (Figure 6).Numbers in brackets represented the percentage of variation of the data explained by each factor.RDA axis 1 explained 89.2% of the total variation and RDA axis 2 explains 7.8% of the total variation in bulk soil, while RDA axis 1 explained 67.3% of the total variation and RDA axis 2 explains 18.9% of the total variation in rhizosphere soil, respectively.The multivariate analysis showed that, apart from total/DTPA-extractable Pb concentrations, soil EC, pH and SOM content greatly influenced bacterial abundance and diversity in bulk soil (Figure 6a).Detailed information indicated that, in rhizosphere soil, the bacterial abundance and community was highly associated with total/DTPA-extractable Pb, soil pH, EC, TP content (Figure 6b).Although environmental factors influencing bacterial community structure are varied, several studies have found that HM pollutions (Pb and Cd) could be the major environmental variables affecting bacterial diversity [64].Soil pH has been proved to be one of the most important environmental factors determining the movement, availability, and toxicity HM, due to its strong effects on solubility of metals in soil [65].Our study provided the evidence that soil pH was negatively corrected with HM availability (p < 0.05) (Table S1).The current study also found that SOM had significantly negative correlations with soil total/DTPA-extractable Pb, Zn and Pb availability in bulk soil (p < 0.05), but no correlation between SOM and HM availability could be detected in rhizosphere soil (p > 0.05) (Table S1).Soil pH and SOM content altered bacterial community diversity and composition probably through influencing HM availability and changing the toxic effects on growth and metabolism of microorganisms in polluted soils.On the other hand, SOM was able to influence the function of the soil microbial community via affecting soil structural properties and availability of nutritional substrates [66].Soil EC has been used as a good indicator of soil stress level and is generally correlated with soil pH.In agreement with our findings, Kim et al. [67] suggested that the bacterial community in intensively cultivated greenhouse soils could be particularly affected by soil pH and EC.Soil pH had also shown one of the most important environmental factors impacting bacterial abundance and diversity in natural and agricultural ecosystems [68].However, it should be noted that the four plots in each HM pollution level might not be appropriate replicates from the view of large scale, and other environmental factors (vegetation, stage of plant growth, etc.) which was not mentioned in the study probably had great effects on soil property, bacterial abundance and community.Therefore, further research is still needed to improve the understanding of detailed correlations among other environmental factors and their specific effects on soil property, enzyme activity, bacterial abundance and community in bulk and rhizosphere soils along a gradient of HM contamination in forest ecosystem.The arrow length and direction correspond to the variance that can be explained by the environmental and response variables.The direction of an arrow reflects the extent to which the given factor is influenced by each RDA variable.The perpendicular distance between the abundance of bacterial taxa and environmental variable axes in the plot indicates their correlations.The smaller the distance, the stronger the correlation.EC, electrical conductivity; SOM, soil organic matter; SOC, soil organic carbon; TN, total nitrogen; AN, available nitrogen; TP, total phosphorus; AP, available phosphorus; TPb, total Pb; TZn, total Zn; DPb, DTPA-extractable Pb; DZn, DTPA-extractable Zn.

Conclusions
It has been widely accepted that the physical, chemical, and biological properties of rhizosphere soil were significantly different from those of the surrounding bulk soil.This study was conducted to illustrate the rhizosphere effect on soil properties, enzyme activities and bacterial community associated with Robinia pseudoacacia in HM contaminated soils, Northwest China.The results showed that the pH of rhizosphere soil was lower than that of bulk soil (p < 0.05).In contrast, soil organic matter (SOM), available nitrogen (AN) and phosphorus (AP) contents in rhizosphere soil were significantly higher than those in bulk soil at HM contaminated sites (p < 0.05).The activities of β-glucosidase, protease, urease and alkaline phosphatase in rhizosphere soil across all study sites increased by 47.5%, 64.1%, 52.9% and 103.8%, respectively.Bacterial 16S rRNA copy number in rhizosphere soil was significantly higher than that in bulk soil (p < 0.05), and was significantly inhibited by total/DTPA-extractable Pb concentrations (p < 0.01).Alphaproteobacteria, Gammaproteobacteria and Firmicutes were the most dominant groups of bacteria at different study sites.The bacterial diversity index of Species richness (S) and Margalef (dMa) were significantly higher in rhizosphere soil compared with that in bulk soil, although no difference could be found in Simpson index (D) between bulk and rhizosphere soils (p > 0.05).Redundancy analysis (RDA) showed that soil pH, EC, SOM and total/DTPA-extractable Pb concentrations were the most important variables affecting bacterial abundance, composition and diversity (p < 0.05).The information and insights obtained on the properties of soil nutrient content, enzyme activities, bacterial abundance and community in bulk and rhizosphere soils will enhance our understanding The direction of an arrow reflects the extent to which the given factor is influenced by each RDA variable.The perpendicular distance between the abundance of bacterial taxa and environmental variable axes in the plot indicates their correlations.The smaller the distance, the stronger the correlation.EC, electrical conductivity; SOM, soil organic matter; SOC, soil organic carbon; TN, total nitrogen; AN, available nitrogen; TP, total phosphorus; AP, available phosphorus; TPb, total Pb; TZn, total Zn; DPb, DTPA-extractable Pb; DZn, DTPA-extractable Zn.

Conclusions
It has been widely accepted that the physical, chemical, and biological properties of rhizosphere soil were significantly different from those of the surrounding bulk soil.This study was conducted to illustrate the rhizosphere effect on soil properties, enzyme activities and bacterial community associated with Robinia pseudoacacia in HM contaminated soils, Northwest China.The results showed that the pH of rhizosphere soil was lower than that of bulk soil (p < 0.05).In contrast, soil organic matter (SOM), available nitrogen (AN) and phosphorus (AP) contents in rhizosphere soil were significantly higher than those in bulk soil at HM contaminated sites (p < 0.05).The activities of β-glucosidase, protease, urease and alkaline phosphatase in rhizosphere soil across all study sites increased by 47.5%, 64.1%, 52.9% and 103.8%, respectively.Bacterial 16S rRNA copy number in rhizosphere soil was significantly higher than that in bulk soil (p < 0.05), and was significantly inhibited by total/DTPA-extractable Pb concentrations (p < 0.01).Alphaproteobacteria, Gammaproteobacteria and Firmicutes were the most dominant groups of bacteria at different study sites.The bacterial diversity index of Species richness (S) and Margalef (d Ma ) were significantly higher in rhizosphere soil compared with that in bulk soil, although no difference could be found in Simpson index (D) between bulk and rhizosphere soils (p > 0.05).Redundancy analysis (RDA) showed that soil pH, EC, SOM and total/DTPA-extractable Pb concentrations were the most important variables affecting bacterial abundance, composition and diversity (p < 0.05).The information and insights obtained on the properties of soil nutrient content, enzyme activities, bacterial abundance and community in bulk and rhizosphere soils will enhance our understanding of rhizosphere effect and improve remediation of HM contaminated soils in forest ecosystem.

Figure 1 .
Figure 1.Activities of β-glucosidase (a), protease (b), urease (c) and alkaline phosphatase (d) in bulk and rhizosphere soils at different HM pollution levels.Potential differences in parameters between bulk and rhizosphere soils were analyzed using t-test (** p < 0.01; ns, not significant).

Figure 2 .
Figure 2.Copy number of bacterial 16S rRNA gene in bulk and rhizosphere soils at different HM pollution levels.Potential differences in parameters between bulk and rhizosphere soils were analyzed using t-test (** p < 0.01; * p < 0.05; ns, not significant).

Figure 2 .
Figure 2.Copy number of bacterial 16S rRNA gene in bulk and rhizosphere soils at different HM pollution levels.Potential differences in parameters between bulk and rhizosphere soils were analyzed using t-test (** p < 0.01; ns, not significant).

Figure 3 .Figure 3 .
Figure 3. Neighbor-joining bacterial phylogenetic tree construction from partial 16S rRNA sequences.Clones from bulk soils are coded as BS; clones from rhizosphere soils are coded as RS; clones from four study sites are coded as S1, S2, S3 and S4, representing light, low, medium and high Figure 3. Neighbor-joining bacterial phylogenetic tree construction from partial 16S rRNA sequences.Clones from bulk soils are coded as BS; clones from rhizosphere soils are coded as RS; clones from four study sites are coded as S1, S2, S3 and S4, representing light, low, medium and high pollution levels.The access number for each division is shown on the tree.Bootstrap values (1000 replicates) higher than 50% are shown.Scale bar represents the 5% substitution percentage.

Figure 4 .
Figure 4. Bacterial composition in bulk soil (BS) and rhizosphere soil (RS) at different HM pollution levels.The relative abundance of each class based on 16S rRNA gene analysis is expressed in percentage.

Figure 4 .
Figure 4. Bacterial composition in bulk soil (BS) and rhizosphere soil (RS) at different HM pollution levels.The relative abundance of each class based on 16S rRNA gene analysis is expressed in percentage.

Figure 5 .
Figure 5.The bacterial diversity index of Species richness (S), Margalef (dMa) and Simpson index (D) in bulk and rhizosphere soils associated with R. pseudoacacia grown at different HM pollution levels.Potential differences in parameters between bulk and rhizosphere soils were analyzed using t-test (** p < 0.01; * p < 0.05; ns, not significant).

Figure 5 .
Figure 5.The bacterial diversity index of Species richness (S), Margalef (d Ma ) and Simpson index (D) in bulk and rhizosphere soils associated with R. pseudoacacia grown at different HM pollution levels.Potential differences in parameters between bulk and rhizosphere soils were analyzed using t-test (* p < 0.05; ns, not significant).

Forests
improve the understanding of detailed correlations among other environmental factors and their specific effects on soil property, enzyme activity, bacterial abundance and community in bulk and rhizosphere soils along a gradient of HM contamination in forest ecosystem.

Figure 6 .
Figure 6.Redundancy analysis (RDA) of the correlations between the bacteria taxa and soil physiochemical properties in bulk soil (a) and rhizosphere soil (b) at different HM pollution levels.The arrow length and direction correspond to the variance that can be explained by the environmental and response variables.The direction of an arrow reflects the extent to which the given factor is influenced by each RDA variable.The perpendicular distance between the abundance of bacterial taxa and environmental variable axes in the plot indicates their correlations.The smaller the distance, the stronger the correlation.EC, electrical conductivity; SOM, soil organic matter; SOC, soil organic carbon; TN, total nitrogen; AN, available nitrogen; TP, total phosphorus; AP, available phosphorus; TPb, total Pb; TZn, total Zn; DPb, DTPA-extractable Pb; DZn, DTPA-extractable Zn.

Figure 6 .
Figure 6.Redundancy analysis (RDA) of the correlations between the bacteria taxa and soil physiochemical properties in bulk soil (a) and rhizosphere soil (b) at different HM pollution levels.The arrow length and direction correspond to the variance that can be explained by the environmental and response variables.The direction of an arrow reflects the extent to which the given factor is influenced by each RDA variable.The perpendicular distance between the abundance of bacterial taxa and environmental variable axes in the plot indicates their correlations.The smaller the distance, the stronger the correlation.EC, electrical conductivity; SOM, soil organic matter; SOC, soil organic carbon; TN, total nitrogen; AN, available nitrogen; TP, total phosphorus; AP, available phosphorus; TPb, total Pb; TZn, total Zn; DPb, DTPA-extractable Pb; DZn, DTPA-extractable Zn.

Table 1 .
Chemical properties of bulk and rhizosphere soils at different study sites.