The evolutionary history of Antirrhinumin the Pyrenees inferred from phylogeographic analyses
© Liberal et al.; licensee BioMed Central Ltd. 2014
Received: 16 January 2014
Accepted: 17 June 2014
Published: 26 June 2014
The origin and colonisation history after the Quaternary ice ages remain largely unresolved for many plant lineages, mainly owing to a lack of fine-scale studies. Here, we present a molecular phylogeny and a phylogeographic analysis of Antirrhinum, an important model system in plant biology, in the Pyrenees range. Our goal was to reconstruct the evolutionary and colonisation history of four taxa endemic to this region (A. majus subsp. majus, A. majus. subsp. striatum, A. molle, and A. sempervirens) by using a dense sampling strategy, with a total of 452 individuals from 99 populations whose collective distribution spans nearly the entirety of the Pyrenees and adjacent mountains.
Phylogenetic and phylogeographic analyses of the sequences of two plastid (trnS-trnG and trnK-matK) regions revealed the following: (i) historical relationship between the Pyrenees and Iberia (but not with the Alps); (ii) the long persistence of populations in the Pyrenees, at least since the Late Pleistocene; (iii) three different colonisation histories for populations from the Western, Central, and Eastern Pyrenees; (iv) the deep phylogeographic separation of the eastern and western populations; and (v) the colonisation of southern France from the Eastern Pyrenees.
The present study underlines the enormous influence of the glacial history of the mountain ranges on the current configuration of intra- and inter-specific genetic diversity in Antirrhinum, as well as the importance of periglacial areas for the survival of species during glacial periods of the Quaternary.
The three mountain systems (the Pyrenees, Alps, and Balkan mountains) located to the north of the major southern European peninsulas (Iberia, Italy, and Balkans) have played an important role in determining the current distributions of European biota, and they harbour exceptionnally high levels of biodiversity [1–3]. Their high elevation greatly influenced the climate-induced range shifts of many species during Quaternary climatic cycles [4–7], and local speciation promoted by isolation of populations in long-standing refugia and potential contact zones between postglacial recolonizing genetic lineages have been invoked to account for such diversity [8, 9].
The Pyrenees are especially renowned as an important plant diversity hotspot because of their high floristic richness (around 3500 species and subspecies of vascular plants) [10, 11], including considerable endemicity (about 4% of the plant species are endemic) and numerous plant associations [12, 13]. Yet we still know relatively little about the processes that led to species accumulation within this region. Their east–west arrangement clearly acted as a latitudinal barrier for postglacial dispersal predominantly out of the Iberian Peninsula. Moreover, they served as an arena in which postglacial lineage contacts occurred . While the valleys in other European mountain ranges were only partially covered by ice even during the coldest periods of the ice age glaciations (e.g. Balkan Mountains and Carpathians), the Alps and the Pyrenees were almost entirely covered by large ice sheets [14–17]. Such similarities in their glacial history may have resulted in analogous patterns of genetic structure of the plant species that inhabit the two east–west mountain ranges. Unfortunately, in contrast to the high number of studies that have focused on the evolutionary history of plant lineages in the Alps [see , few investigations have addressed the origin and colonisation history of Pyrenean plant species in a phylogenetic context at the appropriate spatial scale. Most of the previous phylogeographic studies that included the Pyrenean range have relied upon limited sampling strategies (but see ). Thus, the major patterns related to the origins and colonisation of Pyrenean plant lineages after the Quaternary ice ages remain largely unresolved.
The present study used Antirrhinum species from the Pyrenees to reconstruct the colonisation history of mountain plants in this region. Monophyly of this genus, an important model system in plant biology primarily distributed in the western Mediterranean basin, has been previously proved . The Iberian Peninsula harbours the highest species diversity; however, a few Antirrhinum species are also found in other European areas, such as the Alps (A. latifolium) and Italy (A. siculum), or are distributed widely across the Mediterranean (A. tortuosum) despite the lack of any obvious long-distance dispersal syndrome. Historically, the number of Antirrhinum species has been the subject of taxonomic debate, resulting in different taxa circumscriptions. The difficulties in species delimitation, as well as the discrepancies between phylogenetic relationships and taxonomic classifications, have been interpreted as strong evidence for recent geographic speciation and extensive hybridisation since the Pliocene [21–24].
Four Antirrhinum taxa from two sections are endemic to the Pyrenees and pre-Pyrenees range. The pre-Pyrenees range is formed by a complex system of foothill ranges which stretches from the western side (influenced by the Atlantic regime) to the Mediterranean coast on the eastern end of the Pyrenees (commonly called sub-Pyrenees). This range mostly runs parallel to the Pyrenees along the Iberian side, but it also includes the Corbières Range, towards the eastern end on the French side, where the Pyrenees’s slopes descend rather abruptly. The four endemic taxa occur in different habitats throughout these mountain ranges, ranging from lowland areas (800 m) to high mountain environments (>2500 m). A. molle and A. sempervirens (sect. Kickxiella) inhabit limestone cliffs from medium to high altitudes, whereas A. majus subsp. majus and A. majus ssp. striatum (sect. Antirrhinum) occur mostly in lowland areas of the Eurosiberian woodlands (see Additional file 1 for further clarification of the different taxonomic treatments). Previous phylogeographic studies found that four of the eight main Antirrhinum lineages are distributed in northeast Iberia . Thus, the Pyrenees along with the pre-Pyrenees mountains are one of the areas with the highest numbers of genotypes and lineages. Recent studies have also been conducted to understand the cause and the maintenance of the parapatric distribution of two of these taxa (A. majus subsp. majus and A. majus subsp. striatum), as well as the ecological process that occurs in narrow contact zones [25–27]. Thus, it seems reasonable to assume that the present-day distribution of Antirrhinum in the Pyrenees is the result of historical processes that have taken place across the mountain range and adjacent mountains during Quaternary times. We evaluate this possibility with molecular phylogenetic and phylogeographic analyses.
Based on previous phylogenetic and phylogeographic studies on plants (see Additional file 2), four hypotheses (not mutually exclusive) can be proposed for explaining the origin of Antirrhinum lineages that are found in the Pyrenees: (i) a predominantly southern (Iberian) origin; (ii) a predominantly northern (central European) origin; (iii) a long-standing isolation of Pyrenean populations; and (iv) recurrent secondary contacts between Iberian and central European lineages.
In this study, we first attempt to infer the historical connections between the Pyrenees and other geographic areas. We also investigate the genetic structure within the species and assess the patterns of genetic diversity within the Pyrenees region. By combining the results of these analyses, we explore potential processes that could explain the evolutionary history of Antirrhinum lineages throughout the Pyrenees.
Features of the two sequenced cpDNA regions
Aligned length (bp)
Ungapped length range
Pairwise % identity
Mean % G + C content
Estimation of divergence times
Ancestral area reconstructions
DPA analysis showed uncertainty in the geographic origin of the Pyrenean lineages; since several areas outside the Pyrenees showed similar probabilities (see Figure 2). Nevertheless, stronger support was found for the geographic origin of more recent subclades, leading to a robust reconstruction of dispersal events within and around the Pyrenees during the Pleistocene. Bayes factor (BF) tests (Additional file 6) for significant non-zero rates revealed a historically close connection between the Pyrenees and the Iberian Peninsula (BF = 21.93). In contrast, South-western Alps populations (where A. latifolium occurs) were not linked to the Pyrenees. The Central Pyrenees were the most likely ancestral location of western haplotypes (BF = 36.74). A colonisation episode involving a few populations may have also occurred from Central Pyrenees to the Southern French basin, as indicated by a well-supported connection between these two areas (BF = 10.45). However, the colonisation of the Southern French basin seems to have occurred mainly from the Eastern Pyrenees as indicated by a considerably higher Bayes factor value (BF = 404.17). Interestingly, a less supported link was also observed between Eastern and Central Pyrenees (BF = 4.13). These results suggest: (i) no close-relationship between Pyrenean lineages and those distributed in the South-western Alps; (ii) an origin of Pyrenean lineages during the Pleistocene; and (iii) recent dispersal episodes within the Pyrenees from two main spreading centers (Eastern and Central Pyrenees).
When the number of areas was reduced to four –Pyrenees and surrounding areas, South-western Alps, Iberian Peninsula (excluding the Pyrenees) and rest of Mediterranean areas– the inferred historical biogeographic scenarios from the S-DIVA (Statistical Dispersal-Vicariance Analysis) and DEC (Dispersal-Extinction-Cladogenesis) analyses were mostly congruent, but not very informative about the origin of the Pyrenean lineages (Additional files 7 and 8 respectively). However, the results agreed with a sister relationship between the South-western Alps lineage and an Iberian clade, indicating: (i) a historical connection between two distant areas and (ii) no recent relationship between the Alps and Pyrenean lineages.
Haplotype data analysis
The statistical parsimony analysis of the Antirrhinum matrix found 11 haplotype lineages, of which seven included samples from NE Iberia. In addition, five haplotype lineages from the Pyrenees were formed by 18 haplotypes (haplotypes 2, 9-12 and 52 - 63, see Additional file 9) obtained from the 99 populations collected across the Pyrenees and adjacent areas. Some missing haplotypes (extinct or not found) and one loop were retrieved for the whole genus. Most of the haplotypes were grouped primarily together within lineages I and II of northern Iberia (Additional file 9). The five samples collected in the South-western Alps provided two new haplotypes (haplotypes 64 and 65) with a tip position, which formed part of lineage VI from southeastern Iberia.
Genetic structure and diversity
The geographical distribution of haplotypes and lineages is shown in Figure 3 (information per grids summarized in Additional file 10). The Eastern Pyrenees harbored only haplotypes of lineage E, which appeared to be primarily distributed in this area. The Western Pyrenees harbored exclusively the haplotype 4. The five lineages were distributed across the Central Pyrenees. The majority of 10×10 km quadrats (54) harbored only one haplotype, although some contained two (8 quadrats) and three (3 quadrats) haplotypes. Moreover, few quadrats had haplotypes of two (grid coordinates I14, I25, J21, J22, K19) and three (J17, F8) lineages, and all of them fell into the Central Pyrenees.
Number and frequency of each haplotype (%)
Genetic diversity in the three Pyrenean sections
Pairwise F ST -values
Conventional F-Statistics from haplotype frequencies
Results of AMOVA
Lineage E-rest of lineages
Tajima’s D and Fu’s FS test
Tajima’s D p-value
Origins of Antirrhinumin the Pyrenees
Regarding the geographic origin of the Pyrenean lineages, an important historical connection between the northeast Iberian Peninsula and the Pyrenees is well supported (Figure 2, Additional file 6). It is notable that none of the analyses identified evolutionary relationships between the lineages from the Pyrenees and those from the South-western Alps, which is a connection found in other plant species with similar distributions e.g. Erodium, , Primula, . The lack of a connection between these two areas is also supported by nrITS sequences (Liberal et al. unpublished, see additional file 12). This is a rather unexpected, but essential, result since the yellow-flowered populations that inhabit the Alps and the Pyrenees have recently been circumscribed into a single species (A. latifolium Mill.) by some taxonomists [32, 33]. The phylogenetic and dating analyses provide strong support for the divergence of Pyrenean sublineages in the Quaternary but poor support for sister-group relationships at the most basal nodes. The analyses presented here and those of Vargas et al.  support the first hypothesis herein posit that the Pyrenean lineages are the result of long-standing isolation and geographic differentiation from a widespread ancestor in the Pleistocene. However, additional nuclear data, such as single-copy nuclear markers, would be required to determine whether the unresolved basal relationships are attributable to hard or soft polytomies .
East–west geographical structure of the genetic diversity
The main genetic discontinuity identified by the analyses (SAMOVA, AMOVA, and FST) indicated the differentiation of the Lineage E from the remaining lineages. These two groups basically correspond to the division between the Eastern and Central-West Pyrenees (Figure 3). This cryptic transition area, which approximately follows the abrupt valley of the Segre River, is supported by the phylogeographic gap in the southeastern part of the Pyrenees that was already identified in previous studies of animals and plants e.g. [35–38]. A similar phytogeographic boundary was previously proposed based on species distributions in the Pyrenees [39, 40]. The high level of differentiation between the Eastern Pyrenean lineage (Lineage E) and the remaining lineages, suggests the isolation of eastern populations over more than one glacial cycle.
Interestingly, the eastern and central areas may have experienced different historical processes during the last glaciations. According to the results of the mismatch distribution analyses, the observed genetic structure in the Eastern Pyrenees is congruent with a range expansion process, whereas this pattern was not observed in the Central Pyrenees. By contrast, the strong evidence for the genetic clustering of Central Pyrenean samples, which was detected by the BAPS and SAMOVA analyses, as well as the significant genetic differentiation among the different lineages that occur in the area, support long-standing genetic isolation in different local refugia without an extensive interglacial (or postglacial) range expansion as already inferred for Pinus uncinata [according to . However, the homogeneously distributed genetic diversity indicates that moderate gene flow occurred among Central Pyrenean refugia during interglacial periods. Thus, two different historical processes were inferred: (i) limited extinction/expansion during glacial cycles in the central area, and (ii) population expansion (probably pre-dated by a demographic bottleneck based on the reduced number of lineages and haplotypes) in the eastern area, may explain the observed differences in the genetic structure and diversity of both regions. A similar interpretation was proposed for Ramonda myconi to explain genetic differences between populations from Central and Eastern Pyrenees  where significant evidence was presented to support the existence of an important refuge near the central area. Further support for a large refuge in Central Pyrenees was also found by palynological data [42, 43], the patterns of phylogeographic concordance among different species , and the high diversity of habitats, species richness, and endemic species [63 of the 160 endemic plant species of the Pyrenean flora are exclusively from the Central Pyrenees; . In agreement with this and other studies e.g. , a scenario of multiple refugia in mountains adjacent to the Central Pyrenees (i.e., the southern mountains of the pre-Pyrenees) is gaining more support (see Figure 3).
Similar to the pattern found in Antirrhinum, Ramonda myconi also exhibits considerable differentiation but reduced genetic diversity in the eastern population . Fossil pollen records and molecular studies suggest the presence of smaller refugia in the southeastern Pyrenees, along the Mediterranean coast [47–49]. Recent studies indicate that climatic fluctuations during the last glacial period affected the Eastern Pyrenees (Mediterranean coast) more strongly than the Central Pyrenees, at least during the last 50,000 years. Overall, more pronounced oscillations between arid and humid periods prevailed in the easternmost regions [50, 51]. Thus, severe reductions in population sizes and the isolation of peripheral populations at the edge of the distribution might have led to genetic drift and local adaptations via selective forces see (see [52–55]). Future comparative phylogeographic studies and coalescent demographic analyses could help to elucidate the historical processes responsible for the phylogeographic and floristic configuration of the Pyrenees.
Migration routes across the Pyrenees
The colonisation of the Western Pyrenees occurred from the Central Pyrenees. Indeed, the western populations have a low level of genetic diversity, which is shared by central populations. Such gradual genetic impoverishment observed from the Central to the Western Pyrenees is interpreted as a recent colonisation, which originated in the most central parts of the mountain range see [4, 5, 56, 57]. The results of the DPA analysis also indicate that the postglacial colonisation of the southern French basin occurred primarily from the Eastern Pyrenees (Figure 3). The geographic distribution of the haplotypes and the mismatch distribution analyses (Figure 4) suggest a gradual colonisation of southern France from northeastern Iberia, where the edge of the glacial distribution range may have retreated early see .
Similarities between the Pyrenees and the Alps
The causes of the separation of Antirrhinum lineages and plant species by the abrupt valley in the Pyrenees (Segre River) are unclear. Similar west–east vicariant lineages abutting along the Etsch Valley and the Brenner Pass in the Alps have been reported [18, 59–61]. Indeed, this discontinuity accounts for the most important phytogeographic and biogeographic boundary in the Alps e.g. [62–64]. Similar genetic and floristic discontinuities are found in the Alps and the Pyrenees, which suggest that plant extinctions due to Holocene glaciers and further recolonisation events were not distributed homogeneously across each of the two mountain ranges. The major effect of the glacial history of mountain valleys in the current configuration of intra- and inter-specific genetic diversity requires further investigation.
The present study also underlines the importance of periglacial areas for the survival of species during glacial periods of the Quaternary. The pre-Pyrenees range, is recognised as an important glacial refuge where many animal and plant species survived during cold episodes [30, 38, 65, 66]. Previous studies focused on the Alps have also recognised the importance of marginal and peripheral areas at lower altitudes with warmer climates as refugia for mountain species, where their present-day distribution were almost entirely glaciated reviewed in .
Our results do not support connections between the Antirrhinum populations of the Alps and those of the Pyrenees. Instead, the analyses agree with ancient connection with Iberia followed by persistence of Antirrhinum populations in the Pyrenees, at least since Late Pleistocene times. The three Pyrenean sections (Western, Eastern, and Central) have different colonisation histories. The Western Pyrenean populations are closely related to the central populations, which is here interpreted as evidence for a recent westward colonisation (see Figure 2). It is likely that the Central pre-Pyrenees played an important role in the maintenance of the genetic diversity of Antirrhinum by providing suitable conditions for the establishment of probably separate refugial populations during glacial periods. A significant influx of migrants from adjacent refugia (Central pre-Pyrenees) during the interglacial and postglacial periods would account for the genetic diversity found in the Central Pyrenees. However, the current genetic structure of Antirrhinum in the Eastern Pyrenees appears to have been shaped by range expansion, which was probably preceded by range contraction at the southernmost edge of this area, from where little gene flow may have occurred into the Central Pyrenees over several glacial cycles. According to this, the Eastern pre-Pyrenees also appear to have been a refugial area, although less important than the central pre-Pyrenees, from where Antirrhinum colonisation southern France.
Maternally inherited plastid markers often display stronger differentiation between populations than biparental markers due to their smaller effective size and lack of recombination [56, 67–69]. Because of this, reconstruction of migration histories by means of seed dispersal have been traditionally based on organelle markers once maternal inheritance has been documented [70, 71]. To this end, artificial crossing experiments under greenhouse conditions were conducted between four Antirrhinum taxa (A. mollissimum, A. majus subsp. majus, A. controversum and A. charidemi) (Additional file 13).
Sampling strategy and DNA sequencing
We collected a total of 452 individuals from 99 populations of A. molle, A. sempervirens, A. majus subsp. majus and A. majus subsp. striatum from the Pyrenees and adjacent mountains. In addition, a total of five individuals from five populations of A. latifolium distributed across the Province/southern Maritime Alps (Additional file 3) were analyzed given that this species has historically been considered closely related to A. majus subsp. striatum [32, 33]. The number of populations sampled per species depended on their distribution range: A. molle (3), A. sempervirens (9), A. majus subsp. majus (61) and A. majus subsp. striatum (26) respectively. Particular effort was made to achieve a sampling as complete as possible for the entire mountain range (see Figure 3). All individuals were collected in the field and dried in silica gel. Total genomic DNA was extracted using Dneasy Plant Mini Kit (QIAGEN Inc., California).
Two sequencing strategies were adopted. First, trnK-matK and trnS-trnG intergenic spacer sequences were obtained for one individual from each population, and added on to the sequence data matrix from Vargas et al. . This dataset (hereafter Antirrhinum matrix) was used to infer the colonisation history of Pyrenean lineages by phylogenetic and phylogeographic analyses. Second, trnS-trnG intergenic spacer sequences were obtained for the 452 individuals (1-6 individuals per population) from the Pyrenees. This sequence matrix (hereafter Pyrenees matrix) was used to investigate underlying spatial genetic structure and genetic diversity across the Pyrenees.
The trnK-matK and trnS-trnG regions were amplified as in Vargas et al. . PCR products were outsourced for sequencing to a contract sequencing facility (Macrogen, Seoul, South Korea) on an ABI Prism® 3730xi DNA sequencer, using the same primer sets as for PCR. Resulting sequence data were assembled and edited using Geneious Pro v5 . Both datasets were aligned with the MAFFT v.6.814b alignment tool , as implemented in Geneious Pro with default parameters. Further adjustments were made by visual inspection.
Phylogenetic analyses (Antirrhinum matrix)
Bayesian phylogenetic analyses were performed on the two separate matrices (trnS-trnG/trnK-matK) to examine plastid gene tree congruence. Gambelia speciosa and Misopates orontium were selected as the outgroup based on previous phylogenetic evidence . In addition, phylogenetic analyses for 190 trnS-trnG/trnK-matK concatenated sequences were conducted using Bayesian inference (BI), as implemented in MrBayes v3.1.2 , maximum likelihood (ML), as implemented in PhyML 3.0 , and maximum parsimony (MP), as implemented in TNT 1.1  (see Additional file 14 for details).
Estimation of divergence times (Antirrhinum matrix)
A relaxed molecular-clock approach implemented in BEAST v.1.6.2  was used to estimate the divergence times between Antirrhinum Pyrenean lineages. This software estimates the phylogenetic tree and node ages simultaneously. Since no reliable fossils of Antirrhinum are known to date, we implemented a secondary basal calibration following  but with some modifications. The divergence time between the basalmost Antirrhinum lineages (crown node) was modelled as a normal distribution with mean = 4.5 Ma and standard deviation = 1.8, on the basis of a previous analysis of ndhF sequences of Antirrhineae that included two most distant Antirrhinum species (A. meonantum and A. majus subsp. majus, following results of phylogenetic analyses). The latter analysis incorporated a calibration of 74 Ma for the divergence time between Oleaceae and Antirrhineae , and minimum stem-age constraints for Lamiales families and tribes based on five fossils see  for details on fossils].
Ancestral area reconstructions (Antirrhinum matrix)
A discrete phylogeographic analysis (DPA) implemented in BEAST  was performed to assess the probability distribution of the geographic locations in each node. A total of 14 discrete areas were delimited: (i) the four Iberian quadrants (northeastern Iberia, NE; northwestern Iberia, NW; southeastern Iberia, SE; southwestern Iberia, SW), as divided by the geographical coordinates 40ºN/5ºW see ; (ii) Eastern, Central and Western Pyrenees, as the three recognized biogeographic regions within the Pyrenees (see below); (iii) the other two northern areas sampled nearby the Pyrenees (Southern French basin and South-western Alps); and (iv) the remaining five areas sampled across Mediterranean basin (northern Africa, Sicily, Sardinia, Italian Peninsula and Turkey). Statistical significance for the rates of the dispersal events was assessed via Bayes factor (BF) test as described by Lemey et al. ; see Additional file 14 for details].
Additional ancestral range reconstructions were conducted using the Bayesian time-calibrated molecular phylogeny with the aim to discriminate between northern and southern origin of Pyrenean lineages. For this purpose, only four areas were delimited (i) Iberian Peninsula, (ii) Pyrenees and adjacent areas, (iii) South-western Alps, and (iv) samples from the Mediterranean basin (excluding Iberia). Ancestors were allowed to be present in all of them. Distribution ranges of sequences (haplotypes) instead of species was used see . Two alternative reconstruction methods were used: (a) Statistical Dispersal-Vicariance Analysis (S-DIVA) implemented in the program RASP 1.1 , and (b) dispersal-extinction-cladogenesis analysis (DEC) implemented in the software package Lagrange v2.0.1 ; see Additional file 14 for details].
Haplotype network analyses (Antirrhinum and Pyrenees matrices)
Haplotype networks using the Antirrhinum and Pyrenees matrices (hereafter Antirrhinum and Pyrenees networks, respectively) were constructed under the statistical parsimony framework in TCS 1.2.1 . The maximum number of differences among haplotypes, as a result of single substitutions, was calculated with 95% confidence limits. Indels of varying lengths resulting from the expansion/contraction of polyA and polyT were excluded due to their high levels of homoplasy . One indel of 9 bp was found in several samples of the population Villefranche-de-Conflent of A. majus. subsp. striatum for the trnS-trnG plastid region (Additional file 3) and was re-coded as a single character and treated as a fifth character state in. For the Pyrenees network, one sequence of A. braun-blanquetii (EU673480) was used as the outgroup, since this sequence appears to retain the ancestral haplotype in the Antirrhinum network analysis.
Genetic diversity and geographic structure (Pyrenees matrix)
An analysis of genetic diversity was carried out across the three recognized biogeographic regions in which the Pyrenees range is divided (Eastern, Central and Western Pyrenees) (see Figure 3a). The boundaries of these three biogeographic areas, although with slight differences, have been traditionally established by both geologists [86, 87] and phytogeographers [39, 40, 88–91] on the basis of geologic, climatic and floristic data. The boundaries of these areas are approximately defined by the Gallego River, between Central and Western Pyrenees, and the Segre River, between Central and Eastern Pyrenees. Haplotype frequencies and molecular diversity indices for each biogeographic area were calculated using DnaSP v5 . In addition, to identify potential hotspots of genetic diversity across the Pyrenees, individuals were geographically grouped by means of a 10×10 km grid. Charts representing haplotype frequencies were constructed for each grid cell, which was named by a generic letter–number code (Figure 3).
To identify genetic subdivisions among the Eastern, Central and Western Pyrenees, we performed an analysis of molecular variance (AMOVA) . Pairwise FST statistics were also calculated to estimate genetic distances. Both analyses were performed by using ARLEQUIN . Additionally, an AMOVA was performed in order to assess the partitioning of variance between the Lineage E (primarily distributed in the eastern part of the Pyrenees) and the rest of lineages (see below). Furthermore, to evaluate the optimal grouping of the sampled sites without a priori assumptions, a Bayesian analysis of population structure, implemented in the BAPS software, version 5.3 [95, 96], and a spatial analysis of molecular variance (SAMOVA) implemented in the software package SAMOVA 1.0  were also performed (see Additional file 13 for details on analyses).
Demographic history (Pyrenees matrix)
Past range expansions were assessed for the two partitioned population groups (lineage E and remaining populations, see below) using two methods. First, we tested for range expansion using Tajima’s D  and Fu’s Fs statistics  which, assuming neutrality, are expected to be significantly negative under population expansion. Significance for each test was assessed by generating null distributions from 10000 simulations. Second, we also conducted a mismatch analysis of cpDNA sequence differences, comparing the observed frequencies of pairwise differences of haplotypes with those expected under a sudden expansion (pure demographic expansion) model [29, 100]. Under expansion, the distribution of the observed differences is expected to be unimodal, whereas, under population contraction or genetic subdivision, a multimodal distribution is expected. Statistical significance was determined by 10000 bootstrap replicates. Goodness of fit was assessed by the sum of square deviations (SSD) and the Harpending’s raggedness index  between the observed and the expected mismatch. These analyses were conducted using ARLEQUIN 3.5. .
We would like to thank to Jose L. Blanco Pastor and Mario Fernández-Mazuecos for helping in the analysis and the manuscript writing; Emilio Cano, for laboratory assistance and Jose Luis Benito and Virginia Valcárcel for help with the sample collection. We would like to thank Betina Perales for assistance with map creation by using ArcGIS; and Dr. Jukka Corander and Dr. Andrew Rambaut for help with the BAPS and DPA analyses performance respectively. This research was supported by a project “Evolution of the personate flower” (CGL2009-10031) of the Spanish Ministry of Science and Technology and by the Synthesys program (CS and MB).
Genbank Accession numbers for DNA sequences: trnS-trnG (KJ911400-KJ911861) and trnK-matK (KJ911294-KJ911399).
- Varga ZS, Schmitt T: Types of oreal and oreotundral disjunctions in the western Palearctic. Biol J Linn Soc. 2008, 93 (2): 415-430.Google Scholar
- Barthlott W, Mutke J, Rafiqpoor D, Kier G, Kreft H: Global Centers Of Vascular Plant Diversity. Plant Taxonomy to Evolutionary Biology Nova Acta Leopoldina NF 92, Nr 342. 2005, 61-83.Google Scholar
- EEA: Biogeographical Regions, Europe 2011. 2011, Copenhagen: European Environment AgencyGoogle Scholar
- Hewitt G: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405 (6789): 907-913.PubMedGoogle Scholar
- Taberlet P, Fumagalli L, Wust-Saucy A-G, Cosson J-F: Comparative phylogeography and postglacial colonization routes in Europe. Mol Ecol. 1998, 7 (4): 453-464.PubMedGoogle Scholar
- Schmitt T: Biogeographical and evolutionary importance of the European high mountain systems. Front Zool. 2009, 6: 9-PubMedPubMed CentralGoogle Scholar
- Hewitt G: Genetic consequences of climatic oscillations in the Quaternary. Philos Trans R Soc Lond Ser B Biol Sci. 2004, 359 (1442): 183-195.Google Scholar
- Médail F, Diadema K: Glacial refugia influence plant diversity patterns in the Mediterranean Basin. J Biogeogr. 2009, 36 (7): 1333-1345.Google Scholar
- Hewitt GM: Post-glacial re-colonization of European biota. Biol J Linn Soc. 1999, 68 (1–2): 87-112.Google Scholar
- Dupias G: Végétation Des Pyrénées. Notice Détaillée De La Partie Pyrénéenne De La Carte De La Végétation De La France au 1/200.000e. 1985, Toulouse, FR: CNRSGoogle Scholar
- Gómez D, Sesé JA, Villar L: The vegetation of the Alpine Zone in the Pyrenees. Alpine Biodiversity In Europe. Edited by: Nagy L, Grabherr G, Koörner C, Thompson DBA. 2003, Berlin, DE: Springer, 85-92. [Ecological Studies no.167]Google Scholar
- Villar L, Dendaletche C: Pyrenees. France, Spain and Andorra. Centres Of Plant Diversity A Guide And Strategy For Their Conservation. Edited by: Heywood VR, Davis SD, Hamilton AC. 1994, Cambridge: WWF and IUCN, 61-64.Google Scholar
- Puşcaş M, Choler P: A biogeographic delineation of the European Alpine System based on a cluster analysis of Carex curvula dominated grasslands. Flora Morphol Distrib Funct Ecol Plants. 2012, 207: 168-178.Google Scholar
- Charlesworth JK: The Quaternary Era: With Special Reference To Its Glaciation. 1957Google Scholar
- Frenzel B, Pecsi M, Velichko AA: Atlas of the Paleoclimates and paleoenvironments of the Northern Hemisphere. Late Pleistocene Holocene. Gustav Fisher, Budapest. 1992Google Scholar
- Vuia F: Studiul Reliefului Glaciar oei Periglaciar din România. Thesis raport, Babes-Bolyai University, Cluj-Napoca. 2005Google Scholar
- Messerli B: Die eiszeitliche und die gegenwärtige Vergletscherung im Mittelmeerraum. Geogr Helv. 1967, 22: 105-228.Google Scholar
- Schönswetter P, Stehlik I, Holderegger R, Tribsch A: Molecular evidence for glacial refugia of mountain plants in the European Alps. Mol Ecol. 2005, 14 (11): 3547-3555.PubMedGoogle Scholar
- Dzialuk A, Muchewicz E, Boratyński A, Montserrat JM, Boratyńska K, Burczyk J: Genetic variation of Pinus uncinata (Pinaceae) in the Pyrenees determined with cpSSR markers. Plant Syst Evol. 2009, 277 (3–4): 197-205.Google Scholar
- Vargas P, Rosselló JA, Oyama R, Güemes J: Molecular evidence for naturalness of genera in the tribe Antirrhineae (Scrophulariaceae) and three independent evolutionary lineages from the New World and the Old. Plant Syst Evol. 2004, 249 (3–4): 151-172.Google Scholar
- Jiménez JF, Sánchez-Gómez P, Güemes J, Rosselló JA: Phylogeny of snapdragon species (Antirrhinum; Scrophulariaceae) using non-coding cpDNA sequences. Isr J Plant Sci. 2005, 53 (1): 47-54.Google Scholar
- Vargas P, Carrió E, Guzmán B, Amat E, Güemes J: A geographical pattern of Antirrhinum (Scrophulariaceae) speciation since the Pliocene based on plastid and nuclear DNA polymorphisms. J Biogeogr. 2009, 36 (7): 1297-1312.Google Scholar
- Carrió E, Forrest AD, Güemes J, Vargas P: Evaluating species nonmonophyly as a trait affecting genetic diversity: A case study of three endangered species of Antirrhinum L. (Scrophulariaceae). Plant Syst Evol. 2010, 288 (1): 43-58.Google Scholar
- Wilson Y, Hudson A: The evolutionary history of Antirrhinum suggests that ancestral phenotype combinations survived repeated hybridizations. Plant J. 2011, 66 (6): 1032-1043.PubMedGoogle Scholar
- Andalo C, Cruzan MB, Cazettes C, Pujol B, Burrus M, Thébaud C: Post-pollination barriers do not explain the persistence of two distinct Antirrhinum subspecies with parapatric distribution. Plant Syst Evol. 2010, 286 (3): 223-234.Google Scholar
- Khimoun A, Cornuault J, Burrus M, Pujol B, Thébaud C, Andalo C: Ecology predicts parapatric distributions in two closely related Antirrhinum majus subspecies. Evol Ecol. 2013, 27 (1): 51-64.Google Scholar
- Tastard E, Ferdy JB, Burrus M, Thébaud C, Andalo C: Patterns of floral colour neighbourhood and their effects on female reproductive success in an Antirrhinum hybrid zone. J Evol Biol. 2012, 25 (2): 388-399.PubMedGoogle Scholar
- Rambaut A, Drummond AJ: Tracer 1.5.0. 2007, Available at: http://beast.bio.ed.ac.uk/TracerGoogle Scholar
- Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9 (3): 552-569.PubMedGoogle Scholar
- Alarcón M, Vargas P, Sáez L, Molero J, Aldasoro JJ: Genetic diversity of mountain plants: Two migration episodes of Mediterranean Erodium (Geraniaceae). Mol Phylogenet Evol. 2012, 63 (3): 866-876.PubMedGoogle Scholar
- Zhang LB, Comes HP, Kadereit JW: The temporal course of Quaternary diversification in the European high mountain endemic Primula sect. Auricula (Primulaceae). Int J Plant Sci. 2004, 165 (1): 191-207.Google Scholar
- Güemes J: Antirrhinum L. Flora Iberica. Edited by: Castroviejo S, Herrero A, Benedí C, Rico E, Güemes J. 2009, Madrid: CSIC, vol. 13Google Scholar
- Sutton DA: A revision of the tribe Antirrhineae. 1988, Oxford: Oxford University PressGoogle Scholar
- Walsh H, Kidd M, Moum T, Friesen V: Polytomies and the power of phylogenetic inference. Evolution. 1999, 53: 932-937.Google Scholar
- Magri D, Vendramin GG, Comps B, Dupanloup I, Geburek T, Gömöry D, Latałowa 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 Phytol. 2006, 171 (1): 199-221.PubMedGoogle Scholar
- Albre J, Gers C, Legal L: Molecular phylogeny of the Erebia tyndarus (Lepidoptera, Rhopalocera, Nymphalidae, Satyrinae) species group combining CoxII and ND5 mitochondrial genes: a case study of a recent radiation. Mol Phylogenet Evol. 2008, 47 (1): 196-210.PubMedGoogle Scholar
- Kropf M, Kadereit JW, Comes HP: Differential cycles of range contraction and expansion in European high mountain plants during the late quaternary: insights from Pritzelago alpina (L.) O. Kuntze (Brassicaceae). Mol Ecol. 2003, 12 (4): 931-949.PubMedGoogle Scholar
- Sánchez de Dios R, Benito-Garzón M, Sainz-Ollero H: Hybrid zones between two European oaks: a plant community approach. Plant Ecol. 2006, 187 (1): 109-125.Google Scholar
- Rivas Martinez S: Memoria Del Mapa De Series De Vegetaci6n De Espafia 1: 400000. 1987, Madrid: ICONAGoogle Scholar
- Vigo J, Ninot J: Los Pirineos. La Vegetación de España. Edited by: Peinado Lorea M, Rivas Martinez S. 1987, Alcalá de Henares: Col. Aula Abierta. Publ. Univ, 351-384.Google Scholar
- Dubreuil M, Riba M, Mayol M: Genetic structure and diversity in Ramonda myconi (Gesneriaceae): effects of historical climate change on a preglacial relict species. Am J Bot. 2008, 95 (5): 577-587.PubMedGoogle Scholar
- Jalut G, Amat AE, Riera I, Mora S, Fontugne M, Mook R, Bonnet L, Gauquelin T: Holocene climatic changes in the western mediterranean: installation of the mediterranean climate. C R Acad Sci II Sci Terre Planetes. 1997, 325 (5): 327-334.Google Scholar
- Valero-Garcés BL, González-Sampériz P, Delgado-Huertas A, Navas A, Machín J, Kelts K: Lateglacial and Late Holocene environmental and vegetational change in Salada Mediana, central Ebro Basin, Spain. Quat Int. 2000, 73–74: 29-46.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, NL: Dordrecht, Springer, 155-188.Google Scholar
- Villar L, Sesé JA, Ferrández JV: Atlas de la flora del Pirineo Aragonés. II. (Pyrolaceae-Orchidaceae. Síntesis). 1997, Huesca, Spain: Consejo de Protección de la Naturaleza de Aragón & Instituto de Estudios AltoaragonensesGoogle Scholar
- Segarra-Moragues JG, Palop-Esteban M, González-Candelas F, Catalán P: Nunatak survival vs. tabula rasa in the Central Pyrenees: a study on the endemic plant species Borderea pyrenaica (Dioscoreaceae). J Biogeogr. 2007, 34 (11): 1893-1906.Google Scholar
- Pèrez-Obiol R, Julià R: Climatic change on the iberian peninsula recorded in a 30,000-yr pollen record from lake banyoles. Quat Res. 1994, 41 (1): 91-98.Google Scholar
- Salvador L, Alía R, Agúndez D, Gil L: Genetic variation and migration pathways of maritime pine (Pinus pinaster Ait) in the Iberian peninsula. Theor Appl Genet. 2000, 100 (1): 89-95.Google Scholar
- Lumaret R, Mir C, Michaud H, Raynal V: Phylogeographical variation of chloroplast DNA in holm oak (Quercus ilex L.). Mol Ecol. 2002, 11 (11): 2327-2336.PubMedGoogle Scholar
- Sánchez Goñi M, Cacho I, Turon J, Guiot J, Sierro F, Peypouquet J, Grimalt J, Shackleton N: Synchroneity between marine and terrestrial responses to millennial scale climatic variability during the last glacial period in the Mediterranean region. Clim Dyn. 2002, 19 (1): 95-105.Google Scholar
- Moreno A, Cacho I, Canals M, Grimalt JO, Sánchez-Goñi MF, Shackleton N, Sierro FJ: Links between marine and atmospheric processes oscillating on a millennial time-scale. A multi-proxy study of the last 50,000 yr from the Alboran Sea (Western Mediterranean Sea). Quat Sci Rev. 2005, 24 (14–15): 1623-1636.Google Scholar
- García-Ramos G, Kirkpatrick M: Genetic models of adaptation and gene flow in peripheral populations. Evolution. 1997, 51 (1): 21-28.Google Scholar
- Cassel-Lundhagen A, Tammaru T, Windig JJ, Ryrholm N, Nylin S: Are peripheral populations special? Congruent patterns in two butterfly species. Ecography. 2009, 32 (4): 591-600.Google Scholar
- Lesica P, Allendorf FW: When are peripheral populations valuable for conservation?. Conserv Biol. 1995, 9 (4): 753-760.Google Scholar
- Losos JB, Glor RE: Phylogenetic comparative methods and the geography of speciation. Trends Ecol Evol. 2003, 18 (5): 220-227.Google Scholar
- Petit RJ, Brewer S, Bordács S, Burg K, Cheddadi R, Coart E, Cottrell J, Csaikl UM, van Dam B, Deans JD: Identification of refugia and post-glacial colonisation routes of European white oaks based on chloroplast DNA and fossil pollen evidence. For Ecol Manag. 2002, 156 (1): 49-74.Google Scholar
- Schmitt T: Molecular biogeography of Europe: pleistocene cycles and postglacial trends. Front Zool. 2007, 4 (11): 1-13.Google Scholar
- Templeton AR: Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history. Mol Ecol. 1998, 7 (4): 381-397.PubMedGoogle Scholar
- Ronikier M, Schneeweiss GM, Schönswetter P: The extreme disjunction between Beringia and Europe in Ranunculus glacialis s. l. (Ranunculaceae) does not coincide with the deepest genetic split - a story of the importance of temperate mountain ranges in arctic-alpine phylogeography. Mol Ecol. 2012, 21 (22): 5561-5578.PubMedGoogle Scholar
- Albach DC, Schönswetter P, Tribsch A: Comparative phylogeography of the Veronica alpina complex in Europe and North America. Mol Ecol. 2006, 15 (11): 3269-3286.PubMedGoogle Scholar
- Schönswetter P, Tribsch A, Barfuss M, Niklfeld H: Several Pleistocene refugia detected in the high alpine plant Phyteuma globulariifolium Sternb. & Hoppe (Campanulaceae) in the European Alps. Mol Ecol. 2002, 11 (12): 2637-2647.PubMedGoogle Scholar
- Kerner A: Die natürlichen Floren im Gelände der deutschen Alpen. Die Dtsch Alpen für Einheimische und Fremde Geschildert. 1871, 1: 126-189.Google Scholar
- Aeschimann D, Rasolofo N, Theurillat JP: Analysis of the flora of the Alps. 1: historical account and biodiversity. Analyse de la flore des Alpes 1: Historique et biodiversité. 2011, 66 (1): 27-55.Google Scholar
- Pampanini R: Essai Sur La Géographie Botanique Des Alpes Et En Particulier Des Alpes Sud-Orientales. Imprimerie Fragnière frères. 1903Google Scholar
- Theissinger K, Bálint M, Feldheim KA, Haase P, Johannesen J, Laube I, Pauls SU: Glacial survival and post-glacial recolonization of an arctic-alpine freshwater insect (Arcynopteryx dichroa, Plecoptera, Perlodidae) in Europe. J Biogeogr. 2013, 40 (2): 236-248.Google Scholar
- Schmitt T, Hewitt GM: Molecular biogeography of the arctic-alpine disjunct burnet moth species Zygaena exulans (Zygaenidae, Lepidoptera) in the Pyrenees and Alps. J Biogeogr. 2004, 31 (6): 885-893.Google Scholar
- Petit RJ, Kremer A, Wagner DB: Finite island model for organelle and nuclear genes in plants. Heredity. 1993, 71 (6): 630-641.Google Scholar
- Ennos R, Sinclair W, Hu X, Langdon A: Using organelle markers to elucidate the history, ecology and evolution of plant populations. Systematics association special volume. 1999, 57: 1-19.Google Scholar
- Soltis DE, Gitzendanner MA, Strenge DD, Soltis PS: Chloroplast DNA intraspecific phylogeography of plants from the Pacific Northwest of North America. Plant Syst Evol. 1997, 206 (1): 353-373.Google Scholar
- Hansen AK, Escobar LK, Gilbert LE, Jansen RK: Paternal, maternal, and biparental inheritance of the chloroplast genome in Passiflora (Passifloraceae): implications for phylogenetic studies. Am J Bot. 2007, 94 (1): 42-46.PubMedGoogle Scholar
- Azhagiri AK, Maliga P: Exceptional paternal inheritance of plastids in Arabidopsis suggests that low-frequency leakage of plastids via pollen may be universal in plants. Plant J. 2007, 52 (5): 817-823.PubMedGoogle Scholar
- Drummond A, Ashton B, Buxton S, Cheung M, Cooper A, Heled J, Kearse M, Moir R, Stones-Havas S, Sturrock S: Geneious v5. 1. 2010, http://www.geneious.com,Google Scholar
- Katoh K, Misawa K, Kuma KI, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30 (14): 3059-3066.PubMedPubMed CentralGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574.PubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704.PubMedGoogle Scholar
- Goloboff PA, Farris JS, Nixon KC: TNT, a free program for phylogenetic analysis. Cladistics. 2008, 24 (5): 774-786.Google Scholar
- Drummond AJ, Rambaut A: BEAST: bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7 (1): 214-PubMedPubMed CentralGoogle Scholar
- Fernández-Mazuecos M, Vargas P: Congruence between distribution modelling and phylogeographical analyses reveals Quaternary survival of a toadflax species (Linaria elegans) in oceanic climate areas of a mountain ring range. New Phytol. 2013, 198: 1274-1289.PubMedGoogle Scholar
- Bell CD, Soltis DE, Soltis PS: The age and diversification of the angiosperms re-revisited. Am J Bot. 2010, 97 (8): 1296-1303.PubMedGoogle Scholar
- Fernández-Mazuecos M, Vargas P: Historical isolation versus recent long-distance connections between Europe and Africa in bifid toadflaxes (Linaria sect. Versicolores). PLoS One. 2011, 6 (7): e22234-PubMedPubMed CentralGoogle Scholar
- Lemey P, Rambaut A, Drummond AJ, Suchard MA: Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009, 5 (9): e1000520-PubMedPubMed CentralGoogle Scholar
- Yu Y, Harris A, He X: Rasp (Reconstruct Ancestral State In Phylogenies) 2.0 Beta. 2011, http://mnh.scu.edu.cn/soft/blog/RASP,Google Scholar
- Ree RH, Smith SA: Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst Biol. 2008, 57 (1): 4-14.PubMedGoogle Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9 (10): 1657-1659.PubMedGoogle Scholar
- Kelchner SA: The evolution of non-coding chloroplast DNA and its application in plant systematics. Ann Mo Bot Gard. 2000, 87 (4): 482-498.Google Scholar
- Souquet P, Bilotte M, Canerot J, Debroas E, Peybernés B, Rey J: Nouvelle interprétation de la structure des Pyrénées. C R Acad Sci Paris. 1975, 281: 609-612.Google Scholar
- Souquet P, Mediavilla F: Nouvelle hypothèse sur la formation des Pyrénées. CR Acad Sci. 1976, 282: 2139-2142.Google Scholar
- Rivas Martínez S, Báscones Carretero JC, Díaz González TE, Fernández González F, Loidi Arregui J: Vegetación del Pirineo occidental y Navarra: VI Excursión Internacional de Fitosociología (AEFA). Itinera Geobotanica. 1991, 5: 5-456.Google Scholar
- Montserrat P: L’exploration floristique des Pyrénées occidentales. (The floristic exploration of West Pyrenees). Bol Soc Brot. 1974, 47: 227-241.Google Scholar
- Villar L: La vegetación del Pirineo occidental: estudio de Geobotánica Ecológica. Príncipe de Viana Suplemento de Ciencias. 1982, 263-434.Google Scholar
- Izard M: Le Climat. La Végétation des Pyrénees. Edited by: Dupias G. 1985, Paris: Editions du CNRS, 17-36.Google Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25 (11): 1451-1452.PubMedGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics. 1992, 131 (2): 479-491.PubMedPubMed CentralGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005, 1: 47-Google Scholar
- Corander J, Marttinen P, Sirén J, Tang J: Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations. BMC bioinformatics. 2008, 9: 539-PubMedPubMed CentralGoogle Scholar
- Corander J, Tang J: Bayesian analysis of population structure based on linked molecular information. Math Biosci. 2007, 205 (1): 19-31.PubMedGoogle Scholar
- Dupanloup I, Schneider S, Excoffier L: A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002, 11 (12): 2571-2581.PubMedGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.PubMedPubMed CentralGoogle Scholar
- Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147 (2): 915-925.PubMedPubMed CentralGoogle Scholar
- Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129 (2): 555-562.PubMedPubMed CentralGoogle Scholar
- Harpending HC: Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Hum Biol. 1994, 66 (4): 591-600.PubMedGoogle Scholar
- Excoffier L, Lischer HEL: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010, 10 (3): 564-567.PubMedGoogle 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 credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.