Phylogeography and Ecological Niche Modeling of the Alashan Pit Viper (Gloydius cognatus; Reptilia, Viperidae) in Northwest China and Adjacent Areas.

The joint impacts of historical geological events and Quaternary climatic oscillations in Northwest China on species evolution have been examined extensively in plant under a phylogeographic perspective. However, animal phylogeographic analyses in this region are still limited. The Alashan pit viper, Gloydius cognatus, occurs primarily in arid Northwest China and adjacent areas. Based on variation at two mtDNA genes (ND4 and Cytb) in 27 individuals representing 24 populations, the spatial genetic structure and demographic history of G. cognatus were examined across its geographic range. Phylogenetic analyses revealed two well-supported allopatric clades (each with two distinct subclades/lineages), distributed across the southern (Qaidam Basin, Lanzhou Basin, and Zoige Basin [S1]; Loess Plateau [S2]) and northern (Ily Basin [N1]; Junggar Basin and Mongolian Plateau [N2]) regions. AMOVA analysis demonstrated that over 76% of the observed genetic variation was related to these lineage splits, indicating substantial genetic differentiation among the four lineages. A strong pattern of isolation-by-distance across the sampling populations suggested that geographic distance principally shaped the genetic structure. The four lineages diverged by 0.9-2.2% for the concatenated data, which were estimated to have coalesced ~1.17 million years ago (Mya), suggesting that the expansions of the Badain Jaran, Tengger, and Mu Us deserts during the Xixiabangma glaciation likely interrupted gene flow and triggered the observed divergence in the southern and northern regions. Subsequently, the early Pleistocene integration of the Yellow River and associated deserts expansion promoted the differentiation of S1 and S2 lineages (~0.9 Mya). Both mitochondrial evidence and ecological niche modeling (ENM) reject the signature of demographic and range contractions during the LGM for G. cognatus. In addition, ENM predicts that the suitable habitat of G. cognatus will contract in the future. As such, the conservation and management of ESUs should be a priority. Our findings provide the first insights on the lineage diversification and population dynamics of the Alashan pit viper in relation to geological history and Pleistocene climatic oscillations in arid Northwest China.

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 Phrynocephalus helioscopus, and their associations with geography and past climate change. We obtained the mtDNA 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 the 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, Zeravshan 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.

Does the Dzungarian racerunner ( Eremias dzungarica Orlova, Poyarkov, Chirikova, Nazarov, Munkhbaatar, Munkhbayar & Terbish, 2017) occur in China? Species delimitation and identification with DNA barcoding and morphometric analyses.

The Eremias multiocellata-przewalskii species complex is a viviparous group in the genus Eremias, and a well-known representative of taxonomically complicated taxa. Within this complex, a new species - E. dzungarica (Orlova et al., 2017) - has been described recently from western Mongolia and eastern Kazakhstan, with an apparent distribution gap in northwestern China. In this study, we used an integrative taxonomic framework to address whether E. dzungarica indeed occurs in China. Thirty specimens previously classified as E. multiocellata were collected in eastern Kazakhstan and the adjacent Altay region in China. The cytochrome c oxidase I ( COI) barcodes were sequenced and compiled with those from Orlova et al. (2017) and analyzed with the standard and diverse barcoding techniques. We detected an absence of a barcoding gap in this complex, which indicates potential cryptic species in Eremias sp. 3 with high intraspecific diversity and multiple recently evolved species in Clade A. Both BIN and GMYC suggested an unrealistically large number of species (23 and 26, respectively), while ABGD, mPTP and BPP indicated a more conservative number of species (10, 12, and 15, respectively), largely concordant with the previously defined species-level lineages according to phylogenetic trees. Based on molecular phylogeny and morphological examination, all 30 individuals collected in this study were reliably identified as E. dzungarica - a distinct species - confirming the occurrence of this species in the Altay region, Xinjiang, China. Potentially owing to the larger sample size in this study, our morphological analyses revealed many inconsistencies with the original descriptions of E. dzungarica, which were primarily associated with sexual dimorphism and a broader range of values for various traits.

Complete mitochondrial genome of Stummer's Racerunner (Eremias stummeri) from Kazakhstan.

The whole mitochondrial genome was determined from a viviparous racerunner, Eremias stummeri, which was collected from southeast Kazakhstan. The mitogenome sequence was 19 602 bp in size, containing 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a control region, which is similar to the typical mtDNA of vertebrates. Mitochondrial genomes analyses using maximum parsimony and Bayesian analyses yielded identical phylogenetic trees, indicating a close phylogenetic affinity of the sampled taxa in genus Eremias. Monophyly of both Eremias and its viviparous group is recovered. The complete mitogenome sequence of E. stummeri provides fundamental data for resolving phylogeneitc and genetic problems related to Eremias viviparity.