Formulation Controls the Potential Neuromuscular Toxicity of Polyethylene Photoproducts in Developing Zebrafish.

Sunlight transforms plastic into water-soluble products, the potential toxicity of which remains unresolved, particularly for vertebrate animals. We evaluated acute toxicity and gene expression in developing zebrafish larvae after 5 days of exposure to photoproduced (P) and dark (D) leachates from additive-free polyethylene (PE) film and consumer-grade, additive-containing, conventional, and recycled PE bags. Using a "worst-case" scenario, with plastic concentrations exceeding those found in natural waters, we observed no acute toxicity. However, at the molecular level, RNA sequencing revealed differences in the number of differentially expressed genes (DEGs) for each leachate treatment: thousands of genes (5442 P, 577 D) for the additive-free film, tens of genes for the additive-containing conventional bag (14 P, 7 D), and none for the additive-containing recycled bag. Gene ontology enrichment analyses suggested that the additive-free PE leachates disrupted neuromuscular processes via biophysical signaling; this was most pronounced for the photoproduced leachates. We suggest that the fewer DEGs elicited by the leachates from conventional PE bags (and none from recycled bags) could be due to differences in photoproduced leachate composition caused by titanium dioxide-catalyzed reactions not present in the additive-free PE. This work demonstrates that the potential toxicity of plastic photoproducts can be product formulation-specific.

Identification of environmental factors that promote intestinal inflammation.

Genome-wide association studies have identified risk loci linked to inflammatory bowel disease (IBD)1-a complex chronic inflammatory disorder of the gastrointestinal tract. The increasing prevalence of IBD in industrialized countries and the augmented disease risk observed in migrants who move into areas of higher disease prevalence suggest that environmental factors are also important determinants of IBD susceptibility and severity2. However, the identification of environmental factors relevant to IBD and the mechanisms by which they influence disease has been hampered by the lack of platforms for their systematic investigation. Here we describe an integrated systems approach, combining publicly available databases, zebrafish chemical screens, machine learning and mouse preclinical models to identify environmental factors that control intestinal inflammation. This approach established that the herbicide propyzamide increases inflammation in the small and large intestine. Moreover, we show that an AHR-NF-κB-C/EBPß signalling axis operates in T cells and dendritic cells to promote intestinal inflammation, and is targeted by propyzamide. In conclusion, we developed a pipeline for the identification of environmental factors and mechanisms of pathogenesis in IBD and, potentially, other inflammatory diseases.

An aryl hydrocarbon receptor from the caecilian Gymnopis multiplicata suggests low dioxin affinity in the ancestor of all three amphibian orders.

The aryl hydrocarbon receptor (AHR) plays pleiotropic roles in the development and physiology of vertebrates in conjunction with xenobiotic and endogenous ligands. It is best known for mediating the toxic effects of dioxin-like pollutants such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). While most vertebrates possess at least one AHR that binds TCDD tightly, amphibian AHRs bind TCDD with very low affinity. Previous analyses of AHRs from Xenopus laevis (a frog; order Anura) and Ambystoma mexicanum (a salamander; order Caudata) identified three amino acid residues in the ligand-binding domain (LBD) that underlie low-affinity binding. In X. laevis AHR1ß, these are A354, A370, and N325. Here we extend the analysis of amphibian AHRs to the caecilian Gymnopis multiplicata, representing the remaining extant amphibian order, Gymnophiona. G. multiplicata AHR groups with the monophyletic vertebrate AHR/AHR1 clade. The LBD includes all three signature residues of low TCDD affinity, and a structural homology model suggests that its architecture closely resembles those of other amphibians. In transactivation assays, the EC50 for reporter gene induction by TCDD was 17.17 nM, comparable to X. laevis AhR1ß (26.23 nM) and Ambystoma AHR (34.09 nM) and dramatically higher than mouse AhR (0.13 nM), a trend generally reflected in direct measures of TCDD binding. These shared properties distinguish amphibian AHRs from the high-affinity proteins typical of both vertebrate groups that diverged earlier (teleost fish) and those that appeared more recently (other tetrapods). These findings suggest the hypothesis that AHRs with low TCDD affinity represent a characteristic that evolved in a common ancestor of all three extant amphibian groups.

Molecular and Functional Properties of the Atlantic Cod (Gadus morhua) Aryl Hydrocarbon Receptors Ahr1a and Ahr2a.

The aryl hydrocarbon receptor (Ahr) is a ligand-activated transcription factor that mediates the toxicity of halogenated and polycyclic aromatic hydrocarbons in vertebrates. Atlantic cod (Gadus morhua) has recently emerged as a model organism in environmental toxicology studies, and increased knowledge of Ahr-mediated responses to xenobiotics is imperative. Genome mining and phylogenetic analyses revealed two Ahr-encoding genes in the Atlantic cod genome, gmahr1a and gmahr2a. In vitro binding assays showed that both gmAhr proteins bind to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), but stronger binding to gmAhr1a was observed. Transactivation studies with a reporter gene assay revealed that gmAhr1a is one order of magnitude more sensitive to TCDD than gmAhr2a, but the maximal responses of the receptors were similar. Other well-known Ahr agonists, such as ß-naphthoflavone (BNF), 3,3',4,4',5-pentachlorobiphenyl (PCB126), and 6-formylindolo[3,2-b]carbazole (FICZ), also activated the gmAhr proteins, but gmAhr1a was, in general, the more sensitive receptor and produced the highest efficacies. The induction of cyp1a in exposed precision-cut cod liver slices confirmed the activation of the Ahr signaling pathway ex vivo. In conclusion, the differences in transcriptional activation by gmAhr's with various agonists, the distinct binding properties with TCDD and BNF, and the distinct tissue-specific expression profiles indicate different functional specializations of the Atlantic cod Ahr's.

Aryl hydrocarbon receptor-mediated activity of gas-phase ambient air derived from passive sampling and an in vitro bioassay.

The gaseous fraction of hydrophobic organic contaminants (HOCs) in ambient air appears to be responsible for a significant portion of aryl hydrocarbon receptor (AhR)-mediated activity, but the majority of compounds contributing to this activity remain unidentified. The present study investigated the use of polyethylene passive samplers to isolate gaseous HOCs from ambient air for use in in vitro bioassays and to improve our understanding of the toxicological relevance of the gaseous fraction of ambient air in urban and residential environments. Concentrations of polycyclic aromatic hydrocarbons (PAHs) and organic flame retardants were measured in polyethylene passive sampler extracts. Extracts were also analyzed using an in vitro bioassay to measure AhR-mediated activity. Bioassay-derived benzo[a]pyrene (BaP) equivalents (BaP-Eqbio ), a measure of potency of HOC mixtures, were greatest in the downtown Cleveland area and lowest at rural/residential sites further from the city center. The BaP-Eqbio was weakly correlated with concentrations of 2-ring alkyl/substituted PAHs and one organophosphate flame retardant, ethylhexyl diphenyl phosphate. Potency predicted based on literature-derived induction equivalency factors (IEFs) explained only 2 to 23% of the AhR-mediated potency observed in bioassay experiments. Our results suggests that health risks of gaseous ambient air pollution predicted using data from targeted chemical analysis may underestimate risks of exposure, most likely due to augmentation of potency by unmonitored chemicals in the mixture, and the lack of relevant IEFs for many targeted analytes. Environ Toxicol Chem 2019;38:748-759. © 2019 SETAC.

Regulation of Ahr signaling by Nrf2 during development: Effects of Nrf2a deficiency on PCB126 embryotoxicity in zebrafish (Danio rerio).

The embryotoxicity of co-planar PCBs is regulated by the aryl hydrocarbon receptor (Ahr), and has been reported to involve oxidative stress. Ahr participates in crosstalk with another transcription factor, Nfe2l2, or Nrf2. Nrf2 binds to antioxidant response elements to regulate the adaptive response to oxidative stress. To explore aspects of the crosstalk between Nrf2 and Ahr and its impact on development, we used zebrafish (Danio rerio) with a mutated DNA binding domain in Nrf2a (nrf2a(fh318/fh318)), rendering these embryos more sensitive to oxidative stress. Embryos were exposed to 2 nM or 5 nM PCB126 at 24 h post fertilization (prim-5 stage of pharyngula) and examined for gene expression and morphology at 4 days post fertilization (dpf; protruding - mouth stage). Nrf2a mutant eleutheroembryos were more sensitive to PCB126 toxicity at 4 dpf, and in the absence of treatment also displayed some subtle developmental differences from wildtype embryos, including delayed inflation of the swim bladder and smaller yolk sacs. We used qPCR to measure changes in expression of the nrf gene family, keap1a, keap1b, the ahr gene family, and known target genes. cyp1a induction by PCB126 was enhanced in the Nrf2a mutants (156-fold in wildtypes vs. 228-fold in mutants exposed to 5 nM). Decreased expression of heme oxygenase (decycling) 1 (hmox1) in the Nrf2a mutants was accompanied by increased nrf2b expression. Target genes of Nrf2a and AhR2, NAD(P)H:quinone oxidoreductase 1 (nqo1) and glutathione S-transferase, alpha-like (gsta1), showed a 2-5-fold increase in expression in the Nrf2a mutants as compared to wildtype. This study elucidates the interaction between two important transcription factor pathways in the developmental toxicity of co-planar PCBs.

Naturally occurring marine brominated indoles are aryl hydrocarbon receptor ligands/agonists.

The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that mediates the toxic and biological effects of structurally diverse chemicals, including the environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). As part of a larger effort to identify the full spectrum of chemicals that can bind to and activate the AhR, we have examined the ability of several naturally occurring marine-derived brominated indoles and brominated (methylthio)indoles (collectively referred to as brominated indoles) to bind to the AhR and stimulate AhR-dependent gene expression. Incubation of mouse, rat, and guinea pig recombinant cell lines containing a stably transfected AhR-responsive luciferase reporter gene with eight brominated indoles revealed that all compounds stimulated luciferase reporter gene activity, although some species-specific differences were observed. All compounds induced significantly more luciferase activity when incubated with cells for 4 h as compared to 24 h, demonstrating that these compounds are transient activators of the AhR signaling pathway. Three of the brominated indoles induced CYP1A1 mRNA in human HepG2 cells in vitro and Cyp1a mRNA in zebrafish embryos in vivo. The identification of the brominated indoles as direct ligands and activators/agonists of the AhR was confirmed by their ability to compete with [(3)H]TCDD for binding to the AhR and to stimulate AhR transformation and DNA binding in vitro. Taken together, these results indicate that marine-derived brominated indoles are members of a new class of naturally occurring AhR agonists.

An Aryl Hydrocarbon Receptor from the Salamander Ambystoma mexicanum Exhibits Low Sensitivity to 2,3,7,8-Tetrachlorodibenzo-p-dioxin.

Structural features of the aryl hydrocarbon receptor (AHR) can underlie species- and population-specific differences in its affinity for 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). These differences often explain variations in TCDD toxicity. Frogs are relatively insensitive to dioxin, and Xenopus AHRs bind TCDD with low affinity. Weak TCDD binding results from the combination of three residues in the ligand-binding domain: A354 and A370, and N325. Here we sought to determine whether this mechanism of weak TCDD binding is shared by other amphibian AHRs. We isolated an AHR cDNA from the Mexican axolotl (Ambystoma mexicanum). The encoded polypeptide contains identical residues at positions that confer low TCDD affinity to X. laevis AHRs (A364, A380, and N335), and homology modeling predicts they protrude into the binding cavity. Axolotl AHR bound one-tenth the TCDD of mouse AHR in velocity sedimentation analysis, and in transactivation assays, the EC50 for TCDD was 23 nM, similar to X. laevis AHR1ß (27 nM) and greater than AHR containing the mouse ligand-binding domain (0.08 nM). Sequence, modeled structure, and function indicate that axolotl AHR binds TCDD weakly, predicting that A. mexicanum lacks sensitivity toTCDD toxicity. We hypothesize that this characteristic of axolotl and Xenopus AHRs arose in a common ancestor of the Caudata and Anura.

Developmental exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin alters DNA methyltransferase (dnmt) expression in zebrafish (Danio rerio).

DNA methylation is one of the most important epigenetic modifications involved in the regulation of gene expression. The DNA methylation reaction is catalyzed by DNA methyltransferases (DNMTs). Recent studies have demonstrated that toxicants can affect normal development by altering DNA methylation patterns, but the mechanisms of action are poorly understood. Hence, we tested the hypothesis that developmental exposure to TCDD affects dnmt gene expression patterns. Zebrafish embryos were exposed to 5nM TCDD for 1h from 4 to 5h post-fertilization (hpf) and sampled at 12, 24, 48, 72, and 96 hpf to determine dnmt gene expression and DNA methylation patterns. We performed a detailed analysis of zebrafish dnmt gene expression during development and in adult tissues. Our results demonstrate that dnmt3b genes are highly expressed in early stages of development, and dnmt3a genes are more abundant in later stages. TCDD exposure upregulated dnmt1 and dnmt3b2 expression, whereas dnmt3a1, 3b1, and 3b4 are downregulated following exposure. We did not observe any TCDD-induced differences in global methylation or hydroxymethylation levels, but the promoter methylation of aryl hydrocarbon receptor (AHR) target genes was altered. In TCDD-exposed embryos, AHR repressor a (ahrra) and c-fos promoters were differentially methylated. To characterize the TCDD effects on DNMTs, we cloned the dnmt promoters with xenobiotic response elements and conducted AHR transactivation assays using a luciferase reporter system. Our results suggest that ahr2 can regulate dnmt3a1, dnmt3a2, and dnmt3b2 expression. Overall, we demonstrate that developmental exposure to TCDD alters dnmt expression and DNA methylation patterns.

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.

The transcriptional response to oxidative stress during vertebrate development: effects of tert-butylhydroquinone and 2,3,7,8-tetrachlorodibenzo-p-dioxin.

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.

In silico identification of an aryl hydrocarbon receptor antagonist with biological activity in vitro and in vivo.

The aryl hydrocarbon receptor (AHR) is critically involved in several physiologic processes, including cancer progression and multiple immune system activities. We, and others, have hypothesized that AHR modulators represent an important new class of targeted therapeutics. Here, ligand shape-based virtual modeling techniques were used to identify novel AHR ligands on the basis of previously identified chemotypes. Four structurally unique compounds were identified. One lead compound, 2-((2-(5-bromofuran-2-yl)-4-oxo-4H-chromen-3-yl)oxy)acetamide (CB7993113), was further tested for its ability to block three AHR-dependent biologic activities: triple-negative breast cancer cell invasion or migration in vitro and AHR ligand-induced bone marrow toxicity in vivo. CB7993113 directly bound both murine and human AHR and inhibited polycyclic aromatic hydrocarbon (PAH)- and TCDD-induced reporter activity by 75% and 90% respectively. A novel homology model, comprehensive agonist and inhibitor titration experiments, and AHR localization studies were consistent with competitive antagonism and blockade of nuclear translocation as the primary mechanism of action. CB7993113 (IC50 3.3 × 10(-7) M) effectively reduced invasion of human breast cancer cells in three-dimensional cultures and blocked tumor cell migration in two-dimensional cultures without significantly affecting cell viability or proliferation. Finally, CB7993113 effectively inhibited the bone marrow ablative effects of 7,12-dimethylbenz[a]anthracene in vivo, demonstrating drug absorption and tissue distribution leading to pharmacological efficacy. These experiments suggest that AHR antagonists such as CB7993113 may represent a new class of targeted therapeutics for immunomodulation and/or cancer therapy.

Identification of cinnabarinic acid as a novel endogenous aryl hydrocarbon receptor ligand that drives IL-22 production.

The aryl hydrocarbon receptor (AHR) binds to environmental toxicants including synthetic halogenated aromatic hydrocarbons and is involved in a diverse array of biological processes. Recently, the AHR was shown to control host immunity by affecting the balance between inflammatory T cells that produce IL-17 (Th17) and IL-22 versus regulatory T cells (Treg) involved in tolerance. While environmental AHR ligands can mediate this effect, endogenous ligands are likely to be more relevant in host immune responses. We investigated downstream metabolites of tryptophan as potential AHR ligands because (1) tryptophan metabolites have been implicated in regulating the balance between Th17 and Treg cells and (2) many of the AHR ligands identified thus far are derivatives of tryptophan. We characterized the ability of tryptophan metabolites to bind and activate the AHR and to increase IL-22 production in human T cells. We report that the tryptophan metabolite, cinnabarinic acid (CA), is an AHR ligand that stimulates the differentiation of human and mouse T cells producing IL-22. We compare the IL-22-stimulating activity of CA to that of other tryptophan metabolites and define stimulation conditions that lead to CA production from immune cells. Our findings link tryptophan metabolism to AHR activation and define a novel endogenous AHR agonist with potentially broad biological functions.

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

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.

Aryl hydrocarbon receptor (AHR) in the cnidarian Nematostella vectensis: comparative expression, protein interactions, and ligand binding.

The aryl hydrocarbon receptor (AHR) is a member of the basic helix-loop-helix/Per-ARNT-Sim (bHLH-PAS) family of transcription factors and has diverse roles in development, physiology, and environmental sensing in bilaterian animals. Studying the expression of conserved genes and function of proteins in outgroups to protostomes and deuterostomes assists in understanding the antiquity of gene function and deciphering lineage-specific differences in these bilaterian clades. We describe the developmental expression of AHR from the sea anemone Nematostella vectensis and compare its expression with three other members of the bHLH-PAS family (AHR nuclear translocator (ARNT), Cycle, and a proto-Single-Minded/Trachealess). NvAHR expression was highest early in the larval stage with spatial expression in the basal portion of the ectoderm that became increasingly restricted to the oral pole with concentrated expression in tentacles of the juvenile polyp. The other bHLH-PAS genes showed a divergent expression pattern in later larval stages and polyps, in which gene expression was concentrated in the aboral end, with broader expression in the endoderm later in development. In co-immunoprecipitation assays, we found no evidence for heterodimerization of AHR with ARNT, contrary to the conservation of this specific interaction in all bilaterians studied to date. Similar to results with other invertebrate AHRs but in contrast to vertebrate AHRs, NvAHR failed to bind two prototypical xenobiotic AHR ligands (2,3,7,8-tetrachlorodibenzo-p-dioxin, ß-naphthoflavone). Together, our data suggest that AHR's original function in Eumetazoa likely involved developmental patterning, potentially of neural tissue. The role of heterodimerization in the function of AHR may have arisen after the cnidarian-bilaterian ancestor. The absence of xenobiotic binding to NvAHR further supports a hypothesis for a derived role of this protein in chemical sensing within the chordates.

Specific ligand binding domain residues confer low dioxin responsiveness to AHR1ß of Xenopus laevis.

The aryl hydrocarbon receptor (AHR) is a Per-ARNT-Sim (PAS) family protein that mediates the toxicity of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in vertebrates. Frogs are remarkably insensitive to TCDD, and AHRs from Xenopus laevis bind TCDD with low affinity. We sought to identify structural features of X. laevis AHR1ß associated with low TCDD sensitivity. Substitution of the entire ligand binding domain (LBD) with the corresponding sequence from mouse AHR(b-1) dramatically increased TCDD responsiveness in transactivation assays. To identify the amino acid residues responsible, we constructed a comparative model of the AHR1ß LBD using homologous domains of PAS proteins HIF2α and ARNT. The model revealed an internal cavity with dimensions similar to those of the putative binding cavity of mouse AHR(b-1), suggesting the importance of side chain interactions over cavity size. Of residues with side chains clearly pointing into the cavity, only two differed from the mouse sequence. When A354, located within a conserved ß-strand, was changed to serine, the corresponding mouse residue, the EC50 for TCDD decreased more than 15-fold. When N325 was changed to serine, the EC50 decreased 3-fold. When the mutations were combined, the EC50 decreased from 18.6 to 0.8 nM, the value nearly matching the TCDD sensitivity of mouse AHR. Velocity sedimentation analysis confirmed that mutant frog AHRs exhibited correspondingly increased levels of TCDD binding. We also assayed mutant AHRs for responsiveness to a candidate endogenous ligand, 6-formylindolo[3,2-b]carbazole (FICZ). Mutations that increased sensitivity to TCDD also increased sensitivity to FICZ. This comparative study represents a novel approach to discerning fundamental information about the structure of AHR and its interactions with biologically important agonists.

Comparative analysis of homology models of the AH receptor ligand binding domain: verification of structure-function predictions by site-directed mutagenesis of a nonfunctional receptor.

The aryl hydrocarbon receptor (AHR) is a ligand-dependent transcription factor that mediates the biological and toxic effects of a wide variety of structurally diverse chemicals, including the toxic environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). While significant interspecies differences in AHR ligand binding specificity, selectivity, and response have been observed, the structural determinants responsible for those differences have not been determined, and homology models of the AHR ligand-binding domain (LBD) are available for only a few species. Here we describe the development and comparative analysis of homology models of the LBD of 16 AHRs from 12 mammalian and nonmammalian species and identify the specific residues contained within their ligand binding cavities. The ligand-binding cavity of the fish AHR exhibits differences from those of mammalian and avian AHRs, suggesting a slightly different TCDD binding mode. Comparison of the internal cavity in the LBD model of zebrafish (zf) AHR2, which binds TCDD with high affinity, to that of zfAHR1a, which does not bind TCDD, revealed that the latter has a dramatically shortened binding cavity due to the side chains of three residues (Tyr296, Thr386, and His388) that reduce the amount of internal space available to TCDD. Mutagenesis of two of these residues in zfAHR1a to those present in zfAHR2 (Y296H and T386A) restored the ability of zfAHR1a to bind TCDD and to exhibit TCDD-dependent binding to DNA. These results demonstrate the importance of these two amino acids and highlight the predictive potential of comparative analysis of homology models from diverse species. The availability of these AHR LBD homology models will facilitate in-depth comparative studies of AHR ligand binding and ligand-dependent AHR activation and provide a novel avenue for examining species-specific differences in AHR responsiveness.

Nrf2b, novel zebrafish paralog of oxidant-responsive transcription factor NF-E2-related factor 2 (NRF2).

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.

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.

Mechanistic basis of resistance to PCBs in Atlantic tomcod from the Hudson River.

The mechanistic basis of resistance of vertebrate populations to contaminants, including Atlantic tomcod from the Hudson River (HR) to polychlorinated biphenyls (PCBs), is unknown. HR tomcod exhibited variants in the aryl hydrocarbon receptor 2 (AHR2) that were nearly absent elsewhere. In ligand-binding assays, AHR2-1 protein (common in the HR) was impaired as compared to widespread AHR2-2 in binding TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) and in driving expression in reporter gene assays in AHR-deficient cells treated with TCDD or PCB126. We identified a six-base deletion in AHR2 as the basis of resistance and suggest that the HR population has undergone rapid evolution, probably due to contaminant exposure. This mechanistic basis of resistance in a vertebrate population provides evidence of evolutionary change due to selective pressure at a single locus.