The complete mitochondrial DNA sequence of Crotalus horridus (timber rattlesnake)

The Phylogenetic Relationships of Major Lizard Families Using Mitochondrial Genomes and Selection Pressure Analyses in Anguimorpha

Anguimorpha, within the order Squamata, represents a group with distinct morphological and behavioral characteristics in different ecological niches among lizards. Within Anguimorpha, there is a group characterized by limb loss, occupying lower ecological niches, concentrated within the subfamily Anguinae. Lizards with limbs and those without exhibit distinct locomotor abilities when adapting to their habitats, which in turn necessitate varying degrees of energy expenditure. Mitochondria, known as the metabolic powerhouses of cells, play a crucial role in providing approximately 95% of an organism's energy. Functionally, mitogenomes (mitochondrial genomes) can serve as a valuable tool for investigating potential adaptive evolutionary selection behind limb loss in reptiles. Due to the variation of mitogenome structures among each species, as well as its simple genetic structure, maternal inheritance, and high evolutionary rate, the mitogenome is increasingly utilized to reconstruct phylogenetic relationships of squamate animals. In this study, we sequenced the mitogenomes of two species within Anguimorpha as well as the mitogenomes of two species in Gekkota and four species in Scincoidea. We compared these data with the mitogenome content and evolutionary history of related species. Within Anguimorpha, between the mitogenomes of limbless and limbed lizards, a branch-site model analysis supported the presence of 10 positively selected sites: Cytb protein (at sites 183 and 187), ND2 protein (at sites 90, 155, and 198), ND3 protein (at site 21), ND5 protein (at sites 12 and 267), and ND6 protein (at sites 72 and 119). These findings suggested that positive selection of mitogenome in limbless lizards may be associated with the energy requirements for their locomotion. Additionally, we acquired data from 205 mitogenomes from the NCBI database. Bayesian inference (BI) and Maximum Likelihood (ML) trees were constructed using the 13 mitochondrial protein-coding genes (PCGs) and two rRNAs (12S rRNA and 16S rRNA) from 213 mitogenomes. Our phylogenetic tree and the divergence time estimates for Squamata based on mitogenome data are consistent with results from previous studies. Gekkota was placed at the root of Squamata in both BI and ML trees. However, within the Toxicofera clade, due to long-branch attraction, Anguimorpha and (Pleurodonta + (Serpentes + Acrodonta)) were closely related groupings, which might indicate errors and also demonstrate that mitogenome-based phylogenetic trees may not effectively resolve long-branch attraction issues. Additionally, we reviewed the origin and diversification of Squamata throughout the Mesozoic era, suggesting that Squamata originated in the Late Triassic (206.05 Mya), with the diversification of various superfamilies occurring during the Cretaceous period. Future improvements in constructing squamate phylogenetic relationships using mitogenomes will rely on identifying snake and acrodont species with slower evolutionary rates, ensuring comprehensive taxonomic coverage of squamate diversity, and increasing the number of genes analyzed.

Timber rattlesnake (Crotalus horridus): Biology, conservation, and envenomation in the Upper Mississippi River Valley (1982-2020)

The Timber Rattlesnake (Crotalus horridus) is the largest pit viper in the Northern United States and is the prominent venomous snake species indigenous to the bluff land habitats of the Upper Mississippi River Valley (UMRV). Conservation of C. horridus in this geographic region not only preserves the ecosystem's biodiversity and ecological balance, but also assures the continued study of their biomedically important venoms/toxins. Field studies of C. horridus biology and natural history performed from 1985 to 2015 in southeastern Minnesota and western Wisconsin along the Mississippi River showed populations have declined. Consequently, the implementation of improved conservation measures afforded the species protective status in both states. Historically, accounts of Timber Rattlesnake bites in the UMRV have been sparse, and medical consequences of envenomation have had limited documentation. However, in recent decades cases of envenomation by C. horridus have continued to occur. Retrospective analysis of clinical toxinology consultations documented from 1982 to 2020 on cases of envenomation by C. horridus in the UMRV revealed a very low incidence of bites annually and revealed that their venom can induce a rapid and precipitous decline in platelets.

Identification of Daboia siamensis venome using integrated multi-omics data

Snakebite, classified by World Health Organization as a neglected tropical disease, causes more than 100,000 deaths and 2 million injuries per year. Currently, available antivenoms do not bind with strong specificity to target toxins, which means that severe complications can still occur despite treatment. Moreover, the cost of antivenom is expensive. Knowledge of venom compositions is fundamental for producing a specific antivenom that has high effectiveness, low side effects, and ease of manufacture. With advances in mass spectrometry techniques, venom proteomes can now be analyzed in great depth at high efficiency. However, these techniques require genomic and transcriptomic data for interpreting mass spectrometry data. This study aims to establish and incorporate genomics, transcriptomics, and proteomics data to study venomics of a venomous snake, Daboia siamensis. Multiple proteins that have not been reported as venom components of this snake such as hyaluronidase-1, phospholipase B, and waprin were discovered. Thus, multi-omics data are advantageous for venomics studies. These findings will be valuable not only for antivenom production but also for the development of novel therapeutics.

The rise of genomics in snake venom research: recent advances and future perspectives

Snake venoms represent a danger to human health, but also a gold mine of bioactive proteins that can be harnessed for drug discovery purposes. The evolution of snakes and their venom has been studied for decades, particularly via traditional morphological and basic genetic methods alongside venom proteomics. However, while the field of genomics has matured rapidly over the past 2 decades, owing to the development of next-generation sequencing technologies, snake genomics remains in its infancy. Here, we provide an overview of the state of the art in snake genomics and discuss its potential implications for studying venom evolution and toxinology. On the basis of current knowledge, gene duplication and positive selection are key mechanisms in the neofunctionalization of snake venom proteins. This makes snake venoms important evolutionary drivers that explain the remarkable venom diversification and adaptive variation observed in these reptiles. Gene duplication and neofunctionalization have also generated a large number of repeat sequences in snake genomes that pose a significant challenge to DNA sequencing, resulting in the need for substantial computational resources and longer sequencing read length for high-quality genome assembly. Fortunately, owing to constantly improving sequencing technologies and computational tools, we are now able to explore the molecular mechanisms of snake venom evolution in unprecedented detail. Such novel insights have the potential to affect the design and development of antivenoms and possibly other drugs, as well as provide new fundamental knowledge on snake biology and evolution.

The complete mitochondrial genome of Bothrops jararaca (Reptilia, Serpentes, Viperidae)

The complete mitochondrial genome, containing 17,526 bp, was determined from the pitviper Bothrops jararaca. It is the first mitogenome for the most medically important genus of snake in Latin America. This mitogenome has common snake mitochondrial features such as a duplicated control region that has nearly identical sequences at two different locations of the mitogenome and a translocation of tRNA-Leu (UUR). Besides, we found a translocation of the tRNA-Pro compared to Colubridae snakes. Finally, an unusual possible duplication containing a tRNA-Phe was observed for the first time and may represent a marker of the genus.

Sequence and organization of the complete mitochondrial genome of the Ussuri mamushi (Gloydius ussuriensis)

The mitochondrial genome sequence of Ussuri mamushi is analyzed and presented publically for the first time. The genome is 17,208 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and 2 control regions. The overall base composition is A (32.4%), C (28.4%), T (26.3%) and G (12.9%). The base compositions clearly presented the A-C skew, which is the most obvious in the protein-coding genes. Mitochondrial genome analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees, indicating a close phylogenetic affinity of the 13 Crotalinae species. It appeared that no less than two major phyletic lineages were present in Crotalinae. The main clades within the Crotalinae include Protobothrops. A clade (G. brevicaudus, G. ussuriensis, G. intermedius, and G. sahatilis) with the Ovophis as the sister taxon to Protobothrops was supported by bootstrap values of 88%. The four Gloydius species formed a paraphyletic group with the high bootstrap value (100%) in all examinations.

The complete mitochondrial genome of Gloydius intermedius (Squamata: Viperidae: Crotalinae) from China

The mitochondrial genome sequence of Gloydius intermedius is analyzed and presented for the first time. The genome was 17, 226 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and 2 control region. The overall base composition was A (32.4%), C (28.8%), T (25.9%), and G (12.9%). The base compositions clearly presented the A-C skew, which was most obvious in the protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of G. intermedius. Mitochondrial genome analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees, indicating a close phylogenetic affinity of the 13 Crotalinae species. It appeared that no less than two major phyletic lineages were present in Crotalinae. The main clades within the Crotalinae supported are: A clade including the Protobothrops. A clade (G. brevicaudus, G. ussuriensis, G. intermedius, G. saxatilis) with the Ovophis appeared as the sister taxon to Protobothrops and was supported by bootstrap values of 88%. The four Gloydius species formed a paraphyletic group with the high bootstrap value (100 %) in all examinations.

Sequencing and analysis of the complete mitochondrial genome of Gloydius saxatilis (Squamata: Viperidae: Crotalinae)

The mitochondrial genome sequence of Gloydius saxatilis is analyzed and presented for the public for the first time. The genome was 17,218 bp in length and contained 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and 2 control regions. The overall base composition was A (32.3%), C (28.9%), T (25.8%), and G (13.0%). The base compositions presented clearly the A-C skew, which was most obviously in the protein-coding genes. The extended termination-associated sequence domain, the central conserved domain and the conserved sequence block domain are defined in the mitochondrial genome control region of G. saxatilis. Mitochondrial genomes analyses based on MP, ML, NJ and Bayesian analyses yielded identical phylogenetic trees, indicating a close phylogenetic affinity of the thirteen Crotalinae species. It appeared that no less than two major phyletic lineages were present in Crotalinae. The main clades within the Crotalinae supported are: A clade including the Protobothrops. A clade (G. brevicaudus, G. ussuriensis, G. intermedius, G. saxatilis) with the Ovophis as the sister taxon to Protobothrops and was supported by bootstrap values of 88%. The four Gloydius species formed a paraphyletic group with the high bootstrap value (100%) in all examinations.