Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Activity and functional diversity of microbial communities in long-term hydrocarbon and heavy metal contaminated soils

Activity and functional diversity of microbial communities in long-term hydrocarbon and heavy... The impacts of long-term polycyclic aromatic hydrocarbons (PAHs) and heavy metal pollution on soil microbial communities functioning were studied in soils taken from an old coke plant. The concentrations of PAHs in the tested soils ranged from 171 to 2137 mg kg-1. From the group of tested heavy metals, concentrations of lead were found to be the highest, ranging from 57 to 3478 mg kg-1, while zinc concentrations varied from 247 to 704 mg kg-1 and nickel from 10 to 666 mg kg-1. High dehydrogenase, acid and alkaline phosphatase activities were observed in the most contaminated soil. This may indicate bacterial adaptation to long-term heavy metal and hydrocarbon contamination. However, the Community Level Physiological Profiles (CLPPs) analysis showed that the microbial functional diversity was reduced and influenced to a higher extent by some metals (Pb, Ni), moisture and conductivity than by PAHs. Introduction Due to progressive industrialization polycyclic aromatic hydrocarbons (PAHs) and heavy metals have been excessively released into the environment thus creating a major global concern. PAHs are ubiquitous contaminants known for their carcinogenic, mutagenic and teratogenic properties (Wang et al. 2010, Veses et al. 2013, Rachwal et al. 2015). These hydrophobic compounds are among the persistent organic pollutants, some of them are listed as US-EPA and European priority contaminants. They reveal high solid-water distribution ratio, low aqueous solubility and low bioavailability (Bamforth and Singleton 2005, Lu et al. 2011). In soils, PAHs diffuse into nanopores of soil particles and bind into organic matter. This natural process, called ageing, strongly affects biological toxicity and biodegradability of hydrocarbons (Tang et al. 2012, Cébron et al. 2013). Hydrocarbon contaminants often co-exist with heavy metals, which is the effect of their common source of pollution. Similarly to PAHs heavy metals, present at elevated concentrations, might decrease the number of soil microorganisms, their activity and biodiversity (Klimek and Nikliska 2007, Lu et al. 2013). It has been reported that the combined effects of organic and metallic contaminants on soil biology may be more toxic than those of PAHs alone (Thavamani et al. 2012, Lu et al. 2013). Since the bioremediation relies on the degradation of hydrocarbons by soil microorganism the estimation of the impact of long- term mixed pollution on the activity and functional diversity of microbial communities has considerable benefits for risk assessment and remediation strategies. The understanding of how autochthonic microbial populations respond to mixed contamination present at various concentrations is still limited. Therefore, the aims of the present study were to determine (i) the chemical profile of PAHs and heavy metals in long-term contaminated soils, (ii) the effects of long-term organic and metallic co-contamination on soil microbial activity by assessing the activity of several soil enzymes (dehydrogenase, urease, acid and alkaline phosphatase), and the functional diversity using the Biolog method. Materials and methods Site description and soil sampling The contaminated soil samples used in this study were collected from a 120-year-old coke oven plant located in Upper Silesia, Poland. Coke production at this site began in 1884. In preliminary investigations concentration of PAH and heavy metal in several samples taken from the coke oven territory and the adjacent area (50°19N, 18°47E) were determined. These analyses revealed that PAHs contamination was unevenly distributed across the tested area. On the basis of different levels of PAH and heavy metal contamination four sample subareas were selected. The soil samples were taken from the upper horizon (0­20 cm depth), and each soil sample was were analysed by means of PCA performed on the correlation matrix of the values of area under the color development profile. a composite of 5 mixed and homogenized subsamples collected at the four corners and the centre of an area of 10 × 10 m. The samples were air-dried and sieved (mesh size 2 × 2 mm) to remove the coarse fraction. Determination of the soil properties The physicochemical properties: pH in water, conductivity, soil total N, P and soil particle size distribution (size classification according to PN-R-04033 1998) were determined using the ISO standards (ISO 11263 1994, ISO 11261 1995, ISO 10390 1997, ISO 11265 1997, ISO 11277 1998). Organic matter (OM) content was determined by sulfochromic oxidation of organic carbon followed by titration of the excess K2Cr2O7 with FeSO4(NH4)2SO4·6H2O as described by Schinner et al. (1996). The 16 EPA priority PAHs in contaminated soils were identified and determined according to the ISO (ISO 13877 1998) and the EPA Standard Methods (USEPA 1992, 1995). PAHs were detected by high-performance liquid chromatography coupled with a fluorescence detector (HPLC-FLD Agilent 1200). Heavy metal contents (Cd, Pb, Ni, Zn, Cu) were determined after mineralization with aqua regia using atomic absorption spectrometry (SpectrAA 300 Varian) with graphite furnace. The limit of detection and recovery values of tested parameters are presented in Table 1. Enzyme assays The microbial activities, in terms of soil enzymes, were measured using the following methods. Dehydrogenase activity (DHA) was determined using triphenyltetrazolium chloride as an artificial electron acceptor and expressed as g triphenyl formazone (TPF) g-1 dry soil 24 h-1 (Alef 1995). Soil urease activity was expressed as g N-NH4 g-1 dry soil h-1 (Tabatabai and Bremner 1972). The activities of alkaline and acid phosphatases were expressed as mg p-nitrophenol g-1 dry soil h-1 (Tabatabai and Bremner 1969). Community-level physiological profiles (CLPPs) Community level physiological profiles (CLPPs) were assessed using the Biolog EcoPlateTM system (Biolog Inc., CA, USA). Briefly, microorganisms were extracted from soil suspensions (10 g dry weight of soil in 90 ml sterile 0.85% NaCl solution) and aliquots of 125 l were inoculated on plates, which were incubated at 24°C in the dark. The readings were taken at 590 nm after inoculation and at 12 h intervals for 168 h using a microplate reader (VictorTM X5 Multilabel Plate Readers, PerkinElmer). The absorbance measurements for individual substrates were corrected as described by Azarbad et al. (2013). The corrected absorbance values were used to calculate the following microbial indices: substrate richness (Rs), functional diversity based on the Shannon-Wiener biodiversity index (H') and evenness (I') (Klimek and Nikliska 2007). The calculated indices were used in statistical analyses. Statistical analysis Statistical analysis was performed using STATISTICA 10.0 PL software (StatSoft, Tulsa, USA). Differences among samples were analysed by one-way analysis of variance (ANOVA), followed by a post hoc least significant difference test in order to discriminate among the means. The Pearson's correlation coefficients of the data were calculated to determine the relationships among the parameters at the significance levels p<0.05 and p<0.01. The physicochemical and biological results Results and discussion Soils contamination with PAHs and metals The total concentrations of PAHs in analysed soils varied from 171 to 2137 mg kg-1 dry soil weight (Table 2). The results of PAHs content in tested soils indicated that the contamination was not only unevenly distributed across the site but also very high. The most contaminated soil II exhibited extremely high concentration of naphthalene and the highest level of organic matter (Table 2, 3). This study revealed a particularly strong link (r=0.99, p<0.01) between the OM and naphthalene (Table 4). This fact is in agreement with the findings of Cai et al. (2007) and Nam et al. (2008), who suggest that low molecular weight (LMW) PAHs (with 2 or 3 benzene rings) might more easily approach equilibrium with soil OM over a larger spatial scale (Cai et al. 2007, Nam et al. 2008). Moreover, it has been mentioned that the determination of OM may be overestimated due to the inclusions of hydrocarbons bound tightly to soil particles (Cai et al. 2007, Nam et al. 2008). In the present study, we observed that the high molecular weight (HMW) compounds (with 4 or more benzene rings) were the predominant PAHs found in soils I, III and IV. In the case of soil II they accounted for 47% of all hydrocarbons. The high level of HMW PAHs in comparison with LMW in soils may result from the different physicochemical properties of these hydrocarbon compounds. The HMW PAHs bind more easily to the soil particles and are less water-soluble thereby these pollutants become less bioavailable and biodegradable by microorganisms (Johnsen and Karlson 2007, Lladó et al. 2009). Total concentrations of seven carcinogenic PAHs (Table 2) varied between 21% and 36% and were much lower than previously reported (50% 16WWA) in the Silesia area (Maliszewska-Kordybach et al. 2009). The presence of heavy metals was also observed in the tested soils (Table 3). Soils III and IV had the highest amount of Pb, Ni and Cu whose concentrations reached the values of 3378, 519, 143 and 3478, 664, 169 mg kg-1 dry soil, respectively. In most soil samples, the concentrations of lead and nickel exceeded the Polish regulation guidelines (Dz.U. 2002). It is well known that the bioavailability of metals is dependent on pH (Vig et al. 2003). In the tested soils, the pH was neutral or slightly alkaline (Table 3) which indicates that heavy metals were most likely in an inaccessible form and hence did not affect the microbial communities (Ross 1994, Antoniadis et al. 2008). This is supported by the observation that no negative correlation between enzymatic activity and heavy metals content was noted (Table 4). In our study a strong positive correlation (p < 0.01) between the sum of 16 PAHs and cadmium (r=0.74) and zinc (r=0.82) was observed (Table 4), pointing to emission from common sources. Surprisingly, there were no statistically significant correlations between PAHs and other metals: Pb, Cu and Ni. This fact is difficult to explain, and may suggest other source of contamination than coke production, especially since these metals have been strongly correlated to each other. Nevertheless, PAHs and heavy metals in long-term contaminated soils could have a strong negative impact on microbial diversity and activity (Maliszewska-Kordybach et al. 2009, Thavamani et al. 2012, Pessacq et al. 2015). Table 1. The limit of detection and recovery values of tested parameters Parameter Cd Pb Zn Cu Ni Naphthalene Acenaphthene Fluorene Phenanthrene Anthracene Fluoranthene Pyrene Benzo(a)anthracene Chrysene Benzo(b)fluoranthene Benzo(k)fluoranthene Benzo(a)pyrene Benzo(ghi)perylene Dibenzo(ah)anthracene Indeno(1,2,3-cd)pyrene Limit of detection [mg kg-1] 0.87 3.39 8.73 4.34 4.16 0.066 0.01 0.003 0.008 0.002 0.016 0.011 0.006 0.006 0.007 0.003 0.004 0.004 0.002 0.004 Recovery [%] 84.8 95.7 98.8 95.0 95.3 >70.0 >70.0 >70.0 >70.0 >70.0 96.1 95.4 93.3 95.3 94.4 93.2 88.8 90.4 92.6 90.3 Table 2. The concentrations of individual and total PAHs in tested soilsa PAHs [mg kg-1] Naphthalene Acenaphthene Fluorene Phenanthrene Anthracene Fluoranthene Pyrene Benzo(a)anthracene* Chrysene* Benzo(b)fluoranthene* Benzo(k)fluoranthene* Benzo(a)pyrene* Benzo(ghi)perylene Dibenzo(ah)anthracene* Indeno(1,2,3-cd)pyrene* Acenaphthylene 16 PAHs SOILS I 24.6 ± 4.1 8.2 ± 1.1 4.9 ± 0.6 25.0 ± 2.5 7.7 ± 0.7 29.9 ± 3.0 18.0 ± 1.8 10.4 ± 0.8 6.8 ± 0.6 7.3 ± 0.7 3.9 ± 0.3 9.7 ± 0.8 4.4 ± 0.4 2.2 ± 0.2 8.3 ± 0.8 ­ 171.4 II 476.6 ± 78.6 54.1 ± 7.3 85.6 ± 11.1 392.8 ± 39.3 119.7 ± 11.4 322.0 ± 32.2 189.3 ± 18.9 101.5± 8.1 68.5 ± 5.5 62.9 ± 6.3 37.1 ± 3.2 88.2 ± 7.1 42.1 ± 3.8 20.5 ± 1.7 76.0 ± 6.8 ­ 2137.1 III 14.0 ± 2.3 5.3 ± 0.7 20.9 ± 2.7 92.5 ± 9.3 21.3 ± 2.0 78.0 ± 7.8 58.3 ± 5.8 36.9 ± 3.0 37.5 ± 3.0 25.2 ± 2.5 13.4 ± 1.1 34.3 ± 2.7 14.1 ± 1.3 1.2 ± 0.1 24.2 ± 2.2 4.7 ± 0.5 481.6 IV 37.8 ± 6.2 4.2 ± 0.6 21.7 ± 2.8 95.3 ± 9.5 18.1 ± 1.7 166.8 ± 16.7 63.6 ± 6.4 37.6 ± 3.0 25.7 ± 2.1 24.5 ± 2.5 12.7 ± 1.1 32.9 ± 2.6 18.3 ± 1.6 4.2 ± 0.4 31.8 ± 2.9 9.4 ± 1.0 604.6 * Carcinogenic PAHs a Data are presented as average ± standard deviation, n = 3. Table 3. The physical and chemical properties of tested soilsa Parameter Sand % Silt % Clay % pH Moisture % OM % Total N % Total P % Conductivity S cm Cd mg kg soil -1 -1 SOILS I 95 5 0 7.2 ± 0.1 20.1 ± 1.0 2.19 ± 0.12 5.3 ± 0.1 0.01 ± 0.00 367.3 ± 18.4 2.6 ± 0.4 57.8 ± 2.9 247 ± 12.4 11.8 ± 0.5 10.5 ± 0.4 II 80 2 18 8.3 ± 0.2 21.7 ± 1.1 8.79 ± 0.81 0.09 ± 0.00 0.01 ± 0.00 214.7 ± 10.7 3.1 ± 0.5 1843 ± 92.2 704 ± 35.2 99.7 ± 4.0 184.2 ± 7.4 III 98 2 0 7.5 ± 0.1 14.1 ± 0.7 1.95 ± 0.53 0.1 ± 0.00 1.24 ± 0.03 95.0 ± 4.8 3.1 ± 0.5 3378 ± 168.9 556 ± 27.8 143.8 ± 5.8 519.7 ± 20.8 IV 96 4 0 7.3 ± 0.1 13.4 ± 0.7 2.25 ± 0.3 0.17 ± 0.00 1.16 ± 0.03 248.1 ± 12.4 2.8 ± 0.5 3478 ± 173.9 504 ± 25.2 169.8 ± 6.8 666.4 ± 26.7 Pb mg kg soil -1 Zn mg kg soil -1 Cu mg kg-1 soil Ni mg kg-1 soil Data are presented as average ± standard deviation, n = 3. Enzyme activity The activity of soil enzymes is often used to monitor the effects of various pollutants on microbial functioning (Turgay et al. 2010, Tang et al. 2014, Tomkiel et al. 2015). In the tested soils, the activities of a range of enzymes contributing to the biogeochemical cycle of carbon, nitrogen and phosphorus were determined (Fig. 1A­C). In addition, the activity of dehydrogenase, present only in viable cells was assessed. DHA is linked to microbial respiration and reflects the ability of microorganisms to degrade hydrocarbons (Achuba and Okoh 2014, Wang et al. 2014). In the tested soils, despite the presence of strong PAHs and metal co-contamination, a significant activity of the investigated enzymes was observed. The most contaminated soil exhibited a particularly high DHA as well as significant acid and alkaline phosphatase activity (Fig. 1A­B). Moreover, we observed strong positive correlations (p<0.01) between the activity of these enzymes and the total PAHs (Table 4). Several studies have demonstrated that DHA is among the most sensitive parameters for evaluation of hydrocarbons toxicity and stress (Irha et al. 2003, Gianfreda et al. 2005). Turgay et al. (2010) observed that the decrease in the content of hydrocarbons resulted in an increase of DHA over time during the bioremediation process of soil contaminated with petroleum hydrocarbons. However, Tang et al. (2012) examining the bioavailability of petroleum hydrocarbons under different attenuation conditions noticed an increase in DHA immediately after the contamination of soil and its decline with the duration of the experiment. Observed in our study high dehydrogenase and other enzyme activities in the most contaminated soil suggested significant intensification of microbial processes. This phenomenon may be a microbial response to a constant stress exerted by contaminants. Moreover, such results may be related to the high level of organic matter and moisture found in the most contaminated soil II. The highest urease activity and the lowest dehydrogenase and acid phosphatase activities were found in soil IV. Soil I, which was the least contaminated one, had the lowest activities of urease and alkaline phosphatase. Contrary to some other studies (Gianfreda et al. 2005, Turgay et al. 2010, Wu et al. 2014) no correlation was observed between urease activity and the level of PAHs. Thus, the variations in enzyme activities between soil samples showed that in tested soils pollution inputs have specific site effects, which may be a result of lower sensitivity of local microbial communities to PAHs and heavy metals. For instance, Qasemian et al. (2012) found an increased activity of acid phosphatase, cellulase and -glucosidase in the leaf litter of forest pine (Pinus halepensis) after 3 months of anthracene contamination. The changes in the activity of an individual enzyme did not allow for a precise quantification of soil degradation by long-time mixed contamination. This finding demonstrates that microorganisms still function when subjected to high levels of pollution. Catabolic diversity of bacterial communities The potential activity of bacterial communities that were able to grow on the substrate provided was measured using Biolog method. Very low metabolic activity with respect to 31 substrates in communities from soils I and II indicated the microbial activity at a basal metabolic rate. It is remarkable that soil III has shown the highest value of average well-color development (AWCD) which expressed the catabolic activity of metabolic communities on EcoPlatesTM but a low enzyme activity. The low activity of tested enzymes may be related to the inhibition effects caused by toxic hydrocarbons. Moreover, soil enzyme activity provides information about the breaking of different chemical bonds whereas the Biolog method reflects the metabolic activity of soil bacteria in response to 31 carbon sources, which were present at high concentrations (Palmroth et al. 2005). Furthermore, CLPPs do not measure the functional abilities of the entire soil microbial community but only that of a very limited subset of microbial fraction capable to grow Table 4. Pearson's correlation coefficients between environmental factors, enzyme activities and functional diversity indices Total P Conductivity Cd Pb Zn Cu Ni 16 PAHs NAH DHA Urease AcdP HMW PAHs Cd pH Moisture OM Total N pH Moisture OM 0.95** Total N -0.51 0.61* -0.32 Total P -0.64* 0.19 0.89** 0.17 0.83** 0.94** -0.41 -0.28 -0.57 0.61* -0.52 -0.25 0.88** 0.90** 0.85** -0.68* 0.20 -0.72** 0.23 0.73** 0.68* -0.71** 0.24 0.73** -0.43 0.79** 0.13 0.78** 0.12 0.10 0.08 -0.14 0.63* -0.12 0.68* -0.18 0.11 0.71** 0.30 0.78** -0.50 0.16 0.63* 0.63* 0.58* -0.12 0.65* -0.15 0.70** -0.08 -0.08 0.63* -0.14 0.70** -0.05 -0.35 -0.38 0.80** -0.33 -0.11 0.73** 0.73** 0.68* -0.31 0.78** 0.19 0.88** 0.28 -0.02 -0.23 0.74** 0.05 0.82** 0.14 -0.16 0.99** 0.98** 0.96** 0.03 0.91** 0.84** -0.43 -0.46 -0.46 0.94** 0.92** 0.15 0.87** 0.83** -0.33 -0.37 0.99** -0.11 0.92** 0.81** -0.56 -0.60* -0.59* -0.20 -0.61* 0.30 0.96** 0.37 0.95** -0.72** 0.56 0.99** 0.63* -0.71** 0.96** 0.59* -0.79** 0.56 -0.83** -0.40 -0.89** -0.57 Conductivity 0.59* -0.07 0.81** Cd 0.81** -0.25 0.61* -0.85** Pb -0.85** -0.15 -0.88** Zn 0.83** -0.26 0.68* -0.90** Cu -0.83** -0.07 -0.90** Ni -0.21 -0.89** -0.36 -0.73** 16 PAHs 0.98** 0.97** -0.52 HMW PAHs 0.96** 0.93** -0.63* NAH 0.96** 0.99** DHA 0.97** 0.98** -0.33 Urease -0.12 -0.63* -0.13 AcdP 0.91** 0.89** 0.89** 0.80** -0.52 -0.55 -0.54 -0.15 -0.12 0.28 0.26 0.22 0.90** -0.55 -0.51 -0.45 AlkP 0.88** 0.76** -0.52 Rs -0.78** -0.54 -0.44 H' -0.37 -0.72** -0.43 I' -0.37 -0.64* -0.57 -0.39 Statistically significant correlations: * p < 0.05; ** p < 0.01. NAH ­ naphthalene Fig. 1. Acid and alkaline phosphatase (A), dehydrogenase (B) and urease (C) activity in tested soils on the carbon sources provided (Palmroth et al. 2005, Hueso et al. 2012). To assess the potential functional diversity of the bacterial community H', I' and Rs indices were calculated (Table 5). The functional diversity and evenness indices were the highest in soil III and the lowest in soil I and II. Microbial communities in soil III were characterised by high species richness and an even distribution of functional abilities within the community. Bacterial communities in the most and the least contaminated soils reveal very low functional biodiversity. No correlations were found between the Shannon indices and the sum of 16 PAHs, which is in contradistinction to the findings of Zhang et al. (2010) who reported a correlation between these two factors. However, Muckian et al. (2007) found that the level of PAHs did not alter the number of ribotypes in the contaminated soils. They indicated that the PAH ring structure plays a greater role in influencing the bacterial community than the total PAHs concentration. In our study, a negative correlation between H' and naphthalene, a 2-ring PAH, was observed (Table 4). This can be explained by different properties of these hydrocarbons. Naphthalene is the most soluble among PAHs and thus more mobile and toxic to microorganisms. No correlations were found between functional diversity indices and enzyme activities. The functional diversity of bacterial communities was positively correlated with total phosphorus (r=0.90), lead (r=0.73) and nickel content (r=0.73, p<0.01). There were also observed strong negative correlations between diversity indices and moisture (r=-0.72) and conductivity (r=-0.72, p<0.01). It seems that some metals and physicochemical factors are more important in controlling the soil bacterial diversity than hydrocarbons (Azarbad et al. 2013). The number of utilized carbon sources, which is defined as the richness index, varied depending on the soil, from 26 to 31 out of 31 carbon sources available on the Biolog EcoPlatesTM (Table 5). Two soils exhibited strongly reduced richness indices: the most contaminated (II) one and the least contaminated (I) one. Table 5. Richness (Rs), Shannon­Weaver (H') and evenness (I') indices calculated on carbon substrate used in Biolog EcoPlateTM Index Rs H' I' SOILS I 27 ± 1.2 a 1.30 ± 0.02 a 0.91 ± 0.02 a II 26 ± 0.6 a 1.28 ± 0.01 a 0.90 ± 0.01 a III 31 ± 0.6 b 1.44 ± 0.03 b 0.97 ± 0.01 b IV 29 ± 0.6 c 1.36 ± 0.01 c 0.94 ± 0.01 c Values with the same letter mean indicate that there are no significant differences (p < 0.05) Fig. 2. Principal Component Analysis (PCA) of (A) carbon utilization profiles of bacterial communities (B) environmental factors and enzymatic activities of tested soils. Each point represents the bar center of three replicates Principal component analysis All the obtained data were subjected to principal component analysis (Fig. 2A­B). The first PCA was carried out for carbon utilization profiles. PCA revealed differences in functional diversity between the bacterial communities present in the tested soils. Both the most contaminated (II) and the least contaminated (I) soils exhibited a noticeable separation from the two other soils in the first principal component (PC1), which explained 78.5% of the variation in the data. The reason for these findings may be the selection and adaptation of microorganisms to the pollutants, which cause a lack of species or strains capable to degrade a variety of carbon sources (Wu et al. 2014). The second PCA analysis was carried out for physicochemical data, 16PAHs concentrations and enzymatic activities (Fig 2B). The majority of the variance (94.4%) was explained by the two principal components. Both soils III and IV exhibited a significant separation from the other two soils (I and II) along the first and the second axis. The PC1 was highly negatively correlated with the sum of PAHs, the activity of dehydrogenase, both phosphatases, as well as with pH, whereas the second factor was highly positively correlated with moisture and negatively correlated with total P, Pb, Cu and Ni. The long-term pressure of the mixture of contaminants may have modified the microbial metabolism and/or the relative composition and size of the constituent populations. This is because microbial communities have the ability to adapt, physiologically and evolutionally, to the presence of contaminants in soil. After the period of "stress", there is an increase in respiration intensity and enzyme activity (Baran et al. 2004, Pessacq et al. 2015). This could have occurred only in the presence of acclimated microbial communities, which were capable of adapting and proliferating in the presence of high levels of PAHs and metals. Understanding these complex interactions is of major importance for the bioremediation of contaminated soils (Thavamani et al. 2012). Conclusions Our study has demonstrated that enzyme activities and microbial functional diversity in PAHs and metal contaminated soils is related and characteristic to a specific site. This fact is connected with the high heterogeneity of soil and uneven distribution of hydrocarbons in soil. The obtained results showed also that the evaluation of a single parameter such as the activity of enzyme does not reflect the real effect of contaminants input on the activity and diversity of microbial communities. Moreover, the results of this study clearly showed that microbial communities have enormous adaptation abilities to the mixed inorganic and organic contamination. Acknowledgements We would like to thank the professional team of Analytical Laboratory in IETU for physical and chemical analysis. This work was supported by Grant No. N305 051340 financed by the Polish Ministry of Science and Higher Education. Klimek, B. & Nikliska, M. (2007). Zinc and copper toxicity to soil bacteria and fungi from zinc polluted and unpolluted soils: A comparative study with different types of Biolog plates, Bulletin of Environmental Contamination and Toxicology, 78, pp. 112­117. Lladó, S., Jiménez, N., Viñas, M. & Solanas, A.M. (2009). Microbial populations related to PAH biodegradation in an aged biostimulated creosote-contaminated soil, Biodegradation, 20, pp. 593­601. Lu, M., Xu, K. & Chen, J. (2013). Effect of pyrene and cadmium on microbial activity and community structure in soil, Chemosphere, 91, pp. 491­497. Lu, X-Y., Zhang, T. & Fang, H.H-P. (2011). Bacteria-mediated PAH degradation in soil and sediment, Applied Microbiology and Biotechnology, 89, pp. 1357­1371. Maliszewska-Kordybach, B., Smreczak, B. & Klimkowicz-Pawlas, A. (2009). Effects of anthropopressure and soil properties on the accumulation of polycyclic aromatic hydrocarbons in the upper layer of soils in selected regions of Poland, Applied Geochemistry, 24, pp. 1918­1926. Muckian, L., Grant, R., Doyle, E. & Clipson, N. (2007). Bacterial community structure in soils contaminated by polycyclic aromatic hydrocarbons, Chemosphere, 68, pp. 1535­1541. Nam, J.J., Thomas, G.O., Jaward, F.M., Steinnes, E., Gustafsson, O. & Jones, K.C. (2008). PAHs in background soils from Western Europe: influence of atmospheric deposition and soil organic matter, Chemosphere, 70, pp. 1596­1602. Palmroth, M.R.T., Münster, U., Pichtel, J. & Puhakka, J.A. (2005). Metabolic responses of microbiota to diesel fuel addition in vegetated soil, Biodegradation, 16, pp. 91­101. Pessacq, J., Medina, R., Terada, C., Bianchini, F.E., Morelli, I.S. & Del Panno, M.T. (2015). Assessment of the responsiveness to different stresses of the microbial community from long-term hydrocarbon-contaminated soils, Water, Air, & Soil Pollution, 226, 20. DOI 10.1007/s11270-014-2262-9. PN-R-04033 (1998). Gleby i utwory mineralne. Podzial na frakcje i grupy granulometryczne. Polish Standard. Distribution of the particle sizes in the soil. Polish Committee for Standardization. (in Polish) Qasemian, L., Guiral, D., Ziarelli, F., Van Dang, T.K. & Farnet, A.M. (2012). Effects of anthracene on microbial activities and organic matter decomposition in a Pinus halepensis litter from a Mediterranean coastal area, Soil Biology and Biochemistry, 46, pp. 148­154. Rachwal, M., Magiera, T. & Wawer, M. (2015). Coke industry and steel metallurgy as the source of soil contamination by technogenic magnetic particles, heavy metals and polycyclic aromatic hydrocarbons, Chemosphere, http://dx.doi.org/10.1016/j. chemosphere.2014.11.077 Ross, S. (1994). Retention, transformation and mobility of toxic metals in soils, in: Toxic metals in soil-plant systems, Ross, S. (Ed.), John Woley Sons, Chichester UK 1994. Schinner, F., Öhlinger, R., Kandeler, E. & Margesin, R. (1996). Methods in soil biology. Springer, Heidelberg Berlin, New York 1996. Tabatabai, M.A. & Bremner, J.M. (1969). Use of p-nitrophenyl phosphate in assay of soil phosphatase activity, Soil Biology and Biochemistry, 1, pp. 301­307. Tabatabai, M.A. & Bremner, J.M. (1972). Assay of urease activity in soils, Soil Biology and Biochemistry, 4, pp. 479­487. Tang, X., Hashmi, M.Z., Long, D., Chen, L., Khan, M.I. & Shen, Ch. (2014). Influence of heavy metals and PCBs pollution on the enzyme activity and microbial community of paddy soils around an e-waste recycling workshop, International Journal of Environmental Research and Public Health, 11, pp. 3118­3131. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Archives of Environmental Protection de Gruyter

Activity and functional diversity of microbial communities in long-term hydrocarbon and heavy metal contaminated soils

Loading next page...
 
/lp/de-gruyter/activity-and-functional-diversity-of-microbial-communities-in-long-ZJAewuHd2x
Publisher
de Gruyter
Copyright
Copyright © 2016 by the
ISSN
2083-4810
eISSN
2083-4810
DOI
10.1515/aep-2016-0041
Publisher site
See Article on Publisher Site

Abstract

The impacts of long-term polycyclic aromatic hydrocarbons (PAHs) and heavy metal pollution on soil microbial communities functioning were studied in soils taken from an old coke plant. The concentrations of PAHs in the tested soils ranged from 171 to 2137 mg kg-1. From the group of tested heavy metals, concentrations of lead were found to be the highest, ranging from 57 to 3478 mg kg-1, while zinc concentrations varied from 247 to 704 mg kg-1 and nickel from 10 to 666 mg kg-1. High dehydrogenase, acid and alkaline phosphatase activities were observed in the most contaminated soil. This may indicate bacterial adaptation to long-term heavy metal and hydrocarbon contamination. However, the Community Level Physiological Profiles (CLPPs) analysis showed that the microbial functional diversity was reduced and influenced to a higher extent by some metals (Pb, Ni), moisture and conductivity than by PAHs. Introduction Due to progressive industrialization polycyclic aromatic hydrocarbons (PAHs) and heavy metals have been excessively released into the environment thus creating a major global concern. PAHs are ubiquitous contaminants known for their carcinogenic, mutagenic and teratogenic properties (Wang et al. 2010, Veses et al. 2013, Rachwal et al. 2015). These hydrophobic compounds are among the persistent organic pollutants, some of them are listed as US-EPA and European priority contaminants. They reveal high solid-water distribution ratio, low aqueous solubility and low bioavailability (Bamforth and Singleton 2005, Lu et al. 2011). In soils, PAHs diffuse into nanopores of soil particles and bind into organic matter. This natural process, called ageing, strongly affects biological toxicity and biodegradability of hydrocarbons (Tang et al. 2012, Cébron et al. 2013). Hydrocarbon contaminants often co-exist with heavy metals, which is the effect of their common source of pollution. Similarly to PAHs heavy metals, present at elevated concentrations, might decrease the number of soil microorganisms, their activity and biodiversity (Klimek and Nikliska 2007, Lu et al. 2013). It has been reported that the combined effects of organic and metallic contaminants on soil biology may be more toxic than those of PAHs alone (Thavamani et al. 2012, Lu et al. 2013). Since the bioremediation relies on the degradation of hydrocarbons by soil microorganism the estimation of the impact of long- term mixed pollution on the activity and functional diversity of microbial communities has considerable benefits for risk assessment and remediation strategies. The understanding of how autochthonic microbial populations respond to mixed contamination present at various concentrations is still limited. Therefore, the aims of the present study were to determine (i) the chemical profile of PAHs and heavy metals in long-term contaminated soils, (ii) the effects of long-term organic and metallic co-contamination on soil microbial activity by assessing the activity of several soil enzymes (dehydrogenase, urease, acid and alkaline phosphatase), and the functional diversity using the Biolog method. Materials and methods Site description and soil sampling The contaminated soil samples used in this study were collected from a 120-year-old coke oven plant located in Upper Silesia, Poland. Coke production at this site began in 1884. In preliminary investigations concentration of PAH and heavy metal in several samples taken from the coke oven territory and the adjacent area (50°19N, 18°47E) were determined. These analyses revealed that PAHs contamination was unevenly distributed across the tested area. On the basis of different levels of PAH and heavy metal contamination four sample subareas were selected. The soil samples were taken from the upper horizon (0­20 cm depth), and each soil sample was were analysed by means of PCA performed on the correlation matrix of the values of area under the color development profile. a composite of 5 mixed and homogenized subsamples collected at the four corners and the centre of an area of 10 × 10 m. The samples were air-dried and sieved (mesh size 2 × 2 mm) to remove the coarse fraction. Determination of the soil properties The physicochemical properties: pH in water, conductivity, soil total N, P and soil particle size distribution (size classification according to PN-R-04033 1998) were determined using the ISO standards (ISO 11263 1994, ISO 11261 1995, ISO 10390 1997, ISO 11265 1997, ISO 11277 1998). Organic matter (OM) content was determined by sulfochromic oxidation of organic carbon followed by titration of the excess K2Cr2O7 with FeSO4(NH4)2SO4·6H2O as described by Schinner et al. (1996). The 16 EPA priority PAHs in contaminated soils were identified and determined according to the ISO (ISO 13877 1998) and the EPA Standard Methods (USEPA 1992, 1995). PAHs were detected by high-performance liquid chromatography coupled with a fluorescence detector (HPLC-FLD Agilent 1200). Heavy metal contents (Cd, Pb, Ni, Zn, Cu) were determined after mineralization with aqua regia using atomic absorption spectrometry (SpectrAA 300 Varian) with graphite furnace. The limit of detection and recovery values of tested parameters are presented in Table 1. Enzyme assays The microbial activities, in terms of soil enzymes, were measured using the following methods. Dehydrogenase activity (DHA) was determined using triphenyltetrazolium chloride as an artificial electron acceptor and expressed as g triphenyl formazone (TPF) g-1 dry soil 24 h-1 (Alef 1995). Soil urease activity was expressed as g N-NH4 g-1 dry soil h-1 (Tabatabai and Bremner 1972). The activities of alkaline and acid phosphatases were expressed as mg p-nitrophenol g-1 dry soil h-1 (Tabatabai and Bremner 1969). Community-level physiological profiles (CLPPs) Community level physiological profiles (CLPPs) were assessed using the Biolog EcoPlateTM system (Biolog Inc., CA, USA). Briefly, microorganisms were extracted from soil suspensions (10 g dry weight of soil in 90 ml sterile 0.85% NaCl solution) and aliquots of 125 l were inoculated on plates, which were incubated at 24°C in the dark. The readings were taken at 590 nm after inoculation and at 12 h intervals for 168 h using a microplate reader (VictorTM X5 Multilabel Plate Readers, PerkinElmer). The absorbance measurements for individual substrates were corrected as described by Azarbad et al. (2013). The corrected absorbance values were used to calculate the following microbial indices: substrate richness (Rs), functional diversity based on the Shannon-Wiener biodiversity index (H') and evenness (I') (Klimek and Nikliska 2007). The calculated indices were used in statistical analyses. Statistical analysis Statistical analysis was performed using STATISTICA 10.0 PL software (StatSoft, Tulsa, USA). Differences among samples were analysed by one-way analysis of variance (ANOVA), followed by a post hoc least significant difference test in order to discriminate among the means. The Pearson's correlation coefficients of the data were calculated to determine the relationships among the parameters at the significance levels p<0.05 and p<0.01. The physicochemical and biological results Results and discussion Soils contamination with PAHs and metals The total concentrations of PAHs in analysed soils varied from 171 to 2137 mg kg-1 dry soil weight (Table 2). The results of PAHs content in tested soils indicated that the contamination was not only unevenly distributed across the site but also very high. The most contaminated soil II exhibited extremely high concentration of naphthalene and the highest level of organic matter (Table 2, 3). This study revealed a particularly strong link (r=0.99, p<0.01) between the OM and naphthalene (Table 4). This fact is in agreement with the findings of Cai et al. (2007) and Nam et al. (2008), who suggest that low molecular weight (LMW) PAHs (with 2 or 3 benzene rings) might more easily approach equilibrium with soil OM over a larger spatial scale (Cai et al. 2007, Nam et al. 2008). Moreover, it has been mentioned that the determination of OM may be overestimated due to the inclusions of hydrocarbons bound tightly to soil particles (Cai et al. 2007, Nam et al. 2008). In the present study, we observed that the high molecular weight (HMW) compounds (with 4 or more benzene rings) were the predominant PAHs found in soils I, III and IV. In the case of soil II they accounted for 47% of all hydrocarbons. The high level of HMW PAHs in comparison with LMW in soils may result from the different physicochemical properties of these hydrocarbon compounds. The HMW PAHs bind more easily to the soil particles and are less water-soluble thereby these pollutants become less bioavailable and biodegradable by microorganisms (Johnsen and Karlson 2007, Lladó et al. 2009). Total concentrations of seven carcinogenic PAHs (Table 2) varied between 21% and 36% and were much lower than previously reported (50% 16WWA) in the Silesia area (Maliszewska-Kordybach et al. 2009). The presence of heavy metals was also observed in the tested soils (Table 3). Soils III and IV had the highest amount of Pb, Ni and Cu whose concentrations reached the values of 3378, 519, 143 and 3478, 664, 169 mg kg-1 dry soil, respectively. In most soil samples, the concentrations of lead and nickel exceeded the Polish regulation guidelines (Dz.U. 2002). It is well known that the bioavailability of metals is dependent on pH (Vig et al. 2003). In the tested soils, the pH was neutral or slightly alkaline (Table 3) which indicates that heavy metals were most likely in an inaccessible form and hence did not affect the microbial communities (Ross 1994, Antoniadis et al. 2008). This is supported by the observation that no negative correlation between enzymatic activity and heavy metals content was noted (Table 4). In our study a strong positive correlation (p < 0.01) between the sum of 16 PAHs and cadmium (r=0.74) and zinc (r=0.82) was observed (Table 4), pointing to emission from common sources. Surprisingly, there were no statistically significant correlations between PAHs and other metals: Pb, Cu and Ni. This fact is difficult to explain, and may suggest other source of contamination than coke production, especially since these metals have been strongly correlated to each other. Nevertheless, PAHs and heavy metals in long-term contaminated soils could have a strong negative impact on microbial diversity and activity (Maliszewska-Kordybach et al. 2009, Thavamani et al. 2012, Pessacq et al. 2015). Table 1. The limit of detection and recovery values of tested parameters Parameter Cd Pb Zn Cu Ni Naphthalene Acenaphthene Fluorene Phenanthrene Anthracene Fluoranthene Pyrene Benzo(a)anthracene Chrysene Benzo(b)fluoranthene Benzo(k)fluoranthene Benzo(a)pyrene Benzo(ghi)perylene Dibenzo(ah)anthracene Indeno(1,2,3-cd)pyrene Limit of detection [mg kg-1] 0.87 3.39 8.73 4.34 4.16 0.066 0.01 0.003 0.008 0.002 0.016 0.011 0.006 0.006 0.007 0.003 0.004 0.004 0.002 0.004 Recovery [%] 84.8 95.7 98.8 95.0 95.3 >70.0 >70.0 >70.0 >70.0 >70.0 96.1 95.4 93.3 95.3 94.4 93.2 88.8 90.4 92.6 90.3 Table 2. The concentrations of individual and total PAHs in tested soilsa PAHs [mg kg-1] Naphthalene Acenaphthene Fluorene Phenanthrene Anthracene Fluoranthene Pyrene Benzo(a)anthracene* Chrysene* Benzo(b)fluoranthene* Benzo(k)fluoranthene* Benzo(a)pyrene* Benzo(ghi)perylene Dibenzo(ah)anthracene* Indeno(1,2,3-cd)pyrene* Acenaphthylene 16 PAHs SOILS I 24.6 ± 4.1 8.2 ± 1.1 4.9 ± 0.6 25.0 ± 2.5 7.7 ± 0.7 29.9 ± 3.0 18.0 ± 1.8 10.4 ± 0.8 6.8 ± 0.6 7.3 ± 0.7 3.9 ± 0.3 9.7 ± 0.8 4.4 ± 0.4 2.2 ± 0.2 8.3 ± 0.8 ­ 171.4 II 476.6 ± 78.6 54.1 ± 7.3 85.6 ± 11.1 392.8 ± 39.3 119.7 ± 11.4 322.0 ± 32.2 189.3 ± 18.9 101.5± 8.1 68.5 ± 5.5 62.9 ± 6.3 37.1 ± 3.2 88.2 ± 7.1 42.1 ± 3.8 20.5 ± 1.7 76.0 ± 6.8 ­ 2137.1 III 14.0 ± 2.3 5.3 ± 0.7 20.9 ± 2.7 92.5 ± 9.3 21.3 ± 2.0 78.0 ± 7.8 58.3 ± 5.8 36.9 ± 3.0 37.5 ± 3.0 25.2 ± 2.5 13.4 ± 1.1 34.3 ± 2.7 14.1 ± 1.3 1.2 ± 0.1 24.2 ± 2.2 4.7 ± 0.5 481.6 IV 37.8 ± 6.2 4.2 ± 0.6 21.7 ± 2.8 95.3 ± 9.5 18.1 ± 1.7 166.8 ± 16.7 63.6 ± 6.4 37.6 ± 3.0 25.7 ± 2.1 24.5 ± 2.5 12.7 ± 1.1 32.9 ± 2.6 18.3 ± 1.6 4.2 ± 0.4 31.8 ± 2.9 9.4 ± 1.0 604.6 * Carcinogenic PAHs a Data are presented as average ± standard deviation, n = 3. Table 3. The physical and chemical properties of tested soilsa Parameter Sand % Silt % Clay % pH Moisture % OM % Total N % Total P % Conductivity S cm Cd mg kg soil -1 -1 SOILS I 95 5 0 7.2 ± 0.1 20.1 ± 1.0 2.19 ± 0.12 5.3 ± 0.1 0.01 ± 0.00 367.3 ± 18.4 2.6 ± 0.4 57.8 ± 2.9 247 ± 12.4 11.8 ± 0.5 10.5 ± 0.4 II 80 2 18 8.3 ± 0.2 21.7 ± 1.1 8.79 ± 0.81 0.09 ± 0.00 0.01 ± 0.00 214.7 ± 10.7 3.1 ± 0.5 1843 ± 92.2 704 ± 35.2 99.7 ± 4.0 184.2 ± 7.4 III 98 2 0 7.5 ± 0.1 14.1 ± 0.7 1.95 ± 0.53 0.1 ± 0.00 1.24 ± 0.03 95.0 ± 4.8 3.1 ± 0.5 3378 ± 168.9 556 ± 27.8 143.8 ± 5.8 519.7 ± 20.8 IV 96 4 0 7.3 ± 0.1 13.4 ± 0.7 2.25 ± 0.3 0.17 ± 0.00 1.16 ± 0.03 248.1 ± 12.4 2.8 ± 0.5 3478 ± 173.9 504 ± 25.2 169.8 ± 6.8 666.4 ± 26.7 Pb mg kg soil -1 Zn mg kg soil -1 Cu mg kg-1 soil Ni mg kg-1 soil Data are presented as average ± standard deviation, n = 3. Enzyme activity The activity of soil enzymes is often used to monitor the effects of various pollutants on microbial functioning (Turgay et al. 2010, Tang et al. 2014, Tomkiel et al. 2015). In the tested soils, the activities of a range of enzymes contributing to the biogeochemical cycle of carbon, nitrogen and phosphorus were determined (Fig. 1A­C). In addition, the activity of dehydrogenase, present only in viable cells was assessed. DHA is linked to microbial respiration and reflects the ability of microorganisms to degrade hydrocarbons (Achuba and Okoh 2014, Wang et al. 2014). In the tested soils, despite the presence of strong PAHs and metal co-contamination, a significant activity of the investigated enzymes was observed. The most contaminated soil exhibited a particularly high DHA as well as significant acid and alkaline phosphatase activity (Fig. 1A­B). Moreover, we observed strong positive correlations (p<0.01) between the activity of these enzymes and the total PAHs (Table 4). Several studies have demonstrated that DHA is among the most sensitive parameters for evaluation of hydrocarbons toxicity and stress (Irha et al. 2003, Gianfreda et al. 2005). Turgay et al. (2010) observed that the decrease in the content of hydrocarbons resulted in an increase of DHA over time during the bioremediation process of soil contaminated with petroleum hydrocarbons. However, Tang et al. (2012) examining the bioavailability of petroleum hydrocarbons under different attenuation conditions noticed an increase in DHA immediately after the contamination of soil and its decline with the duration of the experiment. Observed in our study high dehydrogenase and other enzyme activities in the most contaminated soil suggested significant intensification of microbial processes. This phenomenon may be a microbial response to a constant stress exerted by contaminants. Moreover, such results may be related to the high level of organic matter and moisture found in the most contaminated soil II. The highest urease activity and the lowest dehydrogenase and acid phosphatase activities were found in soil IV. Soil I, which was the least contaminated one, had the lowest activities of urease and alkaline phosphatase. Contrary to some other studies (Gianfreda et al. 2005, Turgay et al. 2010, Wu et al. 2014) no correlation was observed between urease activity and the level of PAHs. Thus, the variations in enzyme activities between soil samples showed that in tested soils pollution inputs have specific site effects, which may be a result of lower sensitivity of local microbial communities to PAHs and heavy metals. For instance, Qasemian et al. (2012) found an increased activity of acid phosphatase, cellulase and -glucosidase in the leaf litter of forest pine (Pinus halepensis) after 3 months of anthracene contamination. The changes in the activity of an individual enzyme did not allow for a precise quantification of soil degradation by long-time mixed contamination. This finding demonstrates that microorganisms still function when subjected to high levels of pollution. Catabolic diversity of bacterial communities The potential activity of bacterial communities that were able to grow on the substrate provided was measured using Biolog method. Very low metabolic activity with respect to 31 substrates in communities from soils I and II indicated the microbial activity at a basal metabolic rate. It is remarkable that soil III has shown the highest value of average well-color development (AWCD) which expressed the catabolic activity of metabolic communities on EcoPlatesTM but a low enzyme activity. The low activity of tested enzymes may be related to the inhibition effects caused by toxic hydrocarbons. Moreover, soil enzyme activity provides information about the breaking of different chemical bonds whereas the Biolog method reflects the metabolic activity of soil bacteria in response to 31 carbon sources, which were present at high concentrations (Palmroth et al. 2005). Furthermore, CLPPs do not measure the functional abilities of the entire soil microbial community but only that of a very limited subset of microbial fraction capable to grow Table 4. Pearson's correlation coefficients between environmental factors, enzyme activities and functional diversity indices Total P Conductivity Cd Pb Zn Cu Ni 16 PAHs NAH DHA Urease AcdP HMW PAHs Cd pH Moisture OM Total N pH Moisture OM 0.95** Total N -0.51 0.61* -0.32 Total P -0.64* 0.19 0.89** 0.17 0.83** 0.94** -0.41 -0.28 -0.57 0.61* -0.52 -0.25 0.88** 0.90** 0.85** -0.68* 0.20 -0.72** 0.23 0.73** 0.68* -0.71** 0.24 0.73** -0.43 0.79** 0.13 0.78** 0.12 0.10 0.08 -0.14 0.63* -0.12 0.68* -0.18 0.11 0.71** 0.30 0.78** -0.50 0.16 0.63* 0.63* 0.58* -0.12 0.65* -0.15 0.70** -0.08 -0.08 0.63* -0.14 0.70** -0.05 -0.35 -0.38 0.80** -0.33 -0.11 0.73** 0.73** 0.68* -0.31 0.78** 0.19 0.88** 0.28 -0.02 -0.23 0.74** 0.05 0.82** 0.14 -0.16 0.99** 0.98** 0.96** 0.03 0.91** 0.84** -0.43 -0.46 -0.46 0.94** 0.92** 0.15 0.87** 0.83** -0.33 -0.37 0.99** -0.11 0.92** 0.81** -0.56 -0.60* -0.59* -0.20 -0.61* 0.30 0.96** 0.37 0.95** -0.72** 0.56 0.99** 0.63* -0.71** 0.96** 0.59* -0.79** 0.56 -0.83** -0.40 -0.89** -0.57 Conductivity 0.59* -0.07 0.81** Cd 0.81** -0.25 0.61* -0.85** Pb -0.85** -0.15 -0.88** Zn 0.83** -0.26 0.68* -0.90** Cu -0.83** -0.07 -0.90** Ni -0.21 -0.89** -0.36 -0.73** 16 PAHs 0.98** 0.97** -0.52 HMW PAHs 0.96** 0.93** -0.63* NAH 0.96** 0.99** DHA 0.97** 0.98** -0.33 Urease -0.12 -0.63* -0.13 AcdP 0.91** 0.89** 0.89** 0.80** -0.52 -0.55 -0.54 -0.15 -0.12 0.28 0.26 0.22 0.90** -0.55 -0.51 -0.45 AlkP 0.88** 0.76** -0.52 Rs -0.78** -0.54 -0.44 H' -0.37 -0.72** -0.43 I' -0.37 -0.64* -0.57 -0.39 Statistically significant correlations: * p < 0.05; ** p < 0.01. NAH ­ naphthalene Fig. 1. Acid and alkaline phosphatase (A), dehydrogenase (B) and urease (C) activity in tested soils on the carbon sources provided (Palmroth et al. 2005, Hueso et al. 2012). To assess the potential functional diversity of the bacterial community H', I' and Rs indices were calculated (Table 5). The functional diversity and evenness indices were the highest in soil III and the lowest in soil I and II. Microbial communities in soil III were characterised by high species richness and an even distribution of functional abilities within the community. Bacterial communities in the most and the least contaminated soils reveal very low functional biodiversity. No correlations were found between the Shannon indices and the sum of 16 PAHs, which is in contradistinction to the findings of Zhang et al. (2010) who reported a correlation between these two factors. However, Muckian et al. (2007) found that the level of PAHs did not alter the number of ribotypes in the contaminated soils. They indicated that the PAH ring structure plays a greater role in influencing the bacterial community than the total PAHs concentration. In our study, a negative correlation between H' and naphthalene, a 2-ring PAH, was observed (Table 4). This can be explained by different properties of these hydrocarbons. Naphthalene is the most soluble among PAHs and thus more mobile and toxic to microorganisms. No correlations were found between functional diversity indices and enzyme activities. The functional diversity of bacterial communities was positively correlated with total phosphorus (r=0.90), lead (r=0.73) and nickel content (r=0.73, p<0.01). There were also observed strong negative correlations between diversity indices and moisture (r=-0.72) and conductivity (r=-0.72, p<0.01). It seems that some metals and physicochemical factors are more important in controlling the soil bacterial diversity than hydrocarbons (Azarbad et al. 2013). The number of utilized carbon sources, which is defined as the richness index, varied depending on the soil, from 26 to 31 out of 31 carbon sources available on the Biolog EcoPlatesTM (Table 5). Two soils exhibited strongly reduced richness indices: the most contaminated (II) one and the least contaminated (I) one. Table 5. Richness (Rs), Shannon­Weaver (H') and evenness (I') indices calculated on carbon substrate used in Biolog EcoPlateTM Index Rs H' I' SOILS I 27 ± 1.2 a 1.30 ± 0.02 a 0.91 ± 0.02 a II 26 ± 0.6 a 1.28 ± 0.01 a 0.90 ± 0.01 a III 31 ± 0.6 b 1.44 ± 0.03 b 0.97 ± 0.01 b IV 29 ± 0.6 c 1.36 ± 0.01 c 0.94 ± 0.01 c Values with the same letter mean indicate that there are no significant differences (p < 0.05) Fig. 2. Principal Component Analysis (PCA) of (A) carbon utilization profiles of bacterial communities (B) environmental factors and enzymatic activities of tested soils. Each point represents the bar center of three replicates Principal component analysis All the obtained data were subjected to principal component analysis (Fig. 2A­B). The first PCA was carried out for carbon utilization profiles. PCA revealed differences in functional diversity between the bacterial communities present in the tested soils. Both the most contaminated (II) and the least contaminated (I) soils exhibited a noticeable separation from the two other soils in the first principal component (PC1), which explained 78.5% of the variation in the data. The reason for these findings may be the selection and adaptation of microorganisms to the pollutants, which cause a lack of species or strains capable to degrade a variety of carbon sources (Wu et al. 2014). The second PCA analysis was carried out for physicochemical data, 16PAHs concentrations and enzymatic activities (Fig 2B). The majority of the variance (94.4%) was explained by the two principal components. Both soils III and IV exhibited a significant separation from the other two soils (I and II) along the first and the second axis. The PC1 was highly negatively correlated with the sum of PAHs, the activity of dehydrogenase, both phosphatases, as well as with pH, whereas the second factor was highly positively correlated with moisture and negatively correlated with total P, Pb, Cu and Ni. The long-term pressure of the mixture of contaminants may have modified the microbial metabolism and/or the relative composition and size of the constituent populations. This is because microbial communities have the ability to adapt, physiologically and evolutionally, to the presence of contaminants in soil. After the period of "stress", there is an increase in respiration intensity and enzyme activity (Baran et al. 2004, Pessacq et al. 2015). This could have occurred only in the presence of acclimated microbial communities, which were capable of adapting and proliferating in the presence of high levels of PAHs and metals. Understanding these complex interactions is of major importance for the bioremediation of contaminated soils (Thavamani et al. 2012). Conclusions Our study has demonstrated that enzyme activities and microbial functional diversity in PAHs and metal contaminated soils is related and characteristic to a specific site. This fact is connected with the high heterogeneity of soil and uneven distribution of hydrocarbons in soil. The obtained results showed also that the evaluation of a single parameter such as the activity of enzyme does not reflect the real effect of contaminants input on the activity and diversity of microbial communities. Moreover, the results of this study clearly showed that microbial communities have enormous adaptation abilities to the mixed inorganic and organic contamination. Acknowledgements We would like to thank the professional team of Analytical Laboratory in IETU for physical and chemical analysis. This work was supported by Grant No. N305 051340 financed by the Polish Ministry of Science and Higher Education. Klimek, B. & Nikliska, M. (2007). Zinc and copper toxicity to soil bacteria and fungi from zinc polluted and unpolluted soils: A comparative study with different types of Biolog plates, Bulletin of Environmental Contamination and Toxicology, 78, pp. 112­117. Lladó, S., Jiménez, N., Viñas, M. & Solanas, A.M. (2009). Microbial populations related to PAH biodegradation in an aged biostimulated creosote-contaminated soil, Biodegradation, 20, pp. 593­601. Lu, M., Xu, K. & Chen, J. (2013). Effect of pyrene and cadmium on microbial activity and community structure in soil, Chemosphere, 91, pp. 491­497. Lu, X-Y., Zhang, T. & Fang, H.H-P. (2011). Bacteria-mediated PAH degradation in soil and sediment, Applied Microbiology and Biotechnology, 89, pp. 1357­1371. Maliszewska-Kordybach, B., Smreczak, B. & Klimkowicz-Pawlas, A. (2009). Effects of anthropopressure and soil properties on the accumulation of polycyclic aromatic hydrocarbons in the upper layer of soils in selected regions of Poland, Applied Geochemistry, 24, pp. 1918­1926. Muckian, L., Grant, R., Doyle, E. & Clipson, N. (2007). Bacterial community structure in soils contaminated by polycyclic aromatic hydrocarbons, Chemosphere, 68, pp. 1535­1541. Nam, J.J., Thomas, G.O., Jaward, F.M., Steinnes, E., Gustafsson, O. & Jones, K.C. (2008). PAHs in background soils from Western Europe: influence of atmospheric deposition and soil organic matter, Chemosphere, 70, pp. 1596­1602. Palmroth, M.R.T., Münster, U., Pichtel, J. & Puhakka, J.A. (2005). Metabolic responses of microbiota to diesel fuel addition in vegetated soil, Biodegradation, 16, pp. 91­101. Pessacq, J., Medina, R., Terada, C., Bianchini, F.E., Morelli, I.S. & Del Panno, M.T. (2015). Assessment of the responsiveness to different stresses of the microbial community from long-term hydrocarbon-contaminated soils, Water, Air, & Soil Pollution, 226, 20. DOI 10.1007/s11270-014-2262-9. PN-R-04033 (1998). Gleby i utwory mineralne. Podzial na frakcje i grupy granulometryczne. Polish Standard. Distribution of the particle sizes in the soil. Polish Committee for Standardization. (in Polish) Qasemian, L., Guiral, D., Ziarelli, F., Van Dang, T.K. & Farnet, A.M. (2012). Effects of anthracene on microbial activities and organic matter decomposition in a Pinus halepensis litter from a Mediterranean coastal area, Soil Biology and Biochemistry, 46, pp. 148­154. Rachwal, M., Magiera, T. & Wawer, M. (2015). Coke industry and steel metallurgy as the source of soil contamination by technogenic magnetic particles, heavy metals and polycyclic aromatic hydrocarbons, Chemosphere, http://dx.doi.org/10.1016/j. chemosphere.2014.11.077 Ross, S. (1994). Retention, transformation and mobility of toxic metals in soils, in: Toxic metals in soil-plant systems, Ross, S. (Ed.), John Woley Sons, Chichester UK 1994. Schinner, F., Öhlinger, R., Kandeler, E. & Margesin, R. (1996). Methods in soil biology. Springer, Heidelberg Berlin, New York 1996. Tabatabai, M.A. & Bremner, J.M. (1969). Use of p-nitrophenyl phosphate in assay of soil phosphatase activity, Soil Biology and Biochemistry, 1, pp. 301­307. Tabatabai, M.A. & Bremner, J.M. (1972). Assay of urease activity in soils, Soil Biology and Biochemistry, 4, pp. 479­487. Tang, X., Hashmi, M.Z., Long, D., Chen, L., Khan, M.I. & Shen, Ch. (2014). Influence of heavy metals and PCBs pollution on the enzyme activity and microbial community of paddy soils around an e-waste recycling workshop, International Journal of Environmental Research and Public Health, 11, pp. 3118­3131.

Journal

Archives of Environmental Protectionde Gruyter

Published: Dec 1, 2016

References