Identification of sex-specific SNPS in burbot Lota lota using RAD sequencing: conservation and management applications

The development of sex-specific genetic assays in a species provides both a method for identifying the system of sex determination and a valuable tool to address questions of conservation and management importance. In this study, we focused on the identification of single nucleotide polymorphisms (SNPs) that differentiate genetic sex in burbot Lota lota. Burbot are the only true freshwater representative of the cod family and a species of conservation and management importance throughout Eurasia and North America. To identify sex-specific SNPs, we utilized restriction site-associated DNA sequencing (RADseq) to interrogate thousands of SNPs in burbot samples of known phenotypic sex. We discovered 170,569 biallelic SNPs, none of which fit the pattern expected under female heterogamety. However, we identified 22 SNPs that fit the pattern expected under male heterogamety (males heterozygous XY, females fixed XX) and, from these, developed two genetic assays that robustly (~ 97% genotyping success) and accurately (> 99% correct) sexed burbot samples. These sex-specific genetic assays will benefit growing conservation aquaculture programs for this species and allow future assessments of sex-specific migration, growth, and mortality.