Abstract

Heavy metals have adverse effects on soil ecology. Given the toxicity of heavy metals, there is an urgent need to select an appropriate indicator that will aid in monitoring their biological effects on soil ecosystems. By combining different monitoring techniques for various aspects of microbiology, the effects of heavy metals on soil microorganisms near a smelter were studied. Our goal was to determine whether proteins could be a proper indicator for soil pollution. This study demonstrated that the activities of acid phosphatase and dehydrogenase, as well as the levels of microbial biomass carbon and proteins, were negatively affected by heavy metals. In addition, significantly negative correlations were observed between these microbial indicators and heavy metals. Denaturing gradient gel electrophoresis analysis was used in this study to demonstrate that heavy metals also have a significantly negative effect on soil microbial diversity and community structure. The soil protein expression was similar across different soils, but a large quantity of presumably low molecular weight protein was observed only in contaminated soil. Based on this research, we determined that the soil protein concentration was more sensitive to heavy metals than acid phosphatase, dehydrogenase, or microbial biomass carbon because it was more dramatically decreased in the contaminated soils. Therefore, we concluded that the soil protein level has great potential to be a sensitive indicator of soil contamination. Further research is essential, particularly to identify the low molecular weight protein that only appears in contaminated soil, so that further insight can be gained into the responses of microbes to heavy metals.

1. Introduction

Heavy metals (e.g., Cd, Cr, Cu, Pb, and Zn) can be introduced into soils from several sources, such as waste from mines and smelters, atmospheric deposition, animal manures, sewage sludge, and, in some circumstances, inorganic fertilizer. Once these elements enter the soil, they can remain for extremely long periods and are difficult to remove [1]. Heavy metals are known to be toxic to most organisms when present in excessive concentrations [2]. Given the toxicity of these elements, there is an urgent need to monitor their biological effects on soil ecosystems.

Soil microbes play vital roles in the recycling of plant nutrients, the maintenance of soil structure, and the detoxification of poisonous chemicals [3]. Therefore, the diversity of microbial communities and the activity of microorganisms are important indicators of soil quality. Alterations in the composition of microbial communities have often been considered effective indicators for soil contamination [4]. Many reports have shown that either long-term or short-term exposure of soil to heavy metals results in a reduction of microbial diversity and microorganism activities [3, 5].

Although soil enzyme activity and microbial biomass are the most tested parameters in soil quality monitoring, these parameters still have many limitations [2, 6]. They need to be supplemented by other valid microbial indicators, working together to fully assess soil contamination. The recent emergence of soil proteomics could provide a new perspective on soil microbial activities [7]. Proteins are the chief actors in organism metabolism processes, and they fall into a variety of categories; for example, enzymes are proteins with a catalytic function. Proteins work together to complete various biological functions. Soil proteins primarily come from microorganisms as well as from flora and fauna tissues. The structure of the microbial community, biomass, and microbial status often determine the synthesis levels of soil proteins. Therefore, any change to the microbial composition, which can be caused by environmental changes such as contamination, could reflect in the composition and expressions of soil proteins [8]. Soil proteomics studies aim to investigate the spatial and temporal changes of proteins extracted from soils because proteins are the functional components of the microbial genomic expression products. Such studies more conclusively determine the ecological functions of soil microbes and their roles in soil pollutant transportation and transformation [911]. Therefore, proteomics analyses have great potential in soil pollution assessment.

Currently, soil proteomics are predominantly used in studies of the biogeochemical cycle [1214] and the rhizosphere soil microecosystem [15, 16]. However, there is little literature on the monitoring of contaminated soil, and the existing studies were all carried out under controlled conditions [8, 17]. Unlike a controlled laboratory, natural soil is a complicated ecosystem, and various factors combine to influence the transportation and transformation of pollutants as well as the structure of the microorganism community. Consequently, we were interested in examining the modification of soil proteins under heavy metals stress and determining whether proteins can be a useful indicator of natural soil pollution.

In this study, a typical heavy metal smelter was selected. We intended to study the effects of heavy metals on soil microbes by combining different biochemical analysis techniques. Three compounds were used for heavy metal extraction, CaCl2, EDTA, and DTPA, to assess the bioavailability of trace metals in soil [18, 19]. Soil microbiological and biochemical properties (including enzyme activities, protein contents, and microbial biomass carbon levels) were measured to evaluate the effects of heavy metals. A soil microbial community analysis was carried out using molecular biology techniques. Partial 16S rDNA genes were amplified from soil microorganism community DNA by the polymerase chain reaction (PCR), using primers that bind to evolutionarily conserved regions within the bacteria and actinomycete genes. The diverse PCR amplified products were transformed to genetic fingerprints using denaturing gradient gel electrophoresis (DGGE) [3, 20]. The soil protein composition was analyzed by SDS-polyacrylamide gel electrophoresis (SDS-PAGE) [7].

2. Materials and Methods

2.1. Soil Research Region and Soil Sampling

The soil area selected for study is downwind of a copper-zinc smelter in Fuyang County, Zhejiang Province, China. The soil in the area has been contaminated by heavy metals, primarily Cu and Zn. This area was formerly used to cultivate crops. Soil samples were taken in April 2012 at nine points along a pollution gradient of Cu and Zn (Table 1). At each point, three replicated samples were taken, and each sample was composed of three soil cores 5 cm in diameter and 20 cm in depth (approximately 500 g). Samples were taken randomly from different areas at each site. Each replicate soil sample was homogenized thoroughly. The site 2,000 m away from the smelter was chosen as the reference site, and the soil collected from it was marked “CK.”

2.2. Sample Treatment and Analysis

Fresh soil samples were sieved ( mm) by a nylon sieve to remove stones, large pieces of plant material, and soil animals. Portions of the samples were kept moist in the dark at 4°C to determine the microbial biomass and then stored at −70°C to extract soil DNA. Other portions were freeze-dried and passed through a 0.25 mm nylon sieve prior to soil protein extraction. The remaining soil was air-dried under cool conditions (approximately 20°C), sieved through a 0.25 mm nylon sieve, and stored at 4°C to analyze the enzyme activity, soil pH, organic matter, total heavy metals, and bioavailable heavy metals.

The soil pH was determined with distilled water in a ratio (soil : water) of 1 : 2.5 (w/v) using a pH meter (Orion 5-Star, Thermo). The soil was heated at 105°C for 12 h and then combusted at 550°C overnight to measure the organic matter content [21]. The samples were digested using the method described by Tang et al. [22]. The dry sample (0.2 g) was weighed and digested with a mixture of nitric acid (HNO3), hydrofluoric acid (HF), and perchloric acid (HClO4). The total concentrations of Cd, Cr, Cu, Pb, and Zn were measured by flame atomic absorption spectrometry (MKIL-M6, Thermo).

The soil bioavailable heavy metal fractions were estimated by extraction with 0.01 M CaCl2 (soil : solution, 1 : 5), 0.05 M EDTA (soil : solution, 1 : 5), and DTPA (0.005 M DTPA, 0.01 M CaCl2, and 0.1 M TEA at pH 7.3) (soil : solution, 1 : 2.5). The soil suspensions were shaken at 200 rpm for 2 h, centrifuged at 4,000 ×g for 10 min, and filtered [2325]. The Cd, Cr, Cu, Pb, and Zn contents of the bioavailable fraction in each filtrate were determined by flame atomic absorption spectrometry.

2.3. Enzyme Activity Determination

We chose acid phosphatase and dehydrogenase as the subjects for our investigation of enzyme changes under heavy metal stress because they are closely related to the soil phosphorous cycle and to microbial metabolic activities [26, 27]. To determine acid phosphatase activity, dry soil was incubated in disodium phenyl phosphate for 2 h, and the results are expressed as the micrograms of P2O5 released per 100 g of dry soil using 4-aminoantipyrine colorimetry at 510 nm. To measure dehydrogenase activity, dry soil was incubated in triphenyltetrazolium chloride (TTC) for 24 h, and the results are expressed as the micrograms of triphenyl formazan (TPF) released per 1 g of dry soil. The concentration of TPF in the extract was measured using a colorimeter at 485 nm [27, 28].

2.4. Microbial Biomass Carbon

The amount of microbial biomass carbon (biomass C) in the soil was measured by a modified fumigation extraction (FE) method [29]. Soil samples were divided into two parts: one part was fumigated with ethanol-free chloroform for 24 h at 25°C in the dark, and the second part was stored in the dark for 24 h at 25°C but not fumigated. The organic carbon was extracted from both soils using 0.5 M K2SO4 (soil : solution, 1 : 4). Soil extractions were shaken at 150 rpm for 30 min and filtered, and then the total organic carbon was determined using a TOC analyzer (TOC-102A, Analytik Jena). The microbial biomass carbon (MBC) was calculated using the following equation: The microbial biomass was expressed as milligrams of biomass C per kilogram of soil.

2.5. Soil Protein Extraction

The citrate-SDS sequential extraction method [15, 16] was used to extract the soil proteins. Specifically, 3.0 g of dried soil was mixed with 15 mL of 0.25 M citrate buffer (pH 8.0), and the homogenate was shaken at 1,200 rpm in room temperature for 4 h. Then, the suspension was centrifuged for 15 min at 15,000 ×g and 4°C and filtered through filter paper (approximately 30 to 50 μm). Next, the soil was extracted using 15 mL of SDS buffer, which contained 1.25% (w/v) SDS, 0.1 M Tris-HCl (pH 6.8), and 20 mM dithiothreitol (DTT). This SDS soil mixture was shaken for 1 h at 1,200 rpm and room temperature, and then it was centrifuged for 15 min at 15,000 ×g and 4°C. For protein recovery, both the citrate extract and the SDS extract were mixed with buffered phenol (pH 8.0) at a volume ratio (extract : phenol) of 3 : 1 and centrifuged at 15,000 ×g and 4°C for 30 min. After centrifugation, the phenol phase was precipitated at −20°C overnight with five volumes of cold 0.1 M ammonium acetate dissolved in methanol. The proteins were recovered by centrifugation at 20,000 ×g and 4°C for 20 min. The pellets were washed once with cold methanol and twice with cold acetone, and then they were air-dried for further use.

2.6. Protein Content Determination and Protein Separation

Protein pellets were solubilized in 500 μL of lysis buffer, which contained 9 M urea, 4% w/v CHAPS, 1% w/v DTT, 0.5% ampholyte, and 1 mM PMSF. The concentration of protein in the supernatant was determined by the Bradford method [30], and the protein was stored at −70°C.

Samples of the extracted proteins were added to an equal volume of loading buffer, which contained 100 mM Tris-HCl (pH 6.8), 4% (w/v) SDS, 20% glycerol, 0.5% (w/v) bromophenol blue, and 100 mM DTT, and then they were heated in water at 95°C for 5 min prior to SDS-PAGE [8]. Subsequently, discontinuous SDS-PAGE was performed using the Mini-PROTEAN 3 Electrophoresis Cell (Bio-Rad) with a 4% stacking gel and a 12% separating gel. The process was run at a constant 75 V/gel through the stacking gel and a constant 150 V/gel through the separating gel. A prestained protein ladder (approximately 10 to 170 kDa, Fermentas) was loaded as a molecular weight marker, and each lane was loaded with the same quantity of protein. After separation, the gels were stained by silver staining [31] for further comparisons.

2.7. DNA Extraction

The extraction of the total soil DNA was accomplished using a Sangon DNA isolation kit, according to the manufacturer’s protocol (SK8233).

2.8. PCR-DGGE Microbial Community Analysis

The primers F357 and R518, 5′-CCT ACG GGA GGC AGC AGC-3′ and 5′-ATT ACC GCG GCT GCT GG-3′ [20], respectively, were used in this study for the amplification of bacterial 16S rDNA genes. The primers Com2xf and Ac1186r-(pH), 5′-AAA CTC AAA GGA ATT GAC GG-3′ and 5′-CTT CCT CCG AGT TGA CCC-3′ [32], respectively, were used for the amplification of actinomycete 16S rDNA genes. A GC clamp (CGC CCG CCG CGC CCC GCG CCC GGC CCG CCG CCC CCG CCC C) was added to the forward primers to facilitate the DGGE [20]. The PCR reaction mixture (Takara) contained 5 μL of 10x reaction buffer (0.1 M Tris-HCl at pH 8.3 and 0.5 M KCl), 6 μL of MgCl2 (25 μM), 1 μL of each primer (20 μM), 4 μL of each dNTP (2.5 mM), 0.5 μL of Taq DNA polymerase (5 U·μL−1), and 2 μL of DNA; milli-Q water was added to reach a total reaction volume of 50 μL. The PCR protocol used to amplify the soil bacterial and actinomycete 16S rDNA gene fragments was as follows: a 5 min initial denaturation step at 94°C, followed by 20 cycles of 94°C for 1 min, 65°C for 1 min, and 72°C for 1 min, then, 10 cycles of 94°C for 1 min, 55°C for 1 min, and 72°C for 1 min, with a final extension at 72°C for 7 min (S1000 Thermal Cycler, Bio-Rad). Prior to DGGE, the PCR products were checked by electrophoresis in 1.0% (w/v) agarose gels stained with ethidium bromide.

For the DGGE analysis, the PCR products generated from each sample were separated on an 8% acrylamide gel using the Bio-Rad Dcode System with a linear denaturant gradient range from 30% to 60%. DGGE was performed in 1x TAE buffer at 60°C and 150 V for 8 h. Gels were stained with SYBR-Green I, and the gels were scanned (Bio-Rad).

2.9. Analysis of DGGE Profiles

The digitized DGGE images were analyzed by Quantity One image analysis software (Version 4.62, Bio-Rad). The Shannon index was used to estimate the soil microbial diversity based on the intensity and number of bands using the following equation: where is the number of bands in each lane of the DGGE gel, is the peak height of band , and is the sum of the peak heights in a given lane of the DGGE gel [33].

Principal component analysis (PCA) was performed by SPSS 16.0 [34], and WPGAMA cluster analysis was performed by Quantity One 4.62 to calculate the similarities of the gel patterns.

2.10. Data Analysis

All measurements of soil pH, enzyme activity, and levels of organic matter, heavy metals, microbial biomass carbon, and proteins were performed in triplicate. All of the values reported are the average of three determinations. All of the statistical analyses were performed using SPSS software version 16.0. One-way analysis of variance (ANOVA) was used for statistical comparisons, and the Pearson coefficient was used for correlation analysis. A value of was considered statistically significant.

3. Results

3.1. Soil Properties

Selected soil physicochemical properties, the total heavy metal contents, and the bioavailable fraction of heavy metals are shown in Table 1. Compared with the contaminated area, soil samples from the reference site were mildly more acidic, with a pH range from 5.72 to 5.94, and had a lower organic matter content between 4.69% and 5.36%. In the field contaminated by heavy metals, the soil pH varied from 6.23 to 8.00, and the range of soil organic matter was from 3.23% to 7.65%. The amount of organic matter in each sample was lower the further the sample was from the smelter.

In all of the soil samples, the Cr concentration was either below or only slightly above the natural background concentration of 90 mg·kg−1, as defined by the Chinese Environmental Quality Standards for Soil. In addition, the bioavailable Cr fraction was not detected by the apparatus. Therefore, the effects of Cr stress on soil microorganisms were not considered.

The soils varied greatly in their concentrations of total Cu (approximately 175 to 1,147 mg·kg−1), Zn (approximately 187 to 2,472 mg·kg−1), Pb (approximately 77 to 1,199 mg·kg−1), and Cd (approximately 2.54 to 7.26 mg·kg−1) with their distances from the smelter. The highest concentrations of Cu, Zn, and Pb were found in the samples nearest to smelter (200 m), but the highest Cd level was found in the sample 400 m away (Table 1).

3.2. Concentration of Bioavailable Heavy Metals

The fraction of metal, that is, in its bioavailable form, is crucial to understanding metal ecological toxicity. The bioavailability of heavy metals is related to their chemical forms in soil. Several fractions or compartments of the soil act as reservoirs for bioavailable metals [3]. Extractions by CaCl2, DTPA, and EDTA are widely used methods of soil analysis, providing an operationally defined soil compartment that is characterized by its mobility and bioavailability. CaCl2 is used to estimate the soluble and exchangeable metals. DTPA and EDTA are predominantly used to study the metals that have bonded with organic matter and the overall phytoavailability of heavy metals [18, 19].

The concentrations of bioavailable Cu, Zn, Pb, and Cd in the soil samples varied significantly (Table 1). The amount of bioavailable heavy metals in the soil decreased with increasing distance from the smelter, and the bioavailable fractions of heavy metals were significantly positively correlated with the total heavy metals (). In all of the soil samples, the concentration of CaCl2-extractable metals < DTPA-extractable metals < EDTA-extractable metals. These results correlated with the extraction abilities of the reagents [35]. The concentration of CaCl2-extractable heavy metals was far less than those of the DTPA- and EDTA-extractable heavy metals. This result indicates that heavy metals bonded with organic matter and were predominantly in their bioavailable form, agreeing with the conclusion that heavy metals in soil are primarily associated with soil organic matter [36, 37]. It is this bioavailable fraction of heavy metals that greatly affects soil microorganisms [1].

3.3. Soil Enzyme Activities

The variations in the activities of acid phosphatase and dehydrogenase in the soil samples are shown in Figure 1. Good correlations were observed between enzyme activity and distance from the smelter. Compared with the reference area (CK), acid phosphatase activity was inhibited in the sites contaminated by heavy metals. Soil acid phosphatase activity tended to decrease with increasing heavy metals. A significant decrease was observed in the samples only 200 m away from the smelter. Conversely, dehydrogenase activity was slightly increased in the soil samples 800 m away from the smelter and decreased as the distance from the smelter decreased. This result indicates that a low level of heavy metals will stimulate soil microbes to synthesize dehydrogenase, promoting microbial metabolic activity.

Soil acid phosphatase activity was significantly negatively correlated with EDTA-extractable Cu and Pb; DTPA-extractable Zn and Pb; and total Cu. Similarly, dehydrogenase activity was significantly negatively correlated with CaCl2-extractable Cu and Zn; EDTA-extractable Cu, Zn, and Pb; DTPA-extractable Cu, Zn, and Pb; and total Cu, Zn, and Pb (Table 2).

3.4. Microbial Biomass Carbon in Soil

The amount of microbial biomass carbon (MBC) in the contaminated soil samples ranged from 44.6 to 105.0 mg·kg−1, as measured by the FE method, which are lower values than those of the reference site (Figure 2). A significant decrease was observed in the soil samples only 200 m away from the smelter, with the MBC increasing from the plot nearest to the plot farthest away (800 m). The amount of soil MBC was significantly negatively correlated with CaCl2-extractable Cu and Zn; EDTA-extractable Cu, Zn, Pb, and Cd; DTPA-extractable Cu, Zn, Pb, and Cd; and total Cu, Zn, Pb, and Cd (Table 2). These results indicate that the concentrations of DTPA- and EDTA-extractable heavy metals are better predictors of the effects that heavy metals have on MBC compared with the concentrations of CaCl2-extractable heavy metals.

3.5. Soil Microorganism Community Composition

The idea of using a DGGE analysis of soil microbial communities to investigate soil microbial diversity was first introduced by Muyzer et al. (1993). DGGE proved to be capable of accurately measuring soil microbial diversity, so it has become the predominant technique in soil microbial ecology for profiling soil microbial communities and monitoring shifts in microbial community composition due to environmental changes [38]. The DGGE profiles of the soil bacteria and actinomycetes are shown in Figures 3(a) and 4(a). Although many bands were found in all lanes, obvious differences in the intensities of the soil microbial communities were clearly observed (Figures 3(a) and 4(a)). As shown in Figure 3(a), the number of bands in the soil was significantly increased with an increasing distance from the smelter; in other words, the number of bands in the DGGE profile increased as the heavy metals presence decreased. A similar result was also found in the soil actinomycete community composition. Compared with the soils from the reference site and from 800 m away, the number of bands in the plots situated at 400 m and 200 m were significantly decreased (Figure 4(a)).

Subsequently, the DGGE gels were interpreted by the Shannon index, principal component analysis (PCA), and cluster analysis to examine the number, presence, and relative intensity of the bands. The Shannon index indicated that the diversity of the soil microorganisms (bacteria and actinomycetes) was decreased significantly with an increase in the heavy metals concentration (Table 3). PCA and the cluster analysis are also useful methods for analyzing DGGE profiles. The PCA plots show a clear separation due to the different distances from the smelter, demonstrating the altered structure and diversity of the microbial community. Using the presence of the bands as input data, the first two principal components (PC1 and PC2) were sufficient to explain 79.5% and 78.0% of the variance in the soil bacteria and the actinomycetes, respectively (Figures 3(b) and 4(b)). The similarity between the soil bacterial and actinomycete communities of the reference site and those of the soil taken around the smelter was obtained from a cluster analysis of the DGGE profile and was determined to be less than 50% (Figures 3(c) and 4(c)). The cluster analysis indicated that the soil microorganism community structure was greatly affected by the presence of heavy metals.

3.6. Soil Protein Expression

The amount of protein in the contaminated soil samples ranged from 15.6 to 104 μg·g−1, far less than the amount found at the reference site (Figure 5(a)). The concentration of protein decreased drastically when the soil was slightly contaminated by heavy metals at a distance of 800 m from the smelter. The quantity of soil proteins was significantly negatively correlated with soil CaCl2-extractable Cu; EDTA-extractable Cu, Zn, Pb, and Cd; DTPA-extractable Cu, Zn, Pb, and Cd; and total Cu, Zn, Pb, and Cd (Table 2).

The protein profiles are shown in Figure 5(b). The results indicate that no discernible differences exist between the composition and location of the citrate-extracted protein band of the reference site and those of the soil samples from 800 m and 400 m. The bands from these three sample sets were almost entirely located between approximately 55 and 70 kD. However, some type of protein, that is, 40 kD in size, was only found in the soil 200 m from the smelter. Additionally, although the SDS buffer extracts showed similar protein patterns between different soil samples, the expression of some proteins was significantly different. Compared with the contaminated soil samples, the uncontaminated samples had more abundant proteins overall. In the contaminated soil, with its increased concentration of heavy metals, there were fewer large molecular weight (>35 kD) proteins and more low molecular weight proteins (approximately 15 kD). At the base of each contaminated profile was a large agglomeration of presumably low molecular weight proteinaceous material that produced dark areas on the gels. Therefore, these low molecular weight proteins were largely aggregated.

4. Discussion

The activity and community composition of soil microbes are closely related to soil fertility and environmental quality. Additionally, soil biological parameters may have potential for use as early and sensitive indicators of soil ecological stress and restoration. Currently, the amount of microbial biomass carbon and enzyme activity are the metrics predominantly used to provide information about the biochemical processes occurring in the soil [3, 39]. The soil microbial biomass is the total mass of microorganisms living in the soil, which is the living part of the soil organic matter [1]. Soil phosphatase plays an essential role in the mineralization of organic phosphorus, while dehydrogenase is an intracellular enzyme that is involved in microbial oxidoreductase metabolism [3, 39].

In our study, the amount of microbial biomass C in the soil decreased when the concentration of heavy metals increased, although the difference was not significant when the heavy metals were at a low level. This finding is not surprising. Microorganisms differ in their sensitivity to metal toxicity, and the development of metal-tolerant strains, as well as shifts in the community structure, could compensate for the loss of more sensitive populations [2]. The reduction in soil microbial biomass C might also explain the inhibition of enzyme activity observed in the heavy metal contaminated soil.

It is important to differentiate between extracellular and intracellular enzymes when we study the changes of enzymes under pollutant stress. The dehydrogenase activity, which is only present in viable cells and essentially depends on the metabolic state of the soil biota, may therefore be considered a direct measure of soil microbial activity. In contrast, phosphatase activity can occur extracellularly as well as within a living cell [1, 3, 6, 26, 39]. For example, Brookes et al. [40] reported less dehydrogenase activity in metal-contaminated soil than in similar uncontaminated soil, while the soil phosphatase was unaffected. This might explain why, in our study, the dehydrogenase activity was also more sensitive to the heavy metal stress than the acid phosphatase activity.

Furthermore, the dehydrogenase activity was significantly correlated with soil microbial biomass C (, ). These results suggest that microbial biomass C and dehydrogenase can be useful measures of the level of heavy metal contamination in a soil sample [26, 39].

The microbial biomass and enzyme activity measurements have their limitations in soil pollution studies, however. For example, the changes were not notable when heavy metal concentrations were low and highly soluble Cu generated chromogen interference to dehydrogenase activity measurements [41]. In addition, the changes in the microorganism community structure cannot be evaluated by such techniques. Therefore, the use of soil enzyme and microbial biomass as indicators for soil contaminated by heavy metals still needs further discussions.

The changes in the metabolic profiles indicate the possibility that heavy metal contamination results in a community that is more variable and less stable [3]. In our study, the genetic structure and protein expression of the indigenous soil microbial communities were evaluated by an advanced molecular technique that is based on PCR-DGGE and SDS-PAGE. The PCR-DGGE method has become a predominant molecular technique to study the changes in soil microorganisms due to contamination and agricultural practices [3, 20, 42]. In contrast, the use of SDS-PAGE to monitor contaminated soil is still in its infancy [7, 8, 17]. The results of the DGGE profiles, as analyzed by principal component analysis (PCA) and cluster analysis in this study, make it clear that the number and intensity of DNA bands in the different samples changed significantly. This significant difference between the samples indicates that the soil microbial community structure changed greatly in response to heavy metals contamination (Table 3, Figures 3-4). We also detected an increased relative abundance in the microbial populations with increasing distance from the smelter, demonstrating that heavy metals exerted a major effect on bacterial diversity by promoting changes in species composition (represented by the position of the bands) and in species richness (represented by the number of occurring bands). Namely, contamination inhibited certain bacterial groups and stimulated others, and it changed the overall microbial diversity [2, 3, 38, 43, 44].

Subsequently, the content and composition of the soil proteins were studied to investigate their changes under heavy metals stress. Compared with the reference site, the soil protein levels decreased significantly due to heavy metals contamination, and the percentage decrease in total soil proteins was larger than the decrease detected in soil enzyme activity or microbial biomass carbon (Figure 5(a)). The soil protein concentrations were significantly negatively correlated with heavy metals (Table 2). These results indicate that proteins are more sensitive to the changes caused by pollutants, making protein concentration a great potential indicator of soil pollution. The variation in soil proteins may be due to modifications in the microbial community composition [7, 8]. This hypothesis was confirmed by the DGGE profiles in our study (Figures 3(a) and 4(a)). Although protein has potential to be an indictor for soil contamination, its reliability still needs more studies since Bradford method is easily affected by humic substances present in protein solution [45].

The collected soil proteins were separated by SDS-PAGE, and we found that several bands were common to all soils. We also found more presumably low molecular weight proteins in the contaminated soil than in the soil of the reference site (Figure 5(b)). This protein could be related to the presence of heavy metals and could be a microbial response to the metals [8, 46]. For example, many eukaryotic microbes are known to produce metallothioneins (a low molecular weight protein) in response to metal (particularly Cd) exposure [47, 48]. Another explanation may be that sufficient metal exposure will result in the immediate death of cells due to a disruption of their essential functions; thus, the cells are lysed, and their proteins are released into the environment, where they are quickly degraded into low molecular weight materials [2, 8]. Obviously, more work is needed to identify this low molecular weight protein and understand its role in the microbial responses to heavy metals stress. This protein should also be investigated as an indicator for soil pollution.

Although PCR-DGGE and SDS-PAGE are useful for monitoring the changes in microbial community structures, they still have drawbacks. These approaches are easily interfered with by impurities that are coextracted with the soil DNA and proteins, primarily humic substances [7, 38]. The DNA and proteins that are mixed in with the organic matter are not easily separated by DGGE or SDS-PAGE. In addition, relatively small populations of microbes and low abundances of proteins may not be identified by DGGE and SDS-PAGE, respectively [8, 38]. In this study, for example, both the DGGE and SDS-PAGE profiles have a dark background and smeared dark areas at the bottom of SDS-PAGE gel, which make the accuracy and reliability of the analysis questionable. Consequently, the molecular techniques used in this study have much room for improvement, and more work is required to develop effective extraction methods for soil DNA and proteins.

No single technique can provide a comprehensive depiction of the soil microbial situation. Therefore, soil microbial activity and diversity are difficult to elucidate. To obtain a better understanding of soil microbial ecology, we need to integrate different methods to create a comprehensive analysis.

5. Conclusions

By combining different monitoring techniques for different aspects of microbiology, we can obtain a better understanding of the soil microorganisms that live in soils contaminated by heavy metals. Our study demonstrated that heavy metals have a significant negative effect on the activities of soil acid phosphatase and dehydrogenase as well as on the levels of proteins and microbial biomass carbon. Compared with the other potential contamination markers tested, the protein level showed the most dramatic decrease in slightly contaminated soil, indicating that it may be more sensitive to heavy metals. Subsequently, the denaturing gradient gel electrophoresis analysis used in this study demonstrated that heavy metals had a significant negative effect on soil microbial diversity and community structure. Furthermore, the soil protein expression was similar in different soils, but a large quantity of presumably low molecular weight protein was observed only in heavy metal contaminated soil. Based on the research described in this paper, we can conclude that the soil protein level has great potential to be a sensitive indicator of soil contamination. Further research is essential to identify the low molecular weight protein that only appears in contaminated soil.

Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (11179025), National Natural Science Foundation for the Youth of China (21007055), Science Foundation for Postdoctor of China (20100471714, 201104718), and Program for New Century Excellent Talents in University (NCET-11-0455).