Abstract
Mining activity is a growing environmental concern as it contributes to heavy metals (HMs) pollution in agricultural soils. Microbial communities play an important role in the biogeochemical cycling of HMs and have the potential to be used as bioindicators. Arsenic (As) and lead (Pb) are the most hazardous HMs and are mainly originated from mining activities. However, spatial variation in microbial community in response to As and Pb contamination in paddy soils remains overlooked. In this study, the biological and chemical properties of sixteen soil samples from four sites (N01, N02, N03, and N04) near a Pb-Zn mining site at different As and Pb levels were examined. The results showed that soil pH, total As and Pb, bioavailable As and Pb, nitrate-nitrogen (NO3−-N) and ammonia-nitrogen (NH4+-N) were the most important factors in shaping the bacterial community structure. In addition, significant correlations between various bacterial genera and As and Pb concentrations were observed, indicating their potential roles in As and Pb biogeochemical cycling. These findings provide insights into the variation of paddy soil bacterial community in soils co-contaminated with different levels of As and Pb.
Introduction
Mining mineral ore has both significant positive and negative impacts on mining areas. Mining resources have been endorsed as a mean to stimulate socioeconomic development and reduce poverty (Li et al., 2014). However, mining activities generate waste, such as mine tailings, acid mine drainages and fly ashes that contain harmful trace heavy metals (HMs) (; Liu et al., 2019). HMs can be transported from mining sites to neighboring agricultural soils through surface runoff and atmospheric dust deposition (Mu et al., 2020). Consequently, most neighboring agricultural soils are contaminated with HMs, and the soil thereafter serves as a long-term sink for HMs (; ).
In China, mining activities are prevalent and this has led to serious soil HMs contamination nationwide (Zhong et al., 2020). According to the Chinese National Soil Contamination Survey Report, approximately 33.4% of soil around mining sites were contaminated with HMs (Mee and Mnr, 2014). Among them, arsenic (As) and lead (Pd) commonly co-exist in high concentrations due to mining and smelting of mineral ores (Wang et al., 2020a; ), and threaten soil, plant and human health. In paddy soil, excess As and Pb can reduce soil quality, stunt plant growth, and ultimately affect human health through the food chain (; Rai et al., 2019). For example, Wang et al. (2018) reported that the concentrations of As and Pb (up to 2749 and 624 μg kg−1, respectively) in rice grains grown in fields neighboring mining sites were significantly higher than the threshold value of 200 μg kg−1 (GB2762–2017). Therefore, the investigation of the HMs pollution in paddy soils surrounding mining areas can assess the risk of HMs to human health and provide the scientific basis for the remediation of HMs contaminated soils.
HMs accumulation in soils affects both the soil physiochemical properties and the microbial community structure over short- and long-term periods (Zhu et al., 2019). Zhang et al. (2015) reported that HMs in paddy soil may have adverse effects on the microbial activities and nutrients bioavailability. Pb contaminations have been shown to inhibit soil respiration and reduce microbial biomass (Liu et al., 2017). Similarly, dynamic changes were observed in the soil microbial community as a result of As contamination (Wang et al., 2016). Furthermore, in the soil-rice system, As contamination can markedly influence the iron (Fe) plaque-associated microbial community structure and metabolic potential (). Although many attempts were tried to remediate As and Pb pollution in soils neighboring mining sites, few have reported the combined effect of As and Pb contamination on the soil microbial community (Tian et al., 2020; She et al., 2021).
Molecular techniques such as high-throughput sequencing techniques have been used to evaluate the changes in microbial community due to HMs pollution (Ros et al., 2009; ). Changes in soil microbial community structure and diversity could serve as biomarkers to study the effect of HMs on soil health (Vinhal-Freitas et al., 2017; Li et al., 2020b; ). However, it should be noted that some studies have reported that HMs contamination had limited effect on the microbial community structure (Lorenz et al., 2006). Other soil factors such as pH, organic matter (OM), moisture, clay content and nutrient content can have a profound influence on the microbial activities, which might play more important role in reshaping microbial communities than the HMs concentrations (). However, studies on the difference between microbial communities in paddy soils with different distances from the mining site and the responses of microorganisms to HMs and soil properties are limited. HMs polluted paddy soil with various distances from the mining site may present different soil characteristics and HMs concentration gradients, then the most important factor in shaping the microbial communities remains unknown.
For this purpose, this present study was conducted with soil samples from four sites along a Pb-Zn mining site to 1) explore the influence of the different distance from the mining site on the physiochemical properties of As and Pb polluted paddy soils; 2) evaluate the combined effects of As and Pb pollution on soil bacterial community structure; 3) examine the correlativity between soil physiochemical properties and bacterial communities. Chemical analysis coupled with Miseq sequencing techniques were applied to achieve the goals. The observations and findings in this study could be used to develop more effective bioremediation strategies for the polluted paddy soil in mining area.
Materials and Methods
Study Area and Sampling
The study area is located near a Pb-Zn mining site (29°59′35.00″ N, 120°45′55.20″ E) in Shangyu city, Zhejiang province, China (Figure 1). The study area has a typical subtropical climate, densely surrounded by paddy field. A total of 16 top-soil samples (0–20 cm) were collected from the four sites referred to as N01, N02, N03, and N04. The four sites are all paddy fields and one growth cycle is performed over one year. Soil samples (four in total, approximately 300 g each) were collected from random points at each site and were mixed uniformly to form a composite sample. Four replicate composite samples were collected from each site (denoted as samples N01–1, N01–2, N01–3, N01–4, and so on). Soil samples were properly sealed and transported to the laboratory on ice. The collected soil samples were divided into two parts; one part was stored at −80°C for molecular analyses, and the other part was used for the physicochemical analysis.
FIGURE 1
Chemical Analysis
Soil samples from N01, N02, N03, and N04 site were analyzed for different physicochemical parameters including soil pH, OM, moisture content, dissolved organic carbon (DOC), available nitrate-nitrogen (NO3−-N) and ammonia-nitrogen (NH4+-N) concentrations, total As, Pb and other HMs concentrations, total sulfur (S) and SO42-. Soil pH and OM were measured following methods described in our pervious study (Zou et al., 2018). Moisture content was measured by drying fresh soil samples overnight in an oven at 105°C (O’Kelly, 2004). DOC was extracted with ultrapure water from fresh soil at a ratio of 1:2 (soil: water) and the concentration was analyzed with an automated N/C analyzer (Multi N/C 2100, Analytik Jena, Germany) (Xu et al., 2015). Available NH4+-N and NO3−-N concentrations were determined in 1 M KCl extracts (1:5, w/v) using the colorimetric methods (Stock, 1983). Soil SO42- was extracted by shaking freeze-dried soils with ultrapure water at a ratio 1:25 (soil:water) and analyzed by ion chromatography (Dionex, United States) (Xu et al., 2015). In order to determine the total Fe, S and other HMs, 0.2 g soil was digested using a concentrated acid mixture (5 ml of HNO3, 3 ml of HClO4, and 5 ml of HF) and was analyzed by inductively coupled plasma-optical emission spectrometry (ICP-OES, iCAP 6300, Thermo, United States). For total As analysis, 0.2 g air-dried soil was digested with 10 ml diluted aqua regia (aqua regia: water = 1:1, v/v) at 100°C for 2 h, and As concentration in digested samples were determined by atomic fluorescence spectrometry (AFS-9130, Beijing Jitian Instrument Company, China) (Zou et al., 2018). In order to determine the accuracy of the total metal analysis, a reference soil (GBW 07429, the National Research Center for standards in China) was used. The extract efficiencies ranged from 97.8 to 104.6% of the target values.
Bioavailable As and Pb Fraction in Soil
As and Pb were fractionated according to a 4-step modified European Community Bureau of Reference (BCR) sequential extraction procedure including F1: extracted with 0.11 M acetic acid (exchangeable fraction); F2: extracted with 0.5 M hydroxylammonium chloride (reducible fraction); F3: extracted with 8.8 M hydrogen peroxide followed by 1.0 M ammonium acetate (oxidizable fraction); F4: extracted with hydrofluoric acid followed by Aqua Regia (nitric acid + hydrochloric acid; 1:3 v/v) (residual fraction) (Rauret et al., 1999). The sum of exchangeable and reducible fractions was defined as bioavailable fractions according to a previous study (Roosa et al., 2014). The extracted solutions were analyzed by ICP-OES for different As and Pb fractions. The difference between the sum of each phase concentration and the total concentrations of As and Pb were compared to verify the accuracy of the sequential extraction procedure.
Microbial Biomass Carbon
Microbial biomass carbon (MBC) in soils was estimated by chloroform fumigation-extraction method (Vance et al., 1987; ). Briefly, 6 subsamples of 5 g moist soil were weighed into different glass centrifuge tubes and three of them were fumigated with ethanol-free chloroform. Both fumigated and non-fumigated samples were incubated for 24 h at 25°C in the dark. After incubation, all the samples were extracted with 25 ml of 0.5 M K2SO4 and filtered. Then, the fumigated samples were heated in a 100°C water bath for 60 min and lost moisture was replenished. Total dissolved carbon in the extracts was determined by an automated N/C analyzer (Multi N/C 2100, Analytik Jena, Germany). MBC was calculated according to Eq. 1 ():where Cf = total dissolved carbon extracted from fumigated soil, Cnf = total dissolved carbon extracted from nonfumigated soil and kEC = 0.45.
Soil DNA Extraction, PCR Amplification and Sequencing Analysis
Total DNA was extracted (four replicate samples) for each site using Fast DNA® Spin Kit for Soil (MP Biomedicals, Santa Ana, United States) according to the manufacturer’s protocol. Total genomic DNA was amplified using primers 926 F (5’ - AAA CTY AAA KGA ATT GAC GG - 3′) and 1392 R (5’ -ACG GGC GGTGTG TRC- 3′) that amplifies the V6-V8 region of the 16 S rRNA gene (Luo et al., 2016). The 16 S rRNA tag-encoded high-throughput sequencing was carried out on Illumina Miseq platform at the Tianke Technology Co., Ltd (Hangzhou, China). Sequencing reads were assigned to each sample according to the individual unique barcode and filtered to remove ambiguous reads. Operational taxonomic units (OTUs) were clustered with 97% similarity cutoff using UPARSE (V8.1.1861) (Wang et al., 2007). For each OTU, a representative sequence was selected and used to assign taxonomic composition by using the RDP classifier (V11.4) against the SILVA database (www.arb-silva.de). The 16S rRNA gene sequences reported in this paper have been deposited in the NCBI SRA under the accession No. PRJNA678872.
Statistical Analysis
Statistical analyses for the experimental data were performed using the SPSS® 19.0 (SPSS, United States) software. The significance of the treatment effects was assessed by analysis of variance (ANOVA), followed by comparisons between treatment means using the least significant difference (LSD) at p < 0.05. Alpha diversity index including Chao1, Shannon, Simpson, ACE and Coverage were calculated using QIIME software (V1.9.1) (). The similarity among microbial communities in different soil samples was determined using UniFrac analysis by QIIME software. Principal coordinate analysis (PCoA) plots of the OTU-based weighed UniFrac distances between all samples and heatmap of Spearman’s rank correlations coefficients between geochemical parameters and the top 35 bacterial genera were generated by R software (V3.5.1) (http://www.r-project.org/). Redundancy analysis (RDA) was performed based on abundant bacteria genera (relative abundance top 35) and selected physicochemical parameters by Canoco 5 (Biometrics Wageningen, The Netherlands) ().
Results and Discussion
Physicochemical Characterization of Soil Samples
The long term and continuous use of As and Pb contaminated paddy soil poses a serious health risks, since rice have a high tendency to accumulate As and Pb in its grains (). Paddy soil contaminated by As and Pb in mining areas has been widely documented (; Zheng et al., 2020). However, the impact of mining activities on paddy soil physicochemical and microbial properties at various distances away from the mine is not well understood. HMs contamination and physicochemical properties of paddy soil play an important role in microbial metabolism (Zhao et al., 2013). Physicochemical characteristics of collected soils from the four sites were presented in Table 1. For all four sites, the soil pH was close to neutral (6.56–7.33), and the maximum pH was observed in site N01. The NO3−-N and NH4+-N concentrations were significantly higher in site N01 as compared to the other studied sites (p < 0.05).
TABLE 1
| Samples | pH | OM (%) | Moisture (%) | Total Fe | DOC | NO3−-N | NH4+-N | SO42- | Total S | Cd | Cr | Cu | Zn | Ni |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| mmol g−1 dry soil | mg kg−1 dry soil | g kg−1 dry soil | mg kg−1 dry soil | |||||||||||
| N01 | 7.33 ± 0.08a | 4.95 ± 0.81ab | 28.2 ± 2.03c | 586 ± 27.2b | 28.0 ± 5.22c | 39.2 ± 9.38a | 45.8 ± 11.1a | 3.41 ± 0.55a | 0.39 ± 0.03a | 1.11 ± 0.11a | 42.0 ± 7.63a | 58.1 ± 9.07a | 168 ± 6.25a | 31.6 ± 0.90b |
| N02 | 6.56 ± 0.27b | 3.81 ± 0.54b | 35.6 ± 1.72bc | 456 ± 28.0c | 41.4 ± 7.25b | 0.95 ± 0.56b | 3.21 ± 0.38b | 1.56 ± 0.19b | 0.21 ± 0.03b | 0.50 ± 0.08c | 24.3 ± 6.99b | 13.3 ± 0.37b | 122 ± 17.9b | 22.2 ± 1.14c |
| N03 | 6.16 ± 0.14c | 6.03 ± 0.22a | 42.0 ± 8.73b | 574 ± 3.85b | 37.7 ± 10.3bc | 2.15 ± 0.90b | 2.16 ± 0.72b | 3.27 ± 0.46a | 0.33 ± 0.04a | 0.74 ± 0.05 b | 31.4 ± 9.19b | 15.3 ± 0.93b | 133 ± 10.7ab | 31.4 ± 0.41b |
| N04 | 6.40 ± 0.02bc | 6.11 ± 1.17a | 55.5 ± 4.90a | 661 ± 25.1a | 67.7 ± 12.3a | 2.73 ± 0.59b | 3.45 ± 1.29b | 1.43 ± 0.30b | 0.41 ± 0.16a | 0.47 ± 0.16c | 44.8 ± 10.2a | 15.1 ± 0.81b | 173 ± 38.7a | 36.9 ± 0.48a |
Physicochemical characteristics of the paddy soil samples from each site. Different letters indicate significant differences among the different sites (p < 0.05).
Abbreviations: DOC, dissolved organic carbon; OM, organic matter; Total S, Total sulfur.
As and Pb were found to be the dominant HMs in the selected sites (Table 1 and Table 2). The high As and Pb concentrations can be ascribed to intensive mining activities in the area. Previous studies have shown that As and Pb can easily be translocated from the mining site to neighboring farmlands via surface runoff and/or dust particles deposition (Wang et al., 2011; ; Li et al., 2020c). The total As and Pb concentrations in the studied area ranged between 44.4–388 mg kg−1 and 0.20–2.12 g kg−1, respectively. The results also showed that the As and Pb concentrations in all the sites exceeded the soil environmental quality risk control standard for soil contamination of agricultural land in China (GB15618–2018) (6.5 < pH ≤ 7.5, As 25 mg kg−1, Pb 120 mg kg−1) (Mee, 2018). For example, the As and Pb concentrations in N01 were 15 and 17 times of the soil environmental quality risk control standard (GB15618–2018), respectively. These results were consistent with the previous studies that suggested agricultural soils near mining sites were seriously polluted (; Zhen et al., 2019).
TABLE 2
| Samples | As (mg kg−1) | Pb (g kg−1) | ||
|---|---|---|---|---|
| Bioavailable Conc. | Total Conc. | Bioavailable Conc. | Total Conc. | |
| N01 | 113 ± 10.8a | 388 ± 10.6a | 1.24 ± 0.06a | 2.12 ± 0.03a |
| N02 | 13.7 ± 1.23b | 124 ± 10.4b | 0.36 ± 0.01 b | 0.54 ± 0.02b |
| N03 | 7.17 ± 2.96b | 74.8 ± 2.75c | 0.24 ± 0.01c | 0.34 ± 0.01c |
| N04 | 6.34 ± 2.13b | 44.4 ± 3.74d | 0.15 ± 0.01 d | 0.20 ± 0.02d |
Total and bioavailability concentrations of As and Pb in soil samples collected from four sites. Different letters indicate significant differences among the different sites (p < 0.05).
Abbreviations: Conc., concentration.
A clear As and Pb gradient was observed at all four sites. The As and Pb concentrations decreased with increasing distance from the mining site and vice versa (Table 2). Previous studies have reported that the HMs concentrations in the topsoil usually decrease with the increase in the distance away from the mine (; Shen et al., 2017; Li et al., 2020c). A significant correlation was observed between total As and total Pb (r = 0.970, p < 0.001) (Supplementary Table S1), which indicated that these metal (loid)s had similar pollution sources (Sun et al., 2010). Furthermore, HMs bioavailable concentrations in soils have been reported to have better correlation with their accumulation in plants, which would directly threaten human health via food chain (Zhang et al., 2018). The bioavailable As concentration in N01 was significantly higher than other sites (p < 0.05). However, bioavailable As concentrations in N02, N03, and N04 were much lower compared to their total As concentrations. On the contrary, the bioavailable Pb concentrations were about half of the total concentrations for all four sites. Bioavailable As (r = 0.715, p < 0.01) and Pb (r = 0.667, p < 0.01) were positively corrected with soil pH, while bioavailable As (r = −0.656, p < 0.01) and Pb (r = −0.585, p < 0.05) were negatively corrected with soil OM (Supplementary Table S1). This indicates that pH and OM might be the main factors influencing As and Pb bioavailability in our study. pH may affect the bioavailability of As and Pb by influencing the adsorption of As and Pb on soil minerals. For example, the decreasing positive surface charge of soil minerals with increasing pH facilitated the desorption of As (Masscheleyn et al., 1991). With the increase of OM, the toxicity of HMs could be reduced through chelation and sequestration of metal cations onto OM (Macoustra et al., 2019).
Changes in Soil Microbial Biomass Carbon
Soil microbial biomass carbon (MBC) has been widely used as a bioindicator to evaluate soil quality since the soil microbial community is more sensitive to changes in soil composition compared to soil physiochemical properties (). The MBC concentrations in N03 (1.39 mg g−1) and N04 (1.15 mg g−1) were significantly (p < 0.05) higher compared to those in N01 (0.52 mg g−1) and N02 (0.47 mg g−1) (Figure 2). The MBC concentrations increased with decreasing As and Pb concentrations. Moreover, significant negative correlations between MBC and As and Pb concentrations were observed (Supplementary Table S1). Our results are consistent with those of a previous study (Zhen et al., 2019) reporting that MBC was significantly reduced in soils collected near a mining site in Shaoguan city, China, particularly due to the presence of high HMs concentration. Other studies also noticed that MBC decreased significantly with the increasing metal concentrations (; Xu et al., 2019). This indicated that due to the presence of excess metals in the contaminated soils, the MBC was suppressed and the microbial population was impeded to multiply (Xu et al., 2019). This may have occurred because soil MBC is sensitive to HMs, since microorganisms require additional energy to deal with stress caused by HMs (Yang et al., 2006). This additional energy requirement may result in a reduction in the amount of substrate used for microbial growth (Xiao et al., 2017).
FIGURE 2
Changes in Microbial Community Structure
Soil microbial community structure and diversity have also been used as an index for soil HMs pollution due to their sensitivity to soil physiochemical changes (Song et al., 2018). A total of 882,322 valid reads were obtained from the sixteen paddy soil samples through Illumina MiSeq platform and the average sequence length was 463 bp. A total of 642,574 operational taxonomic units (OTUs) was revealed. The number of OTUs from N01 was significantly lower than those from N02, N03 and N04 (p < 0.05) (Supplementary Table S2). The alpha diversity indexes including ACE and Chao1 from N03 and N04 were significantly higher than the samples from N01 and N02 (p < 0.05) (Supplementary Table S2). Soil microbial diversity in moderately and heavily contaminated soil were observed to be reduced compared to those of healthy soil (). Similar to our study, Shen et al. (2019) also observed that As and Pb contamination could decrease microbial diversity. Furthermore, the PCoA [PC1 vs PC2 (explaining 86.74%)] analysis suggested an obvious separation of the bacterial communities in N01, N02 and N03 while the bacterial communities in N04 was close to N02 and N03 (Supplementary Figure S1).
The microbial community structures differed among the four sampled sites. The relative abundance of the sequences obtained from the different sites at the phylum level are shown in Figure 3. The top 10 dominant phyla in all samples were, Proteobacteria (28.9–39.8%), Acidobacteria (9.3–16.6%), Chloroflexi (6.2–18.7%), Nitrospirae (2.8–12.7%), Bacteroidetes (2.0–8.9%), Gemmatimonadetes (1.9–4.4%), Planctomycetes (2.5–4.0%), Verrucomicrobia (1.9–3.9%), Chlorobi (1.9–3.9%) and Actinobacteria (1.2–2.2%). These phyla have also been reported in previous studies as dominant bacterial groups in paddy soil (Li et al., 2017; Zhen et al., 2019; ). Pollution-related changes in the relative abundance of the phyla were observed. The relative abundances of Proteobacteria, Bacteroidetes, Planctomycetes, and Verrucomicrobia were significantly higher in N01 than those in N02, N03 and N04 (p < 0.05), while the relative abundances of Chloroflexi, Nitrospirae, and Chlorobi were significantly lower in N01 (p < 0.05). These results indicate that soil HMs pollution may suppress some microbial phyla while promote the growth of HMs resistant microbes (Xu et al., 2018). Thus, we can deduce from our results that the phyla Proteobacteria and Bacteroidetes are able to endure high As and Pb stresses. Previous studies have reported that the members of the phylum Proteobacteria are able to tolerate high HMs stress (Schneider et al., 2017; Li et al., 2020c), which was consistent with the results of our study. Moreover, other studies have shown that Proteobacteria dominate in the metal-tolerant cultures, because many species belonging to this phyla can reduce Fe and facilitate soil-adsorbed metals release which inadvertently increase metal stress (; ). Proteobacteria and Bacteroidetes have also been reported to harbor the majority of HMs resistance genes base on metagenome sequencing (Yan et al., 2020). Furthermore, an increase in Bacteroidetes relative abundance have been observed in Pb spiked soils compared to the unspiked soils ().
FIGURE 3
In order to further understand the microbial community response to geochemical conditions and HMs, the OTUs were further classified to the genus level (Figure 4). The four sampled sites were clustered into two distinct groups (one group included N01 and the other group included N02, N03, and N04) which was consistent with the results shown in the PCoA plot (Figure 4 and Supplementary Figure S1). Flavobacterium (2.0%), Haliangium (1.8%), Nitrospira (1.5%), Sphingomonas (1.4%) and Pseudomonas (1.3%) were the most abundant genera in N01, while Geobacter (1.4–1.9%), Thiobacillus (0.8–3.0%), and Bryobacter (0.6–1.5%) were the most abundant genera in N02, N03 and N04. The relative abundance of Haliangium, Nitrospira, Flavobacterium, Opitutus, Sphingomonas, Pseudomonas, Aquicella, Acidibacter, Blastocatella, and Sorangium were significantly higher in N01 than those in N02, N03, and N04 (p < 0.05). Previous studied have shown that the most of the bacterial strains from these genera are commonly involved in As resistance (; Titah et al., 2018; Li et al., 2020c). Haliangium was observed to flourish in HMs stressed environments (Marti et al., 2017). Similarly, the abundances of Flavobacterium were observed to increase under high Pb stress (Wang et al., 2020b). Flavobacterium was also isolated from As contaminated soil as an As-resistant genus and it was able to reduce arsenate (As(V)) under erobic conditions (Wang and Mulligan, 2006). Besides, Flavobacterium related groups are able to methylate As as a detoxification strategy (). Nitrospira are nitrite-oxidizing bacteria that can provide nitrate which takes part in the oxidation of arsenite (As(III)) (Yu et al., 2020). Studies have found that in HMs polluted soil, Sphingomonas was the most abundant genus and had a potential to remove Pb from seawater (Chen, 2012; ). Pseudomonas has been reported to be involved in As redox processes (). Moreover, the As(III) efflux pump protein gene arsB and As transporter protein gene acr3P was also detected within the genus of Pseudomonas (Sanyal et al., 2016). Acidibacter and Blastocatella increased in HMs polluted soil and their significant positive correlation with Pb or As was also observed by . The relative abundance of Geobacter, Thiobacillus, Sideroxydans, Anaerolinea, Anaeromyxobacter, Leptolinea, Spirochaeta_2, Desulfobacca, Syntrophorhabdus, and Ignavibacterium were significantly higher in N02, N03, N04 than those in N01 (p < 0.05). These results indicated that these microorganisms were relatively sensitive to As and Pb toxicity and the abundances were abruptly declined with the increase of HMs concentrations (). Geobacter populations were reported as Fe and sulfate reducing bacteria, while Anaeromyxobacter were Fe reducing bacteria and Desulfobacca were sulfate reducing bacteria (Sun et al., 2015a). Besides, Thiobacillus were found to be involved in Fe and S oxidation and Thiobacillus can oxidize elemental sulfur (Wang et al., 2021). Sideroxydans were involved in Fe oxidation (). Most of these sensitive microorganisms were Fe- and S-related bacteria indicating their important roles in lightly contaminated sites. We speculated that As and Pb contamination in paddy soil might have an inhibition effect on Fe and S cycling in paddy soil.
FIGURE 4
Correlation Between Environmental Parameters and Microbial Community
In microbial ecology, understanding the relationship between environmental variables and microbial community structure is a key goal (Liang et al., 2011). Spearman’s rank correlations coefficients between geochemical parameters and bacterial abundances were determined (Figure 5). The relative abundance of Acidibacter positively correlated with total As, total Pb, bioavailable Pb (p < 0.05) and bioavailable As (p < 0.001) concentrations in soil. Similar results have also been reported by . Nitrospira, Aquicella and Sorangium were also positively correlated with total As, bioavailable As, total Pb and bioavailable Pb concentrations. Other genera like Sphingomonas, Pseudomonas and Massilia also had positive correlations with total or bioavailable As or Pb. These results suggested that As and Pb contamination may facilitate the growth of these bacterial groups. Thioalkalispira (p < 0.001), Desulfobacca (p < 0.01) and Spirochaeta_2 (p < 0.01) were negatively correlated with total As, bioavailable As, total Pb and bioavailable Pb concentrations. Among them, Nitrospira, Aquicella and Sorangium were positively correlated with different forms of As and Pb, indicating their high HMs tolerance while Thioalkalispira, Desulfobacca and Spirochaeta_2 were sensitive to As and Pb pollution. pH, moisture and NH4+-N had significant correlations with most of the genera. Redundancy analysis (RDA) was used to reveal the possible linkages between the microbial communities and geochemical parameters (Figure 6). Although the pollution source of the soil was from a Pb-Zn mining site, no significant correlation between the microbial communities and the Zn concentration could be found. It might be due to the low Zn concentrations in soil. Besides, it has been reported that high levels of Zn had no effect on the species diversity compared to unpolluted soil in mining areas (Xu et al., 2017). Total As, bioavailable As, total Pb and bioavailable Pb concentrations along with pH, NH4+-N, NO3−-N, moisture, SO42- and MBC were the main factors influencing bacterial community compositions (Figure 6). Total As, bioavailable As, total Pb and bioavailable Pb concentrations along with pH, NH4+-N and NO3−-N positively correlated with microbial communities from N01 while negatively correlated with microbial communities from N02, N03, and N04. Soil moisture positively correlated with microbial communities from N02, N03, and N04. It is widely accepted that pH has a significant effect on the overall diversity and composition of microbial communities (). Any significant deviation in pH could impose stress on single-celled organisms because the intracellular pH of most microorganisms are usually within 1 pH unit of neutral (). Soil pH could also influence other environmental parameters (e.g. nutrient availability or cationic metal solubility), that plays an important role in microbial survival (Sun et al., 2015b). Previous studies also reported that pH had a greater impact on bacterial community composition in HMs contaminated agricultural soils when compared with HMs in soil (; ). Two types of nitrogen (NH4+-N and NO3−-N) were significantly correlated with soil genera in this study, suggesting their important role in shaping the soil microbial community (). reported that NH4+-N and NO3−-N showed positive correlations with the shifts in microbial community structure, which was consistent with our results. Paddy soils could provide good habitats for nitrogen cycling, which in turn have the potential to shape the microbial communities (Sun et al., 2015a). Additionally, moisture is a major environmental factor that determines the redox conditions of soils and the metabolism of microbes (). For example, wet-dry transition in soil would result in a microbial community shift of anaerobic to erobic microbes (Pett-Ridge and Firestone., 2005). This is especially true for paddy soil which are intermittently flooded during rice production.
FIGURE 5
FIGURE 6
Conclusion
In this study, we used the high throughput sequencing to examine the effects of physicochemical parameters on the microbial community variation in HMs contaminated soils along an As and Pb gradient. Our results showed that high amounts of HMs were presented near the mining site, and their concentrations decreased with increasing distance from the mining site. As and Pb contamination have a profound impact on the microbial community structure of paddy soils. Different forms of As and Pb (total and bioavailable) and soil physicochemical components (pH, NH4+-N, NO3−-N, moisture, and MBC) had significant influences on the bacterial community compositions. Some genera including Haliangium, Nitrospira, and Flavobacterium were positively correlated with different forms of As and Pb, indicating their high tolerance and their important role in As and Pb cycling while other genera including Geobacter, Thiobacillus, and Sideroxydans were sensitive to HMs contamination. Our finding provided important information that should be taken into consideration when developing bioremediation strategies for the polluted paddy soil in mining area.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
XT and JN performed the study conceptualization and supervision. LZ performed the literature search and wrote the original draft. YL, YD, MK, WG, and LY carried out the writing, reviewing and JS and JX performed the study supervision and validation editing. All authors read and approved the final manuscript.
Acknowledgments
The National Key Technology R&D Program (2018YFD0800202) supported this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fenvs.2021.630668/full#supplementary-material.
References
1
AazamiJ.Esmaili-SariA.AbdoliSohrabiA. H.Van den BrinkP. J. (2015). Monitoring and Assessment of Water Health Quality in the Tajan River, Iran Using Physicochemical, Fish and Macroinvertebrates Indices. J. Environ. Health Sci. Engineer.13, 29. 10.1186/s40201-015-0186-y
2
AndersonJ. A. H.HooperM. J.ZakJ. C.CoxS. B. (2009). Characterization of the Structural and Functional Diversity of Indigenous Soil Microbial Communities in Smelter-Impacted and Nonimpacted Soils. Environ. Toxicol. Chem.28, 534–541. 10.1897/08-281.1
3
BiswasJ. K.WarkeM.DattaR.SarkarD. (2020). Is Arsenic in Rice a Major Human Health Concern?Curr. Pollut. Rep.6, 37–42. 10.1007/s40726-020-00148-2
4
BlotheM.RodenE. E. (2009). Composition and Activity of an Autotrophic Fe(II)-oxidizing, Nitrate-Reducing Enrichment Culture. Appl. Environ. Microbiol.75, 6937–6940. 10.1128/AEM.01742-09
5
BrookesP. C. (1995). The Use of Microbial Parameters in Monitoring Soil Pollution by Heavy Metals. Biol. Fertil. Soils.19, 269–279. 10.1007/bf00336094
6
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al (2010). QIIME Allows Analysis of High-Throughput Community Sequencing Data. Nat. Methods7, 335–336. 10.1038/nmeth.f.303
7
CavalcaL.CorsiniA.ZaccheoP.AndreoniV.MuyzerG. (2013). Microbial Transformations of Arsenic: Perspectives for Biological Removal of Arsenic from Water. Future Microbiol.8, 753–768. 10.2217/fmb.13.38
8
CavaniL.ManiciL. M.CaputoF.PeruzziE.CiavattaC. (2016). Ecological Restoration of a Copper Polluted Vineyard: Long-Term Impact of Farmland Abandonment on Soil Bio-Chemical Properties and Microbial Communities. J. Environ. Manage.182, 37–47. 10.1016/j.jenvman.2016.07.050
9
CesareA. D.PjevacP.EckertE.CurkowN.ŠparicaM. M.CornoG.et al (2020). The Role of Metal Contamination in Shaping Microbial Communities in Heavily Polluted Marine Sediments. Environ. Pollut.265, 114823. 10.1016/j.envpol.2020.114823
10
ChenM.LuG.WuJ.YangC.NiuX.TaoX.et al (2018). Migration and Fate of Metallic Elements in a Waste Mud Impoundment and Affected River Downstream: A Case Study in Dabaoshan Mine, South China. Ecotoxicology Environ. Saf.164, 474–483. 10.1016/j.ecoenv.2018.08.063
11
ChengW.LeiS.BianZ.ZhaoY.LiY.GanY. (2020). Geographic Distribution of Heavy Metals and Identification of Their Sources in Soils Near Large, Open-Pit Coal Mines Using Positive Matrix Factorization. J. Hazard. Mater.387, 121666. 10.1016/j.jhazmat.2019.121666
12
DengS.KeT.LiL.CaiS.ZhouY.LiuY.et al (2018). Impacts of Environmental Factors on the Whole Microbial Communities in the Rhizosphere of a Metal-Tolerant Plant: Elsholtzia Haichowensis Sun. Environ. Pollut.237, 1088–1097. 10.1016/j.envpol.2017.11.037
13
DuF.YangZ.LiuP.WangL. (2018). Accumulation, Translocation, and Assessment of Heavy Metals in the Soil-Rice Systems Near a Mine-Impacted Region. Environ. Sci. Pollut. Res.25, 32221–32230. 10.1007/s11356-018-3184-7
14
Fernandez-MaciasJ. C.Gonzalez-MilleD. J.Garcia-ArreolaM. E.Cruz-SantiagoO.Rivero-PerezN. E.Perez-VazquezF.et al (2020). Integrated Probabilistic Risk Assessment in Sites Contaminated with Arsenic and Lead by Long-Term Mining Liabilities in San Luis Potosi, Mexico. Ecotox. Environ. Safe.197, 8. 10.1016/j.ecoenv.2020.110568
15
FeyziH.ChoromM.BagheriG. (2020). Urease Activity and Microbial Biomass of Carbon in Hydrocarbon Contaminated Soils. A Case Study of Cheshmeh-Khosh Oil Field. Iran. Ecotox. Environ. Safe.199, 11. 10.1016/j.ecoenv.2020.110664
16
FiererN.JacksonR. B. (2006). The Diversity and Biogeography of Soil Bacterial Communities. Proc. Natl. Acad. Sci.103, 626–631. 10.1073/pnas.0507535103
17
FrancisA. J.DodgeC. J. (1990). Anaerobic Microbial Remobilization of Toxic Metals Coprecipitated with Iron Oxide. Environ. Sci. Technol.24, 373–378. 10.1021/es00073a013
18
GillerK. E.WitterE.McgrathS. P. (1998). Toxicity of Heavy Metals to Microorganisms and Microbial Processes in Agricultural Soils: A Review. Soil Biol. Biochem.30, 1389–1414. 10.1016/s0038-0717(97)00270-8
19
GołębiewskiM.Deja-SikoraE.CichoszM.TretynA.WróbelB. (2014). 16S rDNA Pyrosequencing Analysis of Bacterial Community in Heavy Metals Polluted Soils. Microb. Ecol.67, 635–647. 10.1007/s00248-013-0344-7
20
GuY.Van NostrandJ. D.WuL.HeZ.QinY.ZhaoF. J.et al (2017). Bacterial Community and Arsenic Functional Genes Diversity in Arsenic Contaminated Soils from Different Geographic Locations. Plos. One.12. e0189656. 10.1371/journal.pone.0189656
21
GuoH.NasirM.LvJ.DaiY.GaoJ. (2017). Understanding the Variation of Microbial Community in Heavy Metals Contaminated Soil Using High Throughput Sequencing. Ecotoxicology Environ. Saf.144, 300–306. 10.1016/j.ecoenv.2017.06.048
22
HahnJ.MannB.BangeU.KimmelM. (2019). Horizon-specific Effects of Heavy Metal Mobility on Nitrogen Binding Forms in Forest Soils Near a Historic Smelter (Germany). Geoderma355, 8. 10.1016/j.geoderma.2019.113895
23
HonmaT.OhbaH.Kaneko-KadokuraA.MakinoT.NakamuraK.KatouH. (2016). Optimal Soil Eh, pH, and Water Management for Simultaneously Minimizing Arsenic and Cadmium Concentrations in Rice Grains. Environ. Sci. Technol.50, 4178–4185. 10.1021/acs.est.5b05424
24
HonschoppS.BrunkenN.NehrkornA.BreunigH. J. (1996). Isolation and Characterization of a New Arsenic Methylating Bacterium from Soil. Microbiol. Res.151, 37–41. 10.1016/s0944-5013(96)80053-x
25
HouX. L.HanH.TigabuM.LiQ. Y.LiZ. X.ZhuC. L.et al (2020). Lead Contamination Alters Enzyme Activities and Microbial Composition in the Rhizosphere Soil of the Hyperaccumulator Pogonatherum Crinitum. Ecotoxicol Environ. Saf.207, 111308. 10.1016/j.ecoenv.2020.111308
26
HuM.SunW.KruminsV.LiF. (2019). Arsenic Contamination Influences Microbial Community Structure and Putative Arsenic Metabolism Gene Abundance in Iron Plaque on Paddy Rice Root. Sci. Total Environ.649, 405–412. 10.1016/j.scitotenv.2018.08.388
27
JaenickeS.AnderC.BekelT.BisdorfR.DrögeM.GartemannK. H.et al (2011). Comparative and Joint Analysis of Two Metagenomic Datasets from a Biogas Fermenter Obtained by 454-pyrosequencing. Plos. One.6. e14519. 10.1371/journal.pone.0014519
28
JiangB.AdebayoA.JiaJ.XingY.DengS.GuoL.et al (2019). Impacts of Heavy Metals and Soil Properties at a Nigerian E-Waste Site on Soil Microbial Community. J. Hazard. Mater.362, 187–195. 10.1016/j.jhazmat.2018.08.060
29
JiangJ.PanC. H.XiaoA. P.YangX. A.ZhangG. M. (2017). Isolation, Identification, and Environmental Adaptability of Heavy-Metal-Resistant Bacteria from Ramie Rhizosphere Soil Around Mine Refinery. 3 Biotech.7, 6. 10.1007/s13205-017-0603-2
30
JiangL. F.SongM. K.YangL.ZhangD. Y.SunY. T.ShenZ. G.et al (2016). Exploring the Influence of Environmental Factors on Bacterial Communities within the Rhizosphere of the Cu-Tolerant Plant. Elsholtzia Splendens. Sci. Rep.6, 10. 10.1038/srep36302
31
KashyapR.AhmadM.UniyalS. K.VermaK. S. (2019). Dietary Consumption of Metal(loid)s-Contaminated Rice Grown in Croplands Around Industrial Sectors: a Human Health Risk Perspective. Int. J. Environ. Sci. Technol.16, 8505–8516. 10.1007/s13762-019-02258-x
32
KhanS.El-Latif HeshamA.QiaoM.RehmanS.HeJ.-Z. (2010). Effects of Cd and Pb on Soil Microbial Community Structure and Activities. Environ. Sci. Pollut. Res.17, 288–296. 10.1007/s11356-009-0134-4
33
LeungH. M.Duzgoren-AydinN. S.AuC. K.KrupanidhiS.FungK. Y.CheungK. C.et al (2017). Monitoring and Assessment of Heavy Metal Contamination in a Constructed Wetland in Shaoguan (Guangdong Province, China): Bioaccumulation of Pb, Zn, Cu and Cd in Aquatic and Terrestrial Components. Environ. Sci. Pollut. Res.24, 9079–9088. 10.1007/s11356-016-6756-4
34
LiP. F.LiuM.MaX. Y.WuM.JiangC. Y.LiuK.et al (2020a). Responses of Microbial Communities to a Gradient of Pig Manure Amendment in Red Paddy Soils. Sci. Total. Environ.705, 8. 10.1016/j.scitotenv.2019.135884
35
LiP.LinC.ChengH.DuanX.LeiK. (2015). Contamination and Health Risks of Soil Heavy Metals Around a Lead/zinc Smelter in Southwestern China. Ecotoxicology Environ. Saf.113, 391–399. 10.1016/j.ecoenv.2014.12.025
36
LiS.WuJ.HuoY.ZhaoX.XueL. (2020b). Profiling Multiple Heavy Metal Contamination and Bacterial Communities Surrounding an Iron Tailing Pond in Northwest China. Sci. Total. Environ.752, 141827. 10.1016/j.scitotenv.2020.141827
37
LiS.ZhaoB.JinM.HuL.ZhongH.HeZ. (2020c). A Comprehensive Survey on the Horizontal and Vertical Distribution of Heavy Metals and Microorganisms in Soils of a Pb/Zn Smelter. J. Hazard. Mater.400. 123255. 10.1016/j.jhazmat.2020.123255
38
LiX.MengD.LiJ.YinH.LiuH.LiuX.et al (2017). Response of Soil Microbial Communities and Microbial Interactions to Long-Term Heavy Metal Contamination. Environ. Pollut.231, 908–917. 10.1016/j.envpol.2017.08.057
39
LiZ.MaZ.Van Der KuijpT. J.YuanZ.HuangL. (2014). A Review of Soil Heavy Metal Pollution from Mines in China: Pollution and Health Risk Assessment. Sci. Total Environ.468-469, 843–853. 10.1016/j.scitotenv.2013.08.090
40
LiangY.Van NostrandJ. D.DengY.HeZ.WuL.ZhangX.et al (2011). Functional Gene Diversity of Soil Microbial Communities from Five Oil-Contaminated Fields in China. Isme. J.5, 403–413. 10.1038/ismej.2010.142
41
LiuG.LingS.ZhanX.LinZ.ZhangW.LinK. (2017). Interaction Effects and Mechanism of Pb Pollution and Soil Microorganism in the Presence of Earthworm. Chemosphere173, 227–234. 10.1016/j.chemosphere.2017.01.022
42
LiuJ.YinM.ZhangW.TsangD. C. W.WeiX.ZhouY.et al (2019). Response of Microbial Communities and Interactions to Thallium in Contaminated Sediments Near a Pyrite Mining Area. Environ. Pollut.248, 916–928. 10.1016/j.envpol.2019.02.089
43
LorenzN.HintemannT.KramarewaT.KatayamaA.YasutaT.MarschnerP.et al (2006). Response of Microbial Activity and Microbial Community Composition in Soils to Long-Term Arsenic and Cadmium Exposure. Soil Biol. Biochem.38, 1430–1437. 10.1016/j.soilbio.2005.10.020
44
LuoF.DevineC. E.EdwardsE. A. (2016). Cultivating Microbial Dark Matter in Benzene-Degrading Methanogenic Consortia. Environ. Microbiol.18, 2923–2936. 10.1111/1462-2920.13121
45
MacoustraG.HollandA.StauberJ.JolleyD. F. (2019). Effect of Various Natural Dissolved Organic Carbon on Copper Lability and Toxicity to the Tropical Freshwater Microalga Chlorella Sp. Environ. Sci. Technol.53, 2768–2777. 10.1021/acs.est.8b04737
46
MartiR.Becouze-LareureC.RibunS.MarjoletL.SouibguiC. B.AubinJ. B.et al (2017). Bacteriome Genetic Structures of Urban Deposits Are Indicative of Their Origin and Impacted by Chemical Pollutants. Sci. Rep.7, 14. 10.1038/s41598-017-13594-8
47
MasscheleynP. H.DelauneR. D.PatrickW. H. (1991). Effect of Redox Potential and pH on Arsenic Speciation and Solubility in a Contaminated Soil. Environ. Sci. Technol.25, 1414–1419. 10.1021/es00020a008
48
MeeC. (2018). Soil Environmental Quality Risk Control Standard for Soil Contamination of Agricultural Land.
49
MeeC.MnrC. (2014). The Report on the National Soil Contamination Survey.
50
MuX.WangZ.LiuL.GuoX.GuC.XuH.et al (2020). Multiple Exposure Pathways of First-Year University Students to Heavy Metals in China: Serum Sampling and Atmospheric Modeling. Sci. Total Environ.746, 141405. 10.1016/j.scitotenv.2020.141405
51
O’KellyB. C. (2004). Accurate Determination of Moisture Content of Organic Soils Using the Oven Drying Method. Dry. Technol.22, 1767–1776. 10.1081/DRT-200025642
52
Pett-RidgeJ.FirestoneM. K. (2005). Redox Fluctuation Structures Microbial Communities in a Wet Tropical Soil. Aem71, 6998–7007. 10.1128/aem.71.11.6998-7007.2005
53
RaiP. K.LeeS. S.ZhangM.TsangY. F.KimK.-H. (2019). Heavy Metals in Food Crops: Health Risks, Fate, Mechanisms, and Management. Environ. Int.125, 365–385. 10.1016/j.envint.2019.01.067
54
RauretG.López-SánchezJ. F.SahuquilloA.RubioR.DavidsonC.UreA.et al (1999). Improvement of the BCR Three Step Sequential Extraction Procedure Prior to the Certification of New Sediment and Soil Reference Materials. J. Environ. Monitor.1, 57–61. 10.1039/a807854h
55
RoosaS.WattiezR.PrygielE.LesvenL.BillonG.GillanD. C. (2014). Bacterial Metal Resistance Genes and Metal Bioavailability in Contaminated Sediments. Environ. Pollut.189, 143–151. 10.1016/j.envpol.2014.02.031
56
RosM.PascualJ. A.MorenoJ. L.HernandezM. T.GarciaC. (2009). Evaluation of Microbial Community Activity, Abundance and Structure in a Semiarid Soil under Cadmium Pollution at Laboratory Level. Water Air Soil Pollut.203, 229–242. 10.1007/s11270-009-0006-z
57
SanyalS. K.MouT. J.ChakrabartyR. P.HoqueS.HossainM. A.SultanaM. (2016). Diversity of Arsenite Oxidase Gene and Arsenotrophic Bacteria in Arsenic Affected Bangladesh Soils. Amb. Express.6, 11. 10.1186/s13568-016-0193-0
58
SchneiderA. R.GommeauxM.DuclercqJ.FaninN.ConreuxA.AlahmadA.et al (2017). Response of Bacterial Communities to Pb Smelter Pollution in Contrasting Soils. Sci. Total Environ.605-606, 436–444. 10.1016/j.scitotenv.2017.06.159
59
SheJ. Y.WangJ.WeiX. D.ZhangQ.XieZ. Y.BeiyuanJ. Z.et al (2021). Survival Strategies and Dominant Phylotypes of Maize-Rhizosphere Microorganisms under Metal(loid)s Contamination. Sci. Total. Environ.774, 145143. 10.1016/j.scitotenv.2021.145143
60
ShenF.LiaoR.AliA.MaharA.GuoD.LiR.et al (2017). Spatial Distribution and Risk Assessment of Heavy Metals in Soil Near a Pb/Zn Smelter in Feng County, China. Ecotoxicology Environ. Saf.139, 254–262. 10.1016/j.ecoenv.2017.01.044
61
ShenT.LiuL.LiY.WangQ.DaiJ.WangR. (2019). Long-term Effects of Untreated Wastewater on Soil Bacterial Communities. Sci. Total Environ.646, 940–950. 10.1016/j.scitotenv.2018.07.223
62
SongJ.ShenQ.WangL.QiuG.ShiJ.XuJ.et al (2018). Effects of Cd, Cu, Zn and Their Combined Action on Microbial Biomass and Bacterial Community Structure. Environ. Pollut.243, 510–518. 10.1016/j.envpol.2018.09.011
63
StockW. D. (1983). An Evaluation of Some Manual Colorimetric Methods for the Determination of Inorganic Nitrogen in Soil Extracts. Commun. Soil Sci. Plant Anal.14, 925–936. 10.1080/00103628309367420
64
SunM.XiaoT.NingZ.XiaoE.SunW. (2015a). Microbial Community Analysis in Rice Paddy Soils Irrigated by Acid Mine Drainage Contaminated Water. Appl. Microbiol. Biotechnol.99, 2911–2922. 10.1007/s00253-014-6194-5
65
SunW.XiaoT.SunM.DongY.NingZ.XiaoE.et al (2015b). Diversity of the Sediment Microbial Community in the Aha Watershed (Southwest China) in Response to Acid Mine Drainage Pollution Gradients. Appl. Environ. Microbiol.81, 4874–4884. 10.1128/aem.00935-15
66
SunY.ZhouQ.XieX.LiuR. (2010). Spatial, Sources and Risk Assessment of Heavy Metal Contamination of Urban Soils in Typical Regions of Shenyang, China. J. Hazard. Mater.174, 455–462. 10.1016/j.jhazmat.2009.09.074
67
TianH. J.FengJ.ZhangL. M.HeJ. Z.LiuR. Y. (2020). Ecological Drivers of Methanotrophic Communities in Paddy Soils Around Mercury Mining Areas. Sci. Total. Environ.721. 137760. 10.1016/j.scitotenv.2020.137760
68
TitahH. S.AbdullahS. R. S.IdrisM.AnuarN.BasriH.MukhlisinM.et al (2018). Arsenic Resistance and Biosorption by Isolated Rhizobacteria from the Roots of Ludwigia Octovalvis. Int. J. Microbiol.2018, 1–10. 10.1155/2018/3101498
69
VanceE. D.BrookesP. C.JenkinsonD. S. (1987). An Extraction Method for Measuring Soil Microbial Biomass C. Soil Biol. Biochem.19, 703–707. 10.1016/0038-0717(87)90052-6
70
Vinhal-FreitasI. C.CorrêaG. F.WendlingB.BobuľskáL.FerreiraA. S. (2017). Soil Textural Class Plays a Major Role in Evaluating the Effects of Land Use on Soil Quality Indicators. Ecol. Indicators.74, 182–190. 10.1016/j.ecolind.2016.11.020
71
WangL.ChenL.GuoB.TsangD. C. W.HuangL.OkY. S.et al (2020a). Red Mud-Enhanced Magnesium Phosphate Cement for Remediation of Pb and as Contaminated Soil. J. Hazard. Mater.400. 10.1016/j.jhazmat.2020.123317
72
WangL.MaL.YangZ. (2018). Spatial Variation and Risk Assessment of Heavy Metals in Paddy Rice from Hunan Province, Southern China. Int. J. Environ. Sci. Technol.15, 1561–1572. 10.1007/s13762-017-1504-y
73
WangQ.GarrityG. M.TiedjeJ. M.ColeJ. R. (2007). Naïve Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy. AEM73, 5261–5267. 10.1128/aem.00062-07
74
WangR.HouD.ChenJ.LiJ.FuY.WangS.et al (2020b). Distinct Rhizobacterial Functional Assemblies Assist Two Sedum Alfredii Ecotypes to Adopt Different Survival Strategies under Lead Stress. Environ. Int.143, 105912. 10.1016/j.envint.2020.105912
75
WangS.MulliganC. N. (2006). Natural Attenuation Processes for Remediation of Arsenic Contaminated Soils and Groundwater. J. Hazard. Mater.138, 459–470. 10.1016/j.jhazmat.2006.09.048
76
WangY.-J.ChenZ.LiuP.-P.SunG.-X.DingL.-J.ZhuY.-G. (2016). Arsenic Modulates the Composition of Anode-Respiring Bacterial Community during Dry-Wet Cycles in Paddy Soils. J. Soils Sediments.16, 1745–1753. 10.1007/s11368-016-1369-6
77
WangY.ZhaoQ. L.HuY.DuX.GeWei.LiuW. J. (2011). Survey and Contamination Assessment of Heavy Metals in Soil and Plants Around the Pb/Zn Mine in Shangyu, Zhejiang Province (In Chinese). Environ. Chem.30, 1354–1360.
78
WangZ.ZhangB.HeC.ShiJ.WuM.GuoJ. (2021). Sulfur-based Mixotrophic Vanadium (V) Bio-Reduction towards Lower Organic Requirement and Sulfate Accumulation. Water Res.189, 116655. 10.1016/j.watres.2020.116655
79
XiaoX.-Y.WangM.-W.ZhuH.-W.GuoZ.-H.HanX.-Q.ZengP. (2017). Response of Soil Microbial Activities and Microbial Community Structure to Vanadium Stress. Ecotoxicology Environ. Saf.142, 200–206. 10.1016/j.ecoenv.2017.03.047
80
XuX.ZhangZ.HuS.RuanZ.JiangJ.ChenC.et al (2017). Response of Soil Bacterial Communities to Lead and Zinc Pollution Revealed by Illumina MiSeq Sequencing Investigation. Environ. Sci. Pollut. Res.24, 666–675. 10.1007/s11356-016-7826-3
81
XuY.HeY.ZhangQ.XuJ.CrowleyD. (2015). Coupling between Pentachlorophenol Dechlorination and Soil Redox as Revealed by Stable Carbon Isotope, Microbial Community Structure, and Biogeochemical Data. Environ. Sci. Technol.49, 5425–5433. 10.1021/es505040c
82
XuY.SeshadriB.BolanN.SarkarB.OkY. S.ZhangW.et al (2019). Microbial Functional Diversity and Carbon Use Feedback in Soils as Affected by Heavy Metals. Environ. Int.125, 478–488. 10.1016/j.envint.2019.01.071
83
XuY.SeshadriB.SarkarB.WangH.RumpelC.SparksD.et al (2018). Biochar Modulates Heavy Metal Toxicity and Improves Microbial Carbon Use Efficiency in Soil. Sci. Total Environ.621, 148–159. 10.1016/j.scitotenv.2017.11.214
84
YanC. C.WangF.GengH. H.LiuH. J.PuS. Y.TianZ. J.et al (2020). Integrating High-Throughput Sequencing and Metagenome Analysis to Reveal the Characteristic and Resistance Mechanism of Microbial Community in Metal Contaminated Sediments. Sci. Total. Environ.707, 11. 10.1016/j.scitotenv.2019.136116
85
YangY.CampbellC. D.ClarkL.CameronC. M.PatersonE. (2006). Microbial Indicators of Heavy Metal Contamination in Urban and Rural Soils. Chemosphere63, 1942–1952. 10.1016/j.chemosphere.2005.10.009
86
YuZ. J.LiuX. D.ZengX. B.YinH. Q.ZengW. M. (2020). Effect of Arsenic Pollution Extent on Microbial Community in Shimen Long-Term Arsenic-Contaminated Soil. Water Air Soil Poll.231 (7), 340. 10.1007/s11270-020-04716-6
87
ZhangJ.LiH.ZhouY.DouL.CaiL.MoL.et al (2018). Bioavailability and Soil-To-Crop Transfer of Heavy Metals in Farmland Soils: A Case Study in the Pearl River Delta, South China. Environ. Pollut.235, 710–719. 10.1016/j.envpol.2017.12.106
88
ZhangX. Y.ZhongT. Y.LiuL.OuyangX. Y. (2015). Impact of Soil Heavy Metal Pollution on Food Safety in China. Plos. One.10, 14. 10.1371/journal.pone.0135182
89
ZhaoF.-J.HarrisE.YanJ.MaJ.WuL.LiuW.et al (2013). Arsenic Methylation in Soils and its Relationship with MicrobialarsMAbundance and Diversity, and as Speciation in Rice. Environ. Sci. Technol.47, 7147–7154. 10.1021/es304977m
90
ZhenZ.WangS. B.LuoS. W.RenL.LiangY. Q.YangR. C.et al (2019). Significant Impacts of Both Total Amount and Availability of Heavy Metals on the Functions and Assembly of Soil Microbial Communities in Different Land Use Patterns. Front. Microbiol.10, 14. 10.3389/fmicb.2019.02293
91
ZhengS. N.WangQ.YuH. Y.HuangX. Z.LiF. B. (2020). Interactive Effects of Multiple Heavy Metal(loid)s on Their Bioavailability in Cocontaminated Paddy Soils in a Large Region. Sci. Total. Environ.708, 8. 10.1016/j.scitotenv.2019.135126
92
ZhongX.ChenZ. W.LiY. Y.DingK. B.LiuW. S.LiuY.et al (2020). Factors Influencing Heavy Metal Availability and Risk Assessment of Soils at Typical Metal Mines in Eastern China. J. Hazard. Mater.400. 123289. 10.1016/j.jhazmat.2020.123289
93
ZhuN.QiaoJ.YanT. (2019). Arsenic Immobilization through Regulated Ferrolysis in Paddy Field Amendment with Bismuth Impregnated Biochar. Sci. Total Environ.648, 993–1001. 10.1016/j.scitotenv.2018.08.200
94
ZouL.ZhangS.DuanD.LiangX.ShiJ.XuJ.et al (2018). Effects of Ferrous Sulfate Amendment and Water Management on Rice Growth and Metal(loid) Accumulation in Arsenic and Lead Co-contaminated Soil. Environ. Sci. Pollut. Res.25, 8888–8902. 10.1007/s11356-017-1175-8
Summary
Keywords
Pb-Zn mining site, arsenic, lead, microbial community, paddy soil
Citation
Zou L, Lu Y, Dai Y, Khan MI, Gustave W, Nie J, Liao Y, Tang X, Shi J and Xu J (2021) Spatial Variation in Microbial Community in Response to As and Pb Contamination in Paddy Soils Near a Pb-Zn Mining Site. Front. Environ. Sci. 9:630668. doi: 10.3389/fenvs.2021.630668
Received
18 November 2020
Accepted
15 April 2021
Published
28 April 2021
Volume
9 - 2021
Edited by
Jian Chen, Florida International University, United States
Reviewed by
Xiaomin Li, South China Normal University, China
Baogang Zhang, China University of Geosciences, China
Updates
Copyright
© 2021 Zou, Lu, Dai, Khan, Gustave, Nie, Liao, Tang, Shi and Xu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jun Nie, niejun197@163.com; Xianjin Tang, xianjin@zju.edu.cn
This article was submitted to Toxicology, Pollution and the Environment, a section of the journal Frontiers in Environmental Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.