Lineage Diversification and Population Dynamics of the Qinghai Toad-Headed Agama (Phrynocephalus vlangalii) on the Qinghai–Tibet Plateau, with Particular Attention to the Northern Slope of the Kunlun–Arjin Mountains

The Kunlun, Arjin, and Qilian mountain ranges mark the northern edge of the Qinghai&ndash:Tibet Plateau (QTP), where rapid uplift and Quaternary glacial cycles have shaped a unique cold desert ecosystem and species distribution. Despite sampling challenges, phylogeographic studies are crucial for understanding reptile populations such as the Qinghai toad-headed agama (Phrynocephalus vlangalii), a viviparous lizard with limited dispersal and multiple subspecies in the northeastern QTP. Our fieldwork identified populations of P. vlangalii on the northern slope of the Kunlun&ndash:Arjin Mountains, similar to the controversial subspecies P. v. lidskii. We analyzed 130 individuals from the northern slope of the Kunlun&ndash:Arjin&ndash:Qilian Mountains and 253 individuals from GenBank, using three mitochondrial genes and two nuclear genes to assess intraspecific differentiation and demographic history. We found high haplotype diversity and low nucleotide diversity in P. vlangalii, with phylogenetic analyses revealing six distinct clades. Clade VI, confirmed as P. v. lidskii, and Clade IV, a new genetic lineage, were identified alongside three recognized subspecies. Genetic variation was largely attributed to clade splitting, indicating significant divergence. The Mantel test indicated that geographical and environmental factors drove population differentiation. Bayesian molecular clock analysis suggested that the most recent common ancestor of P. vlangalii lived 2.55 million years ago, influenced by the Qinghai&ndash:Tibet Movement and glacial cycles. Demographic history and ecological niche modeling (ENM) indicated no population decline during the Last Glacial Maximum, supporting the glacial maximum expansion model, with ENM predicting future habitat expansion for P. vlangalii. In addition, morphological data from 13 meristic and 15 metric characters confirmed clade differences. Our findings significantly advance our understanding of P. vlangalii diversification, population dynamics and response to geological and climatic changes in the QTP.