Development of the morpholino gene knockdown technique in Fundulus heteroclitus: a tool for studying molecular mechanisms in an established environmental model

Bioenergetic Effects of Polycyclic Aromatic Hydrocarbon Resistance Manifest Later in Life in Offspring of Fundulus heteroclitus from the Elizabeth River

Shifts in key physiological processes can confer resistance to chemical pollutants. However, these adaptations may come with certain trade-offs, such as altered energy metabolic processes, as evident in Atlantic killifish (Fundulus heteroclitus) in Virginia's Elizabeth River (ER) that have evolved resistance to polycyclic aromatic hydrocarbons (PAHs). We seek to understand the bioenergetic costs of PAH resistance among subpopulations of Atlantic killifish with differing contamination levels in order to examine how these changes manifest across multiple life stages and how these costs might be exacerbated by additional stressors. Bioenergetics data revealed differences in metabolic rates between offspring of PAH-resistant fish and reference fish were absent or minimal in both the embryo and larval stages but pronounced at the juvenile life stage, suggesting that bioenergetic changes in pollution-adapted killifish manifest later in life. We also provide evidence that killifish from remediated sites are more sensitive to PAH exposure than killifish from nonremediated sites, suggesting loss of PAH tolerance following relaxed selection. Collectively, our data suggest that the fitness consequences associated with evolved resistance to anthropogenic stressors may manifest differently over time and depend on the magnitude of the selection pressure. This information can be valuable in effective risk and remediation assessments as well as in broadening our understanding of species responses to environmental change.

Heart development in two populations of Atlantic killifish (Fundulus heteroclitus) following exposure to a polycyclic aromatic hydrocarbon mixture

Historic industrial pollution of the Elizabeth River, Virginia resulted in polycyclic aromatic hydrocarbon (PAH) contamination in sediments. Atlantic killifish (Fundulus heteroclitus) inhabiting the Atlantic Wood (AW) industrial site adapted to complex PAH mixture at this Superfund site. Their embryos have proved highly resistant to cardiac abnormalities indicative of PAH toxicity. In this study, embryos spawned from adults collected at AW and King's Creek (KC), a reference site, were exposed at 24 h post fertilization (hpf) to Elizabeth River Sediment Extract (ERSE), a complex PAH mixture, in a range of concentrations (0, 5.04, 50.45, 100.90, 151.35, or 252.25 µg/L total PAHs). Embryos were processed for histology at 144 hpf to enable evaluations of hearts at tissue and cellular levels. Morphometry and severity scoring were used to evaluate the extent of alterations. Unexposed embryos were similar in both populations. ERSE exposure resulted in multiple changes to hearts of KC embryos but not AW. Alterations were particularly evident in KC embryos exposed to concentrations above 1% ERSE (50.45 µg/L), which had thinner ventricular walls and larger pericardial edema. Individuals with moderate pericardial edema maintained arrangement and proximity of heart chambers, but changes were seen in ventricular myocytes. Severe pericardial edema was prevalent in exposed KC embryos and typically resulted in tube heart formation. Ventricles of tube hearts had very thin walls composed of small, basophilic cells and lacked trabeculae. Edematous pericardial fluid contained small amounts of proteinaceous material, as did controls, and was free of cells. This fluid was primarily unstained, suggesting water influx due to increased permeability. The use of histological approaches provided more specific detail for tissue and cellular effects in hearts of embryos exposed to PAHs and enabled understanding of potential links to later life effects of early life exposure.

Embryonic Fundulus heteroclitus responses to sediment extracts from differentially contaminated sites in the Elizabeth River, VA

Sites along the Elizabeth River are contaminated with polycyclic aromatic hydrocarbons (PAHs) from historical creosote production and other industrial processes. Previous studies have demonstrated that Atlantic killifish collected from sites throughout the Elizabeth River display resistance to the teratogenic effects of PAH-exposure in a manner commensurate with sediment PAH concentrations. The current study characterized various chemical pollutants in sediment and investigated the effects of aqueous sediment extracts from sites along the Elizabeth River to the cardiac development of Atlantic killifish embryos from fish collected from an uncontaminated reference site. Embryonic cardiac deformities were more prevalent after exposure to extracts from sites with high PAH loads. However, activation of cytochrome P4501A, a gene up-regulated by PAH-induction of the aryl hydrocarbon receptor and measured using an in ovo EROD assay, did not consistently increase with PAH concentrations. This work further characterizes sediments in the Elizabeth River, as well as provides insight into the evolutionary pressures at each ER site.

Loss-of-function approaches in comparative physiology: is there a future for knockdown experiments in the era of genome editing?

Loss-of-function technologies, such as morpholino- and RNAi-mediated gene knockdown, and TALEN- and CRISPR/Cas9-mediated gene knockout, are widely used to investigate gene function and its physiological significance. Here, we provide a general overview of the various knockdown and knockout technologies commonly used in comparative physiology and discuss the merits and drawbacks of these technologies with a particular focus on research conducted in zebrafish. Despite their widespread use, there is an ongoing debate surrounding the use of knockdown versus knockout approaches and their potential off-target effects. This debate is primarily fueled by the observations that, in some studies, knockout mutants exhibit phenotypes different from those observed in response to knockdown using morpholinos or RNAi. We discuss the current debate and focus on the discrepancies between knockdown and knockout phenotypes, providing literature and primary data to show that the different phenotypes are not necessarily a direct result of the off-target effects of the knockdown agents used. Nevertheless, given the recent evidence of some knockdown phenotypes being recapitulated in knockout mutants lacking the morpholino or RNAi target, we stress that results of knockdown experiments need to be interpreted with caution. We ultimately argue that knockdown experiments should not be discontinued if proper control experiments are performed, and that with careful interpretation, knockdown approaches remain useful to complement the limitations of knockout studies (e.g. lethality of knockout and compensatory responses).

Small noncoding RNA profiles along alternative developmental trajectories in an annual killifish

Embryonic development of Austrofundulus limnaeus can occur along two phenotypic trajectories that are physiologically and biochemically distinct. Phenotype appears to be influenced by maternal provisioning based on the observation that young females produce predominately non-diapausing embryos and older females produce mostly diapausing embryos. Embryonic incubation temperature can override this pattern and alter trajectory. We hypothesized that temperature-induced phenotypic plasticity may be regulated by post-transcriptional modification via noncoding RNAs. As a first step to exploring this possibility, RNA-seq was used to generate transcriptomic profiles of small noncoding RNAs in embryos developing along the two alternative trajectories. We find distinct profiles of mature sequences belonging to the miR-10 family expressed in increasing abundance during development and mature sequences of miR-430 that follow the opposite pattern. Furthermore, miR-430 sequences are enriched in escape trajectory embryos. MiR-430 family members are known to target maternally provisioned mRNAs in zebrafish and may operate similarly in A. limnaeus in the context of normal development, and also by targeting trajectory-specific mRNAs. This expression pattern and function for miR-430 presents a potentially novel model for maternal-embryonic conflict in gene regulation that provides the embryo the ability to override maternal programming in the face of altered environmental conditions.

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.

Later life swimming performance and persistent heart damage following subteratogenic PAH mixture exposure in the Atlantic killifish (Fundulus heteroclitus)

High-level, acute exposures to individual polycyclic aromatic hydrocarbons (PAHs) and complex PAH mixtures result in cardiac abnormalities in developing fish embryos. Whereas acute PAH exposures can be developmentally lethal, little is known about the later life consequences of early life, lower level PAH exposures in survivors. A population of PAH-adapted Fundulus heteroclitus from the PAH-contaminated Superfund site, Atlantic Wood Industries, Elizabeth River, Portsmouth, Virginia, United States, is highly resistant to acute PAH cardiac teratogenicity. We sought to determine and characterize long-term swimming performance and cardiac histological alterations of a subteratogenic PAH mixture exposure in both reference killifish and PAH-adapted Atlantic Wood killifish embryos. Killifish from a relatively uncontaminated reference site, King's Creek, Virginia, United States, and Atlantic Wood killifish were treated with dilutions of Elizabeth River sediment extract at 24 h post fertilization (hpf). Two proven subteratogenic dilutions, 0.1 and 1.0% Elizabeth River sediment extract (total PAH 5.04 and 50.4 µg/L, respectively), were used for embryo exposures. Then, at 5-mo post hatching, killifish were subjected to a swim performance test. A separate subset of these individuals was processed for cardiac histological analysis. Unexposed King's Creek killifish significantly outperformed the unexposed Atlantic Wood killifish in swimming performance as measured by Ucrit (i.e., critical swimming speed). However, King's Creek killifish exposed to Elizabeth River sediment extract (both 0.1 and 1.0%) showed significant declines in Ucrit. Histological analysis revealed the presence of blood in the pericardium of King's Creek killifish. Although Atlantic Wood killifish showed baseline performance deficits relative to King's Creek killifish, their pericardial cavities were nearly free of blood and atrial and ventricular alterations. These findings may explain, in part, the diminished swimming performance of King's Creek fish. Environ Toxicol Chem 2017;36:3246-3253. © 2017 SETAC.

A non-destructive BFCOD assay for in vivo measurement of cytochrome P450 3A (CYP3A) enzyme activity in fish embryos and larvae

There is increasing interest in quantifying the exposure and effects of anthropogenic contaminants in fish. Determination of exposures in wild fish is routinely performed, but methods to investigate potential effects are less established. One of the most relevant approaches would be the use of in vivo assays, but existing assays are often limited to in vitro determination of enzyme activity. Many pharmaceuticals and some persistent pollutants activate, and are metabolized by cytochrome P4503A (CYP3A), which make it a relevant and desirable target for biomarker research. We altered the established 7-benzyloxy-4-trifluoromethylcoumarin-O-debenzylation (BFCOD) in vitro protocol for CYP3A activity determination, developing a rapid and inexpensive method to measure in vivo (and in ovo) CYP3A activity in two fish systems: Gulf killifish (Fundulus grandis) and zebrafish (Danio rerio) early life stages. Even with very low concentrations of 7-benzyloxy-4-trifluoromethyl coumarin (BFC, 0.06 µM or 20 µg/L), we were able to detect significant induction in CYP3A activity in embryos of F. grandis, as well as in larvae of D. rerio in response to benzo[a]pyrene (BaP) and fluoranthene (FL) exposures. Because of concerns regarding the possible contribution of CYP1A to BFCOD activity from previous research, we have used a CYP1A post-translational inhibitor (FL) in order to calculate the contribution of CYP1A to the BFCOD assay. We also dosed with benzo[k]fluoranthene (BkF) and showed significant induction of CYP1A activity, with no concurrent increase in CYP3A activity. In this paper, we have taken an established in vitro CYP3A activity assay, and utilized the reaction in a novel way to allow for the non-destructive determination of CYP3A. In summary, we describe a sensitive, cheap, fast and easy modified BFCOD assay for in ovo and in vivo determination of CYP3A activity for use in moderate throughput early-life-stage fish experiments.

Induced pesticide tolerance results from detoxification pathway priming

Few studies in developmental toxicology have focused on whether early life contaminant exposure affects future susceptibility. Investigations in frogs suggested that early life exposure to a pesticide resulted in higher tolerance to a subsequent challenge. This led to the hypothesis that early-life stage exposures can alter phenotypically plastic traits during development, resulting in induced tolerance. Here, we used Gulf killifish (Fundulus grandis) to test the role of detoxification pathway priming in this inducible tolerance. In frogs, the induced tolerance is present five days after the end of the pre-exposure, but absent after a month. We show that a pre-exposure early in life with carbaryl, induces the activity of cytochrome P450 1A (CYP1A) and increases the ability of pre-exposed groups to metabolize carbaryl, likely because of activation of the aryl hydrocarbon receptor (AHR) pathway. Embryos pre-exposed to carbaryl had a 350-500% increase in CYP1A activity, threefold greater capacity to metabolize carbaryl and were more tolerant to a lethal challenge five days after the end of pre-exposure. However, ten days later the differences in CYP1A activity, metabolic capacity and tolerance between pre-exposed and control groups were no longer present. Thus, we conclude that the increase in tolerance observed in pre-exposed fish embryos was due to the activation of the AHR and other metabolic pathways, resulting in a prolonged increase in biotransformation capacity. This allowed individuals to more efficiently deal with subsequent chemical challenges for a short period after the initial pre-exposure. However, this induced tolerance was only short-lived due to the recycling of biotransformation enzymes in the cells as part of general cellular protein maintenance. These findings suggest that induced tolerance was likely due to induction of defense mechanisms during the duration of response to the original stressor, rather than a more permanent change in their ability to respond to future challenges.

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.

Resistance to polycyclic aromatic hydrocarbon toxicity and associated bioenergetic consequences in a population of Fundulus heteroclitus

Several locations in the Elizabeth River, VA, USA are highly contaminated with polycyclic aromatic hydrocarbons (PAHs) due to the release of creosote mixtures from wood treatment facilities. Interestingly, some populations of Atlantic killifish (Fundulus heteroclitus) inhabiting the Elizabeth River (ER) are resistant to PAH-induced teratogenesis. However, evolutionary resistance to PAHs due to chronic PAH exposure is associated with reduced fitness and increased susceptibility to other environmental stressors in at least one PAH-resistant ER killifish population. More specifically, wild-caught and first generation PAH-resistant juvenile killifish have altered metabolic demands when compared to non-resistant fish. Herein, we investigated this association further by examining a previously under-studied population captured from the creosote-contaminated site Republic Creosoting (Rep). We assessed PAH toxicity and effects on energy metabolism in Rep killifish in comparison with killifish from the reference site Kings Creek (KC). Following exposures to simple and complex PAH mixtures, Rep killifish exhibited several phenotypes associated with PAH resistance including decreased incidences of developmental cardiovascular deformities and recalcitrant cytochrome P450 1A (CYP1A) activity. We evaluated bioenergetics in killifish embryos throughout development and found elevated basal oxygen consumption rates in Rep embryos relative to KC embryos. Furthermore, juvenile F1 Rep fish had significantly lower maximal metabolic rates and aerobic scopes than KC juveniles. These results suggest that populations of killifish that have adapted or evolved to withstand the toxicity associated with PAHs consequently have altered energetic metabolism or demands. Such consequences could result in an enhanced vulnerability to other environmental and anthropogenic stressors in PAH-resistant killifish.

Regulation of Human Cytochrome P4501A1 (hCYP1A1): A Plausible Target for Chemoprevention?

Human cytochrome P450 1A1 (hCYP1A1) has been an object of study due to its role in precarcinogen metabolism; for this reason it is relevant to know more in depth the mechanisms that rule out its expression and activity, which make this enzyme a target for the development of novel chemiopreventive agents. The aim of this work is to review the origin, regulation, and structural and functional characteristics of CYP1A1 letting us understand its role in the bioactivation of precarcinogen and the consequences of its modulation in other physiological processes, as well as guide us in the study of this important protein.

Embryonic cardiotoxicity of weak aryl hydrocarbon receptor agonists and CYP1A inhibitor fluoranthene in the Atlantic killifish (Fundulus heteroclitus)

High affinity aryl hydrocarbon receptor (AHR) ligands, such as certain polychlorinated biphenyls and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), cause severe cardiac teratogenesis in fish embryos. Moderately strong AHR agonists, for example benzo[a]pyrene and β-naphthoflavone, are capable of causing similar cardiotoxic effects, particularly when coupled with cytochrome P450 1A (CYP1A) inhibitors (e.g., fluoranthene (FL). Additionally, some weaker AHR agonists (carbaryl, 2-methylindole, 3-methylindole, and phenanthrene) are known to also cause cardiotoxicity in zebrafish (Danio rerio) embryos when coupled with FL; however, the cardiotoxic effects were not mediated specifically by AHR stimulation. This study was performed to determine if binary exposure to weak AHR agonists and FL were also capable of causing cardiotoxicity in Atlantic killifish Fundulus heteroclitus embryos. Binary exposures were performed in both naïve and PAH-adapted killifish embryos to examine resistance to weak agonists and FL binary exposures. Weak agonists used in this study included the following: carbaryl, phenanthrene, 2-methylindole, 3-methylindole, indigo, and indirubin. Carbaryl, indigo, and indirubin induced the highest CYP1 activity levels in naïve killifish embryos, but no significant CYP1 induction was observed in the PAH-adapted killifish. Embryos were coexposed to subteratogenic levels of each agonist and 500μg/L FL to assess if binary administration could cause cardiotoxicity. Indigo and indirubin coupled with FL caused cardiac teratogenesis in naïve killifish, but coexposures did not produce cardiac chamber abnormalities in the PAH-adapted population. Knockdown of AHR2 in naïve killifish embryos did not prevent cardiac teratogenesis. The data suggest a unique mechanism of cardiotoxicity that is not driven by AHR2 activation.

Silver toxicity across salinity gradients: the role of dissolved silver chloride species (AgCl x ) in Atlantic killifish (Fundulus heteroclitus) and medaka (Oryzias latipes) early life-stage toxicity

The influence of salinity on Ag toxicity was investigated in Atlantic killifish (Fundulus heteroclitus) early life-stages. Embryo mortality was significantly reduced as salinity increased and Ag(+) was converted to AgCl(solid). However, as salinity continued to rise (>5 ‰), toxicity increased to a level at least as high as observed for Ag(+) in deionized water. Rather than correlating with Ag(+), Fundulus embryo toxicity was better explained (R(2) = 0.96) by total dissolved Ag (Ag(+), AgCl2 (-), AgCl3 (2-), AgCl4 (3-)). Complementary experiments were conducted with medaka (Oryzias latipes) embryos to determine if this pattern was consistent among evolutionarily divergent euryhaline species. Contrary to Fundulus data, medaka toxicity data were best explained by Ag(+) concentrations (R(2) = 0.94), suggesting that differing ionoregulatory physiology may drive observed differences. Fundulus larvae were also tested, and toxicity did increase at higher salinities, but did not track predicted silver speciation. Alternatively, toxicity began to increase only at salinities above the isosmotic point, suggesting that shifts in osmoregulatory strategy at higher salinities might be an important factor. Na(+) dysregulation was confirmed as the mechanism of toxicity in Ag-exposed Fundulus larvae at both low and high salinities. While Ag uptake was highest at low salinities for both Fundulus embryos and larvae, uptake was not predictive of toxicity.

Advances in the Study of Heart Development and Disease Using Zebrafish

Animal models of cardiovascular disease are key players in the translational medicine pipeline used to define the conserved genetic and molecular basis of disease. Congenital heart diseases (CHDs) are the most common type of human birth defect and feature structural abnormalities that arise during cardiac development and maturation. The zebrafish, Danio rerio, is a valuable vertebrate model organism, offering advantages over traditional mammalian models. These advantages include the rapid, stereotyped and external development of transparent embryos produced in large numbers from inexpensively housed adults, vast capacity for genetic manipulation, and amenability to high-throughput screening. With the help of modern genetics and a sequenced genome, zebrafish have led to insights in cardiovascular diseases ranging from CHDs to arrhythmia and cardiomyopathy. Here, we discuss the utility of zebrafish as a model system and summarize zebrafish cardiac morphogenesis with emphasis on parallels to human heart diseases. Additionally, we discuss the specific tools and experimental platforms utilized in the zebrafish model including forward screens, functional characterization of candidate genes, and high throughput applications.

Developmental exposure to a complex PAH mixture causes persistent behavioral effects in naive Fundulus heteroclitus (killifish) but not in a population of PAH-adapted killifish

Acute exposures to some individual polycyclic aromatic hydrocarbons (PAHs) and complex PAH mixtures are known to cause cardiac malformations and edema in the developing fish embryo. However, the heart is not the only organ impacted by developmental PAH exposure. The developing brain is also affected, resulting in lasting behavioral dysfunction. While acute exposures to some PAHs are teratogenically lethal in fish, little is known about the later life consequences of early life, lower dose subteratogenic PAH exposures. We sought to determine and characterize the long-term behavioral consequences of subteratogenic developmental PAH mixture exposure in both naive killifish and PAH-adapted killifish using sediment pore water derived from the Atlantic Wood Industries Superfund Site. Killifish offspring were embryonically treated with two low-level PAH mixture dilutions of Elizabeth River sediment extract (ERSE) (TPAH 5.04 μg/L and 50.4 μg/L) at 24h post fertilization. Following exposure, killifish were raised to larval, juvenile, and adult life stages and subjected to a series of behavioral tests including: a locomotor activity test (4 days post-hatch), a sensorimotor response tap/habituation test (3 months post hatch), and a novel tank diving and exploration test (3months post hatch). Killifish were also monitored for survival at 1, 2, and 5 months over 5-month rearing period. Developmental PAH exposure caused short-term as well as persistent behavioral impairments in naive killifish. In contrast, the PAH-adapted killifish did not show behavioral alterations following PAH exposure. PAH mixture exposure caused increased mortality in reference killifish over time; yet, the PAH-adapted killifish, while demonstrating long-term rearing mortality, had no significant changes in mortality associated with ERSE exposure. This study demonstrated that early embryonic exposure to PAH-contaminated sediment pore water caused long-term locomotor and behavioral alterations in killifish, and that locomotor alterations could be observed in early larval stages. Additionally, our study highlights the resistance to behavioral alterations caused by low-level PAH mixture exposure in the adapted killifish population. Furthermore, this is the first longitudinal behavioral study to use killifish, an environmentally important estuarine teleost fish, and this testing framework can be used for future contaminant assessment.

The Elizabeth River Story: A Case Study in Evolutionary Toxicology

The Elizabeth River system is an estuary in southeastern Virginia, surrounded by the towns of Chesapeake, Norfolk, Portsmouth, and Virginia Beach. The river has played important roles in U.S. history and has been the location of various military and industrial activities. These activities have been the source of chemical contamination in this aquatic system. Important industries, until the 1990s, included wood treatment plants that used creosote, an oil-derived product that is rich in polycyclic aromatic hydrocarbons (PAH). These plants left a legacy of PAH pollution in the river, and in particular Atlantic Wood Industries is a designated Superfund site now undergoing remediation. Numerous studies examined the distribution of PAH in the river and impacts on resident fauna. This review focuses on how a small estuarine fish with a limited home range, Fundulus heteroclitus (Atlantic killifish or mummichog), has responded to this pollution. While in certain areas of the river this species has clearly been impacted, as evidenced by elevated rates of liver cancer, some subpopulations, notably the one associated with the Atlantic Wood Industries site, displayed a remarkable ability to resist the marked effects PAH have on the embryonic development of fish. This review provides evidence of how pollutants have acted as evolutionary agents, causing changes in ecosystems potentially lasting longer than the pollutants themselves. Mechanisms underlying this evolved resistance, as well as mechanisms underlying the effects of PAH on embryonic development, are also described. The review concludes with a description of ongoing and promising efforts to restore this historic American river.

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.

Resistance to teratogenesis by F1 and F2 embryos of PAH-adapted Fundulus heteroclitus is strongly inherited despite reduced recalcitrance of the AHR pathway

Atlantic killifish (Fundulus heteroclitus) inhabiting the Atlantic Wood Superfund site on the Elizabeth River (Portsmouth, VA, USA) are exposed to a complex mixture of polycyclic aromatic hydrocarbons (PAHs) from former creosote operations, but are resistant to the acute toxicity and cardiac teratogenesis caused by PAHs. The resistance is associated with a dramatic recalcitrance to induction of cytochrome P450 (CYP1) metabolism enzymes following exposure to aryl hydrocarbon receptor (AHR) agonists, along with an elevated antioxidant response and increased expression of several other xenobiotic metabolism and excretion enzymes. However, the heritability of the resistance in the absence of chemical stressors has been inconsistently demonstrated. Understanding the heritability of this resistance will help clarify the nature of population-level responses to chronic exposure to PAH mixtures and aid in identifying the important mechanistic components of resistance to aryl hydrocarbons. We compared the response of Atlantic Wood F1 and F2 embryos to benzo[k]fluoranthene (BkF), benzo[a]pyrene (BaP), 3,3',4,4',5-pentachlorobiphenyl (PCB-126), and a mixture of BkF and fluoranthene (Fl) to that of F1 embryos of reference site killifish. Resistance to cardiac teratogenesis and induction of CYP mRNA expression and CYP activity was determined. We found that both Atlantic Wood F1 and F2 embryos were highly resistance to cardiac teratogenesis. However, the resistance by Atlantic Wood F2 embryos to induction of CYP mRNA expression and enzyme activity was intermediate between that of Atlantic Wood F1 embryos and reference embryos. These results suggest that resistance to cardiac teratogenesis in Atlantic Wood fish is conferred by multiple factors, not all of which appear to be fully genetically heritable.

AHR2-Mediated transcriptomic responses underlying the synergistic cardiac developmental toxicity of PAHs

Polycyclic aromatic hydrocarbons (PAHs) induce developmental defects including cardiac deformities in fish. The aryl hydrocarbon receptor (AHR) mediates the toxicity of some PAHs. Exposure to a simple PAH mixture during embryo development consisting of an AHR agonist (benzo(a)pyrene-BaP) with fluoranthene (FL), an inhibitor of cytochrome p450 1(CYP1)--a gene induced by AHR activation--results in cardiac deformities. Exposure to BaP or FL alone at similar concentrations alters heart rates, but does not induce morphological deformities. Furthermore, AHR2 knockdown prevents the toxicity of BaP + FL mixture. Here, we used a zebrafish microarray analysis to identify heart-specific transcriptomic changes during early development that might underlie cardiotoxicity of BaP + FL. We used AHR2 morphant embryos to determine the role of this receptor in mediating toxicity. Control and knockdown embryos at 36 h post-fertilization were exposed to DMSO, 100 μg/l BaP, 500 μg/l FL, or 100 μg/l BaP + 500 μg/l FL, and heart tissues for RNA were extracted at 2, 6, 12, and 18 h-post-exposure (hpe), prior to the appearance of cardiac deformities. Data show AHR2-dependent BaP + FL effects on expression of genes involved in protein biosynthesis and neuronal development in addition to signaling molecules and their associated molecular pathways. Ca(2+)-cycling and muscle contraction genes were the most significantly differentially expressed category of transcripts when comparing BaP + FL-treated AHR2 morphant and control embryos. These differences were most prominent at 2 and 6 hpe. Therefore, we postulate that BaP + FL may affect cellular Ca(2+) levels and subsequently cardiac muscle function, potentially underlying BaP + FL cardiotoxicity.

Salinity-dependent copper accumulation in the guppy Poecilia vivipara is associated with CTR1 and ATP7B transcriptional regulation

Copper (Cu) accumulation and regulation of key-genes involved in Cu homeostasis were evaluated in freshwater- and saltwater-acclimated guppies Poecilia vivipara. Fish were exposed (96h) to environmentally relevant concentrations of dissolved Cu (0, 5.0, 9.0 and 20.0μg/L). In freshwater guppies, gill and liver Cu accumulation was dependent on Cu concentration in the exposure medium. In saltwater guppies, this dependence was observed only in the gut. These findings indicate that Cu accumulation was salinity- and tissue-dependent. Key genes involved in Cu metabolism were sequenced for the first time in P. vivipara. Transcripts coding for the high-affinity copper transporter (CTR1) and copper-transporting ATPase (ATP7B) were identified using polymerase chain reaction (PCR) and gene sequencing. The full-length CTR1 open reading frame (1560bp) and a partial ATP7B (690bp) were discovered. Predicted amino acid sequences shared high identities with the CTR1 of Fundulus heteroclitus (81%) and the ATP7B of Sparus aurata (87%). Basal transcriptional levels addressed by RT-qPCR in control fish indicate that CTR1 and ATP7B was highly transcribed in liver of freshwater guppies while CTR1 was highly transcribed in gut of saltwater guppies. This could explain the higher Cu accumulation observed in liver of freshwater guppies and in gut of saltwater guppies, because CTR1 is involved in Cu uptake. Reduced gill mRNA expression of CTR1 was observed in freshwater guppies exposed to 20.0μg/L Cu and in saltwater guppies exposed to 5.0μg/L Cu. In turn, reduced mRNA expression of gut ATP7B was observed in freshwater and salt water guppies exposed to 9.0 and 20.0μg/L Cu. Liver CTR1 and ATP7B transcription were not affected by Cu exposure. These findings suggest that gill CTR1 and gut ATP7B are down-regulated to limit Cu absorption after exposure to dissolved Cu, while liver CTR1 and ATP7B levels are maintained to allow Cu storage and detoxification. In conclusion, findings reported here indicate that Cu accumulation in the euryhaline guppy P. vivipara is tissue specific and dependent on water salinity. They also suggest that Cu homeostasis involves a differential transcriptional regulation of the newly identified Cu transporters, CTR1 and ATP7B.

Evolved resistance to PCB- and PAH-induced cardiac teratogenesis, and reduced CYP1A activity in Gulf killifish (Fundulus grandis) populations from the Houston Ship Channel, Texas

The Houston Ship Channel (HSC), connecting Houston, Texas to Galveston Bay and ultimately the Gulf of Mexico, is heavily industrialized and includes several areas that have historically been identified as containing significant levels of mercury, dioxins, furans, polychlorinated biphenyls (PCBs), and polycyclic aromatic hydrocarbons (PAHs). Gulf killifish, Fundulus grandis, inhabit this entire estuarine system, including the most contaminated areas. F. grandis is the sister species of the well-established estuarine model organism Fundulus heteroclitus, for which heritable resistance to both PCB and PAH toxicity has been documented in several populations. F. grandis collected from two Superfund sites on the HSC and from a reference population were used to establish breeding colonies. F1 embryos from HSC populations were approximately 1000-fold more resistant to PCB126- and 2-5-fold more resistant to coal tar-induced cardiovascular teratogenesis, relative to embryos from the reference population. Reciprocal crosses between reference and contaminated populations exhibit an intermediate level of resistance, confirming that observed protection is genetic and biparentally inherited. Ethoxyresorufin-O-deethylase (EROD) data confirm a reduction in basal and induced cytochrome P4501A (CYP1A) activity in resistant populations of F. grandis. This result is consistent with responses previously described for resistant populations of F. heteroclitus, specifically a recalcitrant aryl hydrocarbon receptor (AHR) pathway. The decreased levels of cardiovascular teratogenesis, and decrease in CYP1A inducibility in response to PCB126 and a PAH mixture, suggest that HSC F. grandis populations have adapted to chronic contaminants exposures via a mechanism similar to that previously described for F. heteroclitus. To the best of our knowledge, this is the first documentation of evolved pollution resistance in F. grandis. Additionally, the mechanistic similarities between the population adaptation observed in this study and previous work in F. heteroclitus suggest that genetic variation predating the evolutionary divergence of these two species may best explain the apparent rapid parallel evolution of pollution resistance in genetically and geographically distinct species and populations.

Compound- and mixture-specific differences in resistance to polycyclic aromatic hydrocarbons and PCB-126 among Fundulus heteroclitus subpopulations throughout the Elizabeth River estuary (Virginia, USA)

Atlantic killifish (Fundulus heteroclitus) inhabiting the Atlantic Wood Industries Superfund Site (Elizabeth River, Portsmouth, VA, USA) are resistant to the acute toxicity and cardiac teratogenesis caused by high levels of polycyclic aromatic hydrocarbons (PAHs) from creosote. The resistance is linked to down regulation of the aryl hydrocarbon receptor (AHR) pathway. We investigated the association between CYP1 activity, as a marker of potential AHR pathway suppression, and contaminant resistance in killifish subpopulations from sites throughout the estuary that varied significantly in PAH contamination level. Adult killifish and sediments were collected from seven sites across approximately 13.7 km in river length within the estuary and from a nearby reference site. Sediment PAH levels were determined using gas chromatography mass spectrometry. Embryos obtained via manual spawning were exposed to individual AHR agonists and PAH mixtures 24 h post fertilization (hpf); CYP1 activity was determined by in ovo ethoxyresorufin-o-deethylase (EROD) at 96 hpf, and cardiac deformity severity was scored at 144 hpf. The total PAH levels measured among the sites varied from approximately 200 to 125,000 ng/g dry sediment. Overall, the resistance to teratogenesis was strongest in the subpopulations from sites in or closest to the major PAH contamination sites, but even embryos from less-contaminated sites within the Elizabeth River demonstrated at least partial resistance to many challenges. Surprisingly, all of the subpopulations tested were highly resistant to PCB-126 (3,3',4,4',5-pentachlorobiphenyl). However, the degree of CYP1 activity response varied significantly among subpopulations and did not always correlate strongly with resistance to teratogenesis; some subpopulations resisted the cardiac teratogenesis caused by the challenges at doses that still elicited strong EROD induction. Our results suggest that there is variation in the adaptive phenotype exhibited by laboratory-spawned embryos from killifish subpopulations throughout the estuary. Furthermore, the results show that contaminants have affected killifish subpopulations throughout the estuary, even in sites with lower levels of PAHs.

A simple and reliable in vivo EROD activity measurement in single Fundulus heteroclitus embryo and larva

Dioxin-like compounds (DLC) induce toxic responses in early life stages of fish through activation of the aryl hydrocarbon receptor (AhR) which is frequently assessed by ethoxyresorufin-O-deethylase (EROD) activity. A novel spectrofluorimetric method was developed to quantitatively assess EROD activity in individual living embryos and prolarvae of a marine model fish species, the mummichog Fundulus heteroclitus. This in vivo method is based on the measurement of the production of resorufin by single live embryos or prolarvae after 5 h incubation with ethoxyresorufin. Freshly fertilized eggs were treated topically from 2.5 to 50 pg egg⁻¹, with 3,3',4,4',5-pentachlorobiphenyl (PCB126), a prototypical DLC. EROD activity was assessed in embryos (7 days post-fertilization) and prolarvae (16 days post-fertilization). Resorufin was measured both in the culture medium (25‰ seawater) and in whole fish homogenates, to assess the percentage retained in the body. Approximately 95% and 17% of the resorufin produced in vivo was retained in embryos and prolarvae respectively. EROD activity in homogenates of embryos and in the culture medium of prolarvae increased linearly with dose. EROD activity measured by the in vivo method was highly correlated to that measured by a traditional in vitro technique using S9 fractions for both embryos and prolarvae. Both in vivo and in vitro EROD activity were higher in prolarvae than in embryos pretreated with PCB126. EROD induction measured in prolarvae by the in vivo and in vitro methods were similar whereas higher induction was measured in vivo than in vitro in embryos. The in vivo method was more sensitive and as reliable as the in vitro technique, and required a lower number of fish (4 compared to 3 pools of 5). This in vivo method is useful to link EROD induction in individual embryos or prolarvae to other organism-level responses. Further studies with other categories of xenobiotics should be performed to assess potential toxic effects on resorufin absorption/excretion processes which could affect in vivo measurements.

Differential sensitivity to pro-oxidant exposure in two populations of killifish (Fundulus heteroclitus)

New Bedford Harbor (MA, U.S.A.; NBH) is a Superfund site inhabited by Atlantic killifish (Fundulus heteroclitus) with altered aryl hydrocarbon receptor (Ahr) signaling, leading to resistance to effects of polychlorinated biphenyls (PCBs) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). The Ahr is a transcription factor that regulates gene expression of many Phase I and II detoxifying enzymes and interacts with Nrf2, a transcription factor that regulates the response to oxidative stress. This study tested the hypothesis that PCB-resistant killifish exhibit altered sensitivity to oxidative stress. Killifish F(1) embryos from NBH and a clean reference site (Scorton Creek, MA, U.S.A.; SC) were exposed to model pro-oxidant and Nrf2-activator, tert-butylhydroquinone (tBHQ). Embryos were exposed at specific embryonic developmental stages (5, 7, and 9 days post fertilization) and toxicity was assessed, using a deformity score, survival, heart rate, and gene expression to compare sensitivity between PCB -resistant and -sensitive (reference) populations. Acute exposure to tBHQ resulted in transient reduction in heart rate in NBH and SC F(1) embryos. However, embryos from NBH were more sensitive to tBHQ, with more frequent and severe deformities, including pericardial edema, tail deformities, small body size, and reduced pigment and erythrocytes. NBH embryos had lower basal expression of antioxidant genes catalase and glutathione-S-transferase alpha (gsta), and upon exposure to tBHQ, exhibited lower levels of expression of catalase, gsta, and superoxide dismutase compared to controls. This result suggests that adaptation to tolerate PCBs has altered the sensitivity of NBH fish to oxidative stress during embryonic development, demonstrating a cost of the PCB resistance adaptation.

Ahr2-dependence of PCB126 effects on the swim bladder in relation to expression of CYP1 and cox-2 genes in developing zebrafish

The teleost swim bladder is assumed a homolog of the tetrapod lung. Both swim bladder and lung are developmental targets of persistent aryl hydrocarbon receptor (AHR(2)) agonists; in zebrafish (Danio rerio) the swim bladder fails to inflate with exposure to 3,3',4,4',5-pentachlorobiphenyl (PCB126). The mechanism for this effect is unknown, but studies have suggested roles of cytochrome P450 1 (CYP1) and cyclooxygenase 2 (Cox-2) in some Ahr-mediated developmental effects in zebrafish. We determined relationships between swim bladder inflation and CYP1 and Cox-2 mRNA expression in PCB126-exposed zebrafish embryos. We also examined effects on β-catenin dependent transcription, histological effects, and Ahr2 dependence of the effect of PCB126 on swim bladder using morpholinos targeting ahr2. One-day-old embryos were exposed to waterborne PCB126 or carrier (DMSO) for 24h and then held in clean water until day 4, a normal time for swim bladder inflation. The effects of PCB126 were concentration-dependent with EC(50) values of 1.4 to 2.0 nM for induction of the CYP1s, 3.7 and 5.1 nM (or higher) for cox-2a and cox-2b induction, and 2.5 nM for inhibition of swim bladder inflation. Histological defects included a compaction of the developing bladder. Ahr2-morpholino treatment rescued the effect of PCB126 (5 nM) on swim bladder inflation and blocked induction of CYP1A, cox-2a, and cox-2b. With 2nM PCB126 approximately 30% of eleutheroembryos(3) failed to inflate the swim bladder, but there was no difference in CYP1 or cox-2 mRNA expression between those embryos and embryos showing inflated swim bladder. Our results indicate that PCB126 blocks swim bladder inflation via an Ahr2-mediated mechanism. This mechanism seems independent of CYP1 or cox-2 mRNA induction but may involve abnormal development of swim bladder cells.

Fundulus heteroclitus adapted to PAHs are cross-resistant to multiple insecticides

Atlantic killifish (Fundulus heteroclitus) from the Atlantic Wood Superfund site on the Elizabeth River (ER), VA are dramatically resistant to the acute toxicity and teratogenesis caused by polycyclic aromatic hydrocarbons (PAHs). To understand the consequences of adaptation to chronic PAH pollution, we have attempted to further define the chemical tolerance associated with this resistance. An important component of the PAH adaptation of ER fish is the dramatic down-regulation of the aryl hydrocarbon receptor (AHR) pathway, resulting in decreased cytochrome p450 (CYP) 1 activity. Herein, we compared the susceptibility to several insecticides of ER fish to that of reference site (King's Creek; KC) fish; use of these chemicals as probes of the resistance will help to demonstrate if the contaminant adaptation exhibited by ER fish is broad or narrow and AHR-focused. We hypothesized that ER fish would be less susceptible to the organophosphate chlorpyrifos (activated by CYP) and more susceptible to the pyrethroid permethrin (detoxified by CYP). Comparison of acute toxicity in 5-day-old larvae supported this hypothesis for chlorpyrifos. As expected, chemical up-regulation of CYP by co-exposure to β-naphthoflavone (BNF) enhanced the susceptibility of KC but it did not affect ER larvae. Unexpectedly, ER larvae were much less susceptible to permethrin than KC larvae. However, co-exposure to BNF greatly decreased the susceptibility of KC larvae, indicating that metabolism of permethrin by CYP was protective. Additionally, fish from each population were compared for susceptibility to the carbamate carbaryl, an acute neurotoxicant and weak AHR agonist that induces teratogenesis similar to that caused by PAHs. ER embryos and larvae were less susceptible than KC fish. These results suggest that the adaptive phenotype of ER fish is multi-faceted and that aspects other than CYP response are likely to greatly affect their response to contaminants.

Biomarkers of organic contamination in the South American fish Poecilia vivipara and Jenynsia multidentata

South American cyprinodontiform fish are potential candidates to be used as model biomarker species of exposure in environmental toxicology. The aim of this study was to identify molecular and biochemical biomarkers of pollution using Poecilia vivipara (Poecilidae) and Jenynsia multidentata (Anablepidae). Partial nucleotide sequences for cytochrome P-450 1A (CYP1A), a classical biomarker of exposure to organic contaminants in fish, were identified in P. vivipara and J. multidentata (approximately 650 nucleotides) using degenerated primers and polymerase chain reaction (PCR). These sequences shared approximately 90% identity in the predicted amino acid sequence with the corresponding CYP1A region of Fundulus heteroclitus. Real-time quantitative PCR (RT-qPCR) analysis confirmed that CYP1A transcription was markedly induced in the liver and gills of J. multidentata (approximately185-fold and 20-fold, respectively) and P. vivipara (122-fold and 739-fold, respectively) 24 h after exposure to 1 μM synthetic CYP1A inducer β-naphthoflavone (BNF). At 24 h after injection with 1 μg/g environmental carcinogenic contaminant benzo[a]pyrene (BaP), a decreased total antioxidant capacity against peroxyl radicals was observed both in liver of J. multidentata and gills of P. vivipara. BaP injection in both fish did not produce changes in lipid peroxide (thiobarbituric acid-reactive substances, TBARS) levels, suggesting an absence of an oxidative stress condition. The newly identified CYP1A may thus serve as general biomarker of exposure to organic contaminant in future studies using P. vivipara and J. multidentata. Data also indicate the importance of species-specific differences in biomarker responses in these South American cyprinodontiform fish, suggesting distinct resistance/susceptibility properties to polycyclic aromatic hydrocarbons.

Morpholino gene knockdown in adult Fundulus heteroclitus: role of SGK1 in seawater acclimation

The Atlantic killifish (Fundulus heteroclitus) is an environmental sentinel organism used extensively for studies on environmental toxicants and salt (NaCl) homeostasis. Previous research in our laboratory has shown that rapid acclimation of killifish to seawater is mediated by trafficking of CFTR chloride channels from intracellular vesicles to the plasma membrane in the opercular membrane within the first hour in seawater, which enhances chloride secretion into seawater, thereby contributing to salt homeostasis. Acute transition to seawater is also marked by an increase in both mRNA and protein levels of serum glucocorticoid kinase 1 (SGK1) within 15 minutes of transfer. Although the rise in SGK1 in gill and its functional analog, the opercular membrane, after seawater transfer precedes the increase in membrane CFTR, a direct role of SGK1 in elevating membrane CFTR has not been established in vivo. To test the hypothesis that SGK1 mediates the increase in plasma membrane CFTR we designed two functionally different vivo-morpholinos to knock down SGK1 in gill, and developed and validated a vivo-morpholino knock down technique for adult killifish. Injection (intraperitoneal, IP) of the splice blocking SGK1 vivo-morpholino reduced SGK1 mRNA in the gill after transition from fresh to seawater by 66%. The IP injection of the translational blocking and splice blocking vivo-morpholinos reduced gill SGK1 protein abundance in fish transferred from fresh to seawater by 64% and 53%, respectively. Moreover, knock down of SGK1 completely eliminated the seawater induced rise in plasma membrane CFTR, demonstrating that the increase in SGK1 protein is required for the trafficking of CFTR from intracellular vesicles in mitochondrion rich cells to the plasma membrane in the gill during acclimation to seawater. This is the first report of the use of vivo-morpholinos in adult killifish and demonstrates that vivo-morpholinos are a valuable genetic tool for this environmentally relevant model organism.

Mechanistic research in aquatic toxicology: perspectives and future directions

On the 30th anniversary of the journal, I provide a perspective on some of the questions and opportunities for new understanding that will interest aquatic toxicologists during the next 30 years. I focus on mechanisms of toxicity involving transcription factors, signalling pathways, and gene networks involved in toxic and adaptive responses in aquatic animals. Prominent questions address the value of a toxicity pathways approach in aquatic systems, issues involving extrapolation among species, identification of susceptibility genes and useful biomarkers of adverse effect, new emerging contaminants, the importance of epigenetic mechanisms, effects of multiple stressors, evolutionary toxicology, and the relative roles of technical and conceptual limitations to our understanding of chemical effects on aquatic systems.

AHR2 knockdown prevents PAH-mediated cardiac toxicity and XRE- and ARE-associated gene induction in zebrafish (Danio rerio)

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants often present in aquatic systems as complex mixtures. Embryonic fish are sensitive to the developmental toxicity of some PAHs, but the exact mechanisms involved in this toxicity are still unknown. This study explored the role of the aryl hydrocarbon receptor (AHR) in the oxidative stress response of zebrafish to the embryotoxicity of select PAHs. Embryos were exposed to two PAHs, benzo[k]fluoranthene (BkF; a strong AHR agonist) and fluoranthene (FL; a cytochrome P4501A (CYP1A) inhibitor), alone and in combination. CYP1A, CYP1B1, CYP1C1, and redox-responsive genes glutathione s-transferase pi 2 (GSTp2), glutathione peroxidase 1 (GPx1), the glutamate-cysteine ligase catalytic subunit (GCLc), MnSOD and CuZnSOD mRNA expression was examined. CYP1 activity was measured via an in vivo ethoxyresorufin-O-deethlyase (EROD) activity assay, and the area of the pericardium was measured as an index of cardiotoxicity. BkF or FL alone caused no deformities whereas BkF+FL resulted in extreme pericardial effusion. BkF induced CYP activity above controls and co-exposure with FL inhibited this activity. BkF induced expression of all three CYPs, GSTp2, and GCLc. BkF+FL caused greater than additive induction of the three CYPs, GSTp2, GPx1, and GCLc but had no effect on MnSOD or CuZnSOD. AHR2 knockdown protected against the cardiac deformities caused by BkF+FL and significantly inhibited the induction of the CYPs, GSTp2, GPx1, and GCLc after BkF+FL compared to non-injected controls. These results further show the protective role of AHR2 knockdown against cardiotoxic PAHs and the role of AHR2 as a mediator of redox-responsive gene induction.

Transcriptomic assessment of resistance to effects of an aryl hydrocarbon receptor (AHR) agonist in embryos of Atlantic killifish (Fundulus heteroclitus) from a marine Superfund site

Background

Populations of Atlantic killifish (Fundulus heteroclitus) have evolved resistance to the embryotoxic effects of polychlorinated biphenyls (PCBs) and other halogenated and nonhalogenated aromatic hydrocarbons that act through an aryl hydrocarbon receptor (AHR)-dependent signaling pathway. The resistance is accompanied by reduced sensitivity to induction of cytochrome P450 1A (CYP1A), a widely used biomarker of aromatic hydrocarbon exposure and effect, but whether the reduced sensitivity is specific to CYP1A or reflects a genome-wide reduction in responsiveness to all AHR-mediated changes in gene expression is unknown. We compared gene expression profiles and the response to 3,3',4,4',5-pentachlorobiphenyl (PCB-126) exposure in embryos (5 and 10 dpf) and larvae (15 dpf) from F. heteroclitus populations inhabiting the New Bedford Harbor, Massachusetts (NBH) Superfund site (PCB-resistant) and a reference site, Scorton Creek, Massachusetts (SC; PCB-sensitive).

Results

Analysis using a 7,000-gene cDNA array revealed striking differences in responsiveness to PCB-126 between the populations; the differences occur at all three stages examined. There was a sizeable set of PCB-responsive genes in the sensitive SC population, a much smaller set of PCB-responsive genes in NBH fish, and few similarities in PCB-responsive genes between the two populations. Most of the array results were confirmed, and additional PCB-regulated genes identified, by RNA-Seq (deep pyrosequencing).

Conclusions

The results suggest that NBH fish possess a gene regulatory defect that is not specific to one target gene such as CYP1A but rather lies in a regulatory pathway that controls the transcriptional response of multiple genes to PCB exposure. The results are consistent with genome-wide disruption of AHR-dependent signaling in NBH fish.

Ultraviolet treatment and biodegradation of dibenzothiophene: Identification and toxicity of products

Biodegradation of pollutants often results in incomplete mineralization and formation of degradation products with unknown chemical and toxicological characteristics. Ultraviolet (UV) irradiation, a common technology used in water and wastewater treatment, may help reduce aqueous concentrations of degradation products produced during biological treatment and their associated hazards. Combined biological and UV transformations may be important in natural systems as well. We investigated the effects of UV irradiation (254 nm) on dibenzothiophene (DBT), a sulfur-containing polyaromatic hydrocarbon, in artificial seawater, and its effects on biodegradation products produced from mixed-community microbial transformations of DBT, including DBT sulfone, DBT sulfoxide, hydroxylated and carboxylated benzothiophenes, thiosalicylic acid, and others. Toxicity of solutions after UV exposure was monitored using bioluminescent bacteria (Vibrio fischeri) and by evaluating cardiac deformities in Fundulus heteroclitus embryos. The highest UV fluence reduced DBT concentration by 28% when DBT was present as the sole organic solute. In postbiodegradation solution, the same fluence reduced the initial concentration of DBT by 81%, and 11 DBT biodegradation products to trace levels. Regardless of UV fluence, DBT by itself produced minimal effects in Fundulus embryos but was moderately toxic to V. fischeri. Postbiodegradation solutions were highly toxic to both test organisms. The highest UV fluence slightly reduced toxicity of postbiodegradation solution to V. fischeri but exacerbated cardiac deformities in Fundulus embryos. Toxicity could not be attributed to specific products and was likely a result of mixture effects. These results emphasize that toxicity can increase during remediation and that multiple assays may be necessary for evaluation. The novel approach of combined biodegradation/UV treatment is promising, although further research is needed to reduce toxicity in the case of DBT.

AHR2 mediates cardiac teratogenesis of polycyclic aromatic hydrocarbons and PCB-126 in Atlantic killifish (Fundulus heteroclitus)

Exposure of developing fish to polycyclic aromatic hydrocarbons (PAHs) and halogenated aromatic hydrocarbons (HAHs) results in a suite of defects including cardiac malformation, pericardial and yolk sac edema, craniofacial defects, and hemorrhaging. Several populations of Atlantic killifish or mummichog (Fundulus heteroclitus) on the Atlantic coast of the United States are resistant to the developmental and acute toxicity caused by PAHs and HAHs; this has made Fundulus a valuable model for studying aryl hydrocarbon sensitivity and adaptation. In order to further increase the utility of Fundulus, better understanding of the components of the molecular pathways governing aryl hydrocarbon response in Fundulus is required. The aryl hydrocarbon receptor (AHR) is known to mediate many of the toxic responses to PAHs and HAHs. A single AHR has been identified in mammals, but Fundulus has two AHRs and their relative roles are not clear. In the current study, translation-blocking and splice-junction morpholino gene knockdown was used to determine the roles of AHR1 and AHR2 in mediating cardiac teratogenesis induced by beta-naphthoflavone (BNF), benzo[k]fluoranthene (BkF), and 3,3',4,4',5-pentachlorobiphenyl (PCB-126). Here we report that AHR2 and not AHR1 knockdown resulted in rescue of teratogenicity induced by BNF, BkF, and PCB-126. These data demonstrate that AHR2 is the primary mediator of cardiac teratogenesis caused by multiple aryl hydrocarbons in Fundulus and suggest that suppression of the AHR pathway through modulation of AHR2 is a plausible mechanism for PAH resistance in adapted fish. Additionally, this is the first reported use of splice-junction morpholinos in Fundulus.

Characterization of the recalcitrant CYP1 phenotype found in Atlantic killifish (Fundulus heteroclitus) inhabiting a Superfund site on the Elizabeth River, VA

Fundulus heteroclitus (Atlantic killifish) found at the Atlantic Wood Industries Superfund site on the Elizabeth River (ER) in Portsmouth, VA (USA), have been shown to be resistant to the teratogenic effects of creosote-contaminated sediments found at this highly contaminated site. Many of the polycyclic aromatic hydrocarbons (PAHs) found at the ER are known to activate the aryl hydrocarbon receptor (AHR), and are thought to mediate their toxic effects through this pathway. Activation of the AHR results in the induction of several Phase I and II metabolic enzymes. It has been previously shown that the AHR of killifish from the ER are refractory to induction by AHR agonists. To more fully characterize this altered AHR response, we exposed embryos from the ER and from a reference site on King's Creek, VA (KC) to two PAHs, benzo[alpha]pyrene (BaP) and benzo[k]fluoranthene (BkF), and to the dioxin-like compound (DLC), 3,3',4,4',5-pentachlorobiphenyl (PCB126). We compared their developmental and molecular responses by screening the embryos for CYP1A enzyme activity, cardiac deformities, and mRNA expression of CYP1A, CYP1B1, CYP1C1, and AHR2. Basal gene expression of both CYP1A and CYP1B1 was 40% higher in the KC control embryos compared to those from the ER, while AHR2 and CYP1C1 were not significantly different between the populations. Exposure of KC embryos to BaP, BkF, and PCB126 induced CYP1A activity and cardiac deformities. In contrast, CYP1A activity was induced in ER embryos only in response to BkF exposure, although this induction in ER embryos was significantly lower than that observed in KC fish at comparable concentrations. ER embryos did not develop cardiac deformities in response to any of the chemicals tested. CYP1A, CYP1B1 and CYP1C1 mRNA were all significantly induced in the KC embryos after exposure to BaP, BkF and PCB126. Exposure to BaP and BkF in ER embryos resulted in a significant induction of CYP1A mRNA, albeit significantly lower than observed in KC fish. Interestingly, BaP exposure resulted in induction of CYP1B1 at comparable levels in embryos from both populations. CYP1s were not induced in ER embryos in response to PCB126, nor was CYP1C1 for any treatment examined. Additionally, AHR2 was not significantly induced for any of the treatment groups. This study further characterizes the AHR response in killifish, and provides greater insight into the adapted ER phenotype. The ER adaptation involves the suppression of normal AHR-inducible gene expression for all three CYP1 genes, and therefore is likely an alteration in AHR signaling or control.

Identification of mitochondrial cytochrome P450 induced in response to polycyclic aromatic hydrocarbons in the mummichog (Fundulus heteroclitus)

Increasing evidence suggests that polycyclic aromatic hydrocarbons (PAHs) such as benzo[a]pyrene (BaP) are localized to the mitochondria. Because the toxic effects of many PAHs are the result of metabolism by cytochrome P4501A (CYP1A), it is important to investigate whether active forms of these enzymes can be identified in the mitochondria. In this study, we identified mitochondrial P450s with a monoclonal antibody against scup (Stenotomus chrysops) CYP1A in the isolated mitochondrial fraction of the liver from adult male mummichog (Fundulus heteroclitus) livers. The size of the protein in the mitochondria was similar to that of microsomal CYP1A. Fish dosed with 10mg/kg BaP had increased EROD activity in the mitochondrial fraction compared to controls. In mummichog larvae dosed with 100 microg/L BaP and 100 microg/L benzo[k]fluoranthene, CYP1A protein levels as well as enzyme activity were elevated. However, fish from a PAH-polluted Superfund site (Elizabeth River, Portsmouth VA) showed recalcitrant mitochondrial CYP1A protein levels and enzyme activity in a similar manner to microsomal CYP1A.

Effects of benzo[a]pyrene on mitochondrial and nuclear DNA damage in Atlantic killifish (Fundulus heteroclitus) from a creosote-contaminated and reference site

Benzo[a]pyrene (BaP) is a known genotoxicant that affects both mitochondrial and nuclear DNA (mtDNA, nDNA). Here, we examined mtDNA and nDNA damage in the Atlantic killifish (Fundulus heteroclitus) from a highly contaminated Superfund site (Elizabeth River, VA, USA) and from a reference site (King's Creek, VA, USA) that were dosed with 10 mg/kg BaP. Using the long amplicon quantitative PCR technique, we observed similar increases in mitochondrial and nuclear DNA damage in King's Creek fish treated with BaP. Killifish from the Elizabeth River showed high levels of basal nDNA and mtDNA damage compared to fish from the reference site, but the level of damage induced due to BaP treatment was much lower in Elizabeth River killifish compared to King's Creek fish. Laboratory-reared offspring from both populations showed increased BaP-induced damage in mtDNA, relative to nDNA. Similar to the adult experiment, the Elizabeth River larvae had higher levels of basal DNA damage than those from the reference site, but were less impacted by BaP exposure. Measurements of oxidative DNA damage (8-oxo-deoxyguanine by LC-MS/MS) showed no differences among treatment groups, suggesting that the majority of DNA damage is from covalent binding of BaP metabolites to DNA. This study shows for the first time that mitochondria can be an important target of BaP toxicity in fish, indicating that BaP exposures could have important energetic consequences. Results also suggest that multi-generational exposures in the wild may lead to adaptations that dampen DNA damage arising from BaP exposure.

New cytochrome P450 1B1, 1C2 and 1D1 genes in the killifish Fundulus heteroclitus: Basal expression and response of five killifish CYP1s to the AHR agonist PCB126

Knowledge of the complement of cytochrome P450 (CYP) genes is essential to understanding detoxification and bioactivation mechanisms for organic contaminants. We cloned three new CYP1 genes, CYP1B1, CYP1C2 and CYP1D1, from the killifish Fundulus heteroclitus, an important model in environmental toxicology. Expression of the new CYP1s along with previously known CYP1A and CYP1C1 was measured by qPCR in eight different organs. Organ distribution was similar for the two CYP1Cs, but otherwise patterns and extent of expression differed among the genes. The AHR agonist 3,3',4,4',5-pentachlorobiphenyl (PCB126) (31 pmol/g fish) induced expression of CYP1A and CYP1B1 in all organs examined, while CYP1C1 was induced in all organs except testis. The largest changes in response to PCB126 were induction of CYP1A in testis (approximately 700-fold) and induction of CYP1C1 in liver (approximately 500-fold). CYP1B1 in liver and gut, CYP1A in brain and CYP1C1 in gill also were induced strongly by PCB126 (> 100-fold). CYP1C1 expression levels were higher than CYP1C2 in almost all tissues and CYP1C2 was much less responsive to PCB126. In contrast to the other genes, CYP1D1 was not induced by PCB126 in any of the organs. The organ-specific response of CYP1s to PCB126 implies differential involvement in effects of halogenated aromatic hydrocarbons in different organs. The suite of inducible CYP1s could enhance the use of F. heteroclitus in assessing aquatic contamination by AHR agonists. Determining basal and induced levels of protein and the substrate specificity for all five CYP1s will be necessary to better understand their roles in chemical effects and physiology.

Effect of CYP1A inhibition on the biotransformation of benzo[a]pyrene in two populations of Fundulus heteroclitus with different exposure histories

Sediment contaminated with polycyclic aromatic hydrocarbons (PAHs) from a Superfund site on the Elizabeth River (ER) in Portsmouth, VA, is teratogenic to embryonic killifish (Fundulus heteroclitus) from reference sites. However, embryos born to a population of ER killifish are resistant to PAH-induced teratogenicity. Mechanisms underlying the resistance are unclear; however, ER killifish are refractory to induction of metabolic enzyme cytochrome P4501A (CYP1A), at the level of mRNA, protein and activity. The contaminated ER sediment comprises a complex mixture of PAHs with different mechanisms of toxicity. While many are inducers of metabolic enzymes involved in both phase I and phase II of biotransformation, some PAHs can also inhibit phase I enzymatic activity. Previous research has shown that co-exposure to PAHs with different modes of action can result in synergistic embryotoxicity (Billiard, S.M., Meyer, J.N.D., Wassenberg, M., Hodson, P.V., Di Giulio, R.T., 2008. Nonadditive effects of PAHs on early vertebrate development: mechanisms and implications for risk assessment. Toxicological Sciences 105, 5-23). Two of the abundant PAHs at the ER are fluoranthene (FL), a CYP1A inhibitor, and benzo[a]pyrene (BaP), a CYP1A inducer. Based on the ER resistant phenotype and the PAH mixture in the ER sediment, we hypothesized that the inhibition of CYP1A activity affects the teratogenicity of PAHs through a biotransformation-mediated mechanism. To examine this hypothesis, we compared the responses of killifish embryos born to parents from the ER and from a reference site (King's Creek (KC), VA) after a water-borne exposure to BaP (0-400 microg/L) in the presence or absence of FL (0-500 microg/L). Embryos were dosed from 24 to 120 h post-fertilization (hpf) and were analyzed for induction of CYP1 enzymatic activity as measured by the ethoxyresorufin-O-deethylase (EROD) assay, cardiac deformities, and BaP metabolic profile. KC embryos showed significant induction of CYP1 protein activity at all BaP concentrations examined. Co-exposure to 500 microg/L of FL significantly decreased CYP1 activity and increased cardiac deformities. ER embryos showed no change in CYP1 activity or cardiac deformities for any treatment. Significantly greater concentrations of BaP and BaP 9,10-dihydrodiol were recovered from ER embryos compared to those from KC. Co-exposure with FL did not significantly alter the amount of BaP or the metabolites recovered in either population. These findings suggest that the teratogenicity observed by co-exposure to BaP and FL cannot fully be explained by alteration in BaP metabolism. This study also indicates that the metabolic adaptation observed in the ER killifish cannot be explained simply by the refractory CYP1 phenotype.

Fluoranthene, but not benzo[a]pyrene, interacts with hypoxia resulting in pericardial effusion and lordosis in developing zebrafish

Previous research has documented several PAHs that interact synergistically, causing severe teratogenicity in developing fish embryos. The coexposure of CYP1A inhibitors (e.g. FL or ANF) with AHR agonists (e.g. BaP or BNF) results in a synergistic increase in toxicity. As with chemical CYP1A inhibitors, it has also been shown that CYP1A morpholinos exacerbate BNF-induced embryotoxicity. We hypothesized that a hypoxia-induced reduction in CYP1A activity in BNF or BaP-exposed zebrafish embryos would similarly enhance pericardial effusion and other developmental abnormalities. BaP, BNF, ANF, and FL exposures, both individually and as BaP+FL or BNF+ANF combinations, were performed under hypoxia and normoxia. CYP1A activity in the BaP+hypoxia and BNF+hypoxia embryos was reduced by approximately 60% relative to normoxia embryos. Although CYP1A activity was significantly reduced, we did not observe any increase in pericardial effusion in either group. An unexpected yet particularly interesting result of these experiments was the observed interaction of both FL and ANF with hypoxia. Relatively high, yet environmentally relevant concentrations of FL (100-500 microg L(-1)) interact with moderate hypoxia (7.3% DO) through an unknown mechanism, resulting in pericardial effusion and severe lordosis. Additionally, ANF exposures (100 microg L(-1)) which are not normally teratogenic caused dramatic pericardial effusion, but not lordosis, when embryos were coexposed to hypoxia. These results suggest that reduced CYP1A activity may not exclusively underlie observed developmental toxicity, and that hypoxia may exacerbate the developmental toxicity of some PAH mixtures.