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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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".

Early Life Exposure to Low Levels of AHR Agonist PCB126 (3,3',4,4',5-Pentachlorobiphenyl) Reprograms Gene Expression in Adult Brain.

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, nonembryotoxic 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 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.

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

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.6kDa (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.

Diversity as Opportunity: Insights from 600 Million Years of AHR Evolution.

The aryl hydrocarbon receptor (AHR) was for many years of interest only to pharmacologists and toxicologists. However, this protein has fundamental roles in biology that are being revealed through studies in diverse animal species. The AHR is an ancient protein. AHR homologs exist in most major groups of modern bilaterian animals, including deuterostomes (chordates, hemichordates, echinoderms) and the two major clades of protostome invertebrates [ecdysozoans (e.g. arthropods and nematodes) and lophotrochozoans (e.g. molluscs and annelids)]. AHR homologs also have been identified in cnidarians such as the sea anemone Nematostella and in the genome of Trichoplax, a placozoan. Bilaterians, cnidarians, and placozoans form the clade Eumetazoa, whose last common ancestor lived approximately 600 million years ago (MYA). The presence of AHR homologs in modern representatives of all these groups indicates that the original eumetazoan animal possessed an AHR homolog. Studies in invertebrates and vertebrates reveal parallel functions of AHR in the development and function of sensory neural systems, suggesting that these may be ancestral roles. Vertebrate animals are characterized by the expansion and diversification of AHRs, via gene and genome duplications, from the ancestral protoAHR into at least five classes of AHR-like proteins: AHR, AHR1, AHR2, AHR3, and AHRR. The evolution of multiple AHRs in vertebrates coincided with the acquisition of high-affinity binding of halogenated and polynuclear aromatic hydrocarbons and the emergence of adaptive functions involving regulation of xenobiotic-metabolizing enzymes and roles in adaptive immunity. The existence of multiple AHRs may have facilitated subfunction partitioning and specialization of specific AHR types in some taxa. Additional research in diverse model and non-model species will continue to enrich our understanding of AHR and its pleiotropic roles in biology and toxicology.

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.

The genomic landscape of rapid repeated evolutionary adaptation to toxic pollution in wild fish.

Atlantic killifish populations have rapidly adapted to normally lethal levels of pollution in four urban estuaries. Through analysis of 384 whole killifish genome sequences and comparative transcriptomics in four pairs of sensitive and tolerant populations, we identify the aryl hydrocarbon receptor-based signaling pathway as a shared target of selection. This suggests evolutionary constraint on adaptive solutions to complex toxicant mixtures at each site. However, distinct molecular variants apparently contribute to adaptive pathway modification among tolerant populations. Selection also targets other toxicity-mediating genes and genes of connected signaling pathways; this indicates complex tolerance phenotypes and potentially compensatory adaptations. Molecular changes are consistent with selection on standing genetic variation. In killifish, high nucleotide diversity has likely been a crucial substrate for selective sweeps to propel rapid adaptation.

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.

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

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.

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.