Single hybrid population but multiple parental individuals at the origin of parthenogenetic rock lizards Darevskia sapphirina and D. bendimahiensis Schmidtler, & Eiselt Darevsky (1994) endemic to the area of Lake Van in East Turkey

Among vertebrates, obligate parthenogenesis is only found in Squamata, where it always has a hybrid origin and a few lizard genera contain most of the known hybridogenous parthenogenetic taxa. Parthenogenesis thus seems to be pre-conditioned at the genus level, but it is not clear how often the encounter between two parental sexually reproducing species can result in the parthenogenetic offspring, nor whether the success of such hybridization event requires certain conditions or the specific time frame. To address this question, we studied the rock lizards of genus Darevskia, where a pair of parental species, D. valentini and D. raddei, as well as their parthenogenetic daughter species D. bendimahiensis and D. sapphirina, are found in close proximity NE of the Lake Van in East Anatolia. Using ddRAD-seq genotyping on 19 parental and 18 hybrid individuals, we found that (i) all parthenogenetic individuals from both D. bendimahiensis and D. sapphirina have a monophyletic origin tracing back to a single initial hybrid population, but their current genetic variation is geographically structured; (ii) unlike the most probable paternal ancestor, the genetically closest extant population of the maternal ancestor is not the geographically nearest one; and (iii) in the parthenogens, about 1% of loci carry multiple haplotypes, frequently differentiated by multiple substitutions. This pattern, in addition to biases in the relative frequency of haplotypes of maternal and paternal origin, does not appear compatible with a scenario of the entire parthenogenic clonal population having descended from a single pair of parental individuals. Instead, the data suggest that multiple parental individual ancestries still persist in the parthenogenetic gene pool. This supports the notion that although hybridization leading to parthenogenesis is generally rare at the level of species, it may be more common at the individual/population level once the right conditions are met.

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 favoured the deviation of reproductive behaviour when species switched to interspecific mates. To date, at least seven 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 two 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.

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.

Parthenogenesis is self-destructive for scaled reptiles

Parthenogenesis is rare in nature. With 39 described true parthenogens, scaled reptiles (Squamata) are the only vertebrates that evolved this reproductive strategy. Parthenogenesis is ecologically advantageous in the short term, but the young age and rarity of parthenogenetic species indicate it is less advantageous in the long term. This suggests parthenogenesis is self-destructive: it arises often but is lost due to increased extinction rates, high rates of reversal or both. However, this role of parthenogenesis as a self-destructive trait remains unknown. We used a phylogeny of Squamata (5388 species), tree metrics, null simulations and macroevolutionary scenarios of trait diversification to address the factors that best explain the rarity of parthenogenetic species. We show that parthenogenesis can be considered as self-destructive, with high extinction rates mainly responsible for its rarity in nature. Since these parthenogenetic species occur, this trait should be ecologically relevant in the short term.

Fine-scale analysis of habitat occupancy by Kura lizard (Darevskia portschinskii) and its daughter parthenogenetic form (Darevskia dahli)

Two species of rock lizards, the parthenogenetic D. dahli and the sexually reproducing D. portschinskii, coexist in a rocky outcrop in an area of ca. 1 ha, in the vicinity of Tbilisi, Georgia; the location has been well-known since the middle 1960s. The population density of the parthenogenetic lizard is five times higher than that of the sexual breeder. We studied the distribution of active lizards in space and time over three consecutive years, during the spring and autumn activity periods, to explore spatial and temporal differences between the species on a fine spatial scale. We studied the influence of temperature, humidity, and quantitative characteristics of the surface and the distance from permanent water source on the spatial distribution of D. dahli and D. portschinskii. Darevskia portschinskii was less dependent on the distance from the water source and more evenly distributed in space and time than D. dahli. Despite potential competitive interactions, the species did not avoid each other on the microhabitat scale, suggesting that the observed ecological differences are not caused by a niche shift. More individuals of the sexual breeder than individuals of the parthenogen were found in suboptimal habitats. This feature may increase the evolutionary success of D. portschinskii in a long-term perspective.

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.

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.

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.

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.

Hybridization and polyploidy in the weeping lizard Liolaemus chiliensis (Squamata: Liolaemidae)

In reptiles, polyploidy is an unusual phenomenon that can originate from interspecific hybridization. In Chile, the lizard Liolaemus chiliensis provides a unique model with which to study the origin of polyploidy because it occurs in populations of diploid individuals and a few populations that also contain triploid and diploid–triploid mosaic lizards. To test whether L. chiliensis polyploids were hybrids between lineages within the species, we analysed the karyotype, the mitochondrial gene Cytb, seven microsatellite loci, and the linear and geometric morphometry of individuals throughout the species distribution. We found polyploidy at different localities throughout the distribution. Hybrids were detected with the microsatellite loci and morphological data. A clear relationship between hybridization and polyploidy could not be established because triploid and mosaic lizards were hybrids and purebreds. However, our results support the occurrence of both phenomena, and further research is needed to clarify how hybridization could be involved in the origin of polyploidy.

Infection of parthenogenetic lizards by blood parasites does not support the "Red Queen hypothesis" but reveals the costs of sex

Sexual organisms should be better suited than asexual ones in a context of continuous evolution in response to opposite organisms in changing environments ("Red Queen" hypothesis of sex). However, sex also carries costs associated with the maintenance of males and mating (sex cost hypothesis). Here, both non-mutually excluding hypotheses are tested by analysing the infestation by haemogregarines of mixed communities of Darevskia rock lizards composed of parthenogens generated by hybridisation and their bisexual relatives. Prevalence and intensity were recorded from 339 adult lizards belonging to six species from five syntopic localities and analysed using Generalized Mixed-Models (GLMM). Both infestation parameters depended on host-size (like due to longer exposure with age), sex and, for intensity, species. Once accounting for locality and species, males were more parasitized than conspecific females with bisexual species, but no signal of reproductive mode itself on parasitization was recovered. Essentially, male-male interactions increased haemogregarine intensity while females either sexual or asexual had similar reproductive costs when in the same conditions. These findings deviate from the predictions from "Red Queen" dynamics while asymmetric gender costs are here confirmed. Thus, increased parasitization pressure on males adds to other costs, such as higher social interactions and lower fecundity, to explain why parthenogenetic lizards apparently prevail in the short-term evolutionary scale. How this is translated in the long-term requires further phylogenetic analysis.

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.

Comparing developmental stability in unisexual and bisexual rock lizards of the genus Darevskia

Parthenogenetic species are usually considered to be short-lived due to the accumulation of adverse mutations, lack of genetic variability, and inability to adapt to changing environment. If so, one may expect that the phenotype of clonal organisms may reflect such genetic and/or environmental stress. To test this hypothesis, we compared the developmental stability of bisexual and parthenogenetic lizards of the genus Darevskia. We assessed asymmetries in three meristic traits: ventral, preanal, and supratemporal scales. Our results suggest that the amount of ventral and preanal asymmetries is significantly higher in clones compared with their maternal, but not paternal, progenitor species. However, it is questionable, whether this is a consequence of clonality, as it may be considered a mild form of outbreeding depression as well. Moreover, most ventral asymmetries were found in the bisexual species Darevskia valentini. We suggest that greater differences in asymmetry levels among bisexuals may be, for instance, a consequence of the population size: the smaller the population, the higher the inbreeding and the developmental instability. On the basis of the traits examined in this study, the parthenogens do not seem to be of significantly poorer quality.

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5359-5) contains supplementary material, which is available to authorized users.

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.