Genome-wide scan reveals signatures of selection related to pollution adaptation in non-model estuarine Atlantic killifish (Fundulus heteroclitus)

Abstract

In many human-altered ecosystems, organisms are increasingly faced with more diverse and complex environmental stressors and pollutant mixtures, to which the adaptations necessary to survive exposure are likely to be numerous and varied. Improving our understanding of the molecular mechanisms that underlie complex polygenic adaptations in natural settings requires significant toxicological, biochemical, physiological, and genomic data rarely available for non-model organisms. Here, we build upon two decades of study of adaptation to anthropogenic pollutants in a population of Atlantic killifish (Fundulus heteroclitus) that inhabits the creosote-contaminated Atlantic Wood Industries Superfund (AW) site on the Elizabeth River, Virginia in the United States. To better understand the genotypes that underlie previously characterized resistance to PCBs and PAHs, we performed Restriction site-Associated DNA sequencing (RADseq) on killifish from AW and two relatively clean reference sites (King’s Creek–KC, and Mains Creek–MC). Across the genome, we analyzed over 83,000 loci and 12,000 single nucleotide polymorphisms (SNPs). Shared across both comparisons of killifish from polluted (AW) and relatively unpolluted (KC and MC) sites, we found eight genomic regions with smoothed FST values significantly (p < 0.001) elevated above background. Using the recently published F. heteroclitus reference genome, we identified candidate genes in these significant regions involved in the AHR pathway (e.g. AIP, ARNT1c), as well as genes relating to cardiac structure and function. These genes represent both previously characterized and potentially novel molecular adaptations involved with various aspects of resistance to these environmental toxins.