Evolution of tolerance to PCBs and susceptibility to a bacterial pathogen (\u3cem\u3eVibrio harveyi\u3c/em\u3e) in Atlantic killifish (\u3cem\u3eFundulus heteroclitus\u3c/em\u3e) from New Bedford (MA, USA) harbor

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

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

Author Posting. © The Author(s), 2014. This is the author's version of the work. It is posted here by permission of Elsevier for personal use, not for redistribution. The definitive version was published in Aquatic Toxicology 159 (2015): 156-166, doi:10.1016/j.aquatox.2014.12.017.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), Ca2+ 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 Ca2+ 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 Ca2+ signaling proteins or potential signaling partners in NBH killifish.Funding was provided through the NIEHS Superfund Research Program UC Davis (INP and EBF; P42-ES004699) and Boston University (JJS and JVG; P42-ES007381). Support was supplied via the UC Davis NHLBI Training Grant (T32-HL086350, EBF). Additional support came from NIEHS 1R01-ES014901, 1R01-ES017425, the UC Davis Center for Children’s Environmental Health (1P01-ES011269, U.S. Environmental Protection Agency Grant 8354320), and an unrestricted JB Johnson Foundation gift grant.2015-12-1

Population genetic diversity and fitness in multiple environments

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

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

Author Posting. © The Author(s), 2014. This is the author's version of the work. It is posted here by permission of Elsevier for personal use, not for redistribution. The definitive version was published in Aquatic Toxicology 159 (2015): 198-207, doi:10.1016/j.aquatox.2014.12.010.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 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 AhR2 or CYP1A in liver and gills of NBH fish. In NBH fish, but not in SC fish, there was increased expression of hepatic PXR, CYP3A and Pgp genes 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 expression of these genes between the two populations. In SC fish, but not in NBH fish, there was increased 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.This study was supported by grants from FORMAS (216-2007-468) and University of Gothenburg to MCC, and by the Superfund Research Program at Boston University, NIH grant P42ES007381 to JJS, MEH, and SIK. Data interpretation was aided by reference to a preliminary draft of the Fundulus heteroclitus genome sequence, which was supported by funding from the National Science Foundation (collaborative research grants DEB-1120512, DEB-1265282, DEB-1120013, DEB-1120263, DEB-1120333, DEB-1120398). This study was also supported by NOAA Grant No. NA16RG2273 (WHOI Sea Grant Project No. R/P-70 to SIK and MEH) and by funding from Adlerbertska Forskningsstiftelsen, Helge Ax:son Johnsons Stiftelse and Wilhelm och Martina Lundgrens Vetenskapsfond to BW and JG

Proestou_etal_Killifish_SNP_sequences

Fasta format file containing the flanking sequences for validated single nucleotide polymorphisms (SNPs) included in the population genetic analysis

Data from: 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

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

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

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

Background: The non-migratory killifish Fundulus heteroclitus inhabits clean and polluted environments interspersed throughout its range along the Atlantic coast of North America. Several populations of this species have successfully adapted to environments contaminated with toxic aromatic hydrocarbon pollutants such as polychlorinated biphenyls (PCBs). Previous studies suggest that the mechanism of resistance to these and other “dioxin-like compounds” (DLCs) may involve reduced signaling through the aryl hydrocarbon receptor (AHR) pathway. Here we investigated gene diversity and evidence for positive selection at three AHR-related loci (AHR1, AHR2, AHRR) in F. heteroclitus by comparing alleles from seven locations ranging over 600 km along the northeastern US, including extremely polluted and reference estuaries, with a focus on New Bedford Harbor (MA, USA), a PCB Superfund site, and nearby reference sites. Results: We identified 98 single nucleotide polymorphisms within three AHR-related loci among all populations, including synonymous and nonsynonymous substitutions. Haplotype distributions were spatially segregated and F-statistics suggested strong population genetic structure at these loci, consistent with previous studies showing strong population genetic structure at other F. heteroclitus loci. Genetic diversity at these three loci was not significantly different in contaminated sites as compared to reference sites. However, for AHR2 the New Bedford Harbor population had significant FST values in comparison to the nearest reference populations. Tests for positive selection revealed ten nonsynonymous polymorphisms in AHR1 and four in AHR2. Four nonsynonymous SNPs in AHR1 and three in AHR2 showed large differences in base frequency between New Bedford Harbor and its reference site. Tests for isolation-by-distance revealed evidence for non-neutral change at the AHR2 locus. Conclusion: Together, these data suggest that F. heteroclitus populations in reference and polluted sites have similar genetic diversity, providing no evidence for strong genetic bottlenecks for populations in polluted locations. However, the data provide evidence for genetic differentiation among sites, selection at specific nucleotides in AHR1 and AHR2, and specific AHR2 SNPs and haplotypes that are associated with the PCB-resistant phenotype in the New Bedford Harbor population. The results suggest that AHRs, and especially AHR2, may be important, recurring targets for selection in local adaptation to dioxin-like aromatic hydrocarbon contaminants

AHRRhaplotypes_final

Sequences of 65 AHRR haplotypes. See published paper for details