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.

The African coelacanth genome provides insights into tetrapod evolution.

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

Molecular adaptation to high pressure in cytochrome P450 1A and aryl hydrocarbon receptor systems of the deep-sea fish Coryphaenoides armatus.

Limited knowledge of the molecular evolution of deep-sea fish proteomes so far suggests that a few widespread residue substitutions in cytosolic proteins binding hydrophilic ligands contribute to resistance to the effects of high hydrostatic pressure (HP). Structure-function studies with additional protein systems, including membrane bound proteins, are essential to provide a more general picture of adaptation in these extremophiles. We explored molecular features of HP adaptation in proteins binding hydrophobic ligands, either in lipid bilayers (cytochrome P450 1A - CYP1A) or in the cytosol (the aryl hydrocarbon receptor - AHR), and their partners P450 oxidoreductase (POR) and AHR nuclear translocator (ARNT), respectively. Cloning studies identified the full-length coding sequence of AHR, CYP1A and POR, and a partial sequence of ARNT from Coryphaenoides armatus, an abyssal gadiform fish thriving down to 5000m depth. Inferred protein sequences were aligned with many non-deep-sea homologs to identify unique amino acid substitutions of possible relevance in HP adaptation. Positionally unique substitutions of various physicochemical properties were found in all four proteins, usually at sites of strong-to-absolute residue conservation. Some were in domains deemed important for protein-protein interaction or ligand binding. In addition, some involved removal or addition of beta-branched residues; local modifications of beta-branched residue patterns could be important to HP adaptation. In silico predictions further suggested that some unique substitutions might substantially modulate the flexibility of the polypeptide segment in which they are found. Repetitive motifs unique to the abyssal fish AHR were predicted to be rich in glycosylation sites, suggesting that post-translational changes could be involved in adaptation as well. Recombinant CYP1A and AHR showed functional properties (spectral characteristics, catalytic activity and ligand binding) that demonstrate proper folding at 1atm, indicating that they could be used as deep-sea fish protein models to further evaluate protein function under pressure. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone".

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

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

Integrating Monitoring and Genetic Methods To Infer Historical Risks of PCBs and DDE to Common and Roseate Terns Nesting Near the New Bedford Harbor Superfund Site (Massachusetts, USA).

Common and roseate terns are migratory piscivorous seabirds with major breeding colonies within feeding range of the polychlorinated biphenyl (PCB)-contaminated New Bedford Harbor (NBH, MA, USA) Superfund site. Our longitudinal study shows that before PCB discharges into NBH ceased (late 1970s), tern eggs had very high but variable PCB concentrations. However, egg concentrations of PCBs as well as DDE (1,1-bis(p-chlorophenyl)-2,2-dichloroethene), the degradation product of the ubiquitous global contaminant DDT (1,1,1-trichloro-2,2-bis(p-chlorophenyl) ethane), have since declined. Rate constants for temporal decline of PCB congeners in tern eggs varied inversely with log10KOW (n-octanol-water partition coefficient), shifting egg congener patterns away from those characterizing NBH sediment. To estimate the toxic effects on tern eggs of PCB dioxin-like congener (DLC) exposures, we extrapolated published laboratory data on common terns to roseate terns by characterizing genetic and functional similarities in species aryl hydrocarbon receptors (AHRs), which mediate DLC sensitivity. Our assessment of contaminant risks suggests that terns breeding near NBH were exposed historically to toxic levels of PCBs and DDE; however, acute effects on tern egg development have become less likely since the 1970s. Our approach demonstrates how comparative genetics at target loci can effectively increase the range of inference for chemical risk assessments from tested to untested and untestable species.

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.

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.

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect.

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts that adaptation is due to a few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and using RAD-seq genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that one to two large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling. One QTL locus was shared across all populations and another was shared across three populations. One QTL locus showed strong signatures of recent natural selection in the corresponding wild population but another QTL locus did not. Some candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

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.

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect.

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts the influence of few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that a few large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling, where some (but not all) of these QTL loci were shared across all populations, and some (but not all) of these loci showed signatures of recent natural selection in the corresponding wild population. Some strong candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

Sequence and in vitro function of chicken, ring-necked pheasant, and Japanese quail AHR1 predict in vivo sensitivity to dioxins.

There are large differences in sensitivity to the toxic and biochemical effects of dioxins and dioxin-like compounds (DLCs) among vertebrates. Previously, we demonstrated that the difference in sensitivity between domestic chicken (Gallus gallus domesticus) and common tern (Sterna hirundo) to aryl hydrocarbon receptor 1 (AHR1)-dependent changes in gene expression following exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is based upon the identities of the amino acids at two sites within the ligand binding domain of AHR1 (chicken--highly sensitive; Ile324_Ser380 vs common tern--250-fold less sensitive than chicken; Val325_Ala381). Here, we tested the hypotheses that (i) the sensitivity of other avian species to TCDD, 2,3,4,7,8-pentachlorodibenzofuran (PeCDF), and 2,3,7,8-tetrachlorodibenzofuran (TCDF) is also determined by the amino acids at sites that are equivalent to sites 324 and 380 in chicken, and (ii) Ile324_Ala380 and Val324_Ser380 genotypes confer intermediate sensitivity to DLCs in birds. We compared ligand-induced transactivation function of full-length AHR1s from chicken, common tern, ring-necked pheasant (Phasianus colchicus; Ile324_Ala380) and Japanese quail (Coturnix japonica; Val324_Ala380), and three Japanese quail AHR1 mutants. The results support our hypothesis that avian species can be grouped into three general classes of sensitivity to DLCs. Both AHR1 genotype and in vitro transactivation assays predict in vivo sensitivity. Contrary to the assumption that TCDD is the most potent DLC, PeCDF was more potent than TCDD at activating Japanese quail (13- to 26-fold) and common tern (23- to 30-fold) AHR1. Our results support and expand previous in vitro and in vivo work that demonstrated ligand-dependent species differences in AHR1 affinity. The findings and methods will be of use for DLC risk assessments.

Nematostella vectensis exhibits an enhanced molecular stress response upon co-exposure to highly weathered oil and surface UV radiation.

Crude oil released into the environment undergoes weathering processes that gradually change its composition and toxicity. Co-exposure to petroleum mixtures and other stressors, including ultraviolet (UV) radiation, may lead to synergistic effects and increased toxicity. Laboratory studies should consider these factors when testing the effects of oil exposure on aquatic organisms. Here, we study transcriptomic responses of the estuarine sea anemone Nematostella vectensis to naturally weathered oil, with or without co-exposure to environmental levels of UV radiation. We find that co-exposure greatly enhances the response. We use bioinformatic analyses to identify molecular pathways implicated in this response, which suggest phototoxicity and oxidative damage as mechanisms for the enhanced stress response. Nematostella's stress response shares similarities with the vertebrate oxidative stress response, implying deep conservation of certain stress pathways in animals. We show that exposure to weathered oil along with surface-level UV exposure has substantial physiological consequences in a model cnidarian.

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.

PCB126 Exposure Revealed Alterations in m6A RNA Modifications in Transcripts Associated With AHR Activation.

Chemical modifications of proteins, DNA, and RNA moieties play critical roles in regulating gene expression. Emerging evidence suggests the RNA modifications (epitranscriptomics) have substantive roles in basic biological processes. One of the most common modifications in mRNA and noncoding RNAs is N6-methyladenosine (m6A). In a subset of mRNAs, m6A sites are preferentially enriched near stop codons, in 3' UTRs, and within exons, suggesting an important role in the regulation of mRNA processing and function including alternative splicing and gene expression. Very little is known about the effect of environmental chemical exposure on m6A modifications. As many of the commonly occurring environmental contaminants alter gene expression profiles and have detrimental effects on physiological processes, it is important to understand the effects of exposure on this important layer of gene regulation. Hence, the objective of this study was to characterize the acute effects of developmental exposure to PCB126, an environmentally relevant dioxin-like PCB, on m6A methylation patterns. We exposed zebrafish embryos to PCB126 for 6 h starting from 72 h post fertilization and profiled m6A RNA using methylated RNA immunoprecipitation followed by sequencing (MeRIP-seq). Our analysis revealed 117 and 217 m6A peaks in the DMSO and PCB126 samples (false discovery rate 5%), respectively. The majority of the peaks were preferentially located around the 3' UTR and stop codons. Statistical analysis revealed 15 m6A marked transcripts to be differentially methylated by PCB126 exposure. These include transcripts that are known to be activated by AHR agonists (eg, ahrra, tiparp, nfe2l2b) as well as others that are important for normal development (vgf, cebpd, sned1). These results suggest that environmental chemicals such as dioxin-like PCBs could affect developmental gene expression patterns by altering m6A levels. Further studies are necessary to understand the functional consequences of exposure-associated alterations in m6A levels.

Developmental Exposure to PCB153 (2,2',4,4',5,5'-Hexachlorobiphenyl) Alters Circadian Rhythms and the Expression of Clock and Metabolic Genes.

Polychlorinated biphenyls (PCBs) are highly persistent and ubiquitously distributed environmental pollutants. Based on their chemical structure, PCBs are classified into non-ortho-substituted and ortho-substituted congeners. Non-ortho-substituted PCBs are structurally similar to dioxin and their toxic effects and mode of action are well-established. In contrast, very little is known about the effects of ortho-substituted PCBs, particularly, during early development. The objective of this study is to investigate the effects of exposure to an environmentally prominent ortho-substituted PCB (2,2',4,4',5,5'-hexachlorobiphenyl; PCB153) on zebrafish embryos. We exposed zebrafish embryos to 3 different concentrations of PCB153 starting from 4 to 120 hours post-fertilization (hpf). We quantified gross morphological changes, behavioral phenotypes, gene expression changes, and circadian behavior in the larvae. There were no developmental defects during the exposure period, but starting at 7 dpf, we observed spinal deformity in the 10 µM PCB153 treated group. A total of 633, 2227, and 3378 differentially expressed genes were observed in 0.1 µM (0.036 µg/ml), 1 µM (0.36 µg/ml), and 10 µM (3.6 µg/ml) PCB153-treated embryos, respectively. Of these, 301 genes were common to all treatment groups. KEGG pathway analysis revealed enrichment of genes related to circadian rhythm, FoxO signaling, and insulin resistance pathways. Behavioral analysis revealed that PCB153 exposure significantly alters circadian behavior. Disruption of circadian rhythms has been associated with the development of metabolic and neurological diseases. Thus, understanding the mechanisms of action of environmental chemicals in disrupting metabolism and other physiological processes is essential.

Developmental Regulation of Nuclear Factor Erythroid-2 Related Factors (nrfs) by AHR1b in Zebrafish (Danio rerio).

Interactions between regulatory pathways allow organisms to adapt to their environment and respond to stress. One interaction that has been recently identified occurs between the aryl hydrocarbon receptor (AHR) and the nuclear factor erythroid-2 related factor (NRF) family. Each transcription factor regulates numerous downstream genes involved in the cellular response to toxicants and oxidative stress; they are also implicated in normal developmental pathways. The zebrafish model was used to explore the role of AHR regulation of nrf genes during development and in response to toxicant exposure. To determine if AHR1b is responsible for transcriptional regulation of 6 nrf genes during development, a loss-of-function experiment using morpholino-modified oligonucleotides was conducted followed by a chromatin immunoprecipitation study at the beginning of the pharyngula period (24 h postfertilization). The expression of nrf1a was AHR1b dependent and its expression was directly regulated through specific XREs in its cis-promoter. However, nrf1a expression was not altered by exposure to 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin (TCDD), a toxicant and prototypic AHR agonist. The expression of nrf1b, nrf2a, and nfe2 was induced by TCDD, and AHR1b directly regulated their expression by binding to cis-XRE promoter elements. Last, nrf2b and nrf3 were neither induced by TCDD nor regulated by AHR1b. These results show that AHR1b transcriptionally regulates nrf genes under toxicant modulation via binding to specific XREs. These data provide a better understanding of how combinatorial molecular signaling potentially protects embryos from embryotoxic events following toxicant exposure.

Altered lipid homeostasis in a PCB-resistant Atlantic killifish (Fundulus heteroclitus) population from New Bedford Harbor, MA, U.S.A.

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