Population genomics of North American northern pike: variation and sex-specific signals from a chromosome-level, long read genome assembly

We present a chromosome-level, long-read genome assembly as a reference for northern pike (Esox lucius) where 97.5% of the genome is chromosome-anchored and N50 falls at 37.5 Mb. Whole-genome resequencing was genotyped using this assembly for 47 northern pike representing six North American populations from Alaska to New Jersey. We discovered that a disproportionate frequency of genetic polymorphism exists among populations east and west of the North American Continental Divide (NACD), indicating reproductive isolation across this barrier. Genome-wide analysis of heterozygous SNP density revealed a remarkable lack of genetic variation with 1 polymorphic site every 6.3kb in the Yukon River drainage and one every 16.5kb east of the NACD. Observed heterozygosity (Ho), nucleotide diversity (π), and Tajima’s D are depressed in populations east of the NACD (east vs. west: Ho: 0.092 vs 0.31; π: 0.092 vs 0.28; Tajima’s D: -1.61 vs -0.47). We confirm the presence of the master sex determining (MSD) gene, amhby, in the Yukon River drainage and in an invasive population in British Columbia and confirm its absence in populations east of the NACD. We also describe an Alaskan population where amhby is present but not associated with male gender determination. Our results support that northern pike originally colonized North America through Beringia, that Alaska provided an unglaciated refugium for northern pike during the last ice age, and southeast of the NACD was colonized by a small founding population(s) that lost amhby.