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.

Phylogeographic and Paleoclimatic Modelling Tools Improve Our Understanding of the Biogeographic History of Hierophis viridiflavus (Colubridae)

Phylogeographic and paleoclimatic modelling studies have been combined to infer the role of Pleistocene climatic oscillations as drivers of the genetic structure and distribution of Mediterranean taxa. For the European whip snake, Hierophis viridiflavus, previous studies based on paleoclimatic modelling have depicted a low reliability in the pattern of past climatic suitability across the central Mediterranean Basin, which barely fits the species' genetic structure. In this study, we combined phylogeographic and paleoclimatic modelling tools to improve our understanding of the biogeographic history of H. viridiflavus, particularly extending the sampling and phylogeographic inferences to previously under-sampled regions. Phylogeographic analyses recovered two major clades that diverged at the beginning of the Pleistocene and had diversified in different ways by the late Pleistocene: the east clade (composed of three subclades) and the west clade (with no further structure). Paleoclimatic models highlighted the temperate character of H. viridiflavus, indicating range contractions during both the last inter-glacial and last glacial maximum periods. Range expansions from southern-located climatic refugia likely occurred in the Bølling-Allerød and Middle Holocene periods, which are supported by signals of demographic growth in the west clade and South-East-North subclade. Overall, this work improves our understanding of the historical biogeography of H. viridiflavus, providing further insights into the evolutionary processes that occurred in the Mediterranean Basin hotspot.

Genome-wide footprints in the carob tree (Ceratonia siliqua) unveil a new domestication pattern of a fruit tree in the Mediterranean

Intense research efforts over the last two decades have renewed our understanding of plant phylogeography and domestication in the Mediterranean basin. Here we aim to investigate the evolutionary history and the origin of domestication of the carob tree (Ceratonia siliqua), which has been cultivated for millennia for food and fodder. We used >1000 microsatellite genotypes to delimit seven carob evolutionary units (CEUs). We investigated genome‐wide diversity and evolutionary patterns of the CEUs with 3557 single nucleotide polymorphisms generated by restriction‐site associated DNA sequencing (RADseq). To address the complex wild vs. cultivated status of sampled trees, we classified 56 sampled populations across the Mediterranean basin as wild, seminatural or cultivated. Nuclear and cytoplasmic loci were identified from RADseq data and separated for analyses. Phylogenetic analyses of these genomic‐wide data allowed us to resolve west‐to‐east expansions from a single long‐term refugium probably located in the foothills of the High Atlas Mountains near the Atlantic coast. Our findings support multiple origins of domestication with a low impact on the genetic diversity at range‐wide level. The carob was mostly domesticated from locally selected wild genotypes and scattered long‐distance westward dispersals of domesticated varieties by humans, concomitant with major historical migrations by Romans, Greeks and Arabs. Ex situ efforts to preserve carob genetic resources should prioritize accessions from both western and eastern populations, with emphasis on the most differentiated CEUs situated in southwest Morocco, south Spain and eastern Mediterranean. Our study highlights the relevance of wild and seminatural habitats in the conservation of genetic resources for cultivated trees.

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.

A Genus-Wide Bioactivity Analysis of Daboia (Viperinae: Viperidae) Viper Venoms Reveals Widespread Variation in Haemotoxic Properties

The snake genus Daboia (Viperidae: Viperinae; Oppel, 1811) contains five species: D. deserti, D. mauritanica, and D. palaestinae, found in Afro-Arabia, and the Russell's vipers D. russelii and D. siamensis, found in Asia. Russell's vipers are responsible for a major proportion of the medically important snakebites that occur in the regions they inhabit, and their venoms are notorious for their coagulopathic effects. While widely documented, the extent of venom variation within the Russell's vipers is poorly characterised, as is the venom activity of other species within the genus. In this study we investigated variation in the haemotoxic activity of Daboia using twelve venoms from all five species, including multiple variants of D. russelii, D. siamensis, and D. palaestinae. We tested the venoms on human plasma using thromboelastography, dose-response coagulometry analyses, and calibrated automated thrombography, and on human fibrinogen by thromboelastography and fibrinogen gels. We assessed activation of blood factors X and prothrombin by the venoms using fluorometry. Variation in venom activity was evident in all experiments. The Asian species D. russelii and D. siamensis and the African species D. mauritanica possessed procoagulant venom, while D. deserti and D. palaestinae were net-anticoagulant. Of the Russell's vipers, the venom of D. siamensis from Myanmar was most toxic and D. russelli of Sri Lanka the least. Activation of both factor X and prothrombin was evident by all venoms, though at differential levels. Fibrinogenolytic activity varied extensively throughout the genus and followed no phylogenetic trends. This venom variability underpins one of the many challenges facing treatment of Daboia snakebite envenoming. Comprehensive analyses of available antivenoms in neutralising these variable venom activities are therefore of utmost importance.

Biogeographic and demographic history of the Mediterranean snakes Malpolon monspessulanus and Hemorrhois hippocrepis across the Strait of Gibraltar

BACKGROUND

The contribution of North Africa to the assembly of biodiversity within the Western Palaearctic is still poorly documented. Since the Miocene, multiple biotic exchanges occurred across the Strait of Gibraltar, underlying the high biogeographic affinity between the western European and African sides of the Mediterranean basin. We investigated the biogeographic and demographic dynamics of two large Mediterranean-adapted snakes across the Strait and assess their relevance to the origin and diversity patterns of current European and North African populations.

RESULTS

We inferred phylogeographic patterns and demographic history of M. monspessulanus and H. hippocrepis, based on range-wide multilocus data, combined with fossil data and species distribution modelling, under present and past bioclimatic envelopes. For both species we identified endemic lineages in the High Atlas Mountains (Morocco) and in eastern Iberia, suggesting their persistence in Europe during the Pleistocene. One lineage is shared between North Africa and southern Iberia and likely spread from the former to the latter during the sea-level low stand of the last glacial stage. During this period M. monspessulanus shows a sudden demographic expansion, associated with increased habitat suitability in North Africa. Lower habitat suitability is predicted for both species during interglacial stages, with suitable areas restricted to coastal and mountain ranges of Iberia and Morocco. Compiled fossil data for M. monspessulanus show a continuous fossil record in Iberia at least since the Pliocene and throughout the Pleistocene.

CONCLUSIONS

The previously proposed hypothesis of Pleistocene glacial extinction of both species in Europe is not supported based on genetic data, bioclimatic envelopes models, and the available fossil record. A model of range retraction to mountain refugia during arid periods and of glacial expansion (demographic and spatial) associated to an increase of Mediterranean habitats during glacial epochs emerges as a general pattern for mesic vertebrates in North Africa. Moreover, the phylogeographic pattern of H. hippocrepis conforms to a well-established biogeographic partition between western and eastern Maghreb.

Impact of model assumptions on demographic inferences: the case study of two sympatric mouse lemurs in northwestern Madagascar

BACKGROUND

Quaternary climate fluctuations have been acknowledged as major drivers of the geographical distribution of the extraordinary biodiversity observed in tropical biomes, including Madagascar. The main existing framework for Pleistocene Malagasy diversification assumes that forest cover was strongly shaped by warmer Interglacials (leading to forest expansion) and by cooler and arid glacials (leading to forest contraction), but predictions derived from this scenario for forest-dwelling animals have rarely been tested with genomic datasets.

RESULTS

We generated genomic data and applied three complementary demographic approaches (Stairway Plot, PSMC and IICR-simulations) to infer population size and connectivity changes for two forest-dependent primate species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar. The analyses suggested major demographic changes in both species that could be interpreted in two ways, depending on underlying model assumptions (i.e., panmixia or population structure). Under panmixia, the two species exhibited larger population sizes across the Last Glacial Maximum (LGM) and towards the African Humid Period (AHP). This peak was followed by a population decline in M. ravelobensis until the present, while M. murinus may have experienced a second population expansion that was followed by a sharp decline starting 3000 years ago. In contrast, simulations under population structure suggested decreasing population connectivity between the Last Interglacial and the LGM for both species, but increased connectivity during the AHP exclusively for M. murinus.

CONCLUSION

Our study shows that closely related species may differ in their responses to climatic events. Assuming that Pleistocene climatic conditions in the lowlands were similar to those in the Malagasy highlands, some demographic dynamics would be better explained by changes in population connectivity than in population size. However, changes in connectivity alone cannot be easily reconciled with a founder effect that was shown for M. murinus during its colonization of the northwestern Madagascar in the late Pleistocene. To decide between the two alternative models, more knowledge about historic forest dynamics in lowland habitats is necessary. Altogether, our study stresses that demographic inferences strongly depend on the underlying model assumptions. Final conclusions should therefore be based on a comparative evaluation of multiple approaches.

Past environmental changes affected lemur population dynamics prior to human impact in Madagascar

Quaternary climatic changes have been invoked as important drivers of species diversification worldwide. However, the impact of such changes on vegetation and animal population dynamics in tropical regions remains debated. To overcome this uncertainty, we integrated high-resolution paleoenvironmental reconstructions from a sedimentary record covering the past 25,000 years with demographic inferences of a forest-dwelling primate species (Microcebus arnholdi), in northern Madagascar. Result comparisons suggest that climate changes through the African Humid Period (15.2 - 5.5 kyr) strongly affected the demographic dynamics of M. arnholdi. We further inferred a population decline in the last millennium which was likely shaped by the combination of climatic and anthropogenic impacts. Our findings demonstrate that population fluctuations in Malagasy wildlife were substantial prior to a significant human impact. This provides a critical knowledge of climatically driven, environmental and ecological changes in the past, which is essential to better understand the dynamics and resilience of current biodiversity.

Macroevolutionary variation and environmental correlates of scalation traits in Eurasian vipers (Serpentes: Viperinae)

Understanding how phenotypic variation across species is shaped by the combination of shared evolutionary history and environmental factors is key to elucidating the processes that underlie biodiversity. In reptiles, morphological traits have traditionally been used to delimit species and make systematic inferences. Recent studies highlight the possibility that phenotypic variation, particularly in scalation traits, might instead be driven by environmental factors and therefore not reflect the phylogenetic relationships among species. In this study, we combined morphological and ecological data in a macroevolutionary framework, in order to describe the morphological variation across species of Eurasian vipers (Serpentes: Viperinae), investigate the phylogenetic structure of scalation traits and test the contribution of environmental factors in shaping morphological patterns. We found considerable variation in all examined traits, which, in most cases, agreed with the phylogenetic relationships among species, reinforcing their usefulness for taxonomic inferences. Interestingly, however, the number of ventral scales exhibited lower phylogenetic signal and a tight association with environmental factors of geographical ranges, suggesting potential adaptive or developmental sources of variation in the trait. This is the first comparative study of macroevolutionary variation in scalation traits in Eurasian vipers, validating the use of most of them for systematic inferences, but also indicating possible environmental factors that might shape phenotypic variation across species.

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.

Thermal melanism explains macroevolutionary variation of dorsal pigmentation in Eurasian vipers

Colouration may endorse thermoregulatory and antipredatory functions in snakes. The thermal melanism hypothesis predicts that dark-coloured individuals are ecologically favoured in cool climates. However, the loss of aposematic and cryptic colourations may imply high predation for melanistic snakes. Here, we used the monophyletic group of Eurasian vipers (subfamily Viperinae) to test whether an increase in the extent of dark area inside the characteristic zigzag dorsal pattern is associated to colder environments. We measured two colouration traits in zigzag-patterned individuals (number of dorsal marks and weighted pigmentation index) and used a phylogenetic comparative approach to explore macroevolutionary patterns of dorsal pigmentation and test whether its extent is associated to ecogeographic characteristics of lineages’ ranges. Phylogenetically-naïve and phylogenetically-informed analyses yielded a significant association between the degree of pigmentation of the zigzag pattern and environmental variables such as solar radiation, elevation and latitude. The degree of pigmentation of the zigzag pattern is highlighted as an adaptive trait that matches range attributes mirroring cold environments irrespective of the phylogeny. These results constitute the first large-scale evidence supporting the thermal melanism hypothesis in snakes, opening new avenues of inquiry for the mechanisms that shape the evolution of colour phenotypes.

Origin, extinction and ancient DNA of a new fossil insular viper: molecular clues of overseas immigration

Viperinae is a subfamily of viperid snakes whose fossil record in the Mediterranean islands is, until now, restricted to 12 palaeontological deposits on seven islands. Revision of the material excavated 30 years ago from the Middle/Late Pleistocene–Holocene deposit of Es Pouàs [Eivissa (= Ibiza), Balearic Islands, western Mediterranean] revealed about 6000 bones of a small-sized viper across different stratigraphic levels. Its morphological characteristics are different enough to known species of Vipera to warrant the description of a new species, but the nearly complete mitochondrial genome obtained from this snake based on a sample dated to 16 130 ± 45 bp, suggested it belonged to a new insular population of Lataste’s viper (Vipera latastei), Vipera latastei ebusitana subsp. nov. Phylogenetic analysis indicates that the dispersal of the ancestors of V. l. ebusitana to Eivissa, most probably from a north-east Iberian population, occurred via overwater colonization < 1.5 Mya, well after the Messinian Salinity Crisis (5.97–5.32 Mya) when land bridges allowed terrestrial colonization of the Balearic Islands by mainland faunas. The morphological differences between V. l. ebusitana and the Iberian populations suggest that it is a new dwarf taxon resulting from insular evolutionary processes, becoming extinct shortly after the first human arrival to this island about 4000 years ago.

Chasing the phantom: biogeography and conservation of Vipera latastei-monticola in the Maghreb (North Africa)

The Maghreb region (North Africa) constitutes a major component of the Mediterranean Basin biodiversity hotspot. During the last centuries, a consistent human population growth has led to an unprecedented rate of habitat transformation and loss in the region and thus, threatening its biodiversity. The Western Mediterranean viper Vipera latastei-monticola inhabits humid and subhumid areas in the main mountain ranges of the Maghreb, facing such threatening factors; however, its elusive character and rarity hindered data collection for distinct biological purposes. Here, we study the biogeographical patterns and conservation status of the Maghrebian V. latastei-monticola resulting from recent sampling campaigns in Morocco and Tunisia. We update species distribution, and integrate phylogeographic and ecological niche modelling analyses at both species and lineage level to identify suitable areas, and to evaluate the impact of anthropogenic transformation and level of protection of their suitable space. We identified four highly divergent mitochondrial lineages, including a new lineage endemic to the Western High Atlas, with allopatric distributions and restricted to mountain ranges, supporting the role of mountains as past climatic refugia. Despite the remoteness of suitable areas, we report widespread habitat degradation and identify the low effectiveness of the current protected areas system in preserving the species and lineages range. Our study shows the urgent need to apply management actions for the long-term conservation of this vulnerable species and suggests a revaluation of the specific status of V. monticola, as these populations likely represent an ecotype of V. latastei.