Genetic architecture drives seasonal onset of hibernation in the 13-lined ground squirrel

Hibernation is a highly dynamic phenotype whose timing, for many mammals, is controlled by a circannual clock and accompanied by rhythms in body mass and food intake. When housed in an animal facility, 13-lined ground squirrels exhibit individual variation in the seasonal onset of hibernation, which is not explained by environmental or biological factors, such as body mass and sex. We hypothesized that underlying genetic architecture instead drives variation in this timing. After first increasing the contiguity of the genome assembly, we therefore employed a genotype-by-sequencing approach to characterize genetic variation in 153 13-lined ground squirrels. Combining this with datalogger records, we estimated high heritability (61-100%) for the seasonal onset of hibernation. After applying a genome-wide scan with 46,996 variants, we also identified 21 loci significantly associated with hibernation immergence, which alone accounted for 54% of the variance in the phenotype. The most significant marker (SNP 15, p=3.81x10-6) was located near prolactin-releasing hormone receptor (PRLHR), a gene that regulates food intake and energy homeostasis. Other significant loci were located near genes functionally related to hibernation physiology, including muscarinic acetylcholine receptor M2 (CHRM2), involved in the control of heart rate, exocyst complex component 4 (EXOC4) and prohormone convertase 2 (PCSK2), both of which are involved in insulin signaling and processing. Finally, we applied an expression quantitative loci (eQTL) analysis using existing transcriptome datasets, and we identified significant (q<0.1) associations for 9/21 variants. Our results highlight the power of applying a genetic mapping strategy to hibernation and present new insight into the genetics driving its seasonal onset.