Knowing whether triploid hybrids resulting from natural hybridization of parthenogenetic and bisexual species are fertile is crucial for understanding the mechanisms of reticulate evolution in rock lizards. Here, using males of the bisexual diploid rock lizard species Darevskia raddei nairensis and Darevskia valentini and a triploid hybrid male Darevskia unisexualis × Darevskia valentini, we performed karyotyping and comparative immunocytochemistry of chromosome synapsis and investigated the distribution of RAD51 and MLH1 foci in spread spermatocyte nuclei in meiotic prophase I. Three chromosome sets were found to occur in cell nuclei in the D. unisexualis × D. valentini hybrid, two originating from a parthenogenetic D. unisexualis female and one from the D. valentini male. Despite this distorted chromosome synapsis and incomplete double-strand breaks repair in meiotic prophase I, the number of mismatch repair foci in the triploid hybrid was enough to pass through both meiotic divisions. The defects in synapsis and repair did not arrest meiosis or spermatogenesis. Numerous abnormal mature spermatids were observed in the testes of the studied hybrid.

Among vertebrates, true parthenogenesis is known only in reptiles. Parthenogenetic lizards of the genus Darevskia emerged as a result of the hybridization of bisexual parental species. However, uncertainty remains about the mechanisms of the co-existence of these forms. The geographical parthenogenesis hypothesis suggests that unisexual forms can co-exist with their parental species in the “marginal” habitats. Our goal is to investigate the influence of environmental factors on the formation of ecological niches and the distribution of lizards. For this reason, we created models of species distribution and ecological niches to predict the potential geographical distribution of the parthenogenetic and its parental species. We also estimated the realized niches breadth, their overlap, similarities, and shifts in the entire space of predictor variables. We found that the centroids of the niches of the three studied lizards were located in the mountain forests. The “maternal” species D. mixta prefers forest habitats located at high elevations, “paternal” species D. portschinskii commonly occurs in arid and shrub habitats of the lower belt of mountain forests, and D. dahli occupies substantially an intermediate or “marginal” position along environmental gradients relative to that of its parental species. Our results evidence that geographical parthenogenesis partially explains the co-existence of the lizards.

Despite numerous works devoted to hybrid origin of parthenogenesis in reptiles, the causes of hybridization between different species, resulting in the origin of parthenogenetic forms, remain uncertain. Recent studies demonstrate that sexual species considered parental to parthenogenetic rock lizards (Darevskia spp.) avoid interspecific mating in the secondary overlap areas. A specific combination of environmental factors during last glaciation period was critical for ectotherms, which led to a change in their distribution and sex ratio. Biased population structure (e.g., male bias) and limited available distributional range favored the deviation of reproductive behavior when species switched to interspecific mates. To date, at least 7 diploid parthenogenetic species of rock lizards (Darevskia, Lacertidae) originated through interspecific hybridization in the past. The cytogenetic specifics of meiosis, in particular the weak checkpoints of prophase I, may have allowed the formation of hybrid karyotypes in rock lizards. Hybridization and polyploidization are 2 important evolutionary forces in the genus Darevskia. At present, throughout backcrossing between parthenogenetic and parental species, the triploid and tetraploid hybrid individuals appear annually, but no triploid species found among Darevskia spp. on current stage of evolution. The speciation by hybridization with the long-term stage of diploid parthenogenetic species, non-distorted meiosis, together with the high ecological plasticity of Caucasian rock lizards provide us with a new model for considering the pathways and persistence of the evolution of parthenogenesis in vertebrates.

The parthenogenetic Caucasian rock lizard Darevskia armeniaca, like most other parthenogenetic vertebrate species, originated through interspecific hybridization between the closely related sexual Darevskia mixta and Darevskia valentini. Darevskia armeniaca was shown to consist of one widespread allozyme clone and a few rare ones, but notwithstanding the origin of clonal diversity remains unclear. We conduct genomic analysis of D. armeniaca and its parental sexual species using microsatellite and SNP markers to identify the origin of parthenogenetic clonal lineages. Four microsatellite-containing loci were genotyped for 111 specimens of D. armeniaca, 17 D. valentini, and four D. mixta. For these species, a total of 47 alleles were isolated and sequenced. Analysis of the data revealed 13 genotypes or presumptive clones in parthenogenetic D. armeniaca, including one widespread clone, two apparently geographically restricted clones, and ten rare clones. Comparisons of genotype-specific markers in D. armeniaca with those of its parental species revealed three founder-events including a common and two rare clones. All other clones appeared to have originated via post-formation microsatellite mutations in the course of evolutionary history of D. armeniaca. Our new approach to microsatellite genotyping reveals allele-specific microsatellite and SNP markers for each locus studied. Interspecies comparison of these markers identifies alleles inherited by parthenospecies from parental species, and provides new information on origin and evolution of clonal diversity in D. armeniaca. SNP analyses reveal at least three interspecific origins of D. armeniaca, and microsatellite mutations in these initial clones give rise to new clones. Thus, we first establish multiple origins of D. armeniaca. Our study identifies the most effective molecular markers for elucidating the origins of clonal diversity in other unisexual species that arose via interspecific hybridization.