# Comparative phylogeography of amphibians and reptiles in Algeria suggests common causes for the east-west phylogeographic breaks in the Maghreb

Menad Beddek, Saliha Zenboudji-Beddek, Philippe Geniez, Raouaa Fathalla, Patricia Sourouille, Véronique Arnal, Boualem Dellaoui, Fatiha Koudache, Salah Telailia, Olivier Peyre, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Menad Beddek, Saliha Zenboudji-Beddek, Philippe Geniez, Raouaa Fathalla, Patricia Sourouille, et al. Comparative phylogeography of amphibians and reptiles in Algeria suggests common causes for the east-west phylogeographic breaks in the Maghreb. PLoS ONE, 2018, 13 (8), e02626320. [https://doi.org/10.1371/journal.pone.0201218](https://doi.org/10.1371/journal.pone.0201218).

HAL Id: hal-02626320

https://hal.inrae.fr/hal-02626320

Submitted on 26 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# Comparative phylogeography of amphibians and reptiles in Algeria suggests common causes for the east-west phylogeographic breaks in the Maghreb

Menad Beddek¹,²*, Saliha Zenboudji-Beddek¹, Philippe Geniez¹, Raouaa Fathalla¹,³, Patricia Sourouille¹, Véronique Arnal¹, Boualem Dellaoui⁴, Fatiha Koudache⁴, Salah Telailia⁵, Olivier Peyre², Pierre-André Crochet¹

1. CEFE, CNRS, Univ. Montpellier, Univ Paul Valéry Montpellier 3, INRA, EPHE, IRD, Montpellier, France,

2. Naturalia Environnement, Site Agroparc, Avignon, France, 3. Unité de recherche de biodiversité et biologie des populations, Faculté des sciences de Tunis, Université de Tunis Elmanar, Tunis, Tunisie, 4. Département de l'Environnement, Faculté des Sciences, Université Djillali Liabes, Sidi Bel Abbes, Algérie, 5. Département des Sciences Agronomiques, Faculté des Sciences de la Nature et de la Vie, Université Chadli Bendjedid, El Tarf, Algérie

* menad-beddek@hotmail.com

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Beddek M, Zenboudji-Beddek S, Geniez P, Fathalla R, Sourouille P, Arnal V, et al. (2018) Comparative phylogeography of amphibians and reptiles in Algeria suggests common causes for the east-west phylogeographic breaks in the Maghreb. PLoS ONE 13(8): e0201218. https://doi.org/10.1371/journal.pone.0201218

**Editor:** Ulrich Joger, State Museum of Natural History, GERMANY

**Received:** October 9, 2017

**Accepted:** July 11, 2018

**Published:** August 29, 2018

**Copyright:** © 2018 Beddek et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All the genetic sequences produced for this paper are available on GenBank with accession numbers available in the Supporting Information files.

**Funding:** Menad BEDDEK received a Ph.D grant from Naturalia-Environnement/ANRT CIFRE convention N° 2013/0145. Naturalia Environnement: http://www.naturalia-environnement.fr/; ANRT CIFRE: http://www.anrt.asso.fr/fr/cifre-7843. The funder "Naturalia

## Abstract

A series of phylogeographic studies in the Maghreb identified a repeated pattern of deep genetic divergence between an eastern (Tunisia) and western (Morocco) lineage for several taxa but lack of sampling in Algeria made it difficult to know if the range limits between the eastern and western lineages were shared among taxa or not. To address this question, we designed a comparative phylogeographic study using 8 reptile and 3 amphibian species with wide distribution in the Maghreb as models. We selected species where previous studies had identified an East-West phylogeographic divide and collected sampled in Algeria to 1) examine whether the simple East-West divergence pattern still holds after filling the sampling gap in Algeria or if more complex diversity patterns emerge; 2) if the E-W pattern still holds, test whether the limits between the E and W clades are shared between species, suggesting that common historical process caused the E-W divergences; 3) if E-W limits are shared between species, use information on the age of the divergence to identify possible geological or climatic events that could have triggered these E-W differentiations. We found that the E-W pattern was generally maintained after additional sampling in Algeria and identified two common disjunction areas, one around the Algeria-Morocco border, the other one in Kabylia (central Algeria), suggesting that common historical mechanisms caused the E-W divergences in the Maghreb. Our estimates for the times to most common recent ancestors to the E and W clades span a wide range between the Messinian salinity crisis and the Plio-Pleistocene limit (except for one older split), suggesting different origins for the initial divergences and subsequent preservation of the E and W lineages in common climatic refugia in the west and the east of the Maghreb.
Environnement" represented by Olivier Peyre has
participated to the study design and manuscript
preparation.

**Competing interests:** The authors have declared
that no competing interests exist.

**Introduction**

Understanding how global biodiversity patterns emerged from historical and ecological pro-
cesses has been a major task of biogeographical studies ever since they started; identifying
areas of higher diversity and endemism is also the basis for applied conservation [1,2]. The
Mediterranean Basin constitutes a global biodiversity hotspot [3] and as a consequence has
been identified as a conservation priority [4] because it contains high rates of endemism and
elevated species richness in several taxa: higher plants [3,5], freshwater fishes [6], mammals
[7], reptiles and amphibians [8], dragonflies [9].

Within the Mediterranean Basin, the Maghreb (here defined as Northern Algeria, Tunisia, and Morocco) is a well-defined ecological unit surrounded by natural barriers: the Sahara to the south and the east, the Mediterranean to the north and the Atlantic Ocean to the west [10]. It contains several plant diversity hotspots [11,12] and is inhabited by many endemic species of reptiles and amphibians [8,13]. There is however a huge discrepancy between the amount of knowledge on biodiversity in the northern and the southern sides of the Mediterranean basin [14,15] and further research in the Maghreb would allow a better understanding of the historical processes that shaped Mediterranean biodiversity at a regional scale.

Many studies have attempted to understand how the Mediterranean biodiversity hotspot was generated (for ex. [16–20]). One of the emerging patterns of these studies is that for many species North African and South European populations are genetically related, demonstrating that the Mediterranean Sea has not always prevented dispersal and gene flow between Europe and North Africa. Much interchange occurred via land connexion before the opening of Gibraltar (see references in [19]) but in some species dispersal occurred after the filing of the Mediterranean Basin and the opening of Gibraltar at the end of the Messinian crisis (e.g. [21,22]). This succession of divergence on either side of the Mediterranean Sea and colonisation across the Mediterranean contributed building biodiversity in the Mediterranean Basin, and the Maghreb has been the source of colonisation of Europe in many thermophilic organisms [23–25].

Another factor explaining the high diversity of the Mediterranean Basin is its complex
topography and history that promoted elevated rates of in situ diversification and allowed per-
sistence of this diversity in multiple refugia during cold periods of the glacial cycles. This phe-
nomenon has been particularly well studied in the north of the region (southern Europe: see
[26]) but the Maghreb has not been as intensively investigated. Like southern Europe, the
Maghreb encompasses a large diversity of landscapes and bioclimatic regions [27] and a series
of phylogeographic studies also identified a high level genetic variability in many species ([25]
and references therein). However, we still lack a comprehensive understanding of the main
drivers of diversification within the Maghreb.

One of the shortcoming that currently impedes understanding of how the biodiversity pat-
terns in North Africa emerged is the lack of a multi-taxa phylogeography approach. As advo-
cated by [28], a comparative phylogeography approach is needed to identify relevant historical
causes explaining observed diversity patterns. Because many concurrent demographic, geologi-
cal and /or climatic events can shape the current distribution of any given taxon, only “the com-
parison of geographical patterns of genetic variation among multiple co-distributed taxa” [29]
can enlighten demographic and historical processes responsible for the building up of genetic
diversity within species and between closely related species [30]. Accordingly, a common
genetic diversity pattern shared in a given area between several taxa is supposed to be shaped by
the same historical events. For example, location of contact zones between divergent lineages in
Europe exhibit common patterns that have been linked to common glacial refugia and post-gla-
cial colonization routes [31–33]. No such studies are available for the Maghreb, but several bio-
diversity hotspots have been suggested to correspond with glacial refugia [20,25].
While we still lack a proper comparative phylogeography approach in the Maghreb, the existing phylogeographic studies identify a seemingly repeated geographical pattern of deep lineages divergence between Tunisia in the east and Morocco in the west of the region. This E-W divergence is verified in the majority of the studied taxa: reptiles (Natrix maura, [34,35], Chalcides ocellatus, [36], Trogonophis wiegmanni [37], Ptyodactylus oudrii, [38,39], Timon tan-gitanus / pater, [40,41], amphibians (Hyla meridionalis, [42,43], Pelophylax saharicus, [44]), scorpions [45], rodents (Jaculus orientalis, [46], Meriones shawii complex, [47]), molluscs: Cornu aspersum [48]. In birds the pattern does not seem to be as general but it has been recov-ered in Dupont's Lark [49] and Crested Lark [50] at least.

This common pattern is currently difficult to interpret because in most cases there are few or no sample from Algeria. As a consequence, the locations of the suture zones between the E and W clades are unknown. The reasons for the lack of sampling in Algeria are mainly the cur-rent political unrest and safety issues, resulting in a large sampling gap that currently precludes our understanding of the origin of the E-W genetic divide in Maghreban lineages.

In order to improve our knowledge on the origin of current biodiversity patterns in the Maghreb, we designed a comparative phylogeography approach in the Maghreb, specifically aiming at reducing the sampling gap in Algeria. We propose that reptiles and amphibians con-stitute a relevant model for this study, firstly because previous studies allowed us to select several taxa showing deep E-W genetic divergence, secondly because their relatively lower dispersal ability relative to birds or large mammals make them more sensitive to landscape barriers or cli-matic oscillations. The aim of this study were thus: 1) to examine whether the simple East-West divergence pattern that divides the populations of many species into two main clades distrib-uted east and west of the Maghreb still holds after filling the sampling gap in Algeria or if more complex diversity patterns with different geographical distributions emerge; 2) if the E-W pat-tern still holds, to test whether the limits between the E and W clades are shared between spe-cies, suggesting that common historical process caused the E-W divergences; 3) if E-W limits are shared between species, to use information on the age of the divergence to identify possible geological or climatic events that could have triggered these E-W differentiations.

To carry out this study, we selected species that fulfil two criteria: 1) they occur throughout most of the Maghreb and 2) they have been the subject of previous phylogeographic studies and exhibit deep divergence between an eastern and western clade. We conducted targeted sampling in Algeria to produce original sequence data that were merged with sequence avail-able in Genbank. We use mitochondrial DNA (mtDNA) as a proxy for evolutionary units because it allowed us to investigate a larger number of species (11, see material and methods) with our available resources. We are fully aware of the risks of single-markers studies [51] and of the benefits of multi-locus approaches but multi-locus analyses would have constrained us to use a reduced number of species. Moreover, most of former studies were mainly or exclu-sively based on mtDNA, so mtDNA is the only way to merge all existing data together. Even if several evolutionary mechanisms can cause discordance in divergence patterns between mtDNA and the nuclear genome [52], the most common pattern is by far general congruence between the two types of markers [53,54], making us confident that our general conclusion will not be affected by this problem.

**Material and methods**

**Ethic statement**

Specimens from Morocco were collected under the authorisation n°15 HCEFLCD / DLCDPN / DPRN / CFF delivered by the « Haut-Commissariat aux Eaux et Forêts et à la Lutte Contre la Désertification ». The list of protected species of reptiles and amphibians in Algeria is has been
determined by the Executive Decree n°12-235 05/24/2012. Among the species we have used, the following are protected under Algerian law: *Acanthodactylus savignyi*, *A. blanci*, *Chalcides ocellatus*, *Timon pater*, *T. tangitanus*. However, the legislative authority supposed to deliver capture licenses for protected species according to the Law n°14-07 relative to Natural Resources is not identified in any official document and Algerian biologists therefore have to collect tissue samples or voucher specimens without legal permits. Most tissue samples used for this study have been obtained by non-invasive processes (tail tips or muscle samples from road-killed specimens). Voucher collection has followed practice of Algerian herpetologists. Field work in protected areas has been approved, and often supported in the field, by local authorities (see list in Acknowledgments).

## Sampling

A total of 323 tissue samples (mainly muscle samples, sometimes buccal swabs) mostly collected from the field in Northern Algeria from 2013 to 2015 and deposited in the CEFE-CNRS tissues collection in Montpellier (BEV collection) were used in this study (S1 Table). The samples belong to 11 taxa, 8 reptiles and 3 amphibians. All these taxa had been subjected to previous phylogeographic studies revealing an E-W divergence in mtDNA: *Acanthodactylus erythrurus* [3], *Chalcides ocellatus* [36], *Hemorrhois hippocrepis* [55], *Natrix maura* [34], *Podarcis vaucheri* [56-58], *Ptyodactylus oudrii*, [38,39], *Timon pater* / *tangitanus* [40,41]. The samples origin and their GenBank accession codes are detailed in S1 Table.

## DNA extraction, amplification and sequencing

DNA was extracted with the Qiagen DNeasy Blood and Tissues kit following the manufacturer's instructions. Standard protocols were used for Polymerase Chain Reaction (PCR) amplifications with the Sigma kit. The choice of the molecular marker sequenced for each species was dependent on previous published studies for this species. We used sequences of ND4 and adjacent tRNA for *Acanthodactylus erythrurus*, *Podarcis vaucheri* and *Natrix maura*, CytB for *Chalcides ocellatus*, *Hemorrhois hippocrepis* and *Discoglossus pictus*, CO1 for *Hyla meridionalis* and *Pelophylax saharicus*, 16S RNA for *Timon pater* / *tangitanus* and 12S RNA for *Ptyodactylus oudrii*, allowing us to merge our sequences with the sequences already published. PCR conditions and primers are detailed in S2 Table. Purification and sequencing of PCR products were carried out by Eurofins MWG Germany. The chromatograms were checked using Codon code aligner version 4.2.7 for Windows (CodonCode Corporation, Dedham, MA, USA). The alignment of the sequences was carried out with KlustalW implemented in Mega version 7.0 [59] and verified by eye.

## Phylogenetic analyses

Rooted trees were produced with two approaches. Firstly, RaxML [60] with partitioning and the online platform PhyML [61] with automatic AIC criterion based substitution model selection and NNI option for tree improvement were used for maximum likelihood reconstructions. The support of the nodes was evaluated by bootstrapping with 1000 replicates. As GTR is the only nucleotide substitution model available in RAxML, GTR + I + G was applied to all partitions. Secondly, MrBayes 3.2 [62] was used for Bayesian inference (BI). Two separate runs of 6 million generations were conducted with a sampling every 500 trees. Four Monte Carlo Markov chains were initiated for which the starting topology was randomly selected. The convergence of the runs was visualized using Tracer 1.4 [63]. The first 10% trees were eliminated as the burn-in period. The majority consensus tree was computed from the remaining trees. For BI methods, best evolution models and partitioning were calculated using Partitionfinder
[64]. Fifty-six different models were tested and the Bayesian information criterion (BIC) was used for the final choice of model. Both reconstruction methods gave similar topologies for all taxa except *Timon* and we retained the BI topology for inference and mapping of the deepest lineages (see below). For *Timon* we did not manage to obtain the same topology with and without our short sequences so we inferred the topology using long sequences from Genbank only and placed the new Algerian lineage (only represented by our short sequences) in this tree in a second step. For *P. oudrii* and *T. wiegmanni*, we added longer Genbank sequences to the main alignment comprising our sequences (12 S for *P. oudrii*, 16 S for *T. wiegmanni*, see S1 Table) to improve the nodal robustness.

## Time calibration

The time to most recent common ancestor (TMRCA) of the eastern and western mitochondrial clades were estimated as follows: time calibrated trees were obtained with Beast v1.8.2 [65] under an uncorrelated lognormal relaxed clock; we used a Yule prior on rates of cladogenesis for almost all taxa since the trees included outgroup species but the constant size prior was used for *P. saharicus*; four runs were conducted for each taxon (parameters of chain's lengths and substitution rates are given in Table 1). For each time calibration tree, we used published substitution rates: [66] for *A. erythrurus*, medium evolution rate (see [67]) for *H. meridionalis* and *D. pictus*, [39] for *P. oudrii*, [68] for *P. vaucheri*, [69] for *T. pater* and *T. tangitanus*, [70] for *P. saharicus*, except *N. maura* for which we used published time calibrated nodes [71]. The four runs were first analysed in Tracer 1.4 [63] to check the convergence of the chains and then combined in LogCombiner v1.8.2 [72]. TreeAnnotator v1.7.0 [73] was used to produce the chronograms. We performed the dating analysis only on 8 of the 11 taxa (see Table 1). For *C. ocellatus* the analysis was already done by [36]. For *H. hippocrepis* and *T. wiegmanni* because of non-convergence of the runs we did not manage to have accurate results even after several runs (>4).

## Lineage delimitation and genetic diversity

To estimate genetic diversity for each species we used Dnasp [74] to compute the following indices: S: number of segregating sites, H: number of haplotypes, Hd: haplotypic diversity. To assess if an E-W divide is the main process shaping mitochondrial diversity in our taxa, we used the deepest split in the phylogeny of every taxon to identify the two main clades and mapped the distribution of these two clades. To explore further significant structure in mitochondrial diversity, we used the automatic barcode gap discovery (ABGD) method [75]. Although this program is used for species delimitation, we use it to objectively describe distinct evolutionary units regardless to their taxonomic rank. The ABGD approach is based on the principle that

<table><caption>Table 1. Substitution models, tree prior, mean substitution rate and the number of Monte Carlo Markov chains used in the reconstruction of time calibrated trees.</caption><thead><tr><td>Taxon</td><td>Substitution Model</td><td>Tree prior</td><td>Mean Rate</td><td>MCMC</td></tr></thead><tbody><tr><td>Acanthodactylus erythrurus</td><td>GTR</td><td>Yule process</td><td>0.0174</td><td>2.5X10<sup>7</sup></td></tr><tr><td>Natrix maura</td><td>TN93+G+I</td><td>Yule process</td><td>-</td><td>2X10<sup>7</sup></td></tr><tr><td>Podarcis vaucheri</td><td>TN93+G+I</td><td>Yule process</td><td>0.0174</td><td>2X10<sup>7</sup></td></tr><tr><td>Ptyodactylus oudrii</td><td>GTR+G</td><td>Yule process</td><td>0.00755</td><td>2X10<sup>7</sup></td></tr><tr><td>Timon spp</td><td>GTR</td><td>Yule process</td><td>0.0016</td><td>2X10<sup>7</sup></td></tr><tr><td>Discoglossus pictus</td><td>TN93+G+I</td><td>Yule process</td><td>0.00825</td><td>2X10<sup>7</sup></td></tr><tr><td>Hyla meridionalis</td><td>TN93+G+I</td><td>Yule process</td><td>0.0085</td><td>2X10<sup>7</sup></td></tr><tr><td>Pelophylax saharicus</td><td>HKY</td><td>Constant size</td><td>0.0195</td><td>2X10<sup>7</sup></td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0201218.t001
divergence between individuals belonging to the same “species” (here lineage) is smaller than divergence between individuals from different “species”. This approach is designed to work with short DNA sequences, hence it is appropriate for our short sequence data. The analyses were run with default settings. However, as recommended by the program’s author when the run with the default gap width gives only one partition, the relative gap width (X = 1.5 by default) was decreased to 1.2 for *Trogonophis wiegmanni* and to 1.0 for *Acanthodactylus erythrurus*. The most stable distribution value of the initial partition was used to decide the most likely number of distinct lineages. We spatially displayed the results using ARC map10 to get an

<figure><img src="image_5.png" /><figcaption>**Fig 1.** Geographic distribution of the reptiles' genetic lineages delimited by ABGD program. The thick black lines correspond to the suture zones of the deepest splits for each taxon. Dashed lines correspond to the geographical limits between lineages delimited by ABGD program. The ABGD lineages are represented by the coloured dots.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g001
explicit map of lineage distribution, then, the geographical limits between linages were plotted using the mid distance between the nearest individuals of adjacent lineages.

# Results

## Mitochondrial diversity and geographical distribution of major clades

Delimitation of mtDNA clades by ABGD was concordant with previous studies but our additional sampling in Algeria revealed previously overlooked lineages in several taxa (*Timon spp*, *H. meridionalis*, *A. erythrurus*, *P. oudrii*, *P. vaucheri*: compare Figs 1 and 2 with previous

<figure><img src="image_6.png" /><figcaption>**Fig 2.** Geographic distribution of the amphibians' genetic lineages delimited by ABGD program. The thick black lines correspond to the suture zones of the deepest splits for each taxon. Dashed lines correspond to the geographical limits between lineages delimited by ABGD program. The ABGD lineages are represented by the coloured dots.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g002
<table><caption>Table 2. Number of groups inferred by ABGD program, the sequences lengths, genetic distances between and within ABGD groups.</caption><thead><tr><th>Taxon</th><th>ABGD</th><th>Length (pb)</th><th>Distance within group</th><th>Distance between groups</th></tr></thead><tbody><tr><td><i>Acanthodactylus erythrurus</i></td><td>15</td><td>776</td><td>0-0.0156-0.0297</td><td>0.068-0.125-0.166</td></tr><tr><td><i>Chalcides ocellatus</i></td><td>2</td><td>425</td><td>0.0286-0.0338-0.039</td><td>0.102-Nc-0.102</td></tr><tr><td><i>Hemorrhois hippocrepis</i></td><td>2</td><td>641</td><td>0.0042-0.0046-0.0049</td><td>0.029-Nc-0.029</td></tr><tr><td><i>Natrix maura</i></td><td>5</td><td>604</td><td>0-0.0023-0.0052</td><td>0.025-0.045-0.057</td></tr><tr><td><i>Podarcis vaucheri</i></td><td>6</td><td>659</td><td>0.001-0.0129-0.0294</td><td>0.064-0.096-0.134</td></tr><tr><td><i>Ptyodactylus oudrii</i></td><td>9</td><td>454</td><td>0.0033-0.0107-0.0268</td><td>0.034-0.117-0.177</td></tr><tr><td><i>Timon pater / tangitanus</i></td><td>5</td><td>477</td><td>0.0009-0.0048-0.123</td><td>0.021-0.042-0.057</td></tr><tr><td><i>Trogonophis wiegmanni</i></td><td>3</td><td>474</td><td>0.0123-0.0148-0.01814</td><td>0.0123-0.0148-0.0446</td></tr><tr><td><i>Discoglossus pictus</i></td><td>2</td><td>964</td><td>0.007-0.008-0.009</td><td>0.021-Nc-0.021</td></tr><tr><td><i>Hyla meridionalis</i></td><td>4</td><td>775</td><td>0.0011-0.0049-0.0096</td><td>0.012-0.047-0.076</td></tr><tr><td><i>Pelophylax saharicus</i></td><td>3</td><td>654</td><td>0.0034-0.005-0.0068</td><td>0.072-0.108-0.13</td></tr></tbody></table>

From left to right the minimum, medium and maximum distances are given within and between groups.

https://doi.org/10.1371/journal.pone.0201218.t002

<figure><img src="image_7.png" /><figcaption>Fig 3. *Acanthodactylus erythrurus* Bayesian inference phylogenetic tree and geographic distribution of deep clades. The green area in the map corresponds to the distribution of the eastern clade populations and the red area corresponds to the distribution of the western clade populations. The blue lineage is missing nodal support, its phylogenetic position is not resolved with the present data set.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g003
<figure><img src="image_8.png" /><figcaption>Fig 4. *Chalcides ocellatus* Bayesian inference phylogenetic tree and geographic distribution of deep clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g004

studies). Patterns of mitochondrial diversity in Algeria can be grouped into three main categories: i) Algeria is inhabited by the same lineage that is found in Morocco and Tunisia (*H. hippocrepis*, *C. ocellatus*, *D. pictus*), ii) a central lineage occupies parts of Algeria and sometimes adjacent territories between the previously identified eastern and western lineages that occur in Tunisia and Morocco respectively (*Pelophylax*, *Timon*, *Natrix*, *Trogonophis*, *Hyla*) or iii) Algeria is inhabited by several lineages with restricted range (*Podarcis*, *Acanthodactylus*, *Ptyodactylus*). Several of the lineages we have identified are endemic from Algeria based on current sampling (see Figs 1 and 2).

Sequence divergence between clades and diversity within clades are given in Table 2. In spite of the diversity of phylogeographic patterns uncovered by sampling in Algeria, the previously identified general pattern of an E-W divergence corresponding with the major break in mtDNA diversity (deepest node within North Africa) still holds for most taxa (see Figs 3–13).

In addition, a more complex situation emerges alongside the E-W general pattern in *Acanthodactylus* (an endemic lineage in the Atlas mountain of Morocco is of uncertain position and might be basal), *Podarcis* (a deeply divergent lineage from the Atlas Mountain of Morocco groups with a high support with the E clade, see REF) and *Ptyodactylus* (a basal lineage found in the pre-Sahara of Algeria is of uncertain position).

All studied taxa showed a high haplotypic diversity (Mean Hd = 0.094 ± 0.05, Table 3). The (E) and (W) clades don't show a significant difference in Hd (Kruskal-Wallis chi-squared = 0.434, df = 1, p = 0.51, see Table 3).
<figure><img src="image_9.png" /><figcaption>Fig 5. *Hemorrhois hippocrepis* Bayesian inference phylogenetic tree and geographic distribution of deep clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g005

## Location of E-W suture zones

Algerian populations belong to the E or W main lineages depending on the taxon (Fig 14). For most species the limits between the E and W lineages span a narrow longitudinal range but for *Timon*, specimens of the E lineages occupy the north of Algeria while populations of the W lineage are found along the Saharan Atlas in Algeria, resulting in a north-south phylogeographic divide within Algeria that contrast with the situation with the other taxa.

Two repeated patterns are identified for the location of the E-W phylogeographic breaks: the Moulouya suture zone around the Algeria-Morocco border and the Kabylian suture zone in the eastern part of Algeria (Fig 14). Six taxa out of 11 obey the “Moulouya pattern”, with suture zones located between 3.4 and 0.1°W around the Algerian-Moroccan boundaries; four species have suture zones in the Kabylia region between 3.6 and 5.3°E. *Hyla meridionalis* is the only species that deviates from these patterns as its disjunction area is located further east around 7.5°E. As can be seen from Fig 13, the limits in the Kabylian suture zone are more tightly spaced than around the Moulouya valley.
<figure><img src="image_10.png" /><figcaption>Fig 6. *Natrix maura* Bayesian inference phylogenetic tree and geographic distribution of deep clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g006

**Divergence time of the E-W clades.** Divergence times estimated for nine of the eleven taxa, as well as the 95% highest posterior densities (HPD), are reported in Table 3. Three main diversification periods can be identified (Table 3). The first period, represented by *P. oudrii*, *N. maura*, and *Timon.spp.*, with an old diversification dated to the Miocene, long before the Messinian crisis. The second diversification event corresponds with the Messinian salinity crisis (*A. erythrurus* and *C. ocellatus*). The last period corresponds with more recent divergence within the Maghreb after the Messinian crisis between 2 and 3.6 MYA (*D. pictus*, *H. meridionalis*, *P. saharicus*, *P. vaucheri*)

**Geographical distribution of internal clades and secondary phylogeographic suture zones**

As can be seen from Figs 1 and 2, the secondary suture zones between the internal clades delimited by ABGD do not follow a general trend as internal clade limits are spread over the whole of Algeria. The only common pattern is that most sutures are E-W breaks, reinforcing
<figure><img src="image_11.png" /><figcaption>**Fig 7.** *Podarcis voucheri* Bayesian inference phylogenetic tree and geographic distribution of deep clades. The red rectangle in the phylogenetic tree highlights the Moroccan lineage falling in the eastern clade.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g007

the idea that for Mediterranean or mesic species, longitudinal fragmentation of ancestral areas was the main force shaping divergence in mtDNA.

## Discussion

### How significant is the persistence of the E-W pattern after inclusion of Algerian samples?

The fact that we identified an E-W split as the most basal divergence within North Africa in most species (and that this E-W split structures most of the mitochondrial diversity in all of them) should not be interpreted as a support for the generality of the E-W divergence pattern in North Africa. We have chosen species where the E-W pattern had been previously recovered and we wanted to know if this pattern was disrupted by inclusion of Algerian samples. However, the persistence of this E-W divergence was not the only possible outcome of our study. In theory, new Algerian lineages could place at the basis of the mtDNA tree, leaving Moroccan and Tunisian samples in a monophyletic clade, if parts of Algeria harboured lineages that are older than the separation of Tunisian and Moroccan populations. The fact that we never recovered this pattern supports the idea that allopatry in E and W refugia is an important factor
<figure><img src="image_12.png" /><figcaption>**Fig 8.** *Ptyodactylus oudrii* Bayesian inference phylogenetic tree and geographic distribution of deep clades. The phylogenetic position of the specimen from Ghardaia (the blue dot on the map and blue highlighted branch in the tree) is not supported.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g008

driving partitioning of genetic mitochondrial diversity in the Maghreb herpetofauna. A proper examination of this hypothesis would require phylogeographic studies of a random sample of species with a large Maghreban distribution, however.

A rapid survey of published phylogeography for Maghreban amphibians and reptiles uncovers one other general patterns in addition to the main E-W divergence: the deepest divergence is located in Morocco, with the rest of the Maghreb colonized by one of the Moroccan lineages (*Mauremys leprosa*: [76]; *Daboia mauritanica*: [77]; *Bufotes boulengeri*: [78]. The other taxa either follow the general E-W divergence pattern: *Testudo graeca*: [79], *Salamandra algira* [80–81], or their phylogeography within North Africa was too poorly resolved to assign them to one of these situations: *Sclerophrys mauritanica*: [82]; *Psammodromus algirus*: [83].

In addition, more complex patterns have also emerged aside the main E-W divergence in North Africa. For example, in *Timon spp*, the disjunction corresponds with the western break area in the northern part of the range but isolated populations along the Saharan atlas to east-central Algeria belong to the western clade. These populations are geographically separated from the northern ones in Algeria carrying the eastern mtDNA haplotypes by the arid areas of the High Plateaux. Isolated populations genetically assigned to the Eastern clade are sometimes found within the range of western populations, for example in *Podarcis*. Last, basal clades
<figure><img src="image_13.png" /><figcaption>**Fig 9.** Bayesian inference phylogenetic tree and geographic distribution of deep clades. *Timon spp.* The reciprocal monophyly of E-W clades is not supported by any reconstruction method. However, the previous studies have shown that the Moroccan populations and the Tunisian ones are reciprocally monophyletic and highly supported (see references in the discussion).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g009

occur in the High Atlas Mountains (Morocco) for *A. erythrurus* or in the pre-Sahara for *P. oudrii*.

**Concordance of suture zones between E and W lineages in several taxa
suggest common historical processes**

The limits of the E and W clades are concentrated in two main areas: a suture zone around the boundaries of Algeria and Morocco, broadly corresponding with the previously suggested suture zone of the Moulouya valley and one newly identified disjunction zone in the Kabylia region (see [34,44]). Both areas are also suture zones between ABGD lineages for several species. This repeated pattern, supported by results of other studies, indicate that common historical processes shaped the current divergence around these two areas. Note that there is no apparent difference in TMRCA of E and W lineages between the two main disjunction zones (see Table 3).

The Kabylian break is located around the Soumam Valley in the centre-east of Algeria. It had not been identified before this study, even if *Podarcis vaucheri* [56] and *Cornu aspersum* [84] seemed to conform to this pattern. We show here that this area also constitutes a major divide for *Chalcides ocellatus*, *Hemorrhois hippocrepis* (Fig 1) and *Pelophylax saharicus* (Fig 2),
<figure><img src="image_14.png" /><figcaption>**Fig 10.** *Trogonophis wiegmanni*. Bayesian inference phylogenetic tree and geographic distribution of deep clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g010

hence for four out of 11 taxa. The Kabylian break is also the second main disjunction area for *Natrix maura* and *Timon pater* (Fig 1).

Six of our 11 taxa conform more or less strictly to the western break pattern (see Fig 1), while this area is also a secondary break for *Hyla meridonalis* and *Pelophylax saharicus* (Fig 2). This suture zone has been identified by several previous studies [34,44,70,85]. These studies suggested that this suture zone was a direct consequence of the Moulouya River valley and its arid corridor, acting as a barrier for Mediterranean species (but see below).

## Geographical barriers and glacial refugia

Our results suggest that the transition zone is wider than the Moulouya River valley itself (see Figs 1 and 2), suggesting that it cannot be the direct consequence of this valley. In addition, the W suture zone affects species with very diverse ecology, some of them occur in the Moulouya River valley itself (for ex. *Pelophylax*, *Natrix maura*, *Discoglossus*, *Acanthodactylus erythrurus*). It is thus obvious that the western suture zone cannot be explained by the Moulouya River itself acting as a barrier to dispersal. The Moulouya River is located in the Rifan corridor, a former marine corridor joining the Mediterranean east of the Rif Mountains with the Atlantic near present-day Rabat [86,87]. The width of this marine corridor was at least 60 km. The corridor opened between ~ 8.0 and 7.6 Ma, and the maximum sea level was reached in the late-Tortonian and decreased at the Tortonian-Messinian limit (7.2 MYA) but the connexion between the Mediterranean and the Atlantic continued during the Messinian. The final closure
<figure><img src="image_15.png" /><figcaption>**Fig 11. Discoglossus pictus.** Bayesian inference phylogenetic tree and geographic distribution of deep clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g011

of the Rif corridor is dated to ~6.1 to 5.6 MYA with a large marine basin persisting until the end of the upper Pliocene 5,6 MYA [86,88–90].

The Rifan corridor was undoubtedly a major geographic barrier that split the Maghreb for a long period (Miocene to mid Pliocene) and several E-W divergences that we have identified can be dated to this period (see Table 3). The fact that the western break zone coincides with the area of the Rifan Corridor and that the molecular dating of these splits corresponds broadly to this period suggests that the distribution of many species were fragmented on either side of this corridor. This major past geographical barrier could explain East-west divergence in the North and North-South divergence for other species in the southern part of the corridor (Middle Atlas).

The Kabylian break is located around the Soumam valley. Although, like the Moulouya, the Soumam is a relatively dry valley compared to adjacent areas, several species that are currently affected by the Kabylian break inhabit the valley itself. We thus suggest that the Kabylian break is not a direct consequence of the presence of the Soumam River either. Unfortunately, the palaeogeography of the area is poorly documented, but no major barrier comparable to the Rifan Corridor has been identified in Kabylia.

As mentioned above (introduction), the role of the glacial / interglacial cycles in generating biodiversity below or around the species level has been thoroughly documented in Europe and North America. Nothing comparable exists for Northern Africa but we show here with a
<figure><img src="image_16.png" /><figcaption>**Fig 12.** *Hyla meridionalis*. Bayesian inference phylogenetic tree and geographic distribution of deep clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g012

comparative phylogeography approach a good congruence between the Kabylia suture zones and glacial refugia identified by [20] on the basis of intra-specific phylogeographies of plant species.

The major paleo barrier of the Rifan Corridor could thus explain a fair proportion of ancient divergence events in Maghreban species but the origin of many cladogenesis events remain elusive. However, the conservation of these old divergent lineages requires a long term
<figure><img src="image_17.png" /><figcaption>**Fig 13.** *Pelophylax saharicus*. Bayesian inference phylogenetic tree and geographic distribution of deep clades. https://doi.org/10.1371/journal.pone.0201218.g013</figcaption></figure>

separation or a repeated expansion / contraction of the ranges of these lineages. Given the con-
gruence in Kabylia between the location of suture zones and the location of plant refugia, the
hypothesis of repeated contraction to glacial refugia seem the most likely scenario to explain
the current diversity at species and lineage scales.

**Evolutionary significance of mitochondrial diversity**

The lack of contemporary geographical barrier and the narrowness of the contact zones sug-
gests the possibility of competitive exclusion between lineages. Such reciprocal exclusion can
be neutral (drift in local demes) or selective (if the mtDNA lineages do not freely interbred and
exclude each other). We need a multilocus approach to examine these hypotheses. Firstly, we
need to verify that the mtDNA lineages that we have identified correspond with real evolution-
ary units, as profound mitochondrial divergence can sometime persist without nuclear diver-
gence (see [91]). We already know that some of the mitochondrial divergence we have studied
are mirrored in nuclear genes: for *Pelophylax*, the central and W lineages differ in allozymes
[92], (no sample from the E lineage were included) and for *Podarcis* E and W lineages are
highly divergent in allozymes [67]. Secondly, we need to assess whether these evolutionary
units persist in the face of gene flow (implying restriction to dispersal) or are associated with
speciation events. Indeed, several of the “species” that we have examined are probably made of
several species, implying that current taxonomy is inadequate.

A better understanding of the role of speciation events in generating the pattern of diversity
that we have exposed here would also have important consequences in conservation. If some
<table><caption>Table 3. Relevant node ages and haplotypic diversity.</caption><thead><tr><th>Taxon</th><th>Suture Zones</th><th>N</th><th>Gene</th><th>Lenght (pb)</th><th>S</th><th>H</th><th>Hd</th><th>Pi</th><th>MRCA (MY)</th><th></th><th>N</th><th>S</th><th>H</th><th>Hd</th><th>Pi</th></tr></thead><tbody><tr><td rowspan="2">Acanthodactylus erythrurus</td><td rowspan="2">M</td><td rowspan="2">92</td><td rowspan="2">ND4 + tRNA</td><td rowspan="2">755</td><td rowspan="2">285</td><td rowspan="2">70</td><td rowspan="2">0.99</td><td rowspan="2">0.096</td><td rowspan="2">4.5 [3-6.25]</td><td>Est</td><td>32</td><td>208</td><td>23</td><td>0.97</td><td>0.08</td></tr><tr><td>West</td><td>60</td><td>224</td><td>47</td><td>0.99</td><td>0.08</td></tr><tr><td rowspan="2">Chalcides ocellatus</td><td rowspan="2">K</td><td rowspan="2">51</td><td rowspan="2">CytB</td><td rowspan="2">425</td><td rowspan="2">113</td><td rowspan="2">44</td><td rowspan="2">0.99</td><td rowspan="2">0.06</td><td rowspan="2">4.57* [4.4-4.74]</td><td>Est</td><td>17</td><td>46</td><td>14</td><td>0.98</td><td>0.024</td></tr><tr><td>West</td><td>34</td><td>98</td><td>30</td><td>0.98</td><td>0.036</td></tr><tr><td rowspan="2">Hemorrhois hippocrepis</td><td rowspan="2">K</td><td rowspan="2">35</td><td rowspan="2">Cytb</td><td rowspan="2">300</td><td rowspan="2">15</td><td rowspan="2">11</td><td rowspan="2">0.83</td><td rowspan="2">0.009</td><td rowspan="2">-</td><td>Est</td><td>9</td><td>5</td><td>3</td><td>0.55</td><td>0.007</td></tr><tr><td>West</td><td>26</td><td>10</td><td>9</td><td>0.74</td><td>0.005</td></tr><tr><td rowspan="2">Natrix maura</td><td rowspan="2">M</td><td rowspan="2">48</td><td rowspan="2">ND4</td><td rowspan="2">604</td><td rowspan="2">49</td><td rowspan="2">17</td><td rowspan="2">0.87</td><td rowspan="2">0.025</td><td rowspan="2">9.2 [5.7-14]</td><td>Est</td><td>34</td><td>44</td><td>13</td><td>0.82</td><td>0.015</td></tr><tr><td>West</td><td>14</td><td>6</td><td>5</td><td>0.59</td><td>0.0017</td></tr><tr><td rowspan="2">Podarcis vaucheri</td><td rowspan="2">K</td><td rowspan="2">54</td><td rowspan="2">ND4 + tRNA</td><td rowspan="2">659</td><td rowspan="2">180</td><td rowspan="2">39</td><td rowspan="2">0.98</td><td rowspan="2">0.075</td><td rowspan="2">2.2 [1.5-3]</td><td>Est</td><td>25</td><td>127</td><td>19</td><td>0.97</td><td>0.06</td></tr><tr><td>West</td><td>29</td><td>69</td><td>20</td><td>0.95</td><td>0.018</td></tr><tr><td rowspan="2">Ptyodactylus oudrii</td><td rowspan="2">M</td><td rowspan="2">51</td><td rowspan="2">12S-RNA</td><td rowspan="2">350</td><td rowspan="2">92</td><td rowspan="2">28</td><td rowspan="2">0.96</td><td rowspan="2">0.093</td><td rowspan="2">10.6 [6.5-15.4]</td><td>Est</td><td>21</td><td>80</td><td>14</td><td>0.93</td><td>0.005</td></tr><tr><td>West</td><td>30</td><td>82</td><td>18</td><td>0.95</td><td>0.008</td></tr><tr><td rowspan="2">Timon spp</td><td rowspan="2"></td><td rowspan="2">69</td><td rowspan="2">16S-RNA</td><td rowspan="2">477</td><td rowspan="2">46</td><td rowspan="2">19</td><td rowspan="2">0.92</td><td rowspan="2">0.33</td><td rowspan="2">6.9 [3.8-11.4]</td><td>Est</td><td>32</td><td>27</td><td>8</td><td>0.8</td><td>0.024</td></tr><tr><td>West</td><td>37</td><td>28</td><td>11</td><td>0.88</td><td>0.019</td></tr><tr><td rowspan="2">Trogonophis wiegmanni</td><td rowspan="2">M</td><td rowspan="2">30</td><td rowspan="2">16S-RNA</td><td rowspan="2">474</td><td rowspan="2">48</td><td rowspan="2">20</td><td rowspan="2">0.97</td><td rowspan="2">0.034</td><td rowspan="2">-</td><td>Est</td><td>12</td><td>19</td><td>7</td><td>0.89</td><td>0.017</td></tr><tr><td>West</td><td>18</td><td>37</td><td>13</td><td>0.96</td><td>0.025</td></tr><tr><td rowspan="2">Discoglossus pictus</td><td rowspan="2">M</td><td rowspan="2">58</td><td rowspan="2">CytB</td><td rowspan="2">936</td><td rowspan="2">68</td><td rowspan="2">36</td><td rowspan="2">0.97</td><td rowspan="2">0.009</td><td rowspan="2">2 [1-3.3]</td><td>Est</td><td>52</td><td>52</td><td>31</td><td>0.96</td><td>0.006</td></tr><tr><td>West</td><td>6</td><td>7</td><td>5</td><td>0.93</td><td>0.002</td></tr><tr><td rowspan="2">Hyla meridionalis</td><td rowspan="2">E</td><td rowspan="2">50</td><td rowspan="2">CO1</td><td rowspan="2">775</td><td rowspan="2">97</td><td rowspan="2">34</td><td rowspan="2">0.96</td><td rowspan="2">0.031</td><td rowspan="2">3 [1.5-5.76]</td><td>Est</td><td>10</td><td>22</td><td>9</td><td>0.97</td><td>0.009</td></tr><tr><td>West</td><td>40</td><td>46</td><td>25</td><td>0.95</td><td>0.014</td></tr><tr><td rowspan="2">Pelophylax saharicus</td><td rowspan="2">K</td><td rowspan="2">63</td><td rowspan="2">CO1</td><td rowspan="2">654</td><td rowspan="2">114</td><td rowspan="2">33</td><td rowspan="2">0.93</td><td rowspan="2">0.066</td><td rowspan="2">3.65 [1.8-5.4]</td><td>Est</td><td>26</td><td>21</td><td>16</td><td>0.85</td><td>0.003</td></tr><tr><td>West</td><td>37</td><td>64</td><td>17</td><td>0.87</td><td>0.028</td></tr></tbody></table>

N: Number of specimens for each taxon, gene: genetic marker, L: sequence length, S: number of segregating sites, H: number of haplotypes in the complete data set of each taxon, Hd: haplotypic diversity calculated by DnaSp, Pi: TMRCA: age of most recent ancestor of E-W clades.

*: calibrations obtained from [36]

https://doi.org/10.1371/journal.pone.0201218.t003

of the lineages with restricted distribution that we identified in Algeria correspond to valid biological species, the level of endemism and conservation status of the Algerian herpetofauna would be significantly changed. Indeed, the western and eastern lineages of *Timon* are already widely recognized as two species (*T. tangitanus* in the west and *T. pater* in the east) even if this might require confirmation. A recent detailed investigation on *Tarentola mauritanica* in Morocco [93] identified several unrecognized species and *Ptyodactylus oudrii* has been suggested to constitute a species complex [39]. Many similar cases probably await discovery in Algeria, where the amount of detailed studies on small vertebrates is much lower and ABGD identified multiple lineages in several “species” that we investigated. While some of these lineages undoubtedly correspond with truly intraspecific diversity, a multilocus study is needed in their contact zones to figure out the level of reproductive isolation between these lineages.

## Impact of poorly resolved phylogenies

The mitochondrial phylogenies we used were sometimes poorly supported. While in several taxa, it was a consequence of us using short sequences and the original studies reported better supported phylogenies, in *Trogonophis*, the obtained phylogeny lacks strong nodal supports (posterior probability = 0.55), however a recent study resulted in the same overall topology with highly supported nodes [94]. In the case of *Timon spp*, the original study [42] reported a
<figure><img src="image_18.png" /><figcaption>**Fig 14.** A: Location of suture zones between deep E-W clades. B: Location of internal suture zones between ABGD groups. A: Coloured lines correspond to the disjunction between E and W clades. B: Coloured lines correspond to the limits between ABGD lineages. Red areas correspond to biogeographical refugia described by [36].</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201218.g014

well-supported phylogeny but we were unable to recover the same topology with our own analyses of the same dataset whatever the method used. These cases highlight the need for improved phylogenetic resolution through longer mitochondrial sequences and addition of multilocus data.

**Supporting information**

S1 Table. The list of used tissue samples, their geographic origin and GenBank accession numbers. Latitudes (LAT) and Longitudes (LON) are given in WGS84 geodetic system. The locus is specified only when there is concatenation of 2 loci for the phylogenetic
reconstruction, otherwise see the material and methods for the used marker.
(DOCX)

S2 Table. List of primers, their references and hybridisation temperatures.
(DOCX)

**Acknowledgments**

Menad BEDDEK was supported by a Naturalia-Environnement/ANRT CIFRE convention N° 2013/0145. Field work was mainly funded by Naturalia Environnement. We thank people who helped us in the field work in Algeria especially Abdelkrim Belhadj, Hibat Allah Loumassin, Benjamin Allegrini, We are grateful to the Algerian authorities from the general forests direction (Direction générale des forêts, DGF) and national parks (NPs) for the facilities, Saïd Abderrahmani (Belezma NP), Houari Djardini (Theniet El Had NP), Nadia Ramdane (Taza NP), Amina Fellous (National agency for Nature Conservation) & Mr. Allali and Mr. Ould Mohamed (DGF) for their help during field work. Sofiane Benharriga and all members of the marine ecological association “Barbarous” helped with field work in Oran. We thank. Special thanks for Mohamed Mebarki and the direction of “Le Galion” Hotel for accommodation, ad to Sana Beddek-Chelbabi and Mohamed Khiar for their precious help. We are grateful to the Moroccan “Haut Commissariat aux Eaux et Forêts et à la Lutte Contre la Désertification” for delivering our capture permit, and to Mr. Fekhaoui and Mr. Hassani from Institut Scientifique de Rabat for their help in organising field work in Morocco. The molecular data were generated at the genetic platform of our lab (Centre d'Ecologie Fonctionnelle et Evolutive): “SMGE” (Service des Marqueurs Génétiques en Ecologie). We thank all the SMGE team for their help.

**Author Contributions**

**Conceptualization:** Menad Beddek, Boualem Dellaoui, Fatiha Koudache, Olivier Peyre, Pierre-André Crochet.

**Data curation:** Menad Beddek, Philippe Geniez, Raouaa Fathalla, Patricia Sourouille, Véro-nique Arnal, Fatiha Koudache, Salah Telailia, Olivier Peyre, Pierre-André Crochet.

**Formal analysis:** Menad Beddek, Saliha Zenboudji-Beddek.

**Funding acquisition:** Olivier Peyre, Pierre-André Crochet.

**Investigation:** Menad Beddek, Pierre-André Crochet.

**Methodology:** Menad Beddek, Pierre-André Crochet.

**Project administration:** Olivier Peyre, Pierre-André Crochet.

**Resources:** Olivier Peyre, Pierre-André Crochet.

**Software:** Menad Beddek, Saliha Zenboudji-Beddek.

**Supervision:** Olivier Peyre, Pierre-André Crochet.

**Validation:** Pierre-André Crochet.

**Visualization:** Pierre-André Crochet.

**Writing – original draft:** Menad Beddek, Pierre-André Crochet.

**Writing – review & editing:** Menad Beddek, Pierre-André Crochet.
References

1. Lomolino M V, Sax DF, Brown JH. Foundations of biogeography: classic papers with commentaries. University of Chicago Press; 2004.

2. Lomolino M V, Riddle BR, Brown JH, Brown JH. Biogeography. Sinauer Associates Sunderland, MA; 2006.

3. Myers N, Mittermeier R a, Mittermeier CG, da Fonseca G a, Kent J. Biodiversity hotspots for conserva-
tion priorities. Nature. 2000; 403: 853–8. https://doi.org/10.1038/35002501 PMID: 10706275

4. Brooks TM, Mittermeier RA, da Fonseca GAB, Gerlach J, Hoffmann M, Lamoreux JF, et al. Global bio-
diversity conservation priorities. Science (80-). American Association for the Advancement of Science;
2006; 313: 58–61.

5. Medail F, Quezel P. Hot-spots analysis for conservation of plant biodiversity in the Mediterranean Basin. Ann Missouri Bot Gard. JSTOR; 1997; 112–127. https://doi.org/10.2307/239995

6. Reyjol Y, Hugueny B, Pont D, Bianco PG, Beier U, Caiola N, et al. Patterns in species richness and endemism of European freshwater fish. Glob Ecol Biogeogr. 2007; 16: 65–75. https://doi.org/10.1111/j. 1466-8238.2006.00264.x

7. Randi E. Phylogeography of South European mammals BT—Phylogeography of Southern European Refugia: Evolutionary perspectives on the origins and conservation of European biodiversity. In: Weiss S, Ferrand N, editors. Dordrecht: Springer Netherlands; 2007. pp. 101–126. https://doi.org/10.1007/1-4020-4904-8_3

8. Cox N, Chanson J, Stuart S. The status and distribution of reptiles and amphibians of the Mediterranean Basin. IUCN; 2006.

9. Cuttelod A, García N, Abdul Malak D, Temple HJ, Katariya V. The Mediterranean: a biodiversity hotspot under threat. Wildl a Chang World–an Anal 2008 IUCN Red List Threat Species. by: IUCN, Gland, Switzerland; 2009; 89.

10. Fonseca MM, Brito JC, Paulo OS, Carretero MA, Harris DJ. Molecular Phylogenetics and Evolution Systematic and phylogeographical assessment of the Acanthodactylus erythrurus group (Reptilia: Lacertidae) based on phylogenetic analyses of mitochondrial and nuclear DNA. Mol Phylogenet Evol. Elsevier Inc.; 2009; 51: 131–142. https://doi.org/10.1016/j.ympev.2008.11.021 PMID: 19070671

11. Véla E, Benhouhou S. [Assessment of a new hotspot for plant biodiversity in the Mediterranean basin (North Africa)]. C R Biol. 2007; 330: 589–605. https://doi.org/10.1016/j.crvi.2007.04.006 PMID: 17637440

12. Yahi N, Vela E, Benhouhou S, Belair G De, Gharzouli R. Identifying Important Plants Areas (Key Biodi- versity Areas for Plants) in northern Algeria. J Threat Taxa. 2012; 4: 2753–2765. http://dx.doi.org/10. 11609/JoTT.o2998.2753-65

13. Busack SD. Biogeographic Analysis of the Herpetofauna. Natl Geogr Res. 1986; 2: 17–3G.

14. James Harris D, Froufe E. Taxonomic inflation: species concept or historical geopolitical bias? Trends Ecol Evol (Personal Ed. 2005; 20: 6-7-9. https://doi.org/10.1016/j.tree.2004.11.004 PMID: 16701332

15. Ficetola GF, Bonardi A, Sindaco R, Padoa-Schioppa E. Estimating patterns of reptile biodiversity in remote regions. *J Biogeogr.* 2013; 40: 1202–1211. https://doi.org/10.1111/jbi.12060

16. Gómez A, Lunt DH. Refugia within Refugia: Patterns of Phylogeographic Concordance in the Iberian Peninsula BT—Phylogeography of Southern European Refugia: Evolutionary perspectives on the origins and conservation of European biodiversity. In: Weiss S, Ferrand N, editors. Dordrecht: Springer Netherlands; 2007. pp. 155–188. https://doi.org/10.1007/1-4020-4904-8_5

17. Schmitt T, Rákosy L, Abadjiev S, Müller P. Multiple differentiation centres of a non-Mediterranean but-
terfly species in south-eastern Europe. J Biogeogr. Wiley Online Library; 2007; 34: 939–950. https://doi.
org/10.1111/j.1365-2699.2006.01684.x

18. Nieto Feliner G. Patterns and processes in plant phylogeography in the Mediterranean Basin. A review. Perspect Plant Ecol Evol Syst. Elsevier GmbH.; 2014; 16: 265–278. https://doi.org/10.1016/j.ppees. 2014.07.002

19. García-alix A, García A, Minwer-barakat R, Mart E. Updating the Europe–Africa small mammal exchange during the late Messinian. J Biogeogr. 2016; 2017: 10–13. https://doi.org/10.1111/jbi.12732

20. Médail F, Diadema K. Glacial refugia influence plant diversity patterns in the Mediterranean Basin. J Biogeogr. 2009; 36: 1333–1345. https://doi.org/10.1111/j.1365-2699.2008.02051.x

21. Veríssimo J, Znari M, Stuckas H, Fritz U, Pereira P, Teixeira J, et al. Pleistocene diversification in Morocco and recent demographic expansion in the Mediterranean pond turtle *Mauremys leprosa*. *Biol J Linn Soc.* 2016; https://doi.org/10.1111/bij.12849
22. Ahmadzadeh F, Flecks M, Carretero MA, Ihlow F, Kapli P, Miraldo A, et al. Separate histories in both sides of the Mediterranean: phylogeny and niche evolution of ocellated lizards. 2016; https://doi.org/10.1111/jbi.12703

23. Cosson JF, Hutterer R, Libois R, Sarà M, Taberlet P, Vogel P. Phylogeographical footprints of the Strait of Gibraltar and Quaternary climatic fluctuations in the western Mediterranean: A case study with the greater white-toothed shrew, Crocidura russula (Mammalia: Soricidae). Mol Ecol. 2005; 14: 1151–1162. https://doi.org/10.1111/j.1365-294X.2005.02476.x PMID: 15773942

24. Habel JC, Lens L, Rödder D, Schmitt T. From Africa to Europe and back: refugia and range shifts cause high genetic differentiation in the Marbled White butterfly Melanargia galathea. BMC Evol Biol. 2011; 11: 215. https://doi.org/10.1186/1471-2148-11-215 PMID: 21777453

25. Husemann M, Schmitt T, Zachos FE, Ulrich W, Habel JC. Palaearctic biogeography revisited: evidence for the existence of a North African refugium for Western Palaearctic biota. Riddle B, editor. J Biogeogr. 2014; 41: 81–94. https://doi.org/10.1111/jbi.12180

26. Stewart JR, Lister AM, Barnes I, Dale L. Refugia revisited: individualistic responses of species in space and time. Proc R Soc B Biol Sci. 2010; 277: 661–671. https://doi.org/10.1098/rspb.2009.1272 PMID: 19864280

27. Allen H. Mediterranean ecogeography. Pearson Ed. 2001.

28. Avise JC. Phylogeography: the history and formation of species. Harvard university press; 2000.

29. Arbogast BS, Kenagy GJ. Comparative phylogeography as an integrative approach to historical biogeography. J Biogeogr. Wiley Online Library; 2001; 28: 819–825. https://doi.org/10.1046/j.1365-2699.2001.00594.x

30. Avise JC. Phylogeography: Retrospect and prospect. J Biogeogr. 2009; 36: 3–15. https://doi.org/10.1111/j.1365-2699.2008.02032.x

31. Hewitt GM. Speciation, hybrid zones and phylogeography—or seeing genes in space and time. Mol Ecol. 2001; 1: 537–549. https://doi.org/10.1046/j.1365-294x.2001.01202.x

32. Hewitt GM. Quaternary phylogeography: the roots of hybrid zones. Genetica. 2011; 139: 617–638. https://doi.org/10.1007/s10709-011-9547-3 PMID: 21234647

33. Taberlet P, Fumagalli L, Wust-Saucy A, Cosson J. Comparative phylogeography and postglacial colonization routes in Europe. Mol Ecol. Wiley Online Library; 1998; 7: 453–464. https://doi.org/10.1046/j.1365-294x.1998.00289.x PMID: 9628000

34. Barata M, Harris DJ, Castilho R. Comparative phylogeography of northwest African Natrixmaura (Serpentes: Colubridae) inferred from mtDNA sequences. African Zool. 2008; 43: 1–7. https://doi.org/10.3377/1562-7020(2008)43[1:CPONAN]2.0.CO;2

35. Guicking D, Joger U, Wink M. Molecular phylogeography of the viperine snake Natrixmaura (Serpentes: Colubridae): Evidence for strong intraspecific differentiation. Org Divers Evol. 2008; 8: 130–145. https://doi.org/10.1016/j.ode.2007.05.001

36. Kornilios P, Kyriazi P, Poulakakis N, Kumlutaş Y, Ilgaz C, Mylonas M, et al. Phylogeography of the ocelated skink Chalcides ocellatus (Squamata, Scincidae), with the use of mtDNA sequences: a hitch-hiker's guide to the Mediterranean. Mol Phylogenet Evol. 2010; 54: 445–56. https://doi.org/10.1016/j.ympev.2009.09.015 PMID: 19765663

37. Mendonça B, Harris DJ. Genetic variation within *Trogonophis wiegmanni* Kaup 1830. Belgian J Zool. 2007; 137: 239–242.

38. Perera A, Harris DJ, Harris J. Genetic variability within the Oudri's fan-footed gecko *Ptyodactylus oudrii* in North Africa assessed using mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol. 2010; 54: 634–639. https://doi.org/10.1016/j.ympev.2009.10.020 PMID: 19850142

39. Metallinou M, Červenka J, Crochet P-A, Kratochvil L, Wilms T, Geniez P, et al. Species on the rocks: Systematics and biogeography of the rock-dwelling Ptyodactylus geckos (Squamata: Phyllodactylidae) in North Africa and Arabia. Mol Phylogenet Evol. 2015; 85: 208–220. https://doi.org/10.1016/j.ympev.2015.02.010 PMID: 25724867

40. Perera A, Harris DJ. Genetic Variability in the Ocellated Lizard Timon tangitanus in Morocco. African Zool. Zoological Society of Southern Africa; 2010; 45: 321–329. https://doi.org/10.3377/004.045.0211

41. Paulo OS, Pinheiro J, Miraldo a., Bruford MW, Jordan WC, Nichols R a. The role of vicariance vs. dispersal in shaping genetic patterns in ocellated lizard species in the western Mediterranean. Mol Ecol. 2008; 17: 1535–1551. https://doi.org/10.1111/j.1365-294X.2008.03706.x PMID: 21928468

42. Recuero E, Iraola A, Rubio X, Machordom A, García-París M. Mitochondrial differentiation and biogeography of Hyla meridionalis (Anura: Hylidae): an unusual phylogeographical pattern. J Biogeogr. 2007; 34: 1207–1219. https://doi.org/10.1111/j.1365-2699.2007.01688.x
43. Stöck M, Dubey S, Klütsch C, Litvinchuk SN, Scheidt U, Perrin N. Mitochondrial and nuclear phylogeny of circum-Mediterranean tree frogs from the *Hyla arborea* group. Mol Phylogenet Evol. 2008; 49: 1019–1024. https://doi.org/10.1016/j.ympev.2008.08.029 PMID: 18814852

44. Nicolas V, Mataame A, Crochet P-A, Geniez P, Ohler A. Phylogeographic patterns in North African water frog Pelophylax saharicus (Anura: Ranidae). J Zool Syst Evol Res. 2015; n/a-n/a. https://doi.org/10.1111/jzs.12110

45. Sousa P, Harris DJ, Froufe E, van der Meijden a. Phylogeographic patterns of Buthus scorpions (Scorpiones: Buthidae) in the Maghreb and South-Western Europe based on CO1 mtDNA sequences. J Zool. 2012; 288: 66–75. https://doi.org/10.1111/j.1469-7998.2012.00925.x

46. Ben Faleh a. R, Granjon L, Tatard C, Othmen a. Ben, Said K, Cosson JF. Phylogeography of the Greater Egyptian Jerboa (Jaculus orientalis) (Rodentia: Dipodidae) in Mediterranean North Africa. J Zool. 2011; 286: n/a-n/a. https://doi.org/10.1111/j.1469-7998.2011.00868.x

47. Lalis A, Leblois R, Stoetzel E, Souttou K, Denys C, Nicolas V. Phylogeography and demographic history of Shaw 's Jird (Meriones shawii complex) in North Africa. 2016; 262–279. https://doi.org/10.1111/bij.12725

48. Pfenninger M, Véla E, Jesse R, Elejalde MA, Liberto F, Magnin F, et al. Temporal speciation pattern in the western Mediterranean genus Tudorella P. Fischer, 1885 (Gastropoda, Pomatiidae) supports the Tyrrhenian vicariance hypothesis. Mol Phylogenet Evol. 2010; 54: 427–436. https://doi.org/10.1016/j.ympev.2009.09.024 PMID: 19772926

49. García JT, Suárez F, Garza V, Calero-Riestra M, Hernández J, Pérez-Tris J. Genetic and phenotypic variation among geographically isolated populations of the globally threatened Dupont's lark Chersophorus duponti. Mol Phylogenet Evol. Elsevier; 2008; 46: 237–251. https://doi.org/10.1016/j.ympev.2007.06.022 PMID: 17719801

50. Guillaumet A, CROCHET P, Godelle B. Phenotypic variation in *Galerida* larks in Morocco: the role of history and natural selection. Mol Ecol. Wiley Online Library; 2005; 14: 3809–3821. https://doi.org/10.1111/j.1365-294X.2005.02696.x PMID: 16202098

51. Funk DJ, Omland KE. Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu Rev Ecol Evol Syst. Annual Reviews 4139 El Camino Way, PO Box 10139, Palo Alto, CA 94303–0139, USA; 2003; 34: 397–423. https://doi.org/10.1146/annurev.ecolsys.34.011802.132421

52. Denton RD, Kenyon LJ, Greenwald KR, Gibbs HL. Evolutionary basis of mitonuclear discordance between sister species of mole salamanders (Ambystoma sp.). Mol Ecol. Wiley Online Library; 2014; 23: 2811–2824. https://doi.org/10.1111/mec.12775 PMID: 24787951

53. Zink RM, Barrowclough GF. Mitochondrial DNA under siege in avian phylogeography. Mol Ecol. 2008; 17: 2107–2121. https://doi.org/10.1111/j.1365-294X.2008.03737.x PMID: 18397219

54. Kerr KCR, Stoeckle MY, Dove CJ, Weigt LA, Francis CM, Hebert PDN. Comprehensive DNA barcode coverage of North American birds. Mol Ecol Resour. Wiley Online Library; 2007; 7: 535–543.

55. Carranza S, Arnold EN, Pleguezuelos JM. Phylogeny, biogeography, and evolution of two Mediterranean snakes, Malpolon monspessulanus and Hemorrhois hippocrepis (Squamata, Colubridae), using mtDNA sequences. Mol Phylogenet Evol. 2006; 40: 532–546. https://doi.org/10.1016/j.ympev.2006.03.028 PMID: 16679033

56. Kaliontzopoulou A, Pinho C, Harris DJ, Carretero M a. When cryptic diversity blurs the picture: a cautionary tale from Iberian and North African *Podarcis* wall lizards. Biol J Linn Soc. 2011; 103: 779–800. https://doi.org/10.1111/j.1095-8312.2011.01703.x

57. Lima A, Pinho C, Larbes S, Carretero M, Brito JC, Harris DJ. Relationships of Podarcis wall lizards from Algeria based on mtDNA data. Amphibia-Reptilia. 2009; 30: 483–492. https://doi.org/10.1163/156853809789647103

58. Renoult JP, Geniez P, Beddek M, Crochet P-A. An isolated population of Podarcis vaucheri (Sauria: Lacertidae) in south-eastern Spain: Genetic data suggest human-mediated range expansion. Amphib Reptil. 2010; 31. https://doi.org/10.1163/156853810791069074

59. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. SMBE; 2016; msw054. https://doi.org/10.1093/molbev/msw054 PMID: 27004904

60. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. Oxford Univ Press; 2014; 30: 1312–1313. https://doi.org/10.1093/bioinformatics/btu033 PMID: 24451623

61. Guindon S, Lethiec F, Duroux P, Gascuel O. PHYML Online—a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. Oxford Univ Press; 2005; 33: W557–W559. https://doi.org/10.1093/nar/gki352 PMID: 15980534
62. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. Oxford University Press; 2012; 61: 539–542. https://doi.org/10.1093/sysbio/sys029 PMID: 22357727

63. Rambaud A, Drummond AJ. Tracer, version 1.4. 2007 [Internet]. Available from http://tree.bio.ed.ac.uk/software/tracer/

64. Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. Oxford University Press; 2012; 29: 1695–1701. https://doi.org/10.1093/molbev/mss020 PMID: 22319168

65. Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. BioMed Central; 2007; 7: 214. https://doi.org/10.1186/1471-2148-7-214 PMID: 17996036

66. Barata M, Carranza S, Harris DJ. Extreme genetic diversity in the lizard Atlantolacerta andreanskyi (Werner, 1929): A montane cryptic species complex. BMC Evol Biol. 2012; https://doi.org/10.1186/1471-2148-12-167 PMID: 22946997

67. Guarnizo CE, Amézquita A, Bermingham E. The relative roles of vicariance versus elevational gradients in the genetic differentiation of the high Andean tree frog, Dendropsophus labialis. Mol Phylogenet Evol. Elsevier Inc.; 2009; 50: 84–92. https://doi.org/10.1016/j.ympev.2008.10.005 PMID: 18957324

68. Pinho C, Harris DJ, Ferrand N. Contrasting patterns of population subdivision and historical demography in three western Mediterranean lizard species inferred from mitochondrial DNA variation. Mol Ecol. 2007; 16: 1191–205. https://doi.org/10.1111/j.1365-294X.2007.03230.x PMID: 17391406

69. Guo X, Dai X, Chen D, Papenfuss TJ, Ananjeva NB, Melnikov DA, et al. Phylogeny and divergence times of some racerunner lizards (Lacertidae: Eremias) inferred from mitochondrial 16S rRNA gene segments. Mol Phylogenet Evol. Elsevier Inc.; 2011; 61: 400–412. https://doi.org/10.1016/j.ympev.2011.06.022 PMID: 21767655

70. Li Z, Yu G, Rao D, Yang J. Phylogeography and Demographic History of Babina pleuraden (Anura, Ranidae) in Southwestern China. PLoS One. 2012; 7: e1–10. https://doi.org/10.1371/journal.pone.0034013 PMID: 22448286

71. Guicking D, Lawson R, Joger U, Wink M. Evolution and phylogeny of the genus Natrix (Serpentes: Colubridae). Biol J Linn Soc. 2006; 87: 127–143. https://doi.org/10.1111/j.1095-8312.2006.00561.x

72. Rambaut A, Drummond AJ. LogCombiner v1. 8.0. 2013. [Internet]. Available from: http://beast community/logcombiner

73. Rambaut, A., & Drummond, A. J. (2013). TreeAnnotator v1. 7.0. Available as part of the BEAST package at http://beast.bio.ed.ac.uk.

74. Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. Oxford Univ Press; 2003; 19: 2496–2497. PMID: 14668244

75. Puillandre N, Lambert A, Brouillet S, Achaz G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol Ecol. 2012; 21: 1864–1877. https://doi.org/10.1111/j.1365-294X.2011.05239.x PMID: 21883587

76. Fritz U, Barata M, Busack SD, Fritzsch G, Castilho R. Impact of mountain chains, sea straits and peripheral populations on genetic and taxonomic structure of a freshwater turtle, Mauremys leprosa (Reptilia, Testudines, Geoemydidae). Zool Scr. Wiley Online Library; 2006; 35: 97–108. https://doi.org/10.1111/j.1463-6409.2005.00218.x

77. Martínez-Freiría F, Crochet P-A, Fahd S, Geniez P, Brito JC, Velo-Antón G. Integrative phylogeographical and ecological analysis reveals multiple Pleistocene refugia for Mediterranean *Daboia* vipers in north-west Africa. Biol J Linn Soc. 2017; https://doi.org/10.1093/biolinnean/blx038

78. Nicolas V, Mataame A, Crochet P, Geniez P, Fahd S, Ohler A. Phylogeography and ecological niche modeling unravel the evolutionary history of the African green toad, Bufotes boulengeri boulengeri (Amphibia: Bufonidae), through the Quaternary. J Zool Syst Evol Res. Wiley Online Library; 2018; 56: 102–116.

79. Graciá EVA, Vargas-ramírez M, Delfino M, Anadón JD, Giménez A, Fahd S, et al. Expansion after expansion: dissecting the phylogeography of the widely distributed spur-thighed tortoise, Testudo graeca (Testudines: Testudinidae). Biol J Linn Soc. 2017; 1–14.

80. Merabet K, Sanchez E, Dahmana A, Bogaerts S, Donaire D, Steinfartz S, et al. Phylogeographic relationships and shallow mitochondrial divergence of Algerian populations of Salamandra algira. Amphibia-Reptilia. Brill; 2016; 37: 1–8.

81. Ben Hassine J, Gutiérrez-Rodríguez J, Escoriza D, Martínez-Solano I. Inferring the roles of vicariance, climate and topography in population differentiation in Salamandra algira (Caudata, Salamandridae). J Zool Syst Evol Res. 2016; n/a-n/a. https://doi.org/10.1111/jzs.12123
82. Harris D, Perera A. Phylogeography and genetic relationships of North African *Bufo mauritanicus* Schlegel, 1841 estimated from mitochondrial DNA sequences. Biologia (Bratisl). 2009; 64: 356–360. https://doi.org/10.2478/s11756-009-0042-6

83. Carranza S, Harris DJ, Arnold EN, Batista V, Gonzalez de la Vega JP. Phylogeography of the lacertid lizard, *Psammodromus algirus*, in Iberia and across the Strait of Gibraltar. J Biogeogr. 2006; 33: 1279–1288. https://doi.org/10.1111/j.1365-2699.2006.01491.x

84. Guiller A, Coutellec-Vreto MA, Madec L, Deunff J. Evolutionary history of the land snail *Helix aspersa* in the western mediterranean: Preliminary results inferred from mitochondrial DNA sequences. Mol Ecol. 2001; 10: 81–87. https://doi.org/10.1046/j.1365-294X.2001.01145.x PMID: 11251789

85. Nicolas V, Ndiaye A, Benazzou T, Souttou K, Delapre A, Denys C. Phylogeography of the North African dipodil (Rodentia: Muridae) based on cytochrome-b sequences. J Mammal. 2014; 95: 241–253. https://doi.org/10.1644/13-MAMM-A-241

86. Kouwenhoven T., Hilgen F., van der Zwaan G. Late Tortonian–early Messinian stepwise disruption of the Mediterranean–Atlantic connections: constraints from benthic foraminiferal and geochemical data. Palaeogeogr Palaeoclimatol Palaeoecol. 2003; 198: 303–319. https://doi.org/10.1016/S0031-0182(03)00472-3

87. Barhoun N, Bachiri TN. Événements biostratigraphiques et environnementaux enregistrés dans le corridor sud rifain (Maroc septentrional) au Miocène supérieur avant la crise de salinité messinienne. Geodiversitas. 2008; 30: 21–40.

88. Kouwenhoven TJ, van der Zwaan GJ. A reconstruction of late Miocene Mediterranean circulation patterns using benthic foraminifera. Palaeogeogr Palaeoclimatol Palaeoecol. 2006; 238: 373–385. http://dx.doi.org/10.1016/j.palaeo.2006.03.035

89. Bachiri TN, Barhoun N, et Suc JP. Les environnements continentaux du corridor rifain (Maroc) au Miocène supérieur d'après la palynologie. Geodiversitas, 2008, vol. 30, no 1, p. 41–58.

90. Achalhi M, Münch P, Cornée J-J, Azdimousa A, Melinte-Dobrinescu M, Quillévéré F, et al. The late Miocene Mediterranean-Atlantic connections through the North Rifian Corridor: New insights from the Bou-dinar and Arbaa Taourirt basins (northeastern Rif, Morocco). Palaeogeogr Palaeoclimatol Palaeoecol. 2016; 459: 131–152. https://doi.org/10.1016/j.palaeo.2016.06.040

91. Lindell J, Mendez-De La Cruz FR, Murphy RW. Deep biogeographical history and cytonuclear discor-dance in the black-tailed brush lizard (Urosaurus nigricaudus) of Baja California. Biol J Linn Soc. Oxford University Press; 2008; 94: 89–104. https://doi.org/10.1111/j.1095-8312.2008.00976.x

92. Arano B, Llorente GA, Montori A, Buckley D, Herrero P. Diversification in north-west African water frogs: molecular and morphological evidence. Herpetol J. British Herpetological Society; 1998; 8: 57–64.

93. Rato C, Harris DJ, Carranza S, Machado L, Perera A. The taxonomy of the Tarentola mauritanica species complex (Gekkota: Phyllodactylidae): Bayesian species delimitation supports six candidate species. Mol Phylogenet Evol. Elsevier Inc.; 2016; 94: 271–278. https://doi.org/10.1016/j.ympev.2015.09.008 PMID: 26391222

94. Salvi D, Perera A, Sampaio FL, Carranza S, Harris DJ. Underground cryptic speciation within the Maghreb: Multilocus phylogeography sheds light on the diversification of the checkerboard worm lizard *Trogonophis wiegmanni*. Mol Phylogenet Evol. Elsevier; 2018; 120: 118–128. https://doi.org/10.1016/j.ympev.2017.11.013 PMID: 29196204