Historical biogeography of the land snail Cornu aspersum: a new scenario inferred from haplotype distribution in the Western Mediterranean basin
© Guiller and Madec; licensee BioMed Central Ltd. 2010
Received: 10 July 2009
Accepted: 20 January 2010
Published: 20 January 2010
Despite its key location between the rest of the continent and Europe, research on the phylogeography of north African species remains very limited compared to European and North American taxa. The Mediterranean land mollusc Cornu aspersum (= Helix aspersa) is part of the few species widely sampled in north Africa for biogeographical analysis. It then provides an excellent biological model to understand phylogeographical patterns across the Mediterranean basin, and to evaluate hypotheses of population differentiation. We investigated here the phylogeography of this land snail to reassess the evolutionary scenario we previously considered for explaining its scattered distribution in the western Mediterranean, and to help to resolve the question of the direction of its range expansion (from north Africa to Europe or vice versa). By analysing simultaneously individuals from 73 sites sampled in its putative native range, the present work provides the first broad-scale screening of mitochondrial variation (cyt b and 16S rRNA genes) of C. aspersum.
Phylogeographical structure mirrored previous patterns inferred from anatomy and nuclear data, since all haplotypes could be ascribed to a B (West) or a C (East) lineage. Alternative migration models tested confirmed that C. aspersum most likely spread from north Africa to Europe. In addition to Kabylia in Algeria, which would have been successively a centre of dispersal and a zone of secondary contacts, we identified an area in Galicia where genetically distinct west and east type populations would have regained contact.
Vicariant and dispersal processes are reviewed and discussed in the light of signatures left in the geographical distribution of the genetic variation. In referring to Mediterranean taxa which show similar phylogeographical patterns, we proposed a parsimonious scenario to account for the "east-west" genetic splitting and the northward expansion of the western (B) clade which roughly involves (i) the dispersal of ancestral (eastern) types through Oligocene terranes in the Western Mediterranean (ii) the Tell Atlas orogenesis as gene flow barrier between future west and east populations, (iii) the impact of recurrent climatic fluctuations from mid-Pliocene to the last ice age, (iv) the loss of the eastern lineage during Pleistocene northwards expansion phases.
The aim of phylogeography is to explain the spatial distribution of genetic lineages within species and highlight the most influential historical factors explaining their distribution . Hypotheses commonly advanced in phylogeography are inherent to dispersal and vicariant events, inferences about which being usually based on estimates of genetic diversity, divergence time and demographic history . Inferences on demographic processes have especially become a central challenge in phylogeography and the recent development of analytical tools based on coalescent theory is very helpful for investigating the demographic history of populations. From information provided by phylogenetic trees and the signature in the pattern of DNA substitutions left in the case of historical change in population size, it is possible to identify the demographic process which has occurred and estimate related parameters [3, 4] The wealth of literature dealing with phylogeography and demography in many invertebrate and vertebrate species emphasizes the importance of combining both lineage and population components to reveal biogeographical histories. Moreover, the comparison of phylogeographic structures across co-distributed taxa allows to find concordant partitions and to infer common historical factors of divergence and colonization. Then, in North America (for review see ) as in Europe , various species have been studied to determine the climatic and geological effects on phylogeographic patterns and population structures. More precisely in Europe, the Mediterranean basin characterized by an exceptional level of biodiversity and a complex geological history , has been the subject of intensive studies carrying out on genetic diversity and phylogeography. Nevertheless, most of these studies provide insights into the evolutionary history of southern European species. Despite its key location between the rest of the continent and Europe, only a few cases concern the north African territory. In addition, investigations including Algerian data sets are relatively rare, and main historical process influencing the differentiation is generally found in Quaternary climatic oscillations, meaning that few studies imply events which occurred over times much longer than the last glacial period.
Two competing biogeographical scenarios based on opposite directions of dispersal were considered to explain the spatial pattern of variation in C. aspersum and especially to determine processes responsible for genetic isolation in the eastern part of north Africa (Fig 2) . Both were based on Pliocene geological events and Quaternary cold periods, but colonization routes taken by the species involved migrations either from north Africa to Europe or from western Europe to north Africa, with likely secondary contacts in the Kabylian area. The doubtful direction of species expansion results from discrepancies of spatial structure between allozyme and partial mitochondrial data, and from an inconstant pattern of gene diversity. The previous scenario we proposed [9, 14], implying a northwards progression of the species (model 1, Fig 2), fits well with published reports that place the origin of C. aspersum in north Africa . However, mitochondrial gene diversity estimated in north Africa does not corroborate the northward colonization hypothesis . The value of nucleotide diversity, twice as weak in the east (clade C) as in the west (clade B) north African mtDNA lineages, together with unresolved relationships between eastern populations, is rather consistent with southward dispersal either via Tyrrhenian/Aegean routes or the Strait of Gibraltar (model 2, Fig 2) . Whilst in the first scenario we ascribed the lower diversity in east to the occurrence of transitory and prolonged bottlenecks inherent in successive migrations that populations experienced from Kabylia to Tunisia, this may also simply indicate a more recent divergence within the eastern region.
These two hypotheses were evaluated here by studying the phylogeography and the demographic history of C. aspersum over all European and African populations sampled. Since parameters relevant to fluctuating population size on the evolutionary time scale were never taken into account in previous phylogeographical studies of the species, we hoped that their estimation would help to shed more light upon the scenario we previously considered to account for the scattered distribution of C. aspersum [9, 10]. The aim of this work is to assess the geographic patterns of mtDNA diversity in order (i) to improve our understanding of the species range expansion and infer the dispersal direction, (ii) to elucidate the role of the Kabylian area suspected until now to be a suture zone in Algeria. Without ruling out the fact that demographic processes and selective sweeps can lead to convergent genetic signatures when only one genetic marker is analysed, we hope that demographic estimates will allow to reinforce one of the hypotheses relating to the sense of dispersion in western Mediterranean, as well as that concerning the likely contact zone in Kabylia. Recent demographic expansion in the western lineage, combined with low nucleotide diversity, should favour a south-north dispersal, whilst a recent colonization from northern locations to the south should indicate demographic stability and high molecular diversity indices within the western lineage (Fig 2).
Genetic diversity estimates and results of demographic tests for cyt b and 16S genes in main geographical subdivisions of C. aspersum: B1, B2 for cyt b and Ba, Bb for 16S are subgroups of the B lineage, C1 and C2 are subgroups of the C lineage; Z includes intermediate zone sequences.
θπ ± sd
0.95 ± 0.01
9.34 ± 4.34
0.98 ± 0.01
8.52 ± 4.01
0.90 ± 0.05
2.23 ± 1.30
0.89 ± 0.03
3.13 ± 1.84
0.91 ± 0.08
6.29 ± 3.69
0.97 ± 0.02
7.18 ± 3.83
0.89 ± 0.06
12.42 ± 6.74
0.98 ± 0.01
13.8 ± 6.23
0.95 ± 0.01
15.27 ± 6.83
0.98 ± 0.01
13.20 ± 6.04
0.95 ± 0.03
4.18 ± 2.4
0.88 ± 0.03
3.70 ± 2.11
0.81 ± 0.12
15.64 ± 8.76
0.98 ± 0.01
21.56 ± 9.57
Phylogenetic and network relationships among haplotypes
The left-skewness test for the cyt b gene indicated that the observed tree length distribution was significantly more skewed than expected from random expectation (g1 = -0.193, p < 0.01). The g1-value obtained for 16S rRNA (g1 = -0.109, p < 0.01) was also indicative of a highly significant structure. The hypothesis of homogeneous rate within alignment of 16S and cyt b genes respectively was not rejected by the molecular clock tests. The likelihoods with no molecular clock enforced were not significantly better than those with enforced molecular clock (cyt b: 2ΔlogL = 15.8, d.f. = 60, p = 1; 16S, 2ΔlogL = 64.3, d.f. = 72, p = 0.73). For both genes, the RRT detected no significant heterogeneity between the defined lineages. Calculation of divergence time from phylogenies was then justified.
Pattern of migration
For both genes analysed separately, results of AIC values comparing four migration models support an asymmetric model of migration from north African to European populations of the B lineage (Fig 2). The model 1 (MAE) had an AIC value of 6.0 and 5.9 for cyt b and 16S respectively, whilst values of the MEA (model 2) were 730.3 for cyt b and 14 92.1 for 16S.
Estimates of TMCRA (and 95% Highest Posterior Density) for both cyt b and 16S genes overall C. aspersum sequences and of the B and C lineages.
divergence rate (%)
% per Myr
% per Myr
95% HPD lower
95% HPD upper
95% HPD lower
95% HPD upper
95% HPD lower
95% HPD upper
Our results based on 16S and cyt b mtDNA variation showed phylogenetic relationships among C. aspersum haplotypes consistent with inferences based on anatomy, nuclear (allozyme) and mitochondrial data previously analysed [9, 11, 12, 15]. The split into two haplogroups is supported since haplotypes of populations newly investigated compared to the previous mtDNA study  scattered either in the B or in the C lineage without really forming a new one. In a wider sense of the nomenclature inferred from these previous works, these B and C lineages were named "west" and "east" respectively, although the origin of some samples is not always in accordance with this spatial identification. Indeed, both lineages can include haplotypes originating from geographically distant populations. These discrepancies between genetic and geography are suggestive of a relatively intricate phylogeographical history of C. aspersum. As we previously showed, the genetic structuring of populations does not seem to be simply a function of dispersal . Beyond vicariance involving Tertiary tectonic events and Pleistocene glaciations, migration due to human actions needs to be mentioned to account for the current distribution of the species outside north Africa. However, the relative role of each process seems to vary greatly from one region to another. Before assessing the contribution of each of the different processes and refining the evolutionary scenario for the emergence of the two distinct lineages, we will recapitulate the main results by focusing particularly on contrasting zones of the distribution area that appear to be major historical entities.
The key areas detected
The tricky question of the dispersal scenario of the species can be solved by the genetic information provided by the northern part of the species range and preliminary outcomes of northward versus southward migration models tested. Within this wide area represented by western sublineages (i.e. B2 for cyt b, and Bb for 16S), four points can argue in favour of a northward movement: the better support of an asymmetric dispersal model towards a north Africa to Europe migration (model 1, Fig 2), the very low nucleotide diversity obtained, the star-like structure of the haplotypes network, the presence of Algerian individuals scattered in all subclades knowing that such splitting concerns exclusively Algerian samples. Although not specific to these sublineages, recent demographic expansions registered in B2 and Bb strengthen this dispersal direction that could have possibly started from Algiers area (represented by A4, A5, A6 populations). Expansion time estimates are purely speculative since no calibration point is available. However, values inferred from divergence rates fitting with previous allozyme estimates (see below), i.e. from 1 to 5% per Myr, would suggest that expansions highlighted by historical demographic reconstructions would have occurred twice. The first would have happened during one of the Pleistocene interglacial periods (i.e. from around 700 000 years ago), and the second at mid-Holocene after the maximum extent of the ice sheets during the last glaciation (between 37 000 and 7 400 years ago).
One area in Algeria seems to have played a central role in the dispersal of the species. Both genes showed that the Kabylies (represented by samples A10-A14 and A25) would have constituted a significant centre for dispersal and differentiation (model 1, Fig 2). However, phylogenies display an unstable position of Kabylian sequences. They occupy a basal position within the western lineage except on the ML 16S tree where they form a sister group of the remaining individuals of C. aspersum. The population of Djemila (A14) is an exception since the cyt b haplotype of this sample lies within an eastern subclade that also includes haplotypes of neighbouring populations (A15, A16, A17). Therefore, the Kabylies could possibly be the ancestral range of either C. aspersum or of the western lineage only. A difficulty in distinguishing between these two hypotheses stems from genetic exchanges that could have then occurred between western and eastern populations, leading to disruption in genetic structuring (model 2, Fig 2). The pattern of allele distribution at 13 polymorphic allozyme loci in intermediate Algerian populations ranging from Azzefoun (A12, Great Kaylia) to El Hedaick (A16, Skikda) reflects the possibility that peripheral populations of both East and West North-Africa would have regained contact after having been isolated. As previously suspected with the occurrence of eastern typical allozymes found in Djemila , introgressive hybridization between snails of the mixed population could also explain the presence of one C (east) typical cyt b haplotype in Djemila. The derived position of this haplotype diverging from A16 would moreover support the assumption of horizontal transfer of the mitochondrial genome through hybridization. One might suppose that both C (east) and B (west) haplotypes coexist in that population but quite obviously, more individuals from Djemila and surrounding sites need to be analysed to refine the gene diversity and reinforce the gene exchanges assumption between lineages in contact. Evidence in favour of this contact model is the west-east clinal decrease in a pair of genital organ lengths, i.e. the bursa tract diverticulum and the flagellum . Suspected hybridizing populations and especially those from Djemila (A14), El Ancer (A15) and El Hedaick (A16) have medium-sized organs compared with eastern and western ones. Therefore, whether the Kabylies is a dispersal centre or not, it is likely a zone where genetically distinct individuals from east versus west met and mated, resulting in at least some offspring of mixed ancestry. It is worth noting that other individuals sampled outside the Kabylies join this putative dispersal centre, especially snails from Morocco, Crete and Turkey.
The C clade equivalent to the eastern lineage we previously defined from anatomy, shell and molecular features covers here a larger geographical range. In addition to original samples from eastern north Africa, western Algeria (A1 and A2) and Italy [9–11, 15], samples more distant geographically join this clade. These are isolated specimens from Corsica, Croatia (Istria) and the western coast of Iberia (Cape Finisterre and Muxia in Galicia, coastal sites in Portugal). Although there is no clear evidence to explain this grouping, we can nevertheless point out that all populations of this lineage are outside the zone defined by the -1°C January isotherm during the last glacial maximum (LGM), i.e. about 18 000 BP [17, 18]. Contrary to populations of the B lineage that might have expanded northwards after glacial episodes from either southern differentiation centres or after introductions, populations of the C clade would have been limited in their post-last glacial maximum dispersal.
The biogeographical scenario reconsidered: new insights into haplotype distribution in north Africa and Europe
Cornu aspersum is a typically anthropochorous species which is nowadays widespread throughout the world in many zones with climates differing from the original Mediterranean one. Its presence is reported on the American continents, as well as in Australia and in Asia. Therefore, the first explanation for resemblances between populations located on either side of the Mediterranean could be passive transport due to human activities. Transfers might have started as nearly as the Neolithic revolution (around 8 500 BP; [9, 12]) and nowadays they continue occurring in giving rise in some cases to catastrophic destruction of habitats. Such man-made dispersal is however insufficient to explain the genetic pattern of variation obtained. Signatures left in the geographical distribution and the genetic diversity of extant populations indicated differences in spatial patterning between north Africa and Europe that could not be interpreted as a result of human-assisted dispersal alone. While a real pattern of regional differentiation exists for African samples , a lack of genetic structure due to a great admixture of populations is detected in northern parts of the species range . Hence, European individuals of the same population might originate from distinct native populations more or less differentiated genetically. Different waves of dispersion could explain marked genetic divergence estimates among the homogeneous groups of populations we defined using spatial contiguity constraint (see  for details). For the B lineage, northern and eastern Spanish areas like French ones should have been twice colonized by recurrent individuals coming from native populations of the same region. By contrast, snails from south and west Iberia that exhibit genetic similarities with north African individuals would have originated from more highly differentiated sources, probably Europe and north Africa. The Galician populations of Cape Finisterre (E9) and Muxia (E13) mirror such findings with an extremely frequently recorded level of mitochondrial diversity. The co-occurrence within a given population of both B and C specific haplotypes contributes to extreme divergences (for cytb: θπ (E9) = 13.8 and θπ (E13) = 6.9; for 16S: θπ (E13) = 11.8 and θπ (E9) = 0.5 with only eastern types present in E9) that consequently leads to discord between gene and population trees. This intermixing of distinct haplotypes is questionable since it occurs only in one part of the species range which moreover is peripheral and not expected to be the differentiation centre of C. aspersum. The overall differentiation between these sympatric divergent haplotypes involves two likely sources of variation: (i) divergence which already existed within the common ancestral populations; (ii) divergent haplotypes originating from previously isolated lineages . These components of variation imply underlying debatable mechanisms such as incomplete lineage sorting of ancestral polymorphisms due to recently divergent populations or recent migration events between allopatric populations [20, 21]. Incomplete lineage sorting from an ancestral polymorphic gene pool could possibly be the cause of the mitochondrial polyphyly detected in E9 and E13. The high gene diversity observed overall C. aspersum sequences compared to the number of individuals analysed prevents reliable differentiation between shared (ancestral) and specific (derived) haplotypes. However, a distinction could be proposed using their position on networks. Then, eastern haplotypes for both genes would be specific to E9 and E13 whereas western mitotypes would be either specific to each population or shared by other samples more or less distant geographically (Iberian and French samples). For East types only, the absence of ancestral mitotype could correspond to the "allotypy" pattern of the intermediate polyphyly progression described by Omland et al.  and documented by other studies [23–25]. Intermixing of Iberian and East North-African haplotypes should then indicate that they would not be completely sorted. In that case, nuclear genes having about four times the effective population size of mitochondria, would however fail to locate populations which are reciprocally monophyletic . Although E9 and E13 are clearly distinct genetically from most other Iberian samples belonging to the same B clade, nuclear variation is effectively consistent with monophyly of both populations with no intermixing of typical alleles of B (west) versus C (east) lineages. Allozyme data indicated few common alleles shared only with western north African populations, especially M1, M2, A1 and A2 , and assignment tests performed on microsatellite genotypes clustered almost all individuals of both samples (92 of 95 individuals analysed) within a single cluster rather than distinct ones . Other arguments such (i) the occurrence of close western types of both Iberian populations in very remote geographical regions, (ii) the great genetic divergence estimated between the aforementioned intermixed B and C haplotypes, (iii) the geographical location of populations having non-monophyletic haplotypes, tend to argue against this hypothesis of incomplete sorting in promoting secondary contacts between previously isolated lineages. The pattern of morphometric distinctiveness between western and eastern typical populations also supports such hypothesis of migration between divergent populations. Only explanations based on secondary contacts between snails of the B and C lineages could clarify the morphometric heterogeneity in space and time registered in Galician populations, in which adult snails exhibit B or C shell and genital features according to sampling dates (1991 and 2004 for E9, 2004 for E13) and location. Comparison with the morphometric profile of snails originating from populations of the suspected contact zone in Kabylia may provide insights into the process responsible for the variation observed in Galicia. As mentioned above, snails from Djemila (A14), which can be considered as the most representative of the intermediate zone , have genitalia of an intermediate size, which is about average for the species [13, 28]. Such an intermediate pattern of covariation between shell and genital organs should imply mechanism leading to an eventual precopulatory isolation between snails from the eastern and western lineages. More than a simple coexistence of two distinct morphotypes within a population, the intermediate size of diverticulum of snails from Cape Finisterre could reflect the evolution of premating reproductive isolation among lineages recently in contact...The contradictory genetic signature of mitochondrial versus nuclear loci may support this scenario in which horizontal transfer of distinct haplotypes would occur through introgressive hybridization between both lineages where reproductive barriers are incomplete. The mitotypes characterizing southern and western populations of the peninsula strengthen the putative contact zone inferred in that area since Portuguese samples (P2, P3, P7) showed C types only, whereas northern Iberian populations (E10, E14-17) have B types. An extensive sampling of this western part of Iberia is however imperative to corroborate this conclusion.
Mediterranean taxa as examples of phylogenetic splits in north Africa (NA) (MSC: Messinian Salinity Crisis; ?: no available data).
Divergence Time (Myr)
the cork oak Quercus suber
2 sub-lineages in Morocco vs 1 in Tunisian
the pine Pinus pinaster
Morocco vs Tunisian lineages
Geographical barrier: Moulouya River basin
the frogs Discoglossus pictus-auritus/scovazzi
2 isolated lineages in contact in NA: D. pictus-auritus in Algeria-Tunisa vs D. scovazzi in Morocco
the scorpion Buthus occitanus
Atlas (Morocco) vs Tell-Atlas (Tunisia) lineages
the frog Hyla meridionalis
Morocco vs Tunisian lineages
Geographical barrier: Atlas Mountains in Morocco
13 - 11
the fresh turtle Mauremys leprosa
2 clades north vs south mountains
Geographical barrier: Atlas and Tell systems
13 - 11
the lizards Lacerta pater/L. tangitana
2 species west vs east mountains
Mid-Miocene split -end of MSC
Marine transgression (Edough peninsula)
the newts Pleurodeles waltl P.nebulosus P. poireti
Differentiation P.waltl/ancestor P.nebulosus-P. poireti
Differentiation P.waltl/ancestor P.nebulosus-P. poireti
Geographical barrier: Atlas uplift
the Agamid lizard Agama impalearis
2 lineages in Morocco: NW vs SE Atlas
the snake Malpolon monspessulanus
M.m. monspessulanus from Morocco/Algeria vs M.m insignitus from Tunisia lineages
Geographical barrier: Moulouya River basin, ?
the riverine snake Natrix maura
3 lineages: Morocco, east Morocco + Algeria, Tunisia+ east Algeria
the riverine snake Natrix maura
2 lineages: Tunisia vs Morocco
the salamander Salamandra algira
2 lineages: west Morocco vs east Morocco + Algeria
Marine transgression (Edough peninsula)
the salamander Salamandra algira
Morocco vs Edough peninsula lineages
Marine transgression (Edough peninsula)
the land snail Turodella sulcata
Turodella sulcata ss vs Turodella sp (Enough peninsula)
the lizard Acanthodactylus erythrurus belli
Morocco vs Algeria+Tunisia (A. blanci ) lineages #
in situ differentiation (with Moulouya River basin as possible geographical barrier) post to colonization from Europe
the wall lizard Podarcis hispanica s.l.
3 lineages: northen Morocco vs southern Morocco vs Tunisia
the shrew Crocidura russula
Moroccan vs Tunisian lineages
the frog Rana saharica
2 clades: R. s. saharica from Algeria. and R. s. riodeoroi from Morocco
the lizard Psammodromus algirus
Minor divergence between Morocco vs Tunisia clades
in situ differentiation
the Crested lark Galerida cristata
3 lineages: cristata (west NA) vs randonii (central NA) vs senegallensis (east NA)
Quaternary glacial refuges
the wall lizard Podarcis vaucheri
8 subsets in Morocco
Quaternary glacial refuges
the annual grass Hordeum marinum
Iberia vs Central Mediterranean lineages
the gecko Tarentola mauritanica
4 lineages: 2 in Morocco. 1 in Algeria. 1 in Tunisa
the winter pine processionary moth Thaumetopoea pityocampa/T. wilkinsoni
Morocco-west Algeria vs east Algeria-Tunisia lineages
Kerdelhué, comm. pers.
Finally, Pliocene and Pleistocene climatic fluctuations should have influenced the distribution and the evolution of C. aspersum in north Africa. However, signatures left provide insufficient information to infer precisely the modes and timing of the genetic B-C split. The influence of Quaternary climatic changes is more evident in Iberia, where the occurrence of both B and C lineages can be explained by the persistence of several separate refugia throughout the Pleistocene Ice Ages. Recent parallel phylogeographical patterns support the idea that the Iberian Peninsula comprised a "refugia-within-refugia" rather than a single survival area . Amongst the seven putative terrestrial refugia suggested by these authors (Fig 1), the Betic range and Picos de Europa coincide quite strikingly with the tectonic hypothesis proposed and the spatial genetic pattern inferred from allozymes and/or the haplotype distribution in C. aspersum. It seems reasonable to assume that each lineage previously differentiated evolved in isolated refugia during glacial periods, such as near the Betic range for the C clade and in the Picos de Europa area for the B lineage. Range expansion during favourable climatic conditions, that is, from the Günz-Mindel stage for B and the Mindel-Riss stage for C, would lead to the contact zone observed in Galicia. Moreover, the occurrence of C haplotypes in the Betic area would agree with the geological scenario based on microplate drift from Oligocene, since the Betic with the Rif Cordillera formed a continuous orogenic belt together with the Kabylies blocks, the Balearic Islands, Calabria, Corsica and Sardinia . However, haplotypes identified in Sardinia and Balearic Islands do not support this scenario but sampling is extremely poor in these areas as in Corsica. For the Tyrrhenian region, populations collected very recently in Sicily but not yet analysed should indicate close relationships with C haplotypes, especially Tunisian ones. Beyond geological assumption of a land connection between Europe and Africa in the continuity of the Maghrebides at the level of the present Sicilian canal until late the Pliocene (2 Myr; ), allozyme and morphometric features indicate similarities between southernmost Italian samples analysed and eastern north African populations [9, 14]. An alternative assumption of secondary contacts versus a differentiation centre in the southern Italian Peninsula could equally lead to resemblance between those peri-Tyrrhenian populations. Although we are very short of sequences from these areas, the derived position from eastern African haplotypes of the unique Italian 16S type would rather favour the hypothesis that Italy would have become a contact zone between African and European populations. Morover, as mentioned above, the most likely north African origin of Mediterranean species supports this scenario.
While keeping in mind the role of anthropogenic pathways in the spread of the species, time estimates and genetic discontinuities lead to suggest that Pliocene and Pleistocene climatic changes would have played a significant role in shaping variation in north Africa. However, the distribution of haplotypes compared to the present-day position of blocks drifted in the western Mediterranean basin during the Tertiary could not rule out the hypothesis referring to much an older vicariance process to account for the whole biogeographical pattern of C. aspersum. Contrary to the tectonic scenario we previously proposed , the circum-Mediterranean range of typical eastern types would favour dispersal through plate migration of an eastern ancestor rather than a western one. Subsequent orogenic events such as the Tell Atlas building that may have impeded gene exchange between populations on the two sides of the mountain chain, combined with climatic fluctuations from mid-Pliocene to the last ice age, could thus explain the east-west subdivision in north Africa. The time to the most recent common ancestor estimated for C. aspersum, as for each of the B (west) and C (east) lineages, supports this parsimonious scenario. The intermediate position of Kabylia, as regards both geography and genetics, would explain unstable phylogenetic relationships of haplotypes of this area. This main region, recently recognised with Numidia and Kroumiria as a hotspot of biodiversity in the Mediterranean basin , would have been first a centre of dispersal and then of secondary contact between isolated lineages (models 1 and 2, Fig 2). The retreat of populations into the Edough Peninsula, that seems to have served as refugia for many species during Pliocene marine transgression, could have accentuated the genetic cleavage between the B and C entities. Subsequent population expansion towards the east when floods ended could partially explain the low nuclear (allozyme) and mitochondrial variability characterizing especially Tunisian samples [9, 15]. Transitory and prolonged bottlenecks inherent to successive migrations that likely continued until the LGM would moreover be compatible with demographic expansion observed in that easternmost north African region. Also, northward expansions with climatic improvements since the Mindel-Riss interglacial stage would indicate a loss of the eastern lineage during expansion phases. Whatever the range of C. aspersum, the typical eastern populations seem to have been limited in their dispersal after the LGM.
To resume, the different steps of the scenario inferred involve a northward colonization of the western lineage that would have diverged form eastern ancestral type. Most data support the migration model 1, which however needs to be refined in considering the dual role of the Kabylia (Fig 2). In relation to vicariant and dispersal events that could have played a role in the colonization process of both lineages, discrepancies in the success of northward expansion could also reflect variation in adaptive potential of populations among B and C lineages. The extinction probability of populations in changing environments being higher in small than in large populations, one might imagine that typical eastern populations experienced severe and prolonged bottlenecks. As well as demographic studies, efforts should be made to research other possible causes of variation in the colonization process between both lineages. Moreover, this biogeographic scenario should be evaluated (i) by searching for concordant phylogeographic patterns of codistributed species to indicate the influence of common historical factors, (ii) by analysing populations with more extensive sampling, especially in the two putative contact zones identified (Kabylia and Galicia) as in southern Italian Peninsula, and (iii) in incorporating calibration points different from the dubious shell fossils recorded as C. aspersum.
Most individuals analysed here originated from populations previously sampled for anatomical and biochemical surveys [9, 12, 14, 16, 71]; other populations were sampled specifically for the current analysis. Overall, 73 localities were considered, 66 and 49 for 16S rRNA and cyt b genes respectively (Fig 1). Individuals of the subspecies maxima coming from a snail farm in Brittany (France) were also analysed and used as outgroups in phylogenetic and network-based analyses. Sample site characteristics and haplotype codes are given as additional file 1.
DNA extraction, amplification and sequencing
The protocol followed here has been described in . Briefly, total genomic DNA was obtained from foot and mantle muscles of either fresh or frozen material stored since previous investigations. Protein extracts of liver or muscle were used when foot and mantle were not available. We used either the CTAB, or the chelex extraction protocols for both frozen and fresh material . We amplified fragments of approximately 380 and 560 bp for the 16S rRNA and cyt b mitochondrial genes respectively. The 16S fragment was amplified using primers 16S-1 (5'-TGACTGTGCAAAGGTAGC-3') and 16S2 (5'-CTGGCTTACGCCGGTCTG-3') . The cyt b region was amplified using Cytb1 (5'-TTATTGAGGCGCTACGGTTAT-3') and Cytb2 (5'-GCAAGCGAAATATAAGGTTCT-3') primers. Amplification of template DNA was carried out in a volume of 25 μl consisting of 10 mM Tris-HCl (ph 9.0), 50 mM KCl, 1.5 mM MgCl2, 200 μm of each dNTP, stabilizers including BSA, 1.5 U of Taq DNA polymerase (Ready to Go®, Amersham Pharmacia Biotech.), 0.20 μm of each primer and 50 ng of DNA. The PCR conditions were for 16S rRNA, an initial denaturation step of 94°C (5 min), followed by 35 cycles of 94°C (30 s), 50°C (30 s), 72°C (40 s) and a final extension phase at 72°C for 3 min; for cyt b, an initial denaturation step of 94°C (5 min), followed by 35 cycles of 94°C (40 s), 50°C (40 s), 72°C (1 min) and a final extension phase at 72°C for 5 min. Double-strand sequences were obtained by direct manual sequencing of the PCR products using the Sequenase v.2 PCR product sequencing kit (Amersham), or using an automated sequencer (PE Applied Biosystems 310 Genetic Analyser, UMR 6553; plate-forme de séquençage et génotypage OUEST-genopole®).
Mitochondrial sequences were aligned using CLUSTALW, v1.7 . For the 16S RNA region, they were manually adjusted, based on the secondary structure of the large ribosomal subunit of C. nemoralis . The sequences were submitted to GenBank (Additional file 1). Analyses of sequence polymorphism were carried out using DNASP v4.10.9  and ARLEQUIN v3.1 .
The degree of phylogenetic signal was first evaluated in both gene data sets by performing the left-skewness (g1) test using 10 000 randomly generated trees under the maximum likelihood criterion, as implemented in PAUP* v4.0b10 . For inferring phylogenetic relationships among individuals, we used maximum likelihood (ML) and bayesian-based inference (BI) methods. The best fit model of nucleotide substitutions was selected prior to ML and BI analyses using the Akaike Information Criterion. The software MrAIC v1.4.2  was used to evaluate the fit of the data to 24 different models of nucleotide substitutions. For both genes, the resulting best fit model was HKY, with two different rates for transitions and transversions, unequal base frequencies, a parameter for invariable sites (I) and a gamma distribution parameter that describes rate variation across variable sites (Γ). This model of nucleotide substitution was incorporated in PHYML V2.4.4 and in MRBAYES v3.1.1-p1  for ML and BI analyses respectively. For ML analysis, the robustness of inferences was assessed by bootstrap resampling using 1000 repetitions. For Bayesian analysis, the posterior probabilities of trees and parameters were approximated with Markov Chain Monte Carlo (MCMC) and Metropolis coupling. For each gene, we ran two independent MCMC analyses with four chains each and a temperature set to 0.2. Each chain was run for 2 000 000 cycles with trees sampled every 100 generations. Posterior probabilities were obtained from the 50% majority rules consensus of trees sampled after discarding the trees saved before chains reached apparent stationarity (i.e. a 'burn-in period" of 50 000 generations for 16S rRNA, 52 000 for cyt b). Due to low divergence among individuals and possible persistence of ancestral nodes along with descendants, network-based approaches are better suited for intraspecific phylogeny than tree-based algorithms . Consequently for both mtDNA regions, the median joining algorithm implemented in Network v188.8.131.52 software  was also used with default settings for constructing networks (weight = 10 and ε = 0). Due to the high gene diversity in both mtDNA regions, each data set was first reduced using the star contraction option which identifies and contracts any starlike phylogenetic cluster into one ancestral type. Median networks that contained all possible equally short trees were simplified by running the maximum parsimony (MP) calculation option to eliminate superfluous nodes and links.
Before evaluating divergence-time estimates between lineages, we checked whether sequences were evolving at a uniform rate along all branches in a phylogeny. We first performed a molecular-clock likelihood-ratio test (LRT)  by comparing the likelihood scores (lnL1 and lnL2) of the same tree constructed under alternative molecular clock assumptions (relaxed versus enforced molecular clock). The statistic Δ = 2(lnL1-lnL2) can be compared with a χ2 distribution with (n-2) degrees of freedom (where n is the number of haplotypes). Secondly we conducted a relative rate test (RRT; ) as implemented in RRTree  for assessing variation in substitution rates between lineages relative to the outgroup C. a. maxima. For the RRT, p values were corrected using the Bonferroni method to account for multiple pairwise comparisons when more than two lineages were identified.
Pattern of migration
Alternate hypotheses of migration were tested by estimating historic migration rates between north Africa and European populations. We compared models employing asymmetric versus symmetric migration rates between European and north African clusters of populations for both genes separately. Because of the small number of European populations in the C (East) haplogroup, estimates were however restricted to the B (West) haplogroup with the Z sequences removed. We used the program MIGRATE-n v. 3.0  to obtain maximum likelihood estimates (MLEs) of theta (θ = Nμ) and Nm (effective number of migrants per generation). We obtained MLEs for four models of migration between both regions (bidirectional rates with asymmetric MNA-E ≠ ME-NA or symmetric MNA-E = ME-NA rates; unidirectional rates with MNA-E = 0 or ME-NA = 0). We used likelihood ratio tests and the Akaike Information Criterion to choose the model best supported by the data (smaller values of AIC indicating better fit). The program was run with 15 short chains with an increment of 50 and 2000 trees recorded, followed by 3 long chains of 25 000 000 generations each, from which 50 000 trees were sampled after a burnin-period of 5000. We also optimized all models in specifying Ti/Tv ratio (2.56 for cyt b and 1.51 for 16S) and gamma distribution shape parameter (α = 0.75 for both genes) from empirical estimates in PAUP* 4.0b .
The demographic history of populations was inferred using tests of population growth. Fu's Fs  and Ramos-Onsins & Rozas' R2  are among the most powerful statistics to detect demographic expansions; they were estimated using DNASP v4.10.9  and ARLEQUIN v3.1 , and their significances were assessed using 1000 coalescent simulated resamplings. We also performed mismatch distribution analyses to evaluate possible historical events of population growth and decline . Indeed, populations that have experienced a rapid expansion in the recent past show unimodal distributions, while those at demographic equilibrium have multimodal distributions. Mismatch distributions were computed for each haplogroup and compared to the expected distributions obtained under a model of sudden expansion.
An approach based on coalescent theory was also used to estimate the demographic history of lineages within C. aspersum. Under a coalescent model, it is possible to infer population parameters from genetic sequence data such as mutation rates and divergence time estimates which themselves provide information about effective population sizes through time. We chose the Bayesian Skyline Plot (BSP, ) model for inferring past population dynamics. The BSP has the advantage over commonly-used demographic models and skyline plot methods  of being independent of a prespecified parametric model of demographic history, and its parameters are directly estimated from sequences, thus limiting the error inherent in phylogenetic reconstruction. The BSP model used MCMC sampling procedures to estimate credibility intervals of the demographic parameters. Parameters were estimated with the software BEAST v1.4.2 . After an optimization step during which parameters of the prior function were changed at each run to reach optimum performance and achieve a reasonable effective sampling size (ESS, number of independent samples of the posterior distribution that the trace is equivalent to) of parameters of interest, we carried out two independent runs of 20 million generations each. Samples of trees and parameters were recorded every 1000 steps. Visualization and diagnosis analysis of each independent BEAST run were done using TRACER v1 . This program leads to redefine the burn-in period (set by default to 10% of the MCMC chain length) and check the convergence of the independent chains. The analysis was repeated after splitting the data into a B versus C haplogroup, motivated by earlier results. To provide estimates of time to the most common ancestor (TMCRA) of C. aspersum populations, along with the TMCRAs of the B and C lineages (see Results), divergence dates were computed using BEAST. Estimates were based on published mutation rates for land snails. Evolutionary rates at mtDNA genes have been estimated in a few terrestrial gastropods though not in C. aspersum. Evidence from Pliocene shell fossils (5.3-1.8 Myr) from several sites around the Mediterranean (Oran, Sicile, Nice; ) indicates that C. aspersum was already quite widespread in the late Tertiary. However such fossil records, even if shell specimens were assigned to C. aspersum species with sufficient confidence, do not allow to go down at the subspecies level (aspersa or maxima) even less down to mtDNA lineage within aspersa. Information from the fossil record, inadequate for calibrating a molecular clock, can only indicate the minimum age of C. aspersum. Estimates in gastropods, intraspecific published estimates vary greatly for 16S rRNA [94–96], as for protein-encoding mitochondrial genes [97–99]. Without more information, we used and tested a range of calibrated rates of change from 0.03 to 5% for cyt b excluding the third codon position, and from 0.5 to 10% for 16S rRNA, to calculate a rough estimate of the timing of divergence events. In all cases we used an HKY model of nucleotide substitution. Results were compared with the estimated times based on allozyme evolution (0.08-0.10 DNei/Myr; ) which we published previously .
We warmly thank P. Jarne and Y. Surget-Groba for their advice and valuable criticisms on the manuscript. The comments of S. Chiba and two anonymous reviewers improved greatly our manuscript. We are also grateful for information on geodynamics and climatic of the Mediterranean basin provided by D. Frizon Lamotte. Thanks also to A. Scaife for the linguistic revision of the manuscript.
- Avise JC, Arnold J, Ball RM, Bermingham E, Lamb T, Neigel JE, Reeb CA, Saunders NC: Intraspecific phylogeography: the mitochondrial DNA bridge between population genetics and systematics. Ann Rev Ecol Syst. 1987, 18: 489-522.View ArticleGoogle Scholar
- Dawson MN, Waples RS, Bernardi G: Phylogeography. Ecology of California Marine Fishes. Edited by: Allen LG, Horn MH, Pondella DJ II. 2004, 26-54.Google Scholar
- Emerson BC, Paradis E, Thebaud C: Revealing the demographic histories of species using DNA sequences. TREE. 2001, 16: 707-716.Google Scholar
- Delport W, Crowe TM, Lloyd P, Bloomer P: Population growth confounds phylogeographic inference in Namaqua sandgrouse. J Hered. 2007, 98: 158-64. 10.1093/jhered/esl058.View ArticlePubMedGoogle Scholar
- Avise J: Molecular Markers, Natural History and Evolution. 1994, Chapman & Hall: New YorkView ArticleGoogle Scholar
- Taberlet P, Fumagalli L, Wust-Saucy AG, Cosson JF: Comparative phylogeography and postglacial colonization routes in Europe. Mol Ecol. 1998, 7: 453-464. 10.1046/j.1365-294x.1998.00289.x.View ArticlePubMedGoogle Scholar
- Blondel J, Aronson J: Biology and Wildlife of the Mediterranean Region. 1999, Oxford: Oxford University PressGoogle Scholar
- Taylor JW: Monograph of the land and freshwater mollusca of the British Isles. Edited by: Taylor Brothers. Leeds. 1913, 236-273.Google Scholar
- Guiller A, Madec L, Daguzan J: Geographical patterns of genetic differentiation in the landsnail Helix aspersa Müller (Gastropoda: Pulmonata). J Mollus Stud. 1994, 60: 205-221. 10.1093/mollus/60.3.205.View ArticleGoogle Scholar
- Guiller A, Bellido A, Coutelle A, Madec L: Spatial Genetic Pattern in the Land Mollusc Helix aspersa inferred from a "Center-based Clustering" procedure. Genet Res. 2006, 88: 1-18. 10.1017/S0016672306008305.View ArticleGoogle Scholar
- Guiller A, Bellido A, Madec L: Ordination and genetic distances: the land snail Helix aspersa in north Africa as a test case. Syst Biol. 1998, 47: 208-227. 10.1080/106351598260888.View ArticlePubMedGoogle Scholar
- Madec L, Bellido A, Guiller A: Statistical and biogeographical significances of patterns of morphological and biochemical variation in the land snail Helix aspersa. C R Acad Sci Paris Ser III. 1996, 319: 225-229.Google Scholar
- Madec L, Bellido A, Guiller A: Shell shape of the land snail Cornu aspersum in north Africa: further evidence of a phylogeographical splitting. Heredity. 2003, 91: 224-231. 10.1038/sj.hdy.6800301.View ArticlePubMedGoogle Scholar
- Madec L, Guiller A: Geographic variation of distal genitalia in the landsnail Helix aspersa (Mollusca: Gastropoda). J Zool. 1994, 233: 215-231.View ArticleGoogle Scholar
- Guiller A, Coutellec MA, Madec L, Deunff J: Evolutionary history of the land snail Helix aspersa in western Mediterranean: preliminary results inferred from mitochondrial DNA sequences. Mol Ecol. 2001, 10: 81-87. 10.1046/j.1365-294X.2001.01145.x.View ArticlePubMedGoogle Scholar
- Guiller A, Coutellec-Vreto MA, Madec L: Genetic relationships among suspected contact zone populations of Helix aspersa (Gastropoda: Pulmonata) in Algeria. Heredity. 1996, 74: 113-129. 10.1038/hdy.1996.116.View ArticleGoogle Scholar
- Van Andel TH: The climate and landscape of the middle part of the Weichselian glaciation in Europe: the stage 3 project. Quatern Research. 2002, 57: 2-8. 10.1006/qres.2001.2294.View ArticleGoogle Scholar
- Jakob SS, Ihlow A, Blattner FR: Combined ecological niche modelling and molecular phylogeography revealed the evolutionary history of Hordeum marinum (Poaceae) - niche differentiation, loss of genetic diversity, and speciation in Mediterranean Quaternary refugia. Mol Ecol. 2007, 16: 1713-1727. 10.1111/j.1365-294X.2007.03228.x.View ArticlePubMedGoogle Scholar
- Rokas A, Melika G, Abe Y, Nieves-Aldrey JL, Cook JM, Stone GN: Lifecycle closure, lineage sorting and hybridization revealed in a phylogenetic analysis of European oak gallwasps (Hymenoptera: Cynipidae: Cynipini) using mitochondrial sequence data. Mol Phylogenet Evol. 2003, 26: 36-45. 10.1016/S1055-7903(02)00329-9.View ArticlePubMedGoogle Scholar
- Pamilo P, Nei M: Relationships between gene trees and species trees. Mol Biol Evol. 1988, 5: 568-583.PubMedGoogle Scholar
- Maddison WP, Knowles LL: Inferring phylogeny despite incomplete lineage sorting. Syst Biol. 2006, 55: 21-30. 10.1080/10635150500354928.View ArticlePubMedGoogle Scholar
- Omland KE, Baker JM, Peters JL: Genetic signatures of intermediate divergence: population history of Old and New World Holarctic ravens (Corvus corax). Mol Ecol. 2006, 15: 795-808. 10.1111/j.1365-294X.2005.02827.x.View ArticlePubMedGoogle Scholar
- Hare MP, Cipriano F, Palumbi SR: Genetic Evidence on the Demography of Speciation in Allopatric Dolphin Species. Evolution. 2002, 56: 804-816.View ArticlePubMedGoogle Scholar
- Peters JL, Gretes W, Omland KE: Late Pleistocene divergence between eastern and Western populations of wood ducks (Aix sponsa) inferred by the 'isolation with migration' coalescent method. Mol Ecol. 2005, 14: 3407-3418. 10.1111/j.1365-294X.2005.02618.x.View ArticlePubMedGoogle Scholar
- Joseph L, Wilke T, Ten Have J, Chesser RT: Implications of mitochondrial DNA polyphyly in two ecologically undifferentiated but morphologically distinct migratory birds, the masked and white-browed woodswallows Artamus spp. of inland Australia. J Avian Biol. 2006, 37: 625-636. 10.1111/j.0908-8857.2006.03767.x.View ArticleGoogle Scholar
- Moore WS: Inferring phylogenies from mtDNA variation: mitochondrial-gene trees versus nuclear-gene trees. Evolution. 1995, 49: 718-726. 10.2307/2410325.View ArticleGoogle Scholar
- Lacroix C: Origine de la différenciation moléculaire et morphologique chez l'escargot Helix aspersa: résultats préliminaires. Master thesis. 2006, Rennes UniversityGoogle Scholar
- Madec L, Guiller A: Observations on genital structure and mating activity in three conchologically distinct forms of the land snail Helix aspersa Müller. J Mollus Stud. 1993, 59: 455-460. 10.1093/mollus/59.4.455.View ArticleGoogle Scholar
- Beerli P, Hotz E, Uzzell T: Geologically dated sea barriers calibrate a protein clock for Aegean waterfrogs. Evolution. 1996, 50: 1676-1687. 10.2307/2410903.View ArticleGoogle Scholar
- Ravelo AC, Andreasen DH, Lyle M, Lyle AO, Wara MW: Regional climate shifts caused by gradual global cooling in the Pliocene epoch. Nature. 2004, 429: 263-267. 10.1038/nature02567.View ArticlePubMedGoogle Scholar
- Rosenbaum G, Lister GS, Duboz C: Reconstruction of the tectonic evolution of the western Mediterranean since the Oligocene. In Reconstruction of the evolution of the Alpine-Himalayan orogeny. J Virtual Explorer. Edited by: Rosenbaum G, Lister GS. 2002, 8: 107-130.Google Scholar
- Schettino A, Turco E: Plate kinematics of the Western Mediterranean region during the Oligocene and Early Miocene. Geophysical Journal International. 2006, 166: 1398-1423. 10.1111/j.1365-246X.2006.02997.x.View ArticleGoogle Scholar
- Magri D, Fineschi S, Bellarosa R, Buonamici A, Sebastiana F, Schirone B, Simeone MC, Vendramin GG: The distribution of Quercus suber chloroplast haplotypes matches the palaeogeographical history of the western Mediterranean. Mol Ecol. 2007, 16: 5259-5266. 10.1111/j.1365-294X.2007.03587.x.View ArticlePubMedGoogle Scholar
- Burban C, Petit RJ: Phylogeography of maritime pine inferred with organelle markers having contrasted inheritance. Mol Ecol. 2003, 12: 1487-1495. 10.1046/j.1365-294X.2003.01817.x.View ArticlePubMedGoogle Scholar
- Fromhage L, Vences M, Veith M: Testing alternative vicariance scenarios in western Mediterranean discoglosside frogs. Mol Phylogenet Evol. 2004, 31: 308-322. 10.1016/j.ympev.2003.07.009.View ArticlePubMedGoogle Scholar
- Zangari F, Cimmaruta R, Nascetti G: Genetic relationships of the western Mediterranean painted frogs based on allozymes and mitochondrial markers: evolutionary and taxonomic inferences (Amphibia, Anura, Discoglossidae). Biol J Linn Soc. 2006, 87: 515-536. 10.1111/j.1095-8312.2006.00585.x.View ArticleGoogle Scholar
- Gantenbein B, Largiader CR: The phylogeographic importance of the Strait of Gibraltar as a gene flow barrier in terrestrial arthropods: a case study with the scorpion Buthus occitanus as model organism. Mol Phylogenet Evol. 2003, 28: 119-130. 10.1016/S1055-7903(03)00031-9.View ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1365-2699.2007.01688.x.View ArticleGoogle Scholar
- 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 Scripta. 2006, 35: 97-108. 10.1111/j.1463-6409.2005.00218.x.View ArticleGoogle Scholar
- Paulo OS, Pinheiro J, Miraldo A, Bruford MW, Jordan WC, Nichols RA: The role of vicariance vs. dispersal in shaping genetic patterns in ocellated lizard species in the western Mediterranean. Mol Ecol. 2008, 17: 1535-1551. 10.1111/j.1365-294X.2008.03706.x.View ArticlePubMedGoogle Scholar
- Veith M, Mayer C, Samraoui B, Donaire Barroso D, Bogaerts S: From Europe to Africa and vice versa: evidence for multiple intercontinental dispersal in ribbed salamanders (genus Pleurodeles). J Biogeogr. 2004, 31: 159-171. 10.1111/j.1365-2699.2004.00957.x.View ArticleGoogle Scholar
- Carranza S, Arnold EN: History of west Mediterranean newts, Pleurodeles (Amphibia: Salamandridae), inferred from old and recent DNA sequences. Syst Biodiv. 2004, 1: 327-337. 10.1017/S1477200003001221.View ArticleGoogle Scholar
- Carranza S, Wade E: Taxonomic revision of Algero-Tunisian Pleurodeles (Caudata: Salamandridae) using molecular and morphological data. Revalidation of the taxon Pleurodeles nebulosus (Ghichenot, 1850). Zootaxa. 2004, 488: 1-24.Google Scholar
- Brown RP, Suárez NM, Pestano J: The Atlas mountains as a biogeographical divide in north-west Africa: evidence from mtDNA evolution in the Agamid lizard Agama impalearis. Mol Phylogenet Evol. 2002, 24: 324-332. 10.1016/S1055-7903(02)00218-X.View ArticlePubMedGoogle Scholar
- Carranza S, Arnold EN, Pleguezuelos JM: Phylogeny, biogeography and evolution of two Mediterranean snakes, Malpolon monspessulanus and Hemorrhois hippocrepis (Squamata: Colubriade), using mtDNA sequences. Mol Phylogenet Evol. 2006, 46: 532-546. 10.1016/j.ympev.2006.03.028.View ArticleGoogle Scholar
- Barata M, Harris DJ, Castilho R: Comparative phylogeography of northwest African Natrix maura (Serpentes: Colubridae) inferred from mtDNA sequences. African Zool. 2008, 43: 1-7. 10.3377/1562-7020(2008)43[1:CPONAN]2.0.CO;2.View ArticleGoogle Scholar
- Guicking D, Joger U, Wink M: Molecular phylogeography of the viperine snake Natrix maura (Serpentes: Colubridae): evidence for strong intraspecific differentiation. Org Divers Evol. 2008, 8: 130-145. 10.1016/j.ode.2007.05.001.View ArticleGoogle Scholar
- Escoriza D, Comas MM, Donaire D, Carranza S: Rediscovery of Salamandra algira Bedriaga, 1833 from the Beni Snassen massif (Morocco) and phylogenetic relationships of north African Salamandra. Amph Rept. 2006, 27: 448-455. 10.1163/156853806778190042.View ArticleGoogle Scholar
- Steinfartz S, Veith M, Tautz D: Mitochondrial sequence analysis of Salamandra taxa suggests old splits of major lineages and postglacial recolonizations of Central Europe from distinct source populations of Salamandra salamandra. Mol Ecol. 2000, 9: 397-410. 10.1046/j.1365-294x.2000.00870.x.View ArticlePubMedGoogle Scholar
- Véla E, Magnin F, Pavon D, Pfenninger M: Phylogénie moléculaire et données paléobiogéographiques sur le gastéropode terrestre Tudorella sulcata (Draparnaud, 1805) en France et en Algérie orientale. Geodiv. 2008, 30: 233-246.Google Scholar
- Pfenninger M, Véla E, Ruth J, Arantzazu Elejalde M, Liberto F, Magnin F, Martinez-Orti A: Temporal speciation in the western Mediterranean genus Tudorella P. Fischer, 1885 (Gastropoda, Pomatiidae) supports the Tyrrhenian vicariance hypothesis. Mol Phylogenet Evol. 2009Google Scholar
- Harris DJ, Batista V, Carretero MA: Assessment of genetic diversity within Acanthodactylus erythrurus (Reptilia: Lacertidae) in Morocco and the Iberian Peninsula using mitochondrial DNA sequence data. Amph Rept. 2004, 25: 227-232. 10.1163/1568538041231229.View ArticleGoogle Scholar
- Harris DJ, Carranza S, Arnold EN, Pinho C, Ferrand N: Complex biogeographical distribution of genetic variation within Podarcis wall lizards across the Strait of Gibraltar. J Biogeogr. 2002, 29: 1257-1262. 10.1046/j.1365-2699.2002.00744.x.View ArticleGoogle Scholar
- Pinho C, Ferrand N, Harris DJ: Reexamination of the Iberian and north African Podarcis (Squamata: Lacertidae) phylogeny based on increased mitochondrial DNA sequencing. Mol Phylogenet Evol. 2006, 38: 266-273. 10.1016/j.ympev.2005.06.012.View ArticlePubMedGoogle Scholar
- 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. 10.1111/j.1365-294X.2005.02476.x.View ArticlePubMedGoogle Scholar
- Arano B, Llorente GA, Montori A, Buckley D, Herrero P: Diversification in north-west African water frogs: molecular and morphological evidence. Herpetol J. 1998, 8: 57-64.Google Scholar
- Carranza S, Harris DJ, Arnold EN, Batista V, Gonzalez JP: Phylogeography of the lacertid lizard, Psammodromus algirus, in Iberia and accross the Strait of Gibraltar. J Biogeogr. 2006, 33: 1279-1288. 10.1111/j.1365-2699.2006.01491.x.View ArticleGoogle Scholar
- Harris DJ, Batista V, Carretero MA: Diversity of 12S mitochondrial DNA sequence in Iberian and northwest African water frogs across predicted geographic barriers. Herpetozoa. 2003, 16: 81-83.Google Scholar
- Guillaumet A, Pons JM, Godelle B, Crochet PA: History of the Crested Lark in the Mediterranean region as revealed by mtDNA sequences and morphology. Mol Phylogenet Evol. 2006, 39: 645-56. 10.1016/j.ympev.2006.01.002.View ArticlePubMedGoogle Scholar
- 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-1205. 10.1111/j.1365-294X.2007.03230.x.View ArticlePubMedGoogle Scholar
- Harris DJ, Batista V, Carretero MA, Ferrand N: Genetic variation in Tarentola mauritanica (Reptilia: Gekkonidae) across the Strait of Gibraltar derived from mitochondrial and nuclear DNA sequences. Amph Rept. 2004, 25: 451-459.View ArticleGoogle Scholar
- Harris DJ, Batista V, Lymberakis P, Carretero MA: Complex estimates of evolutionary relationships in Tarentola mauritanica derived from mitochondrial DNA sequences. Mol Phylogenet Evol. 2004, 30: 855-859. 10.1016/S1055-7903(03)00260-4.View ArticlePubMedGoogle Scholar
- Frizon de Lamotte D, Saint-Bezar B, Bracène R, Mercier E: The two steps of the Atlas montains building and the geodynamics of the western mediterranean region. Tectonics. 2000, 19: 740-761. 10.1029/2000TC900003.View ArticleGoogle Scholar
- De Jong H: In search of historical biogeographic patterns in the western Mediterranean terrestrial fauna. Biol J Linn Soc. 1998, 65: 99-16. 10.1006/bijl.1998.0239.View ArticleGoogle Scholar
- Blondel J: Biogéographie évolutive. 1986, Masson. ParisGoogle Scholar
- Petit RJ, Brewer S, Bordács S, Burg K, Cheddadi R, Coart E, Cottrell J, Csaikl UM, van Dam BC, Deans JD, Fineschi S, Finkeldey R, Glaz I, Goicoechea PG, Jensen JS, König AO, Lowe AJ, Madsen SF, Mátyás G, Munro RC, Popescu F, Slade D, Tabbener H, de Vries SMG, Ziegenhagen B, de Beaulieu JL, Kremer A: Identification of refugia and postglacial colonization routes of European white oaks based on chloroplast DNA and fossil pollen evidence. Forest Ecol Management. 2002, 156: 49-74. 10.1016/S0378-1127(01)00634-X.View ArticleGoogle Scholar
- Aouraghe H: Histoire du peuplement paléolithique de l'Afrique du Nord et dynamique des interactions entre l'homme et son environnement. C R Palevol. 2006, 5: 237-242. 10.1016/j.crpv.2005.09.017.View ArticleGoogle Scholar
- Ray N, Adams JM: A GIS-based Vegetation Map of the World at the Last Glacial Maximum (25,000-15,000 BP). Internet Archaeol. 2001, 11: [http://0-intarch.ac.uk.brum.beds.ac.uk/journal/issue11/rayadams_toc.html]Google Scholar
- Gómez A, Lunt DH: Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. Phylogeography of southern European Refugia. Edited by: Weiss S, Ferrand N. 2007, Springer, 155-188. full_text.View ArticleGoogle Scholar
- Véla E, Benhouhou S: Assessment of a new hotspot for plant biodiversity in the Mediterranean basin (north Africa). C R Acad Sci Paris Ser III. 2007, 330: 589-605.Google Scholar
- Madec L: Genetic divergence in natural populations of the landsnail Helix aspersa Müller, 1774. J Mollus Stud. 1991, 57: 483-487. 10.1093/mollus/57.4.483.View ArticleGoogle Scholar
- Thomaz D, Guiller A, Clarke B: Extreme divergence of mitochondrial DNA within species of pulmonate land snails. Proc R Soc Lond B. 1996, 263: 363-368. 10.1098/rspb.1996.0056.View ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Gutell RR, Fox GE: A compilation of large subunit RNA sequences presented in a structural format. Nucleic Acids Res Suppl. 1988, 16: r175-r269.View ArticleGoogle Scholar
- Rozas J, Sanchez-Delbarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin version 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online. 2005, 1: 47-50.PubMed CentralGoogle Scholar
- Swofford D: PAUP*: phylogenetic analysis using parsimony (and other methods). Version 4. 2002, Sinauer Associates, Sunderland, MAGoogle Scholar
- Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian phylogenetic analysis of combined data. Syst Biol. 2004, 53: 47-67. 10.1080/10635150490264699.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck J: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: Intraspecific gene genealogies: trees grafting into networks. TREE. 2001, 16: 37-45.PubMedGoogle Scholar
- Bandelt HJ, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.View ArticlePubMedGoogle Scholar
- Felsenstein J: Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981, 17: 368-376. 10.1007/BF01734359.View ArticlePubMedGoogle Scholar
- Tajima F: Measurement of DNA polymorphism. Mechanisms of Molecular Evolution. Edited by: Takahata N, Clark AG. 1993, Sinauer Associates, MA, 37-59.Google Scholar
- Robinson-Rechavi M, Huchon D: RRTree: relative-rate tests between groups of sequences on a phylogenetic tree. Bioinformatics. 2000, 16: 296-297. 10.1093/bioinformatics/16.3.296.View ArticlePubMedGoogle Scholar
- Beerli P, Felsenstein J: Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach. Proc Natl Acad Sci USA. 2001, 98: 4563-4568. 10.1073/pnas.081068098.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.PubMed CentralPubMedGoogle Scholar
- Ramos-Onsins SE, Rozas J: Statistical properties of new neutrality tests against population growth. Mol Biol Evol. 2002, 19: 2092-2100.View ArticlePubMedGoogle Scholar
- Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9: 552-569.PubMedGoogle Scholar
- Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005, 22: 1185-1192. 10.1093/molbev/msi103.View ArticlePubMedGoogle Scholar
- Strimmer K, Pybus OG: Exploring the demographic history of DNA sequences using the generalized skyline plot. Mol Biol Evol. 2001, 18: 2298-2305.View ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralView ArticlePubMedGoogle Scholar
- Rambaut A, Drummond AJ: Tracer v1.4. 2007, [http://beast.bio.ed.ac.uk/Tracer]Google Scholar
- Chiba S: Accelerated evolution of land snails Mandarina in the oceanic Bonin islands: Evidence from mitochondrial DNA sequences. Evolution. 1999, 53: 460-471. 10.2307/2640782.View ArticleGoogle Scholar
- Davison A: Land snails as a model to understand the role of history and selection in the origins of biodiversity. Pop Ecol. 2002, 44: 129-136. 10.1007/s101440200016.View ArticleGoogle Scholar
- Van Riel P, Jordaens K, Van Houtte N, Frias Martins AM, Verhagen R, Backeljau T: Molecular systematics of the endemic Leptaxini (Gastropoda: Pulmonata) on the Azores islands. Mol Phylogenet Evol. 2005, 37: 132-143. 10.1016/j.ympev.2005.03.019.View ArticlePubMedGoogle Scholar
- Marko PB: Fossil Calibration of Molecular Clocks and the Divergence Times of Geminate Species Pairs Separated by the Isthmus of Panama. Mol Biol Evol. 2002, 19: 2005-2021.View ArticlePubMedGoogle Scholar
- Neiman M, Lively CM: Pleistocene glaciation is implicated in the phylogeographical structure of Potamopyrgus antipodarum, a New Zealand snail. Mol Ecol. 2004, 13: 3085-3098. 10.1111/j.1365-294X.2004.02292.x.View ArticlePubMedGoogle Scholar
- Johnson SG: Age, phylogeography and population structure of the microendemic banded spring snail, Mexipyrgus churinceanus. Mol Ecol. 2005, 14: 2299-2311. 10.1111/j.1365-294x.2005.02580.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.