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.

Historical Demography of an Endangered Salamander, Ranodon sibiricus (Amphibia, Urodela, Hynobiidae): A Reassessment.

Chinese populations of the endangered Siberian salamander Ranodon sibiricus are reported to have diverged only about 120 years ago, and to have the lowest genetic diversity of any amphibian. However, these conclusions require verification, as the main range of the species is in Kazakhstan. Moreover, the generation time used for estimating divergence time has a weak ground. In order to clarify these problems, we investigated the molecular phylogenetic relationship and historical demography of the species covering its whole distribution range using the mitochondrial DNA region reported for Chinese population (1072 bp sequences of the control region), while conducting skeletochronological analysis to estimate accurate generation time. As a result, the range expansion was estimated at 88,000-50,000 YA, based on the generation time of 6-10 years. Degree of intraspecific genetic differentiation is actually very small, but, as a single species, is not so small as had been reported for Chinese population alone.

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.

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.

A molecular phylogenetic hypothesis for the Asian agamid lizard genus Phrynocephalus reveals discrete biogeographic clades implicated by plate tectonics.

Phylogenetic relationships of the agamid lizard genus Phrynocephalus are described in the context of plate tectonics. A near comprehensive taxon sampling reports three data sets: (1) mitochondrial DNA from ND1 to COI (3' end of ND1, tRNAGln, tRNAIle, tRNAMet, ND2, tRNATrp, tRNAAla, tRNAAsn, tRNACys, tRNATyr, and the 5' end of COI) with 1761 aligned positional sites (1595 included, 839 informative), (2) nuclear RAG-1 DNA with 2760 aligned positional sites (342 informative), and (3) 25 informative allozyme loci with 213 alleles (107 informative when coded as presence/absence). It is hypothesized that Phrynocephalus phyletic patterns and speciation reflect fault lines of ancient plates now in Asia rejuvenated by the more recent Indian and Arabian plate collisions. Molecular estimates of lineage splits are highly congruent with geologic dates from the literature.  A southern origin for the genus in Southwest Asia is resolved in phylogenetic estimates and a northern origin is statistically rejected. On the basis of monophyly and molecular evidence several taxa previously recognized as subspecies are recognized as species: P. hongyuanensis, P. sogdianus, and P. strauchi as "Current Status"; Phrynocephalus bannikovi, Phrynocephalus longicaudatus, Phrynocephalus turcomanus, and Phrynocephalus vindumi are formally "New Status". Phylogenetic evaluation indicates a soft substrate habitat of sand for the shared ancestor of modern Phrynocephalus. Size diversity maximally overlaps in the Caspian Basin and northwestern Iranian Plateau. The greatest species numbers of six in sympatry and regional allopatry are found in the southern Caspian Basin and southern Helmand Basin, both from numerous phylogenetic lineages in close proximity attributed to tectonic induced events.

Serous cutaneous glands in the South American horned frog Ceratophrys ornata (Leptodactylformes, Chthonobatrachia, Ceratophrydae): ultrastructural expression of poison biosynthesis and maturation.

Serous cutaneous glands are described in newly metamorphosed and juvenile specimens of the horned frog Ceratophrys ornata using light and transmission electron microscopy. We report patterns of biosynthesis and maturation of the specific product of the gland secretory unit. The syncytial, secretory compartment possesses a complex of endoplasmic reticulum (predominantly smooth endoplasmic reticulum after metamorphosis) and Golgi stacks. The serous product is weak in density and is contained in vesicles involved in repeating merging processes. During this maturation activity, secondary lysosomes are observed, which derive from autophagic processes (crinophagy) involving the secretory materials. Ceratophrys ornata, a species representative of the type genus of the family Ceratophrydae, belongs to the heterogeneous group of anurans that, possibly as the result of convergence, all produce cutaneous poisons consisting of vesicles or faint density granules.

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.

First record of Limnatispaluda (Hirudinida, Arhynchobdellida, Praobdellidae) from Kazakhstan, with comments on genetic diversity of Limnatis leeches.

BACKGROUND: Sawyer (1986) included three species in the nasal leech genus LimnatisMoquin-Tandon 1827: Limnatisnilotica (Savigny 1822), LimnatisbacescuiManoleli 1972 and Limnatispaluda (Tennent 1859). The first and last species have mainly been identified in Middle Eastern countries (e.g. Kinzelbach and Rückert 1985). The second species has been identified only in Romania Dobruja (Manoleli 1972). Although Limnatis leeches are well known species of endoparasitic leeches, Limnatisnilotica was recorded only once in Kazakhstan (Lukin 1976). NEW INFORMATION: Specimens of the genus Limnatis from Almaty Province, Kazakhstan are identified as Limnatispaluda. This is the first record of Limnatispaluda from Kazakhstan. Mitochondrial COI and 12S data demonstrated that the present specimens are genetically close to an Israeli specimen identified as Limnatisnilotica. In addition, molecular data suggest that some Limnatis specimens whose DNA sequences have been reported were misidentified. According to the observed phylogenetic relationships, the taxonomic status of the known Limnatis species should be revisited.

Evolution of mitochondrial relationships and biogeography of Palearctic green toads (Bufo viridis subgroup) with insights in their genomic plasticity.

Taxa involving three bisexually reproducing ploidy levels make green toads a unique amphibian system. We put a cytogenetic dataset from Central Asia in a molecular framework and apply phylogenetic and demographic methods to data from the entire Palearctic range. We study the mitochondrial relationships of diploids to infer their phylogeography and the maternal ancestry of polyploids. Control regions (and tRNAs between ND1 and ND2 in representatives) characterize a deeply branched assemblage of twelve haplotype groups, diverged since the Lower Miocene. Polyploidy has evolved several times: Central Asian tetraploids (B. oblongus, B. pewzowi) have at least two maternal origins. Intriguingly, the mitochondrial ancestor of morphologically distinctive, sexually reproducing triploid taxa (B. pseudoraddei) from Karakoram and Hindukush represents a different lineage. We report another potential case of bisexual triploid toads (B. zugmayeri). Identical d-loops in diploids and tetraploids from Iran and Turkmenistan, which differ in morphology, karyotypes and calls, suggest multiple origins and retained polymorphism and/or hybridization. A similar system involves diploids, triploids and tetraploids from Kyrgyzstan and Kazakhstan where green toads exemplify vertebrate genomic plasticity. A new form from Sicily and its African sister species (B. boulengeri) allow internal calibration and divergence time estimates for major clades. The subgroup may have originated in Eurasia rather than Africa since the earliest diverged lineages (B. latastii, B. surdus) and earliest fossils occur in Asia. We delineate ranges, contact and hybrid zones. Phylogeography, including one of the first non-avian datasets from Central Asian high mountains, reflects Quaternary climate and glaciation.