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.

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.

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.

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.

New records and geographic distribution of the sympatric zones of unisexual and bisexual rock lizards of the genus Darevskia in Armenia and adjacent territories.

BACKGROUND: Caucasian rock lizards of the genus Darevskia are unique taxa, including both bisexual and parthenogenetic species. The parthenogenetic species have originated as a result of natural hybridisation between females and males of different bisexual species. The species involved in interspecific hybridisation are called parental. However, sympatric zones (SZ) of unisexual and bisexual rock lizards of the Caucasus are still poorly studied, although they are very important for understanding the role of hybrid individuals of different origin in reticulate evolution. This paper presents the location of the SZs of parthenogenetic and their parental bisexual rock lizards of the genus Darevskia in Armenia and adjacent territories of Georgia and Nagorno-Karabakh. We summarised the locations of the SZs identified from 1957 to the present, based on our field survey data gathered in 2018-2019 and records from publications and museum collections. This dataset includes 39 SZs of three types: SZ of parental bisexual species, SZ of parental species with unisexual species and SZ of the parthenogenetic species. For each zone, species composition, geographical and altitudinal distribution are presented. New records expand our knowledge of the geographical and altitudinal distribution of SZs in these species and provide additional data for understanding the mechanisms of reticulate evolution and hybridogeneous speciation in the past, present and future. NEW INFORMATION: The new records, including geographical and altitudinal distributions of three types of SZs, are presented, which expand the previously-known list to 39 locations of contact zones for parthenogenetic and its bisexual parental species of rock lizards of the genus Darevskia in Armenia and the adjacent territories of Georgia and Nagorno-Karabakh.

Cytogenetic mechanisms of unisexuality in rock lizards.

Darevskia rock lizards is a unique complex taxa, including more than thirty species, seven of which are parthenogenetic. In mixed populations of Darevskia lizards, tri- and tetraploid forms can be found. The most important issues in the theory of reticulate evolution of Darevskia lizards are the origin of parthenogenetic species and their taxonomic position. However, there is little data on how meiosis proceeds in these species. The present work reports the complex results of cytogenetics in a diploid parthenogenetic species - D. unisexualis. Here we detail the meiotic prophase I progression and the specific features оf mitotic chromosomes organization. The stages of meiosis prophase I were investigated by immunocytochemical analysis of preparations obtained from isolated primary oocytes of D. unisexualis in comparison with maternal species D. raddei nairensis. It has been shown that in D. unisexualis at the leptotene-zygotene stages the axial elements and the synaptonemal complex (SC) form typical "bouquets". At the pachytene-diplotene stage, 18 autosomal SC-bivalents and thickened asynapted sex Z and w univalents were observed. The presence of SYCP1 protein between the lateral elements of autosomal chromosomes proved the formation of assembled SCs. Comparative genomic hybridization (CGH) on the mitotic metaphase chromosomes of D. unisexualis was carried out using the genomic DNA isolated from the parental species D. raddei nairensis and D. valentini. In the pericentromeric regions of half of the mitotic chromosomes of D. unisexualis, specific regions inherited from maternal species have been found. Following our results, we suggest a model for diploid germ cells formation from diploid oocytes without premeiotic duplication of chromosomes in the oogenesis of diploid parthenogenetic lizards D. unisexualis. Taken as a whole, our findings confirm the hybrid nature of D. unisexualis and shed light on heterozygosity and automixis in diploid parthenogenetic forms.

Extraordinary centromeres: differences in the meiotic chromosomes of two rock lizards species Darevskia portschinskii and Darevskia raddei.

According to the synthesis of 30 years of multidisciplinary studies, parthenogenetic species of rock lizards of genus Darevskia were formed as a result of different combination patterns of interspecific hybridization of the four bisexual parental species: Darevskia raddei, D. mixta, D. valentini, and D. portschinskii. In particular, D. portschinskii and D. raddei are considered as the parental species for the parthenogenetic species D. rostombekowi. Here for the first time, we present the result of comparative immunocytochemical study of primary spermatocyte nuclei spreads from the leptotene to diplotene stages of meiotic prophase I in two species: D. portschinskii and D. raddei. We observed similar chromosome lengths for both synaptonemal complex (SC) karyotypes as well as a similar number of crossing over sites. However, unexpected differences in the number and distribution of anti-centromere antibody (ACA) foci were detected in the SC structure of bivalents of the two species. In all examined D. portschinskii spermatocyte nuclei, one immunostained centromere focus was detected per SC bivalent. In contrast, in almost every studied D. raddei nuclei we identified three to nine SCs with additional immunostained ACA foci per SC bivalent. Thus, the obtained results allow us to identify species-specific karyotype features, previously not been detected using conventional mitotic chromosome analysis. Presumably the additional centromere foci are result of epigenetic chromatin modifications. We assume that this characteristic of the D. raddei karyotype could represent useful marker for the future studies of parthenogenetic species hybrid karyotypes related to D. raddei.

Reticulate Evolution of the Rock Lizards: Meiotic Chromosome Dynamics and Spermatogenesis in Diploid and Triploid Males of the Genus Darevskia.

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.

Genomic variation in parthenogenetic lizard Darevskia armeniaca: evidence from DNA fingerprinting data.

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.