Genetic evidence for male-biased dispersal in the Qinghai toad-headed agamid Phrynocephalus vlangalii and its potential link to individual social interactions

Geography and past climate changes have shaped the evolution of a widespread lizard in arid Central Asia

The complex orogenic history and structure of Central Asia, coupled with Pleistocene glacial cycles have generated its stepwise aridification. Such events would have significantly influenced the evolution of many mid-latitude species in arid Central Asia (ACA). In this study, we employed two mitochondrial genes (CO1 and ND2) and genome-wide SNPs, coupled with ecological niche modeling, to investigate the lineage diversification and historical demography within a widespread lizard Phrynocepahlus helioscopus, and their associations with geography and past climate change. We obtained themtDNA dataset for 300 individuals from 96 localities within the known range of the lizard, among which 51 individuals from 27 localities were selected for generating the SNP dataset via genotyping-by-sequencing approach. Phylogenetic analyses of the concatenated mtDNA dataset revealed eight geographically correlated lineages that diverged by 4.21-10.41% for the CO1 gene, which were estimated to have coalesced ∼4.47 million years ago. However, we observed mito-nuclear discordance pattern regarding the population of Clade V (P. helioscopus sergeevi) from Tajikistan. Ancestral area estimations suggested that P. helioscopus originated from the Fergana Valley and then dispersed into the adjacent areas in ACA along with a history of multiple allopatric divergence processes, suggesting that Fergana may have been the cradle of diversification of P. helioscopus. The intensification of aridification across Central Asia during the Late Pliocene may have facilitated the rapid radiation of this arid-adapted lizard throughout this vast territory. Subsequently, the geological events (e.g., uplift of the Hissar-Alay, transgressions of the Caspian Sea) and geographic barriers (e.g., Amu Darya River, Zerarshan River) during the Pleistocene triggered the progressive diversification of P. helioscopus. Interestingly, Clade VIII (P. helioscopus varius) experienced rapid population growth coupled with range expansion while Clade IV (P. helioscopus cameranoi) underwent drastic population expansion associated with range contraction during the Last Glacial Maximum. In Clade IV, but not in Clade VIII, environmental turnover contributes more to mitochondrial genetic distinctiveness than geographic distance. Overall, the SNP dataset demonstrates that geographic distance plays a greater role than environmental distance. Both the mtDNA dataset and the SNP dataset suggest local-scale genetic differentiation in Clade IV and Clade VIII, revealing potential geographic barriers in the Ili River Valley and the Junggar Basin, respectively. Twenty-seven outlier SNPs associated with environmental factors (precipitation and temperature) were identified, which supports the signature of local adaptation to the arid desert environment. Finally, our finding suggests taxonomic implications, such as support for full species status for P. saidalievi (Clade II) and P. meridionalis (Clade I). Future analyses based on further evidence and increased taxon and geographic sampling should be carried out to corroborate our findings.

Genetic diversity and sex‐biased dispersal in the brown spotted pitviper ( Protobothrops mucrosquamatus ): Evidence from microsatellite markers

Dispersal plays a vital role in the geographical distribution, population genetic structure, quantity dynamics, and evolution of a species. Sex‐biased dispersal is common among vertebrates and many studies have documented a tendency toward male‐biased dispersal in mammals and female‐biased dispersal in birds. However, dispersal patterns in reptiles remain poorly understood. In this study, we explored the genetic diversity and dispersal patterns of the widely distributed Asian pitviper Protobothrops mucrosquamatus. In total, 16 polymorphic microsatellite loci were screened in 150 snakes (48 males, 44 females, 58 samples without sex information) covering most of their distribution. Microsatellite analysis revealed high genetic diversity in P. mucrosquamatus. Bayesian clustering of population assignment identified two major clusters for all populations, somewhat inconsistent with the mitochondrial DNA phylogeny of P. mucrosquamatus reported in previous research. Analyses based on 92 sex‐determined and 37 samples of P. mucrosquamatus from three small sites in Sichuan, China (Mingshan, Yibin, and Zizhong) consistently suggested female‐biased dispersal in P. mucrosquamatus, which is the first example of this pattern in snakes. The female‐biased dispersal patterns in P. mucrosquamatus may be explained by local resource competition.

Lack of conspicuous sex-biased dispersal patterns at different spatial scales in an Asian endemic goose species breeding in unpredictable steppe wetlands

Dispersal affects the spatial distribution and population structure of species. Dispersal is often male-biased in mammals while female-biased in birds, with the notable exception of the Anatidae. In this study, we tested genetic evidence for sex-biased dispersal (SBD) in the Swan Goose Anser cygnoides, an Asian endemic and IUCN vulnerable species, which has been increasingly restricted to breeding on Mongolian steppe wetlands. We analyzed the genotypes of 278 Swan Geese samples from 14 locations at 14 microsatellite loci. Results from assignment indices, analysis of molecular variance, and five other population descriptors all failed to support significant SBD signals for the Swan Goose at the landscape level. Although overall results showed significantly high relatedness within colonies (suggesting high levels of philopatry in both sexes), local male genetic structure at the 1,050 km distance indicated greater dispersal distance for females from the eastern sector of the breeding range. Hence, local dispersal is likely scale-dependent and female-biased within the eastern breeding range. These findings are intriguing considering the prevailing expectation for there to be female fidelity in most goose species. We suggest that while behavior-related traits may have facilitated the local genetic structure for the Swan Goose, several extrinsic factors, including the decreasing availability of the nesting sites and the severe fragmentation of breeding habitats, could have contributed to the absence of SBD at the landscape level. The long-distance molt migration that is typical of goose species such as the Swan Goose may also have hampered our ability to detect SBD. Hence, we urge further genetic sampling from other areas in summer to extend our results, complemented by field observations to confirm our DNA analysis conclusions about sex-specific dispersal patterns at different spatial scales in this species.

Extreme mito-nuclear discordance in a peninsular lizard: the role of drift, selection, and climate

Nuclear and mitochondrial genomes coexist within cells but are subject to different tempos and modes of evolution. Evolutionary forces such as drift, mutation, selection, and migration are expected to play fundamental roles in the origin and maintenance of diverged populations; however, divergence may lag between genomes subject to different modes of inheritance and functional specialization. Herein, we explore whole mitochondrial genome data and thousands of nuclear single nucleotide polymorphisms to evidence extreme mito-nuclear discordance in the small black-tailed brush lizard, Urosaurus nigricaudus, of the Peninsula of Baja California, Mexico and southern California, USA, and discuss potential drivers. Results show three deeply divergent mitochondrial lineages dating back to the later Miocene (ca. 5.5 Ma) and Pliocene (ca. 2.8 Ma) that likely followed geographic isolation due to trans-peninsular seaways. This contrasts with very low levels of genetic differentiation in nuclear loci (F_ST < 0.028) between mtDNA lineages. Analyses of protein-coding genes reveal substantial fixed variation between mitochondrial lineages, of which a significant portion comes from non-synonymous mutations. A mixture of drift and selection is likely responsible for the rise of these mtDNA groups, albeit with little evidence of marked differences in climatic niche space between them. Finally, future investigations can look further into the role that mito-nuclear incompatibilities and mating systems play in explaining contrasting nuclear gene flow.

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.