Pleistocene diversification in Morocco and recent demographic expansion in the Mediterranean pond turtleMauremys leprosa

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

RADseq Data Suggest Occasional Hybridization between Microcebus murinus and M. ravelobensis in Northwestern Madagascar

The occurrence of natural hybridization has been reported in a wide range of organisms, including primates. The present study focuses on the endemic lemurs of Madagascar, primates for which only a few species occur in sympatry or parapatry with congeners, thereby creating limited opportunity for natural hybridization. This study examines RADseq data from 480 individuals to investigate whether the recent expansion of Microcebus murinus towards the northwest and subsequent secondary contact with Microcebus ravelobensis has resulted in the occurrence of hybridization between the two species. Admixture analysis identified one individual with 26% of nuclear admixture, which may correspond to an F2- or F3-hybrid. A composite-likelihood approach was subsequently used to test the fit of alternative phylogeographic scenarios to the genomic data and to date introgression. The simulations yielded support for low levels of gene flow (2Nm0 = 0.063) between the two species starting before the Last Glacial Maximum (between 54 and 142 kyr). Since M. murinus most likely colonized northwestern Madagascar during the Late Pleistocene, the rather recent secondary contact with M. ravelobensis has likely created the opportunity for occasional hybridization. Although reproductive isolation between these distantly related congeners is not complete, it is effective in maintaining species boundaries.

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.

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.

Characterization and identification of haemogregarine hemoparasites (Apicomplexa: Adeleina: Hepatozoidae) in natural populations of Mauremys leprosa leprosa and M. leprosa saharica from Morocco

Distribution, prevalence and parasitaemia of apicomplexan parasites of the genus Haemogregarina were studied in 858 freshwater turtles (735 Mauremys leprosa leprosa and 123 Mauremys leprosa saharica) throughout 30 localities from Morocco. Blood smears were collected from the turtles and partial 18S rRNA sequence data used to infer genetic diversity and phylogenetic relationships. Of the 858 individuals analyzed by microscopy 22.7% were infected, from 16 M. leprosa leprosa localities. Individuals of M. leprosa saharica sampled south of the High Atlas Mountains were not infected, probably due to the absence or rarity of the leech vectors in these localities. Within M. leprosa leprosa, we did not identify any patterns between geography and prevalence, which varied between 10% (Oued Nfiss) to 100% (Oued Tassaout). Five distinct genetic lineages were identified, indicating the likely presence of multiple species of haemogregarines, one of which probably corresponds to Haemogregarina stepanowi. Mixed infections were also detected. Additional studies are needed to better understand the ecology and epidemiology of this parasite in turtles, as well as the host-parasite relationship with their definitive hosts, leeches.

Population structure and gene flow of the syntopic turtles Emys and Mauremys from coastal and inland regions of Anatolia (Turkey): results from mitochondrial and microsatellite data

Revealing the genetic basis of the existence of different species living together in different geographic regions provides clarification of this phylogeographic differentiation. In this study, we investigated the population genetics and evaluated the level of genetic variation of inland and coastal populations of Mauremys and Emys in Turkey. Tissue samples of 196 terrapins were studied which were collected from syntopic coastal (Gölbent-Söke/Aydın; M. rivulata and E. orbicularis) and inland populations (Bahçesaray/Aksaray; M. caspica and E. orbicularis). DNA was isolated using the InnuPREP DNA Mini Kit. Mitochondrial DNA sequences and allelic variation at 13 microsatellite loci for Mauremys and 12 microsatellite loci for Emys were examined.  Three haplotypes were found for Emys orbicularis (Im, Ip and Iw) collected from the coastal region and two haplotypes for Emys orbicularis (Ig and Im) collected from inland. Two haplotypes were identified for M. caspica (Cmt8 and Cmt9) and three haplotypes were identified for M. rivulata (Rmt3, Rmt24 and Rmt26). Using microsatellites and the software STRUCTURE the most probable value for K was revealed as two 2 for both species. The F_ST value between M. rivulata and M. caspica was 0.39, and between the coastal and inland populations of E. orbicularis 0.09. It can be concluded that Emys populations tend to evolve by somehow preserving the allelic richness they have and Mauremys populations continue to differentiate so that new species emerge in the evolutionary process to reach the ideal allelic structure.

Population structure and body size of the Sahara blue-eyed pond turtle Mauremys leprosa saharica, from an isolated pond in the Lower Draa Basin, southern Morocco

The marginal populations of the Sahara blue-eyed pond turtle, Mauremys leprosa saharica (Testudines: Geoemydidae) in the pre-Saharan area in the Lower Draa Valley, are faced to extreme environmental conditions of arid climate and anthropogenic activities including overuse of water and land salinization. In this study, we investigated a small isolated population of the blue-eyed pond turtle at Sidi El Mehdaoui Oasis, Lower Draa Valley, to determine its population structure and morphometric characteristics. Water salinity was about 24% of seawater and dissolved oxygen concentration was less than 3 mg l-1 (sub-hypoxia). The average body size of adults (carapace length) was much lower in comparison to other Moroccan populations. Among adult individuals, the mean carapace lengths were 103 ± 21 mm and 104.4 ± 35 mm for males and females, respectively. The corresponding mean body mass were 173 ± 92.4 g and 135 ± 52 g. There were no significant differences between sexes for size or weight. The studied population is small, with a predominance of adults (75%). A significant sexual size dimorphism was found. In fact, the Sexual Dimorphism Index for size had the lowest level ever described for the species possibly because of the low productivity of this marginal isolated habitat. The increased drought and salinization are major threats to the long-term persistence of the vulnerable populations of the Saharan blue-eyed pond turtle and their habitats. Therefore, conservation measures of these populations through the protection of their habitats are very urgent.

Reconstructing evolution at the community level: A case study on Mediterranean amphibians

Reconstructing reliable timescales for species evolution is an important and indispensable goal of modern biogeography. However, many factors influence the estimation of divergence times, and uncertainty in the inferred time trees remains a major issue that is often insufficiently acknowledged. We here focus on a fundamental problem of time tree analysis: the combination of slow-evolving (nuclear DNA) and fast-evolving (mitochondrial DNA) markers in a single time tree. Both markers differ in their suitability to infer divergences at different time scales (the 'genome-timescale-dilemma'). However, strategies to infer shallow and deep divergences in a single time tree have rarely been compared empirically. Using Mediterranean amphibians as model system that is exceptional in its geographic and taxonomic completeness of available genetic information, we analyze 202 lineages of western Palearctic amphibians across the entire Mediterranean region. We compiled data of four nuclear and five mitochondrial genes and used twelve fossil calibration points widely acknowledged for amphibian evolution. We reconstruct time trees for an extensive lineage-level data set and compare the performances of the different trees: the first tree is based on primary fossil calibration and mitochondrial DNA, while the second tree is based on a combination of primary fossil and on secondary calibrations taken from a nuclear tree using mitochondrial DNA (two-step protocol). Focusing on a set of nodes that are most likely explained by vicariance, we statistically compare the reconstructed alternative time trees by applying a biogeographical plausibility test. Our two-step protocol outperformed the alternative approach in terms of spatial and temporal plausibility. It allows us to infer scenarios for Mediterranean amphibian evolution in eight geographic provinces. We identified several tectonic and climatic events explaining the majority of Mediterranean amphibian divergences, with Plio-Pleistocene climatic fluctuations being the dominant driver for intrageneric evolution. However, often more than one event could be invoked for a specific split. We give recommendations for the use of secondary calibrations in future molecular clock analyses at the community level.

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.