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

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.

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.

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.

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.

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.

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.

SAPling: a Scan-Add-Print barcoding database system to label and track asexual organisms

We have developed a 'Scan-Add-Print' database system, SAPling, to track and monitor asexually reproducing organisms. Using barcodes to uniquely identify each animal, we can record information on the life of the individual in a computerized database containing its entire family tree. SAPling has enabled us to carry out large-scale population dynamics experiments with thousands of planarians and keep track of each individual. The database stores information such as family connections, birth date, division date and generation. We show that SAPling can be easily adapted to other asexually reproducing organisms and has a strong potential for use in large-scale and/or long-term population and senescence studies as well as studies of clonal diversity. The software is platform-independent, designed for reliability and ease of use, and provided open source from our webpage to allow project-specific customization.

INFERRING HAPLOTYPES FROM GENOTYPES ON A PEDIGREE WITH MUTATIONS, GENOTYPING ERRORS AND MISSING ALLELES

Inferring the haplotypes of the members of a pedigree from their genotypes has been extensively studied. However, most studies do not consider genotyping errors and de novo mutations. In this paper, we study how to infer haplotypes from genotype data that may contain genotyping errors, de novo mutations, and missing alleles. We assume that there are no recombinants in the genotype data, which is usually true for tightly linked markers. We introduce a combinatorial optimization problem, called haplotype configuration with mutations and errors (HCME), which calls for haplotype configurations consistent with the given genotypes that incur no recombinants and require the minimum number of mutations and errors. HCME is NP-hard. To solve the problem, we propose a heuristic algorithm, the core of which is an integer linear program (ILP) using the system of linear equations over Galois field GF(2). Our algorithm can detect and locate genotyping errors that cannot be detected by simply checking the Mendelian law of inheritance. The algorithm also offers error correction in genotypes/haplotypes rather than just detecting inconsistencies and deleting the involved loci. Our experimental results show that the algorithm can infer haplotypes with a very high accuracy and recover 65%–94% of genotyping errors depending on the pedigree topology.

A novel nucleo-cytoplasmic hybrid clone formed via androgenesis in polyploid gibel carp

BACKGROUND

Unisexual vertebrates have been demonstrated to reproduce by gynogenesis, hybridogenesis, parthenogenesis, or kleptogenesis, however, it is uncertain how the reproduction mode contributes to the clonal diversity. Recently, polyploid gibel carp has been revealed to possess coexisting dual modes of unisexual gynogenesis and sexual reproduction and to have numerous various clones. Using sexual reproduction mating between clone D female and clone A male and subsequent 7 generation multiplying of unisexual gynogenesis, we have created a novel clone strain with more than several hundred millions of individuals. Here, we attempt to identify genetic background of the novel clone and to explore the significant implication for clonal diversity contribution.

METHODS

Several nuclear genome markers and one cytoplasmic marker, the mitochondrial genome sequence, were used to identify the genetic organization of the randomly sampled individuals from different generations of the novel clone.

RESULTS

Chromosome number, Cot-1 repetitive DNA banded karyotype, microsatellite patterns, AFLP profiles and transferrin alleles uniformly indicated that nuclear genome of the novel clone is identical to that of clone A, and significantly different from that of clone D. However, the cytoplasmic marker, its complete mtDNA genome sequence, is same to that of clone D, and different from that of clone A.

CONCLUSIONS

The present data indicate that the novel clone is a nucleo-cytoplasmic hybrid between the known clones A and D, because it originates from the offspring of gonochoristic sexual reproduction mating between clone D female and clone A male, and contains an entire nuclear genome from the paternal clone A and a mtDNA genome (cytoplasm) from the maternal clone D. It is suggested to arise via androgenesis by a mechanism of ploidy doubling of clone A sperm in clone D ooplasm through inhibiting the first mitotic division. Significantly, the selected nucleo-cytoplasmic hybrid female still maintains its gynogenetic ability. Based on the present and previous findings, we discuss the association of rapid genetic changes and high genetic diversity with various ploidy levels and multiple reproduction modes in several unisexual and sexual complexes of vertebrates and even other invertebrates.