Bioaccumulation and effects of dietary exposure to the alternative flame retardant, bis(2-ethylhexyl) tetrabromophthalate (TBPH), in the Atlantic killifish, Fundulus heteroclitus.

Bis(2-ethylhexyl) tetrabromophthalate (TBPH), a high production volume flame retardant chemical used as a replacement for banned flame retardants, has been detected in media and human and wildlife tissues globally. We describe bioaccumulation and biological effects from dietary exposure of TBPH to an estuarine fish, Atlantic killifish, Fundulus heteroclitus. Briefly, adult fish were fed carrier control or chemically amended diets for 28 d, followed by 14 d of control diet feeding. Diets were amended with TBPH (TBPH_LO diet, 139 µg/g dry wt, or TBPH_HI diet, 4360 µg/g dry wt) or a polychlorinated biphenyl congener (PCB153 diet, 13 µg/g dry wt), which was included as a positive control for bioaccumulation. Although bioaccumulation of either chemical correlated with fish size, only a small proportion of the TBPH offered (<0.5% total TBPH) had bioaccumulated into TBPH-treated fish by 28 d. In contrast, 24.5% of the PCB153 offered was accounted for in 28-d PCB-treated fish. Although 28-d bioaccumulated concentrations of TBPH differed by sex and treatment, sexes did not differ in their rates of TBPH bioaccumulation, and the time to achieve 50% of 28 d concentration (T1/2 ) was estimated to be 13 d. Depuration rates of TBPH did not differ by sex or treatment, and the time after exposure to achieve T1/2 was estimated to be 22 d. Independent of treatment, male fish grew faster than female fish, but for both sexes reproductive condition (gonadal somatic index) declined unexpectedly over the experimental period. Across treatments, only the TBPH_LO treatment affected growth, reducing male but increasing female growth rates by small amounts relative to respective controls. In summary, our study used very high concentrations of dietary TBPH to contaminate fish tissues above the highest levels reported to date in wild biota, yet we observed few adverse biological effects. Environ Toxicol Chem 2018;37:2350-2360. © 2018 SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Genetic basis for rapidly evolved tolerance in the wild: adaptation to toxic pollutants by an estuarine fish species.

Atlantic killifish (Fundulus heteroclitus) residing in some urban and industrialized estuaries of the US eastern seaboard demonstrate recently evolved and extreme tolerance to toxic aryl hydrocarbon pollutants, characterized as dioxin-like compounds (DLCs). Here, we provide an unusually comprehensive accounting (69%) through quantitative trait locus (QTL) analysis of the genetic basis for DLC tolerance in killifish inhabiting an urban estuary contaminated with PCB congeners, the most toxic of which are DLCs. Consistent with mechanistic knowledge of DLC toxicity in fish and other vertebrates, the aryl hydrocarbon receptor (ahr2) region accounts for 17% of trait variation; however, QTL on independent linkage groups and their interactions have even greater explanatory power (44%). QTL interpreted within the context of recently available Fundulus genomic resources and shared synteny among fish species suggest adaptation via interacting components of a complex stress response network. Some QTL were also enriched in other killifish populations characterized as DLC-tolerant and residing in distant urban estuaries contaminated with unique mixtures of pollutants. Together, our results suggest that DLC tolerance in killifish represents an emerging example of parallel contemporary evolution that has been driven by intense human-mediated selection on natural populations.

Adaptive Significance of ERα Splice Variants in Killifish (Fundulus heteroclitus) Resident in an Estrogenic Environment.

The possibility that chronic, multigenerational exposure to environmental estrogens selects for adaptive hormone-response phenotypes is a critical unanswered question. Embryos/larvae of killifish from an estrogenic-polluted environment (New Bedford Harbor, MA [NBH]) compared with those from a reference site overexpress estrogen receptor alpha (ERα) mRNA but are hyporesponsive to estradiol. Analysis of ERα mRNAs in the two populations revealed differences in splicing of the gene encoding ERα (esr1). Here we tested the transactivation functions of four differentially expressed ERα mRNAs and tracked their association with the hyporesponsive phenotype for three generations after transfer of NBH parents to a clean environment. Deletion variants ERαΔ6 and ERαΔ6-8 were specific to NBH killifish, had dominant negative functions in an in vitro reporter assay, and were heritable. Morpholino-mediated induction of ERαΔ6 mRNA in zebrafish embryos verified its role as a dominant negative ER on natural estrogen-responsive promoters. Alternate long (ERαL) and short (ERαS) 5'-variants were similar transcriptionally but differed in estrogen responsiveness (ERαS ≫ ERαL). ERαS accounted for high total ERα expression in first generation (F1) NBH embryos/larvae but this trait was abolished by transfer to clean water. By contrast, the hyporesponsive phenotype of F1 NBH embryos/larvae persisted after long-term laboratory holding but reverted to a normal or hyper-responsive phenotype after two or three generations, suggesting the acquisition of physiological or biochemical traits that compensate for ongoing expression of negative-acting ERαΔ6 and ERαΔ6-8 isoforms. We conclude that a heritable change in the pattern of alternative splicing of ERα pre-mRNA is part of a genetic adaptive response to estrogens in a polluted environment.

Cloning of multiple ERα mRNA variants in killifish (Fundulus heteroclitus), and differential expression by tissue type, stage of reproduction, and estrogen exposure in fish from polluted and unpolluted environments.

To test the hypothesis that alternative splicing could be an adaptive mechanism for populations subject to multi-generational estrogenic exposures, we compared estrogen receptor alpha (ERα) splicing variants in two populations of killifish (Fundulus heteroclitus): one resident in an estrogenic polluted environment (New Bedford Harbor, NBH, MA, USA) and one from a relatively uncontaminated reference site (Scorton Creek, SC, MA, USA). In total we identified 19 ERα variants, each with deletions of one or more coding exons. Four of the variants with potential functional relevance were analyzed by qPCR to test for population differences in expression by tissue type, site, sex, seasonal reproductive status and estrogen treatment. Significantly, a 5'-truncated short form variant (ERαS) was highly expressed in liver and ovary, and was associated with seasonal reproductive activity in SC but not NBH fish. Both ERαS and the full-length long variant (ERαL) were estrogen-inducible (ERαS>ERαL) but the induction response was lower in NBH than in SC fish. In contrast, NBH killifish were hyper-responsive to estrogen as measured by expression of two other estrogen responsive genes: vitellogenin (Vtg) and aromatase B (AroB). Most strikingly, two ERα deletion variants (Δ6 and Δ6-8), lacking ligand binding and activation function domains, were identified in a subset of NBH fish, where they were associated with reduced responsiveness to estrogen treatment. Together, these results support the hypothesis that alternative splicing of the esr1 gene of killifish could be an autoregulatory mechanism by which estrogen modulates the differential expression of ERα, and suggests a novel and adaptive mechanistic response to xenoestrogenic exposure.

Regulation of pregnane-X-receptor, CYP3A and P-glycoprotein genes in the PCB-resistant killifish (Fundulus heteroclitus) population from New Bedford Harbor.

Killifish survive and reproduce in the New Bedford Harbor (NBH) in Massachusetts (MA), USA, a site severely contaminated with polychlorinated biphenyls (PCBs) for decades. Levels of 22 different PCB congeners were analyzed in liver from killifish collected in 2008. Concentrations of dioxin-like PCBs in liver of NBH killifish were ∼400 times higher, and the levels of non-dioxin-like PCBs ∼3000 times higher than in killifish from a reference site, Scorton Creek (SC), MA. The NBH killifish are known to be resistant to the toxicity of dioxin-like compounds and to have a reduced aryl hydrocarbon receptor (AhR) signaling response. Little is known about the responses of these fish to non-dioxin-like PCBs, which are at extraordinarily high levels in NBH fish. In mammals, some non-dioxin-like PCB congeners act through nuclear receptor 1I2, the pregnane-X-receptor (PXR). To explore this pathway in killifish, a PXR cDNA was sequenced and its molecular phylogenetic relationship to other vertebrate PXRs was determined. Killifish were also collected in 2009 from NBH and SC, and after four months in the laboratory they were injected with a single dose of either the dioxin-like PCB 126 (an AhR agonist) or the non-dioxin-like PCB 153 (a mammalian PXR agonist). Gills and liver were sampled three days after injection and transcript levels of genes encoding PXR, cytochrome P450 3A (CYP3A), P-glycoprotein (Pgp), AhR2 and cytochrome P450 1A (CYP1A) were measured by quantitative PCR. As expected, there was little effect of PCB exposure on mRNA expression of AhR2 or CYP1A in liver and gills of NBH fish. In NBH fish, but not in SC fish, there was increased mRNA expression of hepatic PXR, CYP3A and Pgp upon exposure to either of the two PCB congeners. However, basal PXR and Pgp mRNA levels in liver of NBH fish were significantly lower than in SC fish. A different pattern was seen in gills, where there were no differences in basal mRNA expression of these genes between the two populations. In SC fish, but not in NBH fish, there was increased mRNA expression of branchial PXR and CYP3A upon exposure to PCB126 and of CYP3A upon exposure to PCB153. The results suggest a difference between the two populations in non-AhR transcription factor signaling in liver and gills, and that this could involve killifish PXR. It also implies possible cross-regulatory interactions between that factor (presumably PXR) and AhR2 in liver of these fish.

Expression and function of ryanodine receptor related pathways in PCB tolerant Atlantic killifish (Fundulus heteroclitus) from New Bedford Harbor, MA, USA.

Atlantic killifish (Fundulus heteroclitus) thrive in New Bedford Harbor (NBH), MA, highly contaminated with polychlorinated biphenyls (PCBs). Resident killifish have evolved tolerance to dioxin-like (DL) PCBs, whose toxic effects through the aryl hydrocarbon receptor (AhR) are well studied. In NBH, non-dioxin like PCBs (NDL PCBs), which lack activity toward the AhR, vastly exceed levels of DL congeners yet how killifish counter NDL toxic effects has not been explored. In mammals and fish, NDL PCBs are potent activators of ryanodine receptors (RyR), Ca(2+) release channels necessary for a vast array of physiological processes. In the current study we compared the expression and function of RyR related pathways in NBH killifish with killifish from the reference site at Scorton Creek (SC, MA). Relative to the SC fish, adults from NBH displayed increased levels of skeletal muscle RyR1 protein, and increased levels of FK506-binding protein 12 kDa (FKBP12) an accessory protein essential for NDL PCB-triggered changes in RyR channel function. In accordance with increased RyR1 levels, NBH killifish displayed increased maximal ligand binding, increased maximal response to Ca(2+) activation and increased maximal response to activation by the NDL PCB congener PCB 95. Compared to SC, NBH embryos and larvae had increased levels of mtor and ryr2 transcripts at multiple stages of development, and generations, while levels of serca2 were decreased at 9 days post-fertilization in the F1 and F2 generations. These findings suggest that there are compensatory and heritable changes in RyR mediated Ca(2+) signaling proteins or potential signaling partners in NBH killifish.

Targeted mutagenesis of aryl hydrocarbon receptor 2a and 2b genes in Atlantic killifish (Fundulus heteroclitus).

Understanding molecular mechanisms of toxicity is facilitated by experimental manipulations, such as disruption of function by gene targeting, that are especially challenging in non-standard model species with limited genomic resources. While loss-of-function approaches have included gene knock-down using morpholino-modified oligonucleotides and random mutagenesis using mutagens or retroviruses, more recent approaches include targeted mutagenesis using zinc finger nuclease (ZFN), transcription activator-like effector nuclease (TALENs) and clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology. These latter methods provide more accessible opportunities to explore gene function in non-traditional model species. To facilitate evaluation of toxic mechanisms for important categories of aryl hydrocarbon pollutants, whose actions are known to be receptor mediated, we used ZFN and CRISPR-Cas9 approaches to generate aryl hydrocarbon receptor 2a (AHR2a) and AHR2b gene mutations in Atlantic killifish (Fundulus heteroclitus) embryos. This killifish is a particularly valuable non-traditional model, with multiple paralogs of AHR whose functions are not well characterized. In addition, some populations of this species have evolved resistance to toxicants such as halogenated aromatic hydrocarbons. AHR-null killifish will be valuable for characterizing the role of the individual AHR paralogs in evolved resistance, as well as in normal development. We first used five-finger ZFNs targeting exons 1 and 3 of AHR2a. Subsequently, CRISPR-Cas9 guide RNAs were designed to target regions in exon 2 and 3 of AHR2a and AHR2b. We successfully induced frameshift mutations in AHR2a exon 3 with ZFN and CRISPR-Cas9 guide RNAs, with mutation frequencies of 10% and 16%, respectively. In AHR2b, mutations were induced using CRISPR-Cas9 guide RNAs targeting sites in both exon 2 (17%) and exon 3 (63%). We screened AHR2b exon 2 CRISPR-Cas9-injected embryos for off-target effects in AHR paralogs. No mutations were observed in closely related AHR genes (AHR1a, AHR1b, AHR2a, AHRR) in the CRISPR-Cas9-injected embryos. Overall, our results demonstrate that targeted genome-editing methods are efficient in inducing mutations at specific loci in embryos of a non-traditional model species, without detectable off-target effects in paralogous genes.

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.

Targeted approach to identify genetic loci associated with evolved dioxin tolerance in Atlantic killifish (Fundulus heteroclitus).

BACKGROUND: The most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed. RESULTS: Since the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs. CONCLUSIONS: The observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.

Common mechanism underlies repeated evolution of extreme pollution tolerance.

Human alterations to the environment can exert strong evolutionary pressures, yet contemporary adaptation to human-mediated stressors is rarely documented in wildlife populations. A common-garden experimental design was coupled with comparative transcriptomics to discover evolved mechanisms enabling three populations of killifish resident in urban estuaries to survive normally lethal pollution exposure during development, and to test whether mechanisms are unique or common across populations. We show that killifish populations from these polluted sites have independently converged on a common adaptive mechanism, despite variation in contaminant profiles among sites. These populations are united by a similarly profound desensitization of aryl-hydrocarbon receptor-mediated transcriptional activation, which is associated with extreme tolerance to the lethal effects of toxic dioxin-like pollutants. The rapid, repeated, heritable and convergent nature of evolved tolerance suggests that ancestral killifish populations harboured genotypes that enabled adaptation to twentieth-century industrial pollutants.

The effects of mitochondrial genotype on hypoxic survival and gene expression in a hybrid population of the killifish, Fundulus heteroclitus.

The physiological link between oxygen availability and mitochondrial function is well established. However, whether or not fitness variation is associated with mitochondrial genotypes in the field remains a contested topic in evolutionary biology. In this study, we draw on a population of the teleost fish, Fundulus heteroclitus, where functionally distinct subspecies hybridize, likely as a result of past glacial events. We had two specific aims: (i) to determine the effect of mtDNA genotype on survivorship of male and female fish under hypoxic stress and (ii) to determine the effect of hypoxic stress, sex and mtDNA genotype on gene expression. We found an unexpected and highly significant effect of sex on survivorship under hypoxic conditions, but no significant effect of mtDNA genotype. Gene expression analyses revealed hundreds of transcripts differentially regulated by sex and hypoxia. Mitochondrial transcripts and other predicted pathways were among those influenced by hypoxic stress, and a transcript corresponding to the mtDNA control region was the most highly suppressed transcript under the conditions of hypoxia. An RT-PCR experiment on the control region was consistent with microarray results. Effects of mtDNA sequence variation on genome expression were limited; however, a potentially important epistasis between mtDNA sequence and expression of a nuclear-encoded mitochondrial translation protein was discovered. Overall, these results confirm that mitochondrial regulation is a major component of hypoxia tolerance and further suggest that purifying selection has been the predominant selective force on mitochondrial genomes in these two subspecies.

Population genetic diversity and fitness in multiple environments.

BACKGROUND: When a large number of alleles are lost from a population, increases in individual homozygosity may reduce individual fitness through inbreeding depression. Modest losses of allelic diversity may also negatively impact long-term population viability by reducing the capacity of populations to adapt to altered environments. However, it is not clear how much genetic diversity within populations may be lost before populations are put at significant risk. Development of tools to evaluate this relationship would be a valuable contribution to conservation biology. To address these issues, we have created an experimental system that uses laboratory populations of an estuarine crustacean, Americamysis bahia with experimentally manipulated levels of genetic diversity. We created replicate cultures with five distinct levels of genetic diversity and monitored them for 16 weeks in both permissive (ambient seawater) and stressful conditions (diluted seawater). The relationship between molecular genetic diversity at presumptive neutral loci and population vulnerability was assessed by AFLP analysis. RESULTS: Populations with very low genetic diversity demonstrated reduced fitness relative to high diversity populations even under permissive conditions. Population performance decreased in the stressful environment for all levels of genetic diversity relative to performance in the permissive environment. Twenty percent of the lowest diversity populations went extinct before the end of the study in permissive conditions, whereas 73% of the low diversity lines went extinct in the stressful environment. All high genetic diversity populations persisted for the duration of the study, although population sizes and reproduction were reduced under stressful environmental conditions. Levels of fitness varied more among replicate low diversity populations than among replicate populations with high genetic diversity. There was a significant correlation between AFLP diversity and population fitness overall; however, AFLP markers performed poorly at detecting modest but consequential losses of genetic diversity. High diversity lines in the stressful environment showed some evidence of relative improvement as the experiment progressed while the low diversity lines did not. CONCLUSIONS: The combined effects of reduced average fitness and increased variability contributed to increased extinction rates for very low diversity populations. More modest losses of genetic diversity resulted in measurable decreases in population fitness; AFLP markers did not always detect these losses. However when AFLP markers indicated lost genetic diversity, these losses were associated with reduced population fitness.

PCBs and DDE in tree swallow (Tachycineta bicolor) eggs and nestlings from an estuarine PCB superfund site, New Bedford Harbor, MA, U.S.A.

While breeding tree swallows (Tachycineta bicolor) have been used as biomonitors for freshwater sites, we report the first use of this species to assess contaminant bioaccumulation from estuarine breeding grounds into these aerial insectivores. Eggs and nestlings were collected from nest boxes in a polychlorinated biphenyl (PCB) contaminated estuary, the New Bedford Harbor Superfund site (NBH, Massachusetts, USA), and a reference salt marsh, Fox Hill (FH, Jamestown, Rhode Island, USA). Sediments, eggs, and nestlings were compared on a ng g(-1) wet weight basis for total PCBs and DDE (1,1-bis-(4-chlorophenyl)-2,2-dichloroethene), metabolite of DDT (1,1,1-trichloro-2,2-bis-(p-chlorophenyl)ethane). NBH samples contained high concentrations of PCBs compared to FH for sediment (36,500 and 0.2), eggs (11,200 and 323), and nestlings (16,800 and 26). PCB homologue patterns linked tree swallow contamination to NBH sediment. NBH samples were also contaminated with DDE compared to FH for sediment (207 and 0.9) and nestlings (235 and 30) but not for eggs (526 and 488), suggesting both NBH and nonbreeding ground sources for DDE. The relationships between sediment and tree swallow egg and nestling PCBs were similar to those reported for freshwater sites. Like some highly contaminated freshwater sites, NBH PCB bioaccumulation had little apparent effect on reproductive success.

Evolution of tolerance to PCBs and susceptibility to a bacterial pathogen (Vibrio harveyi) in Atlantic killifish (Fundulus heteroclitus) from New Bedford (MA, USA) harbor.

A population of the non-migratory estuarine fish Fundulus heteroclitus (Atlantic killifish) resident to New Bedford (NB), Massachusetts, USA, an urban harbor highly contaminated with polychlorinated biphenyls (PCBs), demonstrates recently evolved tolerance to some aspects of PCB toxicity. PCB toxicology, ecological theory, and some precedence supported expectations of increased susceptibility to pathogens in NB killifish. However, laboratory bacterial challenges of the marine pathogen Vibrio harveyi to wild fish throughout the reproductive season and to their mature laboratory-raised progeny demonstrated comparable survival by NB and reference killifish, and improved survival by NB males. These results are inconsistent with hypothesized trade-offs of adaptation, and suggest that evolved tolerance in NB killifish may include mechanisms that minimize the immunosuppressive effects of PCBs. Compensatory strategies of populations persisting in highly contaminated environments provide a unique perspective for understanding the long-term ecological effects of toxic chemicals.

Population genetic structure of a nonmigratory estuarine fish (Fundulus heteroclitus) across a strong gradient of polychlorinated biphenyl contamination.

Populations of the estuarine fish Fundulus heteroclitus indigenous to contaminated sites exhibit heritable resistance to some of the toxic effects of early life-stage exposure to polychlorinated biphenyls (PCBs). This evolved tolerance provides evidence of strong selection by PCBs, and it suggests other potential genetic effects of these stressors on resident populations. Environmental contaminants have the potential to affect the genetic structure of populations and to reduce genetic diversity, but species life-history traits, particularly patterns of migration and dispersal, also influence the distribution of genetic variation among populations. Therefore, the present was conducted to determine whether genetic diversity or genetic structure is altered in populations of F. heteroclitus indigenous to 18 sites in Massachusetts (USA) and Rhode Island (USA), representing a steep gradient of sediment PCB contamination and culminating in a Superfund site at New Bedford Harbor (NBH; MA, USA). Allele frequencies at enzymatic loci were used to assess genetic structure and diversity. Selection experiments using a highly toxic PCB congener (3,3',4,4',5-pentachlorobiphenyl) were conducted to determine if genetic patterns at field sites could be associated with contaminant exposures. Although allele frequencies clearly reflected a pattern of isolation by distance, the results indicated neither significant loss of genetic diversity nor alteration of allele frequencies for populations of F. heteroclitus in NBH.

Population genetic structure and tolerance to dioxin-like compounds of a migratory marine fish (Menidia menidia) at polychlorinated biphenyl-contaminated and reference sites.

The present study was conducted to evaluate evidence of genetic adaptation to local contaminants in populations of the migratory marine fish Menidia menidia residing seasonally in reference sites or an industrial harbor contaminated with dioxin-like compounds (DLCs). For this purpose, we compared DLC sensitivity and genetic patterns of populations sampled from sites both inside and outside New Bedford Harbor (NBH; MA, USA), a U.S. Environmental Protection Agency Superfund site with extreme polychlorinated biphenyl (PCB) contamination. Offspring of M. menidia collected from NBH were significantly less sensitive regarding embryonic exposure to the dioxin-like PCB congener 3,3',4,4',5-pentachlorobiphenyl (PCB 126) than offspring of M. menidia from a reference site. Analysis of 10 polymorphic enzymatic loci indicated little genetic differentiation among populations in the study area. However, genotype frequencies of juveniles from both NBH and an adjacent site in Massachusetts exhibited significant deviations from Hardy-Weinberg equilibrium expectations at one locus, phosphoglucomutase (PGM*). Genetic analysis of survivors of embryonic laboratory exposure to PCB 126 indicated that genotypes at PGM* were related to survivorship. Although a relationship was identified between DLC tolerance and PGM* genotype, regional mixing of M. menidia populations during migration and absence of multigeneration exposure at contaminated sites may limit localized adaptation.

Isolation and characterization of two cytochrome P450 aromatase forms in killifish (Fundulus heteroclitus): differential expression in fish from polluted and unpolluted environments.

Populations of killifish (Fundulus heteroclitus) persist in many different highly polluted environment indicative of adaptation or tolerance. In this study, we sought to determine whether long term, multigenerational exposures to environmental contaminants has affected reproductively relevant genes and biological processes. A homology cloning strategy was used to isolate the killifish cytochrome P450 aromatase (P450arom, estrogen synthetase) cDNAs. Consistent with previous fish studies, killifish were found to have two P450arom forms, which segregated into A- and B-gene clades and were differentially expressed in brain (B >> A) and ovary (A >> B). Comparison of killifish from highly polluted (New Bedford Harbor, NBH) and unpolluted (Scorton Creek, SC) environments revealed no site-related differences in P450arom coding sequences or in overall tissue distribution patterns. As measured by real-time quantitative PCR (QPCR) analysis, however, P450arormB (a known marker of estrogen effect) was approximately two-fold higher in the brain of NBH than of SC fish, a difference seen in reproductively active and inactive males and females. Providing further evidence of exposure to estrogen-like pollutants or metabolites in NBH, vitellogenin (vtg) mRNA and protein were elevated in seasonally active and inactive males, and in reproductively inactive females, when compared to SC fish. By contrast, during the period of reproductive activity, NBH females had a lower gonadosomatic index, lower plasma estrogen, a decreased hepatosomatic index, and reduced vtg expression as compared to SC females, indicating that the female hypothalamic-pituitary-gonadal (HPG)-liver axis is impaired in the polluted environment. As measured by a decrease in plasma androgen (but not GSI), the male HPG axis was impaired in reproductively active NBH versus SC fish. In agreement with reports that NBH killifish are resistant to dioxin-like chemicals (DLC) that activate arylhydrocarbon receptor (AhR) signaling, ovarian P450aromA (a marker of dioxin-like effect in zebrafish embryos) did not differ in SC and NBH fish. In conclusion, the killifish population at the NBH Superfund site maintains a level of reproductive competence in the face of evidence of exposure to estrogen-like pollutants and endocrine disruption.

Predicting the occurrence of genetic adaptation to dioxinlike compounds in populations of the estuarine fish Fundulus heteroclitus.

A population of the nonmigratory estuarine fish species Fundulus heteroclitus (mummichog) indigenous to a polychlorinated biphenyl (PCB)-contaminated Superfund site (New Bedford Harbor, MA, USA) demonstrated an inherited tolerance to local, dioxinlike contaminants (DLCs). These findings suggest that DLCs have acted as selective agents, allowing the survival of only the most tolerant individuals, forming DLC-adapted populations. We hypothesized that DLC-tolerant mummichog populations would reside where local conditions are toxic to sensitive individuals, and that toxic environmental conditions could be predicted based on responses of sensitive early life stages to laboratory exposures of DLCs. As a measure of DLC tolerance, progeny of field-collected fish were tested in the laboratory with a dioxinlike PCB congener, 3,3',4,4',5-pentachlorobiphenyl (PCB 126). Mummichog populations were collected from sites with sediment PCB concentrations predicted to range from nontoxic to toxic. Consistent with predictions, tolerant populations were indigenous to sites with elevated sediment PCB concentrations. Also, as predicted, DLC-tolerant populations were resident to sites far less contaminated than the Superfund site. These results suggest that exposures to persistent, bioaccumulative, and toxic contaminants may produce evolutionary effects on a geographic scale larger than previously envisioned. This study presents an approach and describes a model system that may improve understanding of the scale of occurrence for these potentially irreversible ecological effects.