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.

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.

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.

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.

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.

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.

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.

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.

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.

Functional characterization and evolutionary history of two aryl hydrocarbon receptor isoforms (AhR1 and AhR2) from avian species.

Dioxins including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) induce various toxic effects through the aryl hydrocarbon receptor (AhR) signaling pathway. Here, we investigated the structural and functional characteristics and molecular evolution of multiple AhRs in black-footed albatross (Phoebastria nigripes) and common cormorant (Phalacrocorax carbo). We report the complementary DNA sequences of two distinct AhRs, designated AhR1 and AhR2, from these species as well as the identification of an AhR2-like gene sequence from the chicken genome database. Phylogenetic analysis reveals that avian AhR1 and AhR2 are orthologous to mammalian AhR1 and fish AhR2, respectively, supporting the hypothesis that an ancestral AhR gene underwent a tandem duplication prior to the divergence of fish and tetrapod lineages. In vitro-expressed AhR1 and AhR2 isoforms from both albatross and cormorant exhibited specific binding to [3H]TCDD, as assessed by velocity sedimentation. An in vitro reporter gene transactivation assay revealed that both AhR1 and AhR2 are transcriptionally active, but AhR2 appears to have reduced transcriptional efficacy. Hepatic messenger RNA expression level of cormorant AhR1 was greater than that of AhR2. Together, these results suggest that AhR1 is the dominant form of avian AhRs, in contrast to fish, in which AhR2 is the major form. Comparative analysis of AhR diversity and gene synteny among chicken, zebrafish, and human suggests that additional, independent AhR duplications have occurred in the fish and tetrapod lineages following the initial tandem duplication on the ancestral chromosome. The identification and characterization of avian AhR1 and AhR2 provide new insight into the evolution of AhR structure and function in vertebrates.

AHR1B, a new functional aryl hydrocarbon receptor in zebrafish: tandem arrangement of ahr1b and ahr2 genes.

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.

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.

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.

Aryl hydrocarbon receptor polymorphisms and dioxin resistance in Atlantic killifish (Fundulus heteroclitus).

The aryl hydrocarbon receptor (AHR) gene encodes a ligand-activated transcription factor through which planar halogenated aromatic hydrocarbons (HAHs) such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) as well as polynuclear aromatic hydrocarbons (PAHs) cause altered gene expression and toxicity. To understand the role of AHR genetic variability in differential sensitivity to HAHs and PAHs, we are currently studying a population of the teleost Fundulus heteroclitus (Atlantic killifish) that has evolved genetic resistance to the toxic and biochemical effects of these compounds. Here, we report that the killifish AHR1 locus is highly polymorphic and that the frequencies of the major allele types differ between dioxin-sensitive and dioxin-resistant populations. Twenty-five single nucleotide polymorphisms (SNPs), nine of which are non-synonymous, were identified in the AHR1 coding sequence. Seven identified alleles were assigned to three groups, designated AHR1*1, AHR1*2 and AHR1*3. AHR1*1 alleles were under-represented in a population of dioxin- and polychlorinated biphenyl (PCB)-resistant fish from a PCB-contaminated Superfund site (New Bedford Harbor, Massachusetts, USA) compared to dioxin-sensitive fish from a less contaminated reference site (Scorton Creek, Massachusetts, USA). To determine the possible role of these AHR1 variants in differential HAH sensitivity, we expressed representative variant proteins from the two most divergent allelic groups (AHR1*1 and AHR1*3) by in-vitro transcription and translation and assessed their functional properties. AHR1*1A and AHR1*3A proteins displayed similar binding capacities and affinities for [H]TCDD. In transient transfection assays using mammalian cells, AHR1*1A and AHR1*3A exhibited similar abilities to support TCDD-dependent transactivation of a luciferase reporter gene under control of AHR-responsive enhancer elements. We discuss the possibility of other functional differences in AHR1 variants or their interaction with other killifish loci (AHR2, AHRR) that may contribute to differences in dioxin sensitivity.