http://orcid.org/0000-0002-1413-294X
Figures
Abstract
Increased use of salting to de-ice roadways, especially in urban areas, is leading to elevated salinity levels in soil as well as surface- and ground water. This salt pollution may cause long-term ecological changes to soil and aquatic microbial communities. In this study, we examined the impact on microbial communities in soils exposed to urban road salt runoff using both culturing and 16S amplicon sequencing. Both methods showed an increase in halophilic Bacteria and Archaea in samples from road salt-exposed areas and suggested that halophiles are becoming persistent members of microbial communities in urban, road salt-impacted soils. Since salt is a pollutant that can accumulate in soils over time, it is critical to begin assessing its impact on the environment immediately. Toward this goal, we have developed a facile semi-quantitative assay utilizing halophilic microbes as biomarkers to evaluate on-going salt pollution of soils.
Citation: Pecher WT, Al Madadha ME, DasSarma P, Ekulona F, Schott EJ, Crowe K, et al. (2019) Effects of road salt on microbial communities: Halophiles as biomarkers of road salt pollution. PLoS ONE 14(9): e0221355. https://doi.org/10.1371/journal.pone.0221355
Editor: Christopher Staley, University of Minnesota Twin Cities, UNITED STATES
Received: February 23, 2019; Accepted: August 5, 2019; Published: September 4, 2019
Copyright: © 2019 Pecher 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: Raw sequence data have been deposited with links to BioProject accession number PRJNA522438 in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/) and 16S partial sequences of the isolates in GenBank (https://www.ncbi.nlm.nih.gov/genbank/) with the accession numbers MK543962 to MK543988, MK543991 to MK543994, and MK543996 to MK544007.
Funding: The research was supported by private donations through the Haloarchaeal Education & Research Development (HERD) fund (https://medschool-umaryland.givecorps.com/projects/3511-school-of-medicine-research-herd-fund) to SD. 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
Salt pollution is a long-term environmental concern, potentially threatening soil, lake, and stream ecosystems, and groundwater supplies, as well as coastal regions [1–7]. In particular, road salt crystals and brines used to de-ice roads have become major sources of salt pollution in cold regions of the world, while their use has increased substantially over the past decades [2, 8, 9] (Fig 1). Even as the frequency of snow storms and total accumulation have diminished due to climate change, larger and more intense winter precipitation events are leading to increased demand for application of road salt. Global climatic trends are predicted to further increase these extreme events, including record-breaking storms [10–12]. This trend will likely continue, especially in north-central USA, where climate change models predict an increase in the intensity of snowfall and ice storm events [12].
Average snow accumulation in Central Park, New York, and salt usage for de-icing for the USA are shown using a 20-year sliding window starting in the winter season of 1940/41. The sliding window was chosen to minimize annual fluctuations which would otherwise mask the long-term trend. Snow accumulation data were obtained from the National Weather Service [13], and salt usage data from the U. S. Geological Survey [14].
Salt pollution causes damage to roadside vegetation, material goods, and infrastructure, including roads and bridges, as well as water systems. Salt removal from drinking water supplies requires installation of expensive desalination systems [15]. In the 2017/18 winter season, nearly 22 million metric tons of road salt was applied for de-icing in the USA [16]. Estimated damage due to road salt application was between $803 and $3,401 per ton applied [15] resulting in annual costs of up to $75 billion in the USA.
Negative impacts of salt pollution on human health primarily include the effects from rising sodium concentrations in fresh-water drinking reservoirs, where elevated sodium levels contribute to hypertension and consequently an increased burden on health care systems [17–19].
Salt pollution also affects biodiversity and the composition of microbial communities [20–22] which can indirectly impact human health. Increasing salt levels in freshwater systems favor salt-tolerant cyanobacteria, including some that can produce harmful toxins [22–26]. Furthermore, rising salt levels have been linked to the increased release of toxins by cyanobacteria into the air and water [24, 25, 27, 28]. Consumption of or exposure to these toxins is known to cause liver damage and possible carcinogenic effects [27, 29–33]. Moreover, harmful blooms can adversely affect the tourism and food industries, causing additional economic hardship [34].
Other negative effects of salt pollution include changes in soil biogeochemistry, composition, and characteristics such as decreased aeration and decreased soil permeability with increased erosion [22]. In aquatic environments, salt pollution may impact stratification of lakes and shift biological communities towards more salt tolerant species and reduction of biodiversity [20–22, 35, 36].
Increased salt concentrations heighten osmotic stress which in turn results in inhibition of non-halophilic microbial species and stimulation of halophilic ones. The most sensitive microbes are non-halophilic species which tolerate less than 0.2 M sodium chloride (NaCl). These microbes become osmotically stressed in the presence of even small increases in salt concentrations and may only survive in pockets of lower salinity or in dormant spore forms. Growth of slight halophiles (tolerating up to 0.85 M NaCl), moderate halophiles (tolerating up to 3.4 M NaCl), and extreme halophiles, which require high salt concentrations for growth and tolerate up to 5.2 M (saturated) NaCl, may be promoted when salt concentrations are substantially [37–39].
Sodium chloride is the primary salt used to de-ice roads in the USA, employing halite crystals from salt mines and solar salterns. These have long been known to be rich sources of halophilic microorganisms, harboring a wide variety of species which includes Bacteria and their spores, Archaea, and some Eukaryotes, such as green algae [37, 40–43]. In this initial study, we address whether any halophilic microorganisms are becoming either seasonal or persistent members of proximal microbial communities and if their presence may be used as indicator and biomarkers for salt pollution.
Materials and methods
Sampling and sample sites
Soil samples were collected from 3 sites in Baltimore, Maryland, USA. The sample sites are located in the Jones Falls watershed (Maryland (MD) 8-digit watershed code: 02130904) which is a sub-watershed of the Patapsco River watershed (MD 6-digit watershed code: 021309). Based on frequency of road salt exposure, the sites were designated as (1) un-impacted, (2) seasonally, and (3) continuously impacted. The control un-impacted site (longitude: -76.655922°, latitude: 39.351340°, elevation: 126 m) is located within a public arboretum, approximately 30 m south of a service road and atop a gentle ridge top. The site is approximately 1 m higher than the service road. The seasonally impacted site (longitude: -76.627504°, latitude: 39.347162°, elevation: 84 m) is located in a public park that receives runoff from streets of a residential neighborhood that are salted in winter. The continuously impacted site (longitude: -76.625879°, latitude: 39.316416°, elevation: 20 m) is next to a street running parallel to the Jones Falls. An open salt storage and distribution facility of the City of Baltimore was located across the street from the continuously impacted site at the time of sampling. Road salt samples were obtained from the salt storage facility as well as from a salt box in the residential neighborhood close to the seasonally impacted site (Fig 2).
Sample sites are located in the Jones Falls watershed (MD 8-digit watershed code: 02130904), a sub-watershed of the Patapsco River watershed (MD 6-digit watershed code: 021309). The control un-impacted site (with an osmolality of 48 mOsm kg-1) is located in an arboretum on a ridge-top away from any road. The seasonally impacted site (with an osmolality of 51 mOsm kg-1) is located in a park receiving runoff from streets in a residential neighborhood, and the continuously impacted site (with an osmolality of 603 mOsm kg-1) is located adjacent to a year-round salt storage and distribution facility. The location of the Jones Falls, the main contributor to the Inner Harbor of Baltimore, Maryland, USA, and one of its tributaries, the Stony Run are also indicated. Maps were created using Generic Mapping Tools v. 6.0.0 [44] and GIMP v. 2.8.22. Map data were obtained from Geofabrik GMBH (download.geofabrik.de). Map data OpenStreetMap contributors are licensed under the Open Data Commons Open Database License (ODbL; https://opendatacommons.org/licenses/odbl/1.0/index.html) by the OpenStreetMap Foundation, and the Creative Commons Attribution-ShareAlike 2.0 License (CC BY-SA; https://creativecommons.org/licenses/by-sa/2.0/legalcode).
At each location, samples were taken using a sterile spatula from the top 5 cm of soil (or salt pile) from an area about 10 cm in diameter. A total of approximately 40 grams of soil or salt was collected in sterile plastic sample bags or 50 ml conical tubes and stored at 4°C until further processing. For the isolation and culturing of halophiles, soil samples from the seasonally and continuously impacted sites, and salt samples were collected in February 2013. For the quantification and assessment of soil microbial communities, using spot testing, colony counting, and amplicon sequencing of the 16S rRNA gene, soil surface samples were collected from the control un-impacted, the seasonally impacted, and continuously impacted sites in November 2014, before salt was applied (about 8 months after the last salt application), in February 2015 at the end of the salting season, and in June 2015, close to 3 months after the last salt application.
Salt content of soils was estimated by measuring osmolality of water in which soil was suspended. Soil samples collected in winter 2014, were dried at 50°C for 16 hours, followed by suspension in an equal weight of distilled water for 1 hour to dissolve salts. A subsample was centrifuged at 10,000 x g for 2 minutes and the osmolality of the supernatant measured in a Wescor VAPRO Vapor Pressure Osmometer model 5520 (ELITech Group, Puteaux, France).
Enrichment cultures and isolation of halophiles
Five hundred mg of soil or salt samples were resuspended in halophile broth (HB; 80 mM MgSO4, 10 mM trisodium citrate, 27 mM KCl, 10 g/L Oxoid peptone, 2.3 μM ZnSO4, 0.2 μM MnSO4, 0.2 μM (NH4)2Fe(SO4)2, 56 nM CuSO4, pH 7.2) supplemented with 2 M NaCl or 4.3 M NaCl and incubated at 42°C under constant agitation at 220 RPM in an Innova 44R shaker (New Brunswick, Connecticut, USA) for up to two weeks. Enriched cultures were subsequently streaked onto 2 M or 4.3 M NaCl HB agar plates and incubated at 42°C for up to one week. Single colonies were picked and partial 16S rDNA PCR amplified as described below.
Quantification of halophiles
We developed a semi-quantitative spot test to estimate the abundance of halophiles in soil samples. To assess the accuracy of the spot test, we also quantified the presence of halophiles by serial dilution plating and counting of colony forming units (CFU). The spot test was intended to be a convenient and relatively rapid semi-quantitative tool. We therefore chose to reduce the incubation time of plates to 48 hours at 37°C, a temperature routinely used by laboratories that culture microorganisms. At sites with larger numbers of halophiles, the faster growing halophiles may not only out-compete the slow growing halophiles, but also potentially overgrow the plates making an accurate colony count difficult.
Five tenths cm3 of soil was resuspended in 1 ml HB medium supplemented with 2 M NaCl, and 10-fold serial dilutions (from 1 x 100 to 1 x 10−3) prepared. For the spot test, 20 μl of the serial dilutions were spotted in triplicate onto 2 M NaCl HB agar plates in a grid. To quantify halophiles,100 μl of the serial dilutions were plated in triplicate onto 2 M NaCl HB agar plates. Plates were incubated up to 48 hours at 37°C, examined for microbial growth, and CFU counted.
Nucleic acid extraction
To obtain genomic DNA from microbial isolates, single colonies were picked from HB plates, and grown in 5 ml HB medium supplemented with 2 M or 4.3 M NaCl at 42°C with continuous shaking at 200 RPM in an Innova 44R shaker (New Brunswick, Connecticut, USA). Cells were harvested and DNA extracted as described previously [45]. Genomic DNA from soil and salt samples was extracted using the MO BIO PowerLyzer PowerSoil DNA isolation kit per manufacturer’s recommendation (MO BIO Laboratories, Inc., Carlsbad, California, USA).
PCR amplification and sequencing
For the identification of isolates, partial 16S rDNA covering the variable regions V4 to V9 was amplified with the primers 515F (5’-GTGCCAGCMGCCGCGGTAA-3’) [46] and 1492R (5’-GGTTACCTTGTTACGACTT-3’) [47]. The PCR reactions (25μl) were carried out with 4% (v/v) dimethyl sulfoxide (Sigma-Aldrich, St. Louis, Missouri, USA), 10 mM Tris-HCl, 50 mM KCl, and 1.5 mM MgCl2, 200 μM of each dNTP, 6.25 pM of each primer, and 1.25 U of Taq DNA polymerase (all chemicals from New England Biolabs, Ipswich, Massachusetts, USA). Thermocycler conditions used were 94°C for 2 min, 35 cycles of 92°C for 1 min, 50°C for 45 sec with an extension of 1 sec per cycle, 72°C for 1.5 min, with a final extension at 72°C for 7 min.
PCR products were size separated by electrophoresis on 1% agarose gels and visualized by staining with ethidium bromide (0.5 μg/ml) using a Fotodyne FOTO/Analyst FX system (Fotodyne Inc., Hartland, Wisconsin, USA). Purified PCR products were sequenced using the primers described above. PCR product purification and sequencing services were provided by the Bioanalytical Services Laboratory at the Institute for Marine and Environmental Technology, Baltimore, Maryland, USA.
Sequences were checked for quality and trimmed using trev from the Staden package v. 2.0.0b10 (http://staden.sourceforge.net/) on a Linux Fedora 20 based computer (Red Hat Inc., Raleigh, North Carolina, USA). Trimmed sequences were queried with the Ribosomal Database Project (RDP) Classifier [48] trained against the RDP database 11.3 (released on Sep. 17, 2014; http://rdp.cme.msu.edu/index.jsp) [49].
Amplicon sequencing of the 16S rRNA gene
Amplicon sequencing of the 16S rRNA gene was performed by the Bioanalytical Services Laboratory at the Institute for Marine and Environmental Technology, Baltimore, Maryland, USA on an Illumina MiSeq platform as per manufacturer’s instructions (Illumina, Inc., San Diego, California, USA). The Illumina protocol uses the primer pair Bakt_341F and Bakt_805R29 [50]. Sequences were processed using the open source software suite mothur (v. 1.39.5) following procedures developed for sequences generated by the Illumina MiSeq platform [51–53]. Briefly, paired-end sequences were assembled with a quality threshold of 20 and contigs aligned against a customized reference database that was created based on the SILVA non-redundant SSU Ref 99 database (SILVA SSU Ref NR 99 v132, released on Dec. 13, 2017, [54]). After removal of poor alignments and chimeras, sequences that were 97% or more identical were binned into operational taxonomic units (OTUs). OTUs were queried against the customized reference database using a pseudobootstrap value of 80. Finally, sequences that could not be classified (i.e., “unknown” sequences), as well as sequences identified as eukaryotes, mitochondria, and chloroplasts were removed prior to further analysis. Subsequent analysis of the sequence dataset was performed in R (v. 3.4.1) [55].
To create a customized reference database that was compatible with mothur, the SILVA SSU Ref NR 99 database was processed using scripts published online [56, 57] with slight modification. Briefly, to retain most archaeal and bacterial SSU sequences, instead of trimming the SILVA SSU Ref NR 99 database to a region covering almost the entire SSU, sequences were trimmed starting at position 201 and ending at position 1000 of the SSU (based on the 16S sequence of Escherichia coli [GenBank Acc.# J01859.1]). The resulting reference database retained all archaeal genera and 99.8% of the bacterial genera that are present in the unprocessed SILVA SSU reference database. To ensure high quality alignments of the contigs, the customized reference database was further trimmed to cover only V3 to V5 regions of the SSU. All sequence analyses were performed on an Arch Linux based computer (4–16.7-1-ARCH, https://www.archlinux.org).
Sequence data availability
Raw sequence data were deposited in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/) with the accession number PRJNA522438 and 16S partial sequences of isolates deposited in GenBank (https://www.ncbi.nlm.nih.gov/genbank/) with accession numbers MK543962 to MK543988, MK543991 to MK543994, and MK543996 to MK544007.
Statistical analysis
Statistical analysis was performed in R v. 3.4.1[55]. To compare colony counts, a Conover-Iman test was performed using the R package conover.test v. 1.1.5 [58–60]. The Conover-Iman test was chosen over other possible post hoc non-parametric tests (such as the Dunn’s test) due to its power [59, 60]. The abundance of selected microbial groups was compared using a Fisher’s exact test (part of the R base package [55]). The Fisher’s exact test was chosen over a χ2 analysis as the Fisher’s exact test does not rely on approximations. When post hoc pairwise comparisons were performed, P-values were adjusted with the Holm-Bonferroni correction [61]. For all statistical analysis, the significance level α was set to 0.05. Graphs were generated within R using the packages lattice v. 0.20.35 [62], latticeExtra v. 0.6.28[63], grid v. 3.4.1[55], gridExtra v. 2.2.1[64], and with the GNU Image Manipulation Program (GIMP) v. 2.8.22 (http://www.gimp.org) and Inkscape v. 0.92 (https://inkscape.org/en).
Results
Isolation and classification of halophilic prokaryotes
Halophiles were isolated by preparing enrichment cultures in HB medium containing either 2 M or 4.3 M NaCl followed by multiple rounds of sub-culturing of single colonies on the same agar plates. Based on differences in colony and cellular morphology, a total of 43 isolates were chosen for classification through 16S rDNA sequence analysis. Twenty isolates came from road salt samples (6 and 14 from the salt storage facility, and the salt box, respectively), and 23 from the soil samples (13 from the continuously and 10 from the seasonally impacted site). DNA prepared from candidate halophilic strains and partial 16S rDNA (covering the variable regions V4 to V9) was amplified by PCR with primers specific to prokaryotes [46, 47] followed by Sanger sequencing. Thirty of the isolates were assigned to 4 bacterial genera in the family Bacillaceae: Thalassobacillus (14 isolates), Halobacillus (8 isolates), Piscibacillus (6 isolates), and Pontibacillus (1 isolate), and 1 isolate was identified as Staphylococcus (family: Staphylococcaceae). The remaining 13 isolates belonged to the archaeal genera in the family Haloarchaeaceae/Halobacteriaceae: Natrinema (10 isolates) and Haloterrigena (3 isolates). The Thalassobacillus, Halobacillus, and Natrinema species were isolated from both road salt and soil samples (Table 1).
Culturing was first utilized to detect members of halophile families shown in the table as markers for salt pollution. 16S amplicon sequencing data were later obtained and were confirmatory of culturing results, with the exception of detection of Haloarchaeaceae/Halobacteriaceae and Bacillaceae families in the salt sample of the salt storage facility. Checkmarks indicate which families were detected at the site.
Spot test for a semi-quantitative assessment of halophiles
We developed and evaluated a semi-quantitative spot test to assess the presence of halophiles in soil samples. Soil was suspended in 2 M NaCl HB broth, spotted in triplicate on 2 M NaCl HB media plates, and incubated for 48 hours at 37°C. No growth was observed from samples collected from the un-impacted site. The salt-impacted sites showed growth even at the higher dilutions. The continuously impacted site was found to harbor at least 5 to 10 times more halophiles than the seasonally impacted sites (Fig 3).
Twenty μl of serially diluted soil suspensions were spotted in triplicate onto HB agar plates supplemented with 2 M NaCl and incubated for 48 hours at 37°C. Note that the brown to light brown spots seen at the lower dilutions (100 to 10−1) are residual insoluble materials from the soil suspensions. Colonies appear as cream-colored dots. Un-impacted soil shows no growth on the high salt agar, while seasonally and continuously impacted sites show halophile growth. Continuously impacted sites show the highest growth levels. Summer season shows the least amount of growth across all samples.
Quantification of halophiles
We quantified halophiles over time in soil samples collected from the un-impacted, seasonally impacted, and continuously impacted sites. Soil was suspended in HB broth supplemented with 2 M NaCl, diluted as appropriate, and plated in triplicate onto 2 M NaCl HB media plates. After a 48 hour incubation period, colonies were counted. The continuously impacted site had the highest colony counts, ranging from 39,600 ± 5,546 (standard deviation [SD]) CFU/cm3 soil in summer to 241,333 ± 47,258 (SD) CFU/cm3 soil in fall, followed by the seasonally impacted site with counts ranging from 3,087 ± 247 (SD) CFU/cm3 soil in fall to 6,033 ± 571 (SD) CFU/cm3 soil in winter. Colony counts at the un-impacted site were several orders of magnitude lower, ranging from 60 ± 87 (SD) CFU/cm3 soil in summer to 200 ± 231 (SD) CFU/cm3 soil in fall (Fig 4). In each season, colony counts were significantly higher at the road salt impacted sites, compared to the un-impacted site (adjusted P-values of the seasonal and un-impacted site for fall and summer = 0.0104, and winter = 0.0069; adjusted P-values of the continuously and un-impacted site for fall and summer = 0.0005, and winter = 0.0003; 2-sided Conover-Iman test with Holm-Bonferroni correction). Similarly, the continuously impacted site had significantly higher colony counts than the seasonally impacted site (adjusted P-values for fall and summer = 0.0052, and winter = 0.0035) (Fig 4).
Soils collected in November 2014, February 2015, and June 2015 were resuspended in HB broth supplemented with 2 M NaCl, serial dilutions plated onto HB agar plates supplemented with 2 M NaCl, and plates incubated for 48 hours at 37°C. Colonies were counted and colony forming units normalized to 1 cm3 of soil. The seasonal salt impacted site (blue) harbors significantly higher numbers of halophiles of at least one order of magnitude compared to the un-impacted site (green) at all seasons. Continuous salt exposure enhances the differences as the site that is continuously impacted (red) has at least one order of magnitude more halophiles than the seasonally impacted site (green). The y-axis of the plot is a logarithmic scale. Standard deviations are indicated by error bars.
Amplicon sequencing of the 16S rRNA gene
We assessed the microbial community compositions of salt and soil samples collected in winter 2013 and soil samples collected in the fall 2014, winter 2015, and summer 2015. A total of 458,725 OTUs were classified at the domain level. The majority of the sequences were assigned to the Bacteria (99.66%), with 0.34% of the sequences assigned to the Archaea. Among the Bacteria, 52 phyla were identified. The most abundant phylum was the Proteobacteria (36.59%), followed by the Bacteroidetes and Verrucomicrobia (17.86% and 10.84% of all reads, respectively). Among the Archaea, 4 phyla were present: Euryarchaeota (0.322%), Nanoarchaeaeota (0.016%), Thaumarchaeota (0.004%), and Crenarchaeota (0.001%) (S1 Table).
If, as we hypothesized, halophiles become persistent members of the soil community, more halophilic organisms should be identified in road salt impacted soils. We therefore compared the abundance of representative halophiles identified in soil collected in fall 2014, as well as winter and summer 2015 at the 3 sites. These include Archaea belonging to the classes Haloarchaeaceae/Halobacteria and Nanohaloarchaea, as well as Bacteria of the order Halanaerobiales, and the family Halomonadaceae [37, 65–69]. Members of the Halanaerobiales were not detected in the 2014/15 soil samples. Haloarchaea/Halobacteria and Halomonadaceae were found at all sites. Nanohaloarchaea were found at the seasonally impacted and at the un-impacted sites. The continuously impacted site had significantly more halophiles than the un-impacted and seasonally impacted sites (adjusted P-values = 8.21 x 10−10 and 2.96 x 10−9, respectively; 2-sided Fisher’s exact test with Holm-Bonferroni correction) (Fig 5). Similarly, we observed that Cyanobacteria were significantly more prevalent at the seasonal and continuously impacted sites (adjusted P-values = 2.32 x 10−15 and 1.45 x 10−16, respectively; 2-sided Fisher’s exact test with Holm-Bonferroni correction) (Fig 6).
Amplicon sequencing of the 16S rRNA gene identified 3 prokaryotic groups comprised of predominately halophilic members that were present at the sites. The abundance of these groups was significantly higher at the continuously impacted site (adjusted P-values = 8.21 x 10−10 and 2.96 x 10−9, respectively; 2-sided Fisher’s exact test with Holm-Bonferroni correction). Note that Nanohaloarchaea were not found at the continuously impacted site.
Based on amplicon sequencing of the 16S rRNA gene Cyanobacteria abundance was significantly higher at the seasonally and continuously impacted sites (adjusted P-values = 2.32 x 10−15 and 1.45 x 10−16, respectively; 2-sided Fisher’s exact test with Holm-Bonferroni correction).
Discussion
Our data show that seasonal exposure to road salt leads to an increased presence of halophilic microorganisms, even months after road salt exposure. This finding is supported by our semi-quantitative spot test analysis (Fig 3), as well as by a quantitative measure of colony forming units retrieved from each site (Fig 4). These observations are corroborated by amplicon sequencing of the 16S rRNA gene depicted in Fig 5 and S1 Table. Altogether, this evidence shows the significant impact of road salt exposure on soil microbial communities, leading to a persistent and increased presence of halophilic and/or halotolerant microorganisms.
The use of road salt for de-icing has increased by orders of magnitude over the past decades [2, 9] (Fig 1). This trend is predicted to continue due to climate change effects [12]. Recent models have shown that the decline of Arctic sea ice is leading to changes in wind stream patterns resulting in more persistent snow and ice storms [70–72]. This includes more frequent polar vortex states leading to cold extremes during mid to late winters in mid-latitudes [73]. In addition, population growth and continuing urbanization are leading to an expanding network of roads and highways. Both factors will contribute to increased use of de-icing agents for the foreseeable future and consequently result in microbial communities with more halotolerant and/or halophilic species most likely beyond what we observed in this study.
Road salt pollution and the associated shift in microbial communities have wide-reaching implications for municipalities as they strive to mitigate nutrient loading and effect pollutant reductions in urban runoff. Among the practices used to reduce nutrient input in streams and lakes from urban (and agricultural) runoff are the installation of retention basins, rain gardens, and the expansion of riparian buffers. All these best management practices will, however, potentially be impacted by road salt runoff and buildup. In the laboratory setting, Endreny et al. [21] observed these impacts in their study of microbial communities in bioretention media. After exposure to artificial storm water supplemented with an equivalent of 26.4 mM NaCl (935 mg/L Cl), they reported a reduction in the ability of microbial communities to retain NO3- or PO4. Mesocosm studies also showed that the ecology of storm water retention ponds is constrained by environmentally-relevant concentrations of road salt [74].
Similar to what was observed in mesocosm studies, it is conceivable that the observed changes in microbial communities in soils may result in a reduction or alteration of their ecosystem functions. Therefore, the reported reduced water retention of impacted soils [22], in concert with the possible reduced capacity of nutrient removal of salt contaminated soils and retention media, may impede the ability of storm water management structures to reduce nutrient input into urban water bodies. Increased salinity in combination with high nutrient levels in aquatic environments has been shown to lead to more frequent occurrence of algae blooms, especially cyanobacterial blooms, with increased release of cyanotoxins and considerable environmental and public health concerns [24–33]. In addition to the indirect effects that microbial communities of salt contaminated soils may have, the increased sodium levels in fresh water reservoirs may pose a direct threat to human health as sodium levels continue to rise in affected urban fresh water reservoirs [7, 17, 18, 75, 76]. In this context, our results raise concerns regarding the possibility of toxic cyanobacteria resulting from increased salinity in soils (Fig 6).
Our results demonstrate that halophilic microorganisms become persistent members of microbial communities as a result of salting roads for de-icing during winter months. In order to quantify increases in salinity and halophiles in the environment, we have developed and evaluated a facile and rapid spot test that utilizes a customized media formulation for capturing and cultivating halophilic microorganisms. Moreover, we show using this assay, that soils impacted by road salt harbor 10-fold more halophiles than un-impacted soils. As such, the spot test has the potential to be a valuable tool to survey and evaluate salt pollution of soils.
Supporting information
S1 Table. Abundance of prokaryotes at the domain and phylum level.
https://doi.org/10.1371/journal.pone.0221355.s001
(XLSX)
Acknowledgments
This work was supported by private donations through the HERD fund (https://medschool-umaryland.givecorps.com/projects/3511-school-of-medicine-research-herd-fund).
References
- 1. Kaushal SS, Groffman PM, Likens GE, Belt KT, Stack WP, Kelly V2R, et al. Increased salinization of fresh water in the northeastern United States. Proc Nat Acad Sci USA. 2005;102:13517–20. pmid:16157871
- 2.
Mullaney JR, Lorenz DL, Arntson AD. Chloride in groundwater and surface water in areas underlain by the glacial aquifer system, northern United States. Reston, VA, USA: U. S. Geological Survey; 2009. 41p. Report No.: U. S. Geological Survey Scientific Investigations Report 2009–5086.
- 3. Corsi SR, Graczyk DJ, Geis SW, Booth NL, Richards KD. A fresh look at road salt: aquatic toxicity and water-quality impacts on local, regional, and national scales. Environ Sci Technol. 2010;44:7376–82. pmid:20806974
- 4. Van Meter RJ, Swan CM, Snodgrass JW. Salinization alters ecosystem structure in urban stormwater detention ponds. Urban Ecosyst. 2011;14:723–36.
- 5. Cooper CA, Mayer P., and Faulkner B. R. Effects of road salts on groundwater and surface water dynamics of sodium and chloride in an urban restored stream. Biogeochemistry. 2014;121:149–66.
- 6. Corsi SR, De Cicco LA, Lutz MA, Hirsch RM. River chloride trends in snow-affected urban watersheds: increasing concentrations outpace urban growth rate and are common among all seasons. Sci Total Environ. 2015;508:488–97. pmid:25514764
- 7. Dugan HA, Bartlett SL, Burke SM, Doubek JP, Krivak-Tetley FE, Skaff NK, et al. Salting our freshwater lakes. Proc Nat Acad Sci USA. 2017;114:4453–8. pmid:28396392
- 8. Jackson RB, Jobbágy EG. From icy roads to salty streams. Proc Nat Acad Sci USA. 2005;102:14487–8. pmid:16203970
- 9. Kim S-y, Koretsky C. Effects of road salt deicers on sediment biogeochemistry. Biogeochemistry. 2013;112:343–58.
- 10.
Gutowski WJJ, Hegerl GC, Holland GJ, Knutson TR, Mearns LO, Stouffer RJ, et al. Causes of observed changes in extremes and projections of future changes. In: Karl TR, Meehl GA, Miller CD, Hassol SJ, Waple AM, Murray WL, editors. Weather and Climate Extremes in a Changing Climate. Regions of Focus: North America, Hawaii, Caribbean, and U. S. Pacific Islands: U. S. Climate Change Science Program and the Subcommittee on Global Change Research, Washington, DC, USA; 2008. p. 81–114.
- 11. López-Moreno JI, Goyette S, Vicente-Serrano SM, Beniston M. Effects of climate change on the intensity and frequency of heavy snowfall events in the Pyrenees. Clim Change. 2011;105:489–508.
- 12. Notaro M, Lorenz D, Hoving C, Schummer M. Twenty-first-century projections of snowfall and winter severity across central-eastern North America. J Clim. 2014;27:6526–50.
- 13. The National Weather Service (NWS). Monthly and Annual Snowfallbat Central Park [Internet]. 2018 [updated 2018 Aug 24; cited 2018 Aug 27]. Available from: https://www.weather.gov/media/okx/Climate/CentralPark/monthlyseasonalsnowfall.pdf.
- 14. U. S. Geological Survey. Salt: Statistics and Information [Internet]. 2018 [updated 2018 Feb 8; cited 2018 Aug 20]. Available from: https://minerals.usgs.gov/minerals/pubs/commodity/salt/index.html#myb.
- 15.
Fortin Consulting Inc. The real cost of salt use for winter maintenance in the twin cities metropolitan area. Hamel, MN, USA: Minnesota Pollution Control Agency; 2014. 21p.
- 16.
U. S. Geological Survey. Mineral commodity summaries 2018. Report. Reston, VA, USA: U. S. Geological Survey; 2018. 200p. https://doi.org/10.3133/70194932
- 17.
U. S. Environmental Protection Agency. Drinking water advisory: consumer acceptability advice and health effects analysis on sodium. Washington, DC, USA: U. S. Environmental Protection Agency Office of Water (4304T), Health and Ecological Criteria Division; 2003. 34p. Report No.: EPA 822-R-03-006.
- 18. Hallenbeck WH, Brenniman GR, Anderson RJ. High sodium in drinking water and its effect on blood pressure. Am J Epidemiol. 1981;114:817–26. pmid:7315830
- 19.
World Health Organization and International Programme on Chemical Safety. Guidelines for Drinking-Water Quality. Vol. 2, Health criteria and other supporting information 2nd ed. Geneva, Switzerland: World Health Organization; 1996. 973p. Available from: http://www.who.int/iris/handle/10665/38551
- 20. Černohlávková J, Hofman J, Bartoš T, Sáňka M, Anděl P. Effects of road deicing salts on soil microogranisms. Plant Soil Environ. 2008;54:479–85.
- 21. Endreny T, Burke DJ, Burchhardt KM, Fabian MW, Kretzer AM. Bioretention column study of bacteria community response to salt-enriched artificial stormwater. J Environ Qual. 2012;41:1951–9. pmid:23128752
- 22. Ramakrishna DM, Viraraghavan T. Environmental impact of chemical deicers—a review. Water Air Soil Pollut. 2005;166:49–63.
- 23.
Briggins D, Walsh J. The environmental implications of road salting in Nova Scotia. Nova Scotia, Canada: Nova Scotia Department of the Environment and Nova Scotia Department of Transportation and Communication; 1989. 80p.
- 24. O’Neil JM, Davis TW, Burford MA, Gobler CJ. The rise of harmful cyanobacteria blooms: the potential roles of eutrophication and climate change. Harmful Algae. 2012;14:313–34.
- 25. Tonk L, Bosch K, Visser PM, Huisman J. Salt tolerance of the harmful cyanobacterium Microcystis aeruginosa. Aquat. Microb. Ecol. 2007;46:117–23.
- 26. Rastogi RP, Madamwar D, Incharoensakdi A. Bloom dynamics of cyanobacteria and their toxins: environmental health impacts and mitigation strategies. Front Microbiol. 2015;6:1254. Available from: pmid:26635737
- 27. Beversdorf JL, Weirich AC, Bartlett LS, Miller RT. Variable cyanobacterial toxin and metabolite profiles across six eutrophic lakes of differing physiochemical characteristics. Toxins. 2017;9:62. Available from: https://doi.org/10.3390/toxins9020062.
- 28. Rosen BH, Loftin KA, Graham JL, Stahlhut KN, Riley JM, Johnston BD, et al. Understanding the effect of salinity tolerance on cyanobacteria associated with a harmful algal bloom in Lake Okeechobee, Florida. U. S. Geological Survey Scientific Investigations Report 2018–5092. 2018:32p.
- 29. U. S. Environmental Protection Agency. Cyanobacteria/Cyanotoxins [Internet]. 2018 [cited 2018 Jan 9]. Available from: https://www.epa.gov/nutrient-policy-data/cyanobacteriacyanotoxins.
- 30. Zhang F, Lee J, Liang S, Shum CK. Cyanobacteria blooms and non-alcoholic liver disease: evidence from a county level ecological study in the United States. Environ Health. 2015;14:41. Available from: pmid:25948281
- 31. Virginia Department of Public Health. Cyanobacteria [Internet]. 2018 [cited 2018 Aug 1]. Available from: http://www.vdh.virginia.gov/environmental-epidemiology/waterborne-hazards-control-program/cyanobacteria/.
- 32. Murby AL, Haney JF. Field and laboratory methods to monitor lake aerosols for cyanobacteria and microcystins. Aerobiologia. 2016;32:395–403.
- 33. U. S. Environmental Protection Agency. Health and ecological effects [Internet]. 2017 [cited 2018 Aug 1]. Available from: https://www.epa.gov/nutrient-policy-data/health-and-ecological-effects.
- 34. Carmichael WW, Boyer GL. Health impacts from cyanobacteria harmful algae blooms: implications for the north american Great Lakes. Harmful Algae. 2016;54:194–212. pmid:28073476
- 35. Findlay SEG, Kelly VR. Emerging indirect and long-term road salt effects on ecosystems. Ann NY Acad Sci. 2011;1223:58–68. pmid:21449965
- 36. Tiquia SM, Davis D, Hadid H, Imail M, Sahly R, Shim J, et al. Halophilic and halotolerant bacteria from river waters and shallow groundwater along the Rouge River of southeastern Michigan. Environ Technol. 2007;28:297–307. pmid:17432382
- 37. DasSarma S, DasSarma P. Halophiles. eLS. 2017. Available from: https://doi.org/10.1002/9780470015902.a0000394.pub4.
- 38. Larsen H. Halophilic and halotolerant microorganisms—an overview and historical perspective. FEMS Microbiol Lett. 1986;39:3–7.
- 39.
Oren A. Halophilic Microorganisms and their Environment. Dordrecht, Netherlands: Kluver Academic Publishers; 2002. 558p. (Seckbach J, editor. Cellular origin, life in extreme habitats and astrobiology; vol. 5).
- 40. Grant WD. Life at low water activity. Phil Trans R Soc Lond B. 2004;359:1249–67.
- 41. Norton CF, McGenity TJ, Grant WD. Archaeal halophiles (Halobacteria) from two British salt mines. J Gen Microbiol. 1993;139:1077–81.
- 42. Gramain A, Díaz GC, Demergasso C, Lowenstein TK, McGenity TJ. Archaeal diversity along a subterranean salt core from the Salar Grande (Chile). Environ Microbiol. 2011;13:2105–21. pmid:21355972
- 43. DasSarma P, Klebahn G, Klebahn H. Translation of Henrich Klebahn's 'Damaging agents of the klippfish—a contribution to the knowledge of the salt-loving organisms'. Saline Systems. 2010;6:7. Available from: pmid:20569477
- 44. Wessel P, Smith WHF, Scharroo R, F. LJ, Wobbe F. Generic Mapping Tools: improved version released. EOS Trans AGU. 2013;94:409–10.
- 45.
DasSarma S, Robb FT, Place AR, Sowers KR, Schreier HJ, Fleischmann EM, editors. Archaea: a laboratory manual—halophiles. Cold Spring Harbor, NY, USA: Cold Spring Harbor Laboratory Press; 1995.
- 46. Blank CE, Cady SL, Pace NR. Microbial composition of near-boiling silica-depositing thermal springs throughout Yellowstone National Park. Appl Environ Microbiol. 2002;68:5123–35. pmid:12324363
- 47. DeLong EF. Archaea in coastal marine environments. Proc Natl Acad Sci USA. 1992;89:5685–9. pmid:1608980
- 48. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian Classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7. pmid:17586664
- 49. Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, et al. Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucl Acids Res. 2014;42:D633–D42. pmid:24288368
- 50. Herlemann DP, Labrenz M, Jürgens K, Bertilsson S, Waniek JJ, Andersson AF. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 2011;5:1571–9. pmid:21472016
- 51. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41. pmid:19801464
- 52. Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20. pmid:23793624
- 53.
Mothur Wiki. MiSeq SOP. [Internet]. [cited 2017 Dec 20]. Available from: https://www.mothur.org/wiki/MiSeq_SOP.
- 54.
SILVA rRNA database project. Release information: SILVA 132 [Internet]. [cited 2017 Dec 19]. Available from: https://www.arb-silva.de/documentation/release-132/.
- 55.
R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. 2017. Available from: https://www.r-project.org.
- 56. Kerckhof FM. README for SILVA v128 reference files [Internet]. 2017 [updated 2017 Feb 9; cited 2017 Dec 20]. Available from: http://users.ugent.be/~fpkerckh/silvav128_README.html.
- 57. Schloss PD. The mothur blog: README for the SILVA v128 reference files [Internet]. 2017 [updated 2017 Mar 22; cited 2017 Dec 20]. Available from: http://blog.mothur.org/2017/03/22/SILVA-v128-reference-files/.
- 58.
Dinno A. conover.test: Conover-Iman test of multiple comparisons using rank sums. R package version 1.1.5. 2017. Available from: https://CRAN.R-project.org/package=conover.test.
- 59.
Conover WJ. Practical Nonparametric Statistics. 3rd ed. Hoboken, NJ: Wiley; 1998. 592p.
- 60.
Conover WJ, Iman RL. On multiple-comparisons procedures. Los Alamos, NM, USA: Los Alamos Scientific Laboratory; 1979. 14p. Report No.: LA-7677-MS.
- 61. Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6:65–70.
- 62.
Sarkar D. Lattice: Multivariate Data Visualization with R. New York, NY, USA: Springer-Verlag; 2008. 273p.
- 63. Sarkar D, Andrews F. latticeExtra: Extra Graphical Utilities Based on Lattice. 2016. Available from: http://CRAN.R-project.org/package=latticeExtra.
- 64. Auguie B. gridExtra: Miscellaneous Functions for "Grid" graphics. R package version 2.2.1. 2016. Available from: http://CRAN.R-project.org/package=gridExtra.
- 65. Oren A. Microbial life at high salt concentrations: phylogenetic and metabolic diversity. Saline Systems. 2008;4:1–13. Available from:
- 66.
Oren A. Diversity of Halophiles. In: Horikoshi K, Antranikian G, Bull AT, Robb FT, Stetter KO, editors. Extremophiles Handbook. Tokyo, Japan: Springer; 2011. p. 309–25.
- 67.
de la Haba RR, Sánchez-Porro C, Marquez MC, Ventosa A. Taxonomy of Halophiles. In: Horikoshi K, Antranikian G, Bull AT, Robb FT, Stetter KO, editors. Extremophiles Handbook. Tokyo, Japan: Springer; 2011. p. 255–308.
- 68. Ghai R, Pašić L, Fernández AB, Martin-Cuadrado A-B, Mizuno CM, McMahon KD, et al. New abundant microbial groups in aquatic hypersaline environments. Sci Rep. 2011;1:135. Available from: pmid:22355652
- 69. Narasingarao P, Podell S, Ugalde JA, Brochier-Armanet C, Emerson JB, Brocks JJ, et al. De novo metagenomic assembly reveals abundant novel major lineage of Archaea in hypersaline microbial communities. ISME J. 2011;6:81–93. pmid:21716304
- 70. Alexander MA, Bhatt US, Walsh JE, Timlin MS, Miller JS, Scott JD. The atmospheric response to realistic Arctic sea ice anomalies in an AGCM during winter. J Clim. 2004;17:890–905.
- 71. Singarayer JS, Bamber JL, Valde PJ. Twenty-first-century climate impacts from a declining Arctic sea ice cover. J Clim. 2006;19:1109–25.
- 72. Liu J, Curry JA, Wang H, Song M, Horton RM. Impact of declining Arctic sea ice on winter snowfall. Proc Natl Acad Sci USA. 2012;109:4074–9. pmid:22371563
- 73. Kretschmer M, Coumou D, Agel L, Barlow M, Tziperman E, Cohen J. More-persistent weak stratospheric polar vortex states linked to cold extremes. Bull Amer Meteor Soc. 2018;99:49–60.
- 74. Meter RJV, Swan CM. Road salts as environmental constraints in urban pond food webs. PLoS ONE. 2014;9(2):e90168. Available from: pmid:24587259
- 75. Kelly WR, Panno SV, Hackley KC. Impacts of road salt runoff on water quality of the Chicago, Illinois, region. Environ Eng Geosci. 2012;18:65–81.
- 76. Kelly WR. Long-term trends in chloride concentrations in shallow aquifers near Chicago. Ground Water. 2008;46:772–81. pmid:18624692