Resistance to contaminants in North American fish populations

RNA-Seq reveals complex genetic response to Deepwater Horizon oil release in Fundulus grandis

Background

The release of oil resulting from the blowout of the Deepwater Horizon (DH) drilling platform was one of the largest in history discharging more than 189 million gallons of oil and subject to widespread application of oil dispersants. This event impacted a wide range of ecological habitats with a complex mix of pollutants whose biological impact is still not yet fully understood. To better understand the effects on a vertebrate genome, we studied gene expression in the salt marsh minnow Fundulus grandis, which is local to the northern coast of the Gulf of Mexico and is a sister species of the ecotoxicological model Fundulus heteroclitus. To assess genomic changes, we quantified mRNA expression using high throughput sequencing technologies (RNA-Seq) in F. grandis populations in the marshes and estuaries impacted by DH oil release. This application of RNA-Seq to a non-model, wild, and ecologically significant organism is an important evaluation of the technology to quickly assess similar events in the future.

Results

Our de novo assembly of RNA-Seq data produced a large set of sequences which included many duplicates and fragments. In many cases several of these could be associated with a common reference sequence using blast to query a reference database. This reduced the set of significant genes to 1,070 down-regulated and 1,251 up-regulated genes. These genes indicate a broad and complex genomic response to DH oil exposure including the expected AHR-mediated response and CYP genes. In addition a response to hypoxic conditions and an immune response are also indicated. Several genes in the choriogenin family were down-regulated in the exposed group; a response that is consistent with AH exposure. These analyses are in agreement with oligonucleotide-based microarray analyses, and describe only a subset of significant genes with aberrant regulation in the exposed set.

Conclusion

RNA-Seq may be successfully applied to feral and extremely polymorphic organisms that do not have an underlying genome sequence assembly to address timely environmental problems. Additionally, the observed changes in a large set of transcript expression levels are indicative of a complex response to the varied petroleum components to which the fish were exposed.

Relationships between body burdens of trace metals (As, Cu, Fe, Hg, Mn, Se, and Zn) and the relative body size of small tooth flounder (Pseudorhombus jenynsii)

Several studies have described strong relationships between body size and the accumulation of trace metals in animal tissues. However, few of these studies have utilized aging techniques to control for age related effects. We utilized relative body size (gy(-1)) of a model flounder species, Pseudorhombus jenynsii, in order to control for age related effects on growth and size measurements. We investigated links between relative body size, concentrations of trace metals in flounder muscle tissue, physico-chemical variables (temperature, salinity, pH, and turbidity), and levels of trace metals in the sediment. Flounder were sampled using an otter trawl net in the inner areas of eight estuaries that were either heavily modified or relatively unmodified by urbanization and industrial activity. Our results indicate that this commonly eaten fish is accumulating significant levels of some trace metals in their muscle tissue, both in relatively unmodified and heavily modified estuaries. Concentrations of Cu, Zn and Fe in muscle tissue, as well as temperature, showed a negative relationship to the relative body size of flounder. In contrast, Se and Hg in muscle showed a positive relationship to relative body size. Observed growth patterns indicate that these effects are not driven by age related differences in metabolic activity. Instead, our results suggest that differences in food supply or toxicological effects may be responsible for the observed relationships between relative body size and concentrations of Cu, Zn, and Se in muscle tissues. The use of otolith aging and growth measurement techniques represents a novel method for assessing the relationships between trace metal accumulation and the relative body size of fish in a field environment.

Acclimation of life-history traits to experimental changes in environmental contaminant concentrations in brown bullhead (Ameiurus nebulosus)

One adaptive mechanism aquatic populations use to facilitate tolerance to environmental contaminants is acclimation. Polychlorinated biphenyls (PCBs) are a globally ubiquitous class of persistent organic contaminants that have been linked to reproductive impairments in fish. The authors used female brown bullhead (Ameiurus nebulosus) to test whether acclimation of reproductive life-history traits occurs in response to changes in sum PCB exposure. They compared egg diameter, gonadosomatic index (GSI), and fecundity of fish directly caught from wild populations exposed to a range of contaminant concentrations (acute), to those collected from the same populations a year before, which were placed in a clean environment to clear their contaminants throughout that year (cleared). Sum PCB concentrations were also determined for each individual. Brown bullhead from acute treatments had significantly greater sum PCB concentrations compared with cleared treatments. Egg diameter and GSI metrics were greater in cleared treatments compared with acute treatments (by 6 and 14%, respectively). Treatment effect (i.e., acute or cleared), as opposed to where the fish were collected from, accounts for 72 to 89% of the variation in the reproductive life-history trait variables. No difference in fecundity was found between acute and cleared treatments. The authors found support that acclimation of reproductive life-history traits occurs to changes in sum PCB concentration. To their knowledge, the present study is the first experimental test of acclimation responses of female life-history traits to contaminants in wild populations.

Fundulus heteroclitus adapted to PAHs are cross-resistant to multiple insecticides

Atlantic killifish (Fundulus heteroclitus) from the Atlantic Wood Superfund site on the Elizabeth River (ER), VA are dramatically resistant to the acute toxicity and teratogenesis caused by polycyclic aromatic hydrocarbons (PAHs). To understand the consequences of adaptation to chronic PAH pollution, we have attempted to further define the chemical tolerance associated with this resistance. An important component of the PAH adaptation of ER fish is the dramatic down-regulation of the aryl hydrocarbon receptor (AHR) pathway, resulting in decreased cytochrome p450 (CYP) 1 activity. Herein, we compared the susceptibility to several insecticides of ER fish to that of reference site (King's Creek; KC) fish; use of these chemicals as probes of the resistance will help to demonstrate if the contaminant adaptation exhibited by ER fish is broad or narrow and AHR-focused. We hypothesized that ER fish would be less susceptible to the organophosphate chlorpyrifos (activated by CYP) and more susceptible to the pyrethroid permethrin (detoxified by CYP). Comparison of acute toxicity in 5-day-old larvae supported this hypothesis for chlorpyrifos. As expected, chemical up-regulation of CYP by co-exposure to β-naphthoflavone (BNF) enhanced the susceptibility of KC but it did not affect ER larvae. Unexpectedly, ER larvae were much less susceptible to permethrin than KC larvae. However, co-exposure to BNF greatly decreased the susceptibility of KC larvae, indicating that metabolism of permethrin by CYP was protective. Additionally, fish from each population were compared for susceptibility to the carbamate carbaryl, an acute neurotoxicant and weak AHR agonist that induces teratogenesis similar to that caused by PAHs. ER embryos and larvae were less susceptible than KC fish. These results suggest that the adaptive phenotype of ER fish is multi-faceted and that aspects other than CYP response are likely to greatly affect their response to contaminants.

Genotoxic and reprotoxic effects of tritium and external gamma irradiation on aquatic animals

Aquatic ecosystems are chronically exposed to natural radioactivity or to artificial radionuclides released by human activities (e.g., nuclear medicine and biology,nuclear industry, military applications). Should the nuclear industry expand in the future, radioactive environmental releases, under normal operating conditions or accidental ones, are expected to increase, which raises public concerns about possible consequences on the environment and human health. Radionuclide exposures may drive macromolecule alterations, and among macromolecules DNA is the major target for ionizing radiations. DNA damage, if not correctly repaired, may induce mutations, teratogenesis, and reproductive effects. As such, damage at the molecular level may have consequences at the population level. In this review, we present an overview of the literature dealing with the effects of radionuclides on DNA, development, and reproduction of aquatic organisms. The review focuses on the main radionuclides that are released by nuclear power plants under normal operating conditions, γ emitters and tritium. Additionally, we fitted nonlinear curves to the dose-response data provided in the reviewed publications and manuscripts, and thus obtained endpoints commonly associated with ecotoxicological studies, such as the EDR(10). These were then used as a common metric for comparing the values and data published in the literature.The effects of tritium on aquatic organisms were reviewed for dose rates that ranged from 29 nGy/day to 29 Gy/day. Although beta emission from tritium decay presents a rather special risk of damage to DNA, genotoxicity-induced by tritium has been scarcely studied. Most of the effects studied have related to reproduction and development. Species sensitivity and the form of tritium present are important factors that drive the ecotoxicity of tritium. We have concluded from this review that invertebrates are more sensitive to the effects of tritium than are vertebrates.Because several calculated EDR10 values are ten times lower than background levels of γ irradiation the results of some studies either markedly call into question the adequacy of the benchmark value of 0.24 mGy/day for aquatic ecosystems that was recommended by Garnier-Laplace et al. (2006), or the dose rate estimates made in the original research, from which our EDR(10) values were derived, were under estimated, or were inadequate. For γ irradiation, the effects of several different dose rates on aquatic organisms were reviewed, and these ranged from 1 mGy/day to 18 Gy/day. DNA damage from exposure to y irradiation was studied more often than for tritium, but the major part of the literature addressed effects on reproduction and development. These data sets support the benchmark value of 0.24 mGy/day, which is recommended to protect aquatic ecosystems. RBEs, that describe the relative effectiveness of different radiation types to produce the same biological effect, were calculated using the available datasets. These RBE values ranged from 0.06 to 14.9, depending on the biological effect studied, and they had a mean of 3.1 ± 3.7 (standard deviation). This value is similar to the RBE factors of 2-3 recommended by international organizations responsible for providing guidance on radiation safety. Many knowledge gaps remain relative to the biological effects produced from exposure to tritium and y emitters. Among these are: Dose calculations: this review highlights several EDR(10) values that are below the normal range of background radiation. One explanation for this result is that dose rates were underestimated from uncertainties linked to the heterogenous distribution of tritium in cells. Therefore, the reliability of the concept of average dose to organisms must be addressed. Mechanisms of DNA DBS repair: very few studies address the most deleterious form of DNA damage, which are DNA DBSs. Future studies should focus on identifying impaired DNA DBS repair pathways and kinetics, in combination with developmental and reproductive effects. The transmission of genetic damage to offspring, which is of primary concern in the human health arena. However, there has been little work undertaken to assess the potential risk from germ cell mutagens in aquatic organisms, although this is one of the means of extrapolating effects from subcellular levels to populations. Reproductive behavior that is linked to alterations of endocrine function. Despite the importance of reproduction for population dynamics, many key endpoints were scarcely addressed within this topic. Hence, there is, to our knowledge,only one study of courtship behavior in fish exposed to γ rays, while no studies of radionuclide effects on fish endocrine function exist. Recent technical advances in the field of endocrine disrupters can be used to assess the direct or indirect effects of radionuclides on endocrine function. Identifying whether resistance to radiation effects in the field result from adaptation or acclimation mechanisms. Organisms may develop resistance to the toxic effects of high concentrations of radionuclides. Adaptation occurs at the population level by genetic selection for more resistant organisms. To date, very few field studies exist in which adaptation has been addressed, despite the fact that it represents an unknown influence on observed biological responses.

Evolutionary and functional analyses of cytochrome P4501A promoter polymorphisms in natural populations

The functional importance of variable, transcriptional regulatory sequences within and among natural populations is largely unexplored. We analysed the cytochrome P4501A (CYP1A) promoter in three populations of the minnow, Fundulus heteroclitus, because two SNPs in the promoter and first intron of CYP1A are under selection in populations adapted to pollutants. To define the importance of these SNPs, 1630 bp of the CYP1A promoter and first intron and exon were sequenced in eight individuals from three populations: a population from a polluted environment resistant to some aromatic pollutants and two flanking reference populations. CYP1A is induced by many aromatic pollutants, but in populations adapted to pollutants, CYP1A has been shown to be refractory to induction. We were interested in understanding whether variation in the CYP1A promoter explains mechanism(s) of adaptation to these aromatic pollutants. The CYP1A promoter was extremely variable (an average of 9.3% of the promoter nucleotides varied among all populations) and exhibited no fixed differences between populations. As CYP1A is poorly inducible in adapted fish, we hypothesized that CYP1A promoter regions might vary functionally between populations. Unexpectedly, in vitro analysis showed significantly greater transcription from CYP1A promoters found in the population from the polluted environment relative to promoters found in both reference populations. Thus, despite extensive variation among populations and lack of fixed differences between populations, individuals from a polluted environment have significantly enhanced promoter activity. These data demonstrate that intraspecific variation, which provides the raw material for natural selection to act on, can occur while maintaining promoter function.

Species-dependent sensitivity to contaminants: an approach using primary hepatocyte cultures with three marine fish species

There is limited knowledge about the sensitivity of different fish species to environmental pollutants. Such information is pivotal in risk assessment and to understand why some species appear to be more tolerant to contaminants than others. The aim of the current study was to evaluate whether primary hepatocyte cultures of three marine fish species could be established in the field and whether their sensitivity to selected contaminants would differ. Primary hepatocyte cultures of three marine fish species (plaice, long rough dab, Atlantic cod) were established and exposed for 24 h to copper (20-2500 mg L⁻¹) and statins (1-200 mg L⁻¹). Endpoints were esterase activity, metabolic activity and reduced glutathione (GSH) content, all using fluorescent probes. Flatfish hepatocytes were more susceptible to copper and statin exposure than hepatocytes from cod. This study has shown that species-dependent differences in contaminant sensitivity can be investigated using primary hepatocyte cultures.

The effects of intramuscular and intraperitoneal injections of benzo[a]pyrene on selected biomarkers in Clarias gariepinus

This study investigated the dose-dependent and time-course effects of intramuscular (i.m.) and intraperitoneal (i.p.) injection of benzo[a]pyrene (BaP) on the biomarkers EROD activity, GST activity, concentrations of BaP metabolites in bile, and visceral fat deposits (Lipid Somatic Index, LSI) in African catfish (Clarias gariepinus). Intraperitoneal injection resulted in 4.5 times higher accumulation of total selected biliary FACs than i.m. injection. Hepatic GST activities were inhibited by BaP via both injection methods. Dose-response relationships between BaP injection and both biliary FAC concentrations and hepatic GST activities were linear in the i.p. injected group but nonlinear in the i.m. injected fish. Hepatic EROD activity and LSI were not significantly affected by BaP exposure by either injection route. We conclude that i.p. is a more effective route of exposure than i.m. for future ecotoxicological studies of PAH exposure in C. gariepinus.

Strong links between metal contamination, habitat modification and estuarine larval fish distributions

Changes to larval fish assemblages may have far reaching ecological impacts. Correlations between habitat modification, contamination and marine larval fish communities have rarely been assessed in situ. We investigated links between the large-scale distribution of stressors and larval fish assemblages in estuarine environments. Larval fish communities were sampled using a benthic sled within the inner and outer zones of three heavily modified and three relatively unmodified estuaries. Larval abundances were significantly greater in modified estuaries, and there were trends towards greater diversity in these systems. Differences in larval community composition were strongly related to sediment metal levels and reduced seagrass cover. The differences observed were driven by two abundant species, Paedogobius kimurai and Ambassis jacksoniensis, which occurred in large numbers almost exclusively in highly contaminated and pristine locations respectively. These findings suggest that contamination and habitat alteration manifest in substantial differences in the composition of estuarine larval fish assemblages.

Transcriptomic assessment of resistance to effects of an aryl hydrocarbon receptor (AHR) agonist in embryos of Atlantic killifish (Fundulus heteroclitus) from a marine Superfund site

Background

Populations of Atlantic killifish (Fundulus heteroclitus) have evolved resistance to the embryotoxic effects of polychlorinated biphenyls (PCBs) and other halogenated and nonhalogenated aromatic hydrocarbons that act through an aryl hydrocarbon receptor (AHR)-dependent signaling pathway. The resistance is accompanied by reduced sensitivity to induction of cytochrome P450 1A (CYP1A), a widely used biomarker of aromatic hydrocarbon exposure and effect, but whether the reduced sensitivity is specific to CYP1A or reflects a genome-wide reduction in responsiveness to all AHR-mediated changes in gene expression is unknown. We compared gene expression profiles and the response to 3,3',4,4',5-pentachlorobiphenyl (PCB-126) exposure in embryos (5 and 10 dpf) and larvae (15 dpf) from F. heteroclitus populations inhabiting the New Bedford Harbor, Massachusetts (NBH) Superfund site (PCB-resistant) and a reference site, Scorton Creek, Massachusetts (SC; PCB-sensitive).

Results

Analysis using a 7,000-gene cDNA array revealed striking differences in responsiveness to PCB-126 between the populations; the differences occur at all three stages examined. There was a sizeable set of PCB-responsive genes in the sensitive SC population, a much smaller set of PCB-responsive genes in NBH fish, and few similarities in PCB-responsive genes between the two populations. Most of the array results were confirmed, and additional PCB-regulated genes identified, by RNA-Seq (deep pyrosequencing).

Conclusions

The results suggest that NBH fish possess a gene regulatory defect that is not specific to one target gene such as CYP1A but rather lies in a regulatory pathway that controls the transcriptional response of multiple genes to PCB exposure. The results are consistent with genome-wide disruption of AHR-dependent signaling in NBH fish.

Directional genetic selection by pulp mill effluent on multiple natural populations of three-spined stickleback (Gasterosteus aculeatus)

Contamination can cause a rapid environmental change which may require populations to respond with evolutionary changes. To evaluate the effects of pulp mill effluents on population genetics, we sampled three-spined sticklebacks (Gasterosteus aculeatus) near four pulp mills and four adjacent reference sites and analyzed Amplified Fragment Length Polymorphism (AFLP) to compare genetic variability. A fine scale genetic structure was detected and samples from polluted sites separated from reference sites in multidimensional scaling plots (P<0.005, 1000 permutations) and locus-by-locus Analysis of Molecular Variance (AMOVA) further confirmed that habitats are significantly separated (F(ST)=0.021, P<0.01, 1023 permutations). The amount of genetic variation between populations did not differ between habitats, and populations from both habitats had similar levels of heterozygosity (polluted sites Nei's Hs=0.11, reference sites Nei's Hs=0.11). Still, pairwise F(ST): s between three, out of four, pairs of polluted-reference sites were significant. A F(ST)-outlier analysis showed that 21 (8.4%) loci were statistically different from a neutral distribution at the P<0.05 level and therefore indicated to be under divergent selection. When removing 13 F(ST)-outlier loci, significant at the P<0.01 level, differentiation between habitats disappeared in a multidimensional scaling plot. In conclusion, pulp mill effluence has acted as a selective agent on natural populations of G. aculeatus, causing a convergence in genotype composition change at multiple sites in an open environment.

The future of Baltic Sea populations: local extinction or evolutionary rescue?

Environmental change challenges local and global survival of populations and species. In a species-poor environment like the Baltic Sea this is particularly critical as major ecosystem functions may be upheld by single species. A complex interplay between demographic and genetic characteristics of species and populations determines risks of local extinction, chances of re-establishment of lost populations, and tolerance to environmental changes by evolution of new adaptations. Recent studies show that Baltic populations of dominant marine species are locally adapted, have lost genetic variation and are relatively isolated. In addition, some have evolved unusually high degrees of clonality and others are representatives of endemic (unique) evolutionary lineages. We here suggest that a consequence of local adaptation, isolation and genetic endemism is an increased risk of failure in restoring extinct Baltic populations. Additionally, restricted availability of genetic variation owing to lost variation and isolation may negatively impact the potential for evolutionary rescue following environmental change.

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

Atlantic killifish (Fundulus heteroclitus) inhabiting the PCB-contaminated Superfund site in New Bedford Harbor (MA, USA) have evolved genetic resistance to the toxic effects of these compounds. They also lack induction of cytochrome P4501A (CYP1A) and other aryl hydrocarbon receptor (AHR)-dependent responses after exposure to AHR agonists, suggesting an overall down-regulation of the AHR signaling pathway. In this study, we hypothesized that the genetic resistance is due to altered AHR expression resulting from hypermethylation of DNA in the promoter region of AHR genes in fish inhabiting New Bedford Harbor. To test this hypothesis, we cloned and sequenced AHR1 and AHR2 promoter regions and employed bisulfite conversion-polymerase chain reaction (BS-PCR) followed by clonal analysis to compare the methylation status of CpG islands of AHR1 and AHR2 in livers of adult killifish collected from New Bedford Harbor and a reference site (Scorton Creek, MA). No significant differences in methylation profiles were observed in either AHR1 or AHR2 promoter regions between NBH and SC fish. However, hypermethylation of the AHR1 promoter correlated with low expression of transcripts in the liver in both populations. In comparison to AHR1, hepatic mRNA expression of AHR2 is high and its promoter is hypomethylated. Taken together, our results suggest that genetic resistance to contaminants in NBH fish is not due to altered methylation of AHR promoter regions, but that promoter methylation may control tissue-specific expression of AHR genes in killifish.

Ecologically and evolutionarily important SNPs identified in natural populations

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.

Predicting chemical impacts on vertebrate endocrine systems

Animals have evolved diverse protective mechanisms for responding to toxic chemicals of both natural and anthropogenic origin. From a governmental regulatory perspective, these protective responses complicate efforts to establish acceptable levels of chemical exposure. To explore this issue, we considered vertebrate endocrine systems as potential targets for environmental contaminants. Using the hypothalamic-pituitary-thyroid (HPT), hypothalamic-pituitary-gonad (HPG), and hypothalamic-pituitary-adrenal (HPA) axes as case examples, we identified features of these systems that allow them to accommodate and recover from chemical insults. In doing so, a distinction was made between effects on adults and those on developing organisms. This distinction was required because endocrine system disruption in early life stages may alter development of organs and organ systems, resulting in permanent changes in phenotypic expression later in life. Risk assessments of chemicals that impact highly regulated systems must consider the dynamics of these systems in relation to complex environmental exposures. A largely unanswered question is whether successful accommodation to a toxic insult exerts a fitness cost on individual animals, resulting in adverse consequences for populations. Mechanistically based mathematical models of endocrine systems provide a means for better understanding accommodation and recovery. In the short term, these models can be used to design experiments and interpret study findings. Over the long term, a set of validated models could be used to extrapolate limited in vitro and in vivo testing data to a broader range of untested chemicals, species, and exposure scenarios. With appropriate modification, Tier 2 assays developed in support of the U.S. Environmental Protection Agency's Endocrine Disruptor Screening Program could be used to assess the potential for accommodation and recovery and inform the development of mechanistically based models.

Adverse outcome pathways and ecological risk assessment: bridging to population-level effects

Maintaining the viability of populations of plants and animals is a key focus for environmental regulation. Population-level responses integrate the cumulative effects of chemical stressors on individuals as those individuals interact with and are affected by their conspecifics, competitors, predators, prey, habitat, and other biotic and abiotic factors. Models of population-level effects of contaminants can integrate information from lower levels of biological organization and feed that information into higher-level community and ecosystem models. As individual-level endpoints are used to predict population responses, this requires that biological responses at lower levels of organization be translated into a form that is usable by the population modeler. In the current study, we describe how mechanistic data, as captured in adverse outcome pathways (AOPs), can be translated into modeling focused on population-level risk assessments. First, we describe the regulatory context surrounding population modeling, risk assessment and the emerging role of AOPs. Then we present a succinct overview of different approaches to population modeling and discuss the types of data needed for these models. We describe how different key biological processes measured at the level of the individual serve as the linkage, or bridge, between AOPs and predictions of population status, including consideration of community-level interactions and genetic adaptation. Several case examples illustrate the potential for use of AOPs in population modeling and predictive ecotoxicology. Finally, we make recommendations for focusing toxicity studies to produce the quantitative data needed to define AOPs and to facilitate their incorporation into population modeling.

Assessing contaminated sediments in the context of multiple stressors

Sediments have a major role in ecosystem functioning but can also act as physical or chemical stressors. Anthropogenic activities may change the chemical constituency of sediments and the rate, frequency, and extent of sediment transport, deposition, and resuspension. The importance of sediments as stressors will depend on site ecosystem attributes and the magnitude and preponderance of co-occurring stressors. Contaminants are usually of greater ecological consequence in human-modified, depositional environments, where other anthropogenic stressors often co-occur. Risk assessments and restoration strategies should better consider the role of chemical contamination in the context of multiple stressors. There have been numerous advances in the temporal and spatial characterization of stressor exposures and quantification of biological responses. Contaminated sediments causing biological impairment tend to be patchy, whereas more pervasive anthropogenic stressors, such as alterations to habitat and flow, physical disturbance, and nutrient addition, may drive large-scale ecosystem responses. A systematic assessment of relevant ecosystem attributes and reference conditions can assist in understanding the importance of sediments in the context of other stressors. Experimental manipulations then allow for the controlled study of dominant stressors and the establishment of causal links. This approach will result in more effective management of watersheds and waterways.

Characterization of the recalcitrant CYP1 phenotype found in Atlantic killifish (Fundulus heteroclitus) inhabiting a Superfund site on the Elizabeth River, VA

Fundulus heteroclitus (Atlantic killifish) found at the Atlantic Wood Industries Superfund site on the Elizabeth River (ER) in Portsmouth, VA (USA), have been shown to be resistant to the teratogenic effects of creosote-contaminated sediments found at this highly contaminated site. Many of the polycyclic aromatic hydrocarbons (PAHs) found at the ER are known to activate the aryl hydrocarbon receptor (AHR), and are thought to mediate their toxic effects through this pathway. Activation of the AHR results in the induction of several Phase I and II metabolic enzymes. It has been previously shown that the AHR of killifish from the ER are refractory to induction by AHR agonists. To more fully characterize this altered AHR response, we exposed embryos from the ER and from a reference site on King's Creek, VA (KC) to two PAHs, benzo[alpha]pyrene (BaP) and benzo[k]fluoranthene (BkF), and to the dioxin-like compound (DLC), 3,3',4,4',5-pentachlorobiphenyl (PCB126). We compared their developmental and molecular responses by screening the embryos for CYP1A enzyme activity, cardiac deformities, and mRNA expression of CYP1A, CYP1B1, CYP1C1, and AHR2. Basal gene expression of both CYP1A and CYP1B1 was 40% higher in the KC control embryos compared to those from the ER, while AHR2 and CYP1C1 were not significantly different between the populations. Exposure of KC embryos to BaP, BkF, and PCB126 induced CYP1A activity and cardiac deformities. In contrast, CYP1A activity was induced in ER embryos only in response to BkF exposure, although this induction in ER embryos was significantly lower than that observed in KC fish at comparable concentrations. ER embryos did not develop cardiac deformities in response to any of the chemicals tested. CYP1A, CYP1B1 and CYP1C1 mRNA were all significantly induced in the KC embryos after exposure to BaP, BkF and PCB126. Exposure to BaP and BkF in ER embryos resulted in a significant induction of CYP1A mRNA, albeit significantly lower than observed in KC fish. Interestingly, BaP exposure resulted in induction of CYP1B1 at comparable levels in embryos from both populations. CYP1s were not induced in ER embryos in response to PCB126, nor was CYP1C1 for any treatment examined. Additionally, AHR2 was not significantly induced for any of the treatment groups. This study further characterizes the AHR response in killifish, and provides greater insight into the adapted ER phenotype. The ER adaptation involves the suppression of normal AHR-inducible gene expression for all three CYP1 genes, and therefore is likely an alteration in AHR signaling or control.

Embryonic gene expression among pollutant resistant and sensitive Fundulus heteroclitus populations

Changes in gene expression, coupled with biochemical, physiological, and behavioral alterations, play a critical role in adaptation to environmental stress. Our goal was to explore ways natural populations may have adapted to local, polluted environments. We took advantage of natural populations of Fundulus heteroclitus, one of the few studied fish species in North America that has established resistant populations in highly contaminated urban estuaries. We analyzed morphology, physiology, and gene expression of developing F. heteroclitus embryos during late organogenesis (stage 31); these embryos were from both resistant and sensitive populations and were raised in a common, unpolluted environment. While cardiac heart rates show significant differences between embryos of parents from clean and heavily contaminated Superfund sites, time-to-stage, embryo morphology, and gene expression profile analyses do not differ significantly between untreated embryos from resistant and sensitive populations. Further evaluation that includes tissue-specific approaches in gene expression analysis and larger sample sizes may be necessary to highlight important phenotypes associated with mechanisms of sensitivity and resistance among natural F. heteroclitus embryo populations. Alternatively, population differences may be masked by developmental canalization, and biologically important differences between sensitive and resistant embryos may only manifest with exposure (e.g., be dependent on gene by environment interactions).

Comparative study about the effects of pollution on glass and yellow eels (Anguilla anguilla) from the estuaries of Minho, Lima and Douro Rivers (NW Portugal)

The health status of eels (Anguilla anguilla) developing in three estuaries of the NW Portuguese coast with different types and levels of pollution was compared in relation to morphometric parameters, Fulton condition index (F index) and several biomarkers. Relatively to the reference population, glass eels from the Lima estuary had lower weight and length, cholinesterase (ChE) and lactate dehydrogenase (LDH) inhibition, and lower levels of some anti-oxidant parameters, while yellow eels also showed a decreased F index, and increased Na(+)/K(+)-ATPase and lipid peroxidation (LPO) levels. Relatively to the reference population, glass eels from the Douro estuary had increased Na(+)/K(+)-ATPase and glutathione-S-transferase activities and LDH inhibition, while yellow eels also had ChE inhibition and increased LPO. Overall, these results indicate that eels from polluted estuaries showed a poor health status than those from a reference estuary, and that adverse effects become more pronounced after spending several years in polluted estuaries.

SNP identification, verification, and utility for population genetics in a non-model genus

BACKGROUND

By targeting SNPs contained in both coding and non-coding areas of the genome, we are able to identify genetic differences and characterize genome-wide patterns of variation among individuals, populations and species. We investigated the utility of 454 sequencing and MassARRAY genotyping for population genetics in natural populations of the teleost, Fundulus heteroclitus as well as closely related Fundulus species (F. grandis, F. majalis and F. similis).

RESULTS

We used 454 pyrosequencing and MassARRAY genotyping technology to identify and type 458 genome-wide SNPs and determine genetic differentiation within and between populations and species of Fundulus. Specifically, pyrosequencing identified 96 putative SNPs across coding and non-coding regions of the F. heteroclitus genome: 88.8% were verified as true SNPs with MassARRAY. Additionally, putative SNPs identified in F. heteroclitus EST sequences were verified in most (86.5%) F. heteroclitus individuals; fewer were genotyped in F. grandis (74.4%), F. majalis (72.9%), and F. similis (60.7%) individuals. SNPs were polymorphic and showed latitudinal clinal variation separating northern and southern populations and established isolation by distance in F. heteroclitus populations. In F. grandis, SNPs were less polymorphic but still established isolation by distance. Markers differentiated species and populations.

CONCLUSIONS

In total, these approaches were used to quickly determine differences within the Fundulus genome and provide markers for population genetic studies.

Polluted-site killifish (Fundulus heteroclitus) embryos are resistant to organic pollutant-mediated induction of CYP1A activity, reactive oxygen species, and heart deformities

Exposure to coplanar polychlorinated biphenyls (PCBs) and polycyclic aromatic hydrocarbons (PAHs) can adversely affect fish embryonic development, induce expression of cytochrome P4501A (CYP1A), and increase reactive oxygen species (ROS) production, effects believed to be mediated by the aryl hydrocarbon receptor (AHR). Killifish (Fundulus heteroclitus) populations in New Bedford Harbor, Massachusetts, USA (NBH) and Newark Bay, New Jersey, USA (NB) are generationally exposed to coplanar PCBs and PAHs and have developed resistance to PCB mediated induction of CYP1A. We hypothesized that fish resistant to CYP1A induction would also exhibit resistance to PCB and PAH induced ROS production and teratogenesis. Killifish embryos from two contaminated (NB, NBH) and two reference-site populations were exposed to vehicle or 3,3'4,4'5-pentachlorobiphenyl (PCB126) or 3-methylcholanthrene (3-MC) and evaluated for in ovo CYP1A activity, heart deformities, and ROS production. Both chemicals significantly increased in ovo ethoxyresorufin-O-deethylase (EROD) and ROS production in reference-site embryos. These chemicals provoked only moderate induction of in ovo EROD in NBH and NB embryos, and neither PCB126 nor 3-MC induced ROS production in these populations. Similarly, heart deformities were significantly induced by PCB126 in reference-site embryos, but had no significant effects on NB and NBH animals. These results indicate that fish resistant to CYP1A induction also exhibit decreased sensitivity to PCB126 and 3-MC-induced ROS production and teratogenesis. These findings further our understanding of toxicant resistance by demonstrating that reduced response to coplanar PCBs and PAHs extends beyond resistance to CYP1A induction to resistance to the physiological and teratogenic effects of these toxicants, responses that undoubtedly contribute to the increased survival of killifish inhabiting contaminated sites.

Identification of mitochondrial cytochrome P450 induced in response to polycyclic aromatic hydrocarbons in the mummichog (Fundulus heteroclitus)

Increasing evidence suggests that polycyclic aromatic hydrocarbons (PAHs) such as benzo[a]pyrene (BaP) are localized to the mitochondria. Because the toxic effects of many PAHs are the result of metabolism by cytochrome P4501A (CYP1A), it is important to investigate whether active forms of these enzymes can be identified in the mitochondria. In this study, we identified mitochondrial P450s with a monoclonal antibody against scup (Stenotomus chrysops) CYP1A in the isolated mitochondrial fraction of the liver from adult male mummichog (Fundulus heteroclitus) livers. The size of the protein in the mitochondria was similar to that of microsomal CYP1A. Fish dosed with 10mg/kg BaP had increased EROD activity in the mitochondrial fraction compared to controls. In mummichog larvae dosed with 100 microg/L BaP and 100 microg/L benzo[k]fluoranthene, CYP1A protein levels as well as enzyme activity were elevated. However, fish from a PAH-polluted Superfund site (Elizabeth River, Portsmouth VA) showed recalcitrant mitochondrial CYP1A protein levels and enzyme activity in a similar manner to microsomal CYP1A.

Lack of constitutive and inducible ethoxyresorufin-O-deethylase activity in the liver of suckermouth armored catfish (Hypostomus affinis and Hypostomus auroguttatus, Loricariidae)

We investigated the presence and inducibility of CYP1A in suckermouth catfish (Hypostomus affinis and Hypostomus auroguttatus, Loricariidae), tilapia (Oreochromis niloticus, Cichlidae) and mice (Mus musculus, Muridae). Alkoxyresorufin-O-dealkylases (EROD, MROD, PROD and BROD) were detected and proved to be inducible (beta-naphthoflavone, BNF or dimethylbenz[a]anthracene, DMBA, 50 mg/kg bw ip) in liver microsomes from tilapia and mice. In loricariids, alkoxyresorufin-O-dealkylases were either undetectable (MROD/EROD) or very low (PROD/BROD), and so they remained after treatment with BNF or DMBA. Ethoxycoumarin-O-deethylase (ECOD) was recorded in all species and proved not to be inducible by BNF or DMBA. In loricariids and tilapia, ECOD was not depressed by a concentration of alpha-naphthoflavone (CYP1A-inhibitor) that markedly depressed EROD in tilapia. A CYP1A-like protein was detected by a monoclonal antibody in rats, mice and tilapia, but not in loricariids. A polyclonal antibody, however, detected a CYP1A-like protein in liver microsomes of loricariids. Suckermouth catfish, rats, mice and tilapia express a protein reactive with a polyclonal antibody against trout CYP3A. Loricariids and tilapia exhibited marked genotoxic responses (enhanced incidence of micronucleated erythrocytes) following treatment DMBA (50 mg/kg bw ip), a promutagen activated by CYP1A/1B. Therefore, although not exhibiting EROD, a CYP1A-mediated activity, loricariids converted DMBA into its genotoxic metabolites. Our findings suggest that the CYP1A-like protein of locariid catfish recognizes DMBA, but not ethoxyresorufin, as a substrate.

The effects of heavy metal mine drainage on population size structure, reproduction, and condition of western mosquitofish, Gambusia affinis

Anthropogenic degradation of aquatic environments worldwide results in disturbed habitats, altered communities, and stressed populations. Surface waters located in an abandoned lead-zinc mining district in northeastern Oklahoma are no exception. This study examines the reproductive and somatic responses of a pollution-tolerant fish, the western mosquitofish, Gambusia affinis (Teleostei: Poeciliidae), living in mine outflow waters contaminated by heavy metals. Populations were sampled from four streams, which were classified into three habitat types. Populations from Tar Creek and an Unnamed Tributary of Tar Creek receive direct input of mine drainage, while populations living in reference creeks are not known to have mining influence. The influence of mine drainage directly or indirectly (via altered competitor and predator regimes or changes in food availability) affects G. affinis at both the population and the individual level. Metal-contaminated sites had reduced proportions of males and reproductively active females and altered male population size structures. Individual-level effects were apparent, as all G. affinis from Tar Creek invested less in liver weights, and mature males and reproductively active females from Tar Creek invested less in gonad weights. Furthermore, males from impacted sites were significantly lighter than those from reference creeks. Gravid females from Tar Creek had smaller clutch sizes, but average embryo weight did not differ among streams.

Genotoxicity of environmentally relevant concentrations of water-soluble oil components in cod (Gadus morhua)

Large discharges from oil and gas production platforms (produced water) have led to concerns for adverse biological effects in marine areas. The aim of this study was to investigate the development of DNA adductformation and related biomarkers in fish after chronic exposure to water-soluble components of oil. Atlantic cod (Gadus morhua) were exposed for up to 44 weeks to three treatments (low, pulsed, high) containing environmentally relevant concentrations of low-molecular-weight polycyclic aromatic hydrocarbons (PAHs) and short-chained alkylphenols (APs). A time- and dose-related pattern of DNA adduct formation (measured using 32P-postlabeling) was observed. The results suggested that an extended exposure period (more than 16 weeks) would be required for the formation of DNA adduct levels above background. Interestingly, fish receiving pulsed high exposure did not develop elevated concentrations of DNA adducts, possibly due to DNA repair processes. No obvious relationship between DNA adduct concentration and cytochrome P4501A activity (EROD) was observed. This study has demonstrated the genotoxic potential of water-soluble oil components, relevant for operational discharges (produced water) and chronic oil spills. The quantification of PAH metabolites in bile and hepatic DNA adduct formation appear to be suitable for environmental monitoring of chronic oil pollution.

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.

Changes in gene expression due to chronic exposure to environmental pollutants

Populations of the teleost fish Fundulus heteroclitus inhabit and have adapted to highly polluted Superfund sites that are contaminated with persistent toxic chemicals. Populations inhabiting different Superfund sites provide independent contrasts for studying mechanisms of toxicity and resistance due to exposure to environmental pollutants. To identify both shared and unique responses to chronic pollutant exposure, liver, metabolic gene expression in F. heteroclitus populations from each of three Superfund sites (New Bedford Harbor, MA; Newark Bay, NJ; and Elizabeth River, VA) were compared to two flanking reference site populations (nine populations in total). In comparisons to their two clean reference sites, the three Superfund sites had 8-32% of genes with altered expression patterns. Between any two Superfund populations, up to nine genes (4%) show a conserved response, yet among all three populations, there was no gene which had a conserved, altered pattern of expression. Across all three Superfund sites in comparison to all six-reference populations, the most significant gene was fatty acid synthase. Fatty acid synthase is involved in the storage of excess energy as fat, and its lesser expression in the polluted populations suggests that the polluted populations may have limited energy stores. In contrast to previous studies of metabolic gene expression in F. heteroclitus, body weight was a significant covariate for many of the genes which could reflect accumulation and different body burdens of pollutants. Overall, the altered gene expression in these populations likely represents both induced and adaptive changes in gene expression.

Signatures of selection in natural populations adapted to chronic pollution

Background

Populations of the teleost fish Fundulus heteroclitus appear to flourish in heavily polluted and geographically separated Superfund sites. Populations from three Superfund sites (New Bedford Harbor, MA, Newark Bay, NJ, and Elizabeth River, VA) have independently evolved adaptive resistance to chemical pollutants. In these polluted populations, natural selection likely has altered allele frequencies of loci that affect fitness or that are linked to these loci. The aim of this study was to identify loci that exhibit non-neutral behavior in the F. heteroclitus genome in polluted populations versus clean reference populations.

Results

To detect signatures of natural selection and thus identify genetic bases for adaptation to anthropogenic stressors, we examined allele frequencies for many hundreds of amplified fragment length polymorphism markers among populations of F. heteroclitus. Specifically, we contrasted populations from three Superfund sites (New Bedford Harbor, MA, Newark Bay, NJ, and Elizabeth River, VA) to clean reference populations flanking the polluted sites. When empirical F_ST values were compared to a simulated distribution of F_ST values, 24 distinct outlier loci were identified among pairwise comparisons of pollutant impacted F. heteroclitus populations and both surrounding reference populations. Upon removal of all outlier loci, there was a strong correlation (R² = 0.79, p < 0.0001) between genetic and geographical distance. This apparently neutral evolutionary pattern was not evident when outlier loci were included (R² = 0.092, p = 0.0721). Two outlier loci were shared between New Bedford Harbor and Elizabeth River populations, and two different loci were shared between Newark Bay and Elizabeth River populations.

Conclusion

In total, 1% to 6% of loci are implicated as being under selection or linked to areas of the genome under selection in three F. heteroclitus populations that reside in polluted estuaries. Shared loci among polluted sites indicate that selection may be acting on multiple loci involved in adaptation, and loci shared between polluted sites potentially are involved in a generalized adaptive response.

Nonadditive effects of PAHs on Early Vertebrate Development: mechanisms and implications for risk assessment

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants. Traditionally, much of the research has focused on the carcinogenic potential of specific PAHs, such as benzo(a)pyrene, but recent studies using sensitive fish models have shown that exposure to PAHs alters normal fish development. Some PAHs can induce a teratogenic phenotype similar to that caused by planar halogenated aromatic hydrocarbons, such as dioxin. Consequently, mechanism of action is often equated between the two classes of compounds. Unlike dioxins, however, the developmental toxicity of PAH mixtures is not necessarily additive. This is likely related to their multiple mechanisms of toxicity and their rapid biotransformation by CYP1 enzymes to metabolites with a wide array of structures and potential toxicities. This has important implications for risk assessment and management as the current approach for complex mixtures of PAHs usually assumes concentration addition. In this review we discuss our current knowledge of teratogenicity caused by single PAH compounds and by mixtures and the importance of these latest findings for adequately assessing risk of PAHs to humans and wildlife. Throughout, we place particular emphasis on research on the early life stages of fish, which has proven to be a sensitive and rapid developmental model to elucidate effects of hydrocarbon mixtures.

Oxidative damage in gill of Mytilus edulis from Merseyside, UK, and reversibility after depuration

Mussels were collected from the urban/industrialized site of New Brighton, Merseyside and the relatively non-industrial site of Llandudno, North Wales. All mussels were identified as Mytilus edulis by PCR amplification of Mefp1. DNA single strand breaks and 8-oxo-7,8-dihydro-2'-deoxyguanosine were measured in gill within 24h of collection, using the COMET assay, both with and without formamidopyrimidine glycosylase. Gill lipid peroxidation was also measured within 24h. No difference between sites was found for frank SSB and malonaldehyde levels, however 8-oxo-dG and 4-hydroxynonenal were significantly greater in New Brighton mussels compared to Llandudno mussels. After 1-month laboratory maintenance, lipid peroxidation and 8-oxo-dG levels were lower. In contrast, frank SSB were higher. This could reflect enhanced DNA repair excision, though we cannot exclude the possibility of other non-oxidative DNA damage. The results suggest that laboratory maintenance allows recovery from environmentally induced oxidative damage, which was more extensive at Merseyside.

Convergence and divergence in gene expression among natural populations exposed to pollution

Background

Natural populations of the teleost fish Fundulus heteroclitus tolerate a broad range of environmental conditions including temperature, salinity, hypoxia and chemical pollutants. Strikingly, populations of Fundulus inhabit and have adapted to highly polluted Superfund sites that are contaminated with persistent toxic chemicals. These natural populations provide a foundation to discover critical gene pathways that have evolved in a complex natural environment in response to environmental stressors.

Results

We used Fundulus cDNA arrays to compare metabolic gene expression patterns in the brains of individuals among nine populations: three independent, polluted Superfund populations and two genetically similar, reference populations for each Superfund population. We found that up to 17% of metabolic genes have evolved adaptive changes in gene expression in these Superfund populations. Among these genes, two (1.2%) show a conserved response among three polluted populations, suggesting common, independently evolved mechanisms for adaptation to environmental pollution in these natural populations.

Conclusion

Significant differences among individuals between polluted and reference populations, statistical analyses indicating shared adaptive changes among the Superfund populations, and lack of reduction in gene expression variation suggest that common mechanisms of adaptive resistance to anthropogenic pollutants have evolved independently in multiple Fundulus populations. Among three independent, Superfund populations, two genes have a common response indicating that high selective pressures may favor specific responses.

Antioxidative stress proteins and their gene expression in brown trout (Salmo trutta) from three rivers with different heavy metal levels

Three populations of brown trout (Salmo trutta) exposed to different metal levels in their natural environments, were studied with respect to antioxidants metallothionein (MT), superoxide dismutase (SOD) and catalase (CAT) as well as for corresponding mRNA levels. In addition, mRNA levels were studied for glutathione peroxidase (GPx) and glutathione reductase (GR). The Cd/Zn-exposed trout (Naustebekken River) had higher accumulated levels of Cd, Cu and Zn in gills, and higher levels of MT (both protein and mRNA) in liver and kidney as well as in gills compared to the Cu-exposed trout (Rugla River) and trout from an uncontaminated reference river (Stribekken River). Less MT found in the Cu-exposed trout may increase susceptibility to oxidative stress, but no higher levels of antioxidant mRNAs were found in gills of these trouts. The data indicated that chronic exposures of brown trout to Cd, Zn and/or Cu did not involve maintenance of high activities of SOD and CAT enzymes in gills, although SOD mRNA levels were higher in the Cd/Zn-exposed trout. In livers, mRNA levels of SOD, CAT and GPx were higher in the metal-exposed trout, but in the case of GR this was only seen in kidneys of Cd/Zn-exposed trout. However, both metal-exposed groups had higher activities of SOD enzyme in liver compared to the unexposed reference trout, and CAT activity was found to be higher in kidneys of Cu-exposed trout. The Cu-exposed trout did not seem to rely on MT production to avoid Cu toxicity in gills, but rather by keeping the Cu uptake at a low level. A coordinated expression of different stress genes may also be important in chronic metal exposure. It may be concluded that the observed metal effects relies on acclimation rather than on genetic adaptation in the metal exposed populations.

Comparison of hepatic and extra hepatic induction of cytochrome P4501A by graded doses of aryl hydrocarbon receptor agonists in Atlantic tomcod from two populations

Atlantic tomcod Microgadus tomcod from the Hudson River, New York, are exposed to high levels of polycyclic aromatic hydrocarbons (PAHs) and bioaccumulate mixtures of polychlorinated biphenyls (PCBs), polychlorinated dibenzo-p-dioxins and polychlorinatedfurans (PCDD/Fs). Previous studies demonstrated that hepatic cytochrome P4501A (CYP1A) mRNA was not inducible in tomcod from the Hudson River treated with single doses of PCB77 or 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), but was inducible with PAHs. In this study, we sought to determine if CYP1A mRNA was inducible with higher doses of these and other halogenated aromatic hydrocarbons (HAHs) in Hudson River tomcod and if decreased sensitivity to gene inducibility occurs across all tissues. Tomcod from the Hudson River and the cleaner Miramichi River, New Brunswick, were treated individually with graded doses of TCDD and coplanar PCBs (PCB77, PCB81, PCB126, PCB169) and profiles of hepatic CYP1A mRNA expression were compared between the two populations. CYP1A mRNA inducibility was also compared in multiple tissues of tomcod from the two rivers that were treated with PCB77. Additionally, hepatic CYP1A mRNA was characterized in Miramichi River tomcod treated with pairs of PCB congeners that included aryl hydrocarbon receptor (AHR) agonists and antagonists. Hepatic CYP1A mRNA was significantly inducible by all agonists in tomcod from the Miramichi River and TCDD and two of four PCBs in tomcod from the Hudson River. CYP1A mRNA was also significantly inducible in four of five tissues of tomcod from the Miramichi River but only in liver of Hudson River tomcod. In summary, CYP1A mRNA inducibility was approximately two orders of magnitude less sensitive in tomcod from the Hudson River than in those from the Miramichi River. But when achieved, maximum levels of CYP1A expression were similar in tomcod from the two populations. Co-administration of PCB126 and PCB77 did not produce significantly greater CYP1A mRNA induction than administration of PCB126 alone and co-administration of mono-ortho-substituted PCB105 significantly decreased CYP1A mRNA inducibility by PCB77. These results indicate that CYP1A mRNA expression is significantly inducible by HAHs in tomcod from the Hudson River and suggest that all components of the AHR pathway are present and functional, but that the pathway is less sensitive to activation than in tomcod from the Miramichi River. Our results also indicate that CYP1A mRNA levels in environmentally exposed fish may not reflect additive tissue burdens of PCB congeners.

Evidence of spatially extensive resistance to PCBs in an anadromous fish of the Hudson River

Populations of organisms that are chronically exposed to high levels of chemical contaminants may not suffer the same sublethal or lethal effects as naive populations, a phenomenon called resistance. Atlantic tomcod (Microgadus tomcod) from the Hudson River, New York, are exposed to high concentrations of polycyclic aromatic hydrocarbons (PAHs) and bioaccumulate polychlorinated biphenyls (PCBs), polychlorinated dibenzo-p-dioxins (PCDDs), and polychlorinated dibenzofurans (PCDFs). They have developed resistance to PCBs and PCDDs but not to PAHs. Resistance is largely heritable and manifests at early-life-stage toxic end points and in inducibility of cytochrome P4501A (CYP1A) mRNA expression. Because CYP1A induction is activated by the aryl hydrocarbon receptor (AHR) pathway, as are most toxic responses to these compounds, we sought to determine the geographic extent of resistance to CYP1A mRNA induction by PCBs in the Hudson River tomcod population. Samples of young-of-the-year tomcod were collected from seven locales in the Hudson River, extending from the Battery at river mile 1 (RM 1) to RM 90, and from the Miramichi River, New Brunswick, Canada. Laboratory-reared offspring of tomcod adults from Newark Bay, in the western portion of the Hudson River estuary, were also used in this study. Fish were partially depurated in clean water and intraperitoneally injected with 10 ppm coplanar PCB-77, 10 ppm benzo[a]pyrene (BaP), or corn oil vehicle, and levels of CYP1A mRNA were determined. CYP1A was significantly inducible by treatment with BaP in tomcod from the Miramichi River, from laboratory-spawned offspring of Newark Bay origin, and from all Hudson River sites spanning 90 miles of river. In contrast, only tomcod from the Miramichi River displayed significantly induced CYP1A mRNA expression when treated with PCB-77. Our results suggest that the population of tomcod from throughout the Hudson River estuary has developed resistance to CYP1A inducibility and probably other toxicities mediated by the AHR pathway. Tomcod from the Hudson River may represent the most geographically expansive population of vertebrates with resistance to chemical pollutants that has been characterized.

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

The aryl hydrocarbon receptor (AHR) and AHR repressor (AHRR) proteins regulate gene expression in response to some halogenated aromatic hydrocarbons and polycyclic aromatic hydrocarbons. The Atlantic killifish is a valuable model of the AHR signaling pathway, but antibodies are not available to fully characterize AHR and AHRR proteins. Using bacterially expressed AHRs, we developed specific and sensitive polyclonal antisera against the killifish AHR1, AHR2, and AHRR. In immunoblots, these antibodies recognized full-length killifish AHR and AHRR proteins synthesized in rabbit reticulocyte lysate, proteins expressed in mammalian cells transfected with killifish AHR and AHRR constructs, and AHR proteins in cytosol preparations from killifish tissues. Killifish AHR1 and AHR2 proteins were detected in brain, gill, kidney, heart, liver, and spleen. Antisera specifically precipitated their respective target proteins in immunoprecipitation experiments with in vitro-expressed proteins. Killifish ARNT2 co-precipitated with AHR1 and AHR2. These sensitive, specific, and versatile antibodies will be valuable to researchers investigating AHR signaling and other physiological processes involving AHR and AHRR proteins.

Biochemical biomarkers in adult female perch (Perca fluviatilis) in a chronically polluted gradient in the Stockholm recipient (Sweden)

A battery of biochemical biomarkers and the SigmaPCB concentration in adult female perch (Perca fluviatilis) verified an aquatic pollution gradient with the city of Stockholm (Sweden) as a point source of anthropogenic substances. The investigation included both an upstream gradient, 46 km westwards through Lake Mälaren, and a downstream gradient, 84 km eastwards through the Stockholm archipelago. Besides the main gradient from Stockholm, there were strong indications of pollution coming from the Baltic Sea. The results indicated a severe pollution situation in central Stockholm, with poor health status of the perch, characterised by increased specific EROD activity in the liver, increased liver EROD somatic index, decreased AChE activity in the muscle, increased amount of DNA adducts in the liver, and a high concentration of biliary 1-pyrenol. In addition, laboratory exposure to common EROD inducers elicited an abnormal response, suggestive of chronic intoxication.

Co-exposure to metals modulates CYP1A mRNA inducibility in Atlantic tomcod Microgadus tomcod from two populations

Populations from urbanized and industrialized sites are often exposed to mixtures of chemical contaminants including aromatic hydrocarbons (AHs) and heavy metals. The effects of mixtures of these contaminants on these populations are largely unknown. The Hudson River Estuary is highly contaminated with a variety of AHs including, PCBs and PAHs, and metals, and its population of Atlantic tomcod Microgadus tomcod bioaccumulates those which are persistent. The Hudson River's tomcod population exhibits resistance to persistent AHs as exemplified by significantly decreased inducibility of hepatic cytochrome P4501A (CYP1A) mRNA. We used hepatic CYP1A mRNA inducibility in tomcod from the Hudson River and a sensitive population to investigate the effects of acute co-exposure to metals on aryl hydrocarbon receptor (AHR)-mediated gene expression. Adult tomcod from the Hudson River and the cleaner Miramichi River were i.p. injected with one dose of benzo[a]pyrene (B[a]P) or coplanar PCB77 and graded doses of four metals, As, Cd, Cr, and Ni, and levels of hepatic CYP1A mRNA and protein were assayed. We observed no effects of metals treatment on basal levels of hepatic CYP1A mRNA expression, but all four metals significantly reduced CYP1A mRNA inducibility in tomcod from one or both populations. The magnitude of the inhibition of CYP1A mRNA inducibility differed among the metals and fish from the two populations. Also, the profile of the metals modulation of induced CYP1A mRNA showed differences that depended on the time after treatment of sacrifice. Our results demonstrate that co-exposure to several metals can impact inducible, but not basal levels of CYP1A expression and perhaps other toxicities mediated by the AHR.

Analysis of CpG methylation in the killifish CYP1A promoter

Fundulus heteroclitus (Atlantic killifish or mummichog) inhabiting a creosote-contaminated Superfund site on the Elizabeth River (VA, USA), exhibit a lack of induction of cytochrome P4501A (CYP1A) mRNA, immunodetectable protein, and catalytic activity after exposure to typical inducers. This "refractory CYP1A phenotype" is not explained by alterations in mRNA expression of known CYP1A transcription factors. Furthermore, the refractory phenotype is lost progressively during development in laboratory-reared F1 generation fish. Thus, while heritable, the refractory CYP1A phenotype does not appear to be genetically based. To test the hypothesis that cytosine methylation at CpG sites in the promoter region of CYP1A underlies the refractory CYP1A phenotype, we employed bisulfite sequencing to compare the methylation status of CpG sites in the CYP1A promoter region of DNA from killifish from the Elizabeth River and a reference site. We examined genomic DNA both from livers of wild-caught adult killifish and from pools of F1 generation embryos raised in the laboratory. In fish from both the contaminated and the reference site, cytosine methylation was not detectable at any of the 34 CpG sites examined, including 3 that are part of putative xenobiotic response elements.

Differential display of hepatic mRNA from killifish (Fundulus heteroclitus) inhabiting a Superfund estuary

Fundulus heteroclitus (Atlantic killifish, mummichog) from a highly contaminated site on the Elizabeth River (VA, USA) are resistant to the toxicity of sediment from the site. However, the mechanistic changes that confer resistance to the toxicity are not yet well understood. We utilized differential display to identify mRNAs that are differentially expressed in hepatic tissue of male and female killifish from the Elizabeth River environment, compared to killifish from a non-contaminated reference site, King's Creek (VA, USA). Seventy-four differentially expressed mRNAs were initially identified (including sex and population-specific differences), and 65 of these were isolated and sequenced. A reverse northern blot array constructed of these cDNAs (plus an additional 15 previously sequenced mRNAs of interest) was used to confirm and quantify expression differences. High interindividual variability was observed in mRNA expression, but multiple differentially expressed mRNAs were identified, including 11 population-specific differences occurring in both sexes, 24 population-specific differences occurring in only one sex, and 22 sex-specific differences. Many of these differentially expressed mRNAs were novel, or not previously hypothesized to play a role in response to contaminant exposure. In addition, the results indicate that the effect of contaminated sediment exposure on the expression of a large proportion of the differentially expressed mRNAs was dependent on the sex of the fish.

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.

Genotoxicological studies in aquatic organisms: an overview

Substantial progress has been made in the lat two decades to evaluate the impact of physical and chemical genotoxins in aquatic organisms. This overview (a) summarises the major high lights in this stimulating area of research, (b) compares the developments in this field with the developments in mammalian genotoxicological studies, where appropriate, (c) introduces 18 different articles presented in this special issue of Mutation Research in the backdrop of main advances and , (d) hypothesises on future directions of research in this exciting field.