The all-female Caucasian rock lizard Darevskia rostombekowi and other unisexual species of this genus reproduce normally via true parthenogenesis. Typically, diploid parthenogenetic reptiles exhibit some amount of clonal diversity. However, allozyme data from D. rostombekowi have suggested that this species consists of a single clone. Herein, we test this hypothesis by evaluating variation at three variable microsatellite loci for 42 specimens of D. rostombekowi from four populations in Armenia. Analyses based on single nucleotide polymorphisms of each locus reveal five genotypes or presumptive clones in this species. All individuals are heterozygous at the loci. The major clone occurs in 24 individuals and involves three populations. Four rare clones involve one or several individuals from one or two populations. Most variation owes to parent-specific single nucleotide polymorphisms, which occur as heterozygotes. This result fails to reject the hypothesis of a single hybridization founder event that resulted in the initial formation of one major clone. The other clones appear to have originated via post-formation microsatellite mutations of the major clone.

The extant reptiles are one of the most diverse clades among terrestrial vertebrates and one of a few groups with instances of parthenogenesis. Due to the hybrid origin of parthenogenetic species, reference genomes of the parental species as well as of the parthenogenetic progeny are indispensable to explore the genetic foundations of parthenogenetic reproduction. Here, we report on the first genome assembly of rock lizard Darevskia valentini, a paternal species for several parthenogenetic lineages. The novel genome was used in the reconstruction of the comprehensive phylogeny of Squamata inferred independently from 7369 trees of single-copy orthologs and a supermatrix of 378 conserved proteins. We also investigated Hox clusters, the loci that are often regarded as playing an important role in the speciation of animal groups with drastically diverse morphology. We demonstrated that Hox clusters of D. valentini are invaded with transposons and contain the HoxC1 gene that has been considered to be lost in the amniote ancestor. This study provides confirmation for previous works and releases new genomic data that will contribute to future discoveries on the mechanisms of parthenogenesis as well as support comparative studies among reptiles.

A key issue in the study of unisexual (parthenogenetic) vertebrate species is the determination of their genetic and clonal diversity. In pursuing this aim, various markers of nuclear and mitochondrial genomes can be used. The most effective genetic markers include microsatellite DNA, characterized by high variability. The development and characterization of such markers is a necessary step in the genetic studies of parthenogenetic species. In the present study, using locus-specific PCR, for the first time, an analysis of allelic polymorphism of four microsatellite loci is performed in the populations of parthenogenetic species Darevskia armeniaca. In the studied populations, allelic variants of each locus are identified, and the nucleotide sequences of each allele are determined. It is demonstrated that allele differences are associated with the variation in the structure of microsatellite clusters and single nucleotide substitutions at fixed distances in flanking DNA regions. Structural allele variations form haplotype markers that are specific to each allele and are inherited from their parental bisexual species. It is established which of the parental alleles of each locus were inherited by the parthenogenetic species. The characteristics of the distribution and frequency of the alleles of microsatellite loci in the populations of D. armeniaca determining specific features of each population are obtained. The observed heterozygosity of the populations at the studied loci and the mutation rates in genome regions, as well as Nei’s genetic distances between the studied populations, are determined, and the phylogenetic relationships between them are established.

На основе микросателлитного локуса Du215 генотипировано 7 предполагаемых триплоидных гибридных ящериц Darevskia. Выявлены и секвенированы аллели данного локуса у этих особей. Установлено, что 5 из 7 особей являются триплоидными гибридами Darevskia unisexualis x D. valentini, содержащими по 2 аллеля партеновида D. unisexualis и по 1 аллелю D. valentini. Остальные 2 особи, гетерозиготные по данному локусу, содержали только аллели двуполого вида D. valentini.

Seven putative triploid hybrid Darevskia lizards were genotyped on the basis of the Du215 microsatellite locus. The alleles at this locus were identified a