- Open Access
Consistent responses of soil microbial taxonomic and functional attributes to mercury pollution across China
© The Author(s). 2018
- Received: 5 July 2018
- Accepted: 8 October 2018
- Published: 18 October 2018
The ecological consequences of mercury (Hg) pollution—one of the major pollutants worldwide—on microbial taxonomic and functional attributes remain poorly understood and largely unexplored. Using soils from two typical Hg-impacted regions across China, here, we evaluated the role of Hg pollution in regulating bacterial abundance, diversity, and co-occurrence network. We also investigated the associations between Hg contents and the relative abundance of microbial functional genes by analyzing the soil metagenomes from a subset of those sites.
We found that soil Hg largely influenced the taxonomic and functional attributes of microbial communities in the two studied regions. In general, Hg pollution was negatively related to bacterial abundance, but positively related to the diversity of bacteria in two separate regions. We also found some consistent associations between soil Hg contents and the community composition of bacteria. For example, soil total Hg content was positively related to the relative abundance of Firmicutes and Bacteroidetes in both paddy and upland soils. In contrast, the methylmercury (MeHg) concentration was negatively correlated to the relative abundance of Nitrospirae in the two types of soils. Increases in soil Hg pollution correlated with drastic changes in the relative abundance of ecological clusters within the co-occurrence network of bacterial communities for the two regions. Using metagenomic data, we were also able to detect the effect of Hg pollution on multiple functional genes relevant to key soil processes such as element cycles and Hg transformations (e.g., methylation and reduction).
Together, our study provides solid evidence that Hg pollution has predictable and significant effects on multiple taxonomic and functional attributes including bacterial abundance, diversity, and the relative abundance of ecological clusters and functional genes. Our results suggest an increase in soil Hg pollution linked to human activities will lead to predictable shifts in the taxonomic and functional attributes in the Hg-impacted areas, with potential implications for sustainable management of agricultural ecosystems and elsewhere.
- Soil microbial community
- Co-occurrence network
- Functional gene
- Mercury pollution
Environmental pollution resulting from human activities has a great impact on the biodiversity and functioning of terrestrial and aquatic ecosystems and is a major threat for human health across the globe [1–4]. Importantly, recent studies suggested that changing climate exacerbates global mercury (Hg) pollution by releasing historically stocked Hg in the permafrost [5, 6]. Such a threat is of global concern as Hg can be transported globally . The elevated Hg inputs into environments could have significant effects on soil biodiversity and their associated ecosystem functions [8, 9]. Studies conducted at the local scale seemed to support the idea that Hg pollution can alter soil microbial communities [10–12]. Moreover, previous work based on short-term incubation experiments demonstrated that Hg amendments can alter the soil microbiome and its ecological functions [13, 14]. However, much less is known on the role of Hg pollution in shaping the taxonomic and functional attributes of microbial communities in natural ecosystems across large spatial scale and different types of croplands.
Recent studies suggested that microbial taxa strongly co-occur within ecological network that often are called ecological clusters or modules [15, 16]. Such ecological clusters are expected to have multiple implications for the maintenance of soil fertility, decomposition, and ecological services in terrestrial environments [17–19]. The reason is that different microbial ecological clusters often follow very specific environmental preferences. For example, taxa within some ecological clusters have been found to strongly correlate with low or high pH, mesic or arid ecosystem, and nutrient availability and processes rates . However, despite the importance of these ecological clusters for the understanding of the soil microbiome, their response to long-term field Hg pollution remains largely unexplored . Increasing Hg concentrations could potentially impact the relative abundance of major ecological clusters with implications for ecosystem functioning. However, empirical evidence for such assumptions is currently lacking.
Moreover, the response of the functional attributes of the soil microbiome to Hg pollution has been rarely addressed. Such a task is challenging as soil microbial communities are highly diverse, and most soil microbial taxa remains uncharacterized [16, 21]. Recent advances in metagenomic sequencing can infer the functional traits of microbial communities [22, 23] and then determine which are sensitive to increased soil Hg pollution. We expect functional traits of soil microbial communities to follow predictable responses to Hg pollution. Such an expectation is based on the large body of literature reporting consistent microbial responses to other global changes such as climate change and nutrient additions [3, 23, 24].
Herein, we aimed to build a predictive understanding of the response of multiple soil taxonomic and functional attributes to increased Hg contents at the regional scale. Because Hg is a well-known pollutant with important implications for life on Earth , we hypothesized that soil Hg pollution will have consistent effects on taxonomic and functional attributes (e.g., relative abundance of functional genes) across different land use types. To do so, we collected three replicated soils at each of 47 sites (24 paddy fields and 23 adjacent upland fields), resulting in a total of 141 samples from Hg-impacted agricultural locations across two provinces in China. Most of the sites are surrounded by historical Hg mining areas with varying soil Hg contents, and their soils contained a wide range of Hg concentrations under similar vegetation types. These sites also included controls that are far away from the mining sites, with Hg contents similar to local background levels. As such, these locations provided a unique opportunity to empirically evaluate the response of microbial taxonomic and functional attributes to a gradient of soil Hg pollution. We characterized soil bacterial community composition via MiSeq Illumina platform. In addition, we investigated potential shifts in the relative abundance of functional genes linked to microbial communities by analyzing the soil metagenomes from a subset of those sites. Finally, we identified the associations between soil Hg concentrations and bacterial diversity, abundance, and ecological clusters (modules) using a combination of random forest analysis and structure equation modeling, after accounting for other critical environmental predictors.
Mercury pollution altered soil bacterial abundance and diversity
Mercury pollution altered the relative abundance of bacterial phyla
Mercury pollution shifted the relative abundance of ecological clusters within the correlation network
Mercury pollution is a significant predictor of soil bacterial community and ecological clusters after controlling for multiple environmental predictors
Random forest analysis suggested that, in general, Hg variables are major predictors of bacterial diversity and abundance (Additional file 1: Figure S5a, b) and the relative abundance of ecological clusters (Additional file 1: Figure S6a, b). As expected, our results indicate that other environmental factors were also important predictors of microbial communities, although the relative importance of these predictors was highly taxa and module dependent. See Additional file 1: Table S2 for correlations between microbial attributes and environmental predictors.
Functional gene relates to soil Hg pollution
Effects of Hg pollution on soil bacterial abundance and diversity in agricultural ecosystems
Both regression analysis and SEM demonstrate that long-term field Hg pollution had a negative effect on the bacterial abundance from the two land use types. Such a result is in contrast with a previous short-term incubation experiment where no effect of Hg pollution on bacterial abundance was found . Both Hg and MeHg are toxic to life, and thus, long-term field Hg pollution may inhibit microbial growth and then cause a decrease in soil microbial abundance across centuries . Importantly, the lack of an effect of Hg pollution on bacterial abundance in short-term experiments could be related to the fact that relic DNA from dead bacteria is also extracted and then quantified by qPCR [26, 27], obscuring the adverse effect of Hg on soil microbial abundance. Such an artifact should be largely diluted in long-term experiment and observational studies. Similarly, MeHg as the most toxic Hg speciation also negatively correlated to the bacterial abundance in the paddy soils, further emphasizing the strong effect of Hg pollution on bacterial abundance in soil. On the contrary, bacterial diversity tended to increase under moderate Hg content (< 53 mg kg−1, Fig. 1c), partially supporting the intermediate disturbance or stress hypothesis [12, 28]. However, our observations are different from those from short-term incubation experiments, in which Hg amendment has been reported to reduce soil bacterial diversity . As such, our approach suggests that large-scale observational studies are needed to understand the response of bacterial diversity to Hg pollution in real-world ecosystems. The soils examined in this study belong to Hg-mining areas which suffered serious Hg pollution for more than 600 years . Consequently, the microbial community from these soils might have had time to develop resistance to Hg stress . Alternatively, increasing Hg pollution might increase bacterial diversity by releasing subordinate microbial taxa via competition. These results are partially supported by our SEM suggesting a positive effect of MeHg on the diversity in upland soils, while the effect became complicated in paddy soils where no direct effect of total Hg or MeHg on the diversity was observed. However, Hg pollution may indirectly influence the bacterial diversity through impacts on soil properties, which were important drivers of soil microbial diversity [17, 21]. All of these observations illustrate the importance of Hg pollution in regulating the abundance and diversity of soil bacteria. Importantly, these results were maintained after accounting for multiple other drivers including soil and spatial properties.
Mercury pollution altered soil bacterial community composition and the distribution of ecological clusters
Importantly, we also found that Hg pollution consistently related to the relative abundance of the dominant bacterial phyla. For example, the relative abundance of Firmicutes and Bacteroidetes increased in response to elevated total Hg, but the relative abundance of Nitrospirae was negatively related to MeHg contents. Firmicutes and Bacteroidetes are known to be fast-growing opportunistic types of organisms that might benefit from environmental disturbance by taking over niches commonly occupied by other bacterial taxa. We also detected the genera of bacteria within these groups influenced by Hg pollutions (Additional file 1: Figure S3; Figure S4). For example, Sphingobacteriaceae Pedobacter and Clostridium sensu stricto have been reported to be resistant to Hg or heavy metals [30, 31]; however, much less is known about the mechanisms linking the responses of Chitinophagaceae Ferruginibacter to Hg pollution. In addition, our large-scale data suggest a decrease in the relative abundance of Nitrospirae with elevated soil MeHg. Previous studies have found that this phylum is highly sensitive to heavy metal stress at the local scale [12, 32]. The reduction in the relative abundance of Nitrospirae was mainly attributed to a relative decrease in the genus Nitrospirales 0319.6A21sp. in response to elevated MeHg. Increased Hg pollution is known to inhibit the processes of nitrogen cycle driven by soil functional assembles [9, 33], which could be an important result of the inhibiting effect of Hg to nitrifying bacteria such as Nitrospirae. These results suggest that the sensitive taxa could be used as potential ecological indicators for Hg pollution in terrestrial ecosystems.
Our random forest analysis suggested that Hg, in general, is a significant predictor of the ecologically preferential modules within the bacterial co-occurrence network after accounting for other key environmental predictors (Additional file 1: Figure S6). Thus, increases in soil Hg content also led to drastic changes in the co-occurrence network of soil bacterial communities in the studied regions. Interestingly, we found significant patterns of the relative abundances in some modules along increased Hg and MeHg pollution (Fig. 3c; Additional file 1: Table S6). For example, the relative abundance of module #4 tended to increase with elevated total Hg and MeHg in both paddy and upland soils. This module mainly consisted of potentially Hg-resistant microbes such as Proteobacteria, Bacteroidetes, and Actinobacteria, which were reported to contain the Hg-resistance gene merA [34, 35]. We also provide a list of winner (positive effects) and loser (negative effects) community assemblies in response to Hg pollution (Additional file 2: Data S1), which can be used to test for similarities in the response of soil microbial communities to Hg pollution worldwide.
Although SEM is quite a conservative procedure, the results support similar associations between Hg and the relative abundance of modules across the two types of land use. For example, soil Hg had direct effects on the relative abundance of module #1 in paddy soils (Fig. 4a), and a similar effect on this module in upland soils could be reflected through indirect impacts of Hg on other modules (i.e., Mod#2, Mod#3, and Mod#4) due to their interactions (Fig. 4b: BOX2). It is not surprising to note the interactions between the modules in the co-occurrence networks, because increases in a given ecological cluster were often followed by declines in the relative abundance of other ecological clusters . These findings suggest that variability of Hg contents in the ecosystem might change the distributions of ecologically preferential clusters. Of course, as expected, other environmental predictors strongly associate with the relative abundance of microbial ecological clusters (Fig. 4, Additional file 1: Table S2; Additional file 1: Figure S6), which had been well recognized in many previous studies [15–17].
Predicting soil Hg pollution from the variability of microbial functional genes
We found that Hg pollution was also associated with important functional traits within microbial communities. In fact, here, we provide a list of functional genes which were strongly correlated to Hg contents across soils. Of course, we acknowledge the current limitations for linking genes to Hg pollution across large spatial scales due to the correlative nature of these analyses (Fig. 5; Additional file 1: Figure S7). Such results are in agreement with previous studies emphasizing significant effects of Hg stress on microbially driven processes (e.g., nitrification potential and dehydrogenase enzyme activity) driven by diverse soil microorganisms [9, 37]. For example, the relative abundance of genes involved in member protein tended to increase with the level of soil Hg pollution, suggesting that moderated Hg stress may stimulate the enzyme activities responsible for Hg transportation through the cell membrane . We also observed a significant increase in the relative abundances of genes associated with dehydrogenase with elevated Hg. These findings are consistent with those of previous studies showing upregulation of dehydrogenase activities exposed to metals . In addition, we identified some genes that are relevant to Hg transport (i.e., iron transport system) and methylation processes (i.e., CoA synthetase) . Previous studies has suggested that inorganic Hg(II) could be transported into microbial cells probably through ion transport system [41, 42], and the cellular Hg(II) is subsequent methylated to highly neurotoxic MeHg by methylating genes hgcAB . Thus, our predicted iron transport system might be relevant to soil MeHg formation. Furthermore, Hg methylation is an enzyme-catalyzed process associated with the reductive CoA pathway , which may explain our results that the CoA synthetase is a significant predictor for MeHg in the soil. Overall, these gene predictors are associated with nutrient metabolisms and Hg transformations, which are also important biomarkers of soil Hg pollution.
Together, our study represents one of the first attempts to empirically assess how the soil microbiomes respond to long-term Hg pollution across large spatial scales and land use type. We provides solid evidence that Hg pollution has predictable and significant effects on multiple taxonomic and functional microbial attributes including bacterial abundance, diversity, and the relative abundance of ecological clusters and functional genes. Such results are maintained after accounting for other important environmental predictors of soil microbial communities. In general, Hg pollution was negatively related to the bacterial abundance. The relative abundance of Bateroidetes and Firmicutes increased with elevated Hg pollution, while the relative abundance of Nitrospirae was negatively correlated to MeHg. Mercury-induced shifts in the ecological clusters of co-occurrence network and functional traits of microbial communities could have important implications for soil biogeochemical cycling and service functioning of the ecosystem. This work moves us towards a more predictable understanding of the response of microbial communities and their potential function to Hg pollution—a worldwide threat derived from both global warming and intensive anthropogenic activities.
Study area and sampling
The soil samples were collected around two typical Hg mining areas (Fenghuang, FH and Wanshan, WS) in southwest China. Fenghuang is located in West Hunan Province, and Wanshan in East Guizhou Province. Both of them are the major grain-producing areas in China, with a long history Hg mining for more than 600 years . Historical discharges from Hg mining operations and ongoing atmospheric deposition has led to heavy Hg pollution in these areas . We collected soil samples from 24 locations (Additional file 1: Figure S8) along downstream the mining sites to obtain samples with a large gradient of soil Hg concentrations (ranges from 0.27 to 52.4 mg kg−1), and most of them were above risk control value for paddy soils in China (0.6 mg kg−1) . Two typical types of land use including paddy fields (24 sub-sites) and adjacent uplands (23 sub-sites, maize planting) were chosen to compare the effects of land use on the microbial responses to long-term soil Hg pollution. We tried to choose locations where both paddy and upland fields are the main crop planting, and sporadic farmlands are excluded. Three replicated soil (0–15 cm depth) samples were collected at each of the 47 sites to account for small scale variation in Hg pollution. Consequently, a total of 141 soil samples were obtained from the 24 paddy (72 samples) and upland (69 samples) fields. These locations also included two control sites (paddy and adjacent upland fields) located in the natural reserve with little pollution (i.e., location no. 1 in Additional file 1: Figure S8). Total Hg concentrations in the control sites are similar the levels of local background values (0.5~0.7 mg kg−1), which are similar to the Hg risk control value for paddy soil in China . We omitted an upland site due to the lack of representative field in the location. Sampling was conducted in August 2016 before harvesting. Collected soil samples were homogenized and sieved (2.0 mm) to get a representative microbial community [20, 47]. The sieved samples were subsequently divided into two sub-samples. One sub-sample was stored at − 20 °C for microbial analysis, while the other sub-sample was air-dried for the analysis of heavy metals and soil properties.
Measurement of soil heavy metals and chemical indexes
For Hg analysis, 0.20 g of each soil was digested with 10 mL mixed solution (2 mol L−1 HNO3 and 4 mol L−1 HCl) in a Teflon tube at 95 °C for 2 h. The total amount of Hg in these extracts was determined via cold vapor atomic fluorescence spectrometry (CVAFS) , and the method detection limit was 0.2 ng L−1. For the analysis of other heavy metals (i.e., Cu, Pb, Cd, Zn, Ni, and As), 0.30 g soil were digested by trace mixed acids (9.0 mL of HNO3 and 3.0 mL of HF) in a MARS microwave digestion system (CEM, USA). The concentrations of the heavy metals in the final solution were measured using a 7700X Inductively Coupled Plasma-Mass Spectrometer (Agilent, USA). For MeHg, 0.40 g soil was used for MeHg extraction using CuSO4-methanol. The amount MeHg was determined using an automated MeHg analyzer (TEKRAN 2700 GC-CVAFS) . Soil pH was determined on a fresh soil to water ratio of 1: 2.5 using a Delta pH-meter. Soil organic carbon (SOC) was measured using the K2CrO7 oxidation titration method . Total carbon (TC) and total nitrogen (TN) in soils were determined on a LECO TureMac Macro CN analyzer (LECO, St. Joseph, MI, USA). Dissolved organic carbon (DOC) in the soil was extracted with 0.5 M K2SO4 at a ratio of 1:5 by shaking at 200 rpm for 1 h and filtering through a 0.45-μm Millipore filter , and the DOC concentration was measured by a TOC analyzer (TOC-L Analyzer, Shimadzu, Japan). The ammonium and nitrate concentrations in the filtered extracts were analyzed within 24 h using a continuous flow analyzer (SAN++, Skalar, Breda, Holland).
Analysis of soil bacterial communities through Illumina MiSeq sequencing
The genomic DNA was isolated from 0.30 g of soil using the MoBio PowerSoil DNA Isolation Kit (QIAGEN Inc., USA) following the manufacturer’s instructions. The abundance of bacteria was estimated by quantifying the 16S rRNA gene copy number on an iCycler iQ5 thermocycler (Bio-Rad, USA) using the primer pairs Eub338F/Eub518R . To evaluate the bacterial community composition, the V4 region of the bacterial 16S rRNA gene was amplified using the primer pairs of 338F/806R . The 50 μl PCR reaction mixtures consisted of 25 μl PremixTaq™ (Takara Biotechnology, Dalian, China), 1 μl of each primer (10 μM), 3 μl of template DNA, and 20 μl of sterilized ddH2O. The resultant PCR products were purified using the Wizard® SV Gel and PCR Clean-Up System (Promega, San Luis Obispo, USA). The purified amplicons were equimolarly mixed, and 2 × 250 bp paired-end sequencing was carried out on an Illumina Miseq sequencer (Illumina Inc., San Diego, USA). Raw reads generated from the MiSeq paired-end sequencing were merged together using the Fast Length Adjustment of Short reads (FLASH). A chimera filtering approach UPARSE was employed as the Operational Taxonomic Unit (OTU or phylotype) picking strategy at 97% sequence similarity. We randomly selected an even number (30,212) of reads from each sample to account for variability in sequencing depth before downstream analysis (Additional file 1: Figure S9). Representative sequences from individual OTUs generated in UPARSE were processed using the Quantitative Insights into Microbial Ecology (QIIME) pipeline. The bacterial diversity index was calculated based on 97% OTU similarity of obtained bacterial sequences. The taxonomic identity of all phylotypes was determined using The SILVA ribosomal RNA gene database project .
Shotgun metagenomic sequencing and gene analysis
Seven paddy soils and three upland soils were selected from our 141 samples for metagenomic analysis. These samples were chosen, after Hg analysis, to produce a wide range of Hg concentrations (ranged from 0.84 to 32.43 mg kg−1). Sequencing was performed using an Illumina PE150 (Illumina Inc.) at Majorbio, Inc., Shanghai, China. Raw reads (150 bp in length) were trimmed to remove low-quality reads. Paired reads of shotgun metagenomic sequences were merged with FLASH using default parameters . Using MBLASTX, merged reads were also mapped against the protein sequence of the KEGG database (E value cutoff 1e−6), and the relative abundance of each KO gene was also calculated. Additional details on methodology are provided in the Additional file 3. To estimate the influence of elevated Hg and MeHg contents on these genes, we focused on the KO genes related to microbial metabolism, Hg transformations, and other sensitive ones.
Correlation network analyses
We established a co-occurrence network to identify ecological clusters of bacterial taxa across the collected 141 soil samples. A single correlation network including all samples was conducted so that the identified ecological clusters are directly comparable across land use types. We kept those taxa accounting for more than 80% of the relative abundance of soil bacteria (1073 bacterial taxa in Additional file 2: Data S1). We then calculated all pairwise Spearman’s rank correlations (ρ) between all soil bacterial taxa using the psych package of the R statistical software (http://cran.r-project.org/) and focused exclusively on positive correlations as they provide information on bacterial taxa that may respond similarly to environmental conditions . We considered a co-occurrence to be robust if the Spearman’s correlation coefficient was > 0.25 and P < 0.01 . The network was visualized with the interactive platform Gephi . Finally, we used default parameters from the interactive platform Gephi to identify ecological clusters (i.e., modules) of soil taxa strongly interacting with each other . The relative abundance of each ecological cluster was computed by averaging the standardized relative abundances (z-score) of the taxa that belong to each ecological cluster .
We first identified the correlations between total Hg and MeHg with soil bacterial abundance, diversity (Shannon), and the relative abundance of ecological clusters in both paddy and upland soil, using linear or cubic models. We used both regression analyses and Pearson correlation to evaluate correlations between total Hg and MeHg and the relative abundance of ecological clusters, and bacterial community composition for two land use types using one ANOVA, with land use type as the fixed factor. We conducted a classification random forest analysis [56, 57] to identify the statistically significant predictors of the bacterial diversity, total abundance, and the relative abundance of ecological clusters. The major aim of these analyses were to test whether total Hg or MeHg pollution are significant predictors of microbial attributes, after accounting for other key environmental predictors including location, soil properties, heavy metals, and nutrient availability. We also used random forest analysis to identify the genus predictors of the phyla/classes that were significantly correlated to soil Hg. The random forest model determined the importance of each predictor variable via evaluating the decrease in prediction accuracy (i.e., increase in the mean square error between observations and OOB predictions) when the data for that predictor are randomly permuted, as previously described . These analyses were conducted using the rfPermute package  of the R statistical software (http://cran.r-project.org/). Additionally, we used Spearman’s correlation analyses to further evaluate the correlations between bacterial diversity (Shannon), total abundance, and the relative abundance of ecological clusters with environmental predictors. We also used random forest analysis to identify the major functional genes predicting the concentrations of Hg and MeHg in our soils. In these analysis, functional genes are used as predictors of Hg concentrations. After this, we used linear regressions to evaluate the direction of the relationships between the relative abundance of selected genes and soil total Hg concentrations.
We then used structural equation modeling (SEM)  to further clarify associations of Hg (total Hg and MeHg) contents with the bacterial abundance, diversity, and the relative abundance of ecological clusters (modules). Unlike the analysis of regression or ANOVA, SEM allows the ability to separate multiple effect pathways and consider them as parts of a system and thus is useful for isolating the complex relationships among environmental factors commonly found in natural ecosystems [60, 61]. The probability that a path coefficient differs from zero was tested using bootstrap resampling. Bootstrapping is preferred to the classical maximum-likelihood estimation in these cases because in bootstrapping, probability assessments are not based on the assumption that the data match a particular theoretical distribution. Our model also includes spatial autocorrelation (latitude and longitude) and soil properties, which were represented by the three major components by performing principal component analysis of soil variables including pH, soil organic carbon (SOC), C: N, and others (Additional file 1: Table S3). We first established an a priori model according to our current knowledge of environmental variable impacts on soil microbiomes (Additional file 1: Figure S10). The data matrix was fitted to the model using the maximum-likelihood estimation method. There is no single universally accepted test of overall goodness of fit for SEM. Thus, we used the chi-square test (χ2; the model has a good fit when 0 ≤ χ2/d.o.f ≤ 2 and 0.05 < P ≤ 1.00) and the root mean square error of approximation (RMSEA; the model has a good fit when RMSEA 0 ≤ RMSEA ≤ 0.05 and 0.10 < P ≤ 1.00 . The SEM analyses were performed using AMOS 21.0 (SPSS Inc., Chicago, IL, USA).
We thank Xinquan Zhou and Wenli Tang for their assistance in mercury analysis, and Pengpeng Xu for technical support in bioinformatics analysis. We gratefully acknowledge Chao Xiong and Anhui Ge for their assistance during field sampling. We thank Dr. PM Chalk (The University of Melbourne) for his language editing of this manuscript, and the anonymous reviewers for their constructive comments on this paper.
This research was supported in part by the National Natural Science Foundation of China (41877120, 41571453). M.D-B. acknowledges the support from the Marie Sklodowska-Curie Actions of the Horizon 2020 Framework Programme H2020-MSCA-IF-2016 under REA grant agreement no. 702057.
Availability of data and materials
The soil metagenomic dataset has been deposited to the NCBI Sequence Read Archive under accession number PRJNA450451.
YRL and MDB planned the research; YRL and LB conducted the laboratory analyses and collected the raw data. MDB and YRL analyzed the output data and performed all the statistical analyses; YRL, MDB, JZ, and JZH wrote the first draft of the manuscript, and all authors contributed to revisions. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Bissett A, Brown MV, Siciliano SD, Thrall PH. Microbial community responses to anthropogenically induced environmental change: towards a systems approach. Ecol Lett. 2013;16:128–39.View ArticleGoogle Scholar
- Griffiths BS, Philippot L. Insights into the resistance and resilience of the soil microbial community. FEMS Microbiol Rev. 2013;37(2):112–29.View ArticleGoogle Scholar
- Jiao S, Zhang Z, Yang F, Lin Y, Chen W, Wei G. Temporal dynamics of microbial communities in microcosms in response to pollutants. Mol Ecol. 2017;26(3):923–36.View ArticleGoogle Scholar
- Alava JJ, Cheung WWL, Ross PS, Sumaila UR. Climate change–contaminant interactions in marine food webs: toward a conceptual framework. Glob Chang Biol. 2017;23(10):3984–4001.View ArticleGoogle Scholar
- Obrist D, Agnan Y, Jiskra M, Olson CL, Colegrove DP, Hueber J, Moore CW, Sonke JE, Helmig D. Tundra uptake of atmospheric elemental mercury drives Arctic mercury pollution. Nature. 2017;547(7662):201–4.View ArticleGoogle Scholar
- Schuster PF, Schaefer KM, Aiken GR, Antweiler RC, Dewild JF, Gryziec JD, Gusmeroli A, Hugelius G, Jafarov E, Krabbenhoft DP, et al. Permafrost stores a globally significant amount of mercury. Geophysl Res Lett. 2018;45(3):1463–71.View ArticleGoogle Scholar
- Krabbenhoft DP, Sunderland EM. Global change and mercury. Science. 2013;341(6153):1457–8.View ArticleGoogle Scholar
- Mahbub KR, Krishnan K, Naidu R, Andrews S, Megharaj M. Mercury toxicity to terrestrial biota. Ecol Indic. 2017;74:451–62.View ArticleGoogle Scholar
- Liu Y-R, Zheng Y-M, Shen J-P, Zhang L-M, He J-Z. Effects of mercury on the activity and community composition of soil ammonia oxidizers. Environ Sci Pollut Res. 2010;17(6):1237–44.View ArticleGoogle Scholar
- Müller AK, Westergaard K, Christensen S, Sørensen SJ. The effect of long-term mercury pollution on the soil microbial community. FEMS Microbiol Ecol. 2001;36(1):11–9.View ArticleGoogle Scholar
- Harris-Hellal J, Vallaeys T, Garnier-Zarli E, Bousserrhine N. Effects of mercury on soil microbial communities in tropical soils of French Guyana. Appl Soil Ecol. 2009;41(1):59–68.View ArticleGoogle Scholar
- Liu Y-R, Wang J-J, Zheng Y-M, Zhang L-M, He J-Z. Patterns of bacterial diversity along a long-term mercury-contaminated gradient in the paddy soils. Microb Ecol. 2014;68(3):575–83.View ArticleGoogle Scholar
- Frossard A, Hartmann M, Frey B. Tolerance of the forest soil microbiome to increasing mercury concentrations. Soil Biol Biochem. 2017;105:162–76.View ArticleGoogle Scholar
- Frey B, Rieder SR. Response of forest soil bacterial communities to mercury chloride application. Soil Biol Biochem. 2013;65:329–37.View ArticleGoogle Scholar
- Menezes AB, Prendergast-Miller MT, Richardson AE, Toscas P, Farrell M, Macdonald LM, Baker G, Wark T, Thrall PH. Network analysis reveals that bacteria and fungi form modules that correlate independently with soil parameters. Environ Microbiol. 2015;17(8):2677–89.View ArticleGoogle Scholar
- Delgado-Baquerizo M, Oliverio AM, Brewer TE, Benavent-González A, Eldridge DJ, Bardgett RD, Maestre FT, Singh BK, Fierer N. A global atlas of the dominant bacteria found in soil. Science. 2018;359(6373):320–5.View ArticleGoogle Scholar
- Delgado-Baquerizo M, Reith F, Dennis PG, Hamonts K, Powell JR, Young A, Singh BK, Bissett A. Ecological drivers of soil microbial diversity and soil biological networks in the Southern Hemisphere. Ecology. 2018;99(3):583–96.View ArticleGoogle Scholar
- Van Der Heijden MGA, Bardgett RD, Van Straalen NM. The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol Lett. 2008;11(3):296–310.View ArticleGoogle Scholar
- Shi Y, Li Y, Xiang X, Sun R, Yang T, He D, Zhang K, Ni Y, Zhu Y-G, Adams JM, et al. Spatial scale affects the relative role of stochasticity versus determinism in soil bacterial communities in wheat fields across the North China Plain. Microbiome. 2018;6(1):27.View ArticleGoogle Scholar
- Frossard A, Donhauser J, Mestrot A, Gygax S, Bååth E, Frey B. Long- and short-term effects of mercury pollution on the soil microbiome. Soil Biol Biochem. 2018;120:191–9.View ArticleGoogle Scholar
- Fierer N. Embracing the unknown: disentangling the complexities of the soil microbiome. Nat Rev Micro. 2017;15(10):579.View ArticleGoogle Scholar
- Zhang X, Johnston ER, Barberán A, Ren Y, Lü X, Han X. Decreased plant productivity resulting from plant group removal experiment constrains soil microbial functional diversity. Glob Chang Biol. 2017;23(10):4318–32.View ArticleGoogle Scholar
- Leff JW, Jones SE, Prober SM, Barberán A, Borer ET, Firn JL, Harpole WS, Hobbie SE, Hofmockel KS, Knops JMH, et al. Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe. Proc Natl Acad Sci U S A. 2015;112(35):10967–72.View ArticleGoogle Scholar
- Zhang X, Johnston ER, Li L, Konstantinidis KT, Han X. Experimental warming reveals positive feedbacks to climate change in the Eurasian Steppe. ISME J. 2017;11(4):885–95.View ArticleGoogle Scholar
- Obrist D, Kirk JL, Zhang L, Sunderland EM, Jiskra M, Selin NE. A review of global environmental mercury processes in response to human and natural perturbations: changes of emissions, climate, and land use. Ambio. 2018;47(2):116–40.View ArticleGoogle Scholar
- Frostegård Å, Courtois S, Ramisse V, Clerc S, Bernillon D, Le Gall F, Jeannin P, Nesme X, Simonet P. Quantification of bias related to the extraction of DNA directly from soils. Appl Environ Microbiol. 1999;65(12):5409–20.PubMedPubMed CentralGoogle Scholar
- Carini P, Marsden PJ, Leff JW, Morgan EE, Strickland MS, Fierer N. Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat Microbiol. 2017;2:16242.View ArticleGoogle Scholar
- Connell JH. Diversity in tropical rain forests and coral reefs. Science. 1978;199(4335):1302–10.View ArticleGoogle Scholar
- Feng X, Li P, Qiu G, Wang S, Li G, Shang L, Meng B, Jiang H, Bai W, Li Z, et al. Human exposure to methylmercury through rice intake in mercury mining areas, Guizhou Province, China. Environ Sci Technol. 2008;42(1):326–32.View ArticleGoogle Scholar
- Lee J, Jung Y-J, Lee HK, Hong SG, Kim O-S. Complete genome sequence of Pedobacter cryoconitis PAMC 27485, a CRISPR-Cas system-containing psychrophile isolated from Antarctica. J Biotechnol. 2016;226:74–5.View ArticleGoogle Scholar
- Wu S, Li Y, Wang P, Zhong L, Qiu L, Chen J. Shifts of microbial community structure in soils of a photovoltaic plant observed using tag-encoded pyrosequencing of 16S rRNA. Appl Microbiol Biotechnol. 2016;100:3735–45.View ArticleGoogle Scholar
- Das S, Jean J-S, Kar S, Liu C-C. Changes in bacterial community structure and abundance in agricultural soils under varying levels of arsenic contamination. Geomicrobiol J. 2013;30(7):635–44.View ArticleGoogle Scholar
- Zhou Z, Zheng Y, Shen J, Zhang L, Liu Y, He J. Responses of activities, abundances and community structures of soil denitrifiers to short-term mercury stress. J Environm Sci. 2012;24(3):369–75.View ArticleGoogle Scholar
- Barkay T, Miller SM, Summers AO. FEMS Microbiol Rev. 2003;27:355–84.View ArticleGoogle Scholar
- Møller AK, Barkay T, Hansen MA, Norman A, Hansen LH, Sørensen SJ, Boyd ES, Kroer N. Mercuric reductase genes (merA) and mercury resistance plasmids in high Arctic snow, freshwater and sea-ice brine. FEMS Microbiol Ecol. 2014;87(1):52–63.View ArticleGoogle Scholar
- Delgado-Baquerizo M, Eldridge David J, Travers Samantha K, Val J, Oliver I, Bissett A. Effects of climate legacies on above- and belowground community assembly. Glob Chang Biol. 2018;24:4330–9.View ArticleGoogle Scholar
- Mahbub KR, Krishnan K, Megharaj M, Naidu R. Mercury inhibits soil enzyme activity in a lower concentration than the guideline value. Bull Environ Contam Toxicol. 2016;96(1):76–82.View ArticleGoogle Scholar
- Barkay T, Wagner-Dobler I. Microbial transformations of mercury: potentials, challenges, and achievements in controlling mercury toxicity in the environment. Adv Appl Microbiol. 2005;57:1):1–52.PubMedGoogle Scholar
- Lim JC, Thevarajoo S, Selvaratnam C, Goh KM, Shamsir MS, Ibrahim Z, Chong CS. Global transcriptomic response of Anoxybacillus sp. SK 3-4 to aluminum exposure. J Basic Microbiol. 2017;57(2):151–61.View ArticleGoogle Scholar
- Hsu-Kim H, Kucharzyk KH, Zhang T, Deshusses MA. Mechanisms regulating mercury bioavailability for methylating microorganisms in the aquatic environment: a critical review. Environ Sci Technol. 2013;47(6):2441–56.View ArticleGoogle Scholar
- Szczuka A, Morel FMM, Schaefer JK. The effect of thiols, zinc, and redox conditions on Hg uptake in Shewanella oneidensis. Environ Sci Technol. 2015;49(12):7432–8.View ArticleGoogle Scholar
- Schaefer JK, Rocks SS, Zheng W, Liang L, Gu B, Morel FMM. Active transport, substrate specificity, and methylation of Hg(II) in anaerobic bacteria. Proc Natl Acad Sci U S A. 2011;108(21):8714–9.View ArticleGoogle Scholar
- Parks JM, Johs A, Podar M, Bridou R, Hurt RA, Smith SD, Tomanicek SJ, Qian Y, Brown SD, Brandt CC, et al. The genetic basis for bacterial mercury methylation. Science. 2013;339(6125):1332–5.View ArticleGoogle Scholar
- Choi S-C, Chase T, Bartha R. Metabolic pathways leading to mercury methylation in Desulfovibrio desulfuricans LS. Appl Environ Microbiol. 1994;60(11):4072–7.PubMedPubMed CentralGoogle Scholar
- Meng M, Li B, Shao J-J, Wang T, He B, Shi J-J, Ye Z-H, Jiang G-B. Accumulation of total mercury and methylmercury in rice plants collected from different mining areas in China. Environ Pollut. 2014;184:179–86.View ArticleGoogle Scholar
- Ministry of Ecology and Environment of the People's Republic of China. Soil environmental quality, risk control standard for soil contamination of agricultural land (GB 15618—2018), 2018. (in Chinese).Google Scholar
- Delgado-Baquerizo M, Maestre FT, Reich PB, Jeffries TC, Gaitan JJ, Encinar D, Berdugo M, Campbell CD, Singh BK. Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat Commun. 2016;7:10541.View ArticleGoogle Scholar
- Liu Y-R, Dong J-X, Zhang Q-G, Wang J-T, Han L-L, Zeng J, He J-Z. Longitudinal occurrence of methylmercury in terrestrial ecosystems of the Tibetan Plateau. Environ Pollut. 2016;218:1342–9.View ArticleGoogle Scholar
- Walkley A, Black IA. An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil Sci. 1934;37(1):29–38.View ArticleGoogle Scholar
- Jones DL, Willett VB. Experimental evaluation of methods to quantify dissolved organic nitrogen (DON) and dissolved organic carbon (DOC) in soil. Soil Biol Biochem. 2006;38(5):991–9.View ArticleGoogle Scholar
- Cregger MA, Schadt CW, McDowell NG, Pockman WT, Classen AT. Response of the soil microbial community to changes in precipitation in a semiarid ecosystem. Appl Environ Microbiol. 2012;78(24):8587–94.View ArticleGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(D1):D590–6.View ArticleGoogle Scholar
- Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.View ArticleGoogle Scholar
- Barberán A, Casamayor EO, Fierer N. The microbial contribution to macroecology. Front Microbiol. 2014;5:203.Google Scholar
- Bastian M, Heymann S, Jacomy M: Gephi: an open source software for exploring and manipulating networks; 2009.Google Scholar
- Breiman L. Random forests. Mach Learn. 2001;45(1):5–32.View ArticleGoogle Scholar
- Trivedi P, Delgado-Baquerizo M, Trivedi C, Hu H-W, Anderson I-C, Jeffries T-C, Zhou J, Singh B. Microbial regulation of the soil carbon cycle: evidence from gene-enzyme relationships. ISME J. 2016;10:2593–604.View ArticleGoogle Scholar
- Delgado-Baquerizo M, Gallardo A, Covelo F, Prado-Comesaña A, Ochoa V, Maestre FT. Differences in thallus chemistry are related to species-specific effects of biocrust-forming lichens on soil nutrients and microbial communities. Funct Ecol. 2015;29:1087–98.View ArticleGoogle Scholar
- Archer E. rfPermute, Estimate permutation p-values for random forest importance metrics. R package version 1.5.2. 2016.Google Scholar
- Grace JB. Structural equation modeling and natural systems. Cambridge: Cambridge University Press; 2016.Google Scholar
- Maestre FT, Quero JL, Gotelli NJ, Escudero A, Ochoa V, Delgado-Baquerizo M, García-Gómez M, Bowker MA, Soliveres S, Escolar C. Plant species richness and ecosystem multifunctionality in global drylands. Science. 2012;335:214–8.View ArticleGoogle Scholar
- Jing X, Sanders NJ, Shi Y, Chu H, Classen AT, Zhao K, Chen L, Shi Y, Jiang Y, He J-S. The links between ecosystem multifunctionality and above- and belowground biodiversity are mediated by climate. Nat Commun. 2015;6:8159.View ArticleGoogle Scholar