Assessing Different Mechanisms of Toxicity in Mountaintop Removal/Valley Fill Coal Mining-Affected Watershed Samples Using Caenorhabditis elegans

Abstract

Mountaintop removal-valley fill coal mining has been associated with a variety of impacts on ecosystem and human health, in particular reductions in the biodiversity of receiving streams. However, effluents emerging from valley fills contain a complex mixture of chemicals including metals, metalloids, and salts, and it is not clear which of these are the most important drivers of toxicity. We found that streamwater and sediment samples collected from mine-impacted streams of the Upper Mud River in West Virginia inhibited the growth of the nematode Caenorhabditis elegans. Next, we took advantage of genetic and transgenic tools available in this model organism to test the hypotheses that the toxicity could be attributed to metals, selenium, oxidative stress, or osmotic stress. Our results indicate that in general, the toxicity of streamwater to C. elegans was attributable to osmotic stress, while the toxicity of sediments resulted mostly from metals or metalloids.

About the Authors:

Elena A. Turner

Affiliation: Nicholas School of the Environment, Duke University, Durham, North Carolina, United States of America

Gretchen L. Kroeger

Affiliation: Nicholas School of the Environment, Duke University, Durham, North Carolina, United States of America

Mariah C. Arnold

Affiliation: Nicholas School of the Environment, Duke University, Durham, North Carolina, United States of America

B. Lila Thornton

Affiliation: Nicholas School of the Environment, Duke University, Durham, North Carolina, United States of America

Richard T. Di Giulio

Affiliation: Nicholas School of the Environment, Duke University, Durham, North Carolina, United States of America

Joel N. Meyer

* E-mail: [email protected]

Affiliation: Nicholas School of the Environment, Duke University, Durham, North Carolina, United States of America

Introduction

Mountaintop removal/valley fill coal mining is a form of surface mining wherein all vegetation and topsoil are removed, overlying rock is dynamited to expose coal seams, and spoilage is moved into adjacent valleys, thus burying streams [1]. While mountaintop removal/valley fill coal mining (herein abbreviated as MTR/VF) offers a lower-cost way to extract coal than traditional shaft mines, full risk-benefit assessments of this practice require an understanding of potential ecological and human health effects [2]. However, despite the fact that a thorough understanding of those health effects is currently lacking, MTR/VF is already occurring on a very large scale. Lindberg et al. recently estimated the total extent of surface mining as 6.8% of Appalachian coal fields, or 4.68 million square hectares [3].

Evaluating the potential toxic effects of MTR/VF is complicated by the facts that the effluent from such sites is highly complex and variable in space and time. The “TRIAD” approach, which utilizes a combination of chemistry, bioassays and in situ studies, has been proposed to assess the potential health effects of such complex environmental pollutant mixtures. This approach takes advantage of the complementary strengths of field surveys, analytical chemistry, and in-laboratory analysis of the toxicity of field-collected samples [4]. Field surveys of the environmental effects of MTR/VF on stream ecology and benthic invertebrate communities are emerging [5-7], and surface mining has also been linked to negative impacts on human health [8,9]. Analytical chemistry approaches have documented high levels of inorganic solutes in waters and sediments downstream from MTR/VF sites [3,5,10,11]. However, the body of laboratory-based research investigating the toxicity of mine-impacted water and stream sediment is smaller [3,11,12], and currently lacks a strong mechanistic component.

To strengthen our understanding of the potential toxicity of MTR/VF-associated water and sediment samples, we utilized Caenorhabditis elegans. This genetically well-characterized model organism is being increasingly used in environmental toxicology research to gain mechanistic insight into environmental samples which may operate via one or several different mechanisms of toxicity [13-15]. C. elegans normally inhabits moist, microbe-rich soil surface environments [16], but is routinely cultured for multiple generations in liquid medium [17]and frequently used for aquatic toxicity tests [12,14,18,19]. Furthermore, while there are many other invertebrate species available for aquatic toxicology testing, none offer the same genetic tools. C. elegans shares stress-response pathways with higher eukaryotes [20,21], and many mutant and reporter strains that facilitate mechanistic analysis are available and widely used. For example, if knockout mutants that do not produce crucial metal-chelating peptides are hyper-sensitive to exposure to a particular environmental mixture, this suggests that this mixture is causing toxicity in part due to metal content. The same approach can be used to diagnose exposures that cause toxicity via DNA damage, oxidative stress, and other mechanisms using the appropriate mutant strains [19,22-24], an approach that has been dubbed “functional toxicology” [25]. Similarly, transgenic strains that carry a “reporter” such as green fluorescent protein (GFP) driven by the promoter region of stress-response genes can provide a stress-response endpoint at lower thresholds than those that cause outright toxic effects [12,26].

We used C. elegans to investigate the toxicity and mechanism of toxicity of environmental samples from the Mud River and several associated tributaries (Boone County, Southwestern West Virginia). The Mud River flows through the Hobet 21 surface mine from which it receives numerous discharges, both directly and via its tributaries; the sampling sites have previously been characterized with respect to water chemistry and show both conductivity in excess of the proposed EPA benchmark of 300µs/cm-1 and selenium levels in excess of EPA water quality criteria of 5µG Se L-1 [3]. Previous work has suggested that increased levels of metals and total dissolved solids (TDS) may be important stressors in MTR/VF systems [6]. We hypothesized that water and sediment samples from putatively mining-impacted streams would affect the growth of wild-type nematodes, and that mutant nematodes would be differentially affected based on the composition of the stream water and sediment pore water. Furthermore, we analyzed mutant strains that would help us distinguish the contributions of metals or metalloids, osmotic stress and oxidative stress to the toxicity of these samples (Table 1).

Strain Gene Knockout Mechanism

gpdh-1(ok1558); gpdh-2(kb33) glycerol-3 phosphate dehydrogenase sensitive to osmotic stress

mtl-2(gk125) metallothionein sensitive to heavy metals

mtl-2::GFP reporter none produces GFP-tagged metallothionein

N2 none wild-type

pcs-1(tm1748) phytochelatin synthase sensitive to metals and metalloids

smf-2(gk133) divalent metal transporter sensitive to manganese exposure

sod-3(gk257) superoxide dismutase sensitive to oxidative stress

Table 1. Strains.

CSV

Download CSV

Stream water and sediment pore water from putatively mine-impacted sites along the Mud River and its tributaries had significant impacts on the growth of wild-type and mutant nematodes. Effects on growth ranged from positive to highly negative and varied between site, strain and season; additionally, stream water and sediment pore water from the same site often produced different effects on growth.

Materials and Methods

Sample Collection and Handling

Sample collection took place June 24 2010, December 8-10 2010, April 12-15 2011, and May 1 2012 at sites shown in Figure 1. Samples from Mud River 7 were taken on private property with the owner’s permission. All other samples were taken from easements of publicly accessible roads in accordance with local regulations. Rivers are designated as public property in West Virginia. Our field sampling did not involve endangered or protected species. Stream water samples were collected in metal-free 50 ml polypropylene tubes; sediment was collected in glass mason jars using a small plastic trowel. Stream water conductivity was measured during sediment and water collection. Water and sediment samples were stored on ice while in the field and transferred to a 4 °C cold room. Sediment was passed through a 2mm sieve and frozen in 5 to 7 ml aliquots at -80 °C until use. Pore water was prepared by thawing sediment, adding an equal volume of ddH 20, shaking vigorously for two minutes, and centrifuging at 6987xg for 10 minutes. The supernatant (“pore water”) was removed and used for growth assays; all pore water was stored at 4 °C until use and replicates were run within 4 to 6 weeks.

[Figure omitted. See PDF.]

Figure 1. Map of the Mud River watershed.

Sampling sites plus age and extent of surface mining activity. Insets show sequentially larger geographic regions.

https://doi.org/10.1371/journal.pone.0075329.g001

Nematode Culture

Nematodes were maintained at 20 or 15 °C on K agar plates with a lawn of OP 50 strain Escherichia coli [27]. Synchronized L1 larvae were obtained by dissolving gravid adults with a solution of 5% sodium hypochlorite and hatching eggs in K medium plus MgSO4, CaSO4 and cholesterol. The pcs-1 (tm1748 allele, outcrossed 6x) mutant strain was obtained from Stephen Sturzenbaum (King’s College, London). The gpdh-1 and gpdh-2 (ok1558 and kb33 alleles, outcrossed 3 times) mutant strains were obtained from Todd Lamitina (University of Pennsylvania). gpdh-1(ok1558) males were generated by heat-shock and crossed with gpdh-2(kb33) hermaphrodites to generate gpdh-1(ok1558);gpdh-2(kb33) double mutants; homozygosity for both deletions was confirmed by PCR using primers designed using Wormbase (Gpdh-1 forward primer: 5’ ccctcttcaggattcttccc 3’, reverse primer 5’ tccgctagatccatttccag 3’; Gpdh-2 forward primer 5’ cgcgaacaaagaatatttaggaa 3’, reverse primer 5’ aattcgtcaaaggaatctgtgaa 3’). The mtl-2 mutant (JF23 strain, gk125 allele, outcrossed 4 times) was obtained from Jonathan Freedman (NIEHS). N2 (Bristol) wild-type and smf-2 (VC171 strain, gk133 allele) and sod-3 (VC433 strain, gk235 allele) mutant strains were obtained from the Caenorhabditis Genetics Center (CGC, Minneapolis, MN, USA). The sod-3 strain was outcrossed 3 times with N2 (Bristol wildtype) males; genotyping was performed using primers designed using Wormbase (sod-3 forward primer: 5’ gtcccgaaatgcattttttc 3’, reverse primer 5’ gggcggacatttttgtactg 3’). The mtl-2::GFP reporter gene strain was obtained from Phillip Williams (University of Georgia).

Growth assays

On day 0 of each growth assay, synchronized L1 larvae were distributed onto 96-well plates (n=50) with a COPAS Biosort (Union Biometrica, Holliston MA) using 50 µL 2x EPA reconstituted moderately hard water [22]. An additional 50 µL of dosing solution composed of stream water, sediment pore water, or sodium selenite at 2x final dosing strength plus cholesterol (final concentration of 0.01 mg/mL) was added to the wells, and plates were incubated at 20 degrees for the course of the assay. Dosing solution was renewed daily. Worms were fed UVC-killed E. coli (UVRA strain) to eliminate the potentially confounding effect of bacterial metabolism on exposures, as previously described [19]. Bacteria were concentrated to approximately 1x1010 cells/mL [28], at 1:100 by volume for the first day, 1:50 for the second day and 1:10 for the final day of the growth assay. Increasing concentrations of bacteria were added to accommodate the increasing food requirements of the growing nematodes. Growth of the worms was measured using a COPAS Biosort, using extinction (optical density of the worm) as a growth endpoint [22]. Extinction data from the third day of exposure was graphed as a percentage of the growth of EPA water exposed controls. For the purposes of this paper we define 3-day exposures as medium-term, because they encompass most of the developmental trajectory of this species.

Reporter strain assay

50 L4/young adult (staged as having minimally reached the L4 stage based on appearance of the vulval crescent, but not having reached reproductive adulthood based on absence of eggs before and after exposure) mtl-2::GFP nematodes were distributed into 24 well plates containing sediment pore water, EPA water, 0.5µM cadmium chloride, 5mM manganese chloride or 0.15µM sodium selenite for a total volume of 500µL, with no food. Plates were incubated at 20 °C for 24 hours, and fluorescence was measured with the COPAS Biosort and graphed as a ratio of GFP expression to extinction (as a proxy for size).

Statistical analysis

Size (growth) and fluorescence measurements were analyzed via Mann Whitney U tests conducted using Statview for Windows (Version 5.0.1, SAS Institute Inc., Cary, NC). Bonferroni corrections were applied to p value cut-offs to determine significance. For example, for Figure 1a, 13 exposures were compared to the EPA water control, so that the cut-off for significance was p < 0.0038 (0.05/13) rather than p < 0.05.

Results

I: Mud River and tributary samples affected the growth of wild-type nematodes

We exposed wild-type nematodes to water and sediment pore water from one non-mining impacted site and multiple impacted sites shown in Figure 1. Both stream water and sediment pore water from mining-impacted sites on the Mud River and its tributaries reduced the growth of wild-type nematodes in comparison to EPA water controls (Figure 2). Lethality was not observed for any of the exposure conditions. Reduction in growth varied between sites but was statistically significant for all water samples from mining-impacted sites and three sediment pore water samples from mining-impacted sites sampled in April 2011 (Figure 2 and Table 2). Water from a reference site on the unmined left fork of the Mud River did not reduce the growth of wild-types nematodes; in fact, there was a slight positive effect on growth in wild-type nematodes exposed to sediment pore water from the reference site in comparison to EPA water-treated controls. Because growth of wild-type nematodes was better in the reference site water than in EPA MHRW, we carried out statistical analyses for all experiments with wild-type nematodes using the both the reference site and EPA MHRW for comparison (p values reported in Table 2 and Table S1, respectively).

[Figure omitted. See PDF.]

Figure 2. Growth effects of stream water and sediment pore water.

The growth of wild-type nematodes was affected by stream water (2A, n=27-359) and sediment pore water (2B, n=30-359) from the Mud River and its tributaries, but not by water or sediment pore water from a reference site on the unmined Left Fork. Statistically significant difference from EPA water control indicated by asterisks. “Size” refers to the optical density (extinction) of each nematode and describes growth since the data shown is for size after three days of exposure. Each exposure was repeated at least two times separately. EPA water was used as a laboratory control medium.

https://doi.org/10.1371/journal.pone.0075329.g002

April 2011

Site Water Sediment Pore Water Water vs. Sediment Pore Water

Left Mud River (reference) − − 0.1775

Mud River 5 <0.0001 <0.0001 0.3939

Stanley Fork <0.0001 <0.0001 <0.0001

Mud River 5a <0.0001 <0.0001 0.0237

Mud River 6 <0.0001 <0.0001 <0.0001

Laurel Branch <0.0001 <0.0001 0.8584

Mud River 7 <0.0001 <0.0001 0.7772

Mud River 8b <0.0001 <0.0001 0.0741

Connelly Branch <0.0001 <0.0001 ----

Mud River 8 <0.0001 <0.0001 0.0042

Mud River 10a <0.0001 <0.0001 <0.0001

Mud River 10 <0.0001 <0.0001 0.0001

Berry Branch <0.0001 <0.0001 0.7691

p-value cutoff for statistical significance 0.0042 0.004545 0.0042

Table 2. p values, Mann Whitney U test in comparison to reference site.

Comparison of optical density (extinction) of wild-type nematodes treated with steam water and sediment pore water to EPA water-treated controls. The Bonferroni-corrected p-value cut-off for significance based on 12 pairwise comparisons is p < o.oo42.

CSV

Download CSV

The fact that C. elegans grew less well in EPA MHRW than in water from our reference site is not surprising. Other researcher have found that the common aquatic toxicology species Hyallela azteca, for example, grows relatively poorly in this medium, perhaps due to insufficient chloride [29,30], although bromide [31] or other elements may also be important. Since selenium is a known contaminant of the Mud River [3] and a micronutrient [32], we next investigated its effects on the growth of wild-type nematodes.

II: Selenium stimulates the growth of C. elegans at low doses, but is toxic at higher doses

Though the toxic effects of the metalloid selenium in aquatic invertebrates are well-documented [6,33-35], there is considerably less published work on the toxicity of selenium to nematodes. The mechanism of toxicity of selenium in nematodes is not well studied, though oxidative stress-induced cholinergic neurodegeneration has been demonstrated [36]. While toxic at moderate to high doses, selenium is also an essential micronutrient and previous research demonstrated positive effects on the growth and development of nematodes exposed to very low concentrations of selenite [32,37]. To test whether the same phenomenon could be occurring under our experimental conditions, we performed growth assays exposing worms to a wide range of concentrations of sodium selenite. Very low concentrations (0.0001-0.00005µM) of sodium selenite caused a small but significant increase in growth in treated worms compared to untreated controls (Figure 3A). Positive effects on growth became less pronounced as concentration increased; concentrations in the 0.15625 to 0.78125µM range caused highly variable growth and concentrations above 0.5µM caused total growth arrest or death (lethal doses not shown on graph).

[Figure omitted. See PDF.]

Figure 3. Growth of wild-type and pcs-1 knockout nematodes was affected by sodium selenite.

The growth of wild-type (3A, n=49-247) and pcs-1(tm1748) knockout (3B, n=30-279) nematodes was impacted both positively and negatively by addition of sodium selenite to EPA water, according to a non-monotonic dose-response curve. Statistically significant difference from EPA water control indicated by asterisks. Negative effects on growth were seen at lower doses in pcs-1(tm1748) mutants than in wild-type worms, and positive effects were not observed in the pcs-1(tm1748) strain. “Size” refers to the optical density (extinction) of each nematode and describes growth since the data shown is for size after three days of exposure. Each exposure was repeated at least two times separately. Please note the difference in Y-axis scaling.

https://doi.org/10.1371/journal.pone.0075329.g003

Since subsequent tests for toxicity were performed using the pcs-1 mutant strain which cannot synthesize phytochelatin, and since phytochelatin has been shown to bind selenium [38], we also tested the effect of selenite exposure in the pcs-1(tm1748) deletion mutants. Positive growth effects were not statistically significant for most low concentrations in pcs-1(tm1748) deficient nematodes treated with sodium selenite. Overall variability in growth was greater at all concentrations in pcs-1(tm1748) than wild-type nematodes, and overt toxicity was observed at lower concentrations in pcs-1(tm1748) than in wild-type worms (Figure 3b). Comparable concentrations of sodium selenate caused no effect on the growth of either wild-type or pcs-1(tm1748) nematodes (data not shown), in keeping with the relative toxicity of these two Se species in most organisms [39]. We note that the nematodes were fed UVC-inactivated bacteria, so that alteration of the selenium to typically more-toxic organic forms [39] is unlikely. Pilot studies with autoclaved samples ruled out an alternate food source as the cause of positive growth effects (data not shown).

III: Tributaries increased the toxicity of the main stem

We hypothesized that the toxicity of the water and sediment pore water in Mud River was influenced not only by discharge points along the main stem but also by contributions from its mining-impacted tributaries. Comparison of growth data from two different sampling trips supported this hypothesis; stream water from April 2011 (n=27-53) and sediment pore water from December 2010 (n=98-163) from Mud River 5a, downstream from the input of the mining-impacted Stanley Fork, caused greater growth inhibition in wild-type nematodes than stream water or sediment pore water, respectively, from the upstream site of Mud River 5 (Figure 4, p=0.0013 for April 2011, p<0.0001 for December 2010. Stream water and sediment pore water collected from the same site often caused different effects on growth in wild-type nematodes (see Table 2). Potential explanations for these differences include different quantities of potential toxicants, different mechanisms of toxicity, or changing concentration of contaminants resulting from rain events or drought conditions. Having quantified growth inhibition in wild-type nematodes, we next conducted investigations into the mechanism(s) of toxicity using GFP reporter strains and knockout mutants.

[Figure omitted. See PDF.]

Figure 4. Mining-impacted tributaries affect the toxicity of the Mud River.

Contributions from mine-affected tributaries have an additive effect on the toxicity of both water (4A, April 2011n=37-251) and sediment (4B, December 2010 n=148-467) from the main branch of the Mud River. Statistically significant difference between upstream and downstream site indicated by asterisks. “Size” refers to the optical density (extinction) of each nematode and describes growth since the data shown is for size after three days of exposure. Each exposure was repeated at least two times separately. Please note the difference in Y-axis scaling.

https://doi.org/10.1371/journal.pone.0075329.g004

IV: MTR/VF samples did not cause induction of mtl-2::GFP

Sediment pore water from the Mud River and two tributaries, sampled in December 2010 and April 2011, did not cause an increase in metallothionein expression in mtl-2::GFP transgenic nematodes (Figure S1). A statistically significant increase in GFP was induced by cadmium chloride treated positive controls (p= 0.0001). While manganese and selenium are not commonly reported to be metallothionein inducers (we found one report of induction by selenite [40]:), we also tested manganese chloride and sodium selenite because they are suspected contaminants in this system. However, we observed no induction in either case (Figure S1).

V: Genetic approaches reveal different mechanisms of toxicity for water and sediment samples

To further investigate both the potential mechanisms of toxicity responsible for both the growth inhibition seen in wild-type nematodes and the differences observed between stream water and sediment pore water exposures, we employed a number of stress-response knockout mutants which we hypothesized would be more sensitive to potential toxicants and thus display greater growth inhibition than wild type worms. We examined strains lacking genes important in protecting this organism against metal toxicity, osmotic stress, and oxidative stress, three potential mechanisms of toxicity associated with MTR/VF samples.

Metal detoxification in nematodes is mediated by a number of different metals-sequestering peptides, including two metallothioneins, mtl-1 and mtl-2, the latter of which is strongly induced in the gut following exposure to cadmium and other heavy metals [41]. Although mtl-2 deletion mutants have often been employed as metals-sensitive biosensors using a variety of different endpoints [19,42,43], growth of mtl-2(gk125) mutants was not significantly different from wild-type when exposed to water and sediment pore water from Mud River 6, 7, 9, 10 and 11 (data not shown). in keeping with the lack of metallothionein induction seen in the sediment-exposed mtl-2::GFP nematodes.

Phytochelatin, a cysteine-rich metal chelating peptide best studied in yeast and plants, also occurs in nematodes. Production of phytochelatin, using glutathione as a precursor, is driven by phytochelatin synthase. A deletion mutant strain lacking the phytochelatin synthase gene (pcs-1(tm1748)) is more sensitive to a number of metals than either the wild-type strain or metallothionein deletion mutants [42], and is also more sensitive to selenium (Figure 3). Growth inhibition in water-exposed pcs-1(tm1748) mutants was comparable to that of wild-type nematodes (Figure 5). However, sediment pore water from multiple river and tributary sites induced greater growth inhibition in pcs-1(tm1748) mutants than their wild-type counterparts (Figure 5, p = 0.0022 or less), indicating the presence of biologically relevant amounts of metals or metalloids in stream sediment. Lethality was not observed for any of the exposure conditions.

[Figure omitted. See PDF.]

Figure 5. Toxicity of sediment pore water is driven by metals or metalloids.

Stream sediment pore water, but not stream water, caused greater growth inhibition in metals-sensitive pcs-1(tm1748) and smf-2 mutants than in wild-type nematodes, demonstrated in samples collected from Mud River sites 5a and 7 in December 2010 (5A p=<0.0001, n=142-533, 5B p=0.0022, n= 42-277), Mud River 5a in April 2011 (5C p=<0.0001 n=55-277) and Berry Branch in April 2011 (5D p=<0.0001, n= 125-482). Statistically significant difference from sediment-treated wild-type nematodes indicated by asterisks. “Size” refers to the optical density (extinction) of each nematode and describes growth since the data shown is for size after three days of exposure. Each exposure was repeated at least two times separately.

https://doi.org/10.1371/journal.pone.0075329.g005

Manganese is another potential contaminant of interest in MTR/VF impacted rivers and streams [11,44]. Uptake and transport of manganese and iron in nematodes is regulated by the NRAMP-type divalent cation transporters smf-1, smf-2 and smf-3. While smf-1 and smf-3 deletion mutants are deficient in the uptake of manganese and display greater tolerance, smf-2(gk133) deletion mutants are more sensitive to manganese exposure than wild-type worms [11,45]. Growth inhibition in smf-2(gk133) knockout mutants was similar to wild-type for most sites. Smf-2(gk133) knockout mutants were exposed to sediment pore water from Mud River 5-10, Stanley Fork, Laurel Branch and Berry Branch from December 2010, water and sediment pore water from Mud River 5-8b, 10, Stanley Fork, Laurel Branch, Berry Branch and water only from Connelly Branch from April 2011. Lethality was not observed for any of the exposure conditions. For site samples that elicited greater growth inhibition in smf-2(gk133) mutants versus wild-type worms, a similar effect was seen in pcs-1 mutants (Figure 5b). There were no sites that produced growth-inhibited smf-2(gk133) worms without a corresponding effect in pcs-1 mutants.

Response to osmotic stress in nematodes is governed by two glycerol 3-phosphate dehydrogenases, gpdh-1 and gpdh-2, which regulate synthesis and accumulation of glycerol in response to increased osmolarity. gpdh-1 is strongly up-regulated in response to osmotic stress, whereas gpdh-2 is constitutively expressed in several tissues and weakly up-regulated in response to osmotic stress. Mutant nematodes lacking either of these genes are slightly more sensitive to osmotic stress, and double knockouts are markedly more so, exhibiting growth inhibition in solutions with higher solute concentration [46]. We found that gpdh-1(ok1558);gpdh-2(kb33) mutants exhibited greater growth inhibition than wild-type nematodes following exposure to stream water, but not sediment pore water, from two sites with elevated conductivity(Figure 6) in comparison to reference site. Stream water from Mud River 6 and Mud River 7 taken in the same sampling trip did not cause strain-specific growth inhibition in keeping with the fact that the samples were taken shortly after flooding and had lower conductivities (growth data not shown, for conductivity data see Table S2).

[Figure omitted. See PDF.]

Figure 6. Toxicity of stream water is driven by osmotic stress.

Stream water, but not sediment pore water, from the Mud River and its tributaries caused greater growth inhibition in osmotic stress-sensitive gpdh-1(ok1558);gpdh-2(kb33) deletion mutants in comparison to wild-type nematodes, demonstrated in samples collected from Sugartree Branch (6A and 6B p=0.0023n=23-55) and Laurel Branch (6C and 6D p=0.0004n=39-51), May 2012. Statistically significant difference from reference site indicated by asterisks. “Size” refers to the optical density (extinction) of each nematode and describes growth since the data shown is for size after three days of exposure. Each exposure was repeated at least two times separately.

https://doi.org/10.1371/journal.pone.0075329.g006

Lethality was not observed for any of the exposure conditions.

Lastly, we used superoxide dismutase deficient nematodes to assess the potential role of oxidative stress in the toxicity of Mud River samples. sod-3(gk257) mutants, deficient in a mitochondria-localized iron/manganese superoxide dismutase and sensitive to oxidative stress [47,48], did not display greater growth inhibition than their wild-type counterparts when exposed to either water or sediment pore water from Mud River 6, 7, 9, 10 and 11 collected in June 2010 (data not shown).

Discussion

The Mud River watershed is a complex and dynamic system. Growth data from wild-type nematodes exposed to samples taken from 3 different sampling trips (December 2010, April 2011, May 2012) reflects variability of water and sediment chemistry dependent on season, water level and influx of sulfates, metals, metalloids and assorted inorganic contaminants from the Hobet 21 mine. Although nematodes are more robust in response to many stressors than the native stream fauna that are typically first affected by mine runoff in mine-impacted ecosystems [49], wild-type nematodes were growth-inhibited upon exposure to both water and sediment pore water from the mining-impacted portions of the Mud River and its tributaries. As has been reported in previous research [50], impaired tributaries contributed to toxicity of water in the main stem of the river.

However, the major advantage of using C. elegans in these studies is the toolkit available for investigating mechanisms of toxicity. Nematodes employ a wide variety of stress-response mechanisms in order to maintain health in the face of adversity; reporter strains and knockout mutants allow the relationships between sensitivity to specific exposures and the mechanism of toxicity to be examined.

Using mtl-2::GFP expression as an endpoint, we detected no statistically significant increase in expression in comparison to either EPA water controls or nematodes treated with sediment pore water from our reference site. This might be explained by the absence of the types of metals that result in metallothionein induction. Metallothionein expression is induced most effectively by cadmium and other heavy metals, which are not typically present at elevated levels in mine-impacted Appalachian watersheds [11,14]. It is also possible that the lack of inductions reflect insufficient sensitivity of the reporter strain. However, mtl-2(gk125) worms grew comparably to wild-type in the presence of mining-impacted streamwater and sediment pore water. While mtl-2(gk125) deletion mutants are commonly employed as a metals-sensitive strain, recent research indicates that phytochelatin is more important for the mitigation of metal toxicity [22,42,51]. Our comparisons of the growth of mtl-2(gk125) and pcs-1 (tm1748) deletion mutants support this data, since unlike the mtl-2 (gk125) mutants, pcs-1(tm1748) mutants exhibited greater growth inhibition than wild-type nematodes when exposed to sediment from the Mud River (Figure S2).

While a genetic approach by itself cannot resolve which metals or metalloids are involved, our results do provide indirect evidence for an involvement of selenium. We found that pcs-1(tm1748) knockouts were more sensitive than wild-type to selenium, which is a known contaminant of the Mud River. While the relationship between phytochelatin and selenium remains unstudied in nematodes, algae [52] produce phytochelatin upon exposure to selenium, and selenium forms trisulfide linkages with phytochelatin in vitro [38]. The greater sensitivity of pcs-1(tm1748) knockout nematodes to sodium selenite in comparison to wild-type worms and a corresponding lack of induction of metallothionein expression after exposure to selenium, in concert with research by other groups demonstrating the greater importance of phytochelatin in comparison to metallothionein in the elimination of metals, suggest that phytochelatin is responsible for mitigating the effects of excess selenium in nematodes. A limitation of our experimental design is that most environmental selenium uptake is believed to occur via ingestion, rather than direct uptake from the medium. Although this is not something that we could test given the samples available to us, we note that the sediment pore water we used was not filtered and thus likely included ingestible particles. Future experiments using biofilms will be important to further our understanding of the potential role of ingestion of selenium.

Just as the impaired growth of pcs-1(tm1748) knockout nematodes in sediment pore water suggests the presence of biologically relevant levels of one or more metals or metalloids, the impaired growth of osmotic stress sensitive gpdh-1(ok1558);gpdh-2(kb33) knockouts implicates osmotic stress as another mechanism of toxicity. Previous research has demonstrated elevated levels of TDS in MTR/VF impacted Appalachian watersheds, including the Mud River [3,53]. Conductivity values exceeding the water quality benchmark of 300µs/cm-1 have been observed in the water of both the main stem and mining-affected tributaries (Table S2). Also relevant is the potential for mixture effects; mine output and the release of minerals via weathering are highly variable, and synergistic toxicity may be a factor [54-57].

Our interrogation via mutant studies of the differential growth inhibition between water and sediment pore water seen in wild-type nematodes reveals that the relevant mechanism of toxicity in the water column is primarily osmotic stress due to elevated TDS whereas metals or metalloids, most likely including selenium, are the most important toxicant seen in the sediment pore water. Positive growth effects caused by sediment pore water from the Left Fork are likely explained by sub-toxic quantities of selenium; while the Left Fork does not receive mine discharge, selenium in potentially biologically relevant quantities has been observed in biofilms collected from the Left Fork (Mariah Arnold, unpublished observation). Selenium in the mined portion of this system has been reported as typically ~1-20 µg/L in water [3]. However, ~90% of the selenium in the water column is present as selenate (Helen Hsu-Kim and Lusia Liu, personal communication) and is therefore unlikely to account for the growth inhibition that we observed, since selenate up to 150 µg/L had no effect on the growth of C. elegans in our experiments. On the other hand, sediment selenium was typically present as multiple species in this system, including selenite (Helen Hsu-Kim and Lusia Liu, personal communication), and might therefore contribute to growth inhibition as seen in our selenite growth assays (Figure 3).

We extend the previous successful use of C. elegans to investigate the environmental impact of metals mining [12] to the highly complex alkaline mine runoff associated with coal mining, illustrating the utility of C. elegans as a model organism for ecotoxicological research using environmental samples. Our work, in combination with chemical analyses and surveys of stream biota [3,6,58,59] supports the presence of biologically relevant levels of multiple toxicants in the mined portion of the Mud River system.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Sediment pore water did not induce mtl-2::GFP. Neither sediment pore water from Mud River and tributaries nor manganese and selenium controls caused a statistically significant increase in GFP expression in mtl-2::GFP transgenic nematodes in comparison to EPA water and Left Fork controls (n=18-147). Photos show uninduced mtl-2::GFP nematode (left) and cadmium-exposed mtl-2::GFP nematode (right).

https://doi.org/10.1371/journal.pone.0075329.s001

(TIF)

Figure S2.

Sediment pore water caused greater growth inhibition in pcs-1 knockout nematodes vs. mtl-2 knockouts. Sediment pore water from Mud River 6 (2A) and Mud River 7 (2B) caused greater growth inhibition in pcs-1(tm1748) deletion mutants (2A p=<0.0001 n=47-406, 2B p=<0.0001 n=27-406) in comparison to wild-type nematodes, but not mtl-2(gk125) deletion mutants (2A p=0.3609 n=47-406, 2B p=0.9487 n=27-406). Statistically significant difference from sediment-treated wild-type nematodes indicated by asterisks. “Size” refers to the optical density (extinction) of each nematode and describes growth since the data shown is for size after three days of exposure. Each exposure was repeated at least two times separately.

https://doi.org/10.1371/journal.pone.0075329.s002

(TIF)

Table S1.

p values, Mann Whitney U test in comparison to EPA water controls.

https://doi.org/10.1371/journal.pone.0075329.s003

(DOCX)

Table S2.

Conductivity values, Mud River and tributaries.

https://doi.org/10.1371/journal.pone.0075329.s004

(DOCX)

Acknowledgments

We thank Ty Lindberg for sample collection assistance, and Zhirui Zhu for outcrossing the sod-3 strain.

Author Contributions

Conceived and designed the experiments: EAT JNM. Performed the experiments: EAT BLT. Analyzed the data: EAT. Contributed reagents/materials/analysis tools: EAT GLK MCA RDG. Wrote the manuscript: EAT.

Citation: Turner EA, Kroeger GL, Arnold MC, Thornton BL, Di Giulio RT, Meyer JN (2013) Assessing Different Mechanisms of Toxicity in Mountaintop Removal/Valley Fill Coal Mining-Affected Watershed Samples Using Caenorhabditis elegans. PLoS ONE 8(9): e75329. https://doi.org/10.1371/journal.pone.0075329

References

1. USEPA (2011) The Effects of Mountaintop Mines and Valley Fills on Aquatic Ecosystems of the Central Appalachian Coalfields. Washington DC: US Environmental Protection Agency. EPA/600/R-609/138F p.

2. Luoma SN, Presser TS (2009) Emerging opportunities in management of selenium contamination. Environ Sci Technol 43: 8483-8487. doi:https://doi.org/10.1021/es900828h. PubMed: 20028041.

3. Lindberg TT, Bernhardt ES, Bier R, Helton AM, Merola RB et al. (2011) Cumulative impacts of mountaintop mining on an Appalachian watershed. Proc Natl Acad Sci U S A 108: 20929-20934. doi:https://doi.org/10.1073/pnas.1112381108. PubMed: 22160676.

4. Chapman PM (1986) Sediment quality criteria from the sediment quality triad- an example. Environ Toxicol Chem 5: 957-964. doi:https://doi.org/10.1002/etc.5620051104.

5. Bernhardt ES, Lutz BD, King RS, Fay JP, Carter CE et al. (2012) How many mountains can we mine? Assessing the regional degradation of central appalachian rivers by surface coal mining. Environ Sci Technol 46: 8115-8122. doi:https://doi.org/10.1021/es301144q. PubMed: 22788537.

6. Pond GJ, Passmore ME, Borsuk FA, Reynolds L, Rose CJ (2008) Downstream effects of mountaintop coal mining: comparing biological conditions using family- and genus-level macroinvertebrate bioassessment tools. J North American Benthological Society 27: 717-737. doi:https://doi.org/10.1899/08-015.1.

7. Buchwalter DB, Cain DJ, Martin CA, Xie L, Luoma SN et al. (2008) Aquatic insect ecophysiological traits reveal phylogenetically based differences in dissolved cadmium susceptibility. Proc Natl Acad Sci U S A 105: 8321-8326. doi:https://doi.org/10.1073/pnas.0801686105. PubMed: 18559853.

8. Hitt NP, Hendryx M (2010) Ecological integrity of streams related to human cancer mortality rates. Ecohealth 7: 91-104. doi:https://doi.org/10.1007/s10393-010-0297-y. PubMed: 20361230.

9. Ghose MK, Majee SR (2007) Characteristics of hazardous airborne dust around an Indian surface coal mining area. Environ Monit Assess 130: 17-25. doi:https://doi.org/10.1007/s10661-006-9448-6. PubMed: 17285255.

10. Palmer MA, Bernhardt ES, Schlesinger WH, Eshleman KN, Foufoula-Georgiou E et al. (2010) Science and regulation. Mountaintop mining consequences. Science 327: 148-149. doi:https://doi.org/10.1126/science.1180543. PubMed: 20056876.

11. Griffith MB, Norton SB, Alexander LC, Pollard AI, LeDuc SD (2012) The effects of mountaintop mines and valley fills on the physicochemical quality of stream ecosystems in the central Appalachians: A review. Sci Total Environ 417: 1-12. PubMed: 22264919.

12. Anbalagan C, Lafayette I, Antoniou-Kourounioti M, Haque M, King J et al. (2012) Transgenic nematodes as biosensors for metal stress in soil pore water samples. Ecotoxicology 21: 439-455. doi:https://doi.org/10.1007/s10646-011-0804-0. PubMed: 22037694.

13. Leung MC, Williams PL, Benedetto A, Au C, Helmcke KJ et al. (2008) Caenorhabditis elegans: an emerging model in biomedical and environmental toxicology. Toxicol Sci 106: 5-28. doi:https://doi.org/10.1093/toxsci/kfn121. PubMed: 18566021.

14. Williams PL, Dusenbery DB (1990) Aquatic toxicity testing using the nematode, Caenorhabditis-elegans. Environ Toxicol Chem 9: 1285-1290. doi:https://doi.org/10.1897/1552-8618(1990)9[1285:ATTUTN]2.0.CO;2.

15. Mutwakil MH, Reader JP, Holdich DM, Smithurst PR, Candido EPM et al. (1997) Use of stress-inducible transgenic nematodes as biomarkers of heavy metal pollution in water samples from an english river system. Arch Environ Contam Toxicol 32: 146-153. doi:https://doi.org/10.1007/s002449900167. PubMed: 9069189.

16. Félix M-A, Braendle C (2010) The natural history of Caenorhabditis elegans. Curr Biol 20: R965-R969. doi:https://doi.org/10.1016/j.cub.2010.09.050. PubMed: 21093785.

17. Stiernagle T (1999) Maintenance of C. elegans. In: IA Hope. C elegans A Practical Approach. Oxford University Press. pp. 59-60.

18. Traunspurger W, Haitzer M, Höss S, Beier S, Ahlf W et al. (1997) Ecotoxicological assessment of aquatic sediments with Caenorhabditis elegans (nematoda) — a method for testing liquid medium and whole-sediment samples. Environ Toxicol Chem 16: 245-250. doi:https://doi.org/10.1897/1551-5028(1997)016.

19. Meyer JN, Lord CA, Yang XY, Turner EA, Badireddy AR et al. (2010) Intracellular uptake and associated toxicity of silver nanoparticles in Caenorhabditis elegans. Aquat Toxicol 100: 140-150. doi:https://doi.org/10.1016/j.aquatox.2010.07.016. PubMed: 20708279.

20. Koga M, Zwaal R, Guan KL, Avery L, Ohshima Y (2000) A Caenorhabditis elegans MAP kinase kinase, MEK-1, is involved in stress responses. EMBO J 19: 5148-5156. doi:https://doi.org/10.1093/emboj/19.19.5148. PubMed: 11013217.

21. Gartner A, Milstein S, Ahmed S, Hodgkin J, Hengartner MO (2000) A Conserved Checkpoint Pathway Mediates DNA Damage–Induced Apoptosis and Cell Cycle Arrest in C. elegans. Mol Cell 5: 435-443. doi:https://doi.org/10.1016/S1097-2765(00)80438-4. PubMed: 10882129.

22. Yang XY, Gondikas AP, Marinakos SM, Auffan M, Liu J et al. (2012) Mechanism of Silver Nanoparticle Toxicity Is Dependent on Dissolved Silver and Surface Coating in Caenorhabditis elegans. Environ Sci Technol 46: 1119-1127. doi:https://doi.org/10.1021/es202417t. PubMed: 22148238.

23. Martinez-Finley EJ, Aschner M (2011) Revelations from the Nematode Caenorhabditis elegans on the Complex Interplay of Metal Toxicological Mechanisms. J of Toxicology, 2011: 2011: 895236. PubMed: 21876692.

24. Ishii N, Takahashi K, Tomita S, Keino T, Honda S et al. (1990) A methyl viologen-sensitive mutant of the nematode Caenorhabditis-elegans. Mutat Res 237: 165-171. doi:https://doi.org/10.1016/0921-8734(90)90022-J. PubMed: 2233820.

25. North M, Vulpe CD (2010) Functional toxicogenomics: mechanism-centered toxicology. Int J Mol Sci 11: 4796-4813. doi:https://doi.org/10.3390/ijms11124796. PubMed: 21614174.

26. Ma H, Glenn TC, Jagoe CH, Jones KL, Williams PL (2009) A transgenic strain of the nematode Caenorhabditis elegans as a biomonitor for heavy metal contamination. Environ Toxicol Chem 28: 1311-1318. doi:https://doi.org/10.1897/08-496.1. PubMed: 19175297.

27. Lewis JA, Fleming JT (1995) Basic culture methods. Methods Cell Biol 48: 3-29. doi:https://doi.org/10.1016/S0091-679X(08)61381-3. PubMed: 8531730.

28. Gruber J, Ng LF, Poovathingal SK, Halliwell B (2009) Deceptively simple but simply deceptive - Caenorhabditis elegans lifespan studies: Considerations for aging and antioxidant effects. FEBS Lett 583: 3377-3387. doi:https://doi.org/10.1016/j.febslet.2009.09.051. PubMed: 19815017.

29. Smith ME, Lazorchak JM, Herrin LE, BrewerSwartz S, Thoeny WT (1997) A reformulated, reconstituted water for testing the freshwater amphipod, Hyalella azteca. Environ Toxicol Chem 16: 1229-1233. doi:https://doi.org/10.1002/etc.5620160618.

30. Soucek DJ, Kennedy AJ (2005) Effects of hardness, chloride, and acclimation on the acute toxicity of sulfate to freshwater invertebrates. Environ Toxicol Chem 24: 1204-1210. doi:https://doi.org/10.1897/04-142.1. PubMed: 16111001.

31. Borgmann U (1996) Systematic analysis of aqueous ion requirements of Hyalella azteca: A standard artificial medium including the essential bromide ion. Arch Environ Contam Toxicol 30: 356-363. doi:https://doi.org/10.1007/BF00212294.

32. Li WH, Hsu FL, Liu JT, Liao VH (2011) The ameliorative and toxic effects of selenite on Caenorhabditis elegans. Food Chem Toxicol 49: 812-819. doi:https://doi.org/10.1016/j.fct.2010.12.002. PubMed: 21145367.

33. Debruyn AMH, Chapman PM (2007) Selenium toxicity to invertebrates: Will proposed thresholds for toxicity to fish and birds also protect their prey? Environ Sci Technol 41: 1766-1770. doi:https://doi.org/10.1021/es062253j. PubMed: 17396671.

34. Conley JM, Funk DH, Buchwalter DB (2009) Selenium bioaccumulation and maternal transfer in the mayfly Centroptilum triangulifer in a life-cycle, periphyton-biofilm trophic assay. Environ Sci Technol 43: 7952-7957. doi:https://doi.org/10.1021/es9016377. PubMed: 19921919.

35. Conley JM, Funk DH, Cariello NJ, Buchwalter DB (2011) Food rationing affects dietary selenium bioaccumulation and life cycle performance in the mayfly Centroptilum triangulifer. Ecotoxicology 20: 1840-1851. doi:https://doi.org/10.1007/s10646-011-0722-1. PubMed: 21701843.

36. Estevez AO, Mueller CL, Morgan KL, Szewczyk NJ, Teece L et al. (2012) Selenium induces cholinergic motor neuron degeneration in Caenorhabditis elegans. Neurotoxicology 33: 1021-1032. doi:https://doi.org/10.1016/j.neuro.2012.04.019. PubMed: 22560997.

37. Eisler R (2000) Selenium. Handbook of Chemical Risk Assessment. Lewis Publishers. pp. 1649-1693.

38. Spain SM, Rabenstein DL (2004) Characterization of the selenotrisulfide formed by reaction of selenite with end-capped phytochelatin-2. Anal Bioanal Chem 378: 1561-1567. doi:https://doi.org/10.1007/s00216-003-2415-y. PubMed: 15214417.

39. Simmons DBWallschläger D (2005) A critical review of the biogeochemistry and ecotoxicology of selenium in lotic and lentic environments. Environ Toxicol Chem 24: 1331-1343. doi:https://doi.org/10.1897/04-176R.1. PubMed: 16117108.

40. Iwai N, Watanabe C, Suzuki T, Suzuki KT, Tohyama C (1988) Metallothionein induction by sodium selenite at two different ambient temperatures in mice. Arch Toxicol 62: 447-451. doi:https://doi.org/10.1007/BF00288348. PubMed: 3255298.

41. Freedman JH, Slice LW, Dixon D, Fire A, Rubin CS (1993) The novel metallothionein genes of Caenorhabditis elegans. Structural organization and inducible, cell-specific expression. J Biol Chem 268: 2554-2564. PubMed: 8428932.

42. Hughes SL, Bundy JG, Want EJ, Kille P, Stürzenbaum SR (2009) The metabolomic responses of Caenorhabditis elegans to cadmium are largely independent of metallothionein status, but dominated by changes in cystathionine and phytochelatins. J Proteome Res 8: 3512-3519. doi:https://doi.org/10.1021/pr9001806. PubMed: 19466807.

43. Swain SC, Keusekotten K, Baumeister R, Stürzenbaum SR (2004) C. elegans metallothioneins: new insights into the phenotypic effects of cadmium toxicosis. J Mol Biol 341: 951-959. doi:https://doi.org/10.1016/j.jmb.2004.06.050. PubMed: 15328611.

44. Dittman EK, Buchwalter DB (2010) Manganese Bioconcentration in Aquatic Insects: Mn Oxide Coatings, Molting Loss, and Mn(II) Thiol Scavenging. Environ Sci Technol 44: 9182-9188. doi:https://doi.org/10.1021/es1022043. PubMed: 21049994.

45. Au C, Benedetto A, Anderson J, Labrousse A, Erikson K et al. (2009) SMF-1, SMF-2 and SMF-3 DMT1 orthologues regulate and are regulated differentially by manganese levels in C. elegans. PLOS ONE 4: e7792. doi:https://doi.org/10.1371/journal.pone.0007792. PubMed: 19924247.

46. Lamitina T, Huang CG, Strange K (2006) Genome-wide RNAi screening identifies protein damage as a regulator of osmoprotective gene expression. Proc Natl Acad Sci U S A 103: 12173-12178. doi:https://doi.org/10.1073/pnas.0602987103. PubMed: 16880390.

47. Dingley S, Polyak E, Lightfoot R, Ostrovsky J, Rao M et al. (2010) Mitochondrial respiratory chain dysfunction variably increases oxidant stress in Caenorhabditis elegans. Mitochondrion 10: 125-136. doi:https://doi.org/10.1016/j.mito.2009.11.003. PubMed: 19900588.

48. Arnold MC, Badireddy AR, Wiesner MR, Di Giulio RT, Meyer JN (2013) Cerium oxide nanoparticles are more toxic than equimolar bulk cerium oxide in Caenorhabditis elegans. Arch Environ Contam Toxicol, 65: 224–33. PubMed: 23619766.

49. Malaj E, Grote M, Schäfer RB, Brack W, von der Ohe PC (2012) Physiological sensitivity of freshwater macroinvertebrates to heavy metals. Environ Toxicol Chem 31: 1754-1764. doi:https://doi.org/10.1002/etc.1868. PubMed: 22553143.

50. Alexander RB, Boyer EW, Smith RA, Schwarz GE, Moore RB (2007) The Role of Headwater Streams in Downstream Water Quality. J Am Water Resour Assoc 43: 41-59. doi:https://doi.org/10.1111/j.1752-1688.2007.00005.x. PubMed: 22457565.

51. Vatamaniuk OK, Bucher EA, Ward JT, Rea PA (2001) A new pathway for heavy metal detoxification in animals. Phytochelatin synthase is required for cadmium tolerance in Caenorhabditis elegans. J Biol Chem 276: 20817-20820. doi:https://doi.org/10.1074/jbc.C100152200. PubMed: 11313333.

52. Simmons DBEmery RJ (2011) Phytochelatin induction by selenate in Chlorella vulgaris, and regulation of effect by sulfate levels. Environ Toxicol Chem 30: 469-476. doi:https://doi.org/10.1002/etc.392. PubMed: 21038439.

53. Cormier SM, Wilkes SP, Zheng L (2013) Relationship of land use and elevated ionic strength in Appalachian watersheds. Environ Toxicol Chem 32: 296-303. doi:https://doi.org/10.1002/etc.2055. PubMed: 23161531.

54. Power RS, de Pomerai DI (1999) Effect of single and paired metal inputs in soil on a stress-inducible transgenic nematode. Arch Environ Contam Toxicol 37: 503-511. doi:https://doi.org/10.1007/s002449900545. PubMed: 10508898.

55. Dennis JL, Mutwakil MHAZ, Lowe KC, de Pomerai DI (1997) Effects of metal ions in combination with a non-ionic surfactant on stress responses in a transgenic nematode. Aquat Toxicol 40: 37-50. doi:https://doi.org/10.1016/S0166-445X(97)00043-X.

56. Ganther HE, Baumann CA (1962) Selenium metabolism. II. Modifying effects of sulfate. J Nutr 77: 408–414&. PubMed: 13896696.

57. Sasakura C, Suzuki KT (1998) Biological interaction between transition metals (Ag, Cd and Hg), selenide/sulfide and selenoprotein P. J Inorg Biochem 71: 159-162. doi:https://doi.org/10.1016/S0162-0134(98)10048-X. PubMed: 9833321.

58. Hartman KJ, Kaller MD, Howell JW, Sweka JA (2005) How much do valley fills influence headwater streams? Hydrobiologia 532: 91-102. doi:https://doi.org/10.1007/s10750-004-9019-1.

59. Cormier SM, Suter GW 2nd, Zheng L, Pond GJ (2013) Assessing causation of the extirpation of stream macroinvertebrates by a mixture of ions. Environ Toxicol Chem 32: 277-287. doi:https://doi.org/10.1002/etc.2059. PubMed: 23147750.