Resistance to contaminants in North American fish populations

Comparative genomic damage among three shark species with different habits: Sublethal impacts of human origin in a protected island environment in the South Atlantic

Elasmobranchs can bioaccumulate and biomagnify pollutants. However, few studies are directed to the effects of pollutants on the health of these animals, and in most cases, are limited to the analysis of biochemical markers. Thus, the incidence of genomic damage among shark species inhabiting a protected ocean island in the South Atlantic was investigated in association with the analysis of pollutants in seawater sample. High levels of genomic damage were identified, especially in Negaprion brevirostris and Galeocerdo cuvier, in addition to interspecific variations that may be related to characteristics such as animal size, metabolism and habits. High concentrations of Surfactants were observed in seawater sample, in addition to low concentrations of Cadmium, Lead, Copper, Chromium, Zinc, Manganese, and Mercury. The results evidenced the potential of shark species as a bioindicator of environmental quality and allowed assessing the anthropic impact on the archipelago, which currently drives its economy through tourism.

Contamination characteristics of alkyl polycyclic aromatic hydrocarbons in dust and topsoil collected from Huaibei Coalfield, China

Alkyl polycyclic aromatic hydrocarbons (APAHs) are more toxic and persistent than their parent compounds. In this study, the concentrations of polycyclic aromatic compounds (PACs) in dust, topsoil and coal gangue from Huaibei Coal mine, China were analyzed by gas chromatography-mass spectrometry, confirming APAHs were the dominant pollutants. The mean concentrations of APAHs were substantially higher than those of 16 PAHs in both dust and topsoil. The mean concentration of APAHs in dust was 9197 µg kg^-1, accounting for 80% of the total mean concentration of PACs. The mean concentration of APAHs in topsoil was 2835 µg kg^-1, accounting for 77% of the mean concentration of PACs. Alkyl naphthalenes and alkyl phenanthrenes were the primary pollutants in APAHs. Their mean concentrations in dust and topsoil were 7782 µg kg^-1 and 2333 µg kg^-1, respectively. This accounted for 85% and 82% of the concentration of APAHs, respectively. Additionally, low-molecular-weight APAHs dominated the PACs of the coal mine, exhibiting petrogenic characteristics; distribution of C1-C4 NAP and C1-C4 PHE exhibited "bell shape" pattern indicated as petrogenic source. Source identification indicated that the PACs were mainly derived from petrogenic sources and vehicle emissions, followed by biomass and coal burning. Fingerprinting information of dust and topsoil were consistent with coal gangue, indicating that PACs are most likely derived from coal gangue. Coalfields comparable to our study area are widely distributed in China. Therefore, investigating PAC pollution derived from coal gangue warrants further attention.

Correlations between oxidative DNA damage and formation of hepatic tumours in two flatfish species from contaminated environments

Many species in aquatic environments face increased exposure to oncogenic pollution due to anthropogenic environmental change which can lead to higher cancer prevalence. The mechanistic relationship connecting environmental pollution and cancer is multi-factorial and poorly understood, and the specific mechanisms are so far still uncharacterized. One potential mediator between pollutant exposure and cancer is oxidative damage to DNA. We conducted a study in the field with two flatfish species, European flounder (Platichthys flesus L.) and common dab (Limanda limanda L.) with overlapping distribution and similar ecological niche, to investigate if the link between oncogenic pollutants and cancer described in ecotoxicological literature could be mediated by oxidative DNA damage. This was not the case for flounders as neither polycyclic aromatic hydrocarbon (PAH) bile metabolites nor metallic trace element concentrations were related to oxidative DNA damage measurements. However, dabs with higher PAH concentrations did exhibit increased oxidative damage. High oxidative DNA damage also did not predict neoplasm occurrence, rather, healthy individuals tended to have higher oxidative damage measurements compared to fishes with pre-neoplastic tumours. Our analyses showed that flounders had lower concentrations of PAH bile metabolites, suggesting that compared to dab this species is less exposed or better at eliminating these contaminants.

Cardiac physiology and metabolic gene expression during late organogenesis among F. heteroclitus embryo families from crosses between pollution-sensitive and -resistant parents

BACKGROUND

The teleost fish Fundulus heteroclitus inhabit estuaries heavily polluted with persistent and bioaccumulative chemicals. While embryos of parents from polluted sites are remarkably resistant to toxic sediment and develop normally, embryos of parents from relatively clean estuaries, when treated with polluted sediment extracts, are developmentally delayed, displaying deformities characteristic of pollution-induced embryotoxicity. To gain insight into parental effects on sensitive and resistant phenotypes during late organogenesis, we established sensitive, resistant, and crossed embryo families using five female and five male parents from relatively clean and predominantly PAH-polluted estuaries each, measured heart rates, and quantified individual embryo expression of 179 metabolic genes.

RESULTS

Pollution-induced embryotoxicity manifested as morphological deformities, significant developmental delays, and altered cardiac physiology was evident among sensitive embryos resulting from crosses between females and males from relatively clean estuaries. Significantly different heart rates among several geographically unrelated populations of sensitive, resistant, and crossed embryo families during late organogenesis and pre-hatching suggest site-specific adaptive cardiac physiology phenotypes relative to pollution exposure. Metabolic gene expression patterns (32 genes, 17.9%, at p < 0.05; 11 genes, 6.1%, at p < 0.01) among the embryo families indicate maternal pollutant deposition in the eggs and parental effects on gene expression and metabolic alterations.

CONCLUSION

Heart rate differences among sensitive, resistant, and crossed embryos is a reliable phenotype for further explorations of adaptive mechanisms. While metabolic gene expression patterns among embryo families are suggestive of parental effects on several differentially expressed genes, a definitive adaptive signature and metabolic cost of resistant phenotypes is unclear and shows unexpected sensitive-resistant crossed embryo expression profiles. Our study highlights physiological and metabolic gene expression differences during a critical embryonic stage among pollution sensitive, resistant, and crossed embryo families, which may contribute to underlying resistance mechanisms observed in natural F. heteroclitus populations living in heavily contaminated estuaries.

Association between Chromosome 4 and mercury accumulation in muscle of the three-spined stickleback (Gasterosteus aculeatus)

Anthropogenic stressors, such as pollutants, act as selective factors that can leave measurable changes in allele frequencies in the genome. Metals are of particular concern among pollutants, because of interference with vital biological pathways. We use the three-spined stickleback as a model for adaptation to mercury pollution in natural populations. We collected sticklebacks from 21 locations in Flanders (Belgium), measured the accumulated levels of mercury in the skeletal muscle tissue, and genotyped the fish by sequencing (GBS). The spread of muscle mercury content across locations was considerable, ranging from 21.5 to 327 ng/g dry weight (DW). We then conducted a genome-wide association study (GWAS) between 28,450 single nucleotide polymorphisms (SNPs) and the accumulated levels of mercury, using different approaches. Based on a linear mixed model analysis, the GWAS yielded multiple hits with a single top hit on Chromosome 4, with eight more SNPs suggestive of association. A second approach, a latent factor mixed model analysis, highlighted one single SNP on Chromosome 11. Finally, an outlier test identified one additional SNP on Chromosome 4 that appeared under selection. Out of all ten SNPs we identified as associated with mercury in muscle, three SNPs all located on Chromosome 4 and positioned within a 2.5 kb distance of an annotated gene. Based on these results and the genome coverage of our SNPs, we conclude that the selective effect of mercury pollution in Flanders causes a significant association with at least one locus on Chromosome 4 in three-spined stickleback.

Resistance to Cyp3a induction by polychlorinated biphenyls, including non-dioxin-like PCB153, in gills of killifish (Fundulus heteroclitus) from New Bedford Harbor

Previous reports suggested that non-dioxin-like (NDL) PCB153 effects on cytochrome P450 3A (Cyp3a) expression in Atlantic killifish (Fundulus heteroclitus) gills differed between F0 generation fish from a PCB site (New Bedford Harbor; NBH) and a reference site (Scorton Creek; SC). Here, we examined effects of PCB153, dioxin-like (DL) PCB126, or a mixture of both, on Cyp3a56 mRNA in killifish generations removed from the wild, without environmental PCB exposures. PCB126 effects in liver and gills differed between populations, as expected. Gill Cyp3a56 was not affected by either congener in NBH F2 generation fish, but was induced by PCB153 in SC F1 fish, with females showing a greater response. PCB153 did not affect Cyp3a56 in liver of either population. Results suggest a heritable resistance to NDL-PCBs in killifish from NBH, in addition to that reported for DL PCBs. Induction of Cyp3a56 in gills may be a biomarker of exposure to NDL PCBs in fish populations that are not resistant to PCBs.

Linking pollution and cancer in aquatic environments: A review

Due to the interconnectedness of aquatic ecosystems through the highly effective marine and atmospheric transport routes, all aquatic ecosystems are potentially vulnerable to pollution. Whilst links between pollution and increased mortality of wild animals have now been firmly established, the next steps should be to focus on specific physiological pathways and pathologies that link pollution to wildlife health deterioration. One of the pollution-induced pathologies that should be at the centre of attention in ecological and evolutionary research is cancer, as anthropogenic contamination has resulted in a rapid increase of oncogenic substances in natural habitats. Whilst wildlife cancer research is an emerging research topic, systematic reviews of the many case studies published over the recent decades are scarce. This research direction would (1) provide a better understanding of the physiological mechanisms connecting anthropogenic pollution to oncogenic processes in non-model organisms (reducing the current bias towards human and lab-animal studies in cancer research), and (2) allow us to better predict the vulnerability of different wild populations to oncogenic contamination. This article combines the information available within the scientific literature about cancer occurrences in aquatic and semi-aquatic species. For the first aim, we use available knowledge from aquatic species to suggest physiological mechanisms that link pollution and cancer, including main metabolic detoxification pathways, oxidative damage effects, infections, and changes to the microbiome. For the second aim, we determine which types of aquatic animals are more vulnerable to pollution-induced cancer, which types of pollution are mainly associated with cancer in aquatic ecosystems, and which types of cancer pollution causes. We also discuss the role of migration in exposing aquatic and semi-aquatic animals to different oncogenic pollutants. Finally, we suggest novel research avenues, including experimental approaches, analysis of the effects of pollutant cocktails and long-term chronic exposure to lower levels of pollutants, and the use of already published databases of gene expression levels in animals from differently polluted habitats.

Heart development in two populations of Atlantic killifish (Fundulus heteroclitus) following exposure to a polycyclic aromatic hydrocarbon mixture

Historic industrial pollution of the Elizabeth River, Virginia resulted in polycyclic aromatic hydrocarbon (PAH) contamination in sediments. Atlantic killifish (Fundulus heteroclitus) inhabiting the Atlantic Wood (AW) industrial site adapted to complex PAH mixture at this Superfund site. Their embryos have proved highly resistant to cardiac abnormalities indicative of PAH toxicity. In this study, embryos spawned from adults collected at AW and King's Creek (KC), a reference site, were exposed at 24 h post fertilization (hpf) to Elizabeth River Sediment Extract (ERSE), a complex PAH mixture, in a range of concentrations (0, 5.04, 50.45, 100.90, 151.35, or 252.25 µg/L total PAHs). Embryos were processed for histology at 144 hpf to enable evaluations of hearts at tissue and cellular levels. Morphometry and severity scoring were used to evaluate the extent of alterations. Unexposed embryos were similar in both populations. ERSE exposure resulted in multiple changes to hearts of KC embryos but not AW. Alterations were particularly evident in KC embryos exposed to concentrations above 1% ERSE (50.45 µg/L), which had thinner ventricular walls and larger pericardial edema. Individuals with moderate pericardial edema maintained arrangement and proximity of heart chambers, but changes were seen in ventricular myocytes. Severe pericardial edema was prevalent in exposed KC embryos and typically resulted in tube heart formation. Ventricles of tube hearts had very thin walls composed of small, basophilic cells and lacked trabeculae. Edematous pericardial fluid contained small amounts of proteinaceous material, as did controls, and was free of cells. This fluid was primarily unstained, suggesting water influx due to increased permeability. The use of histological approaches provided more specific detail for tissue and cellular effects in hearts of embryos exposed to PAHs and enabled understanding of potential links to later life effects of early life exposure.

Investigation into Adaptation in Genes Associated with Response to Estrogenic Pollution in Populations of Roach (Rutilus rutilus) Living in English Rivers

Exposure of male fish to estrogenic substances from wastewater treatment works (WwTWs) results in feminization and reduced reproductive fitness. Nevertheless, self-sustaining populations of roach (Rutilus rutilus) inhabit river stretches polluted with estrogenic WwTW effluents. In this study, we examine whether such roach populations have evolved adaptations to tolerate estrogenic pollution by comparing frequency differences in single-nucleotide polymorphisms (SNPs) between populations sampled from rivers receiving either high- or low-level WwTW discharges. SNPs within 36 "candidate" genes, selected for their involvement in estrogenic responses, and 120 SNPs in reference genes were genotyped in 465 roaches. There was no evidence for selection in highly estrogen-dependent candidate genes, including those for the estrogen receptors, aromatases, and vitellogenins. The androgen receptor (ar) and cytochrome P450 1A genes were associated with large shifts in allele frequencies between catchments and in individual populations, but there is no clear link to estrogen pollution. Selection at ar in the effluent-dominated River Lee may have resulted from historical contamination with endocrine-disrupting pesticides. Critically, although our results suggest population-specific selection including at genes related to endocrine disruption, there was no strong evidence that the selection resulted from exposure to estrogen pollution.

Immunity and health of two wild marine fishes naturally exposed to anthropogenic pollution

There are increasing global concerns of the alarming pollution impacts on marine life, thus it is becoming essential to generate reliable tools to monitor and understand the effects of these impacts on aquatic organisms. We performed a field study assessing how exposure to anthropogenic pollution impacts immunological and health-state parameters and parasite infection of a wild marine fish, the Brazilian sandperch Pinguipes brasilianus. Then we compared this information to previously published data of a sympatric species, the Patagonian rockfish Sebastes oculatus inhabiting the same polluted and pristine areas. The field study revealed that exposed P. brasilianus showed chronic stress, poor immune condition and higher prevalence and abundance of acanthocephalan parasites. By comparing these former results with already published in S. oculatus, we concluded that, although both species exhibited physiological alterations associate to inhabiting sites exposed to pollution, their specific immunological and health-state responses differed. Our results demonstrate that Patagonian reef-fish assemblages inhabiting sites exposed to pollutant are being affected in their immune and heath condition, which could potentially result in higher susceptibility to disease and in turn population decline. These findings highlight the necessity of more studies incorporating interspecific comparisons to assess variation in fish susceptibility in an ecoimmunotoxicological context and get a more profound understanding of anthropogenic impacts on wildlife.

In Silico Site-Directed Mutagenesis Informs Species-Specific Predictions of Chemical Susceptibility Derived From the Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) Tool

Chemical hazard assessment requires extrapolation of information from model organisms to all species of concern. The Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool was developed as a rapid, cost-effective method to aid cross-species extrapolation of susceptibility to chemicals acting on specific protein targets through evaluation of protein structural similarities and differences. The greatest resolution for extrapolation of chemical susceptibility across species involves comparisons of individual amino acid residues at key positions involved in protein-chemical interactions. However, a lack of understanding of whether specific amino acid substitutions among species at key positions in proteins affect interaction with chemicals made manual interpretation of alignments time consuming and potentially inconsistent. Therefore, this study used in silico site-directed mutagenesis coupled with docking simulations of computational models for acetylcholinesterase (AChE) and ecdysone receptor (EcR) to investigate how specific amino acid substitutions impact protein-chemical interaction. This study found that computationally derived substitutions in identities of key amino acids caused no change in protein-chemical interaction if residues share the same side chain functional properties and have comparable molecular dimensions, while differences in these characteristics can change protein-chemical interaction. These findings were considered in the development of capabilities for automatically generated species-specific predictions of chemical susceptibility in SeqAPASS. These predictions for AChE and EcR were shown to agree with SeqAPASS predictions comparing the primary sequence and functional domain sequence of proteins for more than 90% of the investigated species, but also identified dramatic species-specific differences in chemical susceptibility that align with results from standard toxicity tests. These results provide a compelling line of evidence for use of SeqAPASS in deriving screening level, species-specific, susceptibility predictions across broad taxonomic groups for application to human and ecological hazard assessment.

Quantitative Response-Response Relationships Linking Aromatase Inhibition to Decreased Fecundity are Conserved Across Three Fishes with Asynchronous Oocyte Development

Quantitative adverse outcome pathways (qAOPs) describe quantitative response-response relationships that can predict the probability or severity of an adverse outcome for a given magnitude of chemical interaction with a molecular initiating event. However, the taxonomic domain of applicability for these predictions is largely untested. The present study began defining this applicability for a previously described qAOP for aromatase inhibition leading to decreased fecundity developed using data from fathead minnow (Pimephales promelas). This qAOP includes quantitative response-response relationships describing plasma 17β-estradiol (E2) as a function of plasma fadrozole, plasma vitellogenin (VTG) as a function of plasma E2, and fecundity as a function of plasma VTG. These quantitative response-response relationships simulated plasma E2, plasma VTG, and fecundity measured in female zebrafish (Danio rerio) exposed to fadrozole for 21 days but not these responses measured in female Japanese medaka (Oryzias latipes). However, Japanese medaka had different basal levels of plasma E2, plasma VTG, and fecundity. Normalizing basal levels of each measurement to equal those of female fathead minnow enabled the relationships to accurately simulate plasma E2, plasma VTG, and fecundity measured in female Japanese medaka. This suggests that these quantitative response-response relationships are conserved across these three fishes when considering relative change rather than absolute measurements. The present study represents an early step toward defining the appropriate taxonomic domain of applicability and extending the regulatory applications of this qAOP.

Mixed evidence for adaptation to environmental pollution

Adaptation to pollution has been studied since the first observations of heavy metal tolerance in plants decades ago. To document micro-evolutionary responses to pollution, researchers have used phenotypic, molecular genetics, and demographic approaches. We reviewed 258 articles and evaluated the evidence for adaptive responses following exposure to a wide range of pollutants, across multiple taxonomic groups. We also conducted a meta-analysis to calculate the magnitude of phenotypic change in invertebrates in response to metal pollution. The majority of studies that reported differences in responses to pollution were focused on phenotypic responses at the individual level. Most of the studies that used demographic assays in their investigations found that negative effects induced by pollution often worsened over multiple generations. Our meta-analysis did not reveal a significant relationship between metal pollution intensity and changes in the traits studied, and this was probably due to differences in coping responses among different species, the broad array of abiotic and biotic factors, and the weak statistical power of the analysis. We found it difficult to make broad statements about how likely or how common adaptation is in the presence of environmental contamination. Ecological and evolutionary responses to contamination are complex, and difficult to interpret in the context of taxonomic, and methodological biases, and the inconsistent set of approaches that have been used to study adaptation to pollution in the laboratory and in the field. This review emphasizes the need for: (a) long-term monitoring programs on exposed populations that link demography to phenotypic, genetic, and selection assays; (b) the use of standardized protocols across studies especially for similar taxa. Approaches that combine field and laboratory studies offer the greatest opportunity to reveal the complex eco-evolutionary feedback that can occur under selection imposed by pollution.

An affordable and automated imaging approach to acquire highly resolved individual data-an example of copepod growth in response to multiple stressors

Individual trait variation is essential for populations to cope with multiple stressors and continuously changing environments. The immense number of possible stressor combinations and the influence of phenotypic variation makes experimental testing for effects on organisms challenging. The acquisition of such data requires many replicates and is notoriously laborious. It is further complicated when responses occur over short time periods. To overcome such challenges, we developed an automated imaging platform to acquire temporally highly resolved individual data. We tested this platform by exposing copepods to a combination of a biotic stressor (predator cues) and a toxicant (copper) and measured the growth response of individual copepods. We tested the automatically acquired data against published manually acquired data with much lower temporal resolution. We find the same general potentiating effects of predator cues on the adverse effects of copper, and the influence of an individual’s clutch identity on its ability to resist stress, between the data obtained from low and high temporal resolution. However, when using the high temporal resolution, we also uncovered effects of clutch ID on the timing and duration of stage transitions, which highlights the importance of considering phenotypic variation in ecotoxicological testing. Phenotypic variation is usually not acknowledged in ecotoxicological testing. Our approach is scalable, affordable, and adjustable to accommodate both aquatic and terrestrial organisms, and a wide range of visually detectable endpoints. We discuss future extensions that would further widen its applicability.

The relation between polychlorinated biphenyls and population metrics of 4 species of fish from the upper Hudson River, New York, USA

In the upper Hudson River, New York, USA, fish were exposed to polychlorinated biphenyls (PCBs) from the 1940s to 1977, and PCBs still persist in this environment. Yellow perch (Perca flavescens), brown bullhead (Ameiurus nebulosus), smallmouth bass (Micropterus dolomieu), and largemouth bass (M. salmoides) were collected annually from 2003 to 2009 from 1 control site upstream of the PCB discharge locations and from 2 sites downstream from where PCBs were released. Fish PCB concentrations were estimated, and 3 population metrics were examined: 1) relative abundance, 2) weight-to-length ratio, and 3) growth. Normalized lipid-based PCB concentrations at the 2 PCB exposure pools averaged approximately 100 to 600 μg/g. Estimated relative abundances with electrofishing were higher for largemouth bass, smallmouth bass, and brown bullhead at PCB exposure sites compared to the control site; but yellow perch were more abundant at the control site. Weight to length ratios varied among sites and species, but no consistent pattern was evident in relation to PCBs at the population level or for individual fish. Growth rates for yellow perch and brown bullhead were similar among sites. Largemouth bass growth was slightly higher at the control site compared to the 2 PCB sites, but smallmouth bass growth was much higher at the PCB sites compared to the control site. We could not detect any relation or influence of PCBs on the 3 population metrics that we examined. the present results corroborated those of previous investigations concerning the effects of PCBs on fishes. We recommend stronger consideration of the biological impacts of PCBs at the population level when conducting risk assessments. Environ Toxicol Chem 2019;38:329-339. © 2018 SETAC.

Turning natural adaptations to oncogenic factors into an ally in the war against cancer

Both field and experimental evolution studies have demonstrated that organisms naturally or artificially exposed to environmental oncogenic factors can, sometimes rapidly, evolve specific adaptations to cope with pollutants and their adverse effects on fitness. Although numerous pollutants are mutagenic and carcinogenic, little attention has been given to exploring the extent to which adaptations displayed by organisms living in oncogenic environments could inspire novel cancer treatments, through mimicking the processes allowing these organisms to prevent or limit malignant progression. Building on a substantial knowledge base from the literature, we here present and discuss this progressive and promising research direction, advocating closer collaboration between the fields of medicine, ecology, and evolution in the war against cancer.

Impacts of the Deepwater Horizon oil spill evaluated using an end-to-end ecosystem model

We use a spatially explicit biogeochemical end-to-end ecosystem model, Atlantis, to simulate impacts from the Deepwater Horizon oil spill and subsequent recovery of fish guilds. Dose-response relationships with expected oil concentrations were utilized to estimate the impact on fish growth and mortality rates. We also examine the effects of fisheries closures and impacts on recruitment. We validate predictions of the model by comparing population trends and age structure before and after the oil spill with fisheries independent data. The model suggests that recruitment effects and fishery closures had little influence on biomass dynamics. However, at the assumed level of oil concentrations and toxicity, impacts on fish mortality and growth rates were large and commensurate with observations. Sensitivity analysis suggests the biomass of large reef fish decreased by 25% to 50% in areas most affected by the spill, and biomass of large demersal fish decreased even more, by 40% to 70%. Impacts on reef and demersal forage caused starvation mortality in predators and increased reliance on pelagic forage. Impacts on the food web translated effects of the spill far away from the oiled area. Effects on age structure suggest possible delayed impacts on fishery yields. Recovery of high-turnover populations generally is predicted to occur within 10 years, but some slower-growing populations may take 30+ years to fully recover.

Adaptive capabilities and fitness consequences associated with pollution exposure in fish

Many fish populations are exposed to harmful levels of chemical pollution and selection pressures associated with these exposures have led to the evolution of tolerance. Our understanding of the physiological basis for these adaptations is limited, but they are likely to include processes involved with the absorption, distribution, metabolism and/or excretion of the target chemical. Other potential adaptive mechanisms include enhancements in antioxidant responses, an increased capacity for DNA and/or tissue repair and alterations to the life cycle of fish that enable earlier reproduction. Analysis of single-nucleotide polymorphism frequencies has shown that tolerance to hydrocarbon pollutants in both marine and estuarine fish species involves alteration in the expression of the xenobiotic metabolism enzyme CYP1A. In this review, we present novel data showing also that variants of the CYP1A gene have been under selection in guppies living in Trinidadian rivers heavily polluted with crude oil. Potential costs associated with these adaptations could reduce fitness in unpolluted water conditions. Integrating knowledge of local adaptation to pollution is an important future consideration in conservation practices such as for successful restocking, and improving connectivity within river systems.This article is part of the themed issue 'Human influences on evolution, and the ecological and societal consequences'.

Micronuclei and other erythrocyte nuclear abnormalities in fishes from the Great Lakes Basin, USA

Biological markers (biomarkers) sensitive to genotoxic and mutagenic contamination in fishes are widely used to identify exposure effects in aquatic environments. The micronucleus assay was incorporated into a suite of indicators to assess exposure to genotoxic and mutagenic contamination at five Great Lakes Areas of Concern (AOCs), as well as one non-AOC (reference) site. The assay allowed enumeration of micronuclei as well as other nuclear abnormalities for both site and species comparisons. Erythrocyte abnormality data was also compared to skin and liver tumor prevalence and hepatic transcript abundance. Erythrocyte abnormalities were observed at all sites with variable occurrence and severity among sites and species. Benthic-oriented brown bullhead (Ameiurus nebulosus) and white sucker (Catostomus commersonii) expressed lower rates of erythrocyte abnormalities, but higher rates of skin and liver neoplasms, when compared to pelagic-oriented largemouth bass (Micropterus salmoides) or smallmouth bass (Micropterus dolomieu) at the same site. The reduced erythrocyte abnormalities, increased transcript abundance associated with Phase I and II toxicant responsive pathways, and increased neoplastic lesions among benthic-oriented taxa may indicate the development of contaminant resistance of these species to more acute effects. Environ. Mol. Mutagen. 58:570-581, 2017. © 2017 This article is a U.S. Government work and is in the public domain in the USA. Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

RNA seq- and DEG-based comparison of developmental toxicity in fish embryos of two species exposed to Iranian heavy crude oil

To determine and compare the toxic effects of Iranian heavy crude oil (IHCO) on the embryonic development of two fish species, we examined transcriptome profiles using RNA-seq. The assembled contigs were 66,070 unigenes in olive flounder embryos and 76,498 unigenes in spotted seabass embryos. In the differential gene expression (DEG) profiles, olive flounder embryos showed different up- and down-regulated patterns than spotted seabass embryos in response to fresh IHCO (FIHCO) and weathered IHCO (WIHCO). In this work, we categorized DEG profiles into six pathways: ribosome, oxidative phosphorylation, Parkinson's disease, Alzheimer's disease, Huntington's disease, and cardiac muscle contraction, validating the expression patterns of 13 DEGs using real-time quantitative RT-PCR. The expression of the CYP1A, CYP1B1, and CYP1C1 genes in spotted seabass embryos was higher than in olive flounder embryos, whereas genes related to cell processing, development, and the immune system showed the opposite trend. Orthologous gene cluster analysis showed that olive flounder embryos were sensitive (fold change of genes with cutoff P<0.05) to both FIHCO and WIHCO, but spotted seabass embryos exhibited higher sensitivity to WIHCO than FIHCO, indicating that species-specific differences are likely to be reflected in population levels after oil spills. Overall, our study provides new insight on the different embryonic susceptibilities of two marine fish species to FIHCO and WIHCO and a better understanding of the underlying molecular mechanisms via RNA-seq and DEGs.

Hepatic transcriptional responses to copper in the three-spined stickleback are affected by their pollution exposure history

Some fish populations inhabiting contaminated environments show evidence of increased chemical tolerance, however the mechanisms contributing to this tolerance, and whether this is heritable, are poorly understood. We investigated the responses of two populations of wild three-spined stickleback (Gasterosteus aculeatus) with different histories of contaminant exposure to an oestrogen and copper, two widespread aquatic pollutants. Male stickleback originating from two sites, the River Aire, with a history of complex pollution discharges, and Siblyback Lake, with a history of metal contamination, were depurated and then exposed to copper (46μg/L) and the synthetic oestrogen ethinyloestradiol (22ng/L). The hepatic transcriptomic response was compared between the two populations and to a reference population with no known history of exposure (Houghton Springs, Dorset). Gene responses included those typical for both copper and oestrogen, with no discernable difference in response to oestrogen between populations. There was, however, some difference in the magnitude of response to copper between populations. Siblyback fish showed an elevated baseline transcription of genes encoding metallothioneins and a lower level of metallothionein induction following copper exposure, compared to those from the River Aire. Similarly, a further experiment with an F1 generation of Siblyback fish bred in the laboratory found evidence for elevated transcription of genes encoding metallothioneins in unexposed fish, together with an altered transcriptional response to 125μg/L copper, compared with F1 fish originating from the clean reference population exposed to the same copper concentration. These data suggest that the stickleback from Siblyback Lake have a differential response to copper, which is inherited by the F1 generation in laboratory conditions, and for which the underlying mechanism may include an elevation of baseline transcription of genes encoding metallothioneins. The genetic and/or epigenetic mechanisms contributing to this inherited alteration of metallothionein transcription have yet to be established.

Evolutionary toxicology in an omics world

Evolutionary toxicology is a young field that has grown rapidly in the past two decades. The potential of this field comes from the ability to link chemical contamination to multigenerational and population-wide effects in various species. The advancements and rapidly decreasing costs of -omic tools are improving the power and resolution of evolutionary toxicology studies. In this manuscript, we aim to address the trajectories and perspectives for conducting evolutionary toxicology studies with -omic approaches. We discuss the complementarity of using multiple -omic tools (genomics, eDNA, transcriptomics, proteomics, and metabolomics) for utility in understanding the toxicological relevance of adaptive responses in populations. In addition, we discuss phenotypic plasticity and its relevance to transcriptomic studies in toxicology. As evolutionary toxicology grows and expands its capacity to link toxicology with population-wide end points, we emphasize the applications of such studies in answering questions about ecological and population health, as well as future applicability to regulation. Thus, we aim to emphasize the enormous potential for evolutionary toxicology in an -omics world and give perspectives on the directions of future investigations.

Cross-resistance in Gulf killifish (Fundulus grandis) populations resistant to dioxin-like compounds

The Houston Ship Channel (HSC) in Houston, Texas is an aquatic environment with a long history of contamination, including polychlorinated dibenzodioxins (PCDD), polychlorinated dibenzofurans (PCDF), polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), and heavy metals. Populations of Gulf killifish (Fundulus grandis) from the HSC have adapted to resist developmental cardiac deformities caused by dioxin-like compounds (DLCs). Contaminants in the HSC have acted as a strong selective pressure on resident Gulf killifish populations. Rapid adaptation can lead to fitness costs, some as a direct result of the mechanisms involved in the adaptive process, whereas other adaptations may be more general. To explore potential fitness costs, we evaluated two Gulf killifish populations with documented resistance to DLC-induced cardiac teratogenesis (Patrick Bayou and Vince Bayou), and one previously characterized reference population (Gangs Bayou). We also characterized a previously unstudied population from Galveston Bay as an additional reference population (Smith Point). We tested the sensitivity of F1 larvae from these four populations to two classes of pesticides (pyrethroid (permethrin) and carbamate (carbaryl)) and two model pro-oxidants (tert-butyl hydroquinone (tBHQ) and tert-butyl hydroperoxide (tBOOH)). In addition, we explored their responses to hypoxia and measured resting metabolic rates (M.O2). Both adapted populations were cross-resistant to the toxicity of carbaryl and both pro-oxidants tested. There were no population differences in sensitivity to permethrin. On the other hand, one reference population (Gangs Bayou) was less sensitive to hypoxia, and maintained a lower M.O2 . However, there were no differences in hypoxia tolerance or resting metabolic rate between the second reference and the two adapted populations. This investigation emphasizes the importance of including multiple reference populations to clearly link fitness costs or cross-resistance to pollution adaptation, rather than to unrelated environmental or ecological differences. When compared to previous literature on adapted populations of Fundulus heteroclitus, we see a mixture of similarities and differences, suggesting that F. grandis adapted phenotypes likely involve multiple mechanisms, which may not be completely consistent among adapted populations.

Can ecological history influence response to pollutants? Transcriptomic analysis of Manila clam collected in different Venice lagoon areas and exposed to heavy metal

Chronic exposure to environmental pollutants can exert strong selective pressures on natural populations, favoring the transmission over generations of traits that enable individuals to survive and thrive in highly impacted environments. The lagoon of Venice is an ecosystem subject to heavy anthropogenic impact, mainly due to the industrial activities of Porto Marghera (PM), which led to a severe chemical contamination of soil, groundwater, and sediments. Gene expression analysis on wild Manila clams collected in different Venice lagoon areas enabled to identify differences in gene expression profiles between clams collected in PM and those sampled in clean areas, and the definition of molecular signatures of chemical stress. However, it remains largely unexplored to which extent modifications of gene expression patterns persists after removing the source of contamination. It is also relatively unknown whether chronic exposure to xenobiotics affects the response to other chemical pollutants. To start exploring such issues, in the present study a common-garden experiment was coupled with transcriptomic analysis, to compare gene expression profiles of PM clams with those of clams collected in the less impacted area of Chioggia (CH) during a period under the same control conditions. Part of the two experimental groups were also exposed to copper for seven days to assess whether different "ecological history" does influence response to such pollutant. The results obtained suggest that the chronic exposure to chemical pollution generated a response at the transcriptional level that persists after removal for the contaminated site. These transcriptional changes are centered on key biological processes, such as defense against either oxidative stress or tissue/protein damage, and detoxification, suggesting an adaptive strategy for surviving in the deeply impacted environment of Porto Marghera. On the other hand, CH clams appeared to respond more effectively to copper exposure than PM animals, proposing that chronic exposure to chemical toxicants either lowers the sensitivity to additional toxicants or blunts the capacity to respond to novel chemical challenges in PM clams.

Genetic variation and population structure of American mink Neovison vison from PCB-contaminated and non-contaminated locales in eastern North America

American mink Neovison vison may be particularly vulnerable to toxicities of persistent contaminants such as PCBs because of their aquatic-based diet, position near the top of the food web, and small deme sizes. Furthermore, ranched mink are sensitive to reproductive toxicities of fish diets from PCB-polluted sites. The upper Hudson River is highly contaminated with PCBs and previous studies have shown elevated hepatic burdens of total and coplanar PCBs in mink collected near the river compared with those from more distant locales in New York and elsewhere. We hypothesized that bioaccumulation of PCBs in Hudson River mink has reduced their levels of genetic diversity or altered their genetic population structure. To address this, we conducted microsatellite DNA analysis on collections made in proximity to and from more distant locales in the Hudson River watershed, elsewhere in New York State, and at other sites in eastern North America including New Brunswick, four locales in Ontario, multiple drainages in Maine, and two ecoregions in Rhode Island. We did not find reduced genetic diversity at the individual or population levels in mink collected near (<6 km) to PCB hotspots in the Hudson River nor evidence of altered population structure. Consistent with their distribution in small localized and isolated demes, we did find significant genetic population structure among many mink collections in New York State and elsewhere. Depending on the analytical approach used, genetically distinct populations numbered between 16 when using STRUCTURE to 19-20 when using Exact G tests, F ST, or AMOVA analyses. Genetically distinct population units were found among major ecoregions and minor ecoregions in New York State, among different hydrologic subunits within the Hudson River watershed, among spatially separate locales in Ontario, and among most watersheds in Maine. However, despite this localization and potential heightened impact of stressors, genetic diversity and genetic population structure in mink does not seem to be affected by their bioaccumulation of high levels of PCBs of Hudson River origin.

Sensitivity of shovelnose sturgeon (Scaphirhynchus platorynchus) and pallid sturgeon (S. albus) early life stages to 3,3',4,4',5-pentachlorobiphenyl and 2,3,7,8-tetrachlorodibenzo-P-dioxin exposure

Concern exists that polychlorinated biphenyls (PCBs) may be contributing to the current decline of shovelnose sturgeon (Scaphirhynchus platorynchus) and the US federally endangered pallid sturgeon (Scaphirhynchus albus). Waterborne exposures with newly fertilized eggs were used to assess developmental and morphological effects of 2 of the most potent aryl hydrocarbon receptor (AhR) agonists, 3,3',4,4',5-pentachlorobiphenyl (PCB-126) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), on early life stage shovelnose and pallid sturgeon. No dose-related effects of PCB-126 were observed on percent development or hatch in either species at concentrations as high as 1711 ng/g egg. Effects of TCDD on percent development were not assessed in shovelnose sturgeon. However, percent development was not affected by TCDD in pallid sturgeon, and percent hatch was unaffected by TCDD doses as high as 60 ng/g egg to 81 ng/g egg in either species. Morphological pathologies such as yolk sac edema and craniofacial deformities were typical of AhR agonist exposure and were similar in both species. Calculated PCB-126 50% lethal dose (LD50, 95% fiducial limits) values were 196 ng/g egg (188-203 ng/g) for shovelnose and 159 ng/g egg (122-199 ng/g) for pallid sturgeon. Likewise, calculated TCDD LD50 values were 13 ng/g egg (11-15 ng/g) for shovelnose and 12 ng/g egg (10-14 ng/g) for pallid sturgeon. These LD50 values are among the highest recorded in early life stage fish, suggesting that early life stage Scaphirhynchus sturgeon may be comparatively insensitive to AhR agonists.

Variation in toxicity of a current-use insecticide among resurrected Daphnia pulicaria genotypes

This study examined how genotypes of Daphnia pulicaria from a single population, separated by thousands of generations of evolution in the wild, differ in their sensitivity to a novel anthropogenic stressor. These genotypes were resurrected from preserved resting eggs isolated from sediments belonging to three time periods: 2002-2008, 1967-1977, and 1301-1646 A.D. Toxicity of the organophosphate insecticide chlorpyrifos was determined through a series of acute toxicity tests. There was a significant dose-response effect in all genotypes studied. Moreover, significant variation in toxicity among genotypes within each time period was detected. Importantly, a significant effect of time period on sensitivity to chlorpyrifos was found. Analysis of the median effect concentrations (EC50s) for genotypes within each time period indicated that the 1301-1646 genotypes were 2.7 times more sensitive than the 1967-1977 genotypes. This trend may be partially explained by microevolutionary shifts in response to cultural eutrophication.

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.

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.

The Elizabeth River Story: A Case Study in Evolutionary Toxicology

The Elizabeth River system is an estuary in southeastern Virginia, surrounded by the towns of Chesapeake, Norfolk, Portsmouth, and Virginia Beach. The river has played important roles in U.S. history and has been the location of various military and industrial activities. These activities have been the source of chemical contamination in this aquatic system. Important industries, until the 1990s, included wood treatment plants that used creosote, an oil-derived product that is rich in polycyclic aromatic hydrocarbons (PAH). These plants left a legacy of PAH pollution in the river, and in particular Atlantic Wood Industries is a designated Superfund site now undergoing remediation. Numerous studies examined the distribution of PAH in the river and impacts on resident fauna. This review focuses on how a small estuarine fish with a limited home range, Fundulus heteroclitus (Atlantic killifish or mummichog), has responded to this pollution. While in certain areas of the river this species has clearly been impacted, as evidenced by elevated rates of liver cancer, some subpopulations, notably the one associated with the Atlantic Wood Industries site, displayed a remarkable ability to resist the marked effects PAH have on the embryonic development of fish. This review provides evidence of how pollutants have acted as evolutionary agents, causing changes in ecosystems potentially lasting longer than the pollutants themselves. Mechanisms underlying this evolved resistance, as well as mechanisms underlying the effects of PAH on embryonic development, are also described. The review concludes with a description of ongoing and promising efforts to restore this historic American river.

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

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

Resistance to teratogenesis by F1 and F2 embryos of PAH-adapted Fundulus heteroclitus is strongly inherited despite reduced recalcitrance of the AHR pathway

Atlantic killifish (Fundulus heteroclitus) inhabiting the Atlantic Wood Superfund site on the Elizabeth River (Portsmouth, VA, USA) are exposed to a complex mixture of polycyclic aromatic hydrocarbons (PAHs) from former creosote operations, but are resistant to the acute toxicity and cardiac teratogenesis caused by PAHs. The resistance is associated with a dramatic recalcitrance to induction of cytochrome P450 (CYP1) metabolism enzymes following exposure to aryl hydrocarbon receptor (AHR) agonists, along with an elevated antioxidant response and increased expression of several other xenobiotic metabolism and excretion enzymes. However, the heritability of the resistance in the absence of chemical stressors has been inconsistently demonstrated. Understanding the heritability of this resistance will help clarify the nature of population-level responses to chronic exposure to PAH mixtures and aid in identifying the important mechanistic components of resistance to aryl hydrocarbons. We compared the response of Atlantic Wood F1 and F2 embryos to benzo[k]fluoranthene (BkF), benzo[a]pyrene (BaP), 3,3',4,4',5-pentachlorobiphenyl (PCB-126), and a mixture of BkF and fluoranthene (Fl) to that of F1 embryos of reference site killifish. Resistance to cardiac teratogenesis and induction of CYP mRNA expression and CYP activity was determined. We found that both Atlantic Wood F1 and F2 embryos were highly resistance to cardiac teratogenesis. However, the resistance by Atlantic Wood F2 embryos to induction of CYP mRNA expression and enzyme activity was intermediate between that of Atlantic Wood F1 embryos and reference embryos. These results suggest that resistance to cardiac teratogenesis in Atlantic Wood fish is conferred by multiple factors, not all of which appear to be fully genetically heritable.

Benzo(a)pyrene exposure causes adaptive changes in p53 and CYP1A gene expression in Brown bullhead (Ameiurus nebulosus)

The Brown bullhead (Ameiurus nebulosus) is able to survive and reproduce in high levels of environmentally contaminated areas of the Great Lakes. The purpose of this study was to establish whether there are adaptive genetic/molecular changes occurring in these fish that allow for their survival. Expression of a cell cycle regulator, p53 and the toxin metabolizing protein, CYP1A were measured in liver tissue from bullhead caught from either clean or contaminated areas of Lake Erie and surrounding areas. Wild caught fish and F1 raised offspring (whose parents originated from clean and contaminated sites) were used to measure endogenous gene expression levels. Results revealed that endogenous expression of p53 was on average 6.6× higher in contaminated fish than in fish caught from clean sites. Interestingly, when fed benzo(a)pyrene (BaP)-treated food, p53 expression increased 0.2× in clean fish and decreased 2.6× in contaminated fish. Endogenous CYP1A expression was not detectable in clean fish and low in contaminated fish. Upon exposure to BaP-treated food, CYP1A expression increased in both clean and contaminated fish, although at a higher rate in clean fish. Furthermore, when fish were cleared and then re-exposed to BaP, CYP1A expression increased from basal levels at a higher rate in clean versus contaminated fish. CYP1A and p53 expression in F1 offspring was similar to wild caught fish at the endogenous level and when fed BaP treated food. Results suggest that fish in contaminated regions may be implementing an adaptive response to severe environmental stress by maintaining high expression of p53 and low expression of CYP1A; thus lending increased protection to cells and decreasing the potential amount of carcinogens produced by contaminant metabolism.

Internal effect concentrations of organic substances for early life development of egg-exposed fish

The present study investigates the likelihood that early life development of marine fish from contaminated areas is affected by maternally transferred persistent organic substances (POPs). The common sole (Solea solea) was used as model species. Fertilized eggs were exposed via the water until hatching, 6 days post fertilization. The newly hatched larvae were allowed to develop further under unexposed conditions until the end of the metamorphosis. Effects on the larvae were determined for the dioxin-like polychlorinated biphenyl PCB 126, the technical PCB-mixture Arochlor 1254, polybrominated diphenylethers (PBDEs), and hexabromocyclododecane (HBCD), for an artificial mixture of PCBs and PBDEs, and for 'field mixtures' extracted from sole from the North Sea and the contaminated Western Scheldt estuary. Effect levels were expressed as tissue concentrations in the newly hatched larvae at the end of the exposure period. Exposure to PCBs, PBDEs, and the artificial and field mixtures caused mortality that started to occur shortly after the larvae became free-feeding (10 days post fertilization) and continued to increase until the onset of metamorphosis, 15 days later. The effects induced by the field mixtures correlated well with the ΣPCB concentrations in the tissue of the exposed larvae. No indications were found for synergistic effects or for substantial contribution of other (unknown) substances in the field mixtures. HBCD did not induce toxic effects. As lipid normalized POP levels in fish eggs are in general comparable to the levels in the tissue of the female fish, fish tissue concentrations are indicative of the internal exposure of the developing larvae as a result maternally transferred POPs will occur in the field. In sole from the Western Scheldt estuary POP levels are about twenty times lower than the larval tissue concentration that produced 50 percent early life stage mortality. Levels in North Sea sole are an order of a magnitude lower. At more heavily contaminated sites negative effect of PCBs, especially of those with dioxin-like toxicity can be expected.

Genetic variation at aryl hydrocarbon receptor (AHR) loci in populations of Atlantic killifish (Fundulus heteroclitus) inhabiting polluted and reference habitats

Background

The non-migratory killifish Fundulus heteroclitus inhabits clean and polluted environments interspersed throughout its range along the Atlantic coast of North America. Several populations of this species have successfully adapted to environments contaminated with toxic aromatic hydrocarbon pollutants such as polychlorinated biphenyls (PCBs). Previous studies suggest that the mechanism of resistance to these and other “dioxin-like compounds” (DLCs) may involve reduced signaling through the aryl hydrocarbon receptor (AHR) pathway. Here we investigated gene diversity and evidence for positive selection at three AHR-related loci (AHR1, AHR2, AHRR) in F. heteroclitus by comparing alleles from seven locations ranging over 600 km along the northeastern US, including extremely polluted and reference estuaries, with a focus on New Bedford Harbor (MA, USA), a PCB Superfund site, and nearby reference sites.

Results

We identified 98 single nucleotide polymorphisms within three AHR-related loci among all populations, including synonymous and nonsynonymous substitutions. Haplotype distributions were spatially segregated and F-statistics suggested strong population genetic structure at these loci, consistent with previous studies showing strong population genetic structure at other F. heteroclitus loci. Genetic diversity at these three loci was not significantly different in contaminated sites as compared to reference sites. However, for AHR2 the New Bedford Harbor population had significant F_ST values in comparison to the nearest reference populations. Tests for positive selection revealed ten nonsynonymous polymorphisms in AHR1 and four in AHR2. Four nonsynonymous SNPs in AHR1 and three in AHR2 showed large differences in base frequency between New Bedford Harbor and its reference site. Tests for isolation-by-distance revealed evidence for non-neutral change at the AHR2 locus.

Conclusion

Together, these data suggest that F. heteroclitus populations in reference and polluted sites have similar genetic diversity, providing no evidence for strong genetic bottlenecks for populations in polluted locations. However, the data provide evidence for genetic differentiation among sites, selection at specific nucleotides in AHR1 and AHR2, and specific AHR2 SNPs and haplotypes that are associated with the PCB-resistant phenotype in the New Bedford Harbor population. The results suggest that AHRs, and especially AHR2, may be important, recurring targets for selection in local adaptation to dioxin-like aromatic hydrocarbon contaminants.

Genomic and physiological responses to strong selective pressure during late organogenesis: few gene expression changes found despite striking morphological differences

Background

Adaptations to a new environment, such as a polluted one, often involve large modifications of the existing phenotypes. Changes in gene expression and regulation during critical developmental stages may explain these phenotypic changes. Embryos from a population of the teleost fish, Fundulus heteroclitus, inhabiting a clean estuary do not survive when exposed to sediment extract from a site highly contaminated with polycyclic aromatic hydrocarbons (PAHs) while embryos derived from a population inhabiting a PAH polluted estuary are remarkably resistant to the polluted sediment extract. We exposed embryos from these two populations to surrogate model PAHs and analyzed changes in gene expression, morphology, and cardiac physiology in order to better understand sensitivity and adaptive resistance mechanisms mediating PAH exposure during development.

Results

The synergistic effects of two model PAHs, an aryl hydrocarbon receptor (AHR) agonist (β-naphthoflavone) and a cytochrome P4501A (CYP1A) inhibitor (α-naphthoflavone), caused significant developmental delays, impaired cardiac function, severe morphological alterations and failure to hatch, leading to the deaths of reference embryos; resistant embryos were mostly unaffected. Unexpectedly, patterns of gene expression among normal and moderately deformed embryos were similar, and only severely deformed embryos showed a contrasting pattern of gene expression. Given the drastic morphological differences between reference and resistant embryos, a surprisingly low percentage of genes, 2.24% of 6,754 analyzed, show statistically significant differences in transcript levels during late organogenesis between the two embryo populations.

Conclusions

Our study demonstrates important contrasts in responses between reference and resistant natural embryo populations to synergistic effects of surrogate model PAHs that may be important in adaptive mechanisms mediating PAH effects during fish embryo development. These results suggest that statistically significant changes in gene expression of relatively few genes contribute to the phenotypic changes and large morphological differences exhibited by reference and resistant populations upon exposure to PAH pollutants. By correlating cardiac physiology and morphology with changes in gene expression patterns of reference and resistant embryos, we provide additional evidence for acquired resistance among embryos whose parents live at heavily contaminated sites.

Dioxin pollution disrupts reproduction in male Japanese field mice

Dioxins cause various adverse effects in animals including teratogenesis, induction of drug metabolizing enzymes, tumor promotion, and endocrine disruption. Above all, endocrine disruption is known to disturb reproduction in adult animals and may, also seriously impact their offspring. However, most previous studies have quantified the species-specific accumulation of dioxins, whereas few studies have addressed the physiological impacts of dioxins on wildlife, such as reduced reproductive function. Here we clarify an effect of endocrine disruption caused by dioxins on the Japanese field mouse, Apodemus speciosus. Japanese field mice collected from various sites polluted with dioxins accumulated high concentrations of dioxins in their livers. Some dioxin congeners, especially, 1,2,3,4,6,7,8-heptachlorodibenzo-p-dioxin, 3,3',4,4',5-pentachloro biphenyl, 1,2,3,4,6,7,8-heptachlorodibenzofuran, and octachlorodibenzo-p-dioxin, which showed high biota-soil accumulation factors, contributed to concentration of dioxins in mouse livers with an increase of accumulation of total dioxins. As for physiological effects on the Japanese field mouse, high levels of cytochrome P450 1A1 (CYP1A1) mRNA, a drug metabolizing enzyme induced by dioxins, were found in the livers of mice captured at polluted sites. Furthermore, at such sites polluted with dioxins, increased CYP1A1 expression coincided with reduced numbers of active spermatozoa in mice. Thus, disruption in gametogenesis observed in these mice suggests that dioxins not only negatively impact reproduction among Japanese field mice, but might also act as a kind of selection pressure in a chemically polluted environment.

Population bottlenecks, genetic diversity and breeding ability of the three-spined stickleback (Gasterosteus aculeatus) from three polluted English Rivers

Pollution is a significant environmental pressure on fish populations in both freshwater and marine environments. Populations subjected to chronic exposure to pollutants can experience impacts ranging from altered reproductive capacity to changes in population genetic structure. Few studies, however, have examined the reproductive vigor of individuals within populations inhabiting environments characterized by chronic pollution. In this study we undertook an analysis of populations of three-spined sticklebacks (Gasterosteus aculeatus) from polluted sites, to determine levels of genetic diversity, assess for evidence of historic population genetic bottlenecks and determine the reproductive competitiveness of males from these locations. The sites chosen included locations in the River Aire, the River Tees and the River Birket, English rivers that have been impacted by pollution from industrial and/or domestic effluents for over 100 years. Male reproductive competitiveness was determined via competitive breeding experiments with males and females derived from a clean water site, employing DNA microsatellites to determine parentage outcome. Populations of stickleback collected from the three historically polluted sites showed evidence of recent population bottlenecks, although only the River Aire population showed low genetic diversity. In contrast, fish collected from two relatively unpolluted sites within the River Gowy and Houghton Springs showed weak, or no evidence of such bottlenecks. Nevertheless, males derived from polluted sites were able to reproduce successfully in competition with males derived from clean water exposures, indicating that these bottlenecks have not resulted in any substantial loss of reproductive fitness in males.

Compound- and mixture-specific differences in resistance to polycyclic aromatic hydrocarbons and PCB-126 among Fundulus heteroclitus subpopulations throughout the Elizabeth River estuary (Virginia, USA)

Atlantic killifish (Fundulus heteroclitus) inhabiting the Atlantic Wood Industries Superfund Site (Elizabeth River, Portsmouth, VA, USA) are resistant to the acute toxicity and cardiac teratogenesis caused by high levels of polycyclic aromatic hydrocarbons (PAHs) from creosote. The resistance is linked to down regulation of the aryl hydrocarbon receptor (AHR) pathway. We investigated the association between CYP1 activity, as a marker of potential AHR pathway suppression, and contaminant resistance in killifish subpopulations from sites throughout the estuary that varied significantly in PAH contamination level. Adult killifish and sediments were collected from seven sites across approximately 13.7 km in river length within the estuary and from a nearby reference site. Sediment PAH levels were determined using gas chromatography mass spectrometry. Embryos obtained via manual spawning were exposed to individual AHR agonists and PAH mixtures 24 h post fertilization (hpf); CYP1 activity was determined by in ovo ethoxyresorufin-o-deethylase (EROD) at 96 hpf, and cardiac deformity severity was scored at 144 hpf. The total PAH levels measured among the sites varied from approximately 200 to 125,000 ng/g dry sediment. Overall, the resistance to teratogenesis was strongest in the subpopulations from sites in or closest to the major PAH contamination sites, but even embryos from less-contaminated sites within the Elizabeth River demonstrated at least partial resistance to many challenges. Surprisingly, all of the subpopulations tested were highly resistant to PCB-126 (3,3',4,4',5-pentachlorobiphenyl). However, the degree of CYP1 activity response varied significantly among subpopulations and did not always correlate strongly with resistance to teratogenesis; some subpopulations resisted the cardiac teratogenesis caused by the challenges at doses that still elicited strong EROD induction. Our results suggest that there is variation in the adaptive phenotype exhibited by laboratory-spawned embryos from killifish subpopulations throughout the estuary. Furthermore, the results show that contaminants have affected killifish subpopulations throughout the estuary, even in sites with lower levels of PAHs.

Induction of gene responses in St. Lawrence River northern pike (Esox lucius) environmentally exposed to perfluorinated compounds

Municipal waste water effluents (MWWEs) are important sources of chemical contamination for aquatic environments. This study investigated the presence and effects of perfluorinated compounds (PFCs) in environmentally exposed northern pike (Esox lucius) collected upstream and downstream a major municipal waste water treatment plant (WWTP) in the St. Lawrence River, Canada. Twelve PFCs, including the newly detected perfluoroethylcyclohexane sulfonate (PFECHS), were quantified in fish muscle, liver, and plasma. Additionally, the expression of eight genes and the activity of three biomarkers were analyzed in fish tissues at both sites. Mean ∑PFC concentration in fish plasma collected upstream the WWTP was 185ng/g w.w. compared to 545ng/g w.w. downstream the point of release. PFECHS was quantified for the first time in St. Lawrence River fish (mean plasma concentration in MWWE fish: 5.07±4.72ng/g w.w.). Results of transcriptomic responses were tissue-specific and indicated significant up-regulation for metallothionein (MT) in blood and MT, glutathion-S-transferase (GST), superoxide dismutase (SOD), and cytochromes P450 1A1 (CYP1A1) in gill tissue of fish collected in the MWWE suggesting greater stress responses for organisms at this location. Significant relationships were found between vitellogenin (Vtg) gene expression in liver, Vtg activity in plasma and perfluorotridecanoic acid (PFTrA), perfluorotetradecanoic acid (PFTeA), and perfluorodecane sulfonate (PFDS) plasma concentrations. The possible endocrine effects of these PFCs should be further investigated.

Differential tolerance to copper, but no evidence of population-level genetic differences in a widely-dispersing native barnacle

Despite many estuaries having high levels of metal pollution, species are found to persist in these stressful environments. Populations of estuarine invertebrates exposed to toxic concentrations of such metals may be under selection. However, in species with a wide-dispersal potential, any short-term results of localized selection may be counteracted by external recruitment from populations not under selection. The barnacle Amphibalanus variegatus is found in nearshore coastal environments as well as sheltered embayments and estuaries, including metal-impacted estuaries, from New South Wales, Australia to Western Australia. The fertilised eggs of A. variegatus are brooded internally and released as larvae (nauplii), which remain in the water-column for ~14 days before settling. Hence the species has a considerable dispersal capacity. The purpose of this study was to examine whether populations of A. variegatus from metal-impacted sites, displayed a greater tolerance to a toxicant (copper) than reference populations. Adult barnacles where collected from twenty sites within two metal-impacted and fourteen sites within two reference estuaries. Within 24 h, adults were induced to spawn and the offspring were exposed to copper in a laboratory assay. Larvae collected from the metal-impacted estuaries demonstrated a greater tolerance to copper compared to those from reference sites. To determine if selection/localised in the metal impacted sites was occurring, the genetic structure of populations at three sites was examined using an AFLP methodology. No evidence of unique population identity and or selection (outlier loci) was detected suggesting that: (1) the tolerance displayed in the assay was derived from acclimation during development; and/or (2) that the populations are open preventing the fixation of any unique alleles.

Differential sensitivity to pro-oxidant exposure in two populations of killifish (Fundulus heteroclitus)

New Bedford Harbor (MA, U.S.A.; NBH) is a Superfund site inhabited by Atlantic killifish (Fundulus heteroclitus) with altered aryl hydrocarbon receptor (Ahr) signaling, leading to resistance to effects of polychlorinated biphenyls (PCBs) and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). The Ahr is a transcription factor that regulates gene expression of many Phase I and II detoxifying enzymes and interacts with Nrf2, a transcription factor that regulates the response to oxidative stress. This study tested the hypothesis that PCB-resistant killifish exhibit altered sensitivity to oxidative stress. Killifish F(1) embryos from NBH and a clean reference site (Scorton Creek, MA, U.S.A.; SC) were exposed to model pro-oxidant and Nrf2-activator, tert-butylhydroquinone (tBHQ). Embryos were exposed at specific embryonic developmental stages (5, 7, and 9 days post fertilization) and toxicity was assessed, using a deformity score, survival, heart rate, and gene expression to compare sensitivity between PCB -resistant and -sensitive (reference) populations. Acute exposure to tBHQ resulted in transient reduction in heart rate in NBH and SC F(1) embryos. However, embryos from NBH were more sensitive to tBHQ, with more frequent and severe deformities, including pericardial edema, tail deformities, small body size, and reduced pigment and erythrocytes. NBH embryos had lower basal expression of antioxidant genes catalase and glutathione-S-transferase alpha (gsta), and upon exposure to tBHQ, exhibited lower levels of expression of catalase, gsta, and superoxide dismutase compared to controls. This result suggests that adaptation to tolerate PCBs has altered the sensitivity of NBH fish to oxidative stress during embryonic development, demonstrating a cost of the PCB resistance adaptation.

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.