Evolutionary Toxicology-An Informational Tool for Chemical Regulation?

Pollution offsets the rapid evolution of increased heat tolerance in a natural population

Despite the increasing evidence for rapid thermal evolution in natural populations, evolutionary rescue under global warming may be constrained by the presence of other stressors. Highly relevant in our polluted planet, is the largely ignored evolutionary trade-off between heat tolerance and tolerance to pollutants. By using two subpopulations (separated 40 years in time) from a resurrected natural population of the water flea Daphnia magna that experienced a threefold increase in heat wave frequency during this period, we tested whether rapid evolution of heat tolerance resulted in reduced tolerance to the widespread metal zinc and whether this would affect heat tolerance upon exposure to the pollutant. Our results revealed rapid evolution of increased heat tolerance in the recent subpopulation. Notably, the sensitivity to the metal tended to be stronger (reduction in net energy budget) or was only present (reductions in heat tolerance and in sugar content) in the recent subpopulation. As a result, the rapidly evolved higher heat tolerance of the recent subpopulation was fully offset when exposed to zinc. Our results highlight that the many reports of evolutionary rescue to global change stressors may give a too optimistic view as our warming planet is polluted by metals and other pollutants.

Chronic metal contamination shapes the size structure of Gammarus fossarum populations in French headwater rivers

Assessing the effects of multigenerational exposure of aquatic animal populations to chemical contamination is essential for ecological risk assessment. However, beyond rare examples reporting the sporadic emergence of a toxicological tolerance within populations that persist in contaminated environments, conclusive results are even more limited from field studies when it comes to the alteration of life-history traits. Here, we investigated whether long-term exposure to cadmium (Cd) influences size-related life-history traits (i.e., size at puberty, median adult size, maximum size) in Gammarus fossarum, a keystone species of European stream ecosystems. We studied 13 field populations of G. fossarum (cryptic lineage B) living in headwater rivers located in natural areas scattered at a large geographical scale and exposed to contrasted bioavailable Cd contamination levels due to different local geochemical backgrounds. We achieved a detailed description of the physical and physicochemical conditions of the river reaches investigated. Land-use parameters, hydrological characteristics (flow, slope, river width, flow structure, mosaic of substrates), and physicochemical conditions (temperature, conductivity, dissolved oxygen) were measured. Metallic bioavailable contamination was assessed using a standardized active biomonitoring procedure (Gammarus caging). Based on the field demographic census of the 13 populations, our results demonstrated that chronic Cd contamination significantly influences life-history in the G. fossarum species, with a significant reduction in all size traits of populations (size at puberty, median adult size, maximum size). In addition, we confirmed Cd-tolerance in contaminated populations during exposure tests in the laboratory. Various hypotheses can be then put forward to explain the modification of size-related life-history traits: a direct toxic effect of Cd, a cost of Cd-tolerance, or an adaptive evolution of life-history exposed to toxic pressure.

Contamination acts as a genotype-dependent barrier to gene flow, causing genetic erosion and fine-grained population subdivision in Mussels from the Strait of Istanbul

This study provides evidence of fine-grained genetic structuring in Mediterranean mussels (Mytilus galloprovincialis) from the Strait of Istanbul, caused by barriers to gene flow via contaminant-mediated selection. In this study, mitochondrial D-loop sequences were analyzed in mussels from 8 localities, all less than 30 kilometers apart, with differing contaminant loads. The results were: 1) Intra-population genetic differentiation (ΦST) between sites with high and low contaminant loads was high (up to 0.459), even at distances of only a few kilometers. 2) Genetic diversity was negatively correlated with the contaminant load ("genetic erosion"). 3) There was evidence of selection, based on haplotype frequencies and neutrality tests (Tajima's D), with purifying selection at the most contaminated site and balancing selection at the least contaminated. 4) Genetic distance was not correlated with geographic distance (no isolation-by-distance), but was correlated with contaminant load at each site. 5) Population dendrograms and Bayesian estimators of migration indicated that gene flow between sites was affected by contamination. For the dendrograms of the sampling sites, the clades clustered according to contaminant load more than geographic distance. Overall, these results suggest that 1) contamination may serve as a genotype-dependent dispersal barrier (i.e., contamination may not affect total number of migrants, just the relative proportions of the haplotypes in the established immigrants), leading strong population differentiation over short distances, and 2) genetic erosion may occur by a combination of selection and altered patterns of haplotype-specific gene flow. These effects may be more pronounced in the Strait of Istanbul than in other locations because of the riverine nature and strong, uni-directional current of the strait.

Intraspecific genetic variation is critical to robust toxicological predictions of aquatic contaminants

Environmental risk assessment is a critical tool for protecting aquatic life and its effectiveness is predicated on predicting how natural populations respond to contaminants. Yet, routine toxicity testing typically examines only one genotype, which may render risk assessments inaccurate as populations are most often composed of genetically distinct individuals. To determine the importance of intraspecific variation in the translation of toxicity testing to populations, we quantified the magnitude of genetic variation within 20 Daphnia magna clones derived from one lake using whole genome sequencing and phenotypic assays. We repeated these assays across two exposure levels of microcystins, a cosmopolitan and lethal aquatic contaminant produced by harmful algal blooms. We found considerable intraspecific genetic variation in survival, growth, and reproduction, which was amplified by microcystins exposure. Finally, using simulations we demonstrate that the common practice of employing a single genotype to calculate toxicity tolerance failed to produce an estimate within the 95% confidence interval over half of the time. These results illuminate the importance of incorporating intraspecific genetic variation into toxicity testing to reliably predict how natural populations will respond to aquatic contaminants.

Variation of Tolerance to Isothiazolinones Among Daphnia pulex Clones

Isothiazolinones are a family of broad-spectrum biocides widely used in industry and consumer products. Chloro- and methyl-isothiazolinones (CMIT and MIT) are documented as strong irritants, yet they are still used in a wide variety of applications, including cosmetics, cleansers, hygienic products, and various industrial applications. The subsequent substantial release of these molecules from urban sources into freshwater environments, and their potential impacts on aquatic species, have nevertheless received little attention so far, with few studies reporting on the toxicity of either CMIT or MIT to nontarget organisms. The present study addresses this current knowledge gap by evaluating the acute toxicity to Daphnia pulex (Cladocera) of CMIT/MIT (3:1) and MIT, the two formulations most commonly used by manufacturers. In addition, genetic diversity is known to be a major component of variability in phenotypic responses, although it is largely overlooked in typical toxicity tests. Thus the potential range of responses inherent to genetic diversity is rarely considered. Therefore, to account for intraspecific variations in sensitivity, our design involved eight clonal lines of D. pulex stemming from distinct natural populations or commercial strains. Clones exhibited strong variation in their responses, with median lethal concentration (LC50) values ranging from 0.10 to 1.84 mg/L for the mixture CMIT/MIT, and from 0.68 to 2.84 mg/L for MIT alone. These intraspecific ranges of LC50 values challenge the use of single clones of daphnids in standard ecotoxicological tests and the predictions based on their results. The present study brings new evidence that assessing ecological risk of chemicals while ignoring genotype diversity is neither ecologically relevant, nor a representative evaluation of the diversity of potential adverse outcomes. Environ Toxicol Chem 2023;42:805-814. © 2023 SETAC.

Predicting Resistance: Quantifying the Relationship between Urban Development, Agricultural Pesticide Use, and Pesticide Resistance in a Nontarget Amphipod

Resistance alleles within the voltage-gated sodium channel (vgsc) have been correlated with pyrethroid resistance in wild populations of the nontarget amphipod, Hyalella azteca from California (CA), U.S.A. In the present study, we expand upon the relationship between land use and the evolution of pesticide resistance in H. azteca to develop a quantitative methodology to target and screen novel populations for resistance allele genotypes in a previously uninvestigated region of the U.S. (New England: NE). By incorporating urban land development and toxicity-normalized agricultural pesticide use indices into our site selection, we successfully identified three amino acid substitutions associated with pyrethroid resistance. One of the resistance mutations has been described in H. azteca from CA (L925I). We present the remaining two (vgsc I936F and I936V) as novel pyrethroid-resistance alleles in H. azteca based on previous work in insects and elevated cyfluthrin resistance in one NE population. Our results suggest that urban pesticide use is a strong driver in the evolution of resistance alleles in H. azteca. Furthermore, our method for resistance allele screening provides an applied framework for detecting ecosystem impairment on a nationwide scale that can be incorporated into ecological risk assessment decisions.

Enhanced trophic transfer of chlorpyrifos from resistant Hyalella azteca to inland silversides (Menidia beryllina) and effects on acetylcholinesterase activity and swimming performance at varying temperatures

Chlorpyrifos, an organophosphate (OP) insecticide, is prevalent in aquatic systems globally and is often implicated in aquatic toxicity during storm events. Chlorpyrifos induces toxicity by inhibition of acetylcholinesterase (AChE) activity, which has been related to alterations to fish swimming performance. Resistance to organophosphate insecticides, including chlorpyrifos, is prevalent in populations of the epibenthic amphipod Hyalella azteca in areas with known OP exposure. Previous studies have demonstrated an elevated bioaccumulation potential of insecticide-resistant prey items, however the potential for trophic transfer of chlorpyrifos from OP-resistant prey items and associated neurotoxic effects in fish predators has not been studied. Consequently, the present study aimed to determine the potential for trophic transfer of chlorpyrifos from OP-resistant H. azteca to a known predator, the inland silverside, Menidia beryllina at two temperatures (18 and 23 °C) to simulate temperature changes associated with global climate change (GCC). Fish were fed either ¹⁴C-chlorpyrifos-dosed H. azteca or control animals for 7 d, after which total bioaccumulation, percent parent chlorpyrifos, brain AChE activity and swimming performance (ramp-U_crit) were determined. Fish fed chlorpyrifos-dosed H. azteca bioaccumulated chlorpyrifos ranging from 29.9 to 1250 ng/g lipid, demonstrating the potential for trophic transfer. Lower bioaccumulation and greater biotransformation were observed in M. beryllina at 23 °C as compared to 18 °C, though this was not statistically significant. A significant 36.5% reduction in brain AChE activity was observed in fish fed chlorpyrifos-dosed H. azteca at 23 °C only, which may be attributed to increased biotransformation of parent chlorpyrifos to more potent AChE-inhibiting metabolites. Dietary chlorpyrifos exposure had no significant effect on swimming performance in M. beryllina, though ramp-U_crit was significantly increased by 25% at 23 as compared to 18 °C. These findings confirm the potential for trophic transfer of chlorpyrifos from OP-resistant prey to fish predators and the potential for elevated temperatures to exacerbate the neurotoxic effects of chlorpyrifos.

Bioaccumulation potential of chlorpyrifos in resistant Hyalella azteca: Implications for evolutionary toxicology

Given extensive use of pesticides in agriculture, there is concern for unintended consequences to non-target species. The non-target freshwater amphipod, Hyalella azteca has been found to show resistance to the organophosphate (OP) pesticide, chlorpyrifos, resulting from an amino acid substitution in acetylcholinesterase (AChE), suggesting a selective pressure of unintended pesticide exposure. Since resistant organisms can survive in contaminated habitats, there is potential for them to accumulate higher concentrations of insecticides, increasing the risk for trophic transfer. In the present study, we estimated the uptake and elimination of chlorpyrifos in non-resistant US Lab, and resistant Ulatis Creek (ULC Resistant), H. azteca populations by conducting 24-h uptake and 48-h elimination toxicokinetic experiments with ¹⁴C-chlorpyrifos. Our results indicated that non-resistant H. azteca had a larger uptake clearance coefficient (1467 mL g^-1 h^-1) than resistant animals (557 mL g^-1 h^-1). The half-life derived from the toxicokinetic models also estimated that steady state conditions were reached at 13.5 and 32.5 h for US Lab and ULC, respectively. Bioaccumulation was compared between non-resistant and resistant H. azteca by exposing animals to six different environmentally relevant concentrations for 28 h. Detection of chlorpyrifos in animal tissues indicated that resistant animals exposed to high concentrations of chlorpyrifos were capable of accumulating the insecticide up to 10-fold higher compared to non-resistant animals. Metabolite analysis from the 28-h concentration experiments showed that between 20 and 50 % parent compound was detected in H. azteca. These results imply that bioaccumulation potential can be more significant in chlorpyrifos resistant H. azteca and may be an essential factor in assessing the full impacts of toxicants on critical food webs, especially in the face of increasing pesticide and chemical runoff.

Trends in the Application of "Omics" to Ecotoxicology and Stress Ecology

Our ability to predict and assess how environmental changes such as pollution and climate change affect components of the Earth's biome is of paramount importance. This need positioned the fields of ecotoxicology and stress ecology at the center of environmental monitoring efforts. Advances in these interdisciplinary fields depend not only on conceptual leaps but also on technological advances and data integration. High-throughput "omics" technologies enabled the measurement of molecular changes at virtually all levels of an organism's biological organization and thus continue to influence how the impacts of stressors are understood. This bibliometric review describes literature trends (2000-2020) that indicate that more different stressors than species are studied each year but that only a few stressors have been studied in more than two phyla. At the same time, the molecular responses of a diverse set of non-model species have been investigated, but cross-species comparisons are still rare. While transcriptomics studies dominated until 2016, a shift towards proteomics and multiomics studies is apparent. There is now a wealth of data at functional omics levels from many phylogenetically diverse species. This review, therefore, addresses the question of how to integrate omics information across species.

Recessivity of pyrethroid resistance and limited interspecies hybridization across Hyalella clades supports rapid and independent origins of resistance

Several populations of the amphipod, Hyalella azteca, have developed resistance to pyrethroid insecticides due to non-target exposure, but the dominance of the resistance trait is unknown. The current study investigated the dominance level of point mutations in natural populations of insecticide-resistant H. azteca and determined whether H. azteca from different clades with and without resistant alleles can hybridize and produce viable offspring. A parent generation (P₀) of non-resistant homozygous wild type H. azteca was crossbred with pyrethroid-resistant homozygous mutant animals and the tolerance of the filial 1 (F₁) generation to the pyrethroid insecticide, permethrin, was measured. Then the genotypes of the F₁ generation was examined to assure heterozygosity. The resistant parents had permethrin LC₅₀ values that ranged from 52 to 82 times higher than the non-resistant animals and both crossbreeding experiments produced heterozygous hybrid offspring that had LC₅₀ values similar to the non-resistant H. azteca parent. Dominance levels calculated for each of the crosses showed values close to 0, confirming that the L925I and L925V mutations were completely recessive. The lack of reproduction by hybrids of the C x D breeding confirmed that these clades are reproductively isolated and therefore introgression of adaptive alleles across these clades is unlikely. Potential evolutionary consequences of this selection include development of population bottlenecks, which may arise leading to fitness costs and reduced genetic diversity of H. azteca.