Search

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., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Tributyltin (TBT) and dioxin-like polychlorinated biphenyls (PCBs) are environmental contaminants that are highly toxic to fish and co-occur in New Bedford Harbor (NBH), an estuarine Superfund site located in Massachusetts, USA. Atlantic killifish (Fundulus heteroclitus) that reside in NBH (and other highly contaminated sites along the east coast of the United States) have developed resistance to activation of the aryl hydrocarbon receptor (AHR) pathway and the toxicity of dioxin-like chemicals, such as 3,3',4,4',5-pentachlorobiphenyl, PCB126. In many biological systems, TBT disregulates adipose and bone development via the PPARγ-RXR pathway; AHR activation also disrupts adipose and bone homeostasis, potentially through molecular crosstalk between AHR and PPARγ. However, little is known about how co-exposure and the interaction of these pathways modulate the toxicological effects of these contaminants. Here, we tested the hypotheses that TBT would induce teratogenesis in killifish via activation of PPARγ and that PCB126 co-exposure would suppress PPARγ pathway activation in PCB-sensitive killifish from a reference site (Scorton Creek, SC, PCB-sensitive) but not in PCB-tolerant NBH killifish. Killifish embryos from both populations exposed to TBT (50 and 100 nM) displayed caudal fin deformities. TBT did not change the expression of pparg or its target genes related to adipogenesis (fabp11a and fabp1b) in either population. However, expression of osx/sp7, an osteoblast marker gene, and col2a1b, a chondroblast marker gene, was significantly suppressed by TBT only in SC killifish. An RXR-specific agonist, but not a PPARγ-specific agonist, induced caudal fin deformities like those observed in TBT-treated embryos. PCB126 did not induce caudal fin deformities and did not exacerbate TBT-induced fin deformities. Further, PCB126 increased expression of pparg in SC embryos and not NBH embryos, but did not change the expression of fabp1b. Taken together, these results suggest that in killifish embryos the PPARγ pathway is regulated in part by AHR, but is minimally active at least in this early life stage. In killifish, RXR activation, rather than PPARγ activation, appears to be the mechanism by which TBT induces caudal fin teratogenicity, which is not modulated by AHR responsiveness., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Sentinel species such as the Atlantic killifish (Fundulus heteroclitus) living in urban waterways can be used as toxicological models to understand impacts of environmental metabolism disrupting compound (MDC) exposure on both wildlife and humans. Exposure to MDCs is associated with increased risk of metabolic syndrome, including impaired lipid and glucose homeostasis, adipogenesis, appetite control, and basal metabolism. MDCs are ubiquitous in the environment, including in aquatic environments. New Bedford Harbor (NBH), Massachusetts is polluted with polychlorinated biphenyls (PCBs), and, as we show for the first time, tin (Sn). PCBs and organotins are ligands for two receptor systems known to regulate lipid homeostasis, the aryl hydrocarbon receptor (AHR) and the peroxisome proliferator-activated receptors (PPARs), respectively. In the current study, we compared lipid homeostasis in laboratory-reared killifish from NBH (F2) and a reference location (Scorton Creek, Massachusetts; F1 and F2) to evaluate how adaptation to local conditions may influence responses to MDCs. Adult killifish from each population were exposed to 3,3',4,4',5-pentachlorobiphenyl (PCB126, dioxin-like), 2,2',4,4',5,5'-hexachlorobiphenyl (PCB153, non-dioxin-like), or tributyltin (TBT, a PPARγ ligand) by a single intraperitoneal injection and analyzed after 3 days. AHR activation was assessed by measuring cyp1a mRNA expression. Lipid homeostasis was evaluated phenotypically by measuring liver triglycerides and organosomatic indices, and at the molecular level by measuring the mRNA expression of pparg and ppara and a target gene for each receptor. Acute MDC exposure did not affect phenotypic outcomes. However, overall NBH killifish had higher liver triglycerides and adiposomatic indices than SC killifish. Both season and population were significant predictors of the lipid phenotype. Acute MDC exposure altered hepatic gene expression only in male killifish from SC. PCB126 exposure induced cyp1a and pparg, whereas PCB153 exposure induced ppara. TBT exposure did not induce ppar-dependent pathways. Comparison of lipid homeostasis in two killifish populations extends our understanding of how MDCs act on fish and provides a basis to infer adaptive benefits of these differences in the wild., (Copyright © 2019 Elsevier B.V. All rights reserved.)

There is growing evidence that environmental toxicants can affect various physiological processes by altering DNA methylation patterns. However, very little is known about the impact of toxicant-induced DNA methylation changes on gene expression patterns. The objective of this study was to determine the genome-wide changes in DNA methylation concomitant with altered gene expression patterns in response to 3, 3', 4, 4', 5-pentachlorobiphenyl (PCB126) exposure. We used PCB126 as a model environmental chemical because the mechanism of action is well-characterized, involving activation of aryl hydrocarbon receptor, a ligand-activated transcription factor. Adult zebrafish were exposed to 10 nM PCB126 for 24 h (water-borne exposure) and brain and liver tissues were sampled at 7 days post-exposure in order to capture both primary and secondary changes in DNA methylation and gene expression. We used enhanced Reduced Representation Bisulfite Sequencing and RNAseq to quantify DNA methylation and gene expression, respectively. Enhanced reduced representation bisulfite sequencing analysis revealed 573 and 481 differentially methylated regions in the liver and brain, respectively. Most of the differentially methylated regions are located more than 10 kilobases upstream of transcriptional start sites of the nearest neighboring genes. Gene Ontology analysis of these genes showed that they belong to diverse physiological pathways including development, metabolic processes and regeneration. RNAseq results revealed differential expression of genes related to xenobiotic metabolism, oxidative stress and energy metabolism in response to polychlorinated biphenyl exposure. There was very little correlation between differentially methylated regions and differentially expressed genes suggesting that the relationship between methylation and gene expression is dynamic and complex, involving multiple layers of regulation.

Non-dioxin-like polychlorinated biphenyls (NDL PCBs) activate ryanodine receptors (RyR), microsomal Ca 2+ channels of broad significance. Teleost fish may be important models for NDL PCB neurotoxicity, and we used sequencing databases to characterize teleost RyR and FK506 binding protein 12 or 12.6kDa (genes FKBP1A; FKBP1B), which promote NDL PCB-triggered Ca 2+ dysregulation. Particular focus was placed on describing genes in the Atlantic killifish (Fundulus heteroclitus) genome and searching available RNA-sequencing datasets for single nucleotide variants (SNV) between PCB tolerant killifish from New Bedford Harbor (NBH) versus sensitive killifish from Scorton Creek (SC), MA. Consistent with the teleost whole genome duplication (tWGD), killifish have six RyR genes, corresponding to a and b paralogs of mammalian RyR1, 2 and 3. The presence of six RyR genes was consistent in all teleosts investigated including zebrafish. Killifish have four FKBP1; one FKBP1b and three FKBP1a named FKBP1aa, FKBP1ab, likely from the tWGD and a single gene duplicate FKBP1a3 suggested to have arisen in Atherinomorphae. The RyR and FKBP1 genes displayed tissue and developmental stage-specific mRNA expression, and the previously uncharacterized RyR3, herein named RyR3b, and all FKBP1 genes were prominent in brain. We identified a SNV in RyR3b encoding missense mutation E1458D. In NBH killifish, 57% were heterozygous and 28% were homozygous for this SNV, whereas almost all SC killifish (94%) lacked the variant (n≥39 per population). The outlined sequence differences between mammalian and teleost RyR and FKBP1 together with outlined population differences in SNV frequency may contribute to our understanding of NDL PCB neurotoxicity., (Copyright © 2017 Elsevier B.V. All rights reserved.)

The hypoxia-inducible factor (HIF) family of transcription factors plays central roles in the development, physiology, pathology, and environmental adaptation of animals. Because many aquatic habitats are characterized by episodes of low dissolved oxygen, fish represent ideal models to study the roles of HIF in the response to aquatic hypoxia. The estuarine fish Fundulus heteroclitus is found in habitats prone to hypoxia. It responds to low oxygen via behavioral, physiological, and molecular changes, and one member of the HIF family, HIF2α, has been previously described. Herein, cDNA sequencing, phylogenetic analyses, and genomic approaches were used to determine other members of the HIFα family from F. heteroclitus and their relationships to HIFα subunits from other vertebrates. In vitro and cellular approaches demonstrated that full-length forms of HIF1α, HIF2α, and HIF3α independently formed complexes with the β-subunit, aryl hydrocarbon receptor nuclear translocator, to bind to hypoxia response elements and activate reporter gene expression. Quantitative PCR showed that HIFα mRNA abundance varied among organs of normoxic fish in an isoform-specific fashion. Analysis of the F. heteroclitus genome revealed a locus encoding a second HIF2α-HIF2αb-a predicted protein lacking oxygen sensing and transactivation domains. Finally, sequence analyses demonstrated polymorphism in the coding sequence of each F. heteroclitus HIFα subunit, suggesting that genetic variation in these transcription factors may play a role in the variation in hypoxia responses among individuals or populations., (Copyright © 2017 the American Physiological Society.)

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., (Copyright © 2014 Elsevier B.V. All rights reserved.)

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., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Oxidative stress is an important mechanism of chemical toxicity, contributing to teratogenesis and to cardiovascular and neurodegenerative diseases. Developing animals may be especially sensitive to chemicals causing oxidative stress. The developmental expression and inducibility of anti-oxidant defenses through activation of NF-E2-related factor 2 (NRF2) affect susceptibility to oxidants, but the embryonic response to oxidants is not well understood. To assess the response to chemically mediated oxidative stress and how it may vary during development, zebrafish embryos, eleutheroembryos, or larvae at 1, 2, 3, 4, 5, and 6 days post fertilization (dpf) were exposed to DMSO (0.1%), tert-butylhydroquinone (tBHQ; 10 µM) or 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD; 2 nM) for 6 hr. Transcript abundance was assessed by real-time qRT-PCR and microarray. qRT-PCR showed strong (4- to 5-fold) induction of gstp1 by tBHQ as early as 1 dpf. tBHQ also induced gclc (2 dpf), but not sod1, nqo1, or cyp1a. TCDD induced cyp1a but none of the other genes. Microarray analysis showed that 1477 probes were significantly different among the DMSO-, tBHQ-, and TCDD-treated eleutheroembryos at 4 dpf. There was substantial overlap between genes induced in developing zebrafish and a set of marker genes induced by oxidative stress in mammals. Genes induced by tBHQ in 4-dpf zebrafish included those involved in glutathione synthesis and utilization, signal transduction, and DNA damage/stress response. The strong induction of hsp70 determined by microarray was confirmed by qRT-PCR and by use of transgenic zebrafish expressing enhanced green fluorescent protein (EGFP) under control of the hsp70 promoter. Genes strongly down-regulated by tBHQ included mitfa, providing a molecular explanation for the loss of pigmentation in tBHQ-exposed embryos. These data show that zebrafish embryos are responsive to oxidative stress as early as 1 dpf, that responsiveness varies with development in a gene-specific manner, and that the oxidative stress response is substantially conserved in vertebrate animals.

Background: The non-migratory killifish Fundulus heteroclitus inhabits clean and polluted environments interspersed throughout its range along the Atlantic coast of North America. Several populations of this species have successfully adapted to environments contaminated with toxic aromatic hydrocarbon pollutants such as polychlorinated biphenyls (PCBs). Previous studies suggest that the mechanism of resistance to these and other "dioxin-like compounds" (DLCs) may involve reduced signaling through the aryl hydrocarbon receptor (AHR) pathway. Here we investigated gene diversity and evidence for positive selection at three AHR-related loci (AHR1, AHR2, AHRR) in F. heteroclitus by comparing alleles from seven locations ranging over 600 km along the northeastern US, including extremely polluted and reference estuaries, with a focus on New Bedford Harbor (MA, USA), a PCB Superfund site, and nearby reference sites., Results: We identified 98 single nucleotide polymorphisms within three AHR-related loci among all populations, including synonymous and nonsynonymous substitutions. Haplotype distributions were spatially segregated and F-statistics suggested strong population genetic structure at these loci, consistent with previous studies showing strong population genetic structure at other F. heteroclitus loci. Genetic diversity at these three loci was not significantly different in contaminated sites as compared to reference sites. However, for AHR2 the New Bedford Harbor population had significant FST values in comparison to the nearest reference populations. Tests for positive selection revealed ten nonsynonymous polymorphisms in AHR1 and four in AHR2. Four nonsynonymous SNPs in AHR1 and three in AHR2 showed large differences in base frequency between New Bedford Harbor and its reference site. Tests for isolation-by-distance revealed evidence for non-neutral change at the AHR2 locus., Conclusion: Together, these data suggest that F. heteroclitus populations in reference and polluted sites have similar genetic diversity, providing no evidence for strong genetic bottlenecks for populations in polluted locations. However, the data provide evidence for genetic differentiation among sites, selection at specific nucleotides in AHR1 and AHR2, and specific AHR2 SNPs and haplotypes that are associated with the PCB-resistant phenotype in the New Bedford Harbor population. The results suggest that AHRs, and especially AHR2, may be important, recurring targets for selection in local adaptation to dioxin-like aromatic hydrocarbon contaminants.

The pregnane X receptor (PXR) (nuclear receptor NR1I2) is a ligand activated transcription factor, mediating responses to diverse xenobiotic and endogenous chemicals. The properties of PXR in fish are not fully understood. Here we report on cloning and characterization of full-length PXR of zebrafish, Danio rerio, and pxr expression in vivo. Initial efforts gave a cDNA encoding a 430 amino acid protein identified as zebrafish pxr by phylogenetic and synteny analysis. The sequence of the cloned Pxr DNA binding domain (DBD) was highly conserved, with 74% identity to human PXR-DBD, while the ligand-binding domain (LBD) of the cloned sequence was only 44% identical to human PXR-LBD. Sequence variation among clones in the initial effort prompted sequencing of multiple clones from a single fish. There were two prominent variants, one sequence with S183, Y218 and H383 and the other with I183, C218 and N383, which we designate as alleles pxr*1 (nr1i2*1) and pxr*2 (nr1i2*2), respectively. In COS-7 cells co-transfected with a PXR-responsive reporter gene, the full-length Pxr*1 (the more common variant) was activated by known PXR agonists clotrimazole and pregnenolone 16α-carbonitrile but to a lesser extent than the full-length human PXR. Activation of full-length Pxr*1 was only 10% of that with the Pxr*1 LBD. Quantitative real time PCR analysis showed prominent expression of pxr in liver and eye, as well as brain and intestine of adult zebrafish. The pxr was expressed in heart and kidney at levels similar to that in intestine. The expression of pxr in liver was weakly induced by ligands for mammalian PXR or constitutive androstane receptor (NR1I3). The results establish a foundation for PXR studies in this vertebrate model. PXR allelic variation and the differences between the full-length PXR and the LBD in reporter assays have implications for assessing the action of PXR ligands in zebrafish., (Copyright © 2013 Elsevier B.V. All rights reserved.)

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

NF-E2-related factor 2 (NRF2; also called NFE2L2) and related NRF family members regulate antioxidant defenses by activating gene expression via antioxidant response elements (AREs), but their roles in embryonic development are not well understood. We report here that zebrafish (Danio rerio), an important developmental model species, possesses six nrf genes, including duplicated nrf1 and nrf2 genes. We cloned a novel zebrafish nrf2 paralog, nrf2b. The predicted Nrf2b protein sequence shares several domains with the original Nrf2 (now Nrf2a) but lacks the Neh4 transactivation domain. Zebrafish-human comparisons demonstrate conserved synteny involving nrf2 and hox genes, indicating that nrf2a and nrf2b are co-orthologs of human NRF2. nrf2a and nrf2b displayed distinct patterns of expression during embryonic development; nrf2b was more highly expressed at all stages. Embryos in which Nrf2a expression had been knocked down with morpholino oligonucleotides were more sensitive to tert-butylhydroperoxide but not tert-butylhydroquinone, whereas knockdown of Nrf2b did not affect sensitivity of embryos to either chemical. Gene expression profiling by microarray identified a specific role for Nrf2b as a negative regulator of several genes, including p53, cyclin G1, and heme oxygenase 1, in embryos. Nrf2a and Nrf2b exhibited different mechanisms of cross-talk with the Ahr2 signaling pathway. Together, these results demonstrate distinct roles for nrf2a and nrf2b, consistent with subfunction partitioning, and identify a novel negative regulatory role for Nrf2b during development. The identification of zebrafish nrf2 co-orthologs will facilitate new understanding of the multiple roles of NRF2 in protecting vertebrate embryos from oxidative damage.

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.

Atlantic killifish (Fundulus heteroclitus) inhabiting the PCB-contaminated Superfund site in New Bedford Harbor (MA, USA) have evolved genetic resistance to the toxic effects of these compounds. They also lack induction of cytochrome P4501A (CYP1A) and other aryl hydrocarbon receptor (AHR)-dependent responses after exposure to AHR agonists, suggesting an overall down-regulation of the AHR signaling pathway. In this study, we hypothesized that the genetic resistance is due to altered AHR expression resulting from hypermethylation of DNA in the promoter region of AHR genes in fish inhabiting New Bedford Harbor. To test this hypothesis, we cloned and sequenced AHR1 and AHR2 promoter regions and employed bisulfite conversion-polymerase chain reaction (BS-PCR) followed by clonal analysis to compare the methylation status of CpG islands of AHR1 and AHR2 in livers of adult killifish collected from New Bedford Harbor and a reference site (Scorton Creek, MA). No significant differences in methylation profiles were observed in either AHR1 or AHR2 promoter regions between NBH and SC fish. However, hypermethylation of the AHR1 promoter correlated with low expression of transcripts in the liver in both populations. In comparison to AHR1, hepatic mRNA expression of AHR2 is high and its promoter is hypomethylated. Taken together, our results suggest that genetic resistance to contaminants in NBH fish is not due to altered methylation of AHR promoter regions, but that promoter methylation may control tissue-specific expression of AHR genes in killifish., (Copyright © 2010 Elsevier B.V. All rights reserved.)

The aryl hydrocarbon receptor (AHR) mediates the toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and related compounds. In some mammalian cell lines, TCDD induces G1 cell cycle arrest, which depends on an interaction between the AHR and the retinoblastoma tumor suppressor (RB). Mammals possess one AHR, whereas fishes possess two or more AHR paralogs that differ in the domains important for AHR-RB interactions in mammals. To test the hypothesis that fish AHR paralogs differ in their ability to interact with RB, we cloned RB cDNA from Atlantic killifish, Fundulus heteroclitus, and studied the interactions of killifish RB protein with killifish AHR1 and AHR2. In coimmunoprecipitation experiments, in vitro-expressed killifish RB coprecipitated with both AHR1 and AHR2. Consistent with these results, both killifish AHR1 and AHR2 interacted with RB in mammalian two-hybrid assays. These results suggest that both fish AHR1 and AHR2 paralogs may have the potential to influence cell proliferation through interactions with RB.

The aryl hydrocarbon receptor (AHR) repressor (AHRR) inhibits AHR-mediated transcription and has been associated with reproductive dysfunction and tumorigenesis in humans. Previous studies have characterized the repressor function of AHRRs from mice and fish, but the human AHRR ortholog (AHRR(715)) appeared to be nonfunctional in vitro. Here, we report a novel human AHRR cDNA (AHRRDelta8) that lacks exon 8 of AHRR(715). AHRRDelta8 was the predominant AHRR form expressed in human tissues and cell lines. AHRRDelta8 effectively repressed AHR-dependent transactivation, whereas AHRR(715) was much less active. Similarly, AHRRDelta8, but not AHRR(715), formed a complex with AHR nuclear translocator (ARNT). Repression of AHR by AHRRDelta8 was not relieved by overexpression of ARNT or AHR coactivators, suggesting that competition for these cofactors is not the mechanism of repression. AHRRDelta8 interacted weakly with AHR but did not inhibit its nuclear translocation. In a survey of transcription factor specificity, AHRRDelta8 did not repress the nuclear receptor pregnane X receptor or estrogen receptor alpha but did repress hypoxia-inducible factor (HIF)-dependent signaling. AHRRDelta8-Pro(185) and -Ala(185) variants, which have been linked to human reproductive disorders, both were capable of repressing AHR or HIF. Together, these results identify AHRRDelta8 as the active form of human AHRR and reveal novel aspects of its function and specificity as a repressor.

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) and related halogenated aromatic hydrocarbons (HAHs) are highly toxic to most vertebrate animals, but there are dramatic differences in sensitivity among species and strains. Aquatic birds including the common tern (Sterna hirundo) are highly exposed to HAHs in the environment, but are up to 250-fold less sensitive to these compounds than the typical avian model, the domestic chicken (Gallus gallus). The mechanism of HAH toxicity involves altered gene expression subsequent to activation of the aryl hydrocarbon receptor (AHR), a basic helix-loop-helix-PAS transcription factor. AHR polymorphisms underlie mouse strain differences in sensitivity to HAHs and polynuclear aromatic hydrocarbons, but the role of the AHR in species differences in HAH sensitivity is not well understood. Here, we show that although chicken and tern AHRs both exhibit specific binding of [3H]TCDD, the tern AHR has a lower binding affinity and exhibits a reduced ability to support TCDD-dependent transactivation as compared to AHRs from chicken or mouse. We further show through use of chimeric AHR proteins and site-directed mutagenesis that the difference between the chicken and tern AHRs resides in the ligand-binding domain and that two amino acids (Val-325 and Ala-381) are responsible for the reduced activity of the tern AHR. Other avian species with reduced sensitivity to HAHs also possess these residues. These studies provide a molecular understanding of species differences in sensitivity to dioxin-like compounds and suggest an approach to using the AHR as a marker of dioxin susceptibility in wildlife.

The aryl hydrocarbon receptor (AhR) is an environmental carcinogen-activated transcription factor associated with tumorigenesis. High levels of apparently active AhR characterize a variety of tumors, even in the absence of environmental ligands. Despite this association between transformation and AhR upregulation, little is known of the transcriptional consequences of constitutive AhR activation. Here, the effects of constitutively active and environmental ligand-induced AhR on c-myc, an oncogene whose promoter contains six AhR-binding sites (AhREs (aryl hydrocarbon response elements)), were investigated. A reporter containing the human c-myc promoter, with its six AhREs and two NF-kappaB-binding sites, was constructed. This vector, and variants with deletions in the NF-kappaB and/or AhR-binding sites, was transfected into a human breast cancer cell line, Hs578T, which expresses high levels of apparently active, nuclear AhR. Results indicate that: (1) the AhR constitutively binds the c-myc promoter; (2) there is a low but significant baseline level of c-myc promoter activity, which is not regulated by NF-kappaB and is not affected by an environmental AhR ligand; (3) deletion of any one of the AhREs has no effect on constitutive reporter activity, while deletion of all six increases reporter activity approximately fivefold; (4) a similar increase in reporter activity occurs when constitutively active AhR is suppressed by transfection with an AhR repressor plasmid (AhRR); (5) AhRR transfection significantly increases background levels of endogenous c-myc mRNA and c-Myc protein. These results suggest that the AhR influences the expression of c-Myc, a protein critical to malignant transformation.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor that regulates gene expression following activation by TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) or a variety of other synthetic and natural compounds. Previous studies have identified two AHR genes, AHR1 and AHR2, in zebrafish (Danio rerio), a widely used model species for studying vertebrate development and an emerging model in developmental toxicology. Zebrafish AHR2 binds TCDD with high affinity, is transcriptionally active and has a major role in mediating the developmental toxicity of TCDD. Zebrafish AHR1 lacks the ability to bind TCDD and activate transcription, and has no known function. In the present study, we report a new zebrafish AHR, designated AHR1B, which shares 34% amino acid sequence identity with AHR1 (AHR1A). The ahr1b gene resides on chromosome 22, adjacent to ahr2, whereas the ahr1a gene is located on chromosome 16. AHR1B is expressed in embryos as early as 24 hours post-fertilization and increases through the next 2 days, but expression is not inducible by TCDD. In contrast with the previously identified AHR1A, in vitro-expressed AHR1B protein exhibits specific, high-affinity binding of [3H]TCDD. Furthermore, AHR1B is able to activate the transcription of a reporter gene under the control of AHR response elements with an efficacy comparable with that of AHR2, but with a higher EC50. We speculate that AHR1B may have a physiological role, such as in embryonic development, whereas AHR2 mediates the response to xenobiotics.

The effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and related compounds occur via the aryl hydrocarbon receptor (AHR), a member of the basic helix-loop-helix-Per-ARNT-Sim homology (bHLH-PAS) protein superfamily. A single AHR gene has been identified in mammals, whereas many fish species, including the Atlantic killifish (Fundulus heteroclitus) possess two distinct AHR genes (AHR1 and a novel form, AHR2). A mouse bHLH-PAS protein closely related to AHR and designated AHR repressor (AHRR) is induced by 3-methylcholanthrene and represses the transcriptional activity of the AHR. To determine whether AHRR is the mammalian ortholog of fish AHR2 and to investigate the mechanisms by which AHRR regulates AHR function, we cloned an AHRR ortholog in F. heteroclitus with high sequence identity to the mouse and human AHRRs. Killifish AHRR encodes a 680-residue protein with a predicted molecular mass of 75.2 kDa. We show that in vitro expressed AHRR proteins from human, mouse, and killifish all fail to bind [(3)H]TCDD or [(3)H]beta-naphthoflavone. In transient transfection experiments using a luciferase reporter gene under control of AHR response elements, killifish AHRR inhibited the TCDD-dependent transactivation function of both AHR1 and AHR2. AHRR mRNA is widely expressed in killifish tissues and is inducible by TCDD or polychlorinated biphenyls, but its expression is not altered in a population of fish exhibiting genetic resistance to these compounds. The F. heteroclitus AHRR promoter contains three putative AHR response elements. Both AHR1 and AHR2 activated transcription of luciferase driven by the AHRR promoter, and AHRR could repress its own promoter. Thus, AHRR is an evolutionarily conserved, TCDD-inducible repressor of AHR1 and AHR2 function. Phylogenetic analysis shows that AHRR, AHR1, and AHR2 are distinct genes, members of an AHR gene family; these three vertebrate AHR-like genes descended from a single invertebrate AHR.

The screening of liver and heart cDNA libraries from the teleost Fundulus heteroclitus with degenerate oligonucleotide probes to conserved alpha-helical regions in mammalian P450s resulted in the identification of two cDNAs that together represent a novel P450 subfamily, the CYP2Ns. Northern analysis demonstrated that CYP2N1 transcripts are most abundant in liver and intestine, whereas CYP2N2 mRNAs are most abundant in heart and brain. CYP2N1 and CYP2N2 proteins were co-expressed with NADPH-cytochrome P450 oxidoreductase in Sf9 insect cells, and their ability to metabolize arachidonic acid and xenobiotic substrates was examined. Both CYP2N1 and CYP2N2 metabolize arachidonic acid to epoxyeicosatrienoic acids. Epoxidation is highly regio- and enantioselective with preferential formation of (8R,9S)-epoxyeicosatrienoic acid (optical purities are 91 and 90% for CYP2N1 and CYP2N2, respectively) and (11R, 12S)-epoxyeicosatrienoic acid (optical purities are 92 and 70% for CYP2N1 and CYP2N2, respectively). CYP2N1 and CYP2N2 also catalyze the formation of a variety of hydroxyeicosatetraenoic acids. Both P450s have benzphetamine N-demethylase activities but show minimal alkoxyresorufin O-dealkylase activities. To investigate factors affecting CYP2N expression in vivo, CYP2N transcripts were examined following starvation and/or treatment with 12-O-tetradecanoyl phorbol-13-acetate. Intestinal CYP2N1 mRNAs decrease in starved and/or phorbol ester-treated fish, whereas intestinal CYP2N2 transcripts decrease only following phorbol ester treatment. Interestingly, cardiac CYP2N2 expression decreases following phorbol ester treatment but increases following starvation. These results demonstrate that members of this novel P450 subfamily encode early vertebrate forms of arachidonic acid catalysts that are widely expressed and are regulated by environmental factors. Given the wealth of information on the functional role of P450-derived arachidonate metabolites in mammals, we postulate that CYP2N1 and CYP2N2 products have similar biological functions in early vertebrates. The identity of the mammalian orthologue(s) of the CYP2Ns remains unknown.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which 2,3,7, 8-tetrachlorodibenzo-p-dioxin (TCDD) and related compounds cause altered gene expression and toxicity. The AHR belongs to an emerging multigene family of transcription factors possessing basic helix loop helix (bHLH) and Per-ARNT-Sim (PAS) domains. Most bHLH-PAS proteins occur as duplicates or "paralog groups" in mammals, but only a single mammalian AHR has been identified. Here we report the cDNA cloning of two distinct AHRs, designated FhAHR1 and FhAHR2, from a single vertebrate species, the teleost Fundulus heteroclitus (Atlantic killifish). Both Fundulus AHR proteins possess bHLH and PAS domains that are closely related to those of the mammalian AHR. FhAHR1 and FhAHR2 are highly divergent (40% overall amino acid identity; 61% identity in the N-terminal half), suggesting that they arose from a gene duplication predating the divergence of mammals and fish. Photoaffinity labeling with 2-azido-3-[(125)I]iodo-7, 8-dibromodibenzo-p-dioxin and velocity sedimentation analysis using 2,3,7,8-[1,6-(3)H]TCDD showed that both FhAHR1 and FhAHR2 exhibit specific, high-affinity binding of dioxins. Both AHRs also showed specific, TCDD- and ARNT-dependent interactions with a mammalian xenobiotic response element. The two Fundulus AHR genes displayed different tissue-specific patterns of expression; FhAHR1 transcripts were primarily expressed in brain, heart, ovary, and testis, while FhAHR2 transcripts were equally abundant in many tissues. Phylogenetic analysis demonstrated that Fundulus AHR1 is an ortholog of mammalian AHRs, while AHR2 forms in Fundulus and other fish are paralogous to Fundulus AHR1 and the mammalian AHRs and thus represent a novel vertebrate subfamily of ligand-binding AHRs.

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cause altered gene expression and toxicity. The AHR belongs to the basic helix-loop-helix/Per-ARNT-Sim (bHLH-PAS) family of transcriptional regulatory proteins, whose members play key roles in development, circadian rhythmicity, and environmental homeostasis; however, the normal cellular function of the AHR is not yet known. As part of a phylogenetic approach to understanding the function and evolutionary origin of the AHR, we sequenced the PAS homology domain of AHRs from several species of early vertebrates and performed phylogenetic analyses of these AHR amino acid sequences in relation to mammalian AHRs and 24 other members of the PAS family. AHR sequences were identified in a teleost (the killifish Fundulus heteroclitus), two elasmobranch species (the skate Raja erinacea and the dogfish Mustelus canis), and a jawless fish (the lamprey Petromyzon marinus). Two putative AHR genes, designated AHR1 and AHR2, were found both in Fundulus and Mustelus. Phylogenetic analyses indicate that the AHR2 genes in these two species are orthologous, suggesting that an AHR gene duplication occurred early in vertebrate evolution and that multiple AHR genes may be present in other vertebrates. Database searches and phylogenetic analyses identified four putative PAS proteins in the nematode Caenorhabditis elegans, including possible AHR and ARNT homologs. Phylogenetic analysis of the PAS gene family reveals distinct clades containing both invertebrate and vertebrate PAS family members; the latter include paralogous sequences that we propose have arisen by gene duplication early in vertebrate evolution. Overall, our analyses indicate that the AHR is a phylogenetically ancient protein present in all living vertebrate groups (with a possible invertebrate homolog), thus providing an evolutionary perspective to the study of dioxin toxicity and AHR function.

The PAS domain of a teleost Ah receptor was amplified using reverse transcription-PCR with degenerate primers containing inosine. The deduced amino acid sequence of the amplified cDNA fragment was 62-64% identical with the PAS domains of mammalian Ah receptors. These data demonstrate the homology of Ah receptors in mammals and fish, and reveal regions of this protein that are highly conserved between these diverse vertebrate groups.