Chronic exposure of killifish to a highly polluted environment desensitizes estrogen-responsive reproductive and biomarker genes

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts that adaptation is due to a few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and using RAD-seq genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that one to two large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling. One QTL locus was shared across all populations and another was shared across three populations. One QTL locus showed strong signatures of recent natural selection in the corresponding wild population but another QTL locus did not. Some candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts the influence of few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that a few large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling, where some (but not all) of these QTL loci were shared across all populations, and some (but not all) of these loci showed signatures of recent natural selection in the corresponding wild population. Some strong candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

Feminizing effects of ethinylestradiol in roach (Rutilus rutilus) populations with different estrogenic pollution exposure histories

Experimental exposures aimed at assessing the risks posed by estrogens in waste-water treatment work (WwTW) effluents to fish populations have rarely considered whether populations differ in their sensitivity to estrogenic compounds. This is despite evidence that selection at genes involved in the estrogen response has occurred in wild populations, and evidence that genotype can influence estrogen-response. In this study we compare the effects of a two-year exposure to a low measured concentration (1.3 ng/L) of ethinylestradiol (EE2) on the sexual development of roach (Rutilus rutilus) whose parental generation was sampled from two river stretches heavily contaminated with WwTW effluent and from two without any known WwTW effluent contamination. Exposure to EE2 significantly reduced the proportion of genetic males and induced a range of feminized phenotypes in males. Significantly, exposure also increased the proportion of genetic females with vitellogenic oocytes from 51 to 96%, raising the possibility that estrogen pollution could impact populations of annually spawning fish species through advancing female reproduction by at least a year. However, there was no evidence that river origin affected sensitivity to estrogens in either sex. Thus, we conclude that chronic exposure to low level EE2 has reproductive health outcomes for both male and female roach, but we find no evidence that the nature or magnitude of the response is affected by the population origin.

When evolution is the solution to pollution: Key principles, and lessons from rapid repeated adaptation of killifish (Fundulus heteroclitus) populations

For most species, evolutionary adaptation is not expected to be sufficiently rapid to buffer the effects of human‐mediated environmental changes, including environmental pollution. Here we review how key features of populations, the characteristics of environmental pollution, and the genetic architecture underlying adaptive traits, may interact to shape the likelihood of evolutionary rescue from pollution. Large populations of Atlantic killifish (Fundulus heteroclitus) persist in some of the most contaminated estuaries of the United States, and killifish studies have provided some of the first insights into the types of genomic changes that enable rapid evolutionary rescue from complexly degraded environments. We describe how selection by industrial pollutants and other stressors has acted on multiple populations of killifish and posit that extreme nucleotide diversity uniquely positions this species for successful evolutionary adaptation. Mechanistic studies have identified some of the genetic underpinnings of adaptation to a well‐studied class of toxic pollutants; however, multiple genetic regions under selection in wild populations seem to reflect more complex responses to diverse native stressors and/or compensatory responses to primary adaptation. The discovery of these pollution‐adapted killifish populations suggests that the evolutionary influence of anthropogenic stressors as selective agents occurs widely. Yet adaptation to chemical pollution in terrestrial and aquatic vertebrate wildlife may rarely be a successful “solution to pollution” because potentially adaptive phenotypes may be complex and incur fitness costs, and therefore be unlikely to evolve quickly enough, especially in species with small population sizes.

Hepatic transcriptional responses to copper in the three-spined stickleback are affected by their pollution exposure history

Some fish populations inhabiting contaminated environments show evidence of increased chemical tolerance, however the mechanisms contributing to this tolerance, and whether this is heritable, are poorly understood. We investigated the responses of two populations of wild three-spined stickleback (Gasterosteus aculeatus) with different histories of contaminant exposure to an oestrogen and copper, two widespread aquatic pollutants. Male stickleback originating from two sites, the River Aire, with a history of complex pollution discharges, and Siblyback Lake, with a history of metal contamination, were depurated and then exposed to copper (46μg/L) and the synthetic oestrogen ethinyloestradiol (22ng/L). The hepatic transcriptomic response was compared between the two populations and to a reference population with no known history of exposure (Houghton Springs, Dorset). Gene responses included those typical for both copper and oestrogen, with no discernable difference in response to oestrogen between populations. There was, however, some difference in the magnitude of response to copper between populations. Siblyback fish showed an elevated baseline transcription of genes encoding metallothioneins and a lower level of metallothionein induction following copper exposure, compared to those from the River Aire. Similarly, a further experiment with an F1 generation of Siblyback fish bred in the laboratory found evidence for elevated transcription of genes encoding metallothioneins in unexposed fish, together with an altered transcriptional response to 125μg/L copper, compared with F1 fish originating from the clean reference population exposed to the same copper concentration. These data suggest that the stickleback from Siblyback Lake have a differential response to copper, which is inherited by the F1 generation in laboratory conditions, and for which the underlying mechanism may include an elevation of baseline transcription of genes encoding metallothioneins. The genetic and/or epigenetic mechanisms contributing to this inherited alteration of metallothionein transcription have yet to be established.

Cross-resistance in Gulf killifish (Fundulus grandis) populations resistant to dioxin-like compounds

The Houston Ship Channel (HSC) in Houston, Texas is an aquatic environment with a long history of contamination, including polychlorinated dibenzodioxins (PCDD), polychlorinated dibenzofurans (PCDF), polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), and heavy metals. Populations of Gulf killifish (Fundulus grandis) from the HSC have adapted to resist developmental cardiac deformities caused by dioxin-like compounds (DLCs). Contaminants in the HSC have acted as a strong selective pressure on resident Gulf killifish populations. Rapid adaptation can lead to fitness costs, some as a direct result of the mechanisms involved in the adaptive process, whereas other adaptations may be more general. To explore potential fitness costs, we evaluated two Gulf killifish populations with documented resistance to DLC-induced cardiac teratogenesis (Patrick Bayou and Vince Bayou), and one previously characterized reference population (Gangs Bayou). We also characterized a previously unstudied population from Galveston Bay as an additional reference population (Smith Point). We tested the sensitivity of F1 larvae from these four populations to two classes of pesticides (pyrethroid (permethrin) and carbamate (carbaryl)) and two model pro-oxidants (tert-butyl hydroquinone (tBHQ) and tert-butyl hydroperoxide (tBOOH)). In addition, we explored their responses to hypoxia and measured resting metabolic rates (M.O2). Both adapted populations were cross-resistant to the toxicity of carbaryl and both pro-oxidants tested. There were no population differences in sensitivity to permethrin. On the other hand, one reference population (Gangs Bayou) was less sensitive to hypoxia, and maintained a lower M.O2 . However, there were no differences in hypoxia tolerance or resting metabolic rate between the second reference and the two adapted populations. This investigation emphasizes the importance of including multiple reference populations to clearly link fitness costs or cross-resistance to pollution adaptation, rather than to unrelated environmental or ecological differences. When compared to previous literature on adapted populations of Fundulus heteroclitus, we see a mixture of similarities and differences, suggesting that F. grandis adapted phenotypes likely involve multiple mechanisms, which may not be completely consistent among adapted populations.

Cloning of multiple ERα mRNA variants in killifish (Fundulus heteroclitus), and differential expression by tissue type, stage of reproduction, and estrogen exposure in fish from polluted and unpolluted environments

To test the hypothesis that alternative splicing could be an adaptive mechanism for populations subject to multi-generational estrogenic exposures, we compared estrogen receptor alpha (ERα) splicing variants in two populations of killifish (Fundulus heteroclitus): one resident in an estrogenic polluted environment (New Bedford Harbor, NBH, MA, USA) and one from a relatively uncontaminated reference site (Scorton Creek, SC, MA, USA). In total we identified 19 ERα variants, each with deletions of one or more coding exons. Four of the variants with potential functional relevance were analyzed by qPCR to test for population differences in expression by tissue type, site, sex, seasonal reproductive status and estrogen treatment. Significantly, a 5'-truncated short form variant (ERαS) was highly expressed in liver and ovary, and was associated with seasonal reproductive activity in SC but not NBH fish. Both ERαS and the full-length long variant (ERαL) were estrogen-inducible (ERαS>ERαL) but the induction response was lower in NBH than in SC fish. In contrast, NBH killifish were hyper-responsive to estrogen as measured by expression of two other estrogen responsive genes: vitellogenin (Vtg) and aromatase B (AroB). Most strikingly, two ERα deletion variants (Δ6 and Δ6-8), lacking ligand binding and activation function domains, were identified in a subset of NBH fish, where they were associated with reduced responsiveness to estrogen treatment. Together, these results support the hypothesis that alternative splicing of the esr1 gene of killifish could be an autoregulatory mechanism by which estrogen modulates the differential expression of ERα, and suggests a novel and adaptive mechanistic response to xenoestrogenic exposure.