Common mechanism underlies repeated evolution of extreme pollution tolerance

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.

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

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

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.

Functional genomics to assess biological responses to marine pollution at physiological and evolutionary timescales: toward a vision of predictive ecotoxicology

Marine pollution is ubiquitous, and is one of the key factors influencing contemporary marine biodiversity worldwide. To protect marine biodiversity, how do we surveil, document and predict the short- and long-term impacts of pollutants on at-risk species? Modern genomics tools offer high-throughput, information-rich and increasingly cost-effective approaches for characterizing biological responses to environmental stress, and are important tools within an increasing sophisticated kit for surveiling and assessing impacts of pollutants on marine species. Through the lens of recent research in marine killifish, we illustrate how genomics tools may be useful for screening chemicals and pollutants for biological activity and to reveal specific mechanisms of action. The high dimensionality of transcriptomic responses enables their usage as highly specific fingerprints of exposure, and these fingerprints can be used to diagnose environmental problems. We also emphasize that molecular pathways recruited to respond at physiological timescales are the same pathways that may be targets for natural selection during chronic exposure to pollutants. Gene complement and sequence variation in those pathways can be related to variation in sensitivity to environmental pollutants within and among species. Furthermore, allelic variation associated with evolved tolerance in those pathways could be tracked to estimate the pace of environmental health decline and recovery. We finish by integrating these paradigms into a vision of how genomics approaches could anchor a modernized framework for advancing the predictive capacity of environmental and ecotoxicological science.

Early life exposure to PCB126 results in delayed mortality and growth impairment in the zebrafish larvae

The occurrence of chronic or delayed toxicity resulting from the exposure to sublethal chemical concentrations is an increasing concern in environmental risk assessment. The Fish Embryo Toxicity (FET) test with zebrafish provides a reliable prediction of acute toxicity in adult fish, but it cannot yet be applied to predict the occurrence of chronic or delayed toxicity. Identification of sublethal FET endpoints that can assist in predicting the occurrence of chronic or delayed toxicity would be advantageous. The present study characterized the occurrence of delayed toxicity in zebrafish larvae following early exposure to PCB126, previously described to cause delayed effects in the common sole. The first aim was to investigate the occurrence and temporal profiles of delayed toxicity during zebrafish larval development and compare them to those previously described for sole to evaluate the suitability of zebrafish as a model fish species for delayed toxicity assessment. The second aim was to examine the correlation between the sublethal endpoints assessed during embryonal and early larval development and the delayed effects observed during later larval development. After exposure to PCB126 (3-3000ng/L) until 5 days post fertilization (dpf), larvae were reared in clean water until 14 or 28 dpf. Mortality and sublethal morphological and behavioural endpoints were recorded daily, and growth was assessed at 28 dpf. Early life exposure to PCB126 caused delayed mortality (300 ng/L and 3000 ng/L) as well as growth impairment and delayed development (100 ng/L) during the clean water period. Effects on swim bladder inflation and cartilaginous tissues within 5 dpf were the most promising for predicting delayed mortality and sublethal effects, such as decreased standard length, delayed metamorphosis, reduced inflation of swim bladder and column malformations. The EC50 value for swim bladder inflation at 5 dpf (169 ng/L) was similar to the LC50 value at 8 dpf (188 and 202 ng/L in two experiments). Interestingly, the patterns of delayed mortality and delayed effects on growth and development were similar between sole and zebrafish. This indicates the comparability of critical developmental stages across divergent fish species such as a cold water marine flatfish and a tropical freshwater cyprinid. Additionally, sublethal effects in early embryo-larval stages were found promising for predicting delayed lethal and sublethal effects of PCB126. Therefore, the proposed method with zebrafish is expected to provide valuable information on delayed mortality and delayed sublethal effects of chemicals and environmental samples that may be extrapolated to other species.

Ecotoxicology of polychlorinated biphenyls in fish--a critical review

Polychlorinated biphenyls (PCBs) are widespread persistent anthropogenic contaminants that can accumulate in tissues of fish. The toxicity of PCBs and their transformation products has been investigated for nearly 50 years, but there is a lack of consensus regarding the effects of these environmental contaminants on wild fish populations. The objective of this review is to critically examine these investigations and evaluate publicly available databases for evidence of effects of PCBs in wild fish. Biological activity of PCBs is limited to a small proportion of PCB congeners [e.g., dioxin-like PCBs (DL-PCBs)] and occurs at concentrations that are typically orders of magnitude higher than PCB levels detected in wild fish. Induction of biomarkers consistent with PCB exposure (e.g., induction of cytochrome P450 monooxygenase system) has been evaluated frequently and shown to be induced in fish from some environments, but there does not appear to be consistent reports of damage (i.e., biomarkers of effect) to biomolecules (i.e., oxidative injury) in these fish. Numerous investigations of endocrine system dysfunction or effects on other organ systems have been conducted in wild fish, but collectively there is no consistent evidence of PCB effects on these systems in wild fish. Early life stage toxicity of DL-PCBs does not appear to occur at concentrations reported in wild fish embryos, and results do not support an association between PCBs and decreased survival of early life stages of wild fish. Overall, there appears to be little evidence that PCBs have had any widespread effect on the health or survival of wild fish.

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.

Ocean acidification research in the 'post-genomic' era: Roadmaps from the purple sea urchin Strongylocentrotus purpuratus

Advances in nucleic acid sequencing technology are removing obstacles that historically prevented use of genomics within ocean change biology. As one of the first marine calcifiers to have its genome sequenced, purple sea urchins (Strongylocentrotus purpuratus) have been the subject of early research exploring genomic responses to ocean acidification, work that points to future experiments and illustrates the value of expanding genomic resources to other marine organisms in this new 'post-genomic' era. This review presents case studies of S. purpuratus demonstrating the ability of genomic experiments to address major knowledge gaps within ocean acidification. Ocean acidification research has focused largely on species vulnerability, and studies exploring mechanistic bases of tolerance toward low pH seawater are comparatively few. Transcriptomic responses to high pCO₂ seawater in a population of urchins already encountering low pH conditions have cast light on traits required for success in future oceans. Secondly, there is relatively little information on whether marine organisms possess the capacity to adapt to oceans progressively decreasing in pH. Genomics offers powerful methods to investigate evolutionary responses to ocean acidification and recent work in S. purpuratus has identified genes under selection in acidified seawater. Finally, relatively few ocean acidification experiments investigate how shifts in seawater pH combine with other environmental factors to influence organism performance. In S. purpuratus, transcriptomics has provided insight into physiological responses of urchins exposed simultaneously to warmer and more acidic seawater. Collectively, these data support that similar breakthroughs will occur as genomic resources are developed for other marine species.

Paracellular pathway remodeling enhances sodium secretion by teleost fish in hypersaline environments

In vertebrate salt-secreting epithelia, Na(+) moves passively down an electrochemical gradient via a paracellular pathway. We assessed how this pathway is modified to allow Na(+) secretion in hypersaline environments. Mummichogs (Fundulus heteroclitus) acclimated to hypersaline [2× seawater (2SW), 64‰] for 30 days developed invasive projections of accessory cells with an increased area of tight junctions, detected by punctate distribution of CFTR (cystic fibrosis transmembrane conductance regulator) immunofluorescence and transmission electron miscroscopy of the opercular epithelia, which form a gill-like tissue rich in ionocytes. Distribution of CFTR was not explained by membrane raft organization, because chlorpromazine (50 μmol l(-1)) and filipin (1.5 μmol l(-1)) did not affect opercular epithelia electrophysiology. Isolated opercular epithelia bathed in SW on the mucosal side had a transepithelial potential (Vt) of +40.1±0.9 mV (N=24), sufficient for passive Na(+) secretion (Nernst equilibrium voltage≡ENa=+24.11 mV). Opercular epithelia from fish acclimated to 2SW and bathed in 2SW had higher Vt of +45.1±1.2 mV (N=24), sufficient for passive Na(+) secretion (ENa=+40.74 mV), but with diminished net driving force. Bumetanide block of Cl(-) secretion reduced Vt by 45% and 29% in SW and 2SW, respectively, a decrease in the driving force for Na(+) extrusion. Estimates of shunt conductance from epithelial conductance (Gt) versus short-circuit current (Isc) plots (extrapolation to zero Isc) suggested a reduction in total epithelial shunt conductance in 2SW-acclimated fish. In contrast, the morphological elaboration of tight junctions, leading to an increase in accessory-cell-ionocyte contact points, suggests an increase in local paracellular conductance, compensating for the diminished net driving force for Na(+) and allowing salt secretion, even in extreme salinities.

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.

Genomic and genotoxic responses to controlled weathered-oil exposures confirm and extend field studies on impacts of the Deepwater Horizon oil spill on native killifish

To understand the ecotoxicological impacts of the Deepwater Horizon oil spill, field studies provide a context for ecological realism but laboratory-based studies offer power for connecting biological effects with specific causes. As a complement to field studies, we characterized genome-wide gene expression responses of Gulf killifish (Fundulus grandis) to oil-contaminated waters in controlled laboratory exposures. Transcriptional responses to the highest concentrations of oiled water in the laboratory were predictive of field-observed responses that coincided with the timing and location of major oiling. The transcriptional response to the low concentration (∼10-fold lower than the high concentration) was distinct from the high concentration and was not predictive of major oiling in the field. The high concentration response was characterized by activation of the molecular signaling pathway that facilitates oil metabolism and oil toxicity. The high concentration also induced DNA damage. The low concentration invoked expression of genes that may support a compensatory response, including genes associated with regulation of transcription, cell cycle progression, RNA processing, DNA damage, and apoptosis. We conclude that the gene expression response detected in the field was a robust indicator of exposure to the toxic components of contaminating oil, that animals in the field were exposed to relatively high concentrations that are especially damaging to early life stages, and that such exposures can damage DNA.

Chronic exposure of killifish to a highly polluted environment desensitizes estrogen-responsive reproductive and biomarker genes

Reproductive and endocrine disruption is commonly reported in aquatic species exposed to complex contaminant mixtures. We previously reported that Atlantic killifish (Fundulus heteroclitus) from the chronically contaminated Newark Bay, NJ, exhibit multiple endocrine disrupting effects, including inhibition of vitellogenesis (yolk protein synthesis) in females and false negative vitellogenin biomarker responses in males. Here, we characterized the effects on estrogen signaling and the transcriptional regulation of estrogen-responsive genes in this model population. First, a dose-response study tested the hypothesis that reproductive biomarkers (vtg1, vtg2, chg H, chg Hm, chg L) in Newark Bay killifish are relatively less sensitive to 17β-estradiol at the transcriptional level, relative to a reference (Tuckerton, NJ) population. The second study assessed expression for various metabolism (cyp1a, cyp3a30, mdr) and estrogen receptor (ER α, ER βa, ER βb) genes under basal and estrogen treatment conditions in both populations. Hepatic metabolism of 17β-estradiol was also evaluated in vitro as an integrated endpoint for adverse effects on metabolism. In the third study, gene methylation was evaluated for promoters of vtg1 (8 CpGs) and vtg2 (10 CpGs) in both populations, and vtg1 promoter sequences were examined for single nucleotide polymorphism (SNPs). Overall, these studies show that multi-chemical exposures at Newark Bay have desensitized all reproductive biomarkers tested to estrogen. For example, at 10ng/g 17β-estradiol, inhibition of gene induction ranged from 62% to 97% for all genes tested in the Newark Bay population, relative to induction levels in the reference population. The basis for this recalcitrant phenotype could not be explained by a change in 17β-estradiol metabolism, nuclear estrogen receptor expression, promoter methylation (gene silencing) or SNPs, all of which were unaltered and normal in the Newark Bay population. The decreased transcriptional sensitivity of estrogen-responsive genes is suggestive of a broad effect on estrogen receptor pathway signaling, and provides insight into the mechanisms of the endocrine disrupting effects in the Newark Bay population.

Evolved resistance to PCB- and PAH-induced cardiac teratogenesis, and reduced CYP1A activity in Gulf killifish (Fundulus grandis) populations from the Houston Ship Channel, Texas

The Houston Ship Channel (HSC), connecting Houston, Texas to Galveston Bay and ultimately the Gulf of Mexico, is heavily industrialized and includes several areas that have historically been identified as containing significant levels of mercury, dioxins, furans, polychlorinated biphenyls (PCBs), and polycyclic aromatic hydrocarbons (PAHs). Gulf killifish, Fundulus grandis, inhabit this entire estuarine system, including the most contaminated areas. F. grandis is the sister species of the well-established estuarine model organism Fundulus heteroclitus, for which heritable resistance to both PCB and PAH toxicity has been documented in several populations. F. grandis collected from two Superfund sites on the HSC and from a reference population were used to establish breeding colonies. F1 embryos from HSC populations were approximately 1000-fold more resistant to PCB126- and 2-5-fold more resistant to coal tar-induced cardiovascular teratogenesis, relative to embryos from the reference population. Reciprocal crosses between reference and contaminated populations exhibit an intermediate level of resistance, confirming that observed protection is genetic and biparentally inherited. Ethoxyresorufin-O-deethylase (EROD) data confirm a reduction in basal and induced cytochrome P4501A (CYP1A) activity in resistant populations of F. grandis. This result is consistent with responses previously described for resistant populations of F. heteroclitus, specifically a recalcitrant aryl hydrocarbon receptor (AHR) pathway. The decreased levels of cardiovascular teratogenesis, and decrease in CYP1A inducibility in response to PCB126 and a PAH mixture, suggest that HSC F. grandis populations have adapted to chronic contaminants exposures via a mechanism similar to that previously described for F. heteroclitus. To the best of our knowledge, this is the first documentation of evolved pollution resistance in F. grandis. Additionally, the mechanistic similarities between the population adaptation observed in this study and previous work in F. heteroclitus suggest that genetic variation predating the evolutionary divergence of these two species may best explain the apparent rapid parallel evolution of pollution resistance in genetically and geographically distinct species and populations.

Rapid parallel evolution of standing variation in a single, complex, genomic region is associated with life history in steelhead/rainbow trout

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. Such changes may be population-specific or, alternatively, may involve parallel evolution of the same genomic region in multiple populations, if that region contains genes or co-adapted gene complexes affecting the selected trait(s). Both quantitative and population genetic approaches have identified associations between specific genomic regions and the anadromous (steelhead) and resident (rainbow trout) life-history strategies of Oncorhynchus mykiss. Here, we use genotype data from 95 single nucleotide polymorphisms and show that the distribution of variation in a large region of one chromosome, Omy5, is strongly associated with life-history differentiation in multiple above-barrier populations of rainbow trout and their anadromous steelhead ancestors. The associated loci are in strong linkage disequilibrium, suggesting the presence of a chromosomal inversion or other rearrangement limiting recombination. These results provide the first evidence of a common genomic basis for life-history variation in O. mykiss in a geographically diverse set of populations and extend our knowledge of the heritable basis of rapid adaptation of complex traits in novel habitats.

Common functional targets of adaptive micro- and macro-evolutionary divergence in killifish

Environmental salinity presents a key barrier to dispersal for most aquatic organisms, and adaptation to alternate osmotic environments likely enables species diversification. Little is known of the functional basis for derived tolerance to environmental salinity. We integrate comparative physiology and functional genomics to explore the mechanistic underpinnings of evolved variation in osmotic plasticity within and among two species of killifish; Fundulus majalis harbours the ancestral mainly salt-tolerant phenotype, whereas Fundulus heteroclitus harbours a derived physiology that retains extreme salt tolerance but with expanded osmotic plasticity towards the freshwater end of the osmotic continuum. Common-garden comparative hypo-osmotic challenge experiments show that F. heteroclitus is capable of remodelling gill epithelia more quickly and at more extreme osmotic challenge than F. majalis. We detect an unusual pattern of baseline transcriptome divergence, where neutral evolutionary processes appear to govern expression divergence within species, but patterns of divergence for these genes between species do not follow neutral expectations. During acclimation, genome expression profiling identifies mechanisms of acclimation-associated response that are conserved within the genus including regulation of paracellular permeability. In contrast, several responses vary among species including those putatively associated with cell volume regulation, and these same mechanisms are targets for adaptive physiological divergence along osmotic gradients within F. heteroclitus. As such, the genomic and physiological mechanisms that are associated with adaptive fine-tuning within species also contribute to macro-evolutionary divergence as species diversify across osmotic niches.

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.

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

Background

The most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed.

Results

Since the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs.

Conclusions

The observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.

Spatial variability of metal bioaccumulation in estuarine killifish (Fundulus heteroclitus) at the Callahan mine superfund site, Brooksville, ME

The former Callahan Mine Site in Brooksville, ME, is an open-pit, hardrock mine site in an intertidal system, thus providing a unique opportunity to evaluate how metal-enriched sediments and overlying water impact estuarine food webs. Copper, zinc, cadmium, and lead concentrations in sediment, whole water, and Atlantic killifish (Fundulus heteroclitus) were evaluated at sites in Goose Pond (GP; Callahan Mine Site) and at reference sites. The metal concentrations of sediment, water, and fish were spatially distinct and significantly greater at the mine site than in the reference estuary. Sediment concentrations were particularly elevated and were above probable effects levels for all four metals adjacent to the tailings pile. Even in this well-mixed system, water metal concentrations were significantly elevated adjacent to the tailings pile, and concentrations of Cu and Zn were above ambient water-quality criteria for chronic marine exposure. Neither organic matter in the sediment nor salinity or pH of the water explained the metal concentrations. Adjacent to the tailings pile, killifish metal body burdens were elevated and were significantly related to both sediment and aqueous concentrations. In conclusion, (1) the contaminated sediment and seepage from the tailings impoundment and waste rock pile no. 3 create a continual flux of metals into the water column, (2) the metals are bioavailable and bioconcentrating as evident in the killifish tissue concentrations, and (3) Callahan Mine is directly affecting metal bioaccumulation in fauna residing in the GP estuary and, potentially, in Penobscot Bay by the way of “trophic nekton relay.”

Cancer: a missing link in ecosystem functioning?

Cancer is a disease that affects the majority of metazoan species and, before directly causing host death, is likely to influence the competitive abilities of individuals, their susceptibility to pathogens, their vulnerability to predators, and their ability to disperse. Despite the potential importance of these ecological impacts, cancer is rarely incorporated into model ecosystems. We describe here the diversity of ways in which oncogenic phenomena, from precancerous lesions to generalized metastatic cancers, may affect ecological processes that govern biotic interactions. We argue that oncogenic phenomena, despite their complexity, can have significant and sometimes predictable ecological consequences. Our aim is to provide a new perspective on the ecological and evolutionary significance of cancer in wildlife, and to stimulate research on this topic.

Assessing gene network stability and individual variability in the fathead minnow (Pimephales promelas) transcriptome

Transcriptomics is increasingly used to assess biological responses to environmental stimuli and stressors such as aquatic pollutants. However, fundamental studies characterizing individual variability in mRNA levels are lacking, which currently limits the use of transcriptomics in environmental monitoring assessments. To address individual variability in transcript abundance, we performed a meta-analysis on 231 microarrays that were conducted in the fathead minnow (FHM), a widely used toxicological model. The mean variability for gene probes was ranked from most to least variable based upon the coefficient of variation. Transcripts that were the most variable in individual tissues included NADH dehydrogenase flavoprotein 1, GTPase IMAP family member 7-like and v-set domain-containing T-cell activation inhibitor 1-like while genes encoding ribosomal proteins (rpl24 and rpl36), basic transcription factor 3, and nascent polypeptide-associated complex alpha subunit were the least variable in individuals across a range of microarray experiments. Gene networks that showed high variability (based upon the variation in expression of individual members within the network) included cell proliferation, metabolism (steroid, lipids, and glucose), cell adhesion, vascularization, and regeneration while those that showed low variability (more stability) included mRNA and rRNA processing, regulation of translational fidelity, RNA splicing, and ribosome biogenesis. Real-time PCR was conducted on a subset of genes for comparison of variability collected from the microarrays. There was a significant positive relationship between the two methods when measuring individual variability, suggesting that variability detected in microarray data can be used to guide decisions on sample sizes for measuring transcripts in real-time PCR experiments. A power analysis revealed that measuring estrogen receptor ba (esrba) requires fewer biological replicates than that of estrogen receptor bb (esrbb) in the gonad and samples sizes required to detect a 50% change for reproductive-related transcripts is between 12 and 20. Characterizing individual variability at the molecular level will prove necessary as efforts are made toward integrating molecular tools into environmental risk assessments.

Multitissue molecular, genomic, and developmental effects of the Deepwater Horizon oil spill on resident Gulf killifish (Fundulus grandis)

The Deepwater Horizon oil rig disaster resulted in crude oil contamination along the Gulf coast in sensitive estuaries. Toxicity from exposure to crude oil can affect populations of fish that live or breed in oiled habitats as seen following the Exxon Valdez oil spill. In an ongoing study of the effects of Deepwater Horizon crude oil on fish, Gulf killifish ( Fundulus grandis ) were collected from an oiled site (Grande Terre, LA) and two reference locations (coastal MS and AL) and monitored for measures of exposure to crude oil. Killifish collected from Grande Terre had divergent gene expression in the liver and gill tissue coincident with the arrival of contaminating oil and up-regulation of cytochrome P4501A (CYP1A) protein in gill, liver, intestine, and head kidney for over one year following peak landfall of oil (August 2011) compared to fish collected from reference sites. Furthermore, laboratory exposures of Gulf killifish embryos to field-collected sediments from Grande Terre and Barataria Bay, LA, also resulted in increased CYP1A and developmental abnormalities when exposed to sediments collected from oiled sites compared to exposure to sediments collected from a reference site. These data are predictive of population-level impacts in fish exposed to sediments from oiled locations along the Gulf of Mexico coast.

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.

Reproductive and developmental toxicity of dioxin in fish

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD or dioxin) is a global environmental contaminant and the prototypical ligand for investigating aryl hydrocarbon receptor (AHR)-mediated toxicity. Environmental exposure to TCDD results in developmental and reproductive toxicity in fish, birds and mammals. To resolve the ecotoxicological relevance and human health risks posed by exposure to dioxin-like AHR agonists, a vertebrate model is needed that allows for toxicity studies at various levels of biological organization, assesses adverse reproductive and developmental effects and establishes appropriate integrative correlations between different levels of effects. Here we describe the reproductive and developmental toxicity of TCDD in feral fish species and summarize how using the zebrafish model to investigate TCDD toxicity has enabled us to characterize the AHR signaling in fish and to better understand how dioxin-like chemicals induce toxicity. We propose that such studies can be used to predict the risks that AHR ligands pose to feral fish populations and provide a platform for integrating risk assessments for both ecologically relevant organisms and humans.

Genomic and physiological footprint of the Deepwater Horizon oil spill on resident marsh fishes

The biological consequences of the Deepwater Horizon oil spill are unknown, especially for resident organisms. Here, we report results from a field study tracking the effects of contaminating oil across space and time in resident killifish during the first 4 mo of the spill event. Remote sensing and analytical chemistry identified exposures, which were linked to effects in fish characterized by genome expression and associated gill immunohistochemistry, despite very low concentrations of hydrocarbons remaining in water and tissues. Divergence in genome expression coincides with contaminating oil and is consistent with genome responses that are predictive of exposure to hydrocarbon-like chemicals and indicative of physiological and reproductive impairment. Oil-contaminated waters are also associated with aberrant protein expression in gill tissues of larval and adult fish. These data suggest that heavily weathered crude oil from the spill imparts significant biological impacts in sensitive Louisiana marshes, some of which remain for over 2 mo following initial exposures.