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.

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.

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.

According to the synthesis of 30 years of multidisciplinary studies, parthenogenetic species of rock lizards of genusDarevskiawere formed as a result of different combination patterns of interspecific hybridization of the four bisexual parental species:Darevskia raddei,D. mixta,D. valentini, andD. portschinskii. In particular,D. portschinskiiandD. raddeiare considered as the parental species for the parthenogenetic speciesD. 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. portschinskiiandD. 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 examinedD. portschinskiispermatocyte nuclei, one immunostained centromere focus was detected per SC bivalent. In contrast, in almost every studiedD. raddeinuclei 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 theD. raddeikaryotype could represent useful marker for the future studies of parthenogenetic species hybrid karyotypes related toD. raddei.

Although the pericentromeric regions of chromosomes that are enriched in tandemly repeated satellite DNA represent a significant part of eukaryotic genomes, they remain understudied, which is mainly due to interdisciplinary knowledge gaps. Recent studies suggest their important role in genome regulation, karyotype stability, and evolution. Thus, the idea of satellite DNA as a junk part of the genome has been refuted. The integration of data regarding molecular composition, chromosome behaviour, and the details of the in situ organization of pericentromeric regions is of great interest. The objective of this work was a cytogenetic analysis of the interactions between pericentromeric regions from non-homologous chromosomes in mouse spermatocytes using immuno-FISH. We analysed two events: the associations between centromeric regions of the X chromosome and autosomes and the associations between the centromeric regions of the autosomal bivalents that form chromocenters. We concluded that the X chromosome forms temporary synaptic associations with different autosomes in early meiotic prophase I, which can normally be found until the pachytene–diplotene, without signs of pachytene arrest. These associations are formed between the satellite-DNA-rich centromeric regions of the X chromosome and different autosomes but do not involve the satellite-DNA-poor centromeric region of the Y chromosome. We suggest the hypothetical model of X chromosome competitive replacement from such associations during synaptic correction. We showed that the centromeric region of the X chromosome in association remains free of γH2Ax-dependent chromatin inactivation, while the Y chromosome is completely inactivated. This finding highlights the predominant role of associations between satellite DNA-rich regions of different chromosomes, including the X chromosome. We suppose that X-autosomal transient associations are a manifestation of an additional synaptic disorder checkpoint. These associations are normally corrected before the late diplotene stage. We revealed that the intense spreading conditions that were applied to the spermatocyte I nuclei did not lead to the destruction of stretched chromatin fibers of elongated chromocenters enriched in satellite DNA. The tight associations that we revealed between the pericentromeric regions of different autosomal bivalents and the X chromosome may represent the basis for a mechanism for maintaining the repeats stability in the autosomes and in the X chromosome. The consequences of our findings are discussed.