Ryanodine receptor and FK506 binding protein 1 in the Atlantic killifish (Fundulus heteroclitus) : a phylogenetic and population-based comparison

Author Posting. © The Author(s), 2017. This is the author's version of the work. It is posted here under a nonexclusive, irrevocable, paid-up, worldwide license granted to WHOI. It is made available for personal use, not for redistribution. The definitive version was published in Aquatic Toxicology 192 (2017): 105-115, doi:10.1016/j.aquatox.2017.09.002.Non-dioxin-like polychlorinated biphenyls (NDL PCBs) activate ryanodine receptors (RyR), microsomal Ca2+ channels of broad significance. Teleost fish may be important models for NDL PCB neurotoxicity, and we used sequencing databases to characterize teleost RyR and FK506 binding protein 12 or 12.6 kDa (genes FKBP1A; FKBP1B), which promote NDL PCB-triggered Ca2+ dysregulation. Particular focus was placed on describing genes in the Atlantic killifish (Fundulus heteroclitus) genome and searching available RNA-sequencing datasets for single nucleotide variants (SNV) between PCB tolerant killifish from New Bedford Harbor (NBH) versus sensitive killifish from Scorton Creek (SC), MA. Consistent with the teleost whole genome duplication (tWGD), killifish have six RyR genes, corresponding to a and b paralogs of mammalian RyR1, 2 and 3. The presence of six RyR genes was consistent in all teleosts investigated including zebrafish. Killifish have four FKBP1; one FKBP1b and three FKBP1a named FKBP1aa, FKBP1ab, likely from the tWGD and a single gene duplicate FKBP1a3 suggested to have arisen in Atherinomorphae. The RyR and FKBP1 genes displayed tissue and developmental stage-specific mRNA expression, and the previously uncharacterized RyR3, herein named RyR3b, and all FKBP1 genes were prominent in brain. We identified a SNV in RyR3b encoding missense mutation E1458D. In NBH killifish, 57% were heterozygous and 28% were homozygous for this SNV, whereas almost all SC killifish (94%) lacked the variant (n≥39 per population). The outlined sequence differences between mammalian and teleost RyR and FKBP1 together with outlined population differences in SNV frequency may contribute to our understanding of NDL PCB neurotoxicity.This research was supported by the KC Donnelly Research Externship made possible by the National Institute of Environmental Health Sciences’ Superfund Research Program (EBH) and the Superfund Research Programs at UC Davis (INP and EBH; P42ES004699) and Boston University (JJS, JVG, MEH, SIK; P42ES007381). Additional support was provided by the National Institute of Health (INP; R01 ES014901; and P01 AR052354) and by National Science Foundation collaborative research grants (MEH and SIK; DEB-1265282 and DEB-1120263). This research was also supported in part by an appointment (to BC) with the Postdoctoral Research Program at the U.S. Environmental Protection (US EPA) Office of Research and Development administered by the Oak Ridge Institute for Science and Education (ORISE) through Interagency Agreement No. DW92429801 between the U.S. Department of Energy and the US EPA