Range-Wide Phylogeography and Ecological Niche Modeling Provide Insights into the Evolutionary History of the Mongolian Racerunner (Eremias argus) in Northeast Asia.

The Mongolian racerunner, Eremias argus, is a small lizard endemic to Northeast Asia that can serve as an excellent model for investigating how geography and past climate change have jointly influenced the evolution of biodiversity in this region. To elucidate the processes underlying its diversification and demography, we reconstructed the range-wide phylogeographic pattern and evolutionary trajectory, using phylogenetic, population genetic, landscape genetic, Bayesian phylogeographic reconstruction and ecological niche modeling approaches. Phylogenetic analyses of the mtDNA cyt b gene revealed eight lineages that were unbounded by geographic region. The genetic structure of E. argus was mainly determined by geographic distance. Divergence dating indicated that E. argus and E. brenchleyi diverged during the Mid-Pliocene Warm Period. E. argus was estimated to have coalesced at~0.4351 Ma (Marine Isotope Stage 19). Bayesian phylogeographic diffusion analysis revealed out-of-Inner Mongolia and rapid colonization events from the end of the Last Interglacial to the Last Glacial Maximum, which is consistent with the expanded suitable range of the Last Glacial Maximum. Pre-Last Glacial Maximum growth of population is presented for most lineages of E. argus. The Glacial Maximum contraction model and the previous multiple glacial refugia hypotheses are rejected. This may be due to an increase in the amount of climatically favorable habitats in Northeast Asia. Furthermore, E. argus barbouri most likely represents an invalid taxon. The present study is the first to report a range-wide phylogeography of reptiles over such a large region in Northeast Asia. Our results make a significant contribution towards understanding the biogeography of the entire Northeast Asia.

Transcriptomic profiling revealed immune-related signaling pathways in response to experimental infection of Leishmania donovani in two desert lizards from Northwest China.

Little is known about the immune response of lizards to Leishmania parasties. In this study, we conducted the first liver transcriptome analysis of two lizards (Phrynocephalus przewalskii and Eremias multiocellata) challenged with L. donovani, endemic to the steppe desert region of northwestern China. Our results revealed that multiple biological processes and immune-related signaling pathways are closely associated with the immune response to experimental L. donovani infection in the two lizards, and that both lizards show similar changes to mammals in terms of immunity to Leishmania. However, the interspecific divergence of the two lizards leads to different transcriptomic changes. In particular, in contrast to P. przewalskii, the challenged E. mutltiocellata was characterized by the induction of down-regulation of most DEGs. These findings will contribute to the scarce resources on lizard immunity and provide a reference for further research on immune mechanisms in reptiles.

Genetic Structure and Population History of the Zaisan Toad-Headed Agama (Phrynocephalus melanurus) Inferred from Mitochondrial DNA.

The agamid lizard Phrynocephalus melanurus is restricted to Northwest China (Dzungar Basin) and the adjacent Eastern Kazakhstan (Zaisan and Alakol basins). To elucidate the phylogeography of P. melanurus, we obtained the mitochondrial DNA COI segments of 175 sampled lizards from 44 localities across the whole distribution. Phylogenetic analyses revealed two main Clades comprising five geographically structured lineages (I, IIa, IIb1, IIb2, and IIb3) that fit an isolation-by-distance (IBD) model. The divergence from the most recent common ancestor was dated to ~1.87 million years ago (Ma). Demographic analyses demonstrated lineage-specific response to past climate change: stable population for Clade I, Subclade IIb1; past population expansion for IIb3 since 0.18 Ma, respectively. Bayesian phylogeographic diffusion analyses detected initial spreading at the Saur Mount vicinity, approximately 1.8 Ma. Historical species distribution model (SDM) projected expansion of the suitable habitat in the last interglacial and shift and contraction in the last glacial maximum and Holocene epochs. The SDM predicted a drastic reduction in suitable area throughout the range as a response to future climate change. Our findings suggest that the evolution of P. melanurus followed a parapatric divergence with subsequent dispersal and adaptation to cold and dry environments during the Quaternary. Overall, this work improves our understanding of the lineage diversification and population dynamics of P. melanurus, providing further insights into the evolutionary processes that occurred in Northwest China and adjacent Eastern Kazakhstan.

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.

A complete mitogenome of the Przewalski's Wonder Gecko (Teratoscincus przewalskii) from the Junggar Basin in Northwest China with its phylogenetic implications.

A complete mitogenome of the Przewalski's Wonder Gecko (Teratoscincus przewalskii) from the Junggar Basin in Northwest China was determined by using polymerase chain reaction and directly sequenced with the primer walking method. The total length was 17,184 bp, containing 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR). The order and structure of the genes were identical to those of congeners. The 13 PCGs contain four start codons (ATG, GTG, ATA, and ATC), three complete stop codons (TAA, TAG, and AGG), and two incomplete stop codons (T-, TA-). The concatenated PCGs were used to perform Bayesian phylogenetic analyses together with mitogenome data of the family Sphaerodactylidae and related representative taxa available in GenBank. The resulting tree recovered the monophyly of Sphaerodactylidae, and confirmed the sister relationship between T. przewalskii and T. roborowskii with strong support. The newly determined mitogenome will provide fundamental data for understanding the population genetic structure of T. pzrewalskii in particular, and the mitochondrial DNA evolution in Sphaerodactylidae in general.

Sex-Linked Loci on the W Chromosome in the Multi-Ocellated Racerunner (Eremias multiocellata) Confirm Genetic Sex-Determination Stability in Lacertid Lizards.

The multi-ocellated racerunner, Eremias multiocellata, was considered to have temperature-dependent sex determination (TSD), as its sex ratio can be influenced at different temperatures. However, such an observation contrasts with recent findings that suggest TSD is less common than previously thought. Here, a genotyping-by-sequencing (GBS) approach was employed to identify sex-linked markers in the E. multiocellata, for which the mechanism choice of TSD or GSD is still controversial. We preliminarily identified 119 sex-linked markers based on sex-associated sex-specific sequences, 97% of which indicated female heterogamety. After eliminating the false positives, 38 sex-linked markers were recognized, all of which showed the ZW/ZZ system. Then, eight of the novel markers were verified by PCR amplification from 15 populations of E. multiocellata, which support the GSD in E. multiocellata without geographic variation. To test the conservation of sex chromosome in Eremias, the eight markers were further cross-tested by PCR amplification in 10 individuals of the Mongolian racerunner (Eremias argus), two of which exhibited cross-utility. The novel sex-linked markers could be mapped on the W chromosome of the sand lizard (Lacerta agilis). Our finding that the sex-linked markers are shared in closely related species, along with a conserved synteny of the W chromosome, further supports the homology and conservation of sex chromosomes in the lacertid lizards.

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.

Complete Mitogenome Analysis of Five Leafhopper Species of Idiocerini (Hemiptera: Cicadellidae).

Insect mitochondrial genomes (mitogenomes) are of great interest in exploring molecular evolution, phylogenetics, and biogeography. So far, only 12 mitogenomes of the leafhopper tribe Idiocerini have been released in GenBank, although the tribe comprises 488 known species including some agricultural, forestry, and horticultural pests. In order to compare and analyze the mitochondrial genome structure of Idiocerini and even the selective pressure of 13 protein-coding genes (PCGs) of the family Cicadellidae, the complete mitogenomes of five species including Nabicerus dentimus, Sahlbergotettix salicicola, Podulmorinus opacus, Podulmorinus consimilis, and a new species of a new genus were determined by next-generation sequencing. The size of the newly determined mitogenomes ranged from 14,733 bp to 15,044 bp, comprising the standard set of 13 PCGs, 22 transfer RNA genes, two ribosomal RNA genes, and a long non-coding control region (CR). The extent of purifying selection presented different pictures in the tribe and the family. The less pronounced genes (0.5 < dN/dS < 1) were nad5 and nad4l in Idiocerin, whereas in the family Cicadellidae including the sequences of Idiocerin, nad1-nad6 and cox1 genes were less pronounced. The codon encoding leucine was the most common in all species, and the codon encoding serine 1 was the most common in all species except for P. opacus. Interestingly, in P. opacus, another of the most common codons is that encoding serine 2. Among the 17 examined species of the Idiocerini, 14 species contained the tandem repeats, and 11 species of them contained the motif "TTATA". These findings will promote research on the structure and evolution of the mitochondrial genome and highlight the need for more mitogenomes in Cicadellidae.

Complete mitochondrial genome of the Kyrghyz racerunner (Eremias nikolskii Bedriaga, 1905) from Kyrgyzstan.

The complete mitochondrial genome (mitogenome) of the Kyrghyz racerunner (Eremias nikoskii Bedriaga, 1905) from Kyrgyzstan was determined for the first time by next-generation sequencing. The mitogenome was 20,840 bp in length and comprised the standard set of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. The 13 concatenated PCGs were used to implement Bayesian phylogenetic analyses together with some congeners and three representative lacertids retrieved from GenBank. The monophyly of both Eremias and its viviparous group was recovered in the Bayesian phylogenetic tree, while the subgenus Pareremias was paraphyletic with respect to E. nikoskii. The mitogenome of E. nikoskii will faciliate the research on species delimitation, molecular evolution, and phylogenetic inference in the racerunner lizards.

Complete Mitochondrial Genomes of Five Racerunners (Lacertidae: Eremias) and Comparison with Other Lacertids: Insights into the Structure and Evolution of the Control Region.

Comparative studies on mitochondrial genomes (mitogenomes) as well as the structure and evolution of the mitochondrial control region are few in the Lacertidae family. Here, the complete mitogenomes of five individuals of Eremias scripta (2 individuals), Eremias nikolskii, Eremias szczerbaki, and Eremias yarkandensis were determined using next-generation sequencing and were compared with other lacertids available in GenBank. The circular mitogenomes comprised the standard set of 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes and a long non-coding control region (CR). The extent of purifying selection was less pronounced for the COIII and ND2 genes in comparison with the rest of the PCGs. The codons encoding Leucine (CUN), Threonine, and Isolecucine were the three most frequently present. The secondary structure of rRNA of Lacertidae (herein, E. scripta KZL15 as an example) comprised four domains and 28 helices for 12S rRNA, with six domains and 50 helices for 16S rRNA. Five types and twenty-one subtypes of CR in Lacertidae were described by following the criteria of the presence and position of tandem repeats (TR), termination-associated sequence 1 (TAS1), termination-associated sequence 2 (TAS2), conserved sequence block 1 (CBS1), conserved sequence block 2 (CSB2), and conserved sequence block 3 (CSB3). The compositions of conserved structural elements in four genera, Acanthodactylus, Darevskia, Eremias, and Takydromus, were further explored in detail. The base composition of TAS2 - TATACATTAT in Lacertidae was updated. In addition, the motif "TAGCGGCTTTTTTG" of tandem repeats in Eremias and the motif "GCGGCTT" in Takydromus were presented. Nucleotide lengths between CSB2 and CSB3 remained 35 bp in Eremias and Darevskia. The phylogenetic analyses of Lacertidae recovered the higher-level relationships among the three subfamilies and corroborated a hard polytomy in the Lacertinae phylogeny. The phylogenetic position of E. nikolskii challenged the monophyly of the subgenus Pareremias within Eremias. Some mismatches between the types of CR and their phylogeny demonstrated the complicated evolutionary signals of CR such as convergent evolution. These findings will promote research on the structure and evolution of the CR and highlight the need for more mitogenomes in Lacertidae.

The complete mitochondrial genome of Eremias yarkandensis (Reptilia, Squamata, Lacertidae) from Kyrgyzstan.

The Yarkand racerunner, Eremias yarkandensis Blandford, 1875, is only distributed in China and Kyrgyzstan. Its complete mitogenome was determined by next-generation sequencing for the first time. The mitogenome was 18,743 bp in length, including 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes, and 1 control region. Its gene arrangement was similar to the typical mtDNA of vertebrates. The 13 concatenated PCGs were used to perform Bayesian phylogenetic analyses together with several congeners as well as two representative species of Lacerta with mitogenome data in GenBank. The resulting phylogenetic tree recovered the monophyly of both Eremias and its viviparous group, with E. yarkandensis being more closely related to E. przewalskii than to E. dzungarica. The mitogenome of E. yarkandensis will provide fundamental data for the exploration of the mitogenome evolution in racerunners.

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.

The complete mitochondrial genome of Metidiocerus impressifrons (Kirschbaum, 1868) (Hemiptera: Cicadellidae: Idiocerinae).

We determined the complete mitochondrial genome (mitogenome) of the leafhopper Metidiocerus impressifrons by next-generation sequencing. The mitogenome sequence was 16,426 bp in length and consists of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and a control region. Moreover, the nucleotide composition biases toward A and T, which together made up 78.2% of the entirety. The complete mitochondrial genomes of Metidiocerus impressifrons and other 27 species were used for phylogenetic analysis using the Bayesian method. The above results would facilitate our understanding of the evolution of Idiocerinae mitochondrial genome.

The complete mitochondrial genome of Eremias dzungarica (Reptilia, Squamata, Lacertidae) from the Junggar Basin in Northwest China.

The complete mitochondrial genome (mitogenome) of Eremias dzungarica from the Junggar Basin in Northwest China was determined for the first time by next-generation sequencing. The total length is 19,899 bp, containing13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes and one control region. This gene arrangement is similar to the typical mtDNA of vertebrates. The 13 concatenated PCGs were used to perform Bayesian phylogenetic analyses together with several congeners as well as ten representative lacertids with mitogenome data available in GenBank. The resulting phylogenetic tree supports the monophyly of viviparous species in Eremias, with E. dzungarica being more closely related to E. przewalskii than to E. multiocellata. The mitogenome of E. dzungarica will provide fundamental data for the exploration of the mitogenome evolution in lacertids.

A new species of Psammophis (Serpentes: Psammophiidae) from the Turpan Basin in northwest China.

An adult sand snake specimen was collected during a herpetofaunal survey conducted in the Turpan Basin in northwest China. Phylogenetic analyses revealed that this specimen, along with other snake sloughs and skins collected from different localities in the Turpan Basin formed a clade that is sister to Psammophis lineolatus. This taxon exhibited substantial divergence from its congeners (P. lineolatus and P. condanarus) with uncorrelated p-distances ranging from 11.9 ± 0.9% to 15.8 ± 1.6% for the ND4 gene and from 10.2 ± 0.8% to 13.8 ± 1.1% for the Cytb gene. Given the genetic differences along with morphological differences, we describe the specimen from the Turpan Basin as Psammophis turpanensis sp. nov. We provide detailed morphological descriptions, and compare this specimen with five Asian sand snakes and the Afro-Asian Sand Snake, P. schokari. In addition, we provide brief comments on the biogeography of Psammophis in China.

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.

Records of three mammal tick species parasitizing an atypical host, the multi-ocellated racerunner lizard, in arid regions of Xinjiang, China.

BACKGROUND: Ticks are ectoparasites that feed on blood of a broad taxonomical range of terrestrial and flying vertebrates and are distributed across a wide range of environmental settings. To date, the species identity, diversity, and relationships among the ticks on lizards in China have been poorly understood. METHODS: In this study, 30 ticks, collected from the multi-ocellated racerunner (Eremias multiocellata) lizard in the Tarim Basin and adjacent Yanqi Basin of the Xinjiang Uygur Autonomous Region in China, were identified by morphological observation and confirmed by DNA-based techniques. The mitochondrially encoded 12S rRNA, 16S rRNA, and COI gene fragments of ticks were amplified and sequenced. To understand the genetic polymorphisms, 47 ticks collected from hedgehogs and 1 from brushwood in the Tarim Basin were also included. Species identification was based on both morphological and molecular characters. The median-joining network approach was used to evaluate the intraspecific genealogies of the ticks and their relatedness with the geographical origin or hosts. RESULTS: The sequence similarity analysis confirmed that the 30 ticks belong to three genera and three species including 11 individuals of Hyalomma asiaticum, 3 of Rhipicephalus turanicus, and 16 of Haemaphysalis sulcata. A network approach revealed paraphyletic populations of R. turanicus and Hy. asiaticum at the intraspecies level regarding geographical origin and low host specificity. For R. turanicus and Hy. asiaticum, common ancestry was observed between COI sequences from lizards and other sequence types from different hosts and countries. CONCLUSIONS: To our knowledge, our study is the first to conduct a molecular survey of ticks from lizards in the arid regions of Xinjiang, China. Eremias multiocellata is an atypical host of the three tick species. Notably, two species of ticks, Hy. asiaticum and R. turanicus, have been collected and identified from lizards in China for the first time. Star-like networks suggest both of them might have experienced recent population expansion. The discoveries are closely related to the geographical environments in Xinjiang and will provide information for the control of ticks and tick-borne pathogens in Northwest China.

Next-generation sequencing yields a complete mitochondrial genome of the Tartar sand boa (Eryx tataricus) from Junggar Basin.

The viviparous Tartar sand boa, Eryx tataricus, is a widespread species in arid Central Asia. A complete mitochondrial genome of one individual from Junggar Basin in Northwest China was determined by next-generation sequencing. The mitogenome is 17,537 bp in size, comprising 2 ribosomal RNA genes, 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 control regions. The order and structure of the genes are similar to those of other Henophidia snakes. Phylogenetic analysis based on 13 concatenated PCGs recovered the monophyly of Boidae and indicated that E. tataricus is closely related to Boa constrictor plus Eunectes notaeus.

Next-generation sequencing yields a complete mitochondrial genome of the Asian Glass Lizard (Dopasia gracillis) from the Yungui Plateau in Southwest China.

The Asian Glass Lizard, Dopasia gracillis, has wide distribution in North India, Nepal, South China, and Indochina. In this study, a complete mitochondrial genome of D. gracillis from the Yungui Plateau in Southwest China was determined by next-generation sequencing. Similar to the typical mtDNA of vertebrates, the mitogenome was 17,133 bp in length and comprised the standard set of 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and 1 control region. The concatenated PCGs were used to conduct Bayesian phylogenetic analyses together with mitogenome data of Anguidae and related taxa in GenBank. The resulting phylogenetic tree confirmed the monophyly of Anguidae and Aguinae as well as D. gracillis, respectively. The mitogenome reported here will contribute to the examination of phylogeographic structure for D. gracillis and understanding of mitochondrial DNA evolution in Anguidae.