A multidimensional density dependent matrix population model for assessing risk of stressors to fish populations.

Modeling exposure and recovery of fish and wildlife populations after stressor mitigation serves as a basis for evaluating remediation success. Herein, we develop a novel multidimensional density dependent matrix population model that analyzes both size-structure and age class-structure simultaneously. This modeling approach emphasizes application in conjunction with field monitoring efforts (e.g., through effects-based monitoring programs) and/or laboratory analysis to link effects due to stressors to outcomes in populations. We applied the model to investigate Atlantic killifish (Fundulus heteroclitus) exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin with effects on fertility and survival rates. The Atlantic killifish is an important and well-studied model organism for understanding the effects of pollutants and other stressors in estuarine and marine ecosystems. For each exposure concentration, the corresponding plots of total population size, population size structure, and age structure over time were generated. The present study serves as an example of how a multidimensional matrix population model can integrate effects across the life cycle, provide a linkage between endpoints observed in the individual and ecological risk to the population as a whole, and project outcomes for multiple generations.

Developing a High-Quality Linkage Map for the Atlantic Killifish Fundulus heteroclitus.

Killifish (Fundulus heteroclitus) are widely distributed among different aquatic environments where they demonstrate an impressive range of highly-plastic and locally adaptive phenotypes. High-throughput sequencing has begun to unravel the mechanisms and evolutionary history of these interesting features by establishing relationships in the genotype-phenotype map. However, some genotype-phenotype analyses require a higher order of contiguity than what initial scaffolded (fragmented genome assembly where contigs have been assemble into scaffolds) genome assemblies can provide. Here, we used 5,685 high-quality RAD-Seq markers from a single mapping family to order 84% of the scaffolded genome assembly to 24 chromosomes. This serves to: 1) expand the killifish genomic toolkit, 2) estimate genome-wide recombination rates, and 3) compare genome synteny to humans and other fishes. After initially building our map, we found that the selection of thresholds for sequence data filtration highly impacted scaffold placement in the map. We outline each step of the approach that dramatically improved our map to help guide others toward more effective linkage mapping for genome assembly. Our final map supports strong conservation of genomic synteny among closely related fish species and reveals previously described chromosomal rearrangements between more distantly related clades. However, we also commonly found minor scaffold misorientations in F. heteroclitus and in other assemblies, suggesting that further mapping (such as optical mapping) is necessary for finer scale resolution of genome structure. Lastly, we discuss the problems that would be expected from misoriented/unplaced scaffolds and stress the importance of a quality mapped genome as a key feature for further investigating population and comparative genomic questions with F. heteroclitus and other taxa.

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.

Characterization of the Fundulus heteroclitus embryo transcriptional response and development of a gene expression-based fingerprint of exposure for the alternative flame retardant, TBPH (bis (2-ethylhexyl)-tetrabromophthalate).

Although alternative Flame Retardant (FR) chemicals are expected to be safer than the legacy FRs they replace, their risks to human health and the environment are often poorly characterized. This study used a small volume, fish embryo system to reveal potential mechanisms of action and diagnostic exposure patterns for TBPH (bis (2-ethylhexyl)-tetrabromophthalate), a component of several widely-used commercial products. Two different concentration of TBPH were applied to sensitive early life stages of an ecologically important test species, Fundulus heteroclitus (Atlantic killifish), with a well-annotated genome. Exposed fish embryos were sampled for transcriptomics or chemical analysis of parent compound and primary metabolite or observed for development and survival through larval stage. Global transcript profiling using RNA-seq was conducted (n = 16 per treatment) to provide a non-targeted and statistically robust approach to characterize TBPH gene expression patterns. Transcriptomic analysis revealed a dose-response in the expression of genes associated with a surprisingly limited number of biological pathways, but included the aryl hydrocarbon receptor signal transduction pathway, which is known to respond to several toxicologically-important chemical classes. A transcriptional fingerprint using Random Forests was developed that was able to perfectly discriminate exposed vs. non-exposed individuals in test sets. These results suggest that TBPH has a relatively low potential for developmental toxicity (at least in fishes), despite concerns related to its structural similarities to endocrine disrupting chemicals and that the early life stage Fundulus system may provide a convenient test system for exposure characterization. More broadly, this study advances the usefulness of a biological testing and analysis system utilizing non-targeted transcriptomics profiling and early developmental endpoints that complements current screening methods to characterize chemicals of ecological and human health concern.

Stable genetic structure and connectivity in pollution-adapted and nearby pollution-sensitive populations of Fundulus heteroclitus.

Populations of the non-migratory estuarine fish Fundulus heteroclitus inhabiting the heavily polluted New Bedford Harbour (NBH) estuary have shown inherited tolerance to local pollutants introduced to their habitats in the past 100 years. Here we examine two questions: (i) Is there pollution-driven selection on the mitochondrial genome across a fine geographical scale? and (ii) What is the pattern of migration among sites spanning a strong pollution gradient? Whole mitochondrial genomes were analysed for 133 F. heteroclitus from seven nearby collection sites: four sites along the NBH pollution cline (approx. 5 km distance), which had pollution-adapted fish, as well as one site adjacent to the pollution cline and two relatively unpolluted sites about 30 km away, which had pollution-sensitive fish. Additionally, we used microsatellite analyses to quantify genetic variation over three F. heteroclitus generations in both pollution-adapted and sensitive individuals collected from two sites at two different time points (1999/2000 and 2007/2008). Our results show no evidence for a selective sweep of mtDNA in the polluted sites. Moreover, mtDNA analyses revealed that both pollution-adapted and sensitive populations harbour similar levels of genetic diversity. We observed a high level of non-synonymous mutations in the most polluted site. This is probably associated with a reduction in Ne and concomitant weakening of purifying selection, a demographic expansion following a pollution-related bottleneck or increased mutation rates. Our demographic analyses suggest that isolation by distance influences the distribution of mtDNA genetic variation between the pollution cline and the clean populations at broad spatial scales. At finer scales, population structure is patchy, and neither spatial distance, pollution concentration or pollution tolerance is a good predictor of mtDNA variation. Lastly, microsatellite analyses revealed stable population structure over the last decade.

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.

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.

Heritable oxidative phosphorylation differences in a pollutant resistant Fundulus heteroclitus population.

Populations can adapt to stress including recent anthropogenic pollution. Our published data suggests heritable differences in hepatocyte oxidative phosphorylation (OxPhos) metabolism in field-caught killifish (Fundulus heteroclitus) from the highly polluted Elizabeth River, VA, USA, relative to fish from a nearby, relatively unpolluted reference site in King's Creek VA. Consistent with other studies showing that Elizabeth River killifish are resistant to some of the toxic effects of certain contaminants, OxPhos measurements in hepatocytes from field-caught King's Creek but not field-caught Elizabeth River killifish were altered by acute benzo [a] pyrene exposures. To more definitively test whether the enhanced OxPhos metabolism and toxicity resistance are heritable, we measured OxPhos metabolism in a laboratory-reared F3 generation from the Elizabeth River population versus a laboratory-reared F1 generation from the King's Creek population and compared these results to previous data from the field-caught fish. The F3 Elizabeth River fish compared to F1 King's Creek fish had significantly higher State 3 respiration (routine metabolism) and complex II activity, and significantly lower complex I activity. The consistently higher routine metabolism in the F3 and field-caught Elizabeth River fish versus F1 and field-caught King's Creek fish implies a heritable change in OxPhos function. The observation that LEAK, E-State, Complex I and Complex II were different in laboratory bred versus field-caught fish suggests that different physiological mechanisms produce the enhanced OxPhos differences. Finally, similar to field-caught Elizabeth River fish, acute benzo [a] pyrene exposure did not affect OxPhos function of the laboratory-reared F3 generation, supporting the heritability of the toxicity resistance. Overall, these results suggest that the Elizabeth River population has evolved genetic changes in physiological homeostasis that enhance routine metabolism, and we speculate that these genetic changes interact with environmental factors altering the physiological mechanisms (e.g., alter LEAK, Complex I, and electron transfer system capacity) used to achieve this enhanced metabolism.

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.

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.

Effects of the dioxin-like PCB 126 on larval summer flounder (Paralichthys dentatus).

Little is known about the sensitivity of teleost post-embryonic developmental stages (larval and metamorphic) to dioxin-like compounds. Larval and metamorphosing summer flounder (Paralichthys dentatus) were exposed to the dioxin-like polychlorinated biphenyl congener PCB 126, to compare their sensitivity to other fish species early life stages, and to document effects on metamorphic development, including degree of eye migration and gastric maturation. Median lethal doses (LD 50s) ranged between 30 and 220 ng/g wet mass, indicating that pre- and early-metamorphic stages of summer flounder are equally sensitive to the embryos of some of the most vulnerable fish species tested. Consistent with the presence of a functional aryl hydrocarbon receptor pathway, dose-dependent induction of cytochrome P-4501A (CYP1A) at four days post-exposure was observed in liver, stomach, intestine, and kidney of metamorphosing larvae. Stage-dependent differences in the epithelial distribution of CYP1A immunoreactivity were observed in the developing stomach of fish exposed to relatively high PCB 126 doses. A single sublethal dose (15 ng/g) delayed metamorphic progress (determined by the degree of eye migration), and resulted in abnormally high levels of cell proliferation and abnormal gastric gland morphology in late metamorphic stages. These results suggest that the post-embryonic larval and metamorphic stages of summer flounder, and potentially other fish species with complex life histories, are vulnerable to the effects of dioxin-like compounds, including lethality, developmental delay, and malformations.

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.

Fundulus as the premier teleost model in environmental biology: opportunities for new insights using genomics.

A strong foundation of basic and applied research documents that the estuarine fish Fundulus heteroclitus and related species are unique laboratory and field models for understanding how individuals and populations interact with their environment. In this paper we summarize an extensive body of work examining the adaptive responses of Fundulus species to environmental conditions, and describe how this research has contributed importantly to our understanding of physiology, gene regulation, toxicology, and ecological and evolutionary genetics of teleosts and other vertebrates. These explorations have reached a critical juncture at which advancement is hindered by the lack of genomic resources for these species. We suggest that a more complete genomics toolbox for F. heteroclitus and related species will permit researchers to exploit the power of this model organism to rapidly advance our understanding of fundamental biological and pathological mechanisms among vertebrates, as well as ecological strategies and evolutionary processes common to all living organisms.

Genetic diversity and structure of an estuarine fish (Fundulus heteroclitus) indigenous to sites associated with a highly contaminated urban harbor.

Intense selection on isolated populations can cause loss of genetic diversity, which if persistent, reduces adaptive potential and increases extinction probability. Phenotypic evidence of inherited tolerance suggests that polychlorinated biphenyls (PCBs), have acted as strong selective agents on populations of a non-migratory fish, Fundulus heteroclitus, indigenous to heavily contaminated sites. To evaluate population genetic structure and test for effects of intense, multi-generational PCB contamination on genetic diversity, we used AFLP analysis on fish collected from six sites along the east coast of North America that varied widely in PCB contamination. The sites included a heavily contaminated urban harbor (New Bedford, MA), an adjacent moderately contaminated sub-estuary (Buzzards Bay, MA), and an uncontaminated estuary 60 km away (Narragansett, RI). AFLP markers distinguished populations at moderate and small scales, suggesting genetic differentiation at distances of 2 km or less. Genetic diversity did not differ across the study sites. Genome-wide diversity may have been preserved because of large effective population sizes and/or because the mechanism for genetic adaptation to these contaminants affected only a small number of loci. Alternatively, loss in diversity may have been restored with moderate levels of migration and relatively short generation time for this species.

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.

Cytochrome P4501A (CYP1A) in killifish (Fundulus heteroclitus): heritability of altered expression and relationship to survival in contaminated sediments.

Previous research has shown that killifish (Fundulus heteroclitus) inhabiting a creosote-contaminated site on the Elizabeth River in Virginia exhibit little induction of cytochrome P4501A (CYP1A) protein expression and activity upon exposure to typical CYP1A-inducing chemicals. We characterized the CYP1A response of first, second, and third generation laboratory-raised offspring of feral Elizabeth River killifish to exposure to sediments from the contaminated site as well as the prototypical polycyclic aromatic hydrocarbon (PAH)-type CYP1A inducers beta-naphthoflavone (BNF) and 3-methylcholanthrene (3-MC). The Elizabeth River offspring's responses were compared to those of offspring of killifish from two reference sites (King's Creek, Virginia, and Russell Creek, North Carolina). As with feral Elizabeth River killifish, the first generation embryos and larvae were refractory to CYP1A induction. However, the response observed in 3-year-old first generation adults, as well as with second and third generation fish, was much closer to that observed in reference-site fish. We suggest that the pattern of altered CYP1A response in Elizabeth River killifish, while persistent and heritable for one generation, is mostly nongenetically based. Additionally, we investigated the hypothesis that low CYP1A activity (measured as in ovo EROD activity) would correlate to increased survival in Elizabeth River sediment pore water; this hypothesis was not supported by our results.

Effects of benzo[a]pyrene exposure on a fish population resistant to the toxic effects of dioxin-like compounds.

Effects of a model polycyclic aromatic hydrocarbon (PAH) were compared in populations of the estuarine fish Fundulus heteroclitus indigenous to a reference site and one highly contaminated with polychlorinated biphenyls (PCBs) and other compounds. The fish population resident to the PCB-contaminated site is genetically resistant to those PCB congeners categorized as dioxin-like compounds (DLCs) that act through the aryl hydrocarbon receptor (AHR). In response to DLC exposures, these DLC-resistant fish showed poor inducibility for enzymes known to be regulated by the AHR pathway and important for the metabolism of xenobiotics including some PAHs that also act as AHR agonists. Therefore, a laboratory study using the model PAH, benzo[a]pyrene (BaP), was conducted to evaluate how PAHs might affect these wild fish populations that differed in their inherent sensitivities to DLCs and in their tissue concentrations of contaminants. Following BaP treatment, the activities of two xenobiotic metabolizing enzymes and the concentrations of BaP-DNA adducts, as measured using the 32P-postlabeling method, were lower in the livers of DLC-resistant than reference fish. These results suggest that DLC-resistance could provide protection following chronic exposures to PAHs from the long-term consequences of DNA adduct formation, such as cancer. Alternatively, reduced metabolism and elimination of toxic or photo-activated PAHs could have acute consequences to the health and reproduction of these DLC-resistant fish and their progeny. These fish populations provide useful models to evaluate the potential costs and benefits of genetic adaptation in wildlife populations subject to anthropogenic stress.