Phylogeny and divergence times of some racerunner lizards (Lacertidae: Eremias) inferred from mitochondrial 16S rRNA gene segments

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

The Silk roads: phylogeography of Central Asian dice snakes (Serpentes: Natricidae) shaped by rivers in deserts and mountain valleys

Influenced by rapid changes in climate and landscape features since the Miocene, widely distributed species provide suitable models to study the environmental impact on their evolution and current genetic diversity. The dice snake Natrix tessellata, widely distributed in the Western Palearctic is one such species. We aimed to resolve a detailed phylogeography of N. tessellata with a focus on the Central Asian clade with 4 and the Anatolia clade with 3 mitochondrial lineages, trace their origin, and correlate the environmental changes that affected their distribution through time. The expected time of divergence of both clades began at 3.7 Mya in the Pliocene, reaching lineage differentiation approximately 1 million years later. The genetic diversity in both clades is rich, suggesting different ancestral areas, glacial refugia, demographic changes, and colonization routes. The Caspian lineage is the most widespread lineage in Central Asia, distributed around the Caspian Sea and reaching the foothills of the Hindu Kush Mountains in Afghanistan, and Eastern European lowlands in the west. Its distribution is limited by deserts, mountains, and cold steppe environments. Similarly, Kazakhstan and Uzbekistan lineages followed the Amu Darya and the Syr Darya water systems in Central Asia, with ranges delimited by the large Kyzylkum and Karakum deserts. On the western side, there are several lineages within the Anatolia clade that converged in the central part of the peninsula with 2 being endemic to Western Asia. The distribution of both main clades was affected by expansion from their Pleistocene glacial refugia around the Caspian Sea and in the valleys of Central Asia as well as by environmental changes, mostly through aridification.

Climate-driven mitochondrial selection in lacertid lizards

The mitochondrion, which is an intracellular organelle responsible for most of the energy-producing pathways, can have its genome targeted for climate-driven selection. However, climate-driven mitochondrial selection remains a sparsely studied area in reptiles. Here, we reported the complete mitochondrial genome sequence of a lacertid lizard (Takydromus intermedius) and used mitogenomes from 54 species of lacertid lizards to study their phylogenetic relationships and to identify the mitochondrial genes under positive selection by climate. The length of the complete mitochondrial genome sequence of T. intermedius was 17,713 bp, which was within the range of lengths (17,224-18,943) ever reported for Takydromus species. The arrangement of mitochondrial genes in T. intermedius was the same as in other congeneric species. The 54 lacertid species could be divided into three geographically and climatically different clades. We identified three mitochondrial genes (ATP6, ATP8, and ND3) under positive selection by climate, and found that isothermality, temperature seasonality, precipitation of wettest month, and precipitation seasonality were the most important climatic variables contributing to the gene selection.

Taxonomic diversity of racerunners with descriptions of two new Eremias species (Sauria: Lacertidae) from Central Iran

We provide a diversity assessment of Iranian species of the genus Eremias based on the cytochrome oxidase I mtDNA gene fragment. We analyzed 93 genetic samples from the entire distribution of the Eremias fasciata species complex in Iran and surrounding regions, along with morphological data to support the description of two new species from Central Iran. We hypothesize that the diversification of the Eremias fasciata species complex was largely influenced by the fragmentation of sand massifs in the region. This same hypothesis has been used to explain the high level of endemism among the sand-dwelling species of reptiles along the Iranian Plateau in the same area. The two new species described herein can be distinguished from other congeneric species by their phylogenetic position and a combination of morphological characters. We use these data to discuss the taxonomy of Eremias based on morphology, habitat choice, and genetic data.

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.

Integrative Taxonomy within Eremias multiocellata Complex (Sauria, Lacertidae) from the Western Part of Range: Evidence from Historical DNA

The Kokshaal racerunner, Eremias kokshaaliensis Eremchenko et Panfilov, 1999, together with other central Asian racerunner species, is included in the Eremias multiocellata complex. In the present work, for the first time, the results of the analysis of historical mitochondrial DNA (barcode) are presented and the taxonomic status and preliminary phylogenetic relationships within the complex are specified. We present, for the first time, the results of the molecular analysis using historical DNA recovered from specimens of several species of this complex (paratypes of the Kokshaal racerunner and historical collections of the Kashgar racerunner E. buechneri from Kashgaria) using DNA barcoding.

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.

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.

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.

Genetic diversity of major histocompatibility complex class I genes in Zootoca vivipara

The Major Histocompatibility Complex (MHC), as a family of highly polymorphic genes associated with immunity in the genome of the vertebrate, has become an important indicator for assessing the evolutionary potential of wildlife. In order to better protect Zootoca vivipara in the Greater Khingan Range and Lesser Khingan Range, to understand the genetic structure of Z. vivipara, and to explore the mechanism and phylogenetic relationship of the gene polymorphisms, the MHC molecular marker method was used to analyze Z. vivipara population. Forty-seven alleles were obtained from four populations. The four populations were highly polymorphic, rich in genetic information, and had significant genetic diversity. There were certain inbreeding phenomena. There was a high degree of genetic differentiation among populations, which was caused by genetic drift and natural selection. The sequence undergoes genetic duplication and recombination. The existence of trans-species polymorphism was found in the constructed phylogenetic tree. The present study provides a theoretical basis for species protection of Z. vivipara.

Phylogeography of the red-bellied lizard, Darevskia parvula in Turkey

The groups of red-bellied lizards had a small distribution area in the Pontic zone. The several studies performed on these lizard groups are based on taxonomy and systematics. Although there were several taxonomic or systematic researches on some species of this group, the phylogeographical pattern and species disturbing boundaries of this group is still not clear. In the present study, we aimed to resolve the taxonomic and phylogenetic relationships of the red-bellied lizards in Turkey, based on two combined mitochondrial gene fragments and one protein-coding nuclear gene (rag1). Also, we evaluated ecological niches differentiations among red-bellied lizard groups. The mitochondrial DNA genes were found to be highly polymorphic in this group. One hundred and one variable nucleotide sites were detected on the combined gene sequences. According to phylogenetic trees based on the maximum likelihood (ML) and Bayesian inference (BI), the red-bellied lizards group have three species groups; Darevskia parvula, D. adjarica and unnamed Darevskia sp. (candidate species for Darevskia genus). This situation was supported by high bootstrap and posterior probability values in the trees of mitochondrial DNA gene fragments. However, no genetic variation was detected according to nuclear DNA (rag1) sequence. Because the species groups have no overlaps in terms of their ecological niches, ecological niche modelling (ENM) results revealed differences among the groups of D. parvula, D. adjarica, and unnamed Darevskia sp. Besides, we detected no geographical overlaps among three species groups, since there were geographical isolation zones among the species groups of red-bellied lizard.

Comparative phylogeography of amphibians and reptiles in Algeria suggests common causes for the east-west phylogeographic breaks in the Maghreb

A series of phylogeographic studies in the Maghreb identified a repeated pattern of deep genetic divergence between an eastern (Tunisia) and western (Morocco) lineage for several taxa but lack of sampling in Algeria made it difficult to know if the range limits between the eastern and western lineages were shared among taxa or not. To address this question, we designed a comparative phylogeographic study using 8 reptile and 3 amphibian species with wide distribution in the Maghreb as models. We selected species where previous studies had identified an East-West phylogeographic divide and collected sampled in Algeria to 1) examine whether the simple East-West divergence pattern still holds after filling the sampling gap in Algeria or if more complex diversity patterns emerge; 2) if the E-W pattern still holds, test whether the limits between the E and W clades are shared between species, suggesting that common historical process caused the E-W divergences; 3) if E-W limits are shared between species, use information on the age of the divergence to identify possible geological or climatic events that could have triggered these E-W differentiations. We found that the E-W pattern was generally maintained after additional sampling in Algeria and identified two common disjunction areas, one around the Algeria-Morocco border, the other one in Kabylia (central Algeria), suggesting that common historical mechanisms caused the E-W divergences in the Maghreb. Our estimates for the times to most common recent ancestors to the E and W clades span a wide range between the Messinian salinity crisis and the Plio-Pleistocene limit (except for one older split), suggesting different origins for the initial divergences and subsequent preservation of the E and W lineages in common climatic refugia in the west and the east of the Maghreb.

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.

Extinction vs. Rapid Radiation: The Juxtaposed Evolutionary Histories of Coelotine Spiders Support the Eocene-Oligocene Orogenesis of the Tibetan Plateau

Evolutionary biology has long been concerned with how changing environments affect and drive the spatiotemporal development of organisms. Coelotine spiders (Agelenidae: Coelotinae) are common species in the temperate and subtropical areas of the Northern Hemisphere. Their long evolutionary history and the extremely imbalanced distribution of species richness suggest that Eurasian environments, especially since the Cenozoic, are the drivers of their diversification. We use phylogenetics, molecular dating, ancestral area reconstructions, diversity, and ecological niche analyses to investigate the spatiotemporal evolution of 286 coelotine species from throughout the region. Based on eight genes (6.5 kb) and 2323 de novo DNA sequences, analyses suggest an Eocene South China origin for them. Most extant, widespread species belong to the southern (SCG) or northern (NCG) clades. The origin of coelotine spiders appears to associate with either the Paleocene-Eocene Thermal Maximum or the hot period in early Eocene. Tibetan uplifting events influenced the current diversity patterns of coelotines. The origin of SCG lies outside of the Tibetan Plateau. Uplifting in the southeastern area of the plateau blocked dispersal since the Late Eocene. Continuous orogenesis appears to have created localized vicariant events, which drove rapid radiation in SCG. North-central Tibet is the likely location of origin for NCG and many lineages likely experienced extinction owing to uplifting since early Oligocene. Their evolutionary histories correspond with recent geological evidence that high-elevation orographical features existed in the Tibetan region as early as 40-35 Ma. Our discoveries may be the first empirical evidence that links the evolution of organisms to the Eocene-Oligocene uplifting of the Tibetan Plateau. [Tibet; biogeography; ecology; molecular clock; diversification.].

The reanalysis of biogeography of the Asian tree frog, Rhacophorus (Anura: Rhacophoridae): geographic shifts and climatic change influenced the dispersal process and diversification

Rapid uplifts of the Tibetan Plateau and climate change in Asia are thought to have profoundly modulated the diversification of most of the species distributed throughout Asia. The ranoid tree frog genus Rhacophorus, the largest genus in the Rhacophoridae, is widely distributed in Asia and especially speciose in the areas south and east of the Tibetan Plateau. Here, we infer phylogenetic relationships among species and estimate divergence times, asking whether the spatiotemporal characteristics of diversification within Rhacophorus were related to rapid uplifts of the Tibetan Plateau and concomitant climate change. Phylogenetic analysis recovered distinct lineage structures in Rhacophorus, which indicated a clear distribution pattern from Southeast Asia toward East Asia and India. Molecular dating suggests that the first split within the genus date back to the Middle Oligocene (approx. 30 Ma). The Rhacophorus lineage through time (LTT) showed that there were periods of increased speciation rate: 14–12 Ma and 10–4 Ma. In addition, ancestral area reconstructions supported Southeast Asia as the ancestral area of Rhacophorus. According to the results of molecular dating, ancestral area reconstructions and LTT we think the geographic shifts, the staged rapid rises of the Tibetan Plateau with parallel climatic changes and reinforcement of the Asian monsoons (15 Ma, 8 Ma and 4–3 Ma), possibly prompted a burst of diversification in Rhacophorus.

Phylogenetic relationships of Podarcis siculus (Rafinesque-Schmaltz, 1810) and Podarcis tauricus (Pallas, 1814) in Turkey, based on mitochondrial DNA

The Italian wall lizard and the Balkan wall lizard have a series of taxonomic revisions. However, their phylogenetic relationships still remain uncertain in Turkey. In the present study, we have assessed taxonomic relationships, both of Podarcis siculus and Podarcis tauricus through estimation of phylogenetic relationships among 43 and 42 specimens, respectively, using mtDNA (16 S rRNA and cytb) from great main populations in Turkey. The genetic distances among the populations of P. siculus in Turkey were very low and they were ranged from 0.2 to 1.6% in 16 S rRNA while they were ranged from 0.0% to 3.3% in cytb. On the other hand, the p-distances among the populations of P. tauricus were ranged from 0.0 to 0.6% in 16 S rRNA while they were 0.2% cytb in Turkey. Finally, most of the topologically identical trees of phylogenetic analyses and p-distances showed that monophyly was found in extant populations of P. siculus and P. tauricus. The nominate subspecies, P. s. siculus and P. t. tauricus are representatives of these lizards in Turkey.

Phylogenetic relationships of D. rudis (Bedriaga, 1886) and D. bithynica(Mehely, 1909) based on microsatellite and mitochondrial DNA in Turkey

The spiny-tailed lizard of the genus Darevskia have a series of taxonomic revisions, but still their phylogenetic relationships remain uncertain. In the present study, we have assessed taxonomic relationships among Darevskia bithynica and Darevskia rudis populations through estimation of phylogenetic relationships among 96 specimens using microsatellite DNA (Du215, Du281, Du323 and Du418 loci) and 53 specimens using mtDNA (16S rRNA and cytb) from main populations in Turkey. Although D. b. bithynica and D. r. mirabilis were separated based on the PCoA analysis at low level from other D. rudis and D. bithynica populations, the distance values of Nei's genetic distance, Nei's unbiased genetic distance, Fst and Linear Fst were not high among taxa in microsatellite DNA. On the other hand, our phylogenetic analyses (NJ, ML, MP and BI) did not separate D. rudis and D. bithynica populations. Finally, most of the topologically identical trees of phylogenetic analyses and microsatellite results showed that the extant populations of D. rudis and D. bithynica were found to be polytomy. Based on our molecular phylogenetic study, D. rudis complex is still ongoing revisions.

The phylogenetic position and diversity of the enigmatic mongrel frog Nothophryne Poynton, 1963 (Amphibia, Anura)

The phylogenetic relationships of the African mongrel frog genus Nothophryne are poorly understood. We provide the first molecular assessment of the phylogenetic position of, and diversity within, this monotypic genus from across its range-the Afromontane regions of Malawi and Mozambique. Our analysis using a two-tiered phylogenetic approach allowed us to place the genus in Pyxicephalidae. Within the family, Nothophryne grouped with Tomopterna, a hypothesis judged significantly better than alternative hypotheses proposed based on morphology. Our analyses of populations across the range of Nothophryne suggest the presence of several cryptic species, at least one species per mountain. Formal recognition of these species is pending but there is a major conservation concern for these narrowly distributed populations in an area impacted by major habitat change. The phylogenetic tree of pyxicephalids is used to examine evolution of life history, ancestral habitat, and biogeography of this group.

The evolution of colour pattern complexity: selection for conspicuousness favours contrasting within-body colour combinations in lizards

Many animals display complex colour patterns that comprise several adjacent, often contrasting colour patches. Combining patches of complementary colours increases the overall conspicuousness of the complex pattern, enhancing signal detection. Therefore, selection for conspicuousness may act not only on the design of single colour patches, but also on their combination. Contrasting long- and short-wavelength colour patches are located on the ventral and lateral surfaces of many lacertid lizards. As the combination of long- and short-wavelength-based colours generates local chromatic contrast, we hypothesized that selection may favour the co-occurrence of lateral and ventral contrasting patches, resulting in complex colour patterns that maximize the overall conspicuousness of the signal. To test this hypothesis, we performed a comparative phylogenetic study using a categorical colour classification based on spectral data and descriptive information on lacertid coloration collected from the literature. Our results demonstrate that conspicuous ventral (long-wavelength-based) and lateral (short-wavelength-based) colour patches co-occur throughout the lacertid phylogeny more often than expected by chance, especially in the subfamily Lacertini. These results suggest that selection promotes the evolution of the complex pattern rather than the acquisition of a single conspicuous colour patch, possibly due to the increased conspicuousness caused by the combination of colours with contrasting spectral properties.

Sequencing and analysis of the whole mitochondrial genome of a variegated racerunner from Taklamakan Desert

The whole mitochondrial genome of a variegated racerunner (Eremias vermiculata) from the Taklamakan Desert was determined using polymerase chain reaction and directly sequenced with a primer walking method. The mitogenome sequence was 19 796 bp in size, containing 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a control region (D-loop), 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 seven Eremias species. Monophyly of the genus Eremias and E. vermiculata was recovered. The mitogenome presented here will contribute to the examination of genetic differentiation for E. vermiculata and understanding of the mitochondrial DNA evolution in Eremias.

The nearly complete mitochondrial genome of the rapid racerunner Eremias velox (Squamata: Lacertidae)

We report the nearly complete mitochondrial genome of the rapid racerunner, Eremias velox (Lacertidae), which is a circular molecule of 18,033 bp in size and consists of 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs and a partial control region (2627 bp). The A + T content of overall base composition of H-strand is 58.0% (T: 27.7%, C: 28.3%, A: 30.3%, G: 13.7%). Some short microsatellite-like repeat regions (polyA and polyT) are scattered in the control region. All the results provide powerful data for further study of the molecular systematics, species identification and conservation genetics in this species and its congenators.

Polymorphic microsatellite loci in the rapid racerunner Eremias velox (Squamata: Lacertidae)

We isolated and characterizated 12 polymorphic microsatellite loci in the rapid racerunner Eremias velox (Squamata: Lacertidae). The loci were screened in 37 E. velox individuals. The number of alleles ranged from 6 to 16. The observed heterozygosity ranged from 0.432 to 0.919, and the expected heterozygosity ranged from 0.685 to 0.902. These microsatellite markers should prove useful for population genetic studies of E. velox and other Eremias species.