Genetic variation at aryl hydrocarbon receptor (AHR) loci in populations of Atlantic killifish (Fundulus heteroclitus) inhabiting polluted and reference habitats

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.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

Improving hydrocarbon toxicity tolerance in poultry: role of genes and antioxidants

Sustenance of smallholder poultry production as an alternative source of food security and income is imperative in communities exposed to hydrocarbon pollution. Exposure to hydrocarbon pollutants causes disruption of homeostasis, thereby compromising the genetic potential of the birds. Oxidative stress-mediated dysfunction of the cellular membrane is a contributing factor in the mechanism of hydrocarbon toxicity. Epidemiological studies show that tolerance to hydrocarbon exposure may be caused by the activation of genes that control disease defense pathways like aryl hydrocarbon receptor (AhR) and nuclear factor erythroid 2p45-related factor 2 (Nrf2). Disparity in the mechanism and level of tolerance to hydrocarbon fragments among species may exist and may result in variations in gene expression within individuals of the same species upon exposure. Genomic variability is critical for adaptation and serves as a survival mechanism in response to environmental pollutants. Understanding the interplay of diverse genetic mechanisms in relation to environmental influences is important for exploiting the differences in various genetic variants. Protection against pollutant-induced physiological responses using dietary antioxidants can mitigate homeostasis disruptions. Such intervention may initiate epigenetic modulation relevant to gene expression of hydrocarbon tolerance, enhancing productivity, and possibly future development of hydrocarbon-tolerant breeds.

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.

The Ah Receptor from Toxicity to Therapeutics: Report from the 5th AHR Meeting at Penn State University, USA, June 2022

The aryl hydrocarbon receptor (AHR) is a sensor of low-molecular-weight molecule signals that originate from environmental exposures, the microbiome, and host metabolism. Building upon initial studies examining anthropogenic chemical exposures, the list of AHR ligands of microbial, diet, and host metabolism origin continues to grow and has provided important clues as to the function of this enigmatic receptor. The AHR has now been shown to be directly involved in numerous biochemical pathways that influence host homeostasis, chronic disease development, and responses to toxic insults. As this field of study has continued to grow, it has become apparent that the AHR is an important novel target for cancer, metabolic diseases, skin conditions, and autoimmune disease. This meeting attempted to cover the scope of basic and applied research being performed to address possible applications of our basic knowledge of this receptor on therapeutic outcomes.

Proper modulation of AHR signaling is necessary for establishing neural connectivity and oligodendrocyte precursor cell development in the embryonic zebrafish brain

2,3,7,8-tetrachlorodibenzo-[p]-dioxin (TCDD) is a persistent global pollutant that exhibits a high affinity for the aryl hydrocarbon receptor (AHR), a ligand activated transcription factor. Epidemiological studies have associated AHR agonist exposure with multiple human neuropathologies. Consistent with the human data, research studies using laboratory models have linked pollutant-induced AHR activation to disruptions in learning and memory as well as motor impairments. Our understanding of endogenous AHR functions in brain development is limited and, correspondingly, scientists are still determining which cell types and brain regions are sensitive to AHR modulation. To identify novel phenotypes resulting from pollutant-induced AHR activation and ahr2 loss of function, we utilized the optically transparent zebrafish model. Early embryonic TCDD exposure impaired embryonic brain morphogenesis, resulted in ventriculomegaly, and disrupted neural connectivity in the optic tectum, habenula, cerebellum, and olfactory bulb. Altered neural network formation was accompanied by reduced expression of synaptic vesicle 2. Loss of ahr2 function also impaired nascent network development, but did not affect gross brain or ventricular morphology. To determine whether neural AHR activation was sufficient to disrupt connectivity, we used the Gal4/UAS system to express a constitutively active AHR specifically in differentiated neurons and observed disruptions only in the cerebellum; thus, suggesting that the phenotypes resulting from global AHR activation likely involve multiple cell types. Consistent with this hypothesis, we found that TCDD exposure reduced the number of oligodendrocyte precursor cells and their derivatives. Together, our findings indicate that proper modulation of AHR signaling is necessary for the growth and maturation of the embryonic zebrafish brain.

Intraspecific Variation in the Sublethal Effects of Physically and Chemically Dispersed Crude Oil on Early Life Stages of Atlantic Cod (Gadus morhua)

The offshore oil industry in Atlantic Canada necessitates a greater understanding of the potential impacts of oil exposure and spill response measures on cold-water marine species. We used a standardized scoring index to characterize sublethal developmental impacts of physically and chemically dispersed crude oil in early life stages of Atlantic cod (Gadus morhua) and assessed intraspecific variation in the response among cod families. Cod (origin: Scotian Shelf, Canada) were laboratory-crossed to produce embryos from five specific families, which were subsequently exposed prehatch to gradient dilutions of a water-accommodated fraction (WAF) and a chemically enhanced WAF (CEWAF; prepared with Corexit 9500A) for 24 h. Postexposure, live embryos were transferred into filtered seawater and monitored to hatch; then, all live fish had sublethal endpoints assessed using the blue-sac disease (BSD) severity index. In both WAF and CEWAF groups, increasing exposure concentrations (measured as total petroleum hydrocarbons) resulted in an increased incidence of BSD symptoms (impaired swimming ability, increased degree of spinal curvature, yolk-sac edemas) in cod across all families. This positive concentration-dependent increase in BSD was similar between physically (WAF) versus chemically (CEWAF) dispersed oil exposures, indicating that dispersant addition does not exacerbate the effect of crude oil on BSD incidence in cod. Sensitivity varied between families, with some families having less BSD than others with increasing exposure concentrations. To our knowledge, our study is the first to demonstrate the occurrence in fishes of intraspecific variation among families in sublethal responses to oil and dispersant exposure. Our results suggest that sublethal effects of crude oil exposure will not be uniformly observed across cod populations and that sensitivity depends on genetic background. Environ Toxicol Chem 2022;41:1967-1976. © 2022 SETAC.

The aryl hydrocarbon receptor: A predominant mediator for the toxicity of emerging dioxin-like compounds

The aryl hydrocarbon receptor (AHR) is a member of the basic helix-loop-helix/Per-ARNT-Sim (bHLH-PAS) family of transcription factors and has broad biological functions. Early after the identification of the AHR, most studies focused on its roles in regulating the expression of drug-metabolizing enzymes and mediating the toxicity of dioxins and dioxin-like compounds (DLCs). Currently, more diverse functions of AHR have been identified, indicating that AHR is not just a dioxin receptor. Dioxins and DLCs occur ubiquitously and have diverse health/ecological risks. Additional research is required to identify both shared and compound-specific mechanisms, especially for emerging DLCs such as polyhalogenated carbazoles (PHCZs), polychlorinated diphenyl sulfides (PCDPSs), and others, of which only a few investigations have been performed at present. Many of the toxic effects of emerging DLCs were observed to be predominantly mediated by the AHR because of their structural similarity as dioxins, and the in vitro TCDD-relative potencies of certain emerging DLC congeners are comparable to or even greater than the WHO-TEFs of OctaCDD, OctaCDF, and most coplanar PCBs. Due to the close relationship between AHR biology and environmental science, this review begins by providing novel insights into AHR signaling (canonical and non-canonical), AHR's biochemical properties (AHR structure, AHR-ligand interaction, AHR-DNA binding), and the variations during AHR transactivation. Then, AHR ligand classification and the corresponding mechanisms are discussed, especially the shared and compound-specific, AHR-mediated effects and mechanisms of emerging DLCs. Accordingly, a series of in vivo and in vitro toxicity evaluation methods based on the AHR signaling pathway are reviewed. In light of current advances, future research on traditional and emerging DLCs will enhance our understanding of their mechanisms, toxicity, potency, and ecological impacts.

Long-lasting effects of chronic exposure to chemical pollution on the hologenome of the Manila clam

Chronic exposure to pollutants affects natural populations, creating specific molecular and biochemical signatures. In the present study, we tested the hypothesis that chronic exposure to pollutants might have substantial effects on the Manila clam hologenome long after removal from contaminated sites. To reach this goal, a highly integrative approach was implemented, combining transcriptome, genetic and microbiota analyses with the evaluation of biochemical and histological profiles of the edible Manila clam Ruditapes philippinarum, as it was transplanted for 6 months from the polluted area of Porto Marghera (PM) to the clean area of Chioggia (Venice lagoon, Italy). One month post-transplantation, PM clams showed several modifications to its resident microbiota, including an overrepresentation of the opportunistic pathogen Arcobacter spp. This may be related to the upregulation of several immune genes in the PM clams, potentially representing a host response to the increased abundance of deleterious bacteria. Six months after transplantation, PM clams demonstrated a lower ability to respond to environmental/physiological stressors related to the summer season, and the hepatopancreas-associated microbiota still showed different compositions among PM and CH clams. This study confirms that different stressors have predictable effects in clams at different biological levels and demonstrates that chronic exposure to pollutants leads to long-lasting effects on the animal hologenome. In addition, no genetic differentiation between samples from the two areas was detected, confirming that PM and CH clams belong to a single population. Overall, the obtained responses were largely reversible and potentially related to phenotypic plasticity rather than genetic adaptation. The results here presented will be functional for the assessment of the environmental risk imposed by chemicals on an economically important bivalve species.

Resistance to Cyp3a induction by polychlorinated biphenyls, including non-dioxin-like PCB153, in gills of killifish (Fundulus heteroclitus) from New Bedford Harbor

Previous reports suggested that non-dioxin-like (NDL) PCB153 effects on cytochrome P450 3A (Cyp3a) expression in Atlantic killifish (Fundulus heteroclitus) gills differed between F0 generation fish from a PCB site (New Bedford Harbor; NBH) and a reference site (Scorton Creek; SC). Here, we examined effects of PCB153, dioxin-like (DL) PCB126, or a mixture of both, on Cyp3a56 mRNA in killifish generations removed from the wild, without environmental PCB exposures. PCB126 effects in liver and gills differed between populations, as expected. Gill Cyp3a56 was not affected by either congener in NBH F2 generation fish, but was induced by PCB153 in SC F1 fish, with females showing a greater response. PCB153 did not affect Cyp3a56 in liver of either population. Results suggest a heritable resistance to NDL-PCBs in killifish from NBH, in addition to that reported for DL PCBs. Induction of Cyp3a56 in gills may be a biomarker of exposure to NDL PCBs in fish populations that are not resistant to PCBs.

Dioxin exposure associated with fecundability and infertility in mothers and daughters of Seveso, Italy

STUDY QUESTION

Is there an association between 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) exposure and fecundability and infertility among Seveso women and their daughters?

SUMMARY ANSWER

TCDD exposure is associated with a decrease in fecundability and increased risk of infertility in women, as well as their daughters.

WHAT IS KNOWN ALREADY

In animal studies, maternal exposure to TCDD is associated with decreased fertility in offspring. Effects of TCDD are mediated by activation of the aryl hydrocarbon receptor (AHR) pathway.

STUDY DESIGN, SIZE, DURATION

The Seveso Women's Health Study (SWHS) has followed 981 women exposed to TCDD in a 1976 accident since 1996. In 2014, we initiated the Seveso Second Generation Study to follow-up their children.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We obtained information on pregnancy history including time of trying to conceive from SWHS women and their daughters who were 18 years or older. We considered TCDD exposure as initial 1976 serum TCDD concentration and estimated TCDD at pregnancy. We examined relationships of TCDD exposure with time to pregnancy (TTP, the monthly probability of conception within the first 12 months of trying) and infertility (≥12 months of trying to conceive). We also assessed contributions of polymorphisms in the AHR pathway via genetic risk score.

MAIN RESULTS AND THE ROLE OF CHANCE

Among SWHS women (n = 446), median TTP was 3 months and 18% reported taking ≥12 months to conceive. Initial 1976 TCDD (log10) was associated with longer TTP (adjusted fecundability odds ratio = 0.82; 95% CI 0.68-0.98) and increased risk of infertility (adjusted relative risk = 1.35; 95% CI 1.01-1.79). TCDD at pregnancy yielded similar associations. Among SWHS daughters (n = 66), median TTP was 2 months and 11% reported taking ≥12 months to conceive. Daughters showed similar, but non-significant, associations with maternal TCDD exposure.

LIMITATIONS, REASONS FOR CAUTION

A limitation of this study is time to pregnancy was reported retrospectively, although previous studies have found women are able to recall time to conception with a high degree of accuracy many years after the fact. The number of SWHS daughters who had a live birth was small and we were unable to examine fecundability of SWHS sons.

WIDER IMPLICATIONS OF THE FINDINGS

Consistent with previous findings in animal studies, our study found that TCDD exposure may be associated with decreased fertility in Seveso mothers and potentially in their daughters exposed in utero. There may be susceptible genetic subgroups. The literature has largely considered the genetics of the AHR pathway in the context of male fertility but not female fertility, despite strong biological plausibility. These findings should be replicated in larger populations and of different ancestry. Future studies in Seveso should examine the sons and the grandchildren of exposed mothers given the animal literature suggesting potential heritable epigenetic effects.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by grant numbers F06 TW02075-01 from the National Institutes of Health, R01 ES07171 and 2P30-ESO01896-17 from the National Institute of Environmental Health Sciences, R82471 from the U.S. Environmental Protection Agency and #2896 from Regione Lombardia and Fondazione Lombardia Ambiente, Milan, Italy. J.A. was supported by F31ES026488 from the National Institutes of Health. The authors declare they have no actual or potential competing financial interests.

TRIAL REGISTRATION NUMBER

N/A.

Aryl hydrocarbon receptor ligands increase ABC transporter activity and protein expression in killifish (Fundulus heteroclitus) renal proximal tubules

Many widespread and persistent organic pollutants, for example, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and some polychlorinated biphenyls, activate the aryl hydrocarbon receptor (AhR) causing it to translocate to the cell nucleus where it transactivates target genes, increasing expression of a number of xenobiotic metabolizing enzymes as well as some transporters. AhR's ability to target transporters within the kidney is essentially unexplored. We show here that exposing isolated killifish (Fundulus heteroclitus) renal proximal tubules to micromolar β-naphthoflavone (BNF) or nanomolar TCDD roughly doubled the transport activity of Multidrug resistance-associated proteins Mrp2 and Mrp4, P-glycoprotein (P-gp) and Breast cancer resistance protein (Bcrp), all ATP-driven xenobiotic efflux pumps and critical determinants of renal xenobiotic excretion. These effects were abolished by actinomycin D and cycloheximide and by the AhR antagonist, α-naphthoflavone, indicating that increased transport activity was dependent on transcription and translation as well as ligand binding to AhR. Quantitative immunostaining of renal tubules exposed to BNF and TCDD showed increased luminal membrane expression of Mrp2, Mrp4, P-gp and Bcrp. Thus, in these renal tubules, the four ABC transporters are targets of AhR action.

Phase 0 of the Xenobiotic Response: Nuclear Receptors and Other Transcription Factors as a First Step in Protection from Xenobiotics

This mini-review examines the crucial importance of transcription factors as a first line of defense in the detoxication of xenobiotics. Key transcription factors that recognize xenobiotics or xenobiotic-induced stress such as reactive oxygen species (ROS), include AhR, PXR, CAR, MTF, Nrf2, NF-κB, and AP-1. These transcription factors constitute a significant portion of the pathways induced by toxicants as they regulate phase I-III detoxication enzymes and transporters as well as other protective proteins such as heat shock proteins, chaperones, and anti-oxidants. Because they are often the first line of defense and induce phase I-III metabolism, could these transcription factors be considered the phase 0 of xenobiotic response?

Divergence, evolution and adaptation in ray-finned fish genomes

With the rapid development of next-generation sequencing technologies and bioinformatics, over 50 ray-finned fish genomes by far have been sequenced with high quality. The genomic work provides abundant genetic resources for deep understanding of divergence, evolution and adaptation in the fish genomes. They are also instructive for identification of candidate genes for functional verification, molecular breeding, and development of novel marine drugs. As an example of other omics data, the Fish-T1K project generated a big database of fish transcriptomes to integrate with these published fish genomes for potential applications. In this review, we highlight the above-mentioned recent investigations and core topics on the ray-finned fish genome research, with a main goal to obtain a deeper understanding of fish biology for theoretical and practical applications.

Characterization of the Fundulus heteroclitus embryo transcriptional response and development of a gene expression-based fingerprint of exposure for the alternative flame retardant, TBPH (bis (2-ethylhexyl)-tetrabromophthalate)

Although alternative Flame Retardant (FR) chemicals are expected to be safer than the legacy FRs they replace, their risks to human health and the environment are often poorly characterized. This study used a small volume, fish embryo system to reveal potential mechanisms of action and diagnostic exposure patterns for TBPH (bis (2-ethylhexyl)-tetrabromophthalate), a component of several widely-used commercial products. Two different concentration of TBPH were applied to sensitive early life stages of an ecologically important test species, Fundulus heteroclitus (Atlantic killifish), with a well-annotated genome. Exposed fish embryos were sampled for transcriptomics or chemical analysis of parent compound and primary metabolite or observed for development and survival through larval stage. Global transcript profiling using RNA-seq was conducted (n = 16 per treatment) to provide a non-targeted and statistically robust approach to characterize TBPH gene expression patterns. Transcriptomic analysis revealed a dose-response in the expression of genes associated with a surprisingly limited number of biological pathways, but included the aryl hydrocarbon receptor signal transduction pathway, which is known to respond to several toxicologically-important chemical classes. A transcriptional fingerprint using Random Forests was developed that was able to perfectly discriminate exposed vs. non-exposed individuals in test sets. These results suggest that TBPH has a relatively low potential for developmental toxicity (at least in fishes), despite concerns related to its structural similarities to endocrine disrupting chemicals and that the early life stage Fundulus system may provide a convenient test system for exposure characterization. More broadly, this study advances the usefulness of a biological testing and analysis system utilizing non-targeted transcriptomics profiling and early developmental endpoints that complements current screening methods to characterize chemicals of ecological and human health concern.

Evolutionary Toxicogenomics of the Striped Killifish (Fundulus majalis) in the New Bedford Harbor (Massachusetts, USA)

In this paper, we used a Genotyping-by-Sequencing (GBS) approach to find and genotype more than 4000 genome-wide SNPs (Single Nucleotide Polymorphisms) from striped killifish exposed to a variety of polychlorinated biphenyls (PCBs) and other aromatic pollutants in New Bedford Harbor (NBH, Massachusetts, USA). The aims of this study were to identify the genetic consequences of exposure to aquatic pollutants and detect genes that may be under selection. Low genetic diversity (H_E and π) was found in the site exposed to the highest pollution level, but the pattern of genetic diversity did not match the pollution levels. Extensive connectivity was detected among sampling sites, which suggests that balanced gene flow may explain the lack of genetic variation in response to pollution levels. Tests for selection identified 539 candidate outliers, but many of the candidate outliers were not shared among tests. Differences among test results likely reflect different test assumptions and the complex pollutant mixture. Potentially, selectively important loci are associated with 151 SNPs, and enrichment analysis suggests a likely involvement of these genes with pollutants that occur in NBH. This result suggests that selective processes at genes targeted by pollutants may be occurring, even at a small geographical scale, and may allow the local striped killifish to resist the high pollution levels.

Profiling microRNA expression in Atlantic killifish (Fundulus heteroclitus) gill and responses to arsenic and hyperosmotic stress

The Atlantic killifish (Fundulus heteroclitus), native to estuarine areas of the Atlantic coast of the United States, has become a valuable ecotoxicological model as a result of its ability to acclimate to rapid environmental changes and adapt to polluted habitats. MicroRNAs (miRNAs) are highly conserved small RNAs that regulate gene expression and play critical roles in stress responses in a variety of organisms. Global miRNA expression in killifish and the potential roles miRNA have in environmental acclimation have yet to be characterized. Accordingly, we profiled miRNA expression in killifish gill for the first time and identified a small group of highly expressed, well-conserved miRNAs as well as 16 novel miRNAs not yet identified in other organisms. Killifish respond to large fluctuations in salinity with rapid changes in gene expression and protein trafficking to maintain osmotic balance, followed by a secondary phase of gene and protein expression changes that enable remodeling of the gills. Arsenic, a major environmental toxicant, was previously shown to inhibit gene expression responses in killifish gill, as well the ability of killifish to acclimate to a rapid increase in salinity. Thus, we examined the individual and combined effects of salinity and arsenic on miRNA expression in killifish gill. Using small RNA sequencing, we identified 270 miRNAs expressed in killifish, and found that miR-135b was differentially expressed in response to arsenic and at 24 h following transfer to salt water. Predicted targets of miR-135b are involved in ion transport, cell motility and migration, GTPase mediated signal transduction and organelle assembly. Consistent with previous studies of these two environmental stressors, we found a significant interaction (i.e., arsenic dependent salinity effect), whereby killifish exposed to arsenic exhibited an opposite response in miR-135b expression at 24 h post hyperosmotic challenge compared to controls. By examining mRNA expression of predicted miRNA targets during salinity acclimation and arsenic exposure, we found that miR-135b targets were significantly more likely to decrease during salinity acclimation than non-targets. Our identification of a significant interaction effect of arsenic and salinity on miR-135b expression supports the hypothesis that arsenic alters upstream regulators of stress response networks, which may adversely affect the killifish response to osmotic stress. The characterization of miRNAs in this ecotoxicological model will be a valuable resource for future studies investigating the role of miRNAs in response to environmental stress.

Hepatic metabolite profiling of polychlorinated biphenyl (PCB)-resistant and sensitive populations of Atlantic killifish (Fundulus heteroclitus)

Atlantic killifish inhabiting polluted sites along the east coast of the U.S. have evolved resistance to toxic effects of contaminants. One such contaminated site is the Acushnet River estuary, near New Bedford Harbor (NBH), Massachusetts, which is characterized by very high PCB concentrations in the sediments and in the tissues of resident killifish. Though killifish at this site appear to be thriving, the metabolic costs of survival in a highly contaminated environment are not well understood. In this study we compared the hepatic metabolite profiles of resistant (NBH) and sensitive populations (Scorton Creek (SC), Sandwich, MA) using a targeted metabolomics approach in which polar metabolites were extracted from adult fish livers and quantified. Our results revealed differences in the levels of several metabolites between fish from the two sites. The majority of these metabolites are associated with one-carbon metabolism, an important pathway that supports multiple physiological processes including DNA and protein methylation, nucleic acid biosynthesis and amino acid metabolism. We measured the gene expression of DNA methylation (DNA methyltransferase 1, dnmt1) and demethylation genes (Ten-Eleven Translocation (TET) genes) in the two populations, and observed lower levels of dnmt1 and higher levels of TET gene expression in the NBH livers, suggesting possible differences in DNA methylation profiles. Consistent with this, the two populations differed significantly in the levels of 5-methylcytosine and 5-hydroxymethylcytosine nucleotides. Overall, our results suggest that the unique hepatic metabolite signatures observed in NBH and SC reflect the adaptive mechanisms for survival in their respective habitats.

Genome-wide scan reveals signatures of selection related to pollution adaptation in non-model estuarine Atlantic killifish (Fundulus heteroclitus)

In many human-altered ecosystems, organisms are increasingly faced with more diverse and complex environmental stressors and pollutant mixtures, to which the adaptations necessary to survive exposure are likely to be numerous and varied. Improving our understanding of the molecular mechanisms that underlie complex polygenic adaptations in natural settings requires significant toxicological, biochemical, physiological, and genomic data rarely available for non-model organisms. Here, we build upon two decades of study of adaptation to anthropogenic pollutants in a population of Atlantic killifish (Fundulus heteroclitus) that inhabits the creosote-contaminated Atlantic Wood Industries Superfund (AW) site on the Elizabeth River, Virginia in the United States. To better understand the genotypes that underlie previously characterized resistance to PCBs and PAHs, we performed Restriction site-Associated DNA sequencing (RADseq) on killifish from AW and two relatively clean reference sites (King's Creek-KC, and Mains Creek-MC). Across the genome, we analyzed over 83,000 loci and 12,000 single nucleotide polymorphisms (SNPs). Shared across both comparisons of killifish from polluted (AW) and relatively unpolluted (KC and MC) sites, we found eight genomic regions with smoothed FST values significantly (p < 0.001) elevated above background. Using the recently published F. heteroclitus reference genome, we identified candidate genes in these significant regions involved in the AHR pathway (e.g. AIP, ARNT1c), as well as genes relating to cardiac structure and function. These genes represent both previously characterized and potentially novel molecular adaptations involved with various aspects of resistance to these environmental toxins.

Adaptive capabilities and fitness consequences associated with pollution exposure in fish

Many fish populations are exposed to harmful levels of chemical pollution and selection pressures associated with these exposures have led to the evolution of tolerance. Our understanding of the physiological basis for these adaptations is limited, but they are likely to include processes involved with the absorption, distribution, metabolism and/or excretion of the target chemical. Other potential adaptive mechanisms include enhancements in antioxidant responses, an increased capacity for DNA and/or tissue repair and alterations to the life cycle of fish that enable earlier reproduction. Analysis of single-nucleotide polymorphism frequencies has shown that tolerance to hydrocarbon pollutants in both marine and estuarine fish species involves alteration in the expression of the xenobiotic metabolism enzyme CYP1A. In this review, we present novel data showing also that variants of the CYP1A gene have been under selection in guppies living in Trinidadian rivers heavily polluted with crude oil. Potential costs associated with these adaptations could reduce fitness in unpolluted water conditions. Integrating knowledge of local adaptation to pollution is an important future consideration in conservation practices such as for successful restocking, and improving connectivity within river systems.This article is part of the themed issue 'Human influences on evolution, and the ecological and societal consequences'.

AHR and CYP1A expression link historical contamination events to modern day developmental effects in the American alligator

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor that initiates a transcriptional pathway responsible for the expression of CYP1A subfamily members, key to the metabolism of xenobiotic compounds. Toxic planar halogenated aromatic hydrocarbons, including dioxin and PCBs, are capable of activating the AHR, and while dioxin and PCB inputs into the environment have been dramatically curbed following strict regulatory efforts in the United States, they persist in the environment and exposures remain relevant today. Little is known regarding the effects that long-term chronic exposures to dioxin or dioxin-like compounds might have on the development and subsequent health of offspring from exposed individuals, nor is much known regarding AHR expression in reptilians. Here, we characterize AHR and CYP1A gene expression in embryonic and juvenile specimen of a long-lived, apex predator, the American alligator (Alligator mississippiensis), and investigate variation in gene expression profiles in offspring collected from sites conveying differential exposures to environmental contaminants. Both age- and tissue-dependent patterning of AHR isoform expression are detected. We characterize two downstream transcriptional targets of the AHR, CYP1A1 and CYP1A2, and describe conserved elements of their genomic architecture. When comparisons across different sites are made, hepatic expression of CYP1A2, a direct target of the AHR, appears elevated in embryos from a site associated with a dioxin point source and previously characterized PCB contamination. Elevated CYP1A2 expression is not persistent, as site-specific variation was absent in juveniles originating from field-collected eggs but reared under lab conditions. Our results illustrate the patterning of AHR gene expression in a long-lived environmental model species, and indicate a potential contemporary influence of historical contamination. This research presents a novel opportunity to link contamination events to critical genetic pathways during embryonic development, and carries significant potential to inform our understanding of potential health effects in wildlife and humans.

Annotation of the Nuclear Receptors in an Estuarine Fish species, Fundulus heteroclitus

The nuclear receptors (NRs) are ligand-dependent transcription factors that respond to various internal as well as external cues such as nutrients, pheromones, and steroid hormones that play crucial roles in regulation and maintenance of homeostasis and orchestrating the physiological and stress responses of an organism. We annotated the Fundulus heteroclitus (mummichog; Atlantic killifish) nuclear receptors. Mummichog are a non-migratory, estuarine fish with a limited home range often used in environmental research as a field model for studying ecological and evolutionary responses to variable environmental conditions such as salinity, oxygen, temperature, pH, and toxic compounds because of their hardiness. F. heteroclitus have at least 74 NRs spanning all seven gene subfamilies. F. heteroclitus is unique in that no RXRα member was found within the genome. Interestingly, some of the NRs are highly conserved between species, while others show a higher degree of divergence such as PXR, SF1, and ARα. Fundulus like other fish species show expansion of the RAR (NR1B), Rev-erb (NR1D), ROR (NR1F), COUPTF (NR2F), ERR (NR3B), RXR (NR2B), and to a lesser extent the NGF (NR4A), and NR3C steroid receptors (GR/AR). Of particular interest is the co-expansion of opposing NRs, Reverb-ROR, and RAR/RXR-COUPTF.

Living in the city: urban environments shape the evolution of a native annual plant

Urban environments are warmer, have higher levels of atmospheric CO₂ and have altered patterns of disturbance and precipitation than nearby rural areas. These differences can be important for plant growth and are likely to create distinct selective environments. We planted a common garden experiment with seeds collected from natural populations of the native annual plant Lepidium virginicum, growing in five urban and nearby rural areas in the northern United States to determine whether and how urban populations differ from those from surrounding rural areas. When grown in a common environment, plants grown from seeds collected from urban areas bolted sooner, grew larger, had fewer leaves, had an extended time between bolting and flowering, and produced more seeds than plants grown from seeds collected from rural areas. Interestingly, the rural populations exhibited larger phenotypic differences from one another than urban populations. Surprisingly, genomic data revealed that the majority of individuals in each of the urban populations were more closely related to individuals from other urban populations than they were to geographically proximate rural areas - the one exception being urban and rural populations from New York which were nearly identical. Taken together, our results suggest that selection in urban environments favors different traits than selection in rural environments and that these differences can drive adaptation and shape population structure.

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.

Applying evolutionary genetics to developmental toxicology and risk assessment

Evolutionary thinking continues to challenge our views on health and disease. Yet, there is a communication gap between evolutionary biologists and toxicologists in recognizing the connections among developmental pathways, high-throughput screening, and birth defects in humans. To increase our capability in identifying potential developmental toxicants in humans, we propose to apply evolutionary genetics to improve the experimental design and data interpretation with various in vitro and whole-organism models. We review five molecular systems of stress response and update 18 consensual cell-cell signaling pathways that are the hallmark for early development, organogenesis, and differentiation; and revisit the principles of teratology in light of recent advances in high-throughput screening, big data techniques, and systems toxicology. Multiscale systems modeling plays an integral role in the evolutionary approach to cross-species extrapolation. Phylogenetic analysis and comparative bioinformatics are both valuable tools in identifying and validating the molecular initiating events that account for adverse developmental outcomes in humans. The discordance of susceptibility between test species and humans (ontogeny) reflects their differences in evolutionary history (phylogeny). This synthesis not only can lead to novel applications in developmental toxicity and risk assessment, but also can pave the way for applying an evo-devo perspective to the study of developmental origins of health and disease.

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.

Evolutionary toxicology in an omics world

Evolutionary toxicology is a young field that has grown rapidly in the past two decades. The potential of this field comes from the ability to link chemical contamination to multigenerational and population-wide effects in various species. The advancements and rapidly decreasing costs of -omic tools are improving the power and resolution of evolutionary toxicology studies. In this manuscript, we aim to address the trajectories and perspectives for conducting evolutionary toxicology studies with -omic approaches. We discuss the complementarity of using multiple -omic tools (genomics, eDNA, transcriptomics, proteomics, and metabolomics) for utility in understanding the toxicological relevance of adaptive responses in populations. In addition, we discuss phenotypic plasticity and its relevance to transcriptomic studies in toxicology. As evolutionary toxicology grows and expands its capacity to link toxicology with population-wide end points, we emphasize the applications of such studies in answering questions about ecological and population health, as well as future applicability to regulation. Thus, we aim to emphasize the enormous potential for evolutionary toxicology in an -omics world and give perspectives on the directions of future investigations.

Diversity as Opportunity: Insights from 600 Million Years of AHR Evolution

The aryl hydrocarbon receptor (AHR) was for many years of interest only to pharmacologists and toxicologists. However, this protein has fundamental roles in biology that are being revealed through studies in diverse animal species. The AHR is an ancient protein. AHR homologs exist in most major groups of modern bilaterian animals, including deuterostomes (chordates, hemichordates, echinoderms) and the two major clades of protostome invertebrates [ecdysozoans (e.g. arthropods and nematodes) and lophotrochozoans (e.g. molluscs and annelids)]. AHR homologs also have been identified in cnidarians such as the sea anemone Nematostella and in the genome of Trichoplax, a placozoan. Bilaterians, cnidarians, and placozoans form the clade Eumetazoa, whose last common ancestor lived approximately 600 million years ago (MYA). The presence of AHR homologs in modern representatives of all these groups indicates that the original eumetazoan animal possessed an AHR homolog. Studies in invertebrates and vertebrates reveal parallel functions of AHR in the development and function of sensory neural systems, suggesting that these may be ancestral roles. Vertebrate animals are characterized by the expansion and diversification of AHRs, via gene and genome duplications, from the ancestral protoAHR into at least five classes of AHR-like proteins: AHR, AHR1, AHR2, AHR3, and AHRR. The evolution of multiple AHRs in vertebrates coincided with the acquisition of high-affinity binding of halogenated and polynuclear aromatic hydrocarbons and the emergence of adaptive functions involving regulation of xenobiotic-metabolizing enzymes and roles in adaptive immunity. The existence of multiple AHRs may have facilitated subfunction partitioning and specialization of specific AHR types in some taxa. Additional research in diverse model and non-model species will continue to enrich our understanding of AHR and its pleiotropic roles in biology and toxicology.

The genomic landscape of rapid repeated evolutionary adaptation to toxic pollution in wild fish

Atlantic killifish populations have rapidly adapted to normally lethal levels of pollution in four urban estuaries. Through analysis of 384 whole killifish genome sequences and comparative transcriptomics in four pairs of sensitive and tolerant populations, we identify the aryl hydrocarbon receptor-based signaling pathway as a shared target of selection. This suggests evolutionary constraint on adaptive solutions to complex toxicant mixtures at each site. However, distinct molecular variants apparently contribute to adaptive pathway modification among tolerant populations. Selection also targets other toxicity-mediating genes and genes of connected signaling pathways; this indicates complex tolerance phenotypes and potentially compensatory adaptations. Molecular changes are consistent with selection on standing genetic variation. In killifish, high nucleotide diversity has likely been a crucial substrate for selective sweeps to propel rapid adaptation.

AHR-related activities in a creosote-adapted population of adult atlantic killifish, Fundulus heteroclitus, two decades post-EPA superfund status at the Atlantic Wood Site, Portsmouth, VA USA

Atlantic killifish, Fundulus heteroclitus, are adapted to creosote-based PAHs at the US EPA Superfund site known as Atlantic Wood (AW) on the southern branch of the Elizabeth River, VA USA. Subsequent to the discovery of the AW population in the early 1990s, these fish were shown to be recalcitrant to CYP1A induction by PAHs under experimental conditions, and even to the time of this study, killifish embryos collected from the AW site are resistant to developmental deformities typically associated with exposure to PAHs in reference fish. Historically, however, 90 +% of the adult killifish at this site have proliferative hepatic lesions including cancer of varying severity. Several PAHs at this site are known to be ligands for the aryl hydrocarbon receptor (AHR). In this study, AHR-related activities in AW fish collected between 2011 and 2013 were re-examined nearly 2 decades after first discovery. This study shows that CYP1A mRNA expression is three-fold higher in intestines of AW killifish compared to a reference population. Using immunohistochemistry, CYP1A staining in intestines was uniformly positive compared to negative staining in reference fish. Livers of AW killifish were examined by IHC to show that CYP1A and AHR2 protein expression reflect lesions-specific patterns, probably representing differences in intrinsic cellular physiology of the spectrum of proliferative lesions comprising the hepatocarcinogenic process. We also found that COX2 mRNA expression levels were higher in AW fish livers compared to those in the reference population, suggesting a state of chronic inflammation. Overall, these findings suggest that adult AW fish are responsive to AHR signaling, and do express CYP1A and AHR2 proteins in intestines at a level above what was observed in the reference population.

Advancing toxicology research using in vivo high throughput toxicology with small fish models

Small freshwater fish models, especially zebrafish, offer advantages over traditional rodent models, including low maintenance and husbandry costs, high fecundity, genetic diversity, physiology similar to that of traditional biomedical models, and reduced animal welfare concerns. The Collaborative Workshop on Aquatic Models and 21st Century Toxicology was held at North Carolina State University on May 5-6, 2014, in Raleigh, North Carolina, USA. Participants discussed the ways in which small fish are being used as models to screen toxicants and understand mechanisms of toxicity. Workshop participants agreed that the lack of standardized protocols is an impediment to broader acceptance of these models, whereas development of standardized protocols, validation, and subsequent regulatory acceptance would facilitate greater usage. Given the advantages and increasing application of small fish models, there was widespread interest in follow-up workshops to review and discuss developments in their use. In this article, we summarize the recommendations formulated by workshop participants to enhance the utility of small fish species in toxicology studies, as well as many of the advances in the field of toxicology that resulted from using small fish species, including advances in developmental toxicology, cardiovascular toxicology, neurotoxicology, and immunotoxicology. We alsoreview many emerging issues that will benefit from using small fish species, especially zebrafish, and new technologies that will enable using these organisms to yield results unprecedented in their information content to better understand how toxicants affect development and health.

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.

Adaptive Significance of ERα Splice Variants in Killifish (Fundulus heteroclitus) Resident in an Estrogenic Environment

The possibility that chronic, multigenerational exposure to environmental estrogens selects for adaptive hormone-response phenotypes is a critical unanswered question. Embryos/larvae of killifish from an estrogenic-polluted environment (New Bedford Harbor, MA [NBH]) compared with those from a reference site overexpress estrogen receptor alpha (ERα) mRNA but are hyporesponsive to estradiol. Analysis of ERα mRNAs in the two populations revealed differences in splicing of the gene encoding ERα (esr1). Here we tested the transactivation functions of four differentially expressed ERα mRNAs and tracked their association with the hyporesponsive phenotype for three generations after transfer of NBH parents to a clean environment. Deletion variants ERαΔ6 and ERαΔ6-8 were specific to NBH killifish, had dominant negative functions in an in vitro reporter assay, and were heritable. Morpholino-mediated induction of ERαΔ6 mRNA in zebrafish embryos verified its role as a dominant negative ER on natural estrogen-responsive promoters. Alternate long (ERαL) and short (ERαS) 5'-variants were similar transcriptionally but differed in estrogen responsiveness (ERαS ≫ ERαL). ERαS accounted for high total ERα expression in first generation (F1) NBH embryos/larvae but this trait was abolished by transfer to clean water. By contrast, the hyporesponsive phenotype of F1 NBH embryos/larvae persisted after long-term laboratory holding but reverted to a normal or hyper-responsive phenotype after two or three generations, suggesting the acquisition of physiological or biochemical traits that compensate for ongoing expression of negative-acting ERαΔ6 and ERαΔ6-8 isoforms. We conclude that a heritable change in the pattern of alternative splicing of ERα pre-mRNA is part of a genetic adaptive response to estrogens in a polluted environment.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

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.

Regulation of pregnane-X-receptor, CYP3A and P-glycoprotein genes in the PCB-resistant killifish (Fundulus heteroclitus) population from New Bedford Harbor

Killifish survive and reproduce in the New Bedford Harbor (NBH) in Massachusetts (MA), USA, a site severely contaminated with polychlorinated biphenyls (PCBs) for decades. Levels of 22 different PCB congeners were analyzed in liver from killifish collected in 2008. Concentrations of dioxin-like PCBs in liver of NBH killifish were ∼400 times higher, and the levels of non-dioxin-like PCBs ∼3000 times higher than in killifish from a reference site, Scorton Creek (SC), MA. The NBH killifish are known to be resistant to the toxicity of dioxin-like compounds and to have a reduced aryl hydrocarbon receptor (AhR) signaling response. Little is known about the responses of these fish to non-dioxin-like PCBs, which are at extraordinarily high levels in NBH fish. In mammals, some non-dioxin-like PCB congeners act through nuclear receptor 1I2, the pregnane-X-receptor (PXR). To explore this pathway in killifish, a PXR cDNA was sequenced and its molecular phylogenetic relationship to other vertebrate PXRs was determined. Killifish were also collected in 2009 from NBH and SC, and after four months in the laboratory they were injected with a single dose of either the dioxin-like PCB 126 (an AhR agonist) or the non-dioxin-like PCB 153 (a mammalian PXR agonist). Gills and liver were sampled three days after injection and transcript levels of genes encoding PXR, cytochrome P450 3A (CYP3A), P-glycoprotein (Pgp), AhR2 and cytochrome P450 1A (CYP1A) were measured by quantitative PCR. As expected, there was little effect of PCB exposure on mRNA expression of AhR2 or CYP1A in liver and gills of NBH fish. In NBH fish, but not in SC fish, there was increased mRNA expression of hepatic PXR, CYP3A and Pgp upon exposure to either of the two PCB congeners. However, basal PXR and Pgp mRNA levels in liver of NBH fish were significantly lower than in SC fish. A different pattern was seen in gills, where there were no differences in basal mRNA expression of these genes between the two populations. In SC fish, but not in NBH fish, there was increased mRNA expression of branchial PXR and CYP3A upon exposure to PCB126 and of CYP3A upon exposure to PCB153. The results suggest a difference between the two populations in non-AhR transcription factor signaling in liver and gills, and that this could involve killifish PXR. It also implies possible cross-regulatory interactions between that factor (presumably PXR) and AhR2 in liver of these fish.

Targeted mutagenesis of aryl hydrocarbon receptor 2a and 2b genes in Atlantic killifish (Fundulus heteroclitus)

Understanding molecular mechanisms of toxicity is facilitated by experimental manipulations, such as disruption of function by gene targeting, that are especially challenging in non-standard model species with limited genomic resources. While loss-of-function approaches have included gene knock-down using morpholino-modified oligonucleotides and random mutagenesis using mutagens or retroviruses, more recent approaches include targeted mutagenesis using zinc finger nuclease (ZFN), transcription activator-like effector nuclease (TALENs) and clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology. These latter methods provide more accessible opportunities to explore gene function in non-traditional model species. To facilitate evaluation of toxic mechanisms for important categories of aryl hydrocarbon pollutants, whose actions are known to be receptor mediated, we used ZFN and CRISPR-Cas9 approaches to generate aryl hydrocarbon receptor 2a (AHR2a) and AHR2b gene mutations in Atlantic killifish (Fundulus heteroclitus) embryos. This killifish is a particularly valuable non-traditional model, with multiple paralogs of AHR whose functions are not well characterized. In addition, some populations of this species have evolved resistance to toxicants such as halogenated aromatic hydrocarbons. AHR-null killifish will be valuable for characterizing the role of the individual AHR paralogs in evolved resistance, as well as in normal development. We first used five-finger ZFNs targeting exons 1 and 3 of AHR2a. Subsequently, CRISPR-Cas9 guide RNAs were designed to target regions in exon 2 and 3 of AHR2a and AHR2b. We successfully induced frameshift mutations in AHR2a exon 3 with ZFN and CRISPR-Cas9 guide RNAs, with mutation frequencies of 10% and 16%, respectively. In AHR2b, mutations were induced using CRISPR-Cas9 guide RNAs targeting sites in both exon 2 (17%) and exon 3 (63%). We screened AHR2b exon 2 CRISPR-Cas9-injected embryos for off-target effects in AHR paralogs. No mutations were observed in closely related AHR genes (AHR1a, AHR1b, AHR2a, AHRR) in the CRISPR-Cas9-injected embryos. Overall, our results demonstrate that targeted genome-editing methods are efficient in inducing mutations at specific loci in embryos of a non-traditional model species, without detectable off-target effects in paralogous genes.

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

Reproductive and endocrine disruption is commonly reported in aquatic species exposed to complex contaminant mixtures. We previously reported that Atlantic killifish (Fundulus heteroclitus) from the chronically contaminated Newark Bay, NJ, exhibit multiple endocrine disrupting effects, including inhibition of vitellogenesis (yolk protein synthesis) in females and false negative vitellogenin biomarker responses in males. Here, we characterized the effects on estrogen signaling and the transcriptional regulation of estrogen-responsive genes in this model population. First, a dose-response study tested the hypothesis that reproductive biomarkers (vtg1, vtg2, chg H, chg Hm, chg L) in Newark Bay killifish are relatively less sensitive to 17β-estradiol at the transcriptional level, relative to a reference (Tuckerton, NJ) population. The second study assessed expression for various metabolism (cyp1a, cyp3a30, mdr) and estrogen receptor (ER α, ER βa, ER βb) genes under basal and estrogen treatment conditions in both populations. Hepatic metabolism of 17β-estradiol was also evaluated in vitro as an integrated endpoint for adverse effects on metabolism. In the third study, gene methylation was evaluated for promoters of vtg1 (8 CpGs) and vtg2 (10 CpGs) in both populations, and vtg1 promoter sequences were examined for single nucleotide polymorphism (SNPs). Overall, these studies show that multi-chemical exposures at Newark Bay have desensitized all reproductive biomarkers tested to estrogen. For example, at 10ng/g 17β-estradiol, inhibition of gene induction ranged from 62% to 97% for all genes tested in the Newark Bay population, relative to induction levels in the reference population. The basis for this recalcitrant phenotype could not be explained by a change in 17β-estradiol metabolism, nuclear estrogen receptor expression, promoter methylation (gene silencing) or SNPs, all of which were unaltered and normal in the Newark Bay population. The decreased transcriptional sensitivity of estrogen-responsive genes is suggestive of a broad effect on estrogen receptor pathway signaling, and provides insight into the mechanisms of the endocrine disrupting effects in the Newark Bay population.

Mercury reduces avian reproductive success and imposes selection: an experimental study with adult- or lifetime-exposure in zebra finch

Mercury is a global pollutant that biomagnifies in food webs, placing wildlife at risk of reduced reproductive fitness and survival. Songbirds are the most diverse branch of the avian evolutionary tree; many are suffering persistent and serious population declines and we know that songbirds are frequently exposed to mercury pollution. Our objective was to determine the effects of environmentally relevant doses of mercury on reproductive success of songbirds exposed throughout their lives or only as adults. The two modes of exposure simulated philopatric species versus dispersive species, and are particularly relevant because of the heightened mercury-sensitivity of developing nervous systems. We performed a dosing study with dietary methylmercury in a model songbird species, the zebra finch (Taeniopygia guttata), at doses from 0.3 – 2.4 parts per million. Birds were exposed to mercury either as adults only or throughout their lives. All doses of mercury reduced reproductive success, with the lowest dose reducing the number of independent offspring produced in one year by 16% and the highest dose, representing approximately half the lethal dose for this species, causing a 50% reduction. While mercury did not affect clutch size or survivorship, it had the most consistent effect on the proportion of chicks that fledged from the nest, regardless of mode of exposure. Among birds exposed as adults, mercury caused a steep increase in the latency to re-nest after loss of a clutch. Birds exposed for their entire lifetimes, which were necessarily the offspring of dosed parents, had up to 50% lower reproductive success than adult-exposed birds at low doses of methylmercury, but increased reproductive success at high doses, suggesting selection for mercury tolerance at the highest level of exposure. Our results indicate that mercury levels in prey items at contaminated sites pose a significant threat to populations of songbirds through reduced reproductive success.

Targeted approach to identify genetic loci associated with evolved dioxin tolerance in Atlantic killifish (Fundulus heteroclitus)

Background

The most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed.

Results

Since the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs.

Conclusions

The observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.