Search

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., (© 2022. The Author(s).)

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.

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.

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.

One of the main questions-in the study of the unisexual (parthenogenetic) species of vertebrates is the determination of their genetic diversity. Nuclear and mitochondrial genetic markers can be used for this purpose. One of the most effective genetic markers is the microsatellite DNA, which mutates at a high rate. Development and characteristics of such markers are necessary in studies of parthenogenetic species. In this work, the analysis of the allele polymorphism of three microsatellite loci was performed for the first time via locus-specific PCR in the populations of parthenogenetic species Darevskia rostombekovi (n = 42) and bisexual parental species D. raddei (n = 6) and D. portschinskii (n = 6). All examined individuals of the parthenogenetic D. rostombekovi were heterozygous. Two to five alleles, depending on the locus, were found in the studied populations of the parthenogenetic species. It was shown that the differences were due to the varying structure of the microsatellite cluster and to single nucleotide substitutions at fixed distances in the DNA regions adjacent to the cluster. The allele structure variations form haplotype markers specific for each allele and inherited from the parental bisexual species. It was determined which alleles of the parthenogenetic species were inherited from the maternal species and which from the paternal species. Characteristics of distribution, frequency of occurrence, and combination of alleles of microsatellite loci, which determine the distinctive features of each D. rostombekovi population were obtained. The data can be used in the future to determine the clonal diversity and possible ways of its formation in the populations of the parthenogenetic species D. rostombekovi.

The article presents the genetic parameters of the populations of lizards of the Darevskia raddei complex (D. raddei nairensis and D. raddei raddei) and the populations of D. valentini calculated on the basis of the analysis of variability of 50 allelic variants of the three nuclear genome microsatellite-containing loci of 83 individuals. It was demonstrated that the F(st) genetic distances between the populations of D. raddei nairensis and D. raddei raddei were not statistically significantly different from the F(st) genetic distances between the populations of different species, D. raddei and D. valentini. At the same time, these distances were statistically significantly higher than the F(st) distances between the populations belonging to one species within the genus Darevskia. These data suggest deep divergence between the populations of D. raddei raddei and D. raddei nairensis of the D. raddei complex and there arises the question on considering them as separate species.

Cloning and sequencing of the partial reverse transcriptase gene (750 bp) of the Bov-B LINE retrotransposon have been held in parthenogenetic lizards Darevskia unisexualis and its assumed parental bisexual species D. nairensis and D.valentini. It was shown that the percentage of transcriptionally active copies of this gene, which does not contain a stop codon, is almost the same in the three species and is about 75%. The intragenomic divergence level of these sequences is low and was found to be 2.6% in D. unisexualis, 1.9% in D. nairensis, and 1.6% in D. valentini. The phylogenetic analysis shows the distribution of copies of D. unisexualis in each of the two clusters of RT sequences characteristic of D. nairensis and D. valentini. This result supports the view of the hybrid origin of D. unisexualis and does not exclude intraspecific hybridization between D. nairensis and D. valentini.

The all-female Caucasian rock lizard species Darevskia dahli and other parthenogenetic species of this genus reproduce normally via true parthenogenesis. Previously, the genetic diversity of this species was analyzed using allozymes, mitochondrial DNA, and DNA fingerprint markers. In the present study, variation at three microsatellite loci was studied in 111 specimens of D. dahli from five populations from Armenia, and new information regarding clonal diversity and clone formation in D. dahli was obtained that suggests a multiple hybridization origin. All individuals but one were heterozygous at the loci studied. Based on specific allele combinations, 11 genotypes were identified among the individuals studied. Individuals with the same genotypes formed distinct clonal lineages: one major clone was represented by 72 individuals, an intermediate clone was represented by 21 individuals, and nine other clones were rare and represented by one or several individuals. A new approach based on the detection and comparison of genotype-specific markers formed by combinations of parental-specific markers was developed and used to identify at least three hybridization founder events that resulted in the initial formation of one major and two rare clones. All other clones, including the intermediate and seven rare clones, probably arose through postformation microsatellite mutations of the major clone. This approach can be used to identify hybridization founder events and to study clone formation in other unisexual taxa.

The all-female Caucasian rock lizard species Darevskia dahli and other parthenogenetic species of this genus reproduce normally via true parthenogenesis. Previously, the genetic diversity of this species was analyzed using allozymes, mitochondrial DNA, and DNA fingerprint markers. In the present study, variation at three microsatellite loci was studied in 111 specimens of D. dahli from five populations from Armenia, and new information regarding clonal diversity and clone formation in D. dahli was obtained that suggests a multiple hybridization origin. All individuals but one were heterozygous at the loci studied. Based on specific allele combinations, 11 genotypes were identified among the individuals studied. Individuals with the same genotypes formed distinct clonal lineages: one major clone was represented by 72 individuals, an intermediate clone was represented by 21 individuals, and nine other clones were rare and represented by one or several individuals. A new approach based on the detection and comparison of genotype-specific markers formed by combinations of parental-specific markers was developed and used to identify at least three hybridization founder events that resulted in the initial formation of one major and two rare clones. All other clones, including the intermediate and seven rare clones, probably arose through postformation microsatellite mutations of the major clone. This approach can be used to identify hybridization founder events and to study clone formation in other unisexual taxa.

In the genome of unisexual (parthenogenetic) lizard Darevskia armeniaca, highly variable locus Du 161 (arm) was discovered. Analysis of allelic polymorphism was carried out using locus-specific PCR of the lizard DNA specimens from 13 isolated Armenian populations (N = 138). In the sample examined, a total of 12 Du 161(arm) alleles were identified, and their differences at the level of primary DNA structure were determined. Sequence analysis of the Du 161 (arm) alleles showed that their microsatellite clusters contained repeats of one type (GATA repeats). Allelic Du 161 (arm) variants differed in the number of GATA monomers in microsatellite, point mutations of transition and transversion types, located at fixed distances from microsatellite cluster, and by single nucleotide insertions, as well as by longer insertions located within and outside of the microsatellite cluster. Moreover, point mutations formed different combinations (haplotypes), typical of certain alleles. These combinations can be used for the analysis of the origin and inheritance of these alleles in D. armeniaca, as well as for investigation of their interspecific variation in the representatives of the genus Darevskia.

Experimental data on the molecular structure and variability of microsatellite loci in unisexual and bisexual lizard species of the genus Darevskia were analyzed. The allelic variants of Du281 and Du47 were found to differ in the number of monomers, the structure of microsatellite clusters, and point mutations in these clusters and flanking DNA. Interspecific comparison of alleles of these loci revealed both variable regions in the microsatellite clusters and allele-specific evolutionarily conserved nucleotide groups. In general, the results of comparative structural analysis of allelic variants testify to a high genetic similarity of the unisexual and bisexual lizard species studied and reveals the characteristic features of their interspecies variability.

In the present study, the first molecular genetic investigation of dinucleotide (GT)n microsatellite loci in parthenogenetic lizards Darevskia unisexualis was performed. New polymorphic locus, Du214, (GeneBank ac no. EU252542) was identified and characterized in detail. It was demonstrated that allele of this locus differed in the size and structure of microsatellite locus, as well as in point mutations, the combinations of which enabled the isolation of stabile fixed double nucleotide substitutions A-A (alleles 2 and 4) and G-T (alleles 1, 3, 5, and 6). Double nucleotide substitutions described were also identified in the orthlogous loci of the parental species genomes, D. raddei (G-T) and D. valentine (A-A). Based on the analysis of allele distribution pattern at this locus in all populations of parthenospecies D. unisexualis, mathematic model was elaborated and realized. Using this model, frequencies of allelic variants for all populations of the species of interest were calculated and population genetic structure of D. unisexualis was characterized. Genetic contribution of each population to the species gene pool was determined. The data obtained demonstrated that microsatellite variation was one of the factors of clonal and genetic diversity of a parthenospecies.

Using a pair of primers selected for the loci deltau 215, deltau 281, and deltau 323 of Darevskia unisexualis monolocus PCR analysis of orthologous loci was carried out in populations of the related parthenogenetic species D. armeniaca and in populations of bisexual parental species D. valentini and D. mixta. It was shown that the studied loci were polymorphic and represented in populations of D. armeniaca by several allelic variants. We cloned and sequenced PCR amplification products of the allelic variants of deltau 215, deltau 281, and deltau 323 loci. It was found that allelic differences of microsatellite loci were caused by variation in a number of tandem repeats in the microsatellite clusters and point mutations in the flanking regions. Interspecies comparison of the orthologous locus deltau 215 between parthenogenetic species D. armeniaca and parental species showed that two allelic variants of deltau 215 in D. armeniaca were inherited from the parental bisexual species D. mixta and D. valentini. The third allelic variant was not found in parental species and appeared because of mutation processes in genome of parthenospecies. For the first time, the information about the molecular nature of allelic polymorphism of these microsatellite loci of parthenogenetic species D. armeniaca was received in this study.

Populations of parthenogenetic lizards of the genus Darevskia consist of genetically identical animals, and represent a unique model for studying the molecular mechanisms underlying the variability and evolution of hypervariable DNA repeats. As unisexual lineages, parthenogenetic lizards are characterized by some level of genetic diversity at microsatellite loci. We cloned and sequenced a number of (GATA)n microsatellite loci of Darevskia unisexualis. PCR products from these loci were also sequenced and the degree of intraspecific polymorphism was assessed. Among the five (GATA)n loci analysed, two (Du215 and Du281) were polymorphic. Cross-species analysis of Du215 and Du281 indicate that the priming sites at the D. unisexualis loci are conserved in the bisexual parental species, D. raddei and D. valentini. Sequencing the PCR products amplified from Du215 and Du281 and from monomorphic Du323 showed that allelic differences at the polymorphic loci are caused by microsatellite mutations and by point mutations in the flanking regions. The haplotypes identified among the allelic variants of Du281 and among its orthologues in the parental species provide new evidence of the cross-species origin of D. unisexualis. To our knowledge, these data are the first to characterize the nucleotide sequences of allelic variants at microsatellite loci within parthenogenetic vertebrate animals.

Microsatellites, or short tandem repeats, are abundant across genomes of most organisms. It is evident that the most straightforward and conclusive way of studying mutations in microsatellite-containing loci is to use clonally transmitted genomes or DNA sequences inherited in multigeneration pedigrees. At present, little is known about the origin of genetic variation in species that lack effective genetic recombination. DNA fingerprinting in 43 families of the parthenogenetic lizard species Darevskia armeniaca (131 siblings), using (GACA)(4), (GGCA)(4), (GATA)(4), and (CAC)(5) probes, revealed mutant fingerprints in siblings that differed from their mothers in several restriction DNA fragments. In some cases, the mutant fingerprints detected in siblings were also found in population samples. The mutation rate for new restriction fragment length estimated by using multilocus probes varied from 0.8 x 10(-2) to 4.9 x 10(-2) per band/per sibling. Probably, the most variations detected as restriction fragment length polymorphism have germ-line origin, but somatic changes of (CAC)(n) fingerprints in adult lizards were also observed. These results provide new evidence of existing unstable regions in genomes of parthenogenetic vertebrate animals, which provide genetic variation in unisexual populations.

Using monolocus PCR analysis with the pairs of primers designed for the Du215 locus of Darevskia unisexualis, allelic polymorphism at the orthologous locus in the populations of the related parthenospecies D. armeniaca was investigated. It was demonstrated that Du215 (arm) locus was polymorphic and in the populations of parthenospecies D. armeniaca (n = 127) represented by at least three allelic variants, differing from each other by the size and composition of microsatellite cluster, and by single nucleotide substitutions in flanking DNA. Unlike the Du215 locus, Du215 (arm) was shown contain not only GATA, but also (GACA) repeats, which were absent in D. unisexualis. Thus, in this study, the data on the molecular nature of allelic polymorphism at one of the microsatellite loci of the parthenospecies D. armeniaca were reported.

The Bov-B LINE retrotransposon was first discovered in Ruminantia and was long considered to be specific for this order. Later, this mobile element was described in snakes and some lizard species. Analysis of phylogenetic relationships of Bov-B LINE elements from different ruminants, snakes, and lizard species led to the suggestion on horizontal transfer of this retrotransposon from Squamata to Ruminantia. In the Squamata group, Bov-B LINE element was found in all snakes and some lizard species examined. The element was not detected in the genomes of some species of the genera Lacerta and Podarcis. In the present study, using PCR amplification and sequencing of PCR products, Bov-B LINE element was identified in the genomes of parthenogenetic and bisexual species of the genus Darevskia (Lacertidae), as well as in such species as Lacerta agilis and Zootoca vivipara, where this retrotransposon had not been not detected before.

Using multilocus DNA fingerprinting with microsatellite probes (CAC)5, (GACA)4, (GGCA)4 and (GATA)4, intraspecific variation of the Southeast Asian lizards belonging to the genus Leiolepis (bisexual species Leiolepis reevesii and triploid parthenogenetic species Leiolepis guentherpetersi) was first examined. The L. guentherpetersi lizards were characterized by monophyletic DNA fingerprint profiles for the loci detected by the (GACA)4, (GGCA)4, and (CAC)5 probes, in terms of intrapopulation similarity index constituting S = 0.96. This was different from the individual-specific profiles of the lizards from bisexual, presumably parental species, L. reevesii (S = 0.6; P < 0.001). Genetic homogeneity of triploid L. guentherpetersi lizards at the loci examined serves as one of the arguments for the parthenogenetic nature of this species. Genetic variability of triploid parthenogenetic species L. guentherpetersi appeared to be comparable with that reported earlier for the Caucasian rock lizards of the genus Darevskia, namely, D. dahlia, D. armeniaca, and D. unisexualis (P > 0.05). The results of DNA fingerprinting analysis of the same L. guentherpetersi samples with the (GATA)4 hybridization probe were unexpected. Variability of parthenogenetic species L. guentherpetersi at the (GATA)n markers was remarkably higher than that at other DNA markers (S = 0.35; P = 3.08 x 10(-11)), being comparable to the variation of the (GATA)n DNA markers in bisexual species L. reevesii (P = 0.74). The reasons for high polymorphism of the (GATA)n-containing loci in L. guentherpetersi still remain unclear. This polymorhism is probably associated with high instability of the loci, which can be revealed by means of family analysis of parthenogenetic offspring.

During the pre-rRNA cleavage pathway, the excision of ITS2, a eukaryotic specific insertion, remains the most elusive processing step, even in yeast. Comparison of the ITS2 sequences in different organisms permits to reveal conservative, presumably functionally important elements as well as obtain new information about ITS2 divergence in evolution. We have cloned and sequenced the ITS2 of three lizard species, Agama caucasia (Agamidae), Darevskia armeniaca and Lacerta strigata (Lacertidae) and detected in them a set of specific and conservative structural elements employing secondary structure consensus for vertebrate ITS2. Furthermore, we have performed an alignment and comparative analysis of the ITS2 sequences for the two lizards families. It enables us to propose that modern lizard species formation in evolution was accompanied by ITS2 duplication in the rDNA of their common progenitors.

The genesis of mini- and microsatellite loci, which is under extensive study in humans and some other bisexual species, have been virtually overlooked in species with clonal mode of reproduction. Earlier, using multilocus DNA fingerprinting, we have examined variability of some mini- and microsatellite DNA markers in parthenogenetic lizards from the genus Darevskia. In particular, mutant (GATA)n-restrictive DNA fragments were found in Darevskia unisexualis. In the present study, we examined intraspecific polymorphism of three cloned loci of D. unisexualis--Du323, Du215, and Du281--containing (GATA)7GAT(GATA)2, GAT(GATA)9, and (GATA)10TA(GATA) microsatellite clusters, respectively. Different levels of intrapopulation and interpopulation variability of these loci were found. Locus Du281 showed the highest polymorphism--six allelic variants (in the sample of 68 DNA specimens). Three alleles were found for locus Du215. The Du325 locus was electrophoretically invariant. The primers chosen for loci Du323, Du215, and Du281 were also used for PCR analysis of homologous loci in two presumptive parental bisexual species, D. valentini and D. nairensis. The PCR products of the corresponding loci of the parental species had approximately the same size (approximately 200 bp) as their counterparts in D. unisexualis, but the polymorphism levels of the paternal, maternal, and hybrid species were shown to be somewhat different. These data on the structure of the D. unisexualis loci provide a possibility to study genetic diversity in the parthenogenetic species D. unisexualis and other related unisexual and bisexual species of this genus, which can provide new information on the origin of parthenogenetic species and on the phylogenetic relationships in the genus Darevskia. These data can also be used for resolving problems of marking the lizard genome, which is still poorly studied.

Mini- and microsatellites, comprising tandemly repeated short nucleotide sequences, are abundant dispersed repetitive elements that are ubiquitous in eukaryotic genomes. In humans and other bisexual species hypervariable mini- and microsatellite loci provide highly informative systems for monitoring of germline and somatic instability. However, little is known about the mechanisms by which these loci mutate in species that lack effective genetic recombination. Here, multilocus DNA fingerprinting was used to study M13 minisatellite and (GATA)n microsatellite instability in the parthenogenetic Caucasian rock lizard Darevskia unisexualis (Lacertidae). DNA fingerprinting of 25 parthenogenetic families, from six isolated populations in Armenia (comprising a total of 84 siblings), using the oligonucleotide (GATA)4 as a hybridization probe, revealed mutant fingerprinting phenotypes in 13 siblings that differed from their mothers in several restriction DNA fragments. In three families (8 siblings), the mutations were present in the germline. Moreover, the mutant fingerprint phenotypes detected in siblings were also present in population DNA samples. No intrafamily variations in DNA fingerprint patterns were observed with the M13 minisatellite probe. Estimates of the mutation rate for (GATA)n microsatellite loci in D. unisexualis showed that it was as high as that seen in some bisexual species, reaching 15% per sibling or 0.95% per microsatellite band. Furthermore, in one case, a somatic (GATA)n microsatellite mutation was observed in an adult lizard. These findings directly demonstrate that mutations in (GATA)n microsatellite loci comprise an important source of genetic variation in parthenogenetic populations of D. unisexualis.

Methods of estimating within- and between-population gene diversity in parthenogenetic species using mini- and microsatellite DNA markers and modified Wright's FST statistic are presented with special reference to model populations of lizards of the genus Darevskia (D. dahli, D. armeniaca, D. unisexualis). We used DNA fingerprinting data for several populations of these species examined earlier. The effects of variation in M13, minisatellite, (GACA)n and (TCC)n microsatellite loci on the formation of within-population gene diversity in parthenogenetic species D. dahli and D. armeniaca were shown to be different. The equality of the realized gene diversity H and its maximum possible value Hmax in two populations of D. dahli (Hmax = 0.032, H = 0.031, P < < 0.0431; Hmax = 0.024, H = 0.027, P = 0.09) and D. armeniaca (Hmax = 0.05, H = 0.053, P = 0.03; Hmax = 0.054, H = 0.055, P = 0.02) suggests that variation in (GACA)n loci substantially contributes to the maintenance of within-population genetic diversity. Analysis of between-population genetic diversity using loci M13, (GACA)n, and (TCC)n showed differentiation of D. dahli populations from northeastern and northwestern Armenia (FST = 0.0272, P = 3 x 10(-13)) and genetic homogeneity of the Armenian and Introduced to the Ukraine populations of D. armeniaca characteristic of one clone (FST = 0, P = 1).

Using multilocus DNA fingerprinting, we have examined variability of (TCT)n microsatellite and M13 minisatellite DNA repeats in populations, families, and tissues of Caucasian parthenogenetic rock lizards Darevskia unisexualis (Lacertidae). It has been shown for the first time that population and family DNA samples of D. unisexualis (75 samples in total) have individually specific DNA fingerprinting patterns of (TCT)n fragments. Analysis of inheritance of (TCT)n microsatellites in 46 first-generation progeny in 17 parthenogenetic D. unisexualis families revealed their extremely high instability. Mutant TCT fingerprint phenotypes were found in virtually each animal of the progeny. Moreover, varying fragments in the progeny and their original variants in the mothers were shown to simultaneously contain (TCT)n and (TCC)n polypyrimidine clusters. At the same time, no variability of (TCT)n fragments has been detected in the tissues and organs of mature parthenogenetic lizards and in the analogous tissues of the two-week-old progeny of this year. This suggests the absence of somatic mosaicism and methylation of the corresponding loci in the samples. Along with the hyperinstability of (TCT/TCC)n polypyrimidine clusters, we have shown that the population and family DNA fingerprinting patterns of M13 minisatellites were invariable and monomorphic in the same DNA samples of D. unisexualis. Our results indicate that mutations at loci containing polypyrimidine microsatellites significantly contribute to the total genomic variability of parthenogenetic lizards D. unisexualis.

Population and family samples of two morphological forms (mutant and normal with respect to dorsal color) of pathogenetic lizard Darevskia armeniaca were examined by means of DNA fingerprinting using M13 mini- and (GATA)n and (TCC)n microsatellite DNA markers. The morphological forms examined were characterized by clonally inherited, species-specific patterns of the DNA markers, which were different from the species-specific DNA fingerprints of the other parthenogenetic species of the genus Darevskia (D. dahli. D. unisexualis, and D. rostombekovi). The mean index of similarity (S) obtained for a sample of 36 individuals from three isolated populations using three types of DNA markers was 0.966. This was similar to the variability level observed in D. dahli (0.962) (P > 0.05), but higher than that in D. unisexualis (0.950) (P < 0.05) and D. rostombekovi (0.875) (P < 0.01). Inheritance of M13 minisatellite and (TCC)n microsatellite DNA markers in the F1 offspring of parthenogenetic lizards was examined. It was shown that variability and clonal diversity of the fingerprint phenotypes observed in the populations and families of D. armeniaca could be at least partly explained by RFLP mutations in microsatellite repeats.

Variation and clonal diversity in populations of the parthenogenetic rock lizard Darevskia rostombekovi was examined by means of multilocus DNA fingerprinting using mini- and microsatellite DNA markers M13, (GATA)4, and (TCC)50). The animals examined were shown to exhibit a clonally inherited, species-specific pattern of DNA markers (fingerprint profile) that is different from the species-specific patterns of parthenogenetic species D. dahli, D. armeniaca, and D. unisexualis. The mean intraspecific similarity index S was 0.950 (0.003) for a sample of 19 animals from three isolated populations of North Armenia. This significantly differed from the estimate of this parameter for a sample of 21 animals including two individuals from mountainous, relict population from the vicinity of the Sevan Lake, which was equal to 0.875 (0.001). A comparison of DNA fingerprints showed differences between 21 individuals attaining 79 DNA fragments of 1801 mini- and microsatellite markers included in the analysis. The results obtained show that intraspecific variation in D. rostombekovi is higher than that in the previously studied parthenogenetic species D. dahli (S = 0.962) and D. unisexualis (S = 0.950) (P < 0.001). Taking into account that D. rostombekovi is considered monoclonal on the basis of allozyme data, the problem of clonal variability is discussed with regard to the evidence on nuclear DNA markers. It is suggested that the hybrid karyotype of D. rostombekovi, which is more unstable than that of D. dahli and D. unisexualis, generates a series of chromosomal rearrangements (mutations). This may lead to the appearance of a geographically isolated chromosomal race (clone) in the population inhabiting the southeastern coast of the Sevan Lake.

Multilocus DNA fingerprinting has been used to study the variability of some mini- and microsatellite sequences in parthenogenetic species of Caucasian rock lizards of the genus Lacerta (L. dahli, L. armeniaca and L. unisexualis). We demonstrate that these clonally reproducing lizards possess species-specific DNA fingerprints with a low degree of intra- and interpopulation variation. Mean indices of similarity obtained using M13 DNA, (GACA)4 and (TCC)50 as probes were 0.962 and 0.966 in L. dahli and L. armeniaca, respectively. The mean index of similarity obtained using M 13 and GATA probes in L. unisexualis was estimated to be 0.95. However, despite the high degree of band-sharing, variable DNA fragments were revealed in all populations with the microsatellite probes. An particularly high level of variability was observed for (TCC)n microsatellites in populations of L. unisexualis. In fact TCC-derived DNA fingerprints were close to being individual-specific, with a mean index of similarity of 0.824. Fingerprint analysis of parthenogenetic families of L. armeniaca showed that all maternal fragments were inherited together by the progeny, and no differences in fingerprint patterns were observed. On the other hand, while identical DNA fingerprints were obtained from L. unisexualis families with M13 and (GATA)4 probes, use of the (TCC)50 probe revealed remarkable intrafamily variation in this species. It is assumed that the genetic heterogeneity observed in parthenogenetic populations may be explained, at least in part, by the existence of genetically unstable microsatellite loci. Our data serve to illustrate processes of spontaneous mutagenesis and the initial stages of clonal differentiation in natural populations of the lizard species studied.

In four isolated populations of parthenogenetic Caucasian rock lizard Lacerta unisexualis, variability of (TCC)n loci was examined using multilocus DNA fingerprinting. Unexpectedly high variability of (TCC)n microsatellites was found in all four populations. The mean similarity index was 0.825, which is higher than similarity estimates obtained for other mini- and microsatellite loci in L. unisexualis and parthenogenetic species L. dahli and L. armeniaca studied earlier. The high variation level of (TCC)n loci was shown to be at least partially associated with the presence of a diverged (TCC)n sequence fraction in the L. unisexualis genome. Mutations at some other genetically unstable (TCC)n loci may cause their structural diversity in populations of L. unisexualis.

Multilocus DNA fingerprinting was used to analyze the genome variation of mini- and microsatellite DNA regions in parthenogenetic Caucasian rock lizard Lacerta unisexualis. The DNA fingerprints obtained with probe M13 were nearly identical in all populations examined (the average similarity index S = 0.992). The fingerprints obtained with probe (GATA)4 varied (S = 0.862). Polymorphic fragments were assumed to correspond to allelic variants of genetically unstable GATA loci. Comparison of the fingerprints of animals from four geographically isolated populations revealed several population-specific GATA microsatellite markers. Based on their distribution among the populations, the corresponding alleles were assumed to originate from a common ancestral allele.

A specially optimized restriction analysis of highly repetitive DNA elements, called DNA taxonprint, was applied for phylogenetic study of primates and lizards. It was shown that electrophoretic bands of DNA repeats revealed by the taxonprint technique have valuable properties for molecular systematics. Approximately half of taxonprint bands (TB) are invariable and do not disappear from the genomes during evolution or change spontaneously. Presumably these invariable bands are restriction fragments of dispersed DNA repeats. Another group represents variable taxonprint bands that differ even between closely related species. These variable bands are probably represented by tandem DNA repeats and could be used as species-specific markers. It was shown that taxonprint bands are independent characters since the appearance of a new taxonprint band does not change the previous band pattern. Phylogenetic reconstruction carried out on taxonprint data demonstrated that this approach could be of general utility for molecular systematics and species identification.

It was earlier shown that electrophoretic patterns of the highly repetitive segments of DNA's of some animal species studied seem to have species specificity ("DNA taxonoprint method") [2]. In this paper the results of investigation of more than 20 lizard species of the Lacerta genus in comparison with representatives of some other genera and family are presented. It was shown that in taxonoprints there are one band, at least, but usually more bands common for all Lacerta species. Specimens of genus Podarcis have more similarity with genus Lacerta than genus Eremias. The representatives of Teiidae familie have no resemblance with family Lacertidae at all. Among Lacerta species there are 10, which form the separate group being very similar by most of the taxonoprints. This fact goes in line with hypothesis that 5 parthenogenetic species of this group evolve due to interspecies hybridization of another 5 bisexual species of this group. The other species studied have specific peculiarities in various taxonoprints permitting to draw some conclusions on their relations.