Evolution around the Red Sea: Systematics and biogeography of the agamid genus Pseudotrapelus (Squamata: Agamidae) from North Africa and Arabia

Reconstructing Squamate Biogeography in Afro-Arabia Reveals the Influence of a Complex and Dynamic Geologic Past

The geographic distribution of biodiversity is central to understanding evolutionary biology. Paleogeographic and paleoclimatic histories often help to explain how biogeographic patterns unfold through time. However, such patterns are also influenced by a variety of other factors, such as lineage diversification, that may affect the probability of certain types of biogeographic events. The complex and well-known geologic and climatic history of Afro-Arabia, together with the extensive research on reptile systematics in the region, makes Afro-Arabian squamate communities an ideal system to investigate biogeographic patterns and their drivers. Here we reconstruct the phylogenetic relationships and the ancestral geographic distributions of several Afro-Arabian reptile clades (totaling 430 species) to estimate the number of dispersal, vicariance and range contraction events. We then compare the observed biogeographic history to a distribution of simulated biogeographic events based on the empirical phylogeny and the best-fit model. This allows us to identify periods in the past where the observed biogeographic history was likely shaped by forces beyond the ones included in the model. We find an increase in vicariance following the Oligocene, most likely caused by the fragmentation of the Afro-Arabian plate. In contrast, we did not find differences between observed and expected dispersal and range contraction levels. This is consistent with diversification enhanced by environmental processes and with the establishment of a dispersal corridor connecting Africa, Arabia and Eurasia since the middle Miocene. Finally, here we show that our novel approach is useful to pinpoint events in the evolutionary history of lineages that might reflect external forces not predicted by the underlying biogeographic model.

Systematics, biogeography and evolution of the Saharo-Arabian naked-toed geckos genus Tropiocolotes

Plate tectonics constitute one of the main mechanisms of biological diversification on Earth, often being associated with cladogenetic events at different phylogenetic levels, as well as with exchange of faunas and floras across previously isolated biogeographic regions. North Africa and Arabia share a complex geological history that dates back to the break-up of the Arabian plate from the African plate ~30-25 Mya, followed by various geological events, such as the formation of the Red Sea or the connection between the African, Arabian and Eurasian plates. Species with Saharo-Arabian distributions have shown a close association between their evolutionary history and these geological events. In this study, we investigate the systematics, biogeography and evolution of the genus Tropiocolotes, a group of small ground-dwelling geckos, comprised by 12 species distributed from the Atlantic coast of North Africa to southwestern Iran. Species delimitation analyses uncovered the existence of high levels of undescribed diversity, with forms here considered at the species level including Tropiocolotes tripolitanus (Mauritania and southern Morocco), T. nattereri (southern Israel) and T. scorteccii (Yemen and Oman). Phylogenetic and biogeographic analyses recovered two main clades, an exclusively African clade and a Saharo-Arabian clade, that split ~25 Mya following the vicariant event mediated by the separation of the Arabian and African plates. The complex geological activity around the Red Sea is associated with the diversification within the Saharo-Arabian clade, including the colonization of North Africa from a second Tropiocolotes group. Results also provide new insights into the geographic distribution of Tropiocolotes nubicus, previously considered as exclusively associated to the Nile River valley, extending its known distribution further west, up to the Central Mountains of the Sahara. Accordingly, the Nile River seems to act as a major biogeographic barrier, separating Tropiocolotes nubicus and T. steudneri in their western and eastern margins, respectively.

Phylogeography of montane dragons could shed light on the history of forests and diversification processes on Sumatra

Biogeographical evidence, both, for and against the "regional endemism paradigm" hypothesis has been uncovered across the Greater Sunda Region (Sundaland) of Southeast Asia. Additionally, there are competing hypotheses regarding how Pleistocene forests may have impacted biological patterns and processes in Sumatra. Using montane agamid lizards from Sumatra, we derived and analyzed a phylogenetic dataset, genetic divergence estimates, and contemporary distributional patterns among species. We tested whether (1) Sumatra's highland Draconinae diversification fits the regional endemism paradigm hypothesis and (2) Draconinae phylogeography provides biological evidence for Pleistocene forest extent at various points in history. Our results suggest in situ diversification was the main driver behind montane Draconinae lizard diversification in Sumatra, rejecting the "regional endemism paradigm". Contemporary distribution of endemic species and their genetic relationships may potentially provide biologicalevidence for determining more precise elevational lower limits of montane forests during the Pleistocene epoch. Our data suggests montane forests did not retreat more than 700-750 m during glacial maxima because lower retreating forests would have become interconnected, allowing for widespread dispersal, exchange of gene flow and sympatric distributions contemporarily. To the contrary, our divergence estimates show that cloud forest dragons have been isolated for millions of years, suggesting there may have been a continuous disconnect between some areas, predating the Pleistocene. There may also be other ecological and evolutionary factors that impacted Draconinae distributions, such as competition, making this an excellent system for testing questions regarding montane biogeography. Additionally, we provide the first phylogeny for a wide range of Sundaland agamid species and identify some biogeographic pressures that may have triggered montane Draconinae diversification in Sumatra.

Molecular phylogeny and intraspecific differentiation of the Trapelus agilis species complex in Iran (Squamata: Agamidae) inferred from mitochondrial DNA sequences

Background

Trapelus agilis consists of different morphotypes with restricted distributions in the Iranian Plateau. The phylogeny of the species complex has not been resolved so far, but recently Trapelus sanguinolentus were elevated from this complex into a full species. Other populations of the species complex need to be evaluated taxonomically.

Methods

In the present study, several populations of this species complex along with specimens of its closely related taxa in Iran, T. sanguinolentus, T. ruderatus and T. persicus, were examined using partial nucleotide sequences of two mitochondrial genes (cytb and ND2) (total length 1,322 bp).

Result

Populations of T. sanguinolentus clustered within the T. agilis species complex, thus indicating its paraphyly, but T. sanguinolentus was previously determined to be a species based on morphological features. The T. agilis species complex forms two distinct major clades, each of which is represented by several local populations on the Iranian Plateau. At least five distinct taxa can be identified within this traditional group. Our biogeographic evaluation of the molecular dataset suggested that the Trapelus complex originated in the Late Oligocene (30 mya) and subsequently diversified during the early to middle Miocene (22–13 mya). At first, the predominantly western clade of Trapelus ruderatus diverged from the other clades (22 mya). Afterward, Trapelus persicus diverged around 18 mya ago. The broader T. agilis complex started to diverge about 16 mya, forming several clades on the Iranian Plateau and in Central Asia. The different lineages within this species complex appear to be the result of vicariance events and dispersal waives. The corresponding vicariance events are the formation of the Zagros and Kopet Dagh basins (16–14 mya), and consequently, the aridification of the Iranian Plateau in the late Miocene (11–6 Mya).

Using multilocus approach to uncover cryptic diversity within Pseudotrapelus lizards from Saudi Arabia

The agamid Pseudotrapelus lizards inhabit the mountainous areas of the Arabian Peninsula and eastern North Africa. Currently six Pseudotrapelus species are recognised, though diagnostic morphological characters are still lacking, creating great difficulty in describing new species. Recently, two specimens of Pseudotrapelus were collected from the vicinity of Riyadh in central Saudi Arabia, an area that was not sampled in previous phylogenetic studies. In here we used both mitochondrial and nuclear data to investigate the phylogenetic position of the new samples, and assess their phylogenetic relationships with the other recognised species of Pseudotrapelus from across the distribution range of the genus. We used a multilocus approach of haplotype networks, concatenated datasets and species trees, performed mitochondrial and nuclear species delimitation analyses, and estimated divergence times. In general, our results support previous molecular studies and uncover the presence of cryptic diversity within Pseudotrapelus. The phylogenetic structure of the genus is of two major clades and within them seven distinct, delimited phylogenetic groups belonging to the six recognised species and the seventh to the individuals from Riyadh. The Riyadh specimens were distinct in all analyses performed. We suggest that the new specimens from the Riyadh area are a distinct lineage, forming a clade with their phylogenetic relatives, P. sinaitus and P. chlodnickii. The clade formed by these three species diverged during the Late Miocene around 6.4 Ma, with cladogenesis possibly facilitated by vicariance and isolation caused due to climatic fluctuations and the progression of sandy areas. Our results suggest further morphological research is necessary to revise the taxonomic status of this lineage and of the entire genus.

Evolutionary history of burrowing asps (Lamprophiidae: Atractaspidinae) with emphasis on fang evolution and prey selection

Atractaspidines are poorly studied, fossorial snakes that are found throughout Africa and western Asia, including the Middle East. We employed concatenated gene-tree analyses and divergence dating approaches to investigate evolutionary relationships and biogeographic patterns of atractaspidines with a multi-locus data set consisting of three mitochondrial (16S, cyt b, and ND4) and two nuclear genes (c-mos and RAG1). We sampled 91 individuals from both atractaspidine genera (Atractaspis and Homoroselaps). Additionally, we used ancestral-state reconstructions to investigate fang and diet evolution within Atractaspidinae and its sister lineage (Aparallactinae). Our results indicated that current classification of atractaspidines underestimates diversity within the group. Diversification occurred predominantly between the Miocene and Pliocene. Ancestral-state reconstructions suggest that snake dentition in these taxa might be highly plastic within relatively short periods of time to facilitate adaptations to dynamic foraging and life-history strategies.

Ancient diversification, biogeography, and the role of climatic niche evolution in the Old World cat snakes (Colubridae, Telescopus)

The process of species diversification is often associated with niche shifts in the newly arising lineages so that interspecific competition is minimized. However, an opposing force known as niche conservatism causes that related species tend to resemble each other in their niche requirements. Due to the inherent multidimensionality of niche space, some niche components may be subject to divergent evolution while others remain conserved in the process of speciation. One such possible component is the species' climatic niche. Here, we test the role of climatic niche evolution on the diversification of the Old World cat snakes of the genus Telescopus. These slender, nocturnal snakes are distributed in arid and semiarid areas throughout Africa, southwest Asia and adjoining parts of Europe. Because phylogenetic relationships among the Telescopus species are virtually unknown, we generated sequence data for eight genetic markers from ten of the 14 described species and reconstructed a time-calibrated phylogeny of the genus. Phylogenetic analysesindicate that the genus is of considerably old origin that dates back to the Eocene/Oligocene boundary. Biogeographical analyses place the ancestor of the genus in Africa, where it diversified into the species observed today and from where it colonized Arabia and the Levant twice independently. The colonization of Arabia occurred in the Miocene, that of the Levant either in the Late Oligocene or Early Miocene. We then identified temperature and precipitation niche space and breadth of the species included in the phylogeny and examined whether there is phylogenetic signal in these climatic niche characteristics. Despite the vast range of the genus and its complex biogeographic history, most Telescopus species have similar environmental requirements with preference for arid to semiarid conditions. One may thus expect that the genus' climatic niche will be conserved. However, our results suggest that most of the climatic niche axes examined show no phylogenetic signal, being indicative of no evolutionary constraints on the climatic niche position and niche breadth in Telescopus. The only two variables with positive phylogenetic signal (temperature niche position and precipitation niche breadth) evolved under the Brownian motion model, also indicating no directional selection on these traits. As a result, climatic niche evolution does not seem to be the major driver for the diversification in Telescopus.

Effects of climate change on niche shifts of Pseudotrapelus dhofarensis and Pseudotrapelus jensvindumi (Reptilia: Agamidae) in Western Asia

Genus Pseudotrapelus has a wide distribution in North Africa and in the Middle East. In the present study, we modeled the habitat suitability of two Omani species of the genus (Pseudotrapelus dhofarensis and Pseudotrapelus jensvindumi) to evaluate the potential effects of climate change on their distribution. Mean diurnal range and precipitation of wettest quarter are the most highly contributed variables for P. jensvindumi and P. dhofarensis, respectively. The potential distribution for P. dhofarensis in the current time covers the southern coastal regions of Oman, Yemen, the Horn of Africa, and Socotra Island, but the suitable regions were reduced in the future prediction and limited to Yemen, Socotra Island, and Oman. There have not been any records of the species outside of Oman. Analysis of habitat suitability for P. jensvindumi indicated that the species is restricted to the Al Hajar Mountain of Oman and the southeast coastal region of Iran, but there are no records of the species from Iran. Because mean diurnal range will not be influenced by climate change in future, the potential distribution of the species is not expected to be changed in 2050. All predicted models were performed with the highest AUC (more than 0.97) using the Maxent method. Investigation to find unknown populations of these two species in Iran, Yemen, and Socotra Island is essential for developing conservation programs in the future.

Diversity, distribution and conservation of the terrestrial reptiles of Oman (Sauropsida, Squamata)

In the present work, we use an exceptional database including 5,359 records of 101 species of Oman’s terrestrial reptiles together with spatial tools to infer the spatial patterns of species richness and endemicity, to infer the habitat preference of each species and to better define conservation priorities, with especial focus on the effectiveness of the protected areas in preserving this unique arid fauna. Our results indicate that the sampling effort is not only remarkable from a taxonomic point of view, with multiple observations for most species, but also for the spatial coverage achieved. The observations are distributed almost continuously across the two-dimensional climatic space of Oman defined by the mean annual temperature and the total annual precipitation and across the Principal Component Analysis (PCA) of the multivariate climatic space and are well represented within 17 out of the 20 climatic clusters grouping 10% of the explained climatic variance defined by PC1 and PC2. Species richness is highest in the Hajar and Dhofar Mountains, two of the most biodiverse areas of the Arabian Peninsula, and endemic species richness is greatest in the Jebel Akhdar, the highest part of the Hajar Mountains. Oman’s 22 protected areas cover only 3.91% of the country, including within their limits 63.37% of terrestrial reptiles and 50% of all endemics. Our analyses show that large areas of the climatic space of Oman lie outside protected areas and that seven of the 20 climatic clusters are not protected at all. The results of the gap analysis indicate that most of the species are below the conservation target of 17% or even the less restrictive 12% of their total area within a protected area in order to be considered adequately protected. Therefore, an evaluation of the coverage of the current network of protected areas and the identification of priority protected areas for reptiles using reserve design algorithms are urgently needed. Our study also shows that more than half of the species are still pending of a definitive evaluation by the International Union for Conservation of Nature (IUCN).

Aridification driven diversification of fan-throated lizards from the Indian subcontinent

The establishment of monsoon climate and the consequent aridification has been one of the most important climate change episodes in the Indian subcontinent. However, little is known about how these events might have shaped the diversification patterns among the widely distributed taxa. Fan-throated lizards (FTL) (Genus: Sitana, Sarada) are widespread, diurnal and restricted to the semi-arid zones of the Indian subcontinent. We sampled FTL in 107 localities across its range. We used molecular species delimitation method and delineated 15 species including six putative species. Thirteen of them were distinguishable based on morphology but two sister species were indistinguishable and have minor overlaps in distribution. Five fossils were used to calibrate and date the phylogeny. Diversification of fan-throated lizards lineage started ~18 mya and higher lineage diversification was observed after 11 my. The initial diversification corresponds to the time when monsoon climate was established and the latter was a period of intensification of monsoon and initiation of aridification. Thirteen out of the fifteen FTL species delimited are from Peninsular India; this is probably due to the landscape heterogeneity in this region. The species poor sister genus Otocryptis is paraphyletic and probably represents relict lineages which are now confined to forested areas. Thus, the seasonality led changes in habitat, from forests to open habitats appear to have driven diversification of fan-throated lizards.

Cryptic diversity in Ptyodactylus (Reptilia: Gekkonidae) from the northern Hajar Mountains of Oman and the United Arab Emirates uncovered by an integrative taxonomic approach

The Hajar Mountains of south-eastern Arabia form an isolated massif surrounded by the sea to the east and by a large desert to the west. As a result of their old geological origin, geographical isolation, complex topography and local climate, these mountains provide an important refuge for endemic and relict species of plants and animals. With 19 species restricted to the Hajar Mountains, reptiles are the vertebrate group with the highest level of endemicity, becoming an excellent model for understanding the patterns and processes that generate and shape diversity in this arid mountain range. The geckos of the Ptyodactylus hasselquistii species complex are the largest geckos in Arabia and are found widely distributed across the Arabian Mountains, constituting a very important component of the reptile mountain fauna. Preliminary analyses suggested that their diversity in the Hajar Mountains may be higher than expected and that their systematics should be revised. In order to tackle these questions, we inferred a nearly complete calibrated phylogeny of the genus Ptyodactylus to identify the origin of the Hajar Mountains lineages using information from two mitochondrial and four nuclear genes. Genetic variability within the Hajar Mountains was further investigated using 68 specimens of Ptyodactylus from 46 localities distributed across the entire mountain range and sequenced for the same genes as above. The molecular phylogenies and morphological analyses as well as niche comparisons indicate the presence of two very old sister cryptic species living in allopatry: one restricted to the extreme northern Hajar Mountains and described as a new species herein; the other distributed across the rest of the Hajar Mountains that can be confidently assigned to the species P. orlovi. Similar to recent findings in the geckos of the genus Asaccus, the results of the present study uncover more hidden diversity in the northern Hajar Mountains and stress once again the importance of this unique mountain range as a hot spot of biodiversity and a priority focal point for reptile conservation in Arabia.

An integrative systematic revision and biogeography of Rhynchocalamus snakes (Reptilia, Colubridae) with a description of a new species from Israel

BACKGROUND

The colubrid snakes of the genus Rhynchocalamus are seldom studied and knowledge of their ecology and life history is scarce. Three species of Rhynchocalamus are currently recognized, R. satunini (from Turkey eastwards to Iran), R. arabicus (Yemen and Oman), and R. melanocephalus (from the Sinai Peninsula northwards to Turkey). All are slender, secretive, mainly nocturnal and rare fossorial snakes. This comprehensive study is the first to sample all known Rhynchocalamus species in order to review the intra-generic phylogenetic relationships and historical biogeography of the genus.

METHODS

We revised the systematics of Rhynchocalamus using an integrative approach and evaluated its phylogeography. The phylogenetic position within the Colubridae and the phylogenetic relationships within the genus were inferred using 29 individuals belonging to the three known species, with additional sampling of two other closely-related genera, Muhtarophis and Lytorhynchus. We analysed three mitochondrial (12S, 16S, cytb) and one nuclear (c-mos) gene fragments. Phylogenetic trees were reconstructed using maximum likelihood and Bayesian inference methods; the latter method also used to provide the first time-calibrated molecular phylogeny of the genus. We generated a nuclear network and carried out a topology test and species delimitation analysis. Morphological comparisons were used to differentiate among species and to describe a new species from Israel. The studied material was comprised of 108 alcohol-preserved specimens, 15 photographs, and data from the literature for the examination of 17 mensural, 14 meristic, and two categorical characters.

RESULTS

The molecular results support Rhynchocalamus as monophyletic, and as having split from its sister genus Lytorhynchus during the Late Oligocene. The three recognized species of Rhynchocalamus comprise four independently evolving groups. The molecular results reveal that the genus began to diverge during the Middle Miocene. We revealed that the best-studied species, R. melanocephalus, is paraphyletic. A population, formally ascribed to this species, from the Negev Mountain area in southern Israel is phylogenetically closer to R. arabicus from Oman than to the northern populations of the species from Israel, Syria and Turkey. Herein we describe this population as a new species: Rhynchocalamus dayanaesp. nov.

DISCUSSION

We identify four species within Rhynchocalamus: R. satunini, R. arabicus, R. melanocephalus, and R. dayanaesp. nov., the latter, to the best of our knowledge, is endemic to southern Israel. The onset of Rhynchocalamus diversification is very old and estimated to have occurred during the Middle Miocene, possibly originating in the Levant region. Radiation probably resulted from vicariance and dispersal events caused by continuous geological instability, sea-level fluctuations and climatic changes within the Levant region.

Assessing the diversity, host-specificity and infection patterns of apicomplexan parasites in reptiles from Oman, Arabia

Understanding the processes that shape parasite diversification, their distribution and abundance provides valuable information on the dynamics and evolution of disease. In this study, we assessed the diversity, distribution, host-specificity and infection patterns of apicomplexan parasites in amphibians and reptiles from Oman, Arabia. Using a quantitative PCR approach we detected three apicomplexan parasites (haemogregarines, lankesterellids and sarcocystids). A total of 13 haemogregarine haplotypes were identified, which fell into four main clades in a phylogenetic framework. Phylogenetic analysis of six new lankesterellid haplotypes revealed that these parasites were distinct from, but phylogenetically related to, knownLankesterellaspecies and might represent new taxa. The percentage of infected hosts (prevalence) and the number of haemogregarines in the blood (parasitaemia) varied significantly between gecko species. We also found significant differences in parasitaemia between haemogregarine parasite lineages (defined by phylogenetic clustering of haplotypes), suggesting differences in host–parasite compatibility between these lineages. ForPristurus rupestris, we found significant differences in haemogregarine prevalence between geographical areas. Our results suggest that host ecology and host relatedness may influence haemogregarine distributions and, more generally, highlight the importance of screening wild hosts from remote regions to provide new insights into parasite diversity.

Microendemicity in the northern Hajar Mountains of Oman and the United Arab Emirates with the description of two new species of geckos of the genus Asaccus (Squamata: Phyllodactylidae)

BACKGROUND

The Hajar Mountains of Oman and the United Arab Emirates (UAE) is the highest mountain range in Eastern Arabia. As a result of their old geological origin, geographical isolation, complex topography and local climate, these mountains provide an important refuge for endemic and relict species of plants and animals with strong Indo-Iranian affinities. Among vertebrates, the rock climbing nocturnal geckos of the genus Asaccus represent the genus with the highest number of endemic species in the Hajar Mountains. Recent taxonomic studies on the Zagros populations of Asaccus have shown that this genus is much richer than it was previously thought and preliminary morphological and molecular data suggest that its diversity in Arabia may also be underestimated.

METHODS

A total of 83 specimens originally classified as Asaccus caudivolvulus (including specimens of the two new species described herein), six other Asaccus species from the Hajar and the Zagros Mountains and two representatives of the genus Haemodracon were sequenced for up to 2,311 base pairs including the mitochondrial 12S and cytb and the nuclear c-mos, MC1R and ACM4 genes. Phylogenetic relationships were inferred using both Bayesian and maximum-likelihood approaches and the former method was also used to calibrate the phylogenetic tree. Haplotype networks and phylogenetic trees were inferred from the phased nuclear genes only. Sixty-one alcohol-preserved adult specimens originally classified as Asaccus caudivolvulus from the northern Hajar Mountains were examined for 13 morphometric and the five meristic variables using multivariate methods and were also used to diagnose and describe the two new species.

RESULTS

The results of the molecular and morphological analyses indicate that the species originally classified as Asaccus caudivolvulus is, in fact, an assemblage of three different species that started diversifying during the Mid-Miocene. The molecular phylogenies consistently recovered the Hajar endemic A. montanus as sister taxon to all the other Asaccus species included in the analyses, rendering the Arabian species of Asaccus polyphyletic.

DISCUSSION

Using this integrative approach we have uncovered a very old diversification event that has resulted in a case of microendemicity, where three morphologically and ecologically similar medium-sized lizard species coexist in a very short and narrow mountain stretch. Asaccus caudivolvulus is restricted to a small coastal area of the UAE and at risk from heavy development, while the two new species described herein are widely distributed across the northern tip of the Hajar Mountains and seem to segregate in altitude when found in close proximity in the Musandam Peninsula (Oman). Similarly to other integrative analyses of Hajar reptiles, this study highlights the high level of diversity and endemicity of this arid mountain range, underscoring its status as one of the top hotspots of reptile diversity in Arabia.