A de novo genome assembly and annotation of the southern flying squirrel (Glaucomys volans)

Abstract

Northern (Glaucomys sabrinus) and southern (Glaucomys volans) flying squirrels are widespread species distributed across much of North America. Northern flying squirrels are common inhabitants of the boreal forest, also occurring in coniferous forest remnants farther south, whereas the southern flying squirrel range is centered in eastern temperate woodlands. These two flying squirrel species exhibit a hybrid zone across a latitudinal gradient in an area of recent secondary contact. Glaucomys hybrid offspring are viable and can successfully backcross with either parental species, however, the fitness implications of such events are currently unknown. Some populations of G. sabrinus are endangered, and thus, interspecific hybridization is a key conservation concern in flying squirrels. We sequenced and assembled a de novo long-read genome from a G. volans individual sampled in southern Ontario, Canada, while four short-read genomes (2 G. sabrinus and 2 G. volans, all from Ontario) were re-sequenced on Illumina platforms. The final genome assembly consisted of approximately 2.40Gb with a scaffold N50 of 455.26Kb. Benchmarking Universal Single-Copy Orthologs reconstructed 3,742 (91.2%) complete mammalian genes and genome annotation using RNA-seq identified the locations of 19,124 protein-coding genes. The four short-read individuals were aligned to our reference genome to investigate the demographic history of the two species. A Principal Component Analysis clearly separated re-sequenced individuals, while inferring population size history using the Pairwise Sequentially Markovian Coalescent model noted an approximate species split one million years ago, and a single, possibly recently introgressed individual.