Phylogeny of caucasian rock lizards (Darevskia) and other true lizards based on mitogenome analysis: Optimisation of the algorithms and gene selection

Variants of a major DNA satellite discriminate parental subgenomes in a hybrid parthenogenetic lizard Darevskia unisexualis (Darevsky, 1966)

Hybrid parthenogenetic animals are an exceptionally interesting model for studying the mechanisms and evolution of sexual and asexual reproduction. A diploid parthenogenetic lizard Darevskia unisexualis is a result of an ancestral cross between a maternal species Darevskia raddei nairensis and a paternal species Darevskia valentini and presents a unique opportunity for a cytogenetic and computational analysis of a hybrid karyotype. Our previous results demonstrated a significant divergence between the pericentromeric DNA sequences of the parental Darevskia species; however, an in-depth comparative study of their pericentromeres is still lacking. Here, using target sequencing of microdissected pericentromeric regions, we reveal and compare the repertoires of the pericentromeric tandem repeats of the parental Darevskia lizards. We found species-specific sequences of the major pericentromeric tandem repeat CLsat, which allowed computational prediction and experimental validation of fluorescent DNA probes discriminating parental chromosomes within the hybrid karyotype of D. unisexualis. Moreover, we have implemented a generalizable computational method, based on the optimization of the Levenshtein distance between tandem repeat monomers, for finding species-specific fluorescent probes for pericentromere staining. In total, we anticipate that our comparative analysis of Darevskia pericentromeric repeats, the species-specific fluorescent probes that we found and the pipeline that we developed will form a basis for the future detailed cytogenomic studies of a wide range of natural and laboratory hybrids.

Range-Wide Phylogeography and Ecological Niche Modeling Provide Insights into the Evolutionary History of the Mongolian Racerunner (Eremias argus) in Northeast Asia

The Mongolian racerunner, Eremias argus, is a small lizard endemic to Northeast Asia that can serve as an excellent model for investigating how geography and past climate change have jointly influenced the evolution of biodiversity in this region. To elucidate the processes underlying its diversification and demography, we reconstructed the range-wide phylogeographic pattern and evolutionary trajectory, using phylogenetic, population genetic, landscape genetic, Bayesian phylogeographic reconstruction and ecological niche modeling approaches. Phylogenetic analyses of the mtDNA cyt b gene revealed eight lineages that were unbounded by geographic region. The genetic structure of E. argus was mainly determined by geographic distance. Divergence dating indicated that E. argus and E. brenchleyi diverged during the Mid-Pliocene Warm Period. E. argus was estimated to have coalesced at~0.4351 Ma (Marine Isotope Stage 19). Bayesian phylogeographic diffusion analysis revealed out-of-Inner Mongolia and rapid colonization events from the end of the Last Interglacial to the Last Glacial Maximum, which is consistent with the expanded suitable range of the Last Glacial Maximum. Pre-Last Glacial Maximum growth of population is presented for most lineages of E. argus. The Glacial Maximum contraction model and the previous multiple glacial refugia hypotheses are rejected. This may be due to an increase in the amount of climatically favorable habitats in Northeast Asia. Furthermore, E. argus barbouri most likely represents an invalid taxon. The present study is the first to report a range-wide phylogeography of reptiles over such a large region in Northeast Asia. Our results make a significant contribution towards understanding the biogeography of the entire Northeast Asia.

Hidden Diversity in a Narrow Valley: Description of New Endemic Palearctic Rock Lizard Darevskia (Squamata: Lacertidae) Species from Northeastern Turkey

Darevskia is a particularly species-rich radiation of Palearctic rock lizards from the Caucasus region. Thanks to intense systematic and taxonomic research, the knowledge of species-level diversity within this genus has increased over the last quarter century. Here, we described a new species, Darevskia salihae sp. nov. from northeastern Turkey. The new taxon is differentiated from other nearby taxon by the low number of dorsal scales in the middle of the body, the shorter body length, and the absence of blue dots both on the lateral region above the forelimbs and on the margin of the ventral plates. In addition to their morphological differences, the new taxon is phylogenetically different from close groups. It is located in a separate subclade from the rudis-valentini-portschinskii subclade. This distinction is supported by both a high bootstrap value (100) and a high posterior probability value (1.00). These two subclades are separated from each other by a genetic distance of almost 4%. This separation is supported not only genetically and morphologically, but also geographically. Since the habitat of the new taxon is limited to a high mountain and a narrow valley, it does not provide an opportunity for a different Darevskia species to shelter because it creates geographical isolation. However, Darevskia parvula that live closest to the habitat of the new taxon live only at the habitat boundaries and do not enter areas where the new taxon is found. Therefore, it might be possible that while it was separated from the rudis-valentini-portschinskii group during the evolutionary transformation, it remained as a refuge and relict in a narrow area as a result of the collapse of the valleys and the partial uplift of the Kaçkar Mountains.

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.

Complete Mitochondrial Genomes of Five Racerunners (Lacertidae: Eremias) and Comparison with Other Lacertids: Insights into the Structure and Evolution of the Control Region

Comparative studies on mitochondrial genomes (mitogenomes) as well as the structure and evolution of the mitochondrial control region are few in the Lacertidae family. Here, the complete mitogenomes of five individuals of Eremias scripta (2 individuals), Eremias nikolskii, Eremias szczerbaki, and Eremias yarkandensis were determined using next-generation sequencing and were compared with other lacertids available in GenBank. The circular mitogenomes comprised the standard set of 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes and a long non-coding control region (CR). The extent of purifying selection was less pronounced for the COIII and ND2 genes in comparison with the rest of the PCGs. The codons encoding Leucine (CUN), Threonine, and Isolecucine were the three most frequently present. The secondary structure of rRNA of Lacertidae (herein, E. scripta KZL15 as an example) comprised four domains and 28 helices for 12S rRNA, with six domains and 50 helices for 16S rRNA. Five types and twenty-one subtypes of CR in Lacertidae were described by following the criteria of the presence and position of tandem repeats (TR), termination-associated sequence 1 (TAS1), termination-associated sequence 2 (TAS2), conserved sequence block 1 (CBS1), conserved sequence block 2 (CSB2), and conserved sequence block 3 (CSB3). The compositions of conserved structural elements in four genera, Acanthodactylus, Darevskia, Eremias, and Takydromus, were further explored in detail. The base composition of TAS2 – TATACATTAT in Lacertidae was updated. In addition, the motif “TAGCGGCTTTTTTG” of tandem repeats in Eremias and the motif ”GCGGCTT” in Takydromus were presented. Nucleotide lengths between CSB2 and CSB3 remained 35 bp in Eremias and Darevskia. The phylogenetic analyses of Lacertidae recovered the higher-level relationships among the three subfamilies and corroborated a hard polytomy in the Lacertinae phylogeny. The phylogenetic position of E. nikolskii challenged the monophyly of the subgenus Pareremias within Eremias. Some mismatches between the types of CR and their phylogeny demonstrated the complicated evolutionary signals of CR such as convergent evolution. These findings will promote research on the structure and evolution of the CR and highlight the need for more mitogenomes in Lacertidae.

Mitogenome analyses elucidate the evolutionary relationships of a probable Eocene wet tropics relic in the xerophile lizard genus Acanthodactylus

Climate has a large impact on diversity and evolution of the world's biota. The Eocene-Oligocene transition from tropical climate to cooler, drier environments was accompanied by global species turnover. A large number of Old World lacertid lizard lineages have diversified after the Eocene-Oligocene boundary. One of the most speciose reptile genera in the arid Palearctic, Acanthodactylus, contains two sub-Saharan species with unresolved phylogenetic relationship and unknown climatic preferences. We here aim to understand how and when adaptation to arid conditions occurred in Acanthodactylus and when tropical habitats where entered. Using whole mitogenomes from fresh and archival DNA and published sequences we recovered a well-supported Acanthodactylus phylogeny and underpinned the timing of diversification with environmental niche analyses of the sub-Saharan species A. guineensis and A. boueti in comparison to all arid Acanthodactylus. We found that A. guineensis represents an old lineage that splits from a basal node in the Western clade, and A. boueti is a derived lineage and probably not its sister. Their long branches characterize them-and especially A. guineensis-as lineages that may have persisted for a long time without further diversification or have undergone multiple extinctions. Environmental niche models verified the occurrence of A. guineensis and A. boueti in hot humid environments different from the other 42 arid Acanthodactylus species. While A. guineensis probably remained in tropical habitat from periods prior to the Eocene-Oligocene boundary, A. boueti entered tropical environments independently at a later period. Our results provide an important baseline for studying adaptation and the transition from humid to arid environments in Lacertidae.

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.

Phenotypic divergence, convergence and evolution of Caucasian rock lizards (Darevskia)

Phenotypic evolution can cause either divergent or convergent phenotypic patterns. Even adaptation to the same environment may result in divergence of some elements of phenotype, whereas for other morphological traits it could cause phenotypic convergence. We hypothesize that at least some phenotypic characters diverge monotonically, hence they evolve irreversibly even in very closely related species, and this happens in spite of multiple convergent adaptive patterns. We studied the evolution of phenotype in 13 closely related Caucasian rock lizards (Darevskia), whose phylogenetic relationships are well known. We used head shape and the outlines of three important scales, using geometric morphometrics. We studied the association of the overall head shape, individual principal components of head shape and scale outlines with four predictors: phylogeny, habitat, sex and size. The overall head shape was not correlated with any of these predictors, whereas some principal components were correlated with habitat or phylogeny. Habitat type explained the highest fraction of variation in head shape and anal scale area. The relatedness inferred from the components of phenotype not correlated with habitat was congruent with the phylogenetic tree inferred from molecular data. Although adaptation to local environments may obscure the phylogenetic signal present in phenotype, there are components of phenotype whose evolution is irreversible.