Radiation, multiple dispersal and parallelism in the skinks, Chalcides and Sphenops (Squamata: Scincidae), with comments on Scincus and Scincopus and the age of the Sahara Desert

The macroevolutionary singularity of snakes

Snakes and lizards (Squamata) represent a third of terrestrial vertebrates and exhibit spectacular innovations in locomotion, feeding, and sensory processing. However, the evolutionary drivers of this radiation remain poorly known. We infer potential causes and ultimate consequences of squamate macroevolution by combining individual-based natural history observations (>60,000 animals) with a comprehensive time-calibrated phylogeny that we anchored with genomic data (5400 loci) from 1018 species. Due to shifts in the dynamics of speciation and phenotypic evolution, snakes have transformed the trophic structure of animal communities through the recurrent origin and diversification of specialized predatory strategies. Squamate biodiversity reflects a legacy of singular events that occurred during the early history of snakes and reveals the impact of historical contingency on vertebrate biodiversity.

Hidden in the sand: Phylogenomics unravel an unexpected evolutionary history for the desert-adapted vipers of the genus Cerastes

The desert vipers of the genus Cerastes are a small clade of medically important venomous snakes within the family Viperidae. According to published morphological and molecular studies, the group is comprised by four species: two morphologically similar and phylogenetically sister taxa, the African horned viper (Cerastes cerastes) and the Arabian horned viper (Cerastes gasperettii); a more distantly related species, the Saharan sand viper (Cerastes vipera), and the enigmatic Böhme's sand viper (Cerastes boehmei), only known from a single specimen in captivity allegedly captured in Central Tunisia. In this study, we sequenced one mitochondrial marker (COI) as well as genome-wide data (ddRAD sequencing) from 28 and 41 samples, respectively, covering the entire distribution range of the genus to explore the population genomics, phylogenomic relationships and introgression patterns within the genus Cerastes. Additionally, and to provide insights into the mode of diversification of the group, we carried out niche overlap analyses considering climatic and habitat variables. Both nuclear phylogenomic reconstructions and population structure analyses have unveiled an unexpected evolutionary history for the genus Cerastes, which sharply contradicts the morphological similarities and previously published mitochondrial approaches. Cerastes cerastes and C. vipera are recovered as sister taxa whilst C. gasperettii is a sister taxon to the clade formed by these two species. We found a relatively high niche overlap (OI > 0.7) in both climatic and habitat variables between C. cerastes and C. vipera, contradicting a potential scenario of sympatric speciation. These results are in line with the introgression found between the northwestern African populations of C. cerastes and C. vipera. Finally, our genomic data confirms the existence of a lineage of C. cerastes in Arabia. All these results highlight the importance of genome-wide data over few genetic markers to study the evolutionary history of species.

DNA barcode reference library for the West Sahara-Sahel reptiles

DNA barcode reference libraries are now continuously produced for the tree of life, which are essential pillars for the study of biological diversity. Yet, our knowledge about global diversity is largely limited in undersampled regions such as the largest warm desert, the Sahara-Sahel. This dataset provides a DNA barcode reference library for the reptiles of the Western Sahara-Sahel (WSS) and neighbouring countries across this region. It includes 760 barcodes from 133 reptile taxa, distributed in 23 families, and covering the intraspecific diversity of some species. A total of 84 species were collected in the WSS (83% of the total reptile species richness) over 18 overland field expeditions conducted since 2003. DNA barcodes resulted in a high success rate (95%) of species identification and barcoding gap analysis highlighted the effectiveness of the COI fragment as a barcode marker for the WSS reptiles. This dataset represents a comprehensive and reliable DNA reference library for the WSS, filling an important biodiversity gap across a remote and hard-to-sample region.

Measurement(s)	reptile diversity
Technology Type(s)	DNA barcoding
Factor Type(s)	DNA barcode sequencing
Sample Characteristic - Organism	reptiles
Sample Characteristic - Environment	terrestrial
Sample Characteristic - Location	Western Sahara-Sahel

Aliens Coming by Ships: Distribution and Origins of the Ocellated Skink Populations in Peninsular Italy

Simple Summary

Commercial routes are reported as the main cause of biological invasions. Particularly, naval trade may accidentally bring several species to new areas where they are not native. This is particularly evident for coastal areas, where most biological invasions occur. In our work, we reported, for the first time, the presence of the ocellated skink, native to the largest Italian islands (Sardinia, Sicily and surrounding islets in a port area of continental Central Italy). We collected several individuals of this alien population and we sampled them for molecular analyses, comparing them with those naturally occurring in Sardinia, Sicily and the Mediterranean basin, including individuals accidentally introduced to peninsular Southern Italy. Differently from what previously suggested, the nucleus in Portici (Southern Italy) may have originated from Sardinia. The intense cork trade and touristic traffic between Sardinia and Southern Tuscany may have been responsible for the introduction of this lizard also to Central Italy.

Abstract

The ocellated skink (Chalcides ocellatus) is a widespread lizard, naturally distributed between the Maghreb and coastal Pakistan, with few insular populations in the Mediterranean coastal area. Some populations of this species have also been recorded in peninsular Italy, Campania and Southern Tuscany due to accidental introductions via touristic and commercial routes. In this work, we conducted genetic analyses on mitochondrial DNA COXI, cytb and 16S mtDNA genes on a sample of Italian insular and peninsular populations. Differently from what previously suggested, the nucleus in Portici (Southern Italy) may have originated from Sardinia. The intense trade and touristic traffic between Sardinia and Southern Tuscany may have been responsible for the introduction of this lizard also to Central Italy.

Reptile biodiversity in Souss-Massa National Park: an internationally important hotspot in the Mediterranean region

Souss-Massa National Park (SMNP) is Morocco’s first coastal national park, created to preserve the high diversity of its continental and marine environments. Reptiles play an essential role in balancing SMNP ecosystems, yet little work has been done to study this fauna. The present work aims at providing the first reptile inventory of SMNP since its establishment in 1991. During the period 2019 to 2020, several field surveys were carried out at 30 sites using time-constrained visual encounter surveys (TCVES), with a total sampling effort of 300 person-hours. An inventory of 23 reptile species (including four endemic species) was obtained by combining TCVES results with additional data recorded during random encounters or provided by SMNP researchers. Based on TCVES data, both sampling effort and inventory completeness were evaluated by constructing sample-based accumulation curves and calculating non-parametric estimators (Chao 1, Chao 2, Jackknife 1 and Jackknife 2). These species richness estimators suggest that the current inventory is likely to be complete. Despite its small surface area, SMNP contains nearly 20% of all known Moroccan reptile species and constitutes an important biodiversity hotspot for reptiles in the Mediterranean Region. In terms of reptile conservation concern, five species in SMNP are classified as “vulnerable”, while two species are classified as “near threatened” on the IUCN Red List, underscoring the importance of protected areas for those species.

Hidden limbs in the "limbless skink" Brachymeles lukbani: Developmental observations

Reduced limbs and limblessness have evolved independently in many lizard clades. Scincidae exhibit a wide range of limb-reduced morphologies, but only some species have been used to study the embryology of limb reduction (e.g., digit reduction in Chalcides and limb reduction in Scelotes). The genus Brachymeles, a Southeast Asian clade of skinks, includes species with a range of limb morphologies, from pentadactyl to functionally and structurally limbless species. Adults of the small, snake-like species Brachymeles lukbani show no sign of external limbs in the adult except for small depressions where they might be expected to occur. Here, we show that embryos of B. lukbani in early stages of development, on the other hand, show a truncated but well-developed limb with a stylopod and a zeugopod, but no signs of an autopod. As development proceeds, the limb's small size persists even while the embryo elongates. These observations are made based on external morphology. We used florescent whole-mount immunofluorescence to visualize the morphology of skeletal elements and muscles within the embryonic limb of B. lukabni. Early stages have a humerus and separated ulna and radius cartilages; associated with these structures are dorsal and ventral muscle masses as those found in the embryos of other limbed species. While the limb remains small, the pectoral girdle grows in proportion to the rest of the body, with well-developed skeletal elements and their associated muscles. In later stages of development, we find the small limb is still present under the skin, but there are few indications of its presence, save for the morphology of the scale covering it. By use of CT scanning, we find that the adult morphology consists of a well-developed pectoral girdle, small humerus, extremely reduced ulna and radius, and well-developed limb musculature connected to the pectoral girdle. These muscles form in association with a developing limb during embryonic stages, a hint that "limbless" lizards that possess these muscles may have or have had at least transient developing limbs, as we find in B. lukbani. Overall, this newly observed pattern of ontogenetic reduction leads to an externally limbless adult in which a limb rudiment is hidden and covered under the trunk skin, a situation called cryptomelia. The results of this work add to our growing understanding of clade-specific patterns of limb reduction and the convergent evolution of limbless phenotypes through different developmental processes.

Parasites in a hotspot: diversity and specificity patterns of apicomplexans infecting reptiles from the Socotra Archipelago

Although parasites represent a major component of biodiversity, they remain poorly assessed, especially in remote regions. In this study, we screened 461 reptiles from Socotra, the largest and most biologically diverse archipelago in Arabia. Using 18S rRNA primers, we detected various apicomplexan parasites, namely haemogregarines, sarcocystids and eimeriids. Haemogregarines were the most common and genetically diverse, followed by sarcocystids (genus Sarcocystis) and eimeriids (genera Isospora and Lankesterella). All were related to parasites of other reptiles, including species from Arabia, Northern Africa and Asia. Like their 29 endemic reptile hosts, almost all Socotran parasites presented high genetic divergence and ecological differences from those found elsewhere, and probably represent undescribed endemic species. Among hosts, skinks were the most parasitized, which contrasted with similar studies from other areas, probably due to their more generalist diet and habitat use. As expected due to its high species richness, geckos harboured the highest parasite diversity in the archipelago. Parasite diversity also seemed to be correlated to island size, as the largest island harboured most haplotypes. This study emphasizes the importance of screening parasites in wild hosts from remote regions and of considering host ecology to understand disease transmission across taxa.

Swimming through the sands of the Sahara and Arabian deserts: Phylogeny of sandfish skinks (Scincidae, Scincus) reveals a recent and rapid diversification

Large parts of the Sahara Desert and Arabia are covered by sand seas and sand dunes, which are inhabited by specialized animal communities. For example, many lizards have developed adaptations to life in loose sand, including sand-swimming behavior. The best-known sand swimmers of the Saharo-Arabia are the sandfish skinks (genus Scincus). Although there are currently only four Scincus species recognized, their phylogenetic relationships have not yet been addressed in detail. We use eight genetic markers (three mitochondrial, five nuclear) and a complete sampling of species to infer the relationships within the genus. We employ multiple phylogenetic approaches to reconstruct the evolutionary history of these skinks and to assess the level of reticulation at the onset of their radiation. Our results indicate the presence of five strongly supported species-level lineages, four represented by the currently recognized species and the fifth by S. scincus conirostris, which does not form a clade with S. scincus. Based on these results we elevate the Iranian and northern Arabian S. conirostris to the species level. The two Saharan species, S. albifasciatus and S. scincus, are sister in all analyses. Deeper relationships within the genus, however, remained largely unresolved despite the extensive genetic data set. This basal polytomy, together with the fact that we detected no sign of hybridization in the history of the genus, indicates that the diversification of the five Scincus species was rapid, burst-like, and not followed by secondary hybridization events. Divergence time estimations show a Middle Pliocene crown radiation of the genus (3.3 Mya). We hypothesize that the aridification of the Saharo-Arabia that began in the Late Miocene triggered the initial diversification of Scincus, and that the subsequent expansion of sand deserts enabled their dispersal over the large Saharan and Arabian range. We discuss the evolution of body form in sand swimming lizards and ponder how Scincus retained their fully limbed morphology despite being sand swimmers that are typically limbless.

Morphological study of the integument and corporal skeletal muscles of two psammophilous members of Scincidae (Scincus scincus and Eumeces schneideri)

Sand deserts are common biotopes on the earth's surface. Numerous morphological and physiological adaptations have appeared to cope with the peculiar conditions imposed by sandy substrates, such as abrasion, mechanical resistance and the potential low oxygen levels. The psammophilous scincids (Lepidosauria) Scincus scincus and Eumeces schneideri are among those. S. scincus is a species frequently used to study displacement inside a sandy substrate. E. schneideri is a species phylogenetically closely related to S. scincus with a similar lifestyle. The aims of this study focus on the morphology of the integument and the muscular system. Briefly, we describe interspecific differences at the superficial architecture of the scales pattern and the thickness of the integument. We highlight a high cellular turnover rate at the level of the basal germinal layer of the epidermis, which, we suggest, corresponds to an adaptation to cutaneous wear caused by abrasion. We demonstrate the presence of numerous cutaneous holocrine glands whose secretion probably plays a role in the flow of sand along the integument. Several strata of osteoderms strengthen the skin. We characterize the corporal (M. longissimus dorsi and M. rectus abdominus) and caudal muscular fibers using immunohistochemistry, and quantify them using morphometry. The musculature exhibits a high proportion of glycolytic fast fibers that allow rapid burying and are well adapted to this mechanically resistant and oxygen‐poor substrate. Oxidative slow fibers are low in abundance, less than 10% in S. scincus, but a little higher in E. schneideri.

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.

Evolution of fossorial locomotion in the transition from tetrapod to snake-like in lizards

Dramatic evolutionary transitions in morphology are often assumed to be adaptive in a new habitat. However, these assumptions are rarely tested because such tests require intermediate forms, which are often extinct. In vertebrates, the evolution of an elongate, limbless body is generally hypothesized to facilitate locomotion in fossorial and/or cluttered habitats. However, these hypotheses remain untested because few studies examine the locomotion of species ranging in body form from tetrapod to snake-like. Here, we address these functional hypotheses by testing whether trade-offs exist between locomotion in surface, fossorial and cluttered habitats in Australian Lerista lizards, which include multiple intermediate forms. We found that snake-like species penetrated sand substrates faster than more lizard-like species, representing the first direct support of the adaptation to fossoriality hypothesis. By contrast, body form did not affect surface locomotion or locomotion through cluttered leaf litter. Furthermore, all species with hindlimbs used them during both fossorial and surface locomotion. We found no evidence of a trade-off between fossorial and surface locomotion. This may be either because Lerista employed kinematic strategies that took advantage of both axial- and limb-based propulsion. This may have led to the differential occupation of their habitat, facilitating diversification of intermediate forms.

Evolutionary history of two cryptic species of northern African jerboas

BACKGROUND

Climatic variation and geologic change both play significant roles in shaping species distributions, thus affecting their evolutionary history. In Sahara-Sahel, climatic oscillations shifted the desert extent during the Pliocene-Pleistocene interval, triggering the diversification of several species. Here, we investigated how these biogeographical and ecological events have shaped patterns of genetic diversity and divergence in African Jerboas, desert specialist rodents. We focused on two sister and cryptic species, Jaculus jaculus and J. hirtipes, where we (1) evaluated their genetic differentiation, (2) reconstructed their evolutionary and demographic history; (3) tested the level of gene flow between them, and (4) assessed their ecological niche divergence.

RESULTS

The analyses based on 231 individuals sampled throughout North Africa, 8 sequence fragments (one mitochondrial and seven single copy nuclear DNA, including two candidate genes for fur coloration: MC1R and Agouti), 6 microsatellite markers and ecological modelling revealed: (1) two distinct genetic lineages with overlapping distributions, in agreement with their classification as different species, J. jaculus and J. hirtipes, with (2) low levels of gene flow and strong species divergence, (3) high haplotypic diversity without evident geographic structure within species, and (4) a low level of large-scale ecological divergence between the two taxa, suggesting species micro-habitat specialization.

CONCLUSIONS

Overall, our results suggest a speciation event that occurred during the Pliocene-Pleistocene transition. The contemporary distribution of genetic variation suggests ongoing population expansions. Despite the largely overlapping distributions at a macrogeographic scale, our genetic results suggest that the two species remain reproductively isolated, as only negligible levels of gene flow were observed. The overlapping ecological preferences at a macro-geographic scale and the ecological divergence at the micro-habitat scale suggest that local adaptation may have played a crucial role in the speciation process of these species.

The Key to Understanding the European Miocene Chalcides (Squamata, Scincidae) Comes from Asia: The Lizards of the East Siberian Tagay Locality (Baikal Lake) in Russia

The early middle Miocene (MN 5) lizards from the East Siberian Tagay locality (Baikal Lake, Russia) in Asia are described here. The lizard fauna consists of two clades, Lacertidae and Scincidae. The skink material is allocated to Chalcides. While this taxon was previously reported from Europe, it has rarely been observed in the Neogene record with only jaw fragments and frontal bones described. Its taxonomy was therefore enigmatic. The Tagay material is almost identical to the European fossils of Chalcides from Austria and Hungary, but it also contains the parietal bone. While the material is also similar to the extant Ch. ocellatus, it exhibits several morphological differences. A new species is therefore erected-Chalcides augei sp. nov. These findings further support the connection of the Baikal Lake area with central Europe during the first half of the Miocene. The comparative anatomy of the frontals, parietals and lower jaws was evaluated by micro-CT in selected skink taxa. This comparison highlights several important differences, for example, paired frontals are present in Broadleysaurus (an outgroup taxon), in Acontias and all studied members of Scincidae herein. The character optimization in Mesquite supports fused frontals as being the condition at the basal node of the Ateuchosauridae + Sphenomorphidae + Eugongylidae + Lygosomidae + Egerniidae + Mabuyidae clade. While the parapineal foramen is restricted to the parietal in most taxa studied herein, it is absent (or vestigial) in Acontias and Feylinia. In contrast to all other skinks, this foramen is located on the frontal in Ateuchosaurus chinensis. Anat Rec, 2019. © 2019 American Association for Anatomy Anat Rec, 303:1901-1934, 2020. © 2019 American Association for Anatomy.

Shortcomings of Phylogenetic Studies on Recent Radiated Insular Groups: A Meta-Analysis Using Cabo Verde Biodiversity

Over the previous decades, numerous studies focused on how oceanic islands have contributed to determine the phylogenetic relationships and times of origin and diversification of different endemic lineages. The Macaronesian Islands (i.e., Azores, Madeira, Selvagens, Canaries, and Cabo Verde), harbour biotas with exceptionally high levels of endemism. Within the region, the vascular plants and reptiles constitute two of the most important radiations. In this study we compare relevant published phylogenetic data and diversification rates retrieved within Cabo Verde endemic lineages and discuss the importance of choosing appropriate phylogeny-based methods to investigate diversification dynamics on islands. From this selective literature-based review, we summarize the software packages used in Macaronesian studies and discuss their adequacy considering the published data to obtain well-supported phylogenies in the target groups. We further debate the importance of Next Generation Sequencing (NGS), to investigate the evolutionary processes of diversification in the Macaronesian Islands. Analysis of genomic data provides phylogenetic resolution for rapidly evolving species radiations, suggesting a great potential to improve the phylogenetic signal and divergence time estimates in insular lineages. The most important Macaronesian reptile radiations provide good case-studies to compare classical phylogenetic methods with new tools, such as phylogenomics, revealing a high value for research on this hotspot area.

A paradise for parasites? Seven new haemogregarine species infecting lizards from the Canary Islands

Oceanic islands are hotspots of biodiversity due to their high levels of endemism, with the Canary Islands being a notable example. A previous molecular study on the biogeography and host associations of haemogregarines (Apicomplexa: Adeleorina) infecting lizards from this archipelago detected seven parasite haplogroups. These haplogroups exhibited high host-specificity and geographical structure, suggesting that they might correspond to distinct biological identities. In this study, along with sequencing a longer fragment of the 18S rRNA, we further explore the distinctiveness of these parasites by analysing their morphology, effects on host erythrocytes and parasitaemia levels. These lines of evidence together with their genetics, host associations, frequency of occurrence and geographical distribution support them as different biological entities. As such, we describe seven new species: Karyolysus canariensis sp. nov., Karyolysus galloti sp. nov., Karyolysus stehlini sp. nov., Karyolysus gomerensis sp. nov., Karyolysus atlanticus sp. nov., Karyolysus tinerfensis sp. nov. and Karyolysus makariogeckonis sp. nov. These new taxa are further examples of endemic diversity in the Canarian archipelago. They also contribute to clarify the taxonomy within the Apicomplexa, a phylum estimated to have one of the lowest percentages of described species.

Along for the ride or missing it altogether: exploring the host specificity and diversity of haemogregarines in the Canary Islands

Background

Host-parasite relationships are expected to be strongly shaped by host specificity, a crucial factor in parasite adaptability and diversification. Because whole host communities have to be considered to assess host specificity, oceanic islands are ideal study systems given their simplified biotic assemblages. Previous studies on insular parasites suggest host range broadening during colonization. Here, we investigate the association between one parasite group (haemogregarines) and multiple sympatric hosts (of three lizard genera: Gallotia, Chalcides and Tarentola) in the Canary Islands. Given haemogregarine characteristics and insular conditions, we hypothesized low host specificity and/or occurrence of host-switching events.

Methods

A total of 825 samples were collected from the three host taxa inhabiting the seven main islands of the Canarian Archipelago, including locations where the different lizards occurred in sympatry. Blood slides were screened to assess prevalence and parasitaemia, while parasite genetic diversity and phylogenetic relationships were inferred from 18S rRNA gene sequences.

Results

Infection levels and diversity of haplotypes varied geographically and across host groups. Infections were found in all species of Gallotia across the seven islands, in Tarentola from Tenerife, La Gomera and La Palma, and in Chalcides from Tenerife, La Gomera and El Hierro. Gallotia lizards presented the highest parasite prevalence, parasitaemia and diversity (seven haplotypes), while the other two host groups (Chalcides and Tarentola) harbored one haplotype each, with low prevalence and parasitaemia levels, and very restricted geographical ranges. Host-sharing of the same haemogregarine haplotype was only detected twice, but these rare instances likely represent occasional cross-infections.

Conclusions

Our results suggest that: (i) Canarian haemogregarine haplotypes are highly host-specific, which might have restricted parasite host expansion; (ii) haemogregarines most probably reached the Canary Islands in three colonization events with each host genus; and (iii) the high number of parasite haplotypes infecting Gallotia hosts and their restricted geographical distribution suggest co-diversification. These findings contrast with our expectations derived from results on other insular parasites, highlighting how host specificity depends on parasite characteristics and evolutionary history.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2760-5) contains supplementary material, which is available to authorized users.

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.

Historical colonization and dispersal limitation supplement climate and topography in shaping species richness of African lizards (Reptilia: Agaminae)

To what extent deep-time dispersal limitation shapes present-day biodiversity at broad spatial scales remains elusive. Here, we compiled a continental dataset on the distributions of African lizard species in the reptile subfamily Agaminae (a relatively young, Neogene radiation of agamid lizards which ancestors colonized Africa from the Arabian peninsula) and tested to what extent historical colonization and dispersal limitation (i.e. accessibility from areas of geographic origin) can explain present-day species richness relative to current climate, topography, and climate change since the late Miocene (~10 mya), the Pliocene (~3 mya), and the Last Glacial Maximum (LGM, 0.021 mya). Spatial and non-spatial multi-predictor regression models revealed that time-limited dispersal via arid corridors is a key predictor to explain macro-scale patterns of species richness. In addition, current precipitation seasonality, current temperature of the warmest month, paleo-temperature changes since the LGM and late Miocene, and topographic relief emerged as important drivers. These results suggest that deep-time dispersal constraints — in addition to climate and mountain building — strongly shape current species richness of Africa’s arid-adapted taxa. Such historical dispersal limitation might indicate that natural movement rates of species are too slow to respond to rates of ongoing and projected future climate and land use change.

A reinvestigation of phylogeny and divergence times of the Ablepharus kitaibelii species complex (Sauria, Scincidae) based on mtDNA and nuDNA genes

Morphological and DNA data support that the East Mediterranean snake-eyed skink Ablepharus kitaibelii represents a species complex that includes four species A. kitaibelii, A. budaki, A. chernovi, and A. rueppellii, highlighting the need of its taxonomic reevaluation. Here, we used Bayesian and Maximum Likelihood methods to estimate the phylogenetic relationships of all members of the complex based on two mitochondrial (cyt b, 16S rRNA) and two nuclear markers (MC1R, and NKTR) and using Chalcides, Eumeces, and Eutropis as outgroups. The biogeographic history of the complex was also investigated through the application of several phylogeographic (BEAST) and biogeographic (BBM) analyses. Paleogeographic and paleoclimatic data were used to support the inferred phylogeographic patterns. The A. kitaibelli species complex exhibits high genetic diversity, revealing cases of hidden diversity and cases of non-monophyletic species such as A. kitaibelii and A. budaki. Our results indicate that A. pannonicus branches off first and a group that comprises specimens of A. kitaibelli and A. budaki from Kastelorizo Island group (southeast Greece) and southwest Turkey, respectively is differentiated from the rest A. kitaibelli and A. budaki populations and may represent a new species. The estimated divergence times place the origin of the complex in the Middle Miocene (∼16Mya) and the divergence of most currently recognized species in the Late Miocene. The inferred ancestral distribution suggests that the complex originated in Anatolia, supposing that several vicariance and dispersal events that are related with the formation of the Mid-Aegean Trench, the Anatolian Diagonal and the orogenesis of the mountain chains in southern and eastern Anatolia have led to current distribution pattern of A. kitaibelii species complex in the Balkans and Middle East.

Development of 23 microsatellite loci for Boulenger’s agama (Agama boulengeri) with partial cross-amplification in other Agama species

Patterns of biodiversity and evolutionary processes controlling them are still poorly studied in desert biomes. Fine-scale markers could help answer some of the pressing research questions for desert biomes and Sahara in particular. Such markers are available for some large mammals and crocodiles, but not for small vertebrates. Here we present a battery of microsatellite loci developed for Agama boulengeri, a promising model to study evolutionary and demographic processes in the Sahara-Sahel. Loci were selected by sequencing enriched DNA libraries with 454 pyrosequencing. A total of 23 polymorphic loci were successfully amplified in four multiplex reactions. Cross-amplification of the microsatellite loci in A. agama and A. boueti was partially successful. These markers are a promising tool for assessing genetic diversity, gene-flow dynamics and demographic patterns in this group. Given the genus Agama is distributed throughout Africa, results presented here might also facilitate studies in other regions.

Distinct patterns of desynchronized limb regression in malagasy scincine lizards (squamata, scincidae)

Scincine lizards in Madagascar form an endemic clade of about 60 species exhibiting a variety of ecomorphological adaptations. Several subclades have adapted to burrowing and convergently regressed their limbs and eyes, resulting in a variety of partial and completely limbless morphologies among extant taxa. However, patterns of limb regression in these taxa have not been studied in detail. Here we fill this gap in knowledge by providing a phylogenetic analysis of DNA sequences of three mitochondrial and four nuclear gene fragments in an extended sampling of Malagasy skinks, and microtomographic analyses of osteology of various burrowing taxa adapted to sand substrate. Based on our data we propose to (i) consider Sirenoscincus Sakata & Hikida, 2003, as junior synonym of Voeltzkowia Boettger, 1893; (ii) resurrect the genus name Grandidierina Mocquard, 1894, for four species previously included in Voeltzkowia; and (iii) consider Androngo Brygoo, 1982, as junior synonym of Pygomeles Grandidier, 1867. By supporting the clade consisting of the limbless Voeltzkowia mira and the forelimb-only taxa V. mobydick and V. yamagishii, our data indicate that full regression of limbs and eyes occurred in parallel twice in the genus Voeltzkowia (as hitherto defined) that we consider as a sand-swimming ecomorph: in the Voeltzkowia clade sensu stricto the regression first affected the hindlimbs and subsequently the forelimbs, whereas the Grandidierina clade first regressed the forelimbs and subsequently the hindlimbs following the pattern prevalent in squamates. Timetree reconstructions for the Malagasy Scincidae contain a substantial amount of uncertainty due to the absence of suitable primary fossil calibrations. However, our preliminary reconstructions suggest rapid limb regression in Malagasy scincids with an estimated maximal duration of 6 MYr for a complete regression in Paracontias, and 4 and 8 MYr respectively for complete regression of forelimbs in Grandidierina and hindlimbs in Voeltzkowia.

First assessment on the molecular phylogeny of Anatololacerta (Squamata, Lacertidae) distributed in Southern Anatolia: insights from mtDNA and nDNA markers

The genus Anatololacerta (Lacertidae) occurs mainly in Anatolia (western and southern Turkey) and on the Aegean islands Samos, Ikaria, and Rhodos. Although its taxonomy has long been debated and is currently nascent, three morphological species have been attributed to this genus: Anatololacerta anatolica, Anatololacerta oertzeni, and Anatololacerta danfordi. Here, we investigated the evolutionary history of A. oertzeni and Anatololacerta danfordi based on both mitochondrial and nuclear markers (16S rRNA and cmos). In total, 34 Anatololacerta specimens were analyzed using maximum likelihood (ML) and Bayesian inference (BI) methods. Our results supported the presence of four well-supported lineages: two belongs to A. oertzeni and two to A. danfordi. The temporal diversification of these lineages probably started with the divergence of the first A. oertzeni lineage from western Antalya at 7.9 Mya. The other two major splits may have occurred in early Pliocene (4.4 Mya: the divergence of the second A. oertzeni from A. danfordi) and in late Pliocene (2.7 Mya: the divergence of the two lineages of A. danfordi). The phylogeographical scenario suggests that the major diversification events (from late Miocene to late Pliocene) could be related with climatic oscillations (such as the late Miocene aridification and the Messinian Salinity Crisis) and tectonic movements (such as the uplift of the central Taurus mountain).