The complete mitochondrial genome of an agamid lizard from the Afro-Asian subfamily agaminae and the phylogenetic position of Bufoniceps and Xenagama

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.

The story of a rock-star: multilocus phylogeny and species delimitation in the starred or roughtail rock agama, Laudakia stellio (Reptilia: Agamidae)

Situated at the junction of three continents, Europe, Asia and Africa, the Eastern Mediterranean is an ideal region to study the effects of palaeogeography, ecology and long human presence on animal evolution. Laudakia stellio (Squamata: Agamidae) is found across this region and offers an excellent opportunity for such studies. The high morphological variation across their range suggests that these lizards might represent a species complex. This is the first study exploring their evolutionary history, using molecular markers and individuals from all described subspecies. We employed the latest phylogenetic and species-delimitation methods to identify all distinct evolutionary lineages, their genetic variation and divergence times. The phenotypical diversity of L. stellio matches its genetic differentiation: almost all subspecies correspond to well-supported retrieved subclades and additional distinct lineages representing intermediate morphs have been retrieved. ‘Laudakia stellio’ represents three distinct evolutionary entities that diverged during the Plio-Pleistocene transition, which we propose as distinct species. One includes Greek and Turkish populations, as well as cryptic Anatolian lineages. The second comprises all other Near East populations and the third is endemic to Cyprus. Our results indicate a role of humans in shaping present distribution patterns, and highlight the importance of the Aegean, Anatolia and the Levant as glacial refugia and diversity hotspots.

Comparative Analysis of Mitochondrial Genomes in Two Subspecies of the Sunwatcher Toad-Headed Agama (Phrynocephalus helioscopus): Prevalent Intraspecific Gene Rearrangements in Phrynocephalus

Intraspecific rearrangements of mitochondrial genomes are rarely reported in reptiles, even in vertebrates. The sunwatcher toad-headed agama, Phryncoephalus helioscopus, can serve as an excellent model for investigating the dynamic mitogenome structure at intraspecific level. To date, seven subspecies of P. helioscopus are well recognized, but little is known about the mitogenomic evolution among different subspecies. In this study, complete mitogenomes of subspecies P. helioscopus varius II and P. helioscopus cameranoi were determined by next-generation sequencing, and another P. helioscopus varius I retrieved from GenBank was compiled for comparative analysis. The nucleotide composition and the codon usage are similar to those previously published from toad-headed agamas. P. helioscopus varius II and P. helioscopus cameranoi have 23 tRNA genes, including standard 22 tRNA genes and one extra tRNA-Phe (tRNA-Phe duplication). Gene order and phylogenetic analyses in the genus Phrynocephalus support prevalent intraspecific gene rearrangement in P. helioscopus and other congener species including P. erythrurus, P. vlangalii, and P. forsythii. Six different mitochondrial gene arrangements are observed in Phrynocephalus. Overall, the occurrence of rearrangements may result from multiple independent structural dynamic events. The split of the two subspecies in P. helioscopus was dated at approximately 2.34 million years ago (Ma). Two types of gene rearrangements are found in the three mitogenomes of P. helioscopus, and this intraspecific rearrangement phenomenon can be explained by the tandem duplication/random loss (TDRL) model. Post duplication, the alternative loss types can occur in 0.23–0.72 Ma, suggesting that the duplication and fixation of these rearrangements can occur quite quickly. These findings highlight the need for more mitogenomes at the population level in order to better understand the potentially rampant intraspecific mitogenomic reorganization in Phrynocephalus.

Structure and comparative analysis of the mitochondrial genomes of Liolaemus lizards with different modes of reproduction and ploidy levels

Liolaemus is the most specious genus of the Squamata lizards in South America, presenting exceptional evolutionary radiation and speciation patterns. This recent diversification complicates the formal taxonomic treatment and the phylogenetic analyses of this group, causing relationships among species to remain controversial. Here we used Next-Generation Sequencing to do a comparative analysis of the structure and organization of the complete mitochondrial genomes of three differently related species of Liolaemus and with different reproductive strategies and ploidy levels. The annotated mitochondrial genomes of ca. 17 kb are the first for the Liolaemidae family. Despite the high levels of sequence similarity among the three mitochondrial genomes over most of their lengths, the comparative analyses revealed variations at the stop codons of the protein coding genes and the structure of the tRNAs among species. The presence of a non-canonical dihydrouridine loop is a novelty for the pleurodonts iguanians. But the highest level of variability was observed in two repetitive sequences of the control region, which were responsible for most of the length heterogeneity of the mitochondrial genomes. These tandem repeats may be useful markers to analyze relationships of closely related species of Liolaemus and related genera and to conduct population and phylogenetic studies.

Cenozoic aridization in Central Eurasia shaped diversification of toad-headed agamas (Phrynocephalus; Agamidae, Reptilia)

We hypothesize the phylogenetic relationships of the agamid genus Phrynocephalus to assess how past environmental changes shaped the evolutionary and biogeographic history of these lizards and especially the impact of paleogeography and climatic factors. Phrynocephalus is one of the most diverse and taxonomically confusing lizard genera. As a key element of Palearctic deserts, it serves as a promising model for studies of historical biogeography and formation of arid habitats in Eurasia. We used 51 samples representing 33 of 40 recognized species of Phrynocephalus covering all major areas of the genus. Molecular data included four mtDNA (COI, ND2, ND4, Cytb; 2,703 bp) and four nuDNA protein-coding genes (RAG1, BDNF, AKAP9, NKTR; 4,188 bp). AU-tests were implemented to test for significant differences between mtDNA- and nuDNA-based topologies. A time-calibrated phylogeny was estimated using a Bayesian relaxed molecular clock with nine fossil calibrations. We reconstructed the ancestral area of origin, biogeographic scenarios, body size, and the evolution of habitat preference. Phylogenetic analyses of nuDNA genes recovered a well-resolved and supported topology. Analyses detected significant discordance with the less-supported mtDNA genealogy. The position of Phrynocephalus mystaceus conflicted greatly between the two datasets. MtDNA introgression due to ancient hybridization best explained this result. Monophyletic Phrynocephalus contained three main clades: (I) oviparous species from south-western and Middle Asia; (II) viviparous species of Qinghai–Tibetan Plateau (QTP); and (III) oviparous species of the Caspian Basin, Middle and Central Asia. Phrynocephalus originated in late Oligocene (26.9 Ma) and modern species diversified during the middle Miocene (14.8–13.5 Ma). The reconstruction of ancestral areas indicated that Phrynocephalus originated in Middle East–southern Middle Asia. Body size miniaturization likely occurred early in the history of Phrynocephalus. The common ancestor of Phrynocephalus probably preferred sandy substrates with the inclusion of clay or gravel. The time of Agaminae radiation and origin of Phrynocephalus in the late Oligocene significantly precedes the landbridge between Afro-Arabia and Eurasia in the Early Miocene. Diversification of Phrynocephalus coincides well with the mid-Miocene climatic transition when a rapid cooling of climate drove progressing aridification and the Paratethys salinity crisis. These factors likely triggered the spreading of desert habitats in Central Eurasia, which Phrynocephalus occupied. The origin of the viviparous Tibetan clade has been associated traditionally with uplifting of the QTP; however, further studies are needed to confirm this. Progressing late Miocene aridification, the decrease of the Paratethys Basin, orogenesis, and Plio–Pleistocene climate oscillations likely promoted further diversification within Phrynocephalus. We discuss Phrynocephalus taxonomy in scope of the new analyses.

Historical colonization and dispersal limitation supplement climate and topography in shaping species richness of African lizards (Reptilia: Agaminae)

To what extent deep-time dispersal limitation shapes present-day biodiversity at broad spatial scales remains elusive. Here, we compiled a continental dataset on the distributions of African lizard species in the reptile subfamily Agaminae (a relatively young, Neogene radiation of agamid lizards which ancestors colonized Africa from the Arabian peninsula) and tested to what extent historical colonization and dispersal limitation (i.e. accessibility from areas of geographic origin) can explain present-day species richness relative to current climate, topography, and climate change since the late Miocene (~10 mya), the Pliocene (~3 mya), and the Last Glacial Maximum (LGM, 0.021 mya). Spatial and non-spatial multi-predictor regression models revealed that time-limited dispersal via arid corridors is a key predictor to explain macro-scale patterns of species richness. In addition, current precipitation seasonality, current temperature of the warmest month, paleo-temperature changes since the LGM and late Miocene, and topographic relief emerged as important drivers. These results suggest that deep-time dispersal constraints — in addition to climate and mountain building — strongly shape current species richness of Africa’s arid-adapted taxa. Such historical dispersal limitation might indicate that natural movement rates of species are too slow to respond to rates of ongoing and projected future climate and land use change.

A hybrid phylogenetic-phylogenomic approach for species tree estimation in African Agama lizards with applications to biogeography, character evolution, and diversification

Africa is renowned for its biodiversity and endemicity, yet little is known about the factors shaping them across the continent. African Agama lizards (45 species) have a pan-continental distribution, making them an ideal model for investigating biogeography. Many species have evolved conspicuous sexually dimorphic traits, including extravagant breeding coloration in adult males, large adult male body sizes, and variability in social systems among colorful versus drab species. We present a comprehensive time-calibrated species tree for Agama, and their close relatives, using a hybrid phylogenetic-phylogenomic approach that combines traditional Sanger sequence data from five loci for 57 species (146 samples) with anchored phylogenomic data from 215 nuclear genes for 23 species. The Sanger data are analyzed using coalescent-based species tree inference using (*)BEAST, and the resulting posterior distribution of species trees is attenuated using the phylogenomic tree as a backbone constraint. The result is a time-calibrated species tree for Agama that includes 95% of all species, multiple samples for most species, strong support for the major clades, and strong support for most of the initial divergence events. Diversification within Agama began approximately 23 million years ago (Ma), and separate radiations in Southern, East, West, and Northern Africa have been diversifying for >10Myr. A suite of traits (morphological, coloration, and sociality) are tightly correlated and show a strong signal of high morphological disparity within clades, whereby the subsequent evolution of convergent phenotypes has accompanied diversification into new biogeographic areas.

Giant lizards occupied herbivorous mammalian ecospace during the Paleogene greenhouse in Southeast Asia

Mammals dominate modern terrestrial herbivore ecosystems, whereas extant herbivorous reptiles are limited in diversity and body size. The evolution of reptile herbivory and its relationship to mammalian diversification is poorly understood with respect to climate and the roles of predation pressure and competition for food resources. Here, we describe a giant fossil acrodontan lizard recovered with a diverse mammal assemblage from the late middle Eocene Pondaung Formation of Myanmar, which provides a historical test of factors controlling body size in herbivorous squamates. We infer a predominately herbivorous feeding ecology for the new acrodontan based on dental anatomy, phylogenetic relationships and body size. Ranking body masses for Pondaung Formation vertebrates indicates that the lizard occupied a size niche among the larger herbivores and was larger than most carnivorous mammals. Paleotemperature estimates of Pondaung Formation environments based on the body size of the new lizard are approximately 2–5°C higher than modern. These results indicate that competitive exclusion and predation by mammals did not restrict body size evolution in these herbivorous squamates, and elevated temperatures relative to modern climates during the Paleogene greenhouse may have resulted in the evolution of gigantism through elevated poikilothermic metabolic rates and in response to increases in floral productivity.

Mitochondrial genomes of acrodont lizards: timing of gene rearrangements and phylogenetic and biogeographic implications

Background

Acrodonta consists of Agamidae and Chamaeleonidae that have the characteristic acrodont dentition. These two families and Iguanidae sensu lato are members of infraorder Iguania. Phylogenetic relationships and historical biogeography of iguanian lizards still remain to be elucidated in spite of a number of morphological and molecular studies. This issue was addressed by sequencing complete mitochondrial genomes from 10 species that represent major lineages of acrodont lizards. This study also provided a good opportunity to compare molecular evolutionary modes of mitogenomes among different iguanian lineages.

Results

Acrodontan mitogenomes were found to be less conservative than iguanid counterparts with respect to gene arrangement features and rates of sequence evolution. Phylogenetic relationships were constructed with the mitogenomic sequence data and timing of gene rearrangements was inferred on it. The result suggested highly lineage-specific occurrence of several gene rearrangements, except for the translocation of the tRNA^Pro gene from the 5' to 3' side of the control region, which likely occurred independently in both agamine and chamaeleonid lineages. Phylogenetic analyses strongly suggested the monophyly of Agamidae in relation to Chamaeleonidae and the non-monophyly of traditional genus Chamaeleo within Chamaeleonidae. Uromastyx and Brookesia were suggested to be the earliest shoot-off of Agamidae and Chamaeleonidae, respectively. Together with the results of relaxed-clock dating analyses, our molecular phylogeny was used to infer the origin of Acrodonta and historical biogeography of its descendant lineages. Our molecular data favored Gondwanan origin of Acrodonta, vicariant divergence of Agamidae and Chamaeleonidae in the drifting India-Madagascar landmass, and migration of the Agamidae to Eurasia with the Indian subcontinent, although Laurasian origin of Acrodonta was not strictly ruled out.

Conclusions

We detected distinct modes of mitogenomic evolution among iguanian families. Agamidae was highlighted in including a number of lineage-specific mitochondrial gene rearrangements. The mitogenomic data provided a certain level of resolution in reconstructing acrodontan phylogeny, although there still remain ambiguous relationships. Our biogeographic implications shed a light on the previous hypothesis of Gondwanan origin of Acrodonta by adding some new evidence and concreteness.

Partitioned Bayesian analyses, dispersal-vicariance analysis, and the biogeography of Chinese toad-headed lizards (Agamidae: Phrynocephalus): a re-evaluation

The toad-headed lizards of genus Phrynocephalus are distributed from northwestern China to Turkey and are one of the major components of the central Asian desert fauna. To date, published morphological and molecular phylogenetic hypotheses of Phrynocephalus are only partially congruent, and the relationships within the genus are still far from clear. We re-analyzed published mitochondrial gene sequence data (12S, 16S, cyt b, ND4-tRNA(Leu)) by employing partition-specific modeling in a combined DNA analysis to clarify existing gaps in the phylogeny of Chinese Phrynocephalus. Using this phylogenetic framework, we inferred the genus' historical biogeography by using weighted ancestral-area analysis and dispersal-vicariance analysis in combination with a Bayesian relaxed molecular-clock approach and paleogeographical data. The partitioned Bayesian analyses support the monophyly of Phrynocephalus and its sister-group relationship with Laudakia. An earlier finding demonstrating the monophyly of the viviparous group is corroborated. However, our hypothesis of internal relationships of the oviparous group differs from a previous hypothesis as our results do not support monophyly of the oviparous taxa. Instead, the viviparous taxa form a clade with many oviparous taxa exclusive of P. helioscopus and P. mystaceus. Our results also suggest that: (1) P. putjatia is a valid species, comprising populations from Guide, Qinghai Province and Tianzhu, Gansu Province; (2) P. hongyuanensis is not a valid species, synonymized instead with P. vlangalii; (3) P. zetangensis is not a valid species and should be included in P. theobaldi; (4) the population occurring in Kuytun, Xinjiang Uygur Autonomous Region is recognized as P. guttatus instead of P. versicolor; and (5) the Lanzhou population of P. frontalis is part of P. przewalskii. Congruent with previous hypotheses, the uplift of the Tibetan Plateau played a fundamental role in the diversification of Phrynocephalus. An evolutionary scenario combining aspects of vicariance and dispersal is necessary to explain the distribution of Phrynocephalus. Bayesian divergence-time estimation suggests that Phrynocephalus originated at the Middle-Late Miocene boundary (15.16-10.4 Ma), and diversified from Late Miocene to Pleistocene from a center of origin in Central Asia, Tarim Basin, and Junggar Basin temperate desert, followed by several rapid speciation events in a relatively short time. The proposed biogeographic scenarios also indicate that the Tarim Basin desert may be the secondary diversification center, followed by Junggar Basin temperate desert and Alashan Plateau temperate desert. In the viviparous group, the allopatric speciation of P. theobaldi and P. vlangalii may have been caused by the uplifting of Tanggula Mountain Ranges. In addition, the results of this study make an important contribution to understanding the uplift of the Tibetan Plateau and Tian Shan Mountains and the biogeography of the entire region.