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

Abstract

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.