Project description:The complete DNA sequence of the mitochondrial genome of the bullsnake (Pituophis catenifer sayi) is presented and analyzed in this study. The genome is 17,193?bp in length, and contains 22 transfer RNA genes, 2 ribosomal RNA genes, 13 protein-coding genes and 2 control regions. The overall base composition of the H-strand is A (34.5%), T (26.2%), C (25.8%) and G (12.7%). The gene order and orientation of the mitogenome is consistent with other sequenced genomes from colubrid snakes. Phylogenetic analyses using the ML, NJ and MP methods for a set of colubrids, including every sequenced genus within the sub-family Colubrinae, produced identical trees. We show that Pituophis is most closely related to the American ratsnake genus Pantherophis.
Project description:Venomous animals have traditionally been studied from a proteomic (but also transcriptomic) perspective, often overlooking the study of venom from a genomic point of view until recently. The rise of genomics has led to an increase in the number of reference genomes for non-model organisms, including venomous taxa, enabling new questions on venom evolution from a genomic context. Although venomous snakes are the fundamental model system in venom research, the number of high-quality reference genomes in the group remains limited. In this study, we present a high-quality chromosome-level reference genome for the Arabian horned viper (Cerastes gasperettii), a highly venomous snake native to the Arabian Peninsula. Our highly-contiguous genome allowed us to explore macrochromosomal rearrangements within the Viperidae family, as well as across squamate reptile evolution. Furthermore, we identified a total of ten different toxins conforming the venom’s core, in line with our proteomic results. We also compared microsyntenic changes in the main toxin gene clusters with those of other venomous snake species, highlighting the pivotal role of gene duplication and loss in the emergence and diversification of the two main toxin families for Cerastes gasperettii. Using Illumina data, we reconstructed the demographic history and genome-wide diversity of the species, revealing how historical aridity likely drove population expansions. Finally, this study highlights the importance of using long-read sequencing as well as chromosome-level reference genomes to disentangle the origin and diversification of toxin families in venomous species.
Project description:Venomous animals have traditionally been studied from a proteomic (but also transcriptomic) perspective, often overlooking the study of venom from a genomic point of view until recently. The rise of genomics has led to an increase in the number of reference genomes for non-model organisms, including venomous taxa, enabling new questions on venom evolution from a genomic context. Although venomous snakes are the fundamental model system in venom research, the number of high-quality reference genomes in the group remains limited. In this study, we present a high-quality chromosome-level reference genome for the Arabian horned viper (Cerastes gasperettii), a highly venomous snake native to the Arabian Peninsula. Our highly-contiguous genome allowed us to explore macrochromosomal rearrangements within the Viperidae family, as well as across squamate reptile evolution. Furthermore, we identified a total of ten different toxins conforming the venom’s core, in line with our proteomic results. We also compared microsyntenic changes in the main toxin gene clusters with those of other venomous snake species, highlighting the pivotal role of gene duplication and loss in the emergence and diversification of the two main toxin families for Cerastes gasperettii. Using Illumina data, we reconstructed the demographic history and genome-wide diversity of the species, revealing how historical aridity likely drove population expansions. Finally, this study highlights the importance of using long-read sequencing as well as chromosome-level reference genomes to disentangle the origin and diversification of toxin families in venomous species.
Project description:On the North American Great Plains, several snake species reach their northern range limit where they rely on sparsely distributed hibernacula located in major river valleys. Independent colonization histories for the river valleys and barriers to gene flow caused by the lack of suitable habitat between them may have produced genetically differentiated snake populations. To test this hypothesis, we used 10 microsatellite loci to examine the population structure of two species of conservation concern in Canada: the eastern yellow-bellied racer (Coluber constrictor flaviventris) and bullsnake (Pituophis catenifer sayi) in 3 major river valleys in southern Saskatchewan. Fixation indices (FST) showed that populations in river valleys were significantly differentiated for both species (racers, FST = 0.096, P = 0.001; bullsnakes FST = 0.045-0.157, P = 0.001). Bayesian assignment (STRUCTURE) and ordination (DAPC) strongly supported genetically differentiated groups in the geographically distinct river valleys. Finer-scale subdivision of populations within river valleys was not apparent based on our data, but is a topic that should be investigated further. Our findings highlight the importance of major river valleys for snakes at the northern extent of their ranges, and raise the possibility that populations in each river valley may warrant separate management strategies.