A chromosome-level genome assembly for the Eastern Fence Lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology

Abstract

High-quality genomic resources facilitate population-level and species-level comparisons to answer questions about behavioral ecology, morphological and physiological adaptations, as well as the evolution of genomic architecture. Squamate reptiles (lizards and snakes) are particularly diverse in characteristics that have intrigued evolutionary biologists, but high-quality genomic resources for squamates are relatively sparse. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources. We present a high-quality chromosome-level reference genome assembly, SceUnd1.0, (utilizing 10X Genomics Chromium, HiC, and PacBio data) and tissue/developmental stage transcriptomes for the Eastern Fence Lizard, Sceloporus undulatus. We performed synteny analysis with other available squamate chromosome-level assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes in S. undulatus. Using this new S. undulatus genome assembly we conducted reference-based assemblies for 34 other Sceloporus species to improve draft nuclear genomes assemblies from 1% coverage to 43% coverage on average. Across these species, typically >90% of reads mapped for species within 20 million years divergence from S. undulatus, this dropped to 75% reads mapped for species at 35 million years divergence. Finally we use RNAseq and whole genome resequencing data to compare the three assemblies as references, each representing an increased level of sequencing, cost and assembly efforts: Supernova Assembly with data from10X Genomics Chromium library; HiRise Assembly that added data from HiC library; and PBJelly Assembly that added data from PacBio sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq, but the chromosome-level scaffolds provided by the addition of the HiC data allowed the reference to be used for other whole genome analysis, including synteny and whole genome association mapping analyses. The addition of PacBio data provided negligible gains. Overall, these new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.