Molecular characterization of allelic variants of (GATA)n microsatellite loci in parthenogenetic lizards Darevskia unisexualis (Lacertidae)

Genetic Variability and the Potential Range of Darevskia rostombekowi in Transcaucasia

Abstract—

The results of the analysis of the genetic variability of parthenogenetic Darevskia rostombekowi (Darevsky, 1957) species using four microsatellite-containing loci are presented. Based on 118 records with geographical coordinates of the presence of this species in Transcaucasia, the maps of potential range were created. The analysis of the genetic structure of populations demonstrated that despite the established multiclonality (seven clonal lines in four populations), D. rostombekowi was formed as a result of a single act of hybridization between closely related bisexual species. The predicted distribution of D. rostombekowi using the modelling of potential range revealed new suitable habitats, where the presence of the species has not been reported previously. The results of this study and the absence of multiple acts of hybridization during the formation of these clones may indicate a regression of population size of the species. Consequently, the estimation of the conservation status of this parthenogenetic species seems to be justified.

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.

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.

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.

The origin of multiple clones in the parthenogenetic lizard species Darevskia rostombekowi

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.

Phylogenetic relationships of D. rudis (Bedriaga, 1886) and D. bithynica(Mehely, 1909) based on microsatellite and mitochondrial DNA in Turkey

The spiny-tailed lizard of the genus Darevskia have a series of taxonomic revisions, but still their phylogenetic relationships remain uncertain. In the present study, we have assessed taxonomic relationships among Darevskia bithynica and Darevskia rudis populations through estimation of phylogenetic relationships among 96 specimens using microsatellite DNA (Du215, Du281, Du323 and Du418 loci) and 53 specimens using mtDNA (16S rRNA and cytb) from main populations in Turkey. Although D. b. bithynica and D. r. mirabilis were separated based on the PCoA analysis at low level from other D. rudis and D. bithynica populations, the distance values of Nei's genetic distance, Nei's unbiased genetic distance, Fst and Linear Fst were not high among taxa in microsatellite DNA. On the other hand, our phylogenetic analyses (NJ, ML, MP and BI) did not separate D. rudis and D. bithynica populations. Finally, most of the topologically identical trees of phylogenetic analyses and microsatellite results showed that the extant populations of D. rudis and D. bithynica were found to be polytomy. Based on our molecular phylogenetic study, D. rudis complex is still ongoing revisions.

Clonal diversity and clone formation in the parthenogenetic Caucasian rock lizard Darevskia dahlia

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.

Genetic variation and de novo mutations in the parthenogenetic Caucasian rock lizard Darevskia unisexualis

Unisexual all-female lizards of the genus Darevskia that are well adapted to various habitats are known to reproduce normally by true parthenogenesis. Although they consist of unisexual lineages and lack effective genetic recombination, they are characterized by some level of genetic polymorphism. To reveal the mutational contribution to overall genetic variability, the most straightforward and conclusive way is the direct detection of mutation events in pedigree genotyping. Earlier we selected from genomic library of D. unisexualis two polymorphic microsatellite containing loci Du281 and Du215. In this study, these two loci were analyzed to detect possible de novo mutations in 168 parthenogenetic offspring of 49 D. unisexualis mothers and in 147 offspring of 50 D. armeniaca mothers. No mutant alleles were detected in D. armeniaca offspring at both loci, and in D. unisexualis offspring at the Du215 locus. There were a total of seven mutational events in the germ lines of four of the 49 D. unisexualis mothers at the Du281 locus, yielding the mutation rate of 0.1428 events per germ line tissue. Sequencing of the mutant alleles has shown that most mutations occur via deletion or insertion of single microsatellite repeat being identical in all offspring of the family. This indicates that such mutations emerge at the early stages of embryogenesis. In this study we characterized single highly unstable (GATA)(n) containing locus in parthenogenetic lizard species D. unisexualis. Besides, we characterized various types of mutant alleles of this locus found in the D. unisexualis offspring of the first generation. Our data has shown that microsatellite mutations at highly unstable loci can make a significant contribution to population variability of parthenogenetic lizards.