- Research article
- Open Access
From Africa to Europe and back: refugia and range shifts cause high genetic differentiation in the Marbled White butterfly Melanargia galathea
© Habel et al; licensee BioMed Central Ltd. 2011
- Received: 27 January 2011
- Accepted: 21 July 2011
- Published: 21 July 2011
The glacial-interglacial oscillations caused severe range modifications of biota. Thermophilic species became extinct in the North and survived in southern retreats, e.g. the Mediterranean Basin. These repeated extinction and (re)colonisation events led to long-term isolation and intermixing of populations and thus resulted in strong genetic imprints in many European species therefore being composed of several genetic lineages. To better understand these cycles of repeated expansion and retraction, we selected the Marbled White butterfly Melanargia galathea. Fourty-one populations scattered over Europe and the Maghreb and one population of the sibling taxon M. lachesis were analysed using allozyme electrophoresis.
We obtained seven distinct lineages applying neighbour joining and STRUCTURE analyses: (i) Morocco, (ii) Tunisia, (iii) Sicily, (iv) Italy and southern France, (v) eastern Balkans extending to Central Europe, (vi) western Balkans with western Carpathian Basin as well as (vii) south-western Alps. The hierarchy of these splits is well matching the chronology of glacial and interglacial cycles since the Günz ice age starting with an initial split between the galathea group in North Africa and the lachesis group in Iberia. These genetic structures were compared with past distribution patterns during the last glacial stage calculated with distribution models.
Both methods suggest climatically suitable areas in the Maghreb and the southern European peninsulas with distinct refugia during the last glacial period and underpin strong range expansions to the North during the Postglacial. However, the allozyme patterns reveal biogeographical structures not detected by distribution modelling as two distinct refugia in the Maghreb, two or more distinct refugia at the Balkans and a close link between the eastern Maghreb and Sicily. Furthermore, the genetically highly diverse western Maghreb might have acted as source or speciation centre of this taxon, while the eastern, genetically impoverished Maghreb population might result from a relatively recent recolonisation from Europe via Sicily.
- climatic oscillations
- Melanargia galathea
- Melanargia lachesis
- allozyme electrophoresis
- climate envelope modelling
The impacts of climatic oscillations on the earth's biota have been intensively studied . In the western Palaearctic, thermophilic organisms went extinct over major parts of Central and North Europe during cold stages and survived in the lowlands of lower latitudes in often distinct refugia [2–6]. Molecular studies revealed that most of these taxa exclusively survived glacial periods south of the European high mountain chains in the Iberian, Italian and Balkan peninsulas, and some even in additional extra-Mediterranean refugia [7, 8]. The long-term isolation of populations in these retreats over many thousands of years resulted in genetic differentiation . During the warmer interglacial periods, species expanded their distribution ranges northwards and extended their different genetic lineages over more northern areas [9, 10].
In contrast to the three more intensively studied Mediterranean refugia of southern Europe (Iberia, peninsular Italy and the Balkans), little is known about North African refugia and the biogeographical relation between the Maghreb and southern Europe separated by the two narrow sea straits of Gibraltar and Sicily. It has been shown that the Maghreb is often sub-structured following an east-west [e.g. [11–13]] or south-north differentiation pattern [e.g. [5, 14]]; in some cases, genetic continuity was demonstrated between the Maghreb and Sicily [e.g. [15, 16]]. Other studies underline the important role of Sicily as diversification centre for European taxa unravelling deep genetic splits between this island and peninsular Italy (e.g. Erinaceus europaeus: ; Pseudepidalea viridis: [18, 19]). Few molecular analyses also reveal the outstanding importance of North Africa as a refugium for thermophilic species during glacial periods [e.g. [12, 14, 20, 21]]. However, most studies focus either on the Maghreb or the southern European refugia and do not combine the distribution of species all over north-western Africa and throughout Europe.
To study the biogeographical importance of the Maghreb region and its connection with Europe, we selected the Marbled White butterfly species complex Melanargia galathea (Linnaeus, 1758) and Melanargia lachesis (Hübner, 1790) as a model system using two analytical tools (allozyme polymorphisms and distribution modelling). Today, M. galathea is widely distributed from the Maghreb region (mountain ranges of Morocco, Algeria and Tunisia) [22, 23] to the English Midlands , and from the Pyrenees  to the Baltic Sea in Poland . On the Iberian Peninsula, M. galathea is replaced by its sibling species M. lachesis. Thus, the Italian peninsula is the only possible link between North Africa and Europe for M. galathea.
Which refugia are of importance for the glacial survival of the M. galathea /lachesis species complex during the subsequent glacial periods?
Is there any evidence of genetic structuring within the North African and Italian refugia?
Which routes of expansion and retraction followed the butterfly throughout time?
All enzyme loci had banding patterns consistent with known quaternary structures. While most loci were inherited autosomally, 6PGDH and ME were located on the Z chromosome so that hemizygous females (but not males) had a single copy . No general linkage disequilibrium was observed for any locus (all p > 0.05 after Bonferoni correction). A total of 13 analysed loci were polymorphic, but two loci (FUM, GPDH) were monomorphic throughout all samples. Allele frequencies for each enzyme and population are given in an additional file 1.
Sampling location and five parameters of genetic diversity for 41 populations of Melanargia galathe a from its western Palaearctic distribution area and one population of M. lachesis from the Pyrenees: number of individuals analysed (N), mean number of alleles per locus (A), percentage of expected (H e ) and observed (H o ) heterozygosity, percentage of polymorphic loci not exceeding 95% (P95) and total number of polymorphic loci (Ptot)
Date of sampling
H e (%)
H o (%)
T-Table de Yagurta
F-Col de Tende
RO-Porta di Fier Transilvanici
E-Col de Perbes*
Means of sample sizes and genetic diversities of the different genetic groups of Melanargia galathea and M. lachesis; p values of Kruskal Wallis ANOVAs among groups are given
36.0 ± 7.9
2.10 ± 0.19
17.0 ± 2.2
15.3 ± 2.8
34.5 ± 5.4
45.2 ± 9.4
36.0 ± 0.0
2.38 ± 0.10
18.4 ± 1.6
17.6 ± 3.0
36.7 ± 6.7
53.3 ± 5.4
37.4 ± 7.8
2.00 ± 0.17
13.9 ± 2.1
11.6 ± 2.4
33.3 ± 6.7
49.3 ± 8.9
37.7 ± 5.8
2.13 ± 0.13
15.8 ± 1.8
13.4 ± 1.3
33.3 ± 6.6
55.6 ± 10.2
Italy + SE France
35.3 ± 8.3
2.19 ± 0.24
17.9 ± 2.3
15.1 ± 2.6
38.1 ± 8.3
52.4 ± 11.2
eastern Balkans + Central Europe
36.8 ± 7.4
2.04 ± 0.15
17.8 ± 1.7
16.9 ± 2.0
33.7 ± 2.9
38.7 ± 3.6
26.5 ± 15.6
2.01 ± 0.22
15.3 ± 1.3
14.4 ± 1.1
31.6 ± 3.3
40.0 ± 5.5
p (Kruskal Wallis ANOVA)
eastern Balkans + Romania
35.7 ± 8.4
2.05 ± 0.15
17.1 ± 2.2
15.5 ± 1.8
33.3 ± 3.8
39.0 ± 4.6
37.7 ± 7.0
2.03 ± 0.16
18.5 ± 1.1
18.0 ± 1.4
34.0 ± 2.2
38.5 ± 2.9
p (Kruskal Wallis ANOVA)
Analyses of molecular variance for all Melanargia galathea populations and one population of M.lachesis
among individuals within populations
All galathea and lachesis
Sicily and Tunisia
Italy, Balkans, Central Europe
Italy and S France
west Balkan group
east Balkan and Central Europe
Hierarchical variance analyses of Melanargia galathea and M.lachesis among genetic groups
prop. of among groups variance of total variance among pops.
galathea vs lachesis
Tunisia + Sicily vs rest of Europe
Tunisia vs Morocco
Tunisia vs Sicily
Sicily vs Italy
Italy vs S France (Condat)
Italy vs SW Alps (Col de Tende)
east Balkans + Central Europe vs SW Alps (Col de Tende)
Italy vs SW Alps vs west Balkans vs east Balkans + Central Europe
east Balkans + Romania vs Central Europe
Morocco vs continental Europe
The genetic diversities among these genetic lineages showed significant differences (Table 2). Thus, the Morocco group showed the highest values achieved for A, He as well as Ho, and the means for P95 and Ptot were above average. On the other extreme, Tunisia had the lowest means for A, He and Ho, and the mean for P95 was well below average; the genetic diversities of Tunisia were lower (A, He, Ho, Ptot) or equal than in the otherwise rather similar populations from Sicily. The four groups from mainland Europe all have mostly intermediate genetic diversities scatted around the respective mean values.
Species distribution modelling
According to the classification of Swets , we received 'excellent' AUC values in our 100 models (average training AUC = 0.927, average test AUC = 0.902). On average, the 'temperature annual range' had the highest explanatory power (30.3%), followed by the 'minimum temperature of the coldest month' (16.8%), the 'precipitation of the warmest quarter' (14.8%), the 'annual precipitation' (10.2%), the 'maximum temperature of the warmest month' (8.2%) and the 'precipitation of the driest quarter' (7.2%). All other variables contributed less than 5% each. The average minimum training presence was 0.05, and the lowest 10 percentile training omission threshold was 0.36.
The obtained allozyme data displayed in neighbour-joining phenograms, structure plots and hierarchical variance analyses indicate a profound genetic split between the two taxa, M. galathea and M. lachesis. Nazari et al.  supported this pattern by three lines of evidence: (i) differences of the male genitalia between M. lachesis and M. galathea, (ii) a stronger difference in wing patterns between these two taxa than between M. galathea population in Europe and the Maghreb and (iii) remarkable differences in DNA sequences of the nuclear wg gene between M. lachesis and M. galathea, but no major differentiation between M. galathea samples from Europe and the Maghreb. However, the sequences of the two mtDNA genes cox1 and 16S contradict the common pattern of allozymes, genital structures, wing patterns and nuclear DNA sequences: This marker is not well distinguishing M. galathea from Europe and M. lachesis, but shows remarkable differences between Europe and the Maghreb with this split being dated back to the Messinia Salinity Crises more than 5 My ago . Having in mind the differentiation pattern in all known marker systems, we believe that these two mtDNA lineages in the entire species complex might have originated at that time horizon, but were distributed to different geographical regions only much later by lineage sorting, maybe hereby exemplifying one case of the often observed difference between mtDNA on the one hand and nuclear DNA sequences, morphological characteristics and allozyme pattern on the other .
Our allozyme data further show strong differentiation within M. galathea into two major groups with respective subgroups: (i) Sicily - Tunisia with (i-a) Sicily and (i-b) Tunisia as well as (ii) all other M. galathea with (ii-a) Morocco, (ii-b) Italy with parts of southern France, (ii-c) western Balkan including the western Carpathian Basin, (ii-d) eastern Balkans with Romania and Central Europe, and (ii-e) the south-western Alps.
Atlantic-Mediterranean origin of the M. galathea/lachesisspecies complex
From the Maghreb to Europe
The deepest split in the M. galathea populations is between the Sicily - Tunisia group and all the other populations. As this split is about twice the genetic differentiation among their subgroups and less than half of the distance against M. lachesis, the onset of the Riss glaciation (about 310 ky BP)  might be the trigger for vicariance and thus the beginning of this differentiation. As (i) Iberia was continuously blocked for the expansion of M. galathea to Europe by M. lachesis [cf. 27] and (ii) all European M. galathea populations except Sicily are more similar to populations from Morocco than from Tunisia, a scenario with this split taking place in the Maghreb is little likely. This assumption is further supported by SDMs for ice age conditions predicting mostly continuous distributions over North Africa (Figure 3b) thus allowing vicariance in this region only during the relatively short interglacial stages. For these reasons, M. galathea must have reached Europe before the Riss glaciation.
As the region of the eastern Sahara in Egypt apparently always have been too dry for an expansion of M. galathea, this first expansion of M. galathea to Europe must have been from Tunisia to Sicily (Figure 4a), a sea strait known for biogeographical connections for many taxa [e.g. ; and references therein]. As the Strait of Sicily was considerably narrower during glacial periods due to eustatic sea level lowering, the transition from Mindel glaciation to Holstein interglacial with still low sea level but already higher temperatures might have been a suitable time period for this dispersal. After arrival to Sicily, the Holstein interglacial might have given suitable condition for the expansion of M. galathea over most parts of Europe, including the Balkans but excluding Iberia as this peninsula was already populated by M. lachesis (Figure 4a).
With the climatic cooling of the Riss ice age, which was considerably longer than the following Würm glaciation and had longer durations of minimum temperatures [33, 37], M. galathea most probably was nearly extinct in Europe only surviving in the southernmost possible retreats in Sicily and the southern Balkans (Peleponnesos), but also in the Maghreb; M. lachesis could survive in southern Iberia (Figure 4b). This vicariance might be the origin of the two major European lineages of M. galathea with the eastern one by chance evolving similarly in allele frequencies as the Morocco lineage, with this similarity therefore not representing recent biogeographical connection between them. Riss vicariance events most likely have also been responsible for other differentiation processes as e.g. in the Polyommatus coridon /hispana complex [e.g. ].
...and back to the Maghreb
As the time for differentiation between the four M. galathea lineages from continental Europe is assumed to be the result of one glacial cycle (see above) and as the differentiation between populations from Sicily and Tunisia are in the same order of magnitude, we assume that the onset of this differentiation is in the same time frame. As the genetic diversity is significantly higher in Sicily than in Tunisia and the warm and dry interglacial climatic conditions in Tunisia generally unsuitable for the survival of M. galathea, we assume that a colonisation most likely has taken place from Sicily to Tunisia. While the sea level was still considerably lowered at the transition from Riss to Eem thus facilitating dispersal between these two areas, this time period might be the most likely for this expansion event. During the following Eem interglacial, the Balkan refuge of M. galathea most probably could colonise most parts of Europe apart from Iberia and Sicily, which were occupied by other genetic lineages of this species complex (Figure 4b).
The existence of extra-Mediterranean refugia for thermophilic taxa
During the Würm ice age, which was not more severe than the two previous glaciations but with a shorter maximum , the Marbled White butterflies were not that much pushed to the South than in the previous cases. This is well matching the remarkable differentiation of the species in Europe allowing to distinguish five lineages (see above), which most likely are the result of survival of the Würm ice age in a larger number of different refugia.
This pattern implies at least two different refugia at the Balkan Peninsula at the western and the eastern flank; more in detail analyses also support a third Balkan centre in the peninsula's southern parts  (Figure 4c). This pattern of multiple refugia in the Balkans was already erected by Reinig  postulating different centres of survival in the western, southern and eastern Balkans and was later supported by genetic analyses showing genetic divergences between these areas for a variety of different animal species [e.g. [18, 40–42]].
Furthermore, different Würm refugia have to be postulated for Sicily and peninsular Italy, a pattern also repeated by other genetic analyses [e.g. [17, 43]]. Furthermore, other genetic studies show a remarkable genetic differentiation in the southernmost parts of peninsular Italy [e.g. [34, 44, 45]].
The last remaining lineage of M. galathea in the south-western Alps most likely is not representing a Mediterranean refuge of this species, but an extra-Mediterranean refuge area at the southern slopes of the glaciated Alps (Figure 4c). As already shown by Steward and Lister , glacial survival of temperate species in Europe was not only possible in the classical Mediterranean refugia sensu de Lattin , but also in small climatically buffered pockets in more northern regions [8, 48, 49]. Recent works especially highlight the southern and south-eastern parts of the Alps of particular importance for additional Würm ice age refugia for temperate species [e.g. [42, 50, 51]], and also for species formerly thought to be of exclusive Mediterranean origin [e.g. [52, 53]]. This apparently was also the case for the Marbled White.
During the Postglacial, several lineages of M. galathea were mostly blocked in their expansion by other lineages representing the respective leading edges [cf. 54]. In the case of M. galathea in Morocco, their northwards expansion was blocked by M. lachesis distributed in Iberia. The lineage surviving in the eastern Balkans apparently had the most important impact in the recolonisation of more northern parts of Europe as its dispersal was not hampered by any major mountain obstacle [cf. 9] so that this lineage could expand throughout Central Europe to the western parts of Germany (Figure 4c). However, the samples of north-eastern France and southern Germany show an intermediate genetic structure between this lineage and the south-western Alps lineage, making hybrid origin of these populations rather likely and thus expansion of the southern Alps lineage over the chains of the Alps.
Also the Italian lineage could expand beyond its refugium to southern France. Therefore the entire region of northern France and southern Germany might be a zone of mixing between these three lineages. Hybrid zones between different taxa are frequently observed in this region [e.g. [9, 55]]. Furthermore, the southernmost population in Calabria (southern Italy) has an intermediate genetic texture between the Italian and the Sicily group thus speaking for a postglacial contact and intermixing between these two groups in this region.
The hierarchical structure of our allozyme data set on M. galathea and M. lachesis is consistent with the chronology of the last four glacial-interglacial cycles. Based on this consistency, we derive the following scenario, which in our opinion is the most likely one: (i) The beginning of the Günz ice age might have affected the vicariance between the two species. (ii) M. galathea might have crossed from Tunisia to Sicily at the transition from Mindel ice age to Holstein interglacial and (iii) subsequently spread all over Europe, but retreated in the Maghreb to the higher elevations of the Atlas mountains. (iv) The members of this species complex survived the coldest periods of the Riss glaciation only in southern Iberia, Morocco, Sicily and the southern Balkans (Peleponnesos). (v) At the transition from Riss ice age to Eem interglacial, Tunisia was recolonised from Sicily. (vi) The southern Balkan group might have colonised major parts of Europe during the Eem interglacial including Italy and Central Europe. (vii) Populations of this group survived the Würm ice age in Italy, the southern margin of the Alps, the western and eastern flank of the Balkan peninsula; members of other lineages survived in Sicily, Tunisia, Morocco and Iberia. (viii) During the Postglacial, only the eastern Balkan and the Italian lineage showed major northwards range expansion. (ix) Hybridisation between lineages most probably occurred in western Central Europe and southern Calabria.
Alleles were labelled according to their relative mobility, starting with "1" for the slowest. All laboratory results were stored on cellulose acetate plates. These banding patterns were (re)analysed by one person (JCH). Allele frequencies, Nei's standard genetic distances  and parameters of genetic diversity (i.e. mean number of alleles per locus, A, expected heterozygosity, He, and observed heterozygosity, Ho, total percentage of polymorphic loci, Ptot, and percentage of polymorphic loci with the most common allele not exceeding 95%, P95) were computed with G-Stat . As sample sizes do not differ significantly, the calculation of allelic richness correcting for population sizes was not necessary. For detecting differences of means of genetic diversities among genetic lineages and sublineages, we calculated U-tests using STATISTICA. Conventional F statistics, AMOVAs, hierarchical genetic variance analysis, tests of Hardy-Weinberg equilibrium and linkage disequilibrium were calculated with ARLEQUIN 3.1 . Phenograms using the neighbour joining algorithm  were constructed with PHYLIP , including bootstrap-values (calculated based on 1,000 iterations). To define individual based genetic clusters we performed STRUCTURE analyses . As burn-in and simulation lengths we used 100,000 and 300,000 iterations per run based on the admixture model with correlated gene frequencies comparing different groupings (from K = 2 to K = 10).
Species Distribution Modelling
Over the last few decades, Geographic Information System (GIS) based Species Distribution Models (SDMs) have become vital tools used to predict the potential distribution of species under current conditions and climate change scenarios [62–64]. In combination with palaeoclimatological data, SDMs have been suggested as a mean of inferring species' past distributions [65, 66], especially when combined with phylogeographic techniques .
We compiled a set of 3,483 species records of M. galathea from online data bases (Global Biodiversity Information Facility - GBIF; http://www.gbif.org) and our own field surveys. The accuracy of all records was checked in DIVA-GIS 5.4  and only those which could be unambiguously assigned to a single grid cell with a resolution of 2.5 arc min (ca. 4 km in the study area) were used for further processing. Since unequal spatial clumping of species records may cause problems when computing SDMs, the species records were filtered in geographic space, leaving only 1 record per 10 arc min. The final data set comprised 535 records (Figure 3a) scattered all over the known range of the species in Europe and North Africa.
We obtained information on current and past climate as describedby the Community Climate System Model (CCSM; http://www.ccsm.ucar.edu) with a spatial resolution of 2.5 arc min from the Worldclim data base (; http://www.worldclim.org). Original palaeoclimatological data were previously processed as described by Peterson and Nyári (2007) . A total of 19 BIOCLIM variables were previously suggested as suitable for SDM computation [71, 72]. However, inclusion of too many inter-correlated variables or biologically irrelevant predictors may hamper the transferability of SDMs through space and time [73–76]. Therefore, we first computed a pair-wise correlation matrix based on Pearson's correlation coefficients among all 19 predictor variables and excluded those with R 2 > 0.75. Subsequently, we chose a final set of eleven predictors describing biologically relevant climate conditions for the long-term persistence of M. galathea populations (i.e. annual mean temperature, maximum temperature of warmest month, minimum temperature of coldest month, temperature annual range, mean temperature of wettest quarter, mean temperature of driest quarter, annual precipitation, precipitation of wettest quarter, precipitation of driest quarter, precipitation of warmest quarter, precipitation of coldest quarter).
For SDM computation, Maxent 3.3.2 was applied [77, 78] using the default program settings. Random background records were automatically sampled by Maxent within the study area. Species records were split 100 times into 70% used for model training and 30% for model evaluation via the area under the receiver operating characteristic curve (AUC; ). Subsequently, the average predictions for current and past conditions of the logistic output of the 100 models were computed and transformed into presence/absence maps applying the average minimum training presence and average 10% training omission as thresholds.
We acknowledge financial support by the German Academic Exchange Service (PostDoc Programme) and the Musée national d'histoire naturelle Luxembourg to JCH as well as from the Ministry of Education, Science, Youth and Culture of the Rhineland-Palatinate state of Germany to DR (project: 'Implications of global change for biological resources, law and standards'). We thank Claas Damken (Auckland, New Zealand) and Marc Meyer (Luxembourg) for field assistance. We thank the Fonds National de la Recherche Luxembourg for covering the publication fees.
- Lomolino MV: Biogeography. Sinauer Assn. 2005, 465-Google Scholar
- Reinig WF: Die Holarktis. Gustav-Fischer-Verlag. 1937, JenaGoogle Scholar
- De Lattin G: Grundriß der Zoogeographie. Verlag Gustav Fischer. 1967, JenaGoogle Scholar
- Dennis RLH, Williams WR, Shreeve TG: A multivariate approach to the determination of faunal units among European butterfly species (Lepidoptera: Papilionoidea, Hesperioidea). Zoological Journal of the Linnean Society. 1991, 101: 1-49. 10.1111/j.1096-3642.1991.tb00884.x.View ArticleGoogle Scholar
- Hewitt GM: Some genetic consequences of ice ages, and their role in divergence and speciation. Biological Journal of the Linnean Society. 1996, 58: 247-276.View ArticleGoogle Scholar
- Taberlet P, Fumagalli L, Wust-Saucy A-G, Cosson J-F: Comparative phylogeography and postglacial colonization routes in Europe. Molecular Ecology. 1998, 7: 453-464. 10.1046/j.1365-294x.1998.00289.x.View ArticlePubMedGoogle Scholar
- Hewitt GM: Genetic consequences of climatic oscillation in the quaternary. Phil Trans R Soc Lond B. 2004, 359: 183-195. 10.1098/rstb.2003.1388.View ArticleGoogle Scholar
- Schmitt T: Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Frontiers in Zoology. 2007, 4: 11-10.1186/1742-9994-4-11.View ArticlePubMedPubMed CentralGoogle Scholar
- Hewitt GM: Post-glacial re-colonization of European biota. Biological Journal of the Linnean Society. 1999, 68: 87-112. 10.1111/j.1095-8312.1999.tb01160.x.View ArticleGoogle Scholar
- Hewitt GM: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405: 907-913. 10.1038/35016000.View ArticlePubMedGoogle Scholar
- Cosson J-F, Hutterer R, Libois R, Sara M, Taberlet P, Vogel P: Phylogeographical footprint 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). Molecular Ecology. 2005, 14: 1151-1162. 10.1111/j.1365-294X.2005.02476.x.View ArticlePubMedGoogle Scholar
- Carranza S, Arnold EN, Pleguezuelos JM: Phylogeny, biogeography, and evolution of two Mediterranean snakes, Malpolon monspessulanus and Hemorrhois hippocrepis (Squamata, Colubridae), using mtDNA sequences. Molecular Phylogenetics and Evolution. 2006, 40: 532-546. 10.1016/j.ympev.2006.03.028.View ArticlePubMedGoogle Scholar
- 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. Journal of Biogeography. 2006, 33: 1279-1288. 10.1111/j.1365-2699.2006.01491.x.View ArticleGoogle Scholar
- Fritz U, Barata M, Busack SD, Fritsch 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). Zoologica Scripta. 2006, 35: 97-108. 10.1111/j.1463-6409.2005.00218.x.View ArticleGoogle Scholar
- Franck P, Garnery L, Loiseau A, Oldroyd BP, Hepburn HR, Solignac M, Cornuet J-M: Genetic diversity of the honeybee in Africa: microsatellite and mitochondrial data. Heredity. 2001, 86: 420-430. 10.1046/j.1365-2540.2001.00842.x.View ArticlePubMedGoogle Scholar
- Habel JC, Rödder D, Scalercio S, Meyer M, Schmitt T: Strong genetic cohesiveness between Italy and the Maghreb in four butterfly species. Biological Journal of the Linnean Society. 2010, 99: 818-830. 10.1111/j.1095-8312.2010.01394.x.View ArticleGoogle Scholar
- Santucci F, Emerson B, Hewitt GM: Mitochondrial DNA phylogeography of European hedgehogs. Molecular Ecology. 1998, 7: 1163-1172. 10.1046/j.1365-294x.1998.00436.x.View ArticlePubMedGoogle Scholar
- Seddon JM, Santucci F, Reeve NJ, Hewitt GM: DNA footprints of European hedgehogs, Erinaceus europaeus and E. concolor. Pleistocene refugia, postglacial expansion and colonization routes. Molecular Ecology. 2001, 10: 2187-2198. 10.1046/j.0962-1083.2001.01357.x.View ArticlePubMedGoogle Scholar
- Colliard C, Sicilia A, Turrisi GF, Arculeo M, Perrin N, Stöck M: Strong reproductive barriers in a narrow hybrid zone of West-Mediterranean green toads (Bufo viridis subgroup) with Plio-Pleistocene divergence. BMC Evolutionary Biology. 2010, 10: 232-10.1186/1471-2148-10-232.View ArticlePubMedPubMed CentralGoogle Scholar
- Veith M, Mayer C, Samraoui B, Barrosso DD, Bogaerts S: From Europe to Africa and vice versa: evidence for multipe intercontinental dispersal in ribbed salamanders (Genus Pleurodeles). Journal of Biogeography. 2004, 31: 159-171. 10.1111/j.1365-2699.2004.00957.x.View ArticleGoogle Scholar
- Paulo OS, Pinto I, Bruford MW, Jordan WC, Nichols RA: The double origin of Iberian peninsular chameleons. Biological Journal of the Linnean Society. 2008, 75: 1-7.View ArticleGoogle Scholar
- Tennent J: The butterflies of Morocco, Algeria and Tunisia. Gem publishing Company. 1996, WallingfordGoogle Scholar
- Tolman T, Lewington R: Field guide butterflies of Britain and Europe. Harper Collins Publishers. 1997, LondonGoogle Scholar
- Asher J, Warren M, Fox R, Harding P, Jeffcoate G, Jeffcoate S: The millennium atlas of butterflies in Britain and Ireland. 2001, Oxford University Press OxfordGoogle Scholar
- García-Barros E, Munguira ML, Martín Cano J, Romo Benito H, Garcia-Pereira P, Maravalhas ES: Atlas de las mariposas diurnas de la Península Ibérica e islas Baleares (Lepidoptera: Papilionoidea and Hesperioidea). Monografias Sociedad Entomológica. 2004, Aragonesa, 11:Google Scholar
- Buszko J: Atlas rozmieszczenia motyli dziennych w Polsce (Lepidoptera: Papilionidae, Hesperiidae). Edycja Turpress. 1997, TorunGoogle Scholar
- Habel JC, Schmitt T, Müller P: The fourth paradigm pattern of postglacial range expansion of European terrestrial species: The phylogeography of the Marbled White butterfly (Satyrinae, Lepidoptera). Journal of Biogeography. 2005, 32: 1489-1497. 10.1111/j.1365-2699.2005.01273.x.View ArticleGoogle Scholar
- Habel JC, Meyer M, El Mousadik A, Schmitt T: Africa goes Europa: The complete phylogeography of the Marbled White butterfly species complex Melanargia galathea /lachesis. Organisms, Diversity and Evolution. 2008, 8: 121-129. 10.1016/j.ode.2007.04.002.View ArticleGoogle Scholar
- Schmitt T, Habel JC, Zimmermann M, Müller P: Genetic differentiation of the marble white butterfly, Melanargia galathea, accounts for glacial distribution patterns and postglacial range expansion in southeastern Europe. Molecular Ecology. 2006, 15: 1889-1901. 10.1111/j.1365-294X.2006.02900.x.View ArticlePubMedGoogle Scholar
- Swets K: Measuring the accuracy of diagnostic systems. Science. 1988, 240: 1285-1293. 10.1126/science.3287615.View ArticlePubMedGoogle Scholar
- Nazari V, Hagen WT, Bozano GC: Molecular systematics and phylogeny of the Marbled Whites (Lepidoptera: Nymphidae, Satyrinae, Melanargia Meigen). Sysematics Entomology. 2009, 35: 132-147.View ArticleGoogle Scholar
- Hein J, Schierup MH, Winf C: Gene genealogies, variation and evolution - A primer in coalescent theory. 2005, Oxford University Press, OxfordGoogle Scholar
- Gibbard P, van Kolfschoten T: The Pleistocene and Holocene epochs. A geologic time scale. Edited by: Gradstein FM, Ogg JG, Smith AG. 2004, Cambridge University Press, Cambridge, 98: 441-452.Google 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. Molecular Ecology. 2000, 9: 397-410. 10.1046/j.1365-294x.2000.00870.x.View ArticlePubMedGoogle Scholar
- Carranza S, Wade E: Taxonomic revision of Algero-Tunisian Pleurodeles (Caudata: Salamandridae) using molecular and morphological data. Revalidation of the taxon Pleurodeles nebulosus (Guichenot, 1850). Zootaxa. 2004, 488: 1-24.View ArticleGoogle Scholar
- Fromhage L, Vences M, Veith M: Testing alternative vicariance scenarios in Western Mediterranean discoglossid frogs. Molecular Phylogenetics and Evolution. 2003, 31: 308-322.View ArticleGoogle Scholar
- Quante M: The changing climate: future. Relict species - phylogeography and conservation biology. Edited by: Habel JC, Assmann T. 2010, Springer Heidelberg, 9-56.Google Scholar
- Schmitt T, Varga Z, Seitz A: Are Polyommatus hispana and Polyommatus slovacus bivoltine Polyommatus coridon (Lepidoptera: Lycaenidae)? The discriminatory value of genetics in the taxonomy. Organisms, Diversity and Evolution. 2005, 5: 297-307. 10.1016/j.ode.2005.01.001.View ArticleGoogle Scholar
- Reinig WF: Chorologische Voraussetzungen für die Analyse von Formenkreisen. Syllegomena Biologica, Festschrift für O Kleinschmidt. 1950, 346-378.Google Scholar
- Lenk P, Fritz U, Joger U, Winks M: Mitochondrial phylogeography of the European pond turtle, Emys orbicularis (Linnaeus 1758). Molecular Ecology. 1999, 8: 1911-1922. 10.1046/j.1365-294x.1999.00791.x.View ArticlePubMedGoogle Scholar
- Seddon JM, Reeve N, Hewitt GM: Caucasus Mountains divide postulated postglacial colonization routes in the white-breasted hedgehog, Erinaceus concolor. Journal of Evolutionary Biology. 2002, 15: 463-467. 10.1046/j.1420-9101.2002.00408.x.View ArticleGoogle Scholar
- Pinceel J, Jordaens K, Pfenninger M, Backeljau T: Rangewide phylogeography of a terrestrial slug in Europe: evidence for Alpine refugia and rapid colonization after the Pleistocene glaciations. Molecular Ecology. 2005, 14: 1133-1150. 10.1111/j.1365-294X.2005.02479.x.View ArticlePubMedGoogle Scholar
- Cooper SJ, Ibrahim KM, Hewitt GM: Postglacial expansion and genome subdivision in the European Grasshopper Chorthippus parallelus. Molecular Ecology. 1995, 4: 49-60. 10.1111/j.1365-294X.1995.tb00191.x.View ArticlePubMedGoogle Scholar
- Podnar M, Mayer W, Tvrtkovic N: Phylogeography of the Italian wall lizard, Podarcis sicula, as revealed by mitochondrial DNA sequences. Molecular Ecology. 2005, 14: 575-588. 10.1111/j.1365-294X.2005.02427.x.View ArticlePubMedGoogle Scholar
- Canestrelli D, Cimmaruta R, Costantini V, Nascetti G: Genetic diversity and phylogeography of the Apennine yellow-bellied toad Bombina pachypus, with implications for conservation. Molecular Ecology. 2006, 15: 3741-3754. 10.1111/j.1365-294X.2006.03055.x.View ArticlePubMedGoogle Scholar
- Steward JR, Lister AM: Cryptic northern refugia and the origins of the modern biota. Trends in Ecology and Evolution. 2001, 16: 608-613. 10.1016/S0169-5347(01)02338-2.View ArticleGoogle Scholar
- de Lattin G: Beiträge zur Zoogeographie des Mittelmeergebietes. Verhandlungen der deutschen Zoologischen Gesellschaft, Kiel. 1949, 143-151.Google Scholar
- Nève G, Verlaque R: Genetic differentiation between and among refugia. Relict species - phylogeography and conservation biology. Edited by: Habel JC, Assmann T. 2010, Springer Heidelberg, 277-294.Google Scholar
- Habel JC, Augenstein B, Nève G, Rödder D, Assmann T: Population genetics and ecological niche modelling reveal high fragmentation and potential future extinction of the endangered relict butterfly Lycaena helle. Relict species - phylogeography and conservation biology. Edited by: Habel JC, Assmann T. 2010, Springer Heidelberg, 417-440.Google Scholar
- Schmitt T, Seitz A: Allozyme variation in Polyommatus coridon (Lepidoptera: Lycaenidae): identification of ice-age refugia and reconstruction of post-glacial expansion. Journal of Biogeography. 2001, 28: 1129-1136. 10.1046/j.1365-2699.2001.00621.x.View ArticleGoogle Scholar
- Gratton P, Konopinski MK, Sbordoni V: Pleistocene evolutionary history of the Clouded Apollo (Parnassius mnemosyne): genetic signatures of climate cycles and a 'time-dependent' mitochondrial substitution rate. Molecular Ecology. 2008, 17: 4248-4262. 10.1111/j.1365-294X.2008.03901.x.View ArticlePubMedGoogle Scholar
- Magri D: Patterns of post-glacial spread and the extent of glacial refugia of European beech (Fagus sylvatica). Journal of Biogeography. 2008, 35: 450-463. 10.1111/j.1365-2699.2007.01803.x.View ArticleGoogle Scholar
- Magri D, Vendramin GG, Comps B, Dupanloup I, Geburek T, Gomory D, Latalowa M, Litt T, Paule L, Roure JM, Tantau I, van der Knaap WO, Petit RJ, de Beaulieu JL: A new scenario for the Quaternary history of European beech populations: palaeobotanical evidence and genetic consequences. New Phytologist. 2006, 71: 199-221.View ArticleGoogle Scholar
- Hampe A, Petit RJ: Conserving biodiversity under climate change: the rear edge matters. Ecology Letters. 2005, 8: 461-467. 10.1111/j.1461-0248.2005.00739.x.View ArticlePubMedGoogle Scholar
- Habel JC, Dieker P, Schmitt T: Biogeographical connections between the Maghreb and the Mediterranean peninsulas of southern Europe. Biological Journal of the Linnean Society. 2009, 98: 693-703. 10.1111/j.1095-8312.2009.01300.x.View ArticleGoogle Scholar
- Nei M: Genetic distances between populations. The American Naturalist. 1972, 106: 283-291. 10.1086/282771.View ArticleGoogle Scholar
- Siegismund HR, Müller J: Genetic structure of Gammarus fossarum populations. Heredity. 1991, 66: 419-436. 10.1038/hdy.1991.52.View ArticleGoogle Scholar
- Excoffier L, Larval G, Schneider S: Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online. 2005, 1: 47-50.Google Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Biology and Evolution. 1987, 4: 406-425.PubMedGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny Inference Package) Ver. 3.5.c. Department of Genetics, University of Washington. 1993, Seattle, WashingtonGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-955.PubMedPubMed CentralGoogle Scholar
- Guisan A, Zimmermann N: Predictive habitat distribution models in ecology. Ecological Modelling. 2000, 135: 147-186. 10.1016/S0304-3800(00)00354-9.View ArticleGoogle Scholar
- Jeschke JM, Strayer DL: Usefulness of bioclimatic models for studying climate change and invasive species. Annals of the New York Academy of Sciences. 2008, 1134: 1-24. 10.1196/annals.1439.002.View ArticlePubMedGoogle Scholar
- Elith J, Leathwick JR: Species distribution models: Ecological explanation and prediction across space and time. Annual Reviews in Ecology, Evolution and Systematics. 2009, 40: 677-697. 10.1146/annurev.ecolsys.110308.120159.View ArticleGoogle Scholar
- Waltari E, Hijmans RJ, Peterson AT, Nyári AS, Perkins SL, Guralnick RP: Locating pleistocene refugia: Comparing phylogeographic and ecological niche model predictions. PLoS one. 2007, 7: 1-11.Google Scholar
- Hijmans RJ, Guarino L, Jarvis A, O'brien R, Mathur P, Bussink C, Cruz M, Barrantes I, Rojas E: DIVA-GIS version 5.2 manual. 2005Google Scholar
- Nogués-Bravo D: Predicting the past distributions of species climatic niches. Global Ecology and Biogeography. 2009, 18: 521-531. 10.1111/j.1466-8238.2009.00476.x.View ArticleGoogle Scholar
- Rödder D, Weinsheimer F, Lötters S: Molecules meet macroecology - combining species distribution modells and phylogeographic studies. Zootaxa. 2010, 2426: 54-60.Google Scholar
- Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A: Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology. 2005, 25: 1965-1978. 10.1002/joc.1276.View ArticleGoogle Scholar
- Peterson AT, Nyári ÁS: Ecological niche conservatism and pleistocene refugia in the thrush-like mourner, Shiffornis sp., in the neotropics. Evolution. 2007, 62: 173-183.PubMedGoogle Scholar
- Busby JR: Bioclim - a bioclimatic analysis and prediction system. Nature conservation: Cost effective biological surveys and data analysis. Edited by: CR Margules, MP Austin. 1991, CSIRO, Melbourne, 64-68.Google Scholar
- Beaumont LJ, Hughes L, Poulsen M: Predicting species distributions: Use of climatic parameters in bioclim and its impact on predictions of species' current and future distributions. Ecological Modelling. 2005, 186: 250-269.View ArticleGoogle Scholar
- Heikkinen RK, Luoto M, Araújo MB, Virkkala R, Thuiller W, Sykes MT: Methods and uncertainties in bioclimatic envelope modeling under climate change. Progress in Physical Geography. 2006, 30: 751-777. 10.1177/0309133306071957.View ArticleGoogle Scholar
- Rödder D, Lötters S: Niche shift versus niche conservatism? Climatic charactersistics of the native and invasive ranges of the mediterranean house gecko (Hemidactylus turcicus). Global Ecology and Biogeography. 2009, 18: 674-687. 10.1111/j.1466-8238.2009.00477.x.View ArticleGoogle Scholar
- Rödder D, Schmidtlein S, Veith M, Lötters S: Alien invasive slider turtle in unpredicted habitat: A matter of niche shift or of predictors studied?. PLoS one. 2009, 4: e7843-10.1371/journal.pone.0007843.View ArticlePubMedPubMed CentralGoogle Scholar
- Rödder D, Lötters S: Explanative power of variables used in species distribution modelling: An issue of general model transferability or niche shift in the invasive greenhouse frog (Eleutherodactylus planirostris). Naturwissenschaften. 2010, 97: 781-796. 10.1007/s00114-010-0694-7.View ArticlePubMedGoogle Scholar
- Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions. Ecological Modelling. 2006, 190: 231-259. 10.1016/j.ecolmodel.2005.03.026.View ArticleGoogle Scholar
- Phillips SJ, Dudík M: Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2010, 31: 161-175.View ArticleGoogle Scholar
- Pearce J, Ferrier S: An evaluation of alternative algorithms for fitting species distribution models using logistic regression. Ecological Modelling. 2000, 128: 128-147.View ArticleGoogle 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.