Molecular phylogeny of the Eremias persica complex of the Iranian plateau (Reptilia: Lacertidae), based on mtDNA sequences

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.

The Zagros Mountains acting as a natural barrier to gene flow in the Middle East: more evidence from the evolutionary history of spiny-tailed lizards (Uromasticinae: Saara)

We investigated the phylogeographic structure of the genus Saara and studied whether the Zagros Mountain ranges influence the gene flow and dispersal at the landscape scale between the Iranian Saara species, including S. asmussi in the central Iranian plateau and S. loricata in the Mesopotamian lowlands. Phylogenetic analyses clearly show three well-supported species, including S. loricata, S. asmussi and S. hardwickii, that are distinct from Uromastyx species. The S-DIVA and BBM analyses demonstrate that species of Saara originated from an ancestor somewhere in the Iranian Plateau and then dispersed to its current geographical range. Our results indicate that the separation of S. loricata from S. asmussi coincides with the orogenic events of the Zagros Mountains during the Pliocene, and thus it may show a vicariance event. Diversification within populations of S. loricata and S. assmussi are estimated to have occurred during the Pleistocene. The haplotype network indicates one haplogroup for each of the Iranian Saara species. Population genetic analyses shows signals of demographic expansions at the beginning of the Holocene for S. loricata. Our results support the hypothesis that the Zagros Mountains act as a barrier for gene flow and Quaternary climatic oscillations affected intraspecific genetic divergences of S. loricata and S. asmussi.

Reptile species richness associated to ecological and historical variables in Iran

Spatial gradients of species richness can be shaped by the interplay between historical and ecological factors. They might interact in particularly complex ways in heterogeneous mountainous landscapes with strong climatic and geological contrasts. We mapped the distribution of 171 lizard species to investigate species richness patterns for all species (171), diurnal species (101), and nocturnal species (70) separately. We related species richness with the historical (past climate change, mountain uplifting) and ecological variables (climate, topography and vegetation). We found that assemblages in the Western Zagros Mountains, north eastern and north western parts of Central Iranian Plateau have the highest number of lizard species. Among the investigated variables, annual mean temperature explained the largest variance for all species (10%) and nocturnal species (31%). For diurnal species, temperature change velocity shows strongest explained variance in observed richness pattern (26%). Together, our results reveal that areas with annual temperature of 15–20 °C, which receive 400–600 mm precipitation and experienced moderate level of climate change since the Last Glacial Maximum (LGM) have highest number of species. Documented patterns of our study provide a baseline for understanding the potential effect of ongoing climate change on lizard diversity in Iran.

Genetic evidence for occurrence of Macrovipera razii (Squamata, Viperidae) in the central Zagros region, Iran

This study presents the first molecular evidence of Macrovipera razii from central Zagros, more than 300 km north-west of its prior records in southern Iran. Molecular analyses based on mitochondrial cytochrome b sequences identified the individuals from central Zagros as a lineage of M. razii. Specimens from the new localities are separated by a genetic distance of 1.46% from the known populations of M. razii. The results extend the known distribution range of M. razii as an endemic species of Iran.

The phylogenetic relationships and molecular systematics of scincid lizards of the genus Heremites (Sauria, Scincidae) in the Middle East based on mtDNA sequences

The taxonomic status of species included in the genus Heremites in Iran and Iraq is uncertain. Three of these species have been assigned to the genus based on morphology: Heremites auratus transcaucasica, H. vittatus, and H. septemtaeniatus. We examined the phylogenetic relationships and taxonomic status of the Iranian and Iraqi species of Heremites by performing phylogenetic analyses using mitochondrial DNA sequences (cytochrome b and 16S rRNA). Phylogenetic relationships and estimated genetic distances indicated that the Heremites populations of the area (Iran and Iraq) form five distinct clades. Three of these clades are found only in Iran, specifically in: (1) Fars and Hormozgan provinces; (2) Northeastern Khuzestan; and (3) Khorasan and Isfahan provinces. The fourth clade (H. septemtaeniatus) is found in west and Mahshahr in Iran as well as in eastern and northern parts of Iraq. The fifth clade, Heremites vittatus, is found in Iran and Iraq. We also confirm the absence of H. auratus in Iran and Iraq. It also indicated that H. vittatus is sister taxon to the other groups that our analyses estimate the divergence of this clade in the Middle Miocene (15.9 Mya). The clade containing the Fars-Hormozgan and Khuzestan populations diverged at the end of the Miocene (8.5 Mya). The Isfahan and Khorasan populations separated at the Pliocene (4.2 Mya) from the western Iranian group, the group in Mahshahr, Iran and the groups in northern and eastern Iraq.

Pre-Quaternary divergence and subsequent radiation explain longitudinal patterns of genetic and morphological variation in the striped skink, Heremites vittatus

Background

Many animal and plant species in the Middle East and northern Africa have a predominantly longitudinal distribution, extending from Iran and Turkey along the eastern Mediterranean coast into northern Africa. These species are potentially characterized by longitudinal patterns of biological diversity, but little is known about the underlying biogeographic mechanisms and evolutionary timescales. We examined these questions in the striped skink, Heremites vittatus, one such species with a roughly longitudinal distribution across the Middle East and northern Africa, by analyzing range-wide patterns of mitochondrial DNA (mtDNA) sequence and multi-trait morphological variation.

Results

The striped skink exhibits a basic longitudinal organization of mtDNA diversity, with three major mitochondrial lineages inhabiting northern Africa, the eastern Mediterranean coast, and Turkey/Iran. Remarkably, these lineages are of pre-Quaternary origin, and are characterized by p-distances of 9–10%. In addition, within each of these lineages a more recent Quaternary genetic diversification was observed, as evidenced by deep subclades and high haplotype diversity especially in the Turkish/Iranian and eastern Mediterranean lineages. Consistent with the genetic variation, our morphological analysis revealed that the majority of morphological traits show significant mean differences between specimens from northern Africa, the eastern Mediterranean coast, and Turkey/Iran, suggesting lineage-specific trait evolution. In addition, a subset of traits exhibits clinal variation along the eastern Mediterranean coast, potentially indicating selection gradients at the geographic transition from northern Africa to Anatolia. The existence of allopatric, morphologically and genetically divergent lineages suggests that Heremites vittatus might represent a complex with several taxa.

Conclusions

Our work demonstrates that early divergence events in the Pliocene, likely driven by both climatic and geological factors, established the longitudinal patterns and distribution of Heremites vittatus. Subsequent radiation during the Pleistocene generated the genetic and morphological diversity observed today. Our study provides further evidence that longitudinal diversity patterns and species distributions in the Middle East and northern Africa were shaped by complex evolutionary processes, involving the region’s intricate geological history, climatic oscillations, and the presence of the Sahara.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0969-0) contains supplementary material, which is available to authorized users.

Geographic variation in morphology of the genus Agamura Blanford, 1874 in Iran

The genusAgamurawas previously known from four species occurring on the Iranian Plateau, but was recently revised as a monotypic genus that excluded three species. In the present study, we examined different populations ofAgamura persicamorphologically. We found thatA. persicashows geographic variation with respect to two groups, with the eastern population considered asA. cruralisandA. persicaproposed for the western population. Separation between the two populations ofA. persicawas verified based upon ANOVA results for many morphological characters, including Head Height (HH), Interorbital distance (IO), Forelimb length (FLL), Number of scales across widest part of abdomen (NSA), Loreal scales (LOS), and Number of scales that separate two adjacent tubercles (NTV). Principal Component Analysis (PCA) and Discriminant Function Analysis (DFA) based on metric and meristic characters confirmed the revision of the genusAgamura. Multivariate analysis indicated that all studied OTUs were assigned to the correct classification and have significantly different morphological characters.

Phylogenetic relationships of freshwater fishes of the genus Capoeta (Actinopterygii, Cyprinidae) in Iran

The Middle East contains a great diversity of Capoeta species, but their taxonomy remains poorly described. We used mitochondrial history to examine diversity of the algae‐scraping cyprinid Capoeta in Iran, applying the species‐delimiting approaches General Mixed Yule‐Coalescent (GMYC) and Poisson Tree Process (PTP) as well as haplotype network analyses. Using the BEAST program, we also examined temporal divergence patterns of Capoeta. The monophyly of the genus and the existence of three previously described main clades (Mesopotamian, Anatolian‐Iranian, and Aralo‐Caspian) were confirmed. However, the phylogeny proposed novel taxonomic findings within Capoeta. Results of GMYC, bPTP, and phylogenetic analyses were similar and suggested that species diversity in Iran is currently underestimated. At least four candidate species, Capoeta sp4, Capoeta sp5, Capoeta sp6, and Capoeta sp7, are awaiting description. Capoeta capoeta comprises a species complex with distinct genetic lineages. The divergence times of the three main Capoeta clades are estimated to have occurred around 15.6–12.4 Mya, consistent with a Mio‐Pleistocene origin of the diversity of Capoeta in Iran. The changes in Caspian Sea levels associated with climate fluctuations and geomorphological events such as the uplift of the Zagros and Alborz Mountains may account for the complex speciation patterns in Capoeta in Iran.

Experimental Crossing of Two Distinct Species of Leopard Geckos, Eublepharis angramainyu and E. macularius: Viability, Fertility and Phenotypic Variation of the Hybrids

Hybridization between distinct species of animals and subsequent genetic introgression plays a considerable role in the speciation process and the emergence of adaptive characters. Fitness of between-species hybrids usually sharply decreases with the divergence time of the concerned species and the divergence depth, which still allows for a successful crossing differs among principal clades of vertebrates. Recently, a review of hybridization events among distinct lizard species revealed that lizards belong to vertebrates with a highly developed ability to hybridize. In spite of this, reliable reports of experimental hybridizations between genetically fairly divergent species are only exceptional. Here, we show the results of the crossing of two distinct allopatric species of eyelid geckos possessing temperature sex determination and lacking sex chromosomes: Eublepharis macularius distributed in Pakistan/Afghanistan area and E. angramainyu, which inhabits Mesopotamia and adjacent areas. We demonstrated that F₁ hybrids were viable and fertile, and the introgression of E. angramainyu genes into the E. macularius genome can be enabled via a backcrossing. The examined hybrids (except those of the F₂ generation) displayed neither malformations nor a reduced survival. Analyses of morphometric and coloration traits confirmed phenotypic distinctness of both parental species and their F₁ hybrids. These findings contrast with long-term geographic and an evolutionary separation of the studied species. Thus, the occurrence of fertile hybrids of comparably divergent species, such as E. angramainyu and E. macularius, may also be expected in other taxa of squamates. This would violate the current estimates of species diversity in lizards.

A molecular phylogeny of the eastern group of ocellated lizard genus Timon (Sauria: Lacertidae) based on mitochondrial and nuclear DNA sequences

Timon, a small genus of lacertid lizards, includes four species distributed in two separate ranges in the western and eastern part of the Mediterranean Basin. Phylogenetic relationships between the two groups have not been resolved, and the taxonomic situation of the two subspecies of the eastern representative of the genus, Timon princeps, is not clear. To address these questions, partial DNA sequences of two nuclear (β-fibrinogen intron 7 and C-mos) and three mitochondrial (cytochrome b, 12S rRNA and 16S rRNA) genes were analyzed. Based on the high genetic distance between the two subspecies of T. princeps we promote their taxonomic status to full species, Timon princeps and Timon kurdistanicus. Divergence time estimates based on other lacertid species suggest that the separation of the green (Lacerta) and ocellated (Timon) lizards took place around 12 My ago, and that the Eastern group underwent speciation around 4-5 my ago, perhaps associated with the uplifting of the Zagros mountains. As expected given this ancient divergence and complex paleogeography, considerable levels of genetic diversity are recovered within both taxa, with geographically close individuals showing very divergent haplotypes.

Phylogeny and divergence times of some racerunner lizards (Lacertidae: Eremias) inferred from mitochondrial 16S rRNA gene segments

Eremias, or racerunners, is a widespread lacertid genus occurring in China, Mongolia, Korea, Central Asia, Southwest Asia and Southeast Europe. It has been through a series of taxonomic revisions, but the phylogenetic relationships among the species and subgenera remain unclear. In this study, a frequently studied region of the mitochondrial 16S rRNA was used to (i) reassess the phylogenetic relationships of some Eremias species, (ii) test if the viviparous species form a monophyletic group, and (iii) estimate divergence time among lineages using a Bayesian relaxed molecular-clock approach. The resulting phylogeny supports monophyly of Eremias sensu Szczerbak and a clade comprising Eremias, Acanthodactylus and Latastia. An earlier finding demonstrating monophyly of the subgenus Pareremias is corroborated, with Eremias argus being the sister taxon to Eremias brenchleyi. We present the first evidence that viviparous species form a monophyletic group. In addition, Eremias przewalskii is nested within Eremias multiocellata, suggesting that the latter is likely a paraphyletic species or a species complex. Eremias acutirostris and Eremias persica form a clade that is closely related to the subgenus Pareremias. However, the subgenera Aspidorhinus, Scapteira, and Rhabderemias seem not to be monophyletic, respectively. The Bayesian divergence-time estimation suggests that Eremias originated at about 9.9 million years ago (with the 95% confidence interval ranging from 7.6 to 12 Ma), and diversified from Late Miocene to Pleistocene. Specifically, the divergence time of the subgenus Pareremias was dated to about 6.3 million years ago (with the 95% confidence interval ranging from 5.3 to 8.5 Ma), which suggests that the diversification of this subgenus might be correlated with the evolution of an East Asian monsoon climate triggered by the rapid uplift of the Tibetan Plateau approximately 8 Ma.