http://orcid.org/0000-0001-8517-8454
Figures
Abstract
Long-term exposure to polyaromatic hydrocarbons (PAHs) has been connected to chronic human health disorders. It is also well-known that i) PAH contamination alters soil bacterial communities, ii) human microbiome is associated with environmental microbiome, and iii) alteration in the abundance of members in several bacterial phyla is associated with adverse or beneficial human health effects. We hypothesized that soil pollution by PAHs altered soil bacterial communities that had known associations with human health. The rationale behind our study was to increase understanding and potentially facilitate reconsidering factors that lead to health disorders in areas characterized by PAH contamination. Large containers filled with either spruce forest soil, pine forest soil, peat, or glacial sand were left to incubate or contaminated with creosote. Biological degradation of PAHs was monitored using GC-MS, and the bacterial community composition was analyzed using 454 pyrosequencing. Proteobacteria had higher and Actinobacteria and Bacteroidetes had lower relative abundance in creosote contaminated soils than in non-contaminated soils. Earlier studies have demonstrated that an increase in the abundance of Proteobacteria and decreased abundance of the phyla Actinobacteria and Bacteroidetes are particularly associated with adverse health outcomes and immunological disorders. Therefore, we propose that pollution-induced shifts in natural soil bacterial community, like in PAH-polluted areas, can contribute to the prevalence of chronic diseases. We encourage studies that simultaneously address the classic “adverse toxin effect” paradigm and our novel “altered environmental microbiome” hypothesis.
Citation: Parajuli A, Grönroos M, Kauppi S, Płociniczak T, Roslund MI, Galitskaya P, et al. (2017) The abundance of health-associated bacteria is altered in PAH polluted soils—Implications for health in urban areas? PLoS ONE 12(11): e0187852. https://doi.org/10.1371/journal.pone.0187852
Editor: Marie-Joelle Virolle, Universite Paris-Sud, FRANCE
Received: August 12, 2017; Accepted: October 29, 2017; Published: November 16, 2017
Copyright: © 2017 Parajuli et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The bacterial sequence data are available in the Sequence Read Archive at NCBI under accession number SRR5229978. All other relevant data are within the paper and its Supporting Information files.
Funding: The research was funded by Finnish Funding Agency for Technology and InnovationTekes (decision number 40333/14), https://www.tekes.fi/en/. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Pollutants, such as polyaromatic hydrocarbons (PAHs), can lead to shifts in microbial communities [1–4]. PAH degradation typically lasts for decades in polluted environments including soil, water, air, and sediments [5–7]. Direct PAH toxicity, abiotic transformation, and microbial degradation favor some microbial taxa, whereas others become less prevalent in PAH-polluted soils [8]. The bacterial taxa that thrive under PAH exposure vary depending on the types of soil and environmental conditions [9–12]. A search for general trends in microbial shifts in PAH-polluted soils necessitates concurrent studies of microbial communities in different soil types.
Exposure to PAHs is associated with severe human health deficits and ecological impacts, such as process of PAH-DNA adduct formation and carcinogenesis [13], and therefore they have been classified as priority environmental pollutants by the United States Environmental Protection Agency (US-EPA) and the Environmental European Agency (EEA) [14, 15]. The environmental concentration of PAHs and the connection between direct PAH exposure and human health have been well investigated [16–23]. Almost 90% of the PAHs released into the atmosphere accumulate in surface soil layers [14] where they are primarily degraded by soil bacteria or bind to soil particles [22]. The health outcomes of PAH-exposure have been attributed to direct and indirect toxin effects on humans [23, 24].
Several studies have investigated the relationship between human health and the most common bacterial taxa found in the human microbiome. A rich or dominant Proteobacterial and diminished Bacteroidetes communities in the human gut can contribute to a risk of immune system disorders including chronic obstructive pulmonary disease and asthma [25], and such effects can be transgenerational. For example women whose infants developed IgE-associated eczema had lower diversity of Bacteroidetes in their gut during the pregnancy [26]. Alteration in the Proteobacterial abundance in various regions of human body is associated with several health disorders, but its effects are twofold. On the positive side, diverse Proteobacterial community on skin is related to a reduced risk of atopy [27, 28]. On the negative side, overgrowth of Proteobacteria is connected to asthma and chronic obstructive pulmonary disease [29]. Within Proteobacteria, Betaproteobacteria are enriched in the gut microbiome of individuals with type 2 diabetes [30]. Importantly, recent findings indicate that Proteobacterial community of individual’s skin depends on land use in their living environment and that the environment tunes response to allergens [27–31]. Soil is the major reservoir of Actinobacteria, Bacteroidetes, and Proteobacteria [32, 33], and thus the living environment likely tunes individual’s responses to allergens, sugar metabolism, and health disorders. For these reasons, it is possible that adverse health effects connected to living in contaminated environment may partly be attributable to fundamental changes in the environmental microbiome people are exposed to, in addition to the direct physiological responses caused by the toxins. Surprisingly, to the best of our knowledge, no studies have focused on the potential connection between PAH pollution and soil community changes in the abundance of health-associated bacteria.
In the present study, we sampled four different surface soils, contaminated them with PAHs and followed population shifts in bacterial community in contaminated and non-contaminated containers under controlled conditions. Finally, based on previous empirical evidence as well as our new findings and novel ideas, we propose a new “altered environmental microbiome” hypothesis that should be investigated in parallel with the traditional “direct toxin effect” paradigm as a potential explanation for the complex relationship between human health and environmental pollution.
Materials and methods
PAH source
We selected creosote as the contaminant because it is among the most widely used wood preservatives with a history of more than a century in wood impregnation industry [5]. Creosote consists of 85% PAHs which are chemical compounds consisting of at least two aromatic rings fused together [34, 35]. The remaining is less than 10% phenolic compounds and 5–10% heterocyclic aromatic compounds consisting of oxygen, sulphur, and nitrogen. Therefore, creosote is a relevant compound to test environmental changes caused by PAH-pollution since it consists mainly of PAHs and is a common problem in urban areas with a history of sawmills.
Soil collection
Surface soil was collected at four separate sites in Finland. The first site (Peat in Table 1) was at Haapasuo peat production area in Leivonmäki, Finland [36] (61°54’N 26°4’E). The current surface layer was formed hundreds of years ago, and it represents the transition layer between minerotrophic Carex and ombrotrophic Sphagnum dominated peat bogs. The original bog was drained, and the water table has been kept 30–40 cm below the peat surface for more than three decades. The second site, referred to as pine forest soil, is a boreal pine (Pinus sylvestris L.) forest in Hollola, located in southern Finland (61°0’N 25°29’E) moraine ridge (the undisturbed ecosystem has been described before [36]). The site is characterized by a thin organic soil layer on top of mineral soil [37]. The third site (Spruce forest soil in Table 1) is a spruce forest in Vierumäki, southern Finland (60°52’7N 25°41’E). Spruce (Picea abies) is the dominant tree species, and the field layer consisted of Sphagnum and Pleurozium schreberi mosses as well as dwarf shrubs, mainly Vaccinium species. The fourth site (sand) was next to (distance 50 m) the site 3 (Spruce forest soil), but it consisted of bare glacial sand and scattered ruderal herbs as original surface soil was removed years earlier when the site became a storage field for lumber.
At all four sites, soil was collected from three separate (distance >5m) ca. 1 m2 plots as described in earlier work [38, 39]. In short, at each plot, live vegetation and plant debris were removed and 15–20 L of surface soil (depth 2–15 cm) was collected and mixed thoroughly. The soil was then randomly divided into two 10 L polyethylene containers, and the procedure was repeated at each 1 m2 plot. Filled containers were covered with polyethylene lids with two 5 mm diameter holes sealed loosely with cotton wool to facilitate passive aeration. Soil weight in containers varied between 3–12 kg depending on soil type, the heaviest being the mineral soil-types.
Basic physicochemical analyses.
Moisture content was measured by drying samples in an oven (+ 90°C) for 24 hours. Organic matter (OM) was determined as loss of ignition at 550°C for 4h. Nutrient contents of the different soil types were determined with QuikChem 8000 flow injection analysis system (LACHAT Instruments Inc., USA) [38]. To measure the pH, 10g (fresh weight) of each soil type was mixed in 50 mL of 1M CaCl2, shaken for 5 min and allowed to settle for 2–24 hours [38].
PAH degradation experiment.
For each soil type, three pairs of 10L polyethene containers were included. In the beginning of the experiment, soil was mixed thoroughly within pairs in large polyethene containers and thereafter divided again into two 10L polyethene containers. Of each pair, one container was randomly chosen to be spiked with creosote, whereas another received 100g glacial sand (Lohja Rudus Oy, Lahti, Finland). Initial (day 0) samples for basic chemical analyses were acquired immediately after spiking. Creosote spiking was done as described earlier [38]. In short, creosote (6–24 g) was mixed thoroughly with 100g glacial sand, and the mixture was incorporated into soil in the containers selected for contamination. This resulted in the concentration of ~ 1% PAHs similar to those in creosote contaminated sites. After creosote addition, containers were left to incubate for four weeks (28 days) at 16 ± 1°C and sampled at days 28, 91, and 189 for chemical and bacterial community analyses. Each sample (total weight 10 g) consisted of five subsamples from four depths (2, 5, 10, and 15 cm) collected from each container. Community sequencing analyses were performed on day 0 and when the total PAH concentration had decreased by more than 20% from week four (day 28) value. Therefore, week 13 (day 91) samples were analyzed in most cases. However, week 27 (day 189) samples were utilized for each of control and contaminated mineral soil containers because their PAH concentrations at week 13 were more than 95% of week 4 concentrations. To avoid the transfer of microorganisms between containers, the sampling equipment was carefully flame sterilized with 70% ethanol right before taking each sample.
PAH analyses
PAH concentrations were determined using toluene extraction as described in Nordic Guidelines for Chemical Analysis of Contaminated Soil Samples [40] and analyzed as described earlier [41] with the exception that aqueous sodium pyrophosphate decahydrate solution (0.05 M) instead of hexane was used. The extracts were analyzed with Shimadzu GC–MS-QP5000 system equipped with AOC-20i autoinjector and 30-m ZB-5MS column (0.25 mm i.d., 0.25 μm film thickness). The oven program was set as follows: 80°C for 1 min, 10°C/min to 250°C, 7°C/min to 280°C, 20°C/min to 320°C with a hold of 10 min for a total run time of 34.29 min. PAH-mix 9 (16 PAHs included) was used for GC-MS and PAH-Mix 31 of five deuterated PAHs for soil samples (Dr. Ehrenstorfer, GmbH Germany) as standards and Anthracene-D10 (Dr. Ehrenstorfer, GmbH Germany) as a recovery standard.
DNA extraction, amplification, and sequencing
Total DNA was extracted from each soil sample using FastDNA SPIN Kit for Soil (MP Biomedicals, Illkirch, France) according to the manufacturer’s standard protocol. The highly hypervariable region V3 of bacterial 16S rRNA gene was amplified using primer constructs that incorporated the pyrosequencing adapters (A), sample-specific DNA tags, and MF341 5’ CTA CGG GAG GCA GCA G 3’ or R518 5′ ATT ACC GCG GCT GCT GG 3′ [42]. The PCR were conducted under the following conditions: 200 nM of each forward and reverse primers, 5 ng template DNA, 200 μM of each dNTP, 2.5 mM MgCl2, 1 U GoTaq Hot Start DNA polymerase (Promega, Madison, WI), and 2.5 μl PCR buffer. The PCR cycle parameters consisted of an initial denaturation at 94°C for 3 min, then 25 cycles of denaturation at 94°C for 1 min, annealing at 54°C for 1 min, and extension at 72°C for 2 min, followed by a final extension step at 72°C for 10 min. The PCR products were purified using Ampure XP Magnetic Clean-up (Agencourt Bioscience Corporation, Beverly, MA, USA), and quantified using Nanodrop (Thermo Scientific, Rockford, IL, USA) and Bioanalyzer 2100 with DNA 1000 chips (Agilent Technologies Inc., Santa Clara, CA, USA). The sequencing was performed using the 454 GS FLX protocol and the GS FLX Titanium Rapid Library Preparation Kit (454 Life Sciences, Roche Diagnostics, CT, USA). The bacterial sequence data are available in the Sequence Read Archive at NCBI under accession number SRR5229978.
Sequence analysis.
The sequence data were analyzed using MOTHUR (v1.35.0, 64-bit for Linux) according to the standard operating protocol as described earlier [43]. Briefly, the raw sequence data were quality controlled, and reads with ambiguous bases or homopolymers longer than 8 bp (722 sequences) were removed. The UCHIME algorithm [44] identified 943 sequences as chimeric, and these were subsequently omitted. Sequences were assigned into Operational Taxonomic Units (OTU) at 97% similarity and OTUs assigned to taxon affinities using Naïve Bayesian Classifier [45] against the RDP training set (version 10). Rare OTUs occurring three or fewer times across all samples were omitted to avoid problems caused by uncertainty in origin.
Statistical analyses.
The relative abundance was calculated as the number of sequences in a taxon divided by total number of sequences in a sample. To compare contaminated and pristine treatments, t-test on relative abundances of bacterial phyla was conducted. Significant results are indicated by asterisks (p<0.05). False Detection Rate (FDR) was used for the p-values correction in t-tests. The difference in the pH and the organic matter content across the four soil types was calculated using ANOVA in JMP (v.11.0 64-bit; SAS Institute, Cary, North Carolina). To visualize bacterial community compositions of whole bacterial communities as well as the major bacterial phyla, non-metric multidimensional scaling (NMDS) analyses was performed based on relative abundance of OTUs using vegan package in R (v3.2.2, R Development Core Team 2015). The Bray-Curtis distance was chosen in the NMDS analysis since there were null values between samples in the data [46]. We performed the NMDS analyses at 99, 97, 95, 93, 91% OTU similarity levels but only present the data at 97% because OTU threshold had no impact on the overall conclusions.
Differences in the bacterial community composition between creosote contaminated and control groups were tested using permutative analysis of variation (PERMANOVA, function adonis in R-package vegan and Bray-Curtis metric). The difference in the relative abundances of bacterial phyla between creosote contaminated and control soils was calculated using the T-test in JMP.
Results
Soil chemical characteristics
The four soil types used in this study differed in pH (F = 164.4, df = 3, 20, p < 0.001) and organic matter content (F = 183.2, df = 3, 20, p < 0.001). Sandy soil had the highest and spruce forest soil the lowest pH (Table 1). Organic matter content decreased in the following order: peat > pine forest soil > spruce forest soil > sandy soil (Table 1). Thus, the four soil types represent a range of habitats with highly variable bacterial communities.
Characteristics of bacterial operational units (OTUs) and relative abundances of bacterial phyla
The pyrosequenced bacterial 16S rRNA gene fragment datasets were analysed from 24 soil samples taken on days 0, 91, and 189. Day 0 samples for microbiological analyses were taken before contamination. Then the time point at which 20% or more of PAH concentration on day 28 had disappeared was determined. In sand, 20% degradation was reached on day 189 and in other soils on day 91 (S1 Table). Day 28 was used as a reference day because the reduction in PAH concentration until day 28 is typically caused by the evaporation of naphthalene and other low molecular weight compounds [47, 48].
In total 3626 OTUs were obtained, which represented more than 20 bacterial phyla from different soil types with and without PAH contamination. Bacterial OTUs representing Proteobacteria were the most abundant group accounting for 50.9%, 31.5%, and 35.8% of the total sequences in contaminated, pristine, and day 0 samples, respectively (Table 2). Likewise, Proteobacteria was the predominant phylum when compared across the different soil types with 46.20% (peat soil), 37.7% (pine forest soil), 36.0% (sand), and 37.8% (spruce forest soil) of all sequences. Other dominant phyla in all soils were Acidobacteria, Actinobacteria, Verrucomicrobia, and Bacteroidetes (Tables 2 and 3).
The relative abundances of different phyla between contaminated and pristine treatments were compared. Proteobacteria (df = 1, p = 0.046) were more abundant in contaminated soils, whereas the relative abundances of Actinobacteria (df = 1, p = 0.043) and Bacteroidetes (df = 1, p = 0.005) declined in contaminated soils (Table 2). In addition, it was observed that the relative abundances of Parcubacteria and Candidus Saccharibacteria were lower in the contaminated soils, although the relative abundances of these two were very low (Table 2). As taxa belonging to Proteobacteria represented more than 35% of microbial community in all soil types, we compared the relative abundances of the classes under the phylum Proteobacteria and found that the relative abundance of Betaproteobacteria was higher (df = 1, p = 0.013) in contaminated soils.
Bacterial community composition
The total bacterial community composition did not differ between creosote contaminated and pristine soils at the OTU level as revealed by the nonmetric multidimensional scaling (NMDS) ordination (Fig 1a) and permutational multivariate analysis of variance PERMANOVA (Table 3). The NMDS at the genus level, however, revealed distinctly different bacterial communities in the creosote contaminated and pristine soils (Fig 1b). PERMANOVA indicated that contamination and soil type had distinct main effects on microbial community at the genus level and without a strong interaction effect (Table 3). These findings underline that creosote contamination changed microbial community at the genus level regardless of the soil type.
(a) The bacterial community composition in the two treatment groups overlap at the OTU level. (b) The communities are distinctively different at the genus level. For peat (Peat in the figure), pine forest soil (Pine.F) and spruce forest soil (Spruce.F), samples taken on week 31 (day91) were utilized, whereas week 27 (day 189) samples were used in the case of mineral soil (Sand).
Community composition analysis of major bacterial phyla at OTU level
The major bacterial phyla detected from the soil samples were examined for their differences in community composition at the OTU level between the creosote contaminated and pristine soil samples. Since OTUs representing Bacteroidetes were not detected in two of the creosote contaminated samples, community composition analysis was not performed for Bacteroidetes. Consequently, the analyses were performed only for Proteobacteria and Actinobacteria.
Proteobacteria.
Proteobacterial communities were distinctly different in contaminated than in control soils in NMDS (Fig 2). PERMANOVA revealed that the soil type and creosote contamination had strong individual effects on the Proteobacterial community composition, while the significance of the interaction term was ten times lower (Table 4).
The community composition is noticeably different in the two treatment groups. The ordination is based on Bray-Curtis dissimilarity metric.
Actinobacteria.
Although the relative abundance of Actinobacteria decreased in creosote contaminated soil samples, it was observed that its community composition did not differ between the creosote contaminated and pristine soil samples (Fig 3). Instead, the soil type affected the Actinobacterial community composition found in that soil type (Table 5).
The community composition across the two treatment groups overlap partially. The ordination is based on Bray-Curtis dissimilarity metric.
Temporal variation.
The total bacterial, Proteobacterial and Actinobacterial communities sampled on day 0 were compared with those sampled later. Day 0 samples were compared with both creosote contaminated and pristine samples to determine the temporal community dynamics. Only Proteobacterial community on day 0 distinctly differed from the communities sampled at the later time points (Fig 4a). Neither total bacterial community nor Actinobacterial community changed during the experiment significantly (Fig 4b and 4c).
(a) Proteobacterial community on day 0 is distinctively different from those sampled at later time points (day 91 for peat, pine forest, and spruce forest soils, and day 189 for sand). Neither Actinobacterial community (b) nor the total bacterial community at genus level (c) varied with time.
Discussion
In the present study, we compared the responses of bacterial community to creosote contamination in four different types of soil. Our data indicates shifts in bacterial community composition and alterations in the relative abundances of three major bacterial phyla as a response to creosote contamination. The pollution-induced shift was attributable to an increase in the relative abundance of Proteobacteria, particularly Betaproteobacteria, and a decrease in the relative abundance of Actinobacteria and Bacteroidetes. These phyla are the major bacterial phyla present in soil as well as on human skin and other biological samples [49–51]. Because they have been recently linked to human immune function and chronic diseases [25–27, 52], a major change in their community composition in the environment can potentially affect human exposure to them and, eventually, lead to modulated human immune system functions [27, 31]. As far as we know, this idea that we call the altered environmental microbiome hypothesis has never been proposed nor tested. The first step to evaluate the hypothesis is therefore to perform a literature review on the potential connections between health-associated and pollution-induced changes.
While several studies have investigated and pointed out the alteration of the entire bacterial community as a result of PAH pollution, only a handful of studies have focused on pollution-induced changes in Proteobacterial, Bacteroidetes, and Actinobacterial communities in surface soil [1, 12, 53–58]. Those studies have revealed that Proteobacteria were more abundant in PAH polluted soils compared to pristine soils. Phylum Bacteroidetes, in contrast, either decreased in abundance [55] or was not studied. The relative abundance of Actinobacteria was usually lower in contaminated soil than in non-contaminated soil [55], but its abundance is also dependent on the PAHs involved[57] (Table 6). The increase in the abundance of Proteobacteria and the decrease in Actinobacteria and Bacteroidetes are consistent with our findings. As the studies listed in Table 6 cover a wide range of geographic areas and land use histories, it seems plausible to assume that shifts in the abundance of Proteobacteria, Actinobacteria, and Bacteroidetes are a typical consequence of PAH pollution in soil.
The most common method to compare soil remediation in different treatments is to take samples at selected time intervals and compare the measured variables, such as microbial communities, between sampling dates [38, 59–61]. As the rate of soil remediation is dependent on soil type, it is conceivable that only the time-dependent changes are reported or noticed, and the microbial community changes connected to actual degradation are ignored. Since chemical analyses were done at certain intervals and the results were then compared to 20% degradation level, the abovementioned concern was overcome in this study. Therefore, changes in the relative abundance of Bacteroidetes, Actinobacteria, Proteobacteria, and particularly Betaproteobacteria were detected.
The observed higher relative abundance of the subphylum Betaproteobacteria as a result of 20% PAHs degradation may turn out to be an important observation for several reasons. Earlier findings suggest that Betaproteobacteria are indicators of the endpoint of fast bioremediation in soil [1]; in this study, soil remediation ceased before the last sampling in organic soils (data not shown). Interestingly, an increase in Betaproteobacterial community has been associated with degradation of PAHs [1] [55]. Betaproteobacteria are also enriched in the faeces of diabetic people compared to non-diabetic individuals [62] meaning that Betaproteobacteria have been associated with both soil pollution and health. Based on the current and earlier findings [1], it is worthwhile to study how common and permanent changes in the abundance of Betaproteobacteria are in contaminated environments. If the changes persist, health-perspective should be further investigated, e.g., whether the same species are enriched in polluted soil and diabetics.[55] At the level of the whole phylum Proteobacteria, a similar but less sharp increase was reminiscent of the increase in the abundance of Proteobacteria in the household dust and biological samples that have been associated with several diseases [62–67]. Hence, it should be considered whether pollution-induced changes in Proteobacterial community composition in areas prone to anthropogenic PAHs could be linked to chronic diseases [62–67].
In parallel with Proteobacteria, Bacteroidetes is another dominant phylum in the environment as well as in the human gut [68–70]. Several members of Bacteroidetes have been reported to be sensitive to hydrocarbons and are thought to be outcompeted by other hydrocarbon-degrading bacteria leading to their reduced abundance [71], which could have led to their decreased relative abundance in creosote contaminated soils in this study. Interestingly, in human body, Bacteroidetes are involved in the maintenance of immune system by promoting various regulatory components [72, 73], and therefore their reduced abundance has been implicated to have a role in immune disorders and other inflammatory diseases [26, 67, 74–76].
We also observed that there was a relative decrease in Actinobacterial abundance in creosote contaminated compared to pristine soils. This finding is in contrast with studies that pointed out that Actinobacteria contains efficient degraders of PAH and their abundance is increased in PAH polluted soils [54, 77, 78]. However, two other studies have reported a decreased abundance of Actinobacteria in creosote contaminated soil [53, 79], which is in accordance with findings of this study. A plausible explanation for the varying results is that the PAH-degrading taxa within Actinobacteria were not ubiquitous in this study. Nevertheless, the diminished abundance of Actinobacteria has been observed to have clinical significance in non-communicable diseases [62, 80]. Therefore, a pollution-induced reduction in the supply of Bacteroidetes and Actinobacteria from the environment may be a part of a novel mechanism how environmental PAH pollution is associated with health, and more specifically, with immune disorders. As PAH compounds are a major ingredient of urban smog, the adverse effects of altered soil microbial communities may be common in densely populated areas.
Recently, it has become increasingly more evident that environmental biodiversity contributes in shaping the human microbiota, and that differences in bacterial abundance and diversity are associated with health outcomes [27, 81–86]. Consequently, shifts in environmental microbial community in polluted areas may lead to similar effects in human microbial community and subsequent health consequences. We are well aware of the fact that not all members belonging to the phyla we studied have associations with human health and therefore the shifts in abundance and alteration in the community composition have to be established at a finer taxonomic level to draw definitive conclusions about the association between pollution and shifts in health associated environmental microbiome. In this study, however, we did not search for shifts in abundance and community composition beyond the phylum level (except for Proteobacteria) to avoid erroneous conclusions drawn from the shallow sequencing depth of 454-pyrosequencing technology. Therefore, future studies must utilize more powerful next generation sequencing methods to explore the relationship between pollution induced shifts in the alteration of bacterial communities and human health.
We hypothesized that PAH pollution will cause shift in the abundance and community composition of environmental bacterial phyla that have known associations with human health. Our findings that the relative abundance of Proteobacteria is increased and that of Actinobacteria and Bacteroidetes is decreased after PAH exposure supports our altered environmental microbiome hypothesis. The next step is to study the differences in skin and gut microbiome of people living in severely PAH-exposed versus mildly PAH-polluted areas. If differences in Proteobacterial, Actinobacterial, and Bacteroidetes community on human skin and in stool samples are associated with PAHs or other environmental pollutants, this would further strengthen our altered microbiome hypothesis, and the changes could well turn out to be important factors in explaining the increase in the frequency of immune mediated and other non-communicable diseases among city dwellers.
In parallel to studying the direct and indirect impacts of PAH exposure, it is crucial to know whether shifts in health-associated microbial community are restricted to PAHs and if the changes are irreversible. As this has not been the focus of previous research, a comprehensive answer requires novel research. However, one of our previous studies made with diesel oil indicates that changes in Proteobacterial community are persistent after soil remediation; Gammaproteobacteria were less frequent in naturally remediated diesel contaminated soil as compared to similar pristine soil samples [38]. As the relative abundance of Proteobacteria in general and Gammaproteobacteria in particular has been associated with the quality of living environment in healthy but not atopic individuals [27, 52], the finding of our previous study [38] is in accordance with our novel altered environmental microbiome hypothesis.
Our findings indicate that PAH pollution alters the community composition and changes the relative abundance of soil Proteobacterial, Actinobacterial, and Bacteroidetes communities. As environment can tune bacterial community on human skin [27], alterations in environmental microbiome can be an important component explaining why PAH pollution is connected to several health deficits. It may be worthwhile to consider if other pollutants change health-associated environmental microbiome similarly and also if exposure to PAHs alters human microbial community on skin and other tissues. For these reasons, the altered environmental microbiome hypothesis proposed in this study can supplement the currently prevailing direct toxin effect hypothesis and encourage further studies to explore a potential new dimension in interrelations between environmental pollution, microbiota, and human health.
Supporting information
S1 Table. Sum of total PAHs concentration (μg/kg) before and after natural attenuation.
PAHs were measured after 28, 91 and 189 days.
https://doi.org/10.1371/journal.pone.0187852.s001
(DOCX)
Acknowledgments
Heidi Jääskeläinen provided chemical analytical assistance. The research was funded by Tekes (decision number 40333/14). We thank the reviewers for their constructive feedbacks on the manuscript.
References
- 1. Lors C, Ryngaert A, Périé F, Diels L, Damidot D. Evolution of bacterial community during bioremediation of PAHs in a coal tar contaminated soil. Chemosphere. 2010;81(10):1263–71. pmid:20943246
- 2. Płociniczak T, Sinkkonen A, Romantschuk M, Piotrowska-Seget Z. Characterization of Enterobacter intermedius MH8b and its use for the enhancement of heavy metals uptake by Sinapis alba L. Applied soil ecology. 2013;63:1–7.
- 3. Wang Y, Li P, Jiang Z, Sinkkonen A, Wang S, Tu J, et al. Microbial community of high arsenic groundwater in agricultural irrigation area of Hetao Frontiers in microbiology. 7. 2016.
- 4. Sinkkonen A, Kauppi S, Pukkila V, Nan H, Płociniczak T, Kontro M, et al. Previous exposure advances the degradation of an anthropogenic s-triazine regardless of soil origin. Journal of soils and sediments. 2013;13(8):1430–8.
- 5. Mueller JG, Chapman PJ, Pritchard PH. Creosote-contaminated sites. Their potential for bioremediation. Environmental science & technology. 1989;23(10):1197–201.
- 6. Rivas FJ. Polycyclic aromatic hydrocarbons sorbed on soils: a short review of chemical oxidation based treatments. Journal of Hazardous Materials. 2006;138(2):234–51. pmid:16973263
- 7. Ferrarese E, Andreottola G, Oprea IA. Remediation of PAH-contaminated sediments by chemical oxidation. Journal of Hazardous Materials. 2008;152(1):128–39. pmid:17689010
- 8. Cerniglia CE. Biodegradation of polycyclic aromatic hydrocarbons. Current opinion in biotechnology. 1993;4(3):331–8.
- 9. Juhasz AL, Naidu R. Bioremediation of high molecular weight polycyclic aromatic hydrocarbons: a review of the microbial degradation of benzo [a] pyrene. International biodeterioration & biodegradation. 2000;45(1 SRC—GoogleScholar):57–88.
- 10. Kasai Y, Kishira H, Harayama S. Bacteria belonging to the genus Cycloclasticus play a primary role in the degradation of aromatic hydrocarbons released in a marine environment. Applied and environmental microbiology. 2002;68(11 SRC—GoogleScholar):5625–33.
- 11. Lotfabad S, Gray M. Kinetics of biodegradation of mixtures of polycyclic aromatic hydrocarbons. Applied microbiology and biotechnology. 2002;60(3):361–6. pmid:12436320
- 12. Singleton DR, Sangaiah R, Gold A, Ball LM, Aitken MD. Identification and quantification of uncultivated Proteobacteria associated with pyrene degradation in a bioreactor treating PAH‐contaminated soil. Environmental Microbiology. 2006;8(10):1736–45. pmid:16958754
- 13. Kaushik C, Sangwan P, Haritash A. Association of polycyclic aromatic hydrocarbons (PAHS) with different sizes of atmospheric particulate in hisar city and its health aspects. Polycyclic Aromatic Compounds. 2012;32(5):626–42.
- 14. Wild SR, Jones KC. Polynuclear aromatic hydrocarbons in the United Kingdom environment: a preliminary source inventory and budget. Environmental pollution. 1995;88(1 SRC—GoogleScholar):91–108.
- 15. Sun J-H, Wang G-L, Chai Y, Zhang G, Li J, Feng J. Distribution of polycyclic aromatic hydrocarbons (PAHs) in Henan reach of the Yellow River, Middle China. Ecotoxicology and Environmental Safety. 2009;72(5):1614–24. pmid:18597847
- 16. Boström E, Gerde P, Hanberg A, Jernström B, Johansson C, Kyrklund T, et al. Cancer risk assessment, indicators, and guidelines for polycyclic aromatic hydrocarbons in the ambient air. Environmental health perspectives. 2002;110(Suppl 3):451.
- 17. Unwin J, Cocker J, Scobbie E, Chambers H. An assessment of occupational exposure to polycyclic aromatic hydrocarbons in the UK. Annals of Occupational Hygiene. 2006;50(4 SRC—GoogleScholar):395–403.
- 18. Boffetta P, Jourenkova N, Gustavsson P. Cancer risk from occupational and environmental exposure to polycyclic aromatic hydrocarbons. Cancer Causes Control. 1997;8(3 SRC—GoogleScholar):444–72.
- 19. Diggs DL, Huderson AC, Harris KL, Myers JN, Banks LD, Rekhadevi PV, et al. Polycyclic aromatic hydrocarbons and digestive tract cancers: a perspective. Journal of Environmental Science and Health. 2011;29(4 SRC—GoogleScholar):324–57.
- 20. Olsson AC, Fevotte J, Fletcher T, Cassidy A, Mannetje At, Zaridze D, et al. Occupational exposure to polycyclic aromatic hydrocarbons and lung cancer risk: a multicenter study in Europe. Occupational and environmental medicine. 2010;67(2):98–103. pmid:19773276
- 21. Masih A, Taneja A. Polycyclic aromatic hydrocarbons (PAHs) concentrations and related carcinogenic potencies in soil at a semi-arid region of India. Chemosphere. 2006;65(3):449–56. pmid:16527332
- 22. Cachada A, Pato P, Rocha-Santos T, da Silva EF, Duarte A. Levels, sources and potential human health risks of organic pollutants in urban soils. Science of the Total Environment. 2012;430:184–92. pmid:22652008
- 23. Kim K-H, Jahan SA, Kabir E, Brown RJ. A review of airborne polycyclic aromatic hydrocarbons (PAHs) and their human health effects. Environment international. 2013;60:71–80. pmid:24013021
- 24. Abdel-Shafy HI, Mansour MS. A review on polycyclic aromatic hydrocarbons: source, environmental impact, effect on human health and remediation. Egyptian Journal of Petroleum. 2016;25(1):107–23.
- 25. Park H, Shin JW, Park S-G, Kim W. Microbial communities in the upper respiratory tract of patients with asthma and chronic obstructive pulmonary disease. PloS one. 2014;9(10):e109710. pmid:25329665
- 26. West CE, Renz H, Jenmalm MC, Kozyrskyj AL, Allen KJ, Vuillermin P, et al. The gut microbiota and inflammatory noncommunicable diseases: associations and potentials for gut microbiota therapies. Journal of Allergy and Clinical Immunology. 2015;135(1):3–13. pmid:25567038
- 27. Hanski I, von Hertzen L, Fyhrquist N, Koskinen K, Torppa K, Laatikainen T, et al. Environmental biodiversity, human microbiota, and allergy are interrelated. Proceedings of the National Academy of Sciences. 2012;109(21):8334–9.
- 28. Ruokolainen L, Hertzen L, Fyhrquist N, Laatikainen T, Lehtomäki J, Auvinen P, et al. Green areas around homes reduce atopic sensitization in children. Allergy. 2015;70(2):195–202. pmid:25388016
- 29. Marsland BJ, Gollwitzer ES. Host-microorganism interactions in lung diseases. Nature reviews Immunology. 2014;14(12):827. pmid:25421702
- 30. Zhang H, DiBaise JK, Zuccolo A, Kudrna D, Braidotti M, Yu Y. Human gut microbiota in obesity and after gastric bypass. Proceedings of the National Academy of Sciences. 2009;106(7 SRC—GoogleScholar):2365–70.
- 31. Haahtela T, Laatikainen T, Alenius H, Auvinen P, Fyhrquist N, Hanski I, et al. Hunt for the origin of allergy–comparing the Finnish and Russian Karelia. Clinical & Experimental Allergy. 2015;45(5):891–901.
- 32. Hui N, Liu X-X, Kurola J, Mikola J, Romantschuk M. Lead (Pb) contamination alters richness and diversity of the fungal, but not the bacterial community in pine forest soil. Boreal environment research. 2012.
- 33. Mehta C, Yu D, Srivastava R, Sinkkonen A, Kurola J, Gupta V, et al. Microbial Diversity and Bioactive Substances in Disease Suppressive Composts from India. Compost Science & Utilization. 2016;24(2):105–16.
- 34. Kulik N, Goi A, Trapido M, Tuhkanen T. Degradation of polycyclic aromatic hydrocarbons by combined chemical pre-oxidation and bioremediation in creosote contaminated soil. Journal of Environmental Management. 2006;78(4):382–91. pmid:16154683
- 35. Valderrama C, Gamisans X, De las Heras X, Farran A, Cortina J. Sorption kinetics of polycyclic aromatic hydrocarbons removal using granular activated carbon: intraparticle diffusion coefficients. Journal of Hazardous Materials. 2008;157(2):386–96.
- 36. Rantalainen M-L, Torkkeli M, Strömmer R, Setälä H. Lead contamination of an old shooting range affecting the local ecosystem—a case study with a holistic approach. Science of the total environment. 2006;369(1):99–108.
- 37. Sinkkonen A, Kauppi S, Simpanen S, Rantalainen A-L, Strömmer R, Romantschuk M. Layer of organic pine forest soil on top of chlorophenol-contaminated mineral soil enhances contaminant degradation. Environmental Science and Pollution Research. 2013;20(3):1737–45. pmid:22752813
- 38. Kauppi S, Romantschuk M, Strömmer R, Sinkkonen A. Natural attenuation is enhanced in previously contaminated and coniferous forest soils. Environmental Science and Pollution Research. 2012;19(1):53–63. pmid:21660637
- 39. Sinkkonen A, Ollila S, Romantschuk M. Changes in TcpA gene frequency explain 2, 4, 6-trichlorophenol degradation in mesocosms. Journal of Environmental Science and Health, Part B. 2014;49(10):756–9.
- 40.
Karstensen KH. Nordic guidelines for chemical analysis of contaminated soil samples: Nordtest; 1996.
- 41. Honkonen O, Rantalainen A-L. Impact of urbanization on the concentrations and distribution of organic contaminants in boreal lake sediments. Environmental monitoring and assessment. 2013;185(2):1437–49. pmid:22527470
- 42. Muyzer G, De Waal EC, Uitterlinden AG. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Applied and environmental microbiology. 1993;59(3):695–700. pmid:7683183
- 43. Schloss PD, Gevers D, Westcott SL. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PloS one. 2011;6(12):e27310. pmid:22194782
- 44. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27(16):2194–200. pmid:21700674
- 45. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and environmental microbiology. 2007;73(16):5261–7. pmid:17586664
- 46.
Borcard D, Gillet F, Legendre P. Numerical ecology with R: Springer Science & Business Media; 2011.
- 47. Cousins IT, Beck AJ, Jones KC. A review of the processes involved in the exchange of semi-volatile organic compounds (SVOC) across the air–soil interface. Science of the Total Environment. 1999;228(1):5–24.
- 48. Jacobson M, Hansson HC, Noone K, Charlson R. Organic atmospheric aerosols: Review and state of the science. Reviews of Geophysics. 2000;38(2):267–94.
- 49. Cho I, Blaser MJ. The human microbiome: at the interface of health and disease. Nature reviews Genetics. 2012;13(4):260. pmid:22411464
- 50. Zoetendal E, Rajilić-Stojanović M, De Vos W. High-throughput diversity and functionality analysis of the gastrointestinal tract microbiota. Gut. 2008;57(11):1605–15. pmid:18941009
- 51. Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science. 2008;320(5883):1647–51. pmid:18497261
- 52. Hanski I. Biodiversity, microbes and human well-being. Ethics in Science and Environmental Politics. 2014;14(1):19–25.
- 53. Mukherjee S, Juottonen H, Siivonen P, Quesada CL, Tuomi P, Pulkkinen P, et al. Spatial patterns of microbial diversity and activity in an aged creosote-contaminated site. The ISME journal. 2014;8(10):2131. pmid:25105905
- 54. Peng J, Zhang Y, Su J, Qiu Q, Jia Z, Zhu Y-G. Bacterial communities predominant in the degradation of 13 C 4–4, 5, 9, 10-pyrene during composting. Bioresource technology. 2013;143:608–14. pmid:23845709
- 55. Martin F, Torelli S, Le Paslier D, Barbance A, Martin-Laurent F, Bru D, et al. Betaproteobacteria dominance and diversity shifts in the bacterial community of a PAH-contaminated soil exposed to phenanthrene. Environmental Pollution. 2012;162:345–53. pmid:22243884
- 56. Sawulski P, Clipson N, Doyle E. Effects of polycyclic aromatic hydrocarbons on microbial community structure and PAH ring hydroxylating dioxygenase gene abundance in soil. Biodegradation. 2014;25(6):835–47. pmid:25095739
- 57. Lladó S, Covino S, Solanas A, Petruccioli M, D’annibale A, Viñas M. Pyrosequencing reveals the effect of mobilizing agents and lignocellulosic substrate amendment on microbial community composition in a real industrial PAH-polluted soil. Journal of hazardous materials. 2015;283:35–43. pmid:25261758
- 58. Tejeda-Agredano M, Gallego S, Vila J, Grifoll M, Ortega-Calvo J, Cantos M. Influence of the sunflower rhizosphere on the biodegradation of PAHs in soil. Soil Biology and Biochemistry. 2013;57:830–40.
- 59. Hansi M, Weidenhamer JD, Sinkkonen A. Plant growth responses to inorganic environmental contaminants are density-dependent: Experiments with copper sulfate, barley and lettuce. Environmental Pollution. 2014;184:443–8. pmid:24121419
- 60. Belz RG, Sinkkonen A. Selective toxin effects on faster and slower growing individuals in the formation of hormesis at the population level—A case study with Lactuca sativa and PCIB. Science of the Total Environment. 2016;566:1205–14. pmid:27267716
- 61. Płociniczak T, Sinkkonen A, Romantschuk M, Sułowicz S, Piotrowska-Seget Z. Rhizospheric bacterial strain Brevibacterium casei MH8a colonizes plant tissues and enhances Cd, Zn, Cu phytoextraction by white mustard. Frontiers in plant science. 2016;7.
- 62. Larsen N, Vogensen FK, van den Berg FW, Nielsen DS, Andreasen AS, Pedersen BK, et al. Gut microbiota in human adults with type 2 diabetes differs from non-diabetic adults. PloS one. 2010;5(2):e9085. pmid:20140211
- 63. Ciaccio CE, Barnes C, Kennedy K, Chan M, Portnoy J, Rosenwasser L. Home dust microbiota is disordered in homes of low-income asthmatic children. Journal of Asthma. 2015;52(9):873–80. pmid:26512904
- 64. Lepage P, Häsler R, Spehlmann ME, Rehman A, Zvirbliene A, Begun A, et al. Twin study indicates loss of interaction between microbiota and mucosa of patients with ulcerative colitis. Gastroenterology. 2011;141(1):227–36. pmid:21621540
- 65. Mondot S, Kang S, Furet J-P, Aguirre de Carcer D, McSweeney C, Morrison M, et al. Highlighting new phylogenetic specificities of Crohn's disease microbiota. Inflammatory bowel diseases. 2011;17(1):185–92. pmid:20722058
- 66. Chen Y, Yang F, Lu H, Wang B, Chen Y, Lei D, et al. Characterization of fecal microbial communities in patients with liver cirrhosis. Hepatology. 2011;54(2):562–72. pmid:21574172
- 67. Panzer AR, Lynch SV. Influence and effect of the human microbiome in allergy and asthma. Current opinion in rheumatology. 2015;27(4):373–80. pmid:26002029
- 68. Thomas F, Hehemann J-H, Rebuffet E, Czjzek M, Michel G. Environmental and gut bacteroidetes: the food connection. Frontiers in microbiology. 2011;2.
- 69. Yu D, Sinkkonen A, Hui N, Kurola J, Kukkonen S, Parikka P, et al. Molecular profile of microbiota of Finnish commercial compost suppressive against Pythium disease on cucumber plants. Applied Soil Ecology. 2015;92:47–53.
- 70. Zhang H, DiBaise JK, Zuccolo A, Kudrna D, Braidotti M, Yu Y, et al. Human gut microbiota in obesity and after gastric bypass. Proceedings of the National Academy of Sciences. 2009;106(7):2365–70.
- 71.
Timmis KN, McGenity T, Van Der Meer JR, de Lorenzo V. Handbook of hydrocarbon and lipid microbiology: Springer Berlin; 2010.
- 72. Mazmanian SK, Round JL, Kasper DL. A microbial symbiosis factor prevents intestinal inflammatory disease. Nature. 2008;453(7195):620. pmid:18509436
- 73. Round JL, Mazmanian SK. Inducible Foxp3+ regulatory T-cell development by a commensal bacterium of the intestinal microbiota. Proceedings of the National Academy of Sciences. 2010;107(27):12204–9.
- 74. Kostic AD, Gevers D, Pedamallu CS, Michaud M, Duke F, Earl AM, et al. Genomic analysis identifies association of Fusobacterium with colorectal carcinoma. Genome research. 2012;22(2):292–8. pmid:22009990
- 75. Jensen A, Fagö-Olsen H, Sørensen CH, Kilian M. Molecular mapping to species level of the tonsillar crypt microbiota associated with health and recurrent tonsillitis. PLoS One. 2013;8(2):e56418. pmid:23437130
- 76. de Steenhuijsen Piters WA, Huijskens EG, Wyllie AL, Biesbroek G, Van Den Bergh MR, Veenhoven RH, et al. Dysbiosis of upper respiratory tract microbiota in elderly pneumonia patients. The ISME journal. 2016;10(1):97–108. pmid:26151645
- 77. Antizar-Ladislao B, Spanova K, Beck AJ, Russell NJ. Microbial community structure changes during bioremediation of PAHs in an aged coal-tar contaminated soil by in-vessel composting. International Biodeterioration Biodegradation. 2008;61(4 SRC—GoogleScholar):357–64.
- 78. Johnsen AR, Schmidt S, Hybholt TK, Henriksen S, Jacobsen CS, Andersen O. Strong impact on the polycyclic aromatic hydrocarbon (PAH)-degrading community of a PAH-polluted soil but marginal effect on PAH degradation when priming with bioremediated soil dominated by mycobacteria. Applied and environmental microbiology. 2007;73(5 SRC—GoogleScholar):1474–80.
- 79. Törneman N, Yang X, Bååth E, Bengtsson G. Spatial covariation of microbial community composition and polycyclic aromatic hydrocarbon concentration in a creosote‐polluted soil. Environmental toxicology and chemistry. 2008;27(5):1039–46. pmid:18419193
- 80. Marri PR, Stern DA, Wright AL, Billheimer D, Martinez FD. Asthma-associated differences in microbial composition of induced sputum. Journal of Allergy and Clinical Immunology. 2013;131(2):346–52. e3. pmid:23265859
- 81. Schnorr SL, Candela M, Rampelli S, Centanni M, Consolandi C, Basaglia G, et al. Gut microbiome of the Hadza hunter-gatherers. Nature communications. 2014;5:3654. pmid:24736369
- 82. Garrels W, Talluri TR, Apfelbaum R, Carratalá YP, Bosch P, Pötzsch K, et al. One-step multiplex transgenesis via sleeping beauty transposition in cattle. Scientific reports. 2016;6:21953. pmid:26905416
- 83. Mah KW, Sangsupawanich P, Tunyapanit W, van Bever H, Shek LP, Chua KY, et al. Gut microbiota of children living in rural south Thailand and urban Singapore. Allergology International. 2008;57(1):65–71. pmid:18089942
- 84. De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proceedings of the National Academy of Sciences. 2010;107(33):14691–6.
- 85. Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiome viewed across age and geography. nature. 2012;486(7402):222. pmid:22699611
- 86. Ying S, Zeng D-N, Chi L, Tan Y, Galzote C, Cardona C, et al. The influence of age and gender on skin-associated microbial communities in urban and rural human populations. PLoS One. 2015;10(10):e0141842. pmid:26510185