Variants of a major DNA satellite discriminate parental subgenomes in a hybrid parthenogenetic lizard Darevskia unisexualis (Darevsky, 1966).

Hybrid parthenogenetic animals are an exceptionally interesting model for studying the mechanisms and evolution of sexual and asexual reproduction. A diploid parthenogenetic lizard Darevskia unisexualis is a result of an ancestral cross between a maternal species Darevskia raddei nairensis and a paternal species Darevskia valentini and presents a unique opportunity for a cytogenetic and computational analysis of a hybrid karyotype. Our previous results demonstrated a significant divergence between the pericentromeric DNA sequences of the parental Darevskia species; however, an in-depth comparative study of their pericentromeres is still lacking. Here, using target sequencing of microdissected pericentromeric regions, we reveal and compare the repertoires of the pericentromeric tandem repeats of the parental Darevskia lizards. We found species-specific sequences of the major pericentromeric tandem repeat CLsat, which allowed computational prediction and experimental validation of fluorescent DNA probes discriminating parental chromosomes within the hybrid karyotype of D. unisexualis. Moreover, we have implemented a generalizable computational method, based on the optimization of the Levenshtein distance between tandem repeat monomers, for finding species-specific fluorescent probes for pericentromere staining. In total, we anticipate that our comparative analysis of Darevskia pericentromeric repeats, the species-specific fluorescent probes that we found and the pipeline that we developed will form a basis for the future detailed cytogenomic studies of a wide range of natural and laboratory hybrids.

Premeiotic endoreplication is the mechanism of obligate parthenogenesis in rock lizards of the genus Darevskia.

Among vertebrates, obligate parthenogenesis occurs exclusively in squamate reptiles. Premeiotic endoreplication in a small subset of developing oocytes has been documented as the mechanism of production of unreduced eggs in minutely explored obligate parthenogenetic lineages, namely in teiids and geckos. The situation in the lacertid genus Darevskia has been discussed for decades. Certain observations suggested that the ploidy level is restored during egg formation through a fusion of egg and polar body nuclei in Darevskia unisexualis and D. armeniaca. In this study, we re-evaluated the fusion hypothesis by studying diplotene chromosomes in adult females of sexual species D. raddei nairensis and obligate parthenogens D. armeniaca, D. dahli and D. unisexualis. We revealed 19 bivalents in the sexual species and 38 bivalents in the diploid obligate parthenogens, which uncovers premeiotic endoreplication as the mechanism of the production of non-reduced eggs in parthenogenetic females. The earlier contradicting reports can likely be attributed to the difficulty in identifying mispairing of chromosomes in pachytene, and the fact that in parthenogenetic reptiles relying on premeiotic endoreplication only a small subset of developing oocytes undergo genome doubling and overcome the pachytene checkpoint. This study highlights co-option of premeiotic endoreplication for escape from sexual reproduction in all independent hybrid origins of obligate parthenogenesis in vertebrates studied to date.

Evolution of the parthenogenetic rock lizard hybrid karyotype: Robertsonian translocation between two maternal chromosomes in Darevskia rostombekowi.

Darevskia rostombekowi, the most outstanding of the seven known parthenogenetic species in the genus Darevskia, is the result of an ancestral cross between two bisexual species Darevskia raddei and Darevskia portschinskii. The chromosomal set of this species includes a unique submetacentric autosomal chromosome; the origin of this chromosome was unresolved as only acrocentric chromosomes are described in the karyotypes of Darevskia genus normally. Here, we applied a suite of molecular cytogenetic techniques, including the mapping of telomeric (TTAGGG) n repeats using fluorescence in situ hybridization (FISH), comparative genomic hybridization (CGH), and whole-chromosome painting (WCP) in both D. rostombekowi and parental (D. portschinskii and D. raddei) species. The obtained results in total suggest that a de novo chromosomal rearrangement via Robertsonian translocation (centric fusion) between two maternal (D. raddei) acrocentric chromosomes of different size was involved in the formation of this unique submetacentric chromosome present in the parthenogenetic species D. rostombekowi. Our findings provide new data in specific and rapid evolutional processes of a unisexual reptile species karyotype.

Meiotic synapsis of homeologous chromosomes and mismatch repair protein detection in the parthenogenetic rock lizard Darevskia unisexualis.

Parthenogenetic species of Caucasian rock lizards of the genus Darevksia are important evidence for reticulate evolution and speciation by hybridization in vertebrates. Female-only lineages formed through interspecific hybridization have been discovered in many groups. Nevertheless, critical mechanisms of oogenesis and specifics of meiosis that provide long-term stability of parthenogenetic species are still unknown. Here we report cytogenetic characteristics of somatic karyotypes and meiotic prophase I nuclei in the diploid parthenogenetic species Darevskia unisexualis from the new population "Keti" in Armenia which contains an odd number of chromosomes 2n = 37, instead of the usual 2n = 38. We revealed 36 acrocentric chromosomes and a single metacentric autosomal chromosome, resulting from Robertsonian translocation. Comparative genomic hybridization revealed that chromosome fusion occurred between two chromosomes inherited from the maternal species, similar to another parthenogenetic species D. rostombekowi. To trace the chromosome behaviour in meiosis, we performed an immunocytochemical study of primary oocytes' spread nuclei and studied chromosome synapsis during meiotic prophase I in D. unisexualis based on analysis of synaptonemal complexes (SCs). We found meiotic SC-trivalent composed of one metacentric and two acrocentric chromosomes. We confirmed that the SC was assembled between homeologous chromosomes inherited from two parental species. Immunostaining of the pachytene and diplotene nuclei revealed a mismatch repair protein MLH1 loaded to all autosomal SC bivalents. Possible mechanisms of meiotic recombination between homeologous chromosomes are discussed.

Genotypic similarities among the parthenogenetic Darevskia rock lizards with different hybrid origins.

BACKGROUND: The majority of parthenogenetic vertebrates derive from hybridization between sexually reproducing species, but the exact number of hybridization events ancestral to currently extant clonal lineages is difficult to determine. Usually, we do not know whether the parental species are able to contribute their genes to the parthenogenetic vertebrate lineages after the initial hybridization. In this paper, we address the hypothesis, whether some genotypes of seven phenotypically distinct parthenogenetic rock lizards (genus Darevskia) could have resulted from back-crosses of parthenogens with their presumed parental species. We also tried to identify, as precise as possible, the ancestral populations of all seven parthenogens. RESULTS: We analysed partial mtDNA sequences and microsatellite genotypes of all seven parthenogens and their presumed ansectral species, sampled across the entire geographic range of parthenogenesis in this group. Our results confirm the previous designation of the parental species, but further specify the maternal populations that are likely ancestral to different parthenogenetic lineages. Contrary to the expectation of independent hybrid origins of the unisexual taxa, we found that genotypes at multiple loci were shared frequently between different parthenogenetic species. The highest proportions of shared genotypes were detected between (i) D. sapphirina and D. bendimahiensis and (ii) D. dahli and D. armeniaca, and less often between other parthenogens. In case (ii), genotypes at the remaining loci were notably distinct. CONCLUSIONS: We suggest that both observations (i-ii) can be explained by two parthenogenetic forms tracing their origin to a single initial hybridization event. In case (ii), however, occasional gene exchange between the unisexual and the parental bisexual species could have taken place after the onset of parthenogenetic reproduction. Indeed, backcrossed polyploid hybrids are relatively frequent in Darevskia, although no direct evidence of recent gene flow has been previously documented. Our results further suggest that parthenogens are losing heterozygosity as a result of allelic conversion, hence their fitness is expected to decline over time as genetic diversity declines. Backcrosses with the parental species could be a rescue mechanism which might prevent this decline, and therefore increase the persistance of unisexual forms.

Origin, clonal diversity, and evolution of the parthenogenetic lizard Darevskia unisexualis.

BACKGROUND: The hybridization of female D. raddei and male D. valentini gave rise to the parthenogenetic Caucasian rock lizard Darevskia unisexualis. A previously identified genetic polymorphism in the species consisted of one common and two allozyme clones. Analysis of microsatellites and single nucleotide polymorphisms (SNPs) from the three species yields estimates of clonal diversity and tests the hypothesis of a single origin for D. unisexualis. RESULTS: Genotyping and sequencing of four microsatellite-containing loci for 109 specimens of D. unisexualis, 17 D. valentini, and 45 D. raddei nairensis identified 12 presumptive clones, including one widespread and 11 rare clones. Most individuals in some localities had a rare clone. Clone-specific alleles in D. unisexualis were compared with those of the parental species. The results inferred a single hybridization event. Post-formation mutations best explain the less common clones. CONCLUSIONS: Interspecific analyses identify alleles inherited by D. unisexualis from its bisexual ancestors. SNP analyses fail to reject the hypothesis of a single interspecific origin of D. unisexualis, followed by microsatellite mutations in this initial clone. Microsatellites detect higher clonal diversity in D. unisexualis compared to allozymes and identify the likely origins of clones. Our approach may be applicable to other unisexual species whose origins involve interspecific hybridization.

The role of hybridisation in the origin and evolutionary persistence of vertebrate parthenogens: a case study of Darevskia lizards.

Obligate parthenogenesis is found in only 0.1% of the vertebrate species, is thought to be relatively short lived and is typically of hybrid origin. However, neither the evolutionary persistence of asexuality in vertebrates, nor the conditions that allow the generation of new parthenogenetic lineages are currently well understood. It has been proposed that vertebrate parthenogenetic lineages arise from hybridisation between two divergent taxa within a specific range of phylogenetic distances (the 'Balance Hypothesis'). Moreover, parthenogenetic species often maintain a certain level of hybridisation with their closest sexual relatives, potentially generating new polyploid hybrid lineages. Here we address the role of hybridisation in the origin and evolutionary lifespan of vertebrate parthenogens. We use a set of microsatellite markers to characterise the origins of parthenogens in the lizard genus Darevskia, to study the distinctiveness of sexual and asexual taxa currently in sympatry, and to analyse the evolutionary consequences of interspecific hybridisation between asexual females and sexual males. We find that parthenogens result from multiple past hybridisation events between species from specific lineages over a range of phylogenetic distances. This suggests that the Balance Hypothesis needs to allow for lineage-specific effects, as envisaged in the Phylogenetic Constraint Hypothesis. Our results show recurrent backcrossing between sexual and parthenogenic Darevskia but neither gene flow nor formation of new asexual lineages. We suggest that, along with their demographic advantage, parthenogens gain additional leverage to outcompete sexuals in nature when the retention of sexual reproductive machinery allows backcrossing with their sexual ancestors.

Correction to: The role of hybridisation in the origin and evolutionary persistence of vertebrate parthenogens: a case study of Darevskia lizards.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Multiple interspecific hybridization and microsatellite mutations provide clonal diversity in the parthenogenetic rock lizard Darevskia armeniaca.

BACKGROUND: 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. RESULTS: 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. CONCLUSION: 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.

A standard photomap of the ovarian nurse cell chromosomes for the dominant malaria vector in Europe and Middle East Anopheles sacharovi.

BACKGROUND: Anopheles sacharovi is a dominant malaria vector species in South Europe and the Middle East which has a highly plastic behaviour at both adult and larval stages. Such plasticity has prevented this species from eradication by several anti-vector campaigns. The development of new genome-based strategies for vector control will benefit from genome sequencing and physical chromosome mapping of this mosquito. Although a cytogenetic photomap for chromosomes from salivary glands of An. sacharovi has been developed, no cytogenetic map suitable for physical genome mapping is available. METHODS: Mosquitoes for this study were collected at adult stage in animal shelters in Armenia. Polytene chromosome preparations were prepared from ovarian nurse cells. Fluorescent in situ hybridization (FISH) was performed using PCR amplified probes. RESULTS: This study constructed a high-quality standard photomap for polytene chromosomes from ovarian nurse cells of An. sacharovi. Following the previous nomenclature, chromosomes were sub-divided into 39 numbered and 119 lettered sub-divisions. Chromosomal landmarks for the chromosome recognition were described. Using FISH, 4 PCR-amplified genic probes were mapped to the chromosomes. The positions of the probes demonstrated gene order reshuffling between An. sacharovi and Anopheles atroparvus which has not been seen cytologically. In addition, this study described specific chromosomal landmarks that can be used for the cytotaxonomic diagnostics of An. sacharovi based on the banding pattern of its polytene chromosomes. CONCLUSIONS: This study constructed a high-quality standard photomap for ovarian nurse cell chromosomes of An. sacharovi and validated its utility for physical genome mapping. Based on the map, cytotaxonomic features for identification of An. sacharovi have been described. The cytogenetic map constructed in this study will assist in creating a chromosome-based genome assembly for this mosquito and in developing cytotaxonomic tools for identification of other species from the Maculipennis group.

DNA damage and micronuclei in parthenogenetic and bisexual Darevskia rock lizards from the areas with different levels of soil pollution.

Natural species are widely used as indicator organisms to estimate of the impact of environmental pollution. Here we present the results of first study of a reliability of parthenogenetic Darevskia аrmeniaca and bisexual Darevskia raddei rock lizards as sentinels for monitoring of environmental genotoxicity. The comet assay and micronucleus test were applied to the lizards sampled in six areas in Armenia and Artsakh with different levels of soil contamination. The results obtained showed a clear relationship between the pollution level of lizards' habitats and the frequency of DNA damage in the comet assay. Low baseline frequency of micronuclei in D. аrmeniaca and D. raddei, however, makes this parameter ineffective for environmental genotoxicity evaluation. The parthenogenetic lizards D. аrmeniaca showed higher sensitivity toward genotoxic pollutions compared with bisexual D. raddei living in the same environment. The correlations between soil content of heavy metals Cr, Cu, Zn, Mo, Pb and DNA damage in D. аrmeniaca and between Cu, As, Mo, Pb and DNA damage in D. raddei were revealed. Overall, the lizards D. raddei and D. аrmeniaca appeared to be sensitive species in detecting soil pollution in natural environment. The application of the comet assay in Darevskia lizard species can be considered as a more appropriate method than a micronucleus test. The use of parthenogenetic lizards D. аrmeniaca as bioindicator will permit to assess the environmental genotoxicity independent of the genetic polymorphism of bisexual species.

Parthenogenesis through the ice ages: A biogeographic analysis of Caucasian rock lizards (genus Darevskia).

Darevskia rock lizards include both sexual and parthenogenetic species, mostly distributed in the heterogeneous and ecologically diverse Caucasus. The parthenogenetic species originated via directional hybridogenesis, with only some of the sexual species known to serve as parentals. However, it remains unclear when and where these events happened and how many parental lineages were involved. A multilocus phylogeographic analysis was performed on the parthenogens D. unisexualis, D. bendimahiensis and D. uzzeli, and their putative maternal species D. raddei. Results show the parthenogenetic species all have relatively recent origins, approximately 200-70kyr ago, and at least three hybridization events were involved in their formation. Ecological niche models identify the region where hybridization events leading to the formation of D. unisexualis took place, namely in the northeast of the current distribution. Models also suggest that the sexual D. raddei might have undergone a habitat shift between the Last Interglacial and the Last Glacial Maximum.

First record of Arionvulgaris Moquin-Tandon, 1855 (Arionidae) from Armenia.

Background: Arionvulgaris Moquin-Tandon, 1855 is amongst the fastest-spreading terrestrial slugs Europe-wide. In recent years, it has been recorded in Canada, Mexico and continues to expand eastwards into Eurasia. Renowned for its high invasiveness, combatting its swift spread creates significant challenges in organising effective preventative measures. New information: This study presents the first record of Arionvulgaris from Armenia, which is the second record of this species' invasion of the Caucasus. In 2022, a substantial population of A.vulgaris was observed close to the City of Stepanavan, which is also the first record in Armenia of the family Arionidae. How the species was introduced to Armenia remains unknown. Identification of Arionvulgaris was conducted, based on external and genital morphology and mitochondrial CO1 (cytochrome c oxidase subunit 1) gene sequencing, revealing notable similarities with Central European clades. Our results confirm the introduction and distribution of A.vulgaris to Armenia. Invasion of such species into Armenia will require additional monitoring and would be aided by further research on Armenia's mollusc fauna in the future.

Evolution of parthenogenetic reproduction in Caucasian rock lizards: A review.

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.

Application of land snail Helix lucorum for evaluation of genotoxicity of soil pollution.

Application of native species as sentinels allows environmental scientists to determine real genotoxic impact of environmental pollutants. The present study aims at investigating the DNA damage in the land snail Helix lucorum as a biomarker of soil pollution. For this reason, the genotoxic impact of contaminated soil on H. lucorum, collected from different polluted areas, was investigated using the comet assay in haemocytes and digestive gland cells. An increase in DNA damage was found in the snails sampled from polluted sites compared with the reference one. Strong correlations between DNA damage in haemocytes and digestive gland cells with the level of contamination indicate pollution-induced genotoxic effects in both tissues. At the same time, the digestive gland was more sensitive towards pollutants compared with haemolymph. A direct relationship between concentrations of Cu, As and Mo in soil and the number of damaged cells for hаemolymph and digestive gland tissue was found. However, the data obtained reflect the total genotoxicity of all pollutants in the studied areas. Significant correlations between the DNA damage measured by the comet assay and metal contents in soil indicate that it is a suitable biomarker in ecotoxicological studies. Our results indicate the effectiveness of H. lucorum in biomonitoring of environmental pollution.

Parthenogenesis in Darevskia lizards: A rare outcome of common hybridization, not a common outcome of rare hybridization.

Hybridization is a common evolutionary process with multiple possible outcomes. In vertebrates, interspecific hybridization has repeatedly generated parthenogenetic hybrid species. However, it is unknown whether the generation of parthenogenetic hybrids is a rare outcome of frequent hybridization between sexual species within a genus or the typical outcome of rare hybridization events. Darevskia is a genus of rock lizards with both hybrid parthenogenetic and sexual species. Using capture sequencing, we estimate phylogenetic relationships and gene flow among the sexual species, to determine how introgressive hybridization relates to the origins of parthenogenetic hybrids. We find evidence for widespread hybridization with gene flow, both between recently diverged species and deep branches. Surprisingly, we find no signal of gene flow between parental species of the parthenogenetic hybrids, suggesting that the parental pairs were either reproductively or geographically isolated early in their divergence. The generation of parthenogenetic hybrids in Darevskia is, then, a rare outcome of the total occurrence of hybridization within the genus, but the typical outcome when specific species pairs hybridize. Our results question the conventional view that parthenogenetic lineages are generated by hybridization in a window of divergence. Instead, they suggest that some lineages possess specific properties that underpin successful parthenogenetic reproduction.

First Genome of Rock Lizard Darevskia valentini Involved in Formation of Several Parthenogenetic Species.

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.

The mixed liver and kidney transcriptome dataset of Darevskia valentini rock lizard.

OBJECTIVES: This study is performed in the frame of a bigger study dedicated to genomics and transcriptomics of parthenogenesis in vertebrates. Among vertebrates, obligate parthenogenesis was first described in the lizards of the genus Darevskia. In this genus, all found parthenogenetic species originated via interspecific hybridization. It remains unknown which genetic or genomic factors play a key role in the generation of parthenogenetic organisms. Comparative genomic and transcriptomic analysis of parthenogens and their parental species may elucidate this problem. Darevskia valentini is a paternal species for four (of seven) parthenogens of this genus, which we promote as a particularly important species for the generation of parthenogenetic forms. DATA DESCRIPTION: Total cellular RNA was isolated from kidney and liver tissues using the standard Trizol Tissue RNA Extraction protocol. Sequencing of transcriptome libraries prepared by random fragmentation of cDNA samples was performed on an Illumina HiSeq2500. Obtained raw sequences contained 117,6 million reads with the GC content of 47%. After preprocessing, raw data was assembled by Trinity and produced 491,482 contigs.

Does reproductive mode affect sexually-selected coloration? Evaluating UV-blue spots in parthenogenetic and bisexual lizards of the genus Darevskia.

Sexual selection often leads to evolution of conspicuous signals, raising the chances of attracting not only potential mates, but also predators. In lacertid lizards, ultraviolet (UV)-blue spots on flanks and shoulders represent such a trait. Some level of correlation between male and female ornamentation is also known to exist. Therefore, the phenotype of females may change in the absence of sexual selection. We tested this hypothesis on a complex of parthenogenetic and bisexual lizards of the genus Darevskia. We evaluated area, counts, and chromatic properties (UV opponency, saturation) of UV-blue spots and compared the values between the clones and their bisexual progenitor species. We found a fair heterogeneity between the parthenogenetic species, but no general tendency toward higher crypsis or conspicuousness. Values of the parthenogens were not significantly different from the values of sexual females. A possible explanation is that the changes in selective forces associated with parthenogenetic reproduction are too small to affect the resulting pattern of selective pressures on the studied traits, or that the phenotypes of the parthenogens result from the unique combination of parental genomes and are conserved by clonal reproduction.

De novo transcriptome assembly and annotation of parthenogenetic lizard Darevskia unisexualis and its parental ancestors Darevskia valentini and Darevskia raddei nairensis.

Darevskia rock lizards include 29 sexual and seven parthenogenetic species of hybrid origin distributed in the Caucasus. All seven parthenogenetic species of the genus Darevskia were formed as a result of interspecific hybridization of only four sexual species. It remains unknown what are the main advantages of interspecific hybridization along with switching on parthenogenetic reproduction in evolution of reptiles. Data on whole transcriptome sequencing of parthenogens and their parental ancestors can provide value impact in solving this problem. Here we have sequenced ovary tissue transcriptomes from unisexual parthenogenetic lizard D. unisexualis and its parental bisexual ancestors to facilitate the subsequent annotation and to obtain the collinear characteristics for comparison with other lizard species. Here we report generated RNAseq data from total mRNA of ovary tissues of D. unisexualis, D. valentini and D. raddei with 58932755, 51634041 and 62788216 reads. Obtained RNA reads were assembled by Trinity assembler and 95141, 62123, 61836 contigs were identified with N50 values of 2409, 2801 and 2827 respectively. For further analysis top Gene Ontology terms were annotated for all species and transcript number was calculated. The raw data were deposited in the NCBI SRA database (BioProject PRJNA773939). The assemblies are available in Mendeley Data and can be accessed via doi:10.17632/rtd8cx7zc3.1.