Intraspecific phylogeography of Lacerta vivipara and the evolution of viviparity

The lacertid lizard Lacerta vivipara is one of the few squamate species with two reproductive modes. We present the intraspecific phylogeny obtained from neighbor-joining and maximum-parsimony analyses of the mtDNA cytochrome b sequences for 15 individuals from Slovenian oviparous populations, 34 individuals from western oviparous populations of southern France and northern Spain, 92 specimens from European and Russian viviparous populations, and 3 specimens of the viviparous subspecies L. v. pannonica. The phylogeny indicates that the evolutionary transition from oviparity to viviparity probably occurred once in L. vivipara. The western oviparous group from Spain and southern France is phylogenetically most closely related to the viviparous clade. However, the biarmed W chromosome characterizing the western viviparous populations is an apomorphic character, whereas the uniarmed W chromosome, existing both in the western oviparous populations and in the geographically distant eastern viviparous populations, is a plesiomorphic character. This suggests an eastern origin of viviparity. Various estimates suggest that the oviparous and viviparous clades of L. vivipara split during the Pleistocene. Our results are discussed in the framework of general evolutionary models: the concept of an oviparity-viviparity continuum in squamates, the cold climate model of selection for viviparity in squamates, and the contraction-expansion of ranges in the Pleistocene resulting in allopatric differentiation.     

A phylogenetic analysis of variation in reproductive mode within an Australian lizard (Saiphos equalis, Scincidae)

Saiphos equalis , a semi-fossorial scincid lizard from south-eastern Australia, is one of only three reptile species world-wide that are known to display geographic variation in reproductive mode. Uniquely, Saiphos equalis includes populations with three reproductive modes: oviparous with long (15-day) incubation periods; oviparous with short (5-day) incubation periods; and viviparous (0-day incubation periods). No Saiphos populations show 'normal' scincid oviparity (> 30-day incubation period). We used mitochondrial nucleotide sequences ( ND2 and cytochrome b ) to reconstruct relationships among populations from throughout the species' distribution in New South Wales, Australia. Under the phylogenetic species concept, phylogenetic analyses are consistent with the oviparous and viviparous populations of S. equalis being conspecific. Phylogenetic analyses suggest that the long incubation period oviparous lineage is the sister group to all other populations; and that the viviparous populations belong to a cluster of weakly supported clades basal to the short-incubation-period oviparous clade. These clades correspond to variation in reproductive mode and geographic location.     

Intraspecific Phylogeography of Lacerta vivipara and the Evolution of Viviparity

The lacertid lizard Lacerta vivipara is one of the few squamate species with two reproductive modes. We present the intraspecific phylogeny obtained from neighbor-joining and maximum-parsimony analyses of the mtDNA cytochrome b sequences for 15 individuals from Slovenian oviparous populations, 34 individuals from western oviparous populations of southern France and northern Spain, 92 specimens from European and Russian viviparous populations, and 3 specimens of the viviparous subspecies L. v. pannonica. The phylogeny indicates that the evolutionary transition from oviparity to viviparity probably occurred once in L. vivipara. The western oviparous group from Spain and southern France is phylogenetically most closely related to the viviparous clade. However, the biarmed W chromosome characterizing the western viviparous populations is an apomorphic character, whereas the uniarmed W chromosome, existing both in the western oviparous populations and in the geographically distant eastern viviparous populations, is a plesiomorphic character. This suggests an eastern origin of viviparity. Various estimates suggest that the oviparous and viviparous clades of L. vivipara split during the Pleistocene. Our results are discussed in the framework of general evolutionary models: the concept of an oviparity–viviparity continuum in squamates, the cold climate model of selection for viviparity in squamates, and the contraction–expansion of ranges in the Pleistocene resulting in allopatric differentiation.     

Evolutionary and biogeographical implications of the karyological variations in the oviparous and viviparous forms of the lizard Lacerta (Zootoca) vivipara

The lizard Lacerta ( Zootoca ) vivipara has two modes of reproduction and is variable karyologically. We describe its karyological variation from literature data and from new data on two viviparous populations from France, on two oviparous populations from the Pyrenees in south-western France and on three oviparous populations recently discovered in Slovenia. Males have 36 chromosomes, whereas females have only 35 chromosomes in all viviparous populations and in the Pyrenean oviparous populations. This karyotype has been interpreted to result from a fusion of an ancestral sexual W chromosome with an autosome from the Zl or from the Z2 pair. The karyotype formula is 32 autosomes + ZIZ2W for the female and 32 autosomes + Z1Z1Z2Z2 for the male. The karyotype of the Slovenian oviparous populations, 34 autosomes + ZW in the male and 34 autosomes +ZW in the female, represents an evolutionary stage that preceded the chromosomal fusion. There is minor karyological variation, mainly concerning the W and Z2 chromosomes, within the Pyrenean oviparous populations. This parallels the geographic variation of the W-linked alleles of the MPI enzyme and suggests that allopatric differentiation of these oviparous populations might have occurred in the vicinity of the Pyrenees during the Pleistocene.
The viviparous populations from western Europe carry a metacentric W chromosome, whereas oviparous populations from south-western Europe and eastern viviparous populations both show an acrocentrie, or a subtelocentrie. W chromosome. This suggests that the acrocentric-subtelocentric W is a primitive character and that viviparity probably arose in an eastern lineage of the species.     

Further evidence of the existence of oviparous populations of Lacerta (Zootoca) vivipara in the NW of the Balkan Peninsula

The lizard Lacerta (Zootoca) vivipara, which is viviparous in the greatest part of its distribution range, has however some oviparous populations on the southern margin of its range. The present study aimed at determining the reproductive mode and the ATA (aspartate transaminase) enzyme characteristics of four populations in Slovenia and one population in Croatia. The Slovenian females studied here presented an oviparous reproductive mode which strongly resembled those observed in the oviparous populations of south-western France and north-western Spain. Our electrophoresis analyses revealed the existence of two distinct alleles, ATA-150 and ATA-200, in the oviparous populations of Slovenia. These alleles were identical to those observed in the French and Spanish oviparous group and were distinct from the allele ATA-100 characterizing the viviparous populations that we had previously studied. Although we did not study the reproductive mode of Croatian females, the allele ATA-200 observed in one population of Croatia strongly suggested that this population might also be oviparous.     

Fundamentals of viviparity: comparison of seasonal changes in the uterine epithelium of oviparous and viviparous Lerista bougainvillii (Squamata: Scincidae)

Distinct differences in epithelial response between oviparous and viviparous species of skinks led us to investigate morphological differences in the uterus of a species that exhibits bi-modal reproduction and that may indicate specialities for the different requirements of viviparity and oviparity. The uteri of females from oviparous and viviparous populations of the Australian scincid lizard, Lerista bougainvillii, are described in detail to determine whether the occurrence of uterodomes and the plasma membrane transformation, found in other viviparous species but not oviparous species, are indeed features characteristic of viviparity. Oviductal tissue was dissected at three different stages of reproduction from lizards from both populations: 1) vitellogenic, 2) gravid or pregnant, and 3) non-reproductive or quiescent. Tissue was observed using both scanning and transmission electron microscopy. Lerista bougainvillii has a simple placental morphology with simple squamous epithelium. In contrast to mammals and other viviparous skinks, L. bougainvillii does not undergo a plasma membrane transformation, but early signs of placentation in viviparous individuals are indicated by changes in the uterine surface that occur largely after embryonic stage 30. There are no obvious cellular differences between the uteri of oviparous and viviparous L. bougainvillii at the non-reproductive and vitellogenic phase of the reproductive cycle but throughout gestation/gravidity, the cellular differences that could be related to the changing functional requirements with the retention of the viviparous embryo, became apparent. A plasma membrane transformation with ensuing uterodome formation does not occur, which suggests that these more sophisticated changes are a feature of advanced placental development in reptiles.     

MITOCHONDRIAL DNA SEQUENCE-BASED PHYLOGENY AND THE EVOLUTION OF VIVIPARITY IN THE SCELOPORUS SCALARIS GROUP (REPTILIA,SQUAMATA)

The lizard genus Sceloporus contains both oviparous and viviparous species. The scalaris complex is the only monophyletic group within the genus that includes both reproductive modes, thus it is particularly well suited for studies of the evolution of viviparity. Approximately 874 nucleotides of mtDNA sequence data, collected from 38 specimens, comprising 25 populations of all five recognized species within the group, were used in a phylogenetic analysis of the origin of viviparity. Viviparity appears to have evolved twice in this group: once in S. goldmani, included in a clade formed by a northern group consisting of S. scalaris, S. chaneyi, and S. goldmani, and one more time in S. bicanthalis, included in the southern group formed by S. bicanthalis and S. aeneus. An oviparous population of S. bicanthalis nested within that viviparous clade, indicates that reversal from viviparity to oviparity may be possible. Degree of sequence divergence among several S. bicanthalis individuals pertaining to a population in which both parity modes occur, was no larger between oviparous and viviparous lizards than among viviparous lizards. This suggests that this population is a single species, and it may represent a transition from oviparity to viviparity or vice-versa.     

Comparison of the cold hardiness capacities of the oviparous and viviparous forms of Lacerta vivipara

The lizard Lacerta vivipara has allopatric oviparous and viviparous populations. The cold hardiness strategy of L. vivipara has previously been studied in viviparous populations, but never in oviparous ones. The present study reveals that both the oviparous and viviparous individuals of this species are able to survive in a supercooled state at -3 degrees C for at least one week when kept on dry substrates. The mean crystallisation temperatures of the body, around -4 degrees C on dry substrata and -2 degrees C on wet substrata, do not differ between oviparous and viviparous individuals. All the individuals are able to tolerate up to 48-50% of their body fluid converted into ice, but only viviparous individuals were able to stabilize their body ice content at 48%, and hence were able to survive even when frozen at -3 degrees C for times of up 24 hours. Ice contents higher than 51% have been constantly found lethal for oviparous individuals. This suggests that, in L. vivipara, the evolution towards a higher degree of freezing tolerance could parallel the evolution of the viviparous reproductive mode, a feature believed to be strongly selected under cold climatic conditions. This is the first report, among reptiles, of an intraspecific variation regarding the freeze tolerance capacities.     

The evolution of viviparity within the Australian scincid lizardLerista bougainvillii

The semi-fossorial scincid lizard, Lerista bougainvillii , is oviparous throughout its extensive range in south-eastern mainland Australia. However, two widely separated (by approximately 1000 km) island populations are viviparous; in these populations the eggshell is lost and females retain their offspring in utero until embryogenesis complete. One mainland population in south-eastern Victoria shows an intermediate condition, in which the eggshells are incomplete and uterine embryogenesis is prolonged.
Morphological and electrophoretic analyses confirm a high degree of morphological and genetic similarity between populations (i.e. there is no evidence for the presence of more than one species), and phenetic analyses of these data show that each of the two disjunct viviparous populations more closely resembles adjacent oviparous populations than the other viviparous group. Hence, we infer that viviparity may have arisen twice within L. bougainvillii , in both cases on offshore islands with a cold climate.     

Biogeography of Lacerta (Zootoca) vivipara: reproductive mode and enzyme phenotypes in Bulgaria

Lacerta (Zootoca) vivipara , has allopatric oviparous and viviparous populations viviparity is observed from central France and the British Isles to Scandinavia and Russia, while oviparity is restricted to northern Spain and southwestern France, i e the extreme southwestern part of the range Recent observations in the Rila, Balkan, Vitocha, Pirin and Rhodopes mountains indicate that Bulgarian populations of Lacerta (Zootoca) vivipara are indeed viviparous The electrophoretic study of allozymes and the estimation of genetic distances indicate that viviparous lizards from northwest and central France are more closely related to those of Bulgaria, than to the oviparous lizards of southwest France and northwest Spain Variations in reproductive mode and allozymes are not directly related to geographic distances between populations, nor to their latitude populations located at the southwest limit of distribution are oviparous and exhibit alleles ATA^-150 or ATA^-200, whereas, at a comparable latitude, the Bulgarian populations are viviparous and exhibit allele ATA^-100 characteristic of other distant viviparous populations These findings underline the orginality of the oviparous southwestern populations They do not contradict our previous biogeographic scenario     

Refuge areas and suture zones in the Pyrenean and Cantabrian regions: geographic variation of the female MPI sex-linked alleles among oviparous populations of the lizard Lacerta (Zootoca) vivipara

The viviparous lizard Lacerta (Zootoca) vivipara exhibits several alleles of the mannose-6-phosphate isomerase (MPI) enzyme that are carried exclusively on the female W sex-chromosome. Previous studies showed that both the oviparous and viviparous forms of L. (Zootoca) vivipara have these female sex-linked alleles. We document the existence of geographic variation of these alleles among the oviparous populations of southwestern France and northwestern Spain. Two oviparous subgroups were identified: all females from the eastern and central Pyrenees and most females from Aquitaine and from the northern slope of the western Pyrenees exhibited the fast migrating alleles MPI¹¹⁰ or MPI¹²⁰, whereas all females from the Cantabric mountains, Spanish Basque country, and from the south slope of the western Pyrenees exhibited the slow migrating allele MPI⁹⁰. Populations with both fast and slow migrating alleles occurred at some stations in the upper Ossau valley (western Pyrenees) and also at a lowland station of south Aquitaine. The hypothesis that several oviparous forms could have retreated to different places of the Pyreneo-Iberian refugia during the Quaternary glaciations could explain the conservation or the evolution of the polymorphism of the MPI alleles, and that is consistent with the phylogeographic scenario previously proposed to account for the reproductive and cytogenetical variation observed in this species.     

Extraembryonic membrane development in a reproductively bimodal lizard, Lacerta (Zootoca) vivipara

Reproductive mode has been remarkably labile among squamate reptiles and the evolutionary transition from oviparity to viviparity commonly has been accompanied by a shift in the pattern of embryonic nutrition. Structural specializations for placental transfer of nutrients during intrauterine gestation are highly diverse and many features of the extraembryonic membranes of viviparous species differ markedly from those of oviparous species. However, because of a high degree of evolutionary divergence between the species used for comparisons it is likely that the observed differences arose secondarily to the evolution of viviparity. We studied development of the extraembryonic membranes and placentation in the reproductively bimodal lizard Lacerta vivipara because the influence of reproductive mode on the structural/functional relationship between mothers and embryos can best be understood by studying the most recent evolutionary events. Lecithotrophic viviparity has evolved recently within this species and, although populations with different reproductive modes are allopatric, oviparous and viviparous forms interbreed in the laboratory and share many life history characteristics. In contrast to prior comparisons between oviparous and viviparous species, we found no differences in ontogeny or structure of the extraembryonic membranes between populations with different reproductive modes within L. vivipara. However, we did confirm conclusions from previous studies that the tertiary envelope of the egg, the eggshell, is much reduced in the viviparous population. These conclusions support a widely accepted model for the evolution of squamate placentation. We also found support for work published nearly 80 years ago that the pattern of development of the yolk sac of L. vivipara is unusual and that a function of a unique structure of squamate development, the yolk cleft, is hematopoiesis. The structure of the yolk sac splanchnopleure of L. vivipara is inconsistent with a commonly accepted model for amniote yolk sac function and we suggest that a long standing hypothesis that cells from the yolk cleft participate in yolk digestion requires further study.     

Variation of Reproductive Traits and Female Body Size in the Most Widely-Ranging Terrestrial Reptile: Testing the Effects of Reproductive Mode, Lineage, and Climate

The European common lizard, Zootoca vivipara, is the most widespread terrestrial reptile in the world. It occupies almost the entire Northern Eurasia and includes four viviparous and two oviparous lineages. We analysed how female snout-vent length (SVL), clutch size (CS), hatchling mass (HM), and relative clutch mass (RCM) is associated with the reproductive mode and climate throughout the species range and across the evolutionary lineages within Z. vivipara. The studied variables were scored for 1,280 females and over 3,000 hatchlings from 44 geographically distinct study samples. Across the species range, SVL of reproductive females tends to decrease in less continental climates, whereas CS corrected for female SVL and RCM tend to decrease in climates with cool summer. Both relationships are likely to indicate direct phenotypic responses to climate. For viviparous lineages, the pattern of co-variation between female SVL, CS and HM among populations is similar to that between individual females within populations. Consistent with the hypothesis that female reproductive output is constrained by her body volume, the oviparous clade with shortest retention of eggs in utero showed highest HM, the oviparous clade with longer egg retention showed lower HM, and clades with the longest egg retention (viviparous forms) had lowest HM. Viviparous populations exhibited distinctly lower HM than the other European lacertids of similar female SVL, many of them also displaying unusually high RCM. This pattern is consistent with Winkler and Wallin’s model predicting a negative evolutionary link between the total reproductive investment and allocation to individual offspring.     

Length of activity season drives geographic variation in body size of a widely distributed lizard

Understanding the factors that drive geographic variation in life history is an important challenge in evolutionary ecology. Here, we analyze what predicts geographic variation in life‐history traits of the common lizard, Zootoca vivipara, which has the globally largest distribution range of all terrestrial reptile species. Variation in body size was predicted by differences in the length of activity season, while we found no effects of environmental temperature per se. Females experiencing relatively short activity season mature at a larger size and remain larger on average than females in populations with relatively long activity seasons. Interpopulation variation in fecundity was largely explained by mean body size of females and reproductive mode, with viviparous populations having larger clutch size than oviparous populations. Finally, body size‐fecundity relationship differs between viviparous and oviparous populations, with relatively lower reproductive investment for a given body size in oviparous populations. While the phylogenetic signal was weak overall, the patterns of variation showed spatial effects, perhaps reflecting genetic divergence or geographic variation in additional biotic and abiotic factors. Our findings emphasize that time constraints imposed by the environment rather than ambient temperature play a major role in shaping life histories in the common lizard. This might be attributed to the fact that lizards can attain their preferred body temperature via behavioral thermoregulation across different thermal environments. Length of activity season, defining the maximum time available for lizards to maintain optimal performance, is thus the main environmental factor constraining growth rate and annual rates of mortality. Our results suggest that this factor may partly explain variation in the extent to which different taxa follow ecogeographic rules.     

Lerista bougainvillii,a case study for the evolution of viviparity in reptiles

Many factors, both environmental and biotic, have been suggested to facilitate or hinder the evolution of viviparity (live-bearing) in reptiles. Viviparity has evolved recently within the Australian scincid lizard Lerista bougainvillii and the species includes oviparous, viviparous, and reproductively intermediate (with prolonged egg retention) populations; thus, it offers an exceptional opportunity to evaluate the validity of these hypotheses. We carried out such tests by (i) comparing environmental conditions over the geographic ranges occupied by oviparous, viviparous, and intermediate populations (to identify possible selective forces for the evolution of viviparity), and (ii) comparing morphological, reproductive and ecological traits of L. bougainvillii with those of other sympatric scincid species (to identify traits that may have predisposed this taxon to the evolution of viviparity). The areas occupied by viviparous L. bougainvillii are significantly colder than those occupied by both their intermediate and oviparous conspecifics, in accord with the “cold-climate” hypothesis for reptilian viviparity. Rainfall is similar over the ranges of the three forms. Climatic unpredictability (as assessed by the magnitude of year-to-year thermal variation) is lower for viviparous animals, in contradiction to published speculations. Comparison with 31 sympatric scincid species showed that L. bougainvillii is not atypical for most of the traits we measured (e.g., body size, clutch size, thermal preferenda and tolerances). However, oviparous L. bougainvillii do display several traits that have been suggested to facilitate the evolution of viviparity. For example, pregnancy does not reduce locomotor ability of females; the lizards are semi-fossorial; even the oviparous females produce only a single clutch of eggs per year; and they ovulate relatively late in summer, so that the time available for incubation is limited.     

Mitochondrial and Nuclear DNA Survey of Zootoca vivipara across the Eastern Italian Alps: Evolutionary Relationships,Historical Demography and Conservation Implications

The European common lizard Zootoca vivipara exhibits reproductive bimodality, with populations being either viviparous or oviparous. In the central-eastern Italian Alps oviparous populations (Z. v. carniolica) and viviparous populations (Z. v. vivipara) partly overlap geographically. Studying the evolutionary relationship between these taxa presents an interesting opportunity to gain insight into the evolution of this trait. We aim to: i) test whether Z. v. carniolica, which is endangered, constitutes an ESU (Evolutionary Significant Unity); ii) infer mtDNA divergence time between the Z. v. carniolica clade and all the other Z. vivipara subspecies with the aid of an external calibration point; and iii) describe the phylogeographical and demographic scenarios in the area. To do so we sequenced about 200 individuals for mitochondrial variation; 64 of them were also analysed for three nuclear genes. Furthermore, we analysed the same nuclear markers in 17 individuals from the other oviparous subspecies Z. v. louislantzi and 11 individuals of Z. v. vivipara from widespread geographical origins. The mtDNA and nDNA loci that we examined supported the monophyly of Z. v. carniolica. The mtDNA-based estimate of divergence time between Z. v. carniolica and all the other subspecies indicated a separation at 4.5 Mya (95% CI 6.1–2.6), with about 5% of sequence divergence. Considering that Z. v. carniolica harbours higher genetic diversity, while Z. v. vivipara from central-eastern Alps shows a signature of recent population and spatial expansion, we argue that Z. v. carniolica represents a distinct evolutionary unit, with a presumably long-term evolutionary history of separation. Z. v. carniolica populations, occurring at higher latitudes and altitudes than insofar supposed, live in peat bogs, a seriously threatened habitat: taking into account also its evolutionary distinctness, specific conservation measures should be considered.     

Experimental evidence of early costs of reproduction in conspecific viviparous and oviparous lizards

Reproduction entails costs, and disentangling the relative importance of each stage of the reproductive cycle may be important to assess the costs and benefits of different reproductive strategies. We studied the early costs of reproduction in oviparous and viviparous lizard females of the bimodal reproductive species Zootoca vivipara. Egg retention time in oviparous females is approximately one-third of the time in viviparous females. We compared the vitellogenesis and egg retention stages that are common to both reproductive modes. Precisely, we monitored the thermoregulatory behaviour, the weight gain and the immunocompetence of the females. Moreover, we injected an antigen in half of the females (immune challenge) to study the trade-offs between reproductive mode and immune performance and between different components of the immune system. Finally, we experimentally induced parturition in viviparous females at the time of egg laying in oviparous females. Oviparous and viviparous females did not show strong differences in response to the immune challenge. However, viviparous females spent more time thermoregulating while partially hidden and gained more weight than oviparous females. The greater weight gain indicates that the initial period of egg retention is less costly for viviparous than for oviparous females or that viviparous females are able to save and accumulate energy at this period. This energy may be used by viviparous females to cope with the subsequent costs of the last two-third of the gestation. Such an ability to compensate the higher costs of a longer egg retention period may account for the frequent evolution of viviparity in squamate reptiles.     

Karyological investigations of populations of the lizard Zootoca vivipara (Juaquin, 1787) from the Baltic Sea basin (western region of Russia)

Populations of the lizard Zootoca vivipara (Juaquin, 1787) (Lacertidae, Sauria) from Western and Central Europe are poorly differentiated in morphology. However, they differ from each other by several karyotype markers and their reproductive mode (oviparous and viviparous types). Here we report for the first time the results of combined investigations of Z. vivipara from Eastern Europe (the Baltic Sea basin of Western Russia). The karyotype and reproductive mode of the specimens from four previously not examined populations in Kaliningrad region have been studied. It has been shown that these lizards have viviparous mode of reproduction and their karyotypes vary on morphology of W-sex chromosomes. Females possess 2n = 35 : 32 A + Z1Z2W, with W--SV and 2n = 35 : 32 A + Z1Z2W, with W--A/ST. We have identified these populations karyologically and have first revealed that they belong to two different viviparous forms of nominative subspecies Z. v. vivipara, namely to its "western" and "Russian (eastern)" forms. These data have confirmed that the Baltic Sea basin in Eastern Europe is a zone of secondary contact of these chromosomal forms of Z. v. vivipara. The results obtained do not now allow us to speak about their sympatry, parapatry and (or) hybridization events. The results suggest that during postglasial time "western" form might colonize this region from western Europe whereas "Russian (eastern)" form might do it from Southern-Eastern Europe. It has been recommended to save the biodiversity of Z. vivipara in Kaliningrad region, namely the peripheral and very rare on vast territory of Russia "western" form and unique for western region of Russia "Russian (eastern)" form.     

A phylogeny of Chinese species in the genus Phrynocephalus (Agamidae) inferred from mitochondrial DNA sequences

We investigated the phylogenetic relationships among most Chinese species of lizards in the genus Phrynocephalus (118 individuals collected from 56 populations of 14 well-defined species and several unidentified specimens) using four mitochondrial gene fragments (12S rRNA, 16S rRNA, cytochrome b, and ND4-tRNA(LEU)). The partition-homogeneity tests indicated that the combined dataset was homogeneous, and maximum-parsimony (MP), neighbor-joining (NJ), maximum-likelihood (ML) and Bayesian (BI) analyses were performed on this combined dataset (49 haplotypes including outgroups for 2058bp in total). The maximum-parsimony analysis resulted in 24 equally parsimonious trees, and their strict consensus tree shows that there are two major clades representing the Chinese Phrynocephalus species: the viviparous group (Clade A) and the oviparous group (Clade B). The trees derived from Bayesian, ML, and NJ analyses were topologically identical to the MP analysis except for the position of P. mystaceus. All analyses left the nodes for the oviparous group, the most basal clade within the oviparous group, and P. mystaceus unresolved. The phylogenies further suggest that the monophyly of the viviparous species may have resulted from vicariance, while recent dispersal may have been important in generating the pattern of variation among the oviparous species.