Early life exposure to low levels of AHR agonist PCB126 (3,3’,4,4’,5- pentachlorobiphenyl) reprograms gene expression in adult brain

Author Posting. © The Author(s), 2017. This is the author's version of the work. It is posted here under a nonexclusive, irrevocable, paid-up, worldwide license granted to WHOI. It is made available for personal use, not for redistribution. The definitive version was published in Toxicological Sciences 160 (2017): 386-397, doi:10.1093/toxsci/kfx192.Early life exposure to environmental chemicals can have long-term consequences that are not always apparent until later in life. We recently demonstrated that developmental exposure of zebrafish to low, non-embryotoxic levels of 3,3’,4,4’,5-pentachlorobiphenyl (PCB126) did not affect larval behavior, but caused changes in adult behavior. The objective of this study was to investigate the underlying molecular basis for adult behavioral phenotypes resulting from early life exposure to PCB126. We exposed zebrafish embryos to PCB126 during early development and measured transcriptional profiles in whole embryos, larvae and adult male brains using RNA-sequencing. Early life exposure to 0.3 nM PCB126 induced cyp1a transcript levels in 2-dpf embryos, but not in 5-dpf larvae, suggesting transient activation of aryl hydrocarbon receptor with this treatment. No significant induction of cyp1a was observed in the brains of adults exposed as embryos to PCB126. However, a total of 2209 and 1628 genes were differentially expressed in 0.3 nM and 1.2 nM PCB126-exposed groups, respectively. KEGG pathway analyses of upregulated genes in the brain suggest enrichment of calcium signaling, MAPK and notch signaling, and lysine degradation pathways. Calcium is an important signaling molecule in the brain and altered calcium homeostasis could affect neurobehavior. The downregulated genes in the brain were enriched with oxidative phosphorylation and various metabolic pathways, suggesting that the metabolic capacity of the brain is impaired. Overall, our results suggest that PCB exposure during sensitive periods of early development alters normal development of the brain by reprogramming gene expression patterns, which may result in alterations in adult behavior

Role of DNA methylation of AHR1 and AHR2 promoters in differential sensitivity to PCBs in Atlantic Killifish, Fundulus heteroclitus

Author Posting. © The Author(s), 2010. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Aquatic Toxicology 101 (2011): 288-294, doi:10.1016/j.aquatox.2010.10.010.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.This work is funded in part by the Superfund Basic Research Program at Boston University to MEH (NIH Grant P42ES007381) and the postdoctoral scholar program at WHOI, with funding provided by the Dr. George D. Grice Postdoctoral Scholarship Fund to NA

Role of DNA methylation in altered gene expression patterns in adult zebrafish (Danio rerio) exposed to 3, 3’, 4, 4’, 5-pentachlorobiphenyl (PCB 126)

© The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Environmental Epigenetics 4 (2018): dvy005, doi:10.1093/eep/dvy005.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.This work was supported by the National Institute of Health Outstanding New Environmental Scientist Award to NA (NIH R01ES024915) and Woods Hole Center for Oceans and Human Health [National Institutes of Health (NIH) grant P01ES021923 and National Science Foundation Grant OCE-1314642 to M. Hahn, J. Stegeman, NA and SK]

Interaction of Fish Aryl Hydrocarbon Receptor Paralogs (AHR1 and AHR2) with the Retinoblastoma Protein

The aryl hydrocarbon receptor (AHR) mediates the toxic effects of 2,3,7,8-tetrachlorodibenzo-pdioxin (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

Nrf2 and Nrf2-related proteins in development and developmental toxicity : insights from studies in zebrafish (Danio rerio)

© The Author(s), 2015. This is the author's version of the work and is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Free Radical Biology and Medicine 88B (2015): 275-289, doi:10.1016/j.freeradbiomed.2015.06.022.Oxidative stress is an important mechanism of chemical toxicity, contributing to developmental toxicity and teratogenesis as well as to cardiovascular and neurodegenerative diseases and diabetic embryopathy. Developing animals are especially sensitive to effects of chemicals that disrupt the balance of processes generating reactive species and oxidative stress, and those anti-oxidant defenses that protect against oxidative stress. The expression and inducibility of anti-oxidant defenses through activation of NFE2-related factor 2 (Nrf2) and related proteins is an essential process affecting the susceptibility to oxidants, but the complex interactions of Nrf2 in determining embryonic response to oxidants and oxidative stress are only beginning to be understood. The zebrafish (Danio rerio) is an established model in developmental biology and now also in developmental toxicology and redox signaling. Here we review the regulation of genes involved in protection against oxidative stress in developing vertebrates, with a focus on Nrf2 and related cap’n’collar (CNC)-basic-leucine zipper (bZIP) transcription factors. Vertebrate animals including zebrafish share Nfe2, Nrf1, Nrf2, and Nrf3 as well as a core set of genes that respond to oxidative stress, contributing to the value of zebrafish as a model system with which to investigate the mechanisms involved in regulation of redox signaling and the response to oxidative stress during embryolarval development. Moreover, studies in zebrafish have revealed nrf and keap1 gene duplications that provide an opportunity to dissect multiple functions of vertebrate NRF genes, including multiple sensing mechanisms involved in chemical-specific effects.This work was supported in part by National Institutes of Health grants R01ES016366 (MEH), R01ES015912 (JJS), and F32ES017585 (ART-L).2016-06-2

Sequence and functional characterization of hypoxia inducible factors, HIF1α, HIF2αa, and HIF3α, from the estuarine fish, Fundulus heteroclitus

Author Posting. © The Author(s), 2016. This is the author's version of the work. It is posted here under a nonexclusive, irrevocable, paid-up, worldwide license granted to WHOI. It is made available for personal use, not for redistribution. The definitive version was published in American Journal of Physiology-Regulatory Integrative and Comparative Physiology, 312 (2017): R412-R425, doi:10.1152/ajpregu.00402.2016.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 occurs 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α, 2α, and 3α independently formed complexes with the β subunit (ARNT) 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.This research was supported in part by the National Science Foundation (IBN-0236494 and DEB-1120263) and by National Institute of Environmental Health Sciences (NIEHS) grant P42ES007381 (Superfund Basic Research Program at Boston University). Data interpretation was aided by reference to a preliminary draft of the F. heteroclitus genome sequence, which was supported by funding from the National Science Foundation (collaborative research grants DEB-1120512, DEB-1265282, DEB-1120013, DEB-1120263, DEB-1120333, DEB-1120398)

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

© The Author(s), 2011. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in BMC Genomics 12 (2011): 263, doi:10.1186/1471-2164-12-263.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). 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). 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.This work was supported in part by National Institutes of Health grants P42ES007381 (Superfund Basic Research Program at Boston University) and by a grant from the WHOI Ocean Life Institute

Development and characterization of polyclonal antibodies against the aryl hydrocarbon receptor protein family (AHR1, AHR2, and AHR repressor) of Atlantic killifish Fundulus heteroclitus

Author Posting. © The Authors, 2005. This is the author's version of the work. It is posted here by permission of Elsevier B. V. for personal use, not for redistribution. The definitive version was published in Comparative Biochemistry and Physiology Part C: Toxicology & Pharmacology 142 (2006): 85-94, doi:10.1016/j.cbpc.2005.10.013.The aryl hydrocarbon receptor (AHR) and AHR repressor (AHRR) proteins regulate gene expression in response to some halogenated aromatic hydrocarbons and polycyclic aromatic hydrocarbons. The Atlantic killifish is a valuable model of the AHR signaling pathway, but antibodies are not available to fully characterize AHR and AHRR proteins. Using bacterially expressed AHRs, we developed specific and sensitive polyclonal antisera against the killifish AHR1, AHR2, and AHRR. In immunoblots, these antibodies recognized full-length killifish AHR and AHRR proteins synthesized in rabbit reticulocyte lysate, proteins expressed in mammalian cells transfected with killifish AHR and AHRR constructs, and AHR proteins in cytosol preparations from killifish tissues. Killifish AHR1 and AHR2 proteins were detected in brain, gill, kidney, heart, liver, and spleen. Antisera specifically precipitated their respective target proteins in immunoprecipitation experiments with in vitro-expressed proteins. Killifish ARNT2 co-precipitated with AHR1 and AHR2. These sensitive, specific, and versatile antibodies will be valuable to researchers investigating AHR signaling and other physiological processes involving AHR and AHRR proteins.Funding for this research was provided by National Institutes of Health, National Research Service Award (F32 ES05935) from the National Institute of Environmental Health Sciences (RRM), NIEHS Superfund Basic Research Program Grant P42 ES007381 at Boston University (MEH), and the Oliver S. and Jennie R. Donaldson Charitable Trust (MEH and RRM)

Ryanodine receptor and FK506 binding protein 1 in the Atlantic killifish (Fundulus heteroclitus) : a phylogenetic and population-based comparison

Author Posting. © The Author(s), 2017. This is the author's version of the work. It is posted here under a nonexclusive, irrevocable, paid-up, worldwide license granted to WHOI. It is made available for personal use, not for redistribution. The definitive version was published in Aquatic Toxicology 192 (2017): 105-115, doi:10.1016/j.aquatox.2017.09.002.Non-dioxin-like polychlorinated biphenyls (NDL PCBs) activate ryanodine receptors (RyR), microsomal Ca2+ 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.6 kDa (genes FKBP1A; FKBP1B), which promote NDL PCB-triggered Ca2+ 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.This research was supported by the KC Donnelly Research Externship made possible by the National Institute of Environmental Health Sciences’ Superfund Research Program (EBH) and the Superfund Research Programs at UC Davis (INP and EBH; P42ES004699) and Boston University (JJS, JVG, MEH, SIK; P42ES007381). Additional support was provided by the National Institute of Health (INP; R01 ES014901; and P01 AR052354) and by National Science Foundation collaborative research grants (MEH and SIK; DEB-1265282 and DEB-1120263). This research was also supported in part by an appointment (to BC) with the Postdoctoral Research Program at the U.S. Environmental Protection (US EPA) Office of Research and Development administered by the Oak Ridge Institute for Science and Education (ORISE) through Interagency Agreement No. DW92429801 between the U.S. Department of Energy and the US EPA