Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression
© Milá et al.; licensee BioMed Central Ltd. 2013
Received: 4 April 2013
Accepted: 6 September 2013
Published: 10 September 2013
The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure.
The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France.
The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.
KeywordsCytonuclear incongruence Gene flow Phylogeography Secondary contact Speciation Vicariance
The geographical distribution of intraspecific lineages and the pattern of gene flow among them provides valuable insight into the process of lineage divergence and speciation . Of particular interest are phylogeographic studies at large spatial scales that encompass major geographic barriers to gene flow as well as areas where divergent lineages come into secondary contact. Patterns of introgression at these contact zones can shed light on the degree of reproductive isolation among lineages and the historical factors driving divergence . In Europe, glacial dynamics over the last several million years had a major impact on the phylogeography of animal and plant lineages, which underwent repeated cycles of rapid expansions from southern refugia as species recolonized northern latitudes following glacial maxima [3–5]. These range expansions gave rise to a series of suture zones across the continent as lineages that diverged during glacial periods came into secondary contact during interglacials, providing us with unique natural experimental areas in which to study the process of lineage divergence and the evolution of reproductive isolation .
One of these regions of secondary contact is the Ibero-Pyrenean suture zone, an important biogeographic region where both inter and intra-specific lineages are known to come into contact in a number of taxa [3–6]. In this area of southwestern Europe, the Pyrenees stand as a major barrier separating central Europe from the Iberian peninsula, one of the major glacial refugia in southern Europe, and also a geologically and ecologically complex area that has suffered itself severe climatic changes, giving rise to a complex phylogeographic pattern that consists of “refugia within refugia” . Here we examine the phylogeography and patterns of gene flow among lineages of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Ibero-Pyrenean region. The species has a broad Eurasian distribution composed largely of viviparous lineages, yet individuals in this region belong to an oviparous lineage isolated from the nearest viviparous populations in the French Massif Central by a gap of unsuitable habitat . Recent studies of Z. vivipara in this region have documented a sharp contact zone between two divergent maternal lineages corresponding generally to the French and Spanish sides of the Pyrenees, and showing consistent results with mtDNA 16S and cytochrome b gene sequences [8, 9] and a maternally inherited marker on the female W chromosome . However, because of limited regional sampling in previous phylogeographic studies, and the lack of biparentally inherited markers examined in this system to date, our understanding of region-wide genetic structure, patterns of introgression among lineages and the relative importance of the Pyrenees and other mountain ranges as barriers to gene flow remain general and tentative.
Inferring the demographic and evolutionary history of lineages at various spatial and temporal scales often requires the use of molecular markers from genomic regions with different modes of inheritance and rates of evolution . Mitochondrial DNA sequence markers have proven useful for studying intraspecific genetic variation because of their relatively high mutation rate, small effective population size, short coalescence times, and matrilineal mode of inheritance [1, 12, 13]. However, mt-DNA genomes are susceptible to the effect of historical factors such as population expansions [3, 14], introgressive hybridization [15, 16], incomplete lineage sorting  and even selection [18, 19], all of which can effectively mask real patterns of gene flow among populations and lineages, potentially confounding inferences of evolutionary and demographic events . Complementing mtDNA datasets with biparentally inherited nuclear DNA (nDNA) marker data is often essential to obtain a more complete picture of lineage history [11, 21].
Here we undertake a multilocus phylogeographic survey across the region using both mtDNA sequences and genome-wide AFLP loci and an improved geographic sampling in order to document the presence and distribution of major genetic lineages and test the role of the Pyrenees and the Cantabrian Mountains in driving lineage divergence in Z. vivipara. We also use this large-scale phylogeographic context to assess patterns of introgression and test hypotheses related to the history of contact zones.
Mitochondrial DNA structure and diversity
Sampling localities, sample sizes and mtDNA haplotype frequencies
GH(1), IF(2), JH(1), KH(1)
FH(1), GH(1), HH(11)
GH(5), AA(2), BC(1), CC(2)
AA(6), AB(3), EA(2)
Col de Bergout
Plateau de Bious
Gave du Baralet
Col de Bouesou
AA(4), TL(1), WL(5)
Lac de Fabrege
AA(5), SL(1), WL(3)
Col de Marie-Blanque
Moura de Montrol
Col de l’Ours
Col des Palomieres
Genetic divergence among Zootoca vivipara mtDNA lineages
NW Spain (yellow)
NC Spain (red)
NE Spain (blue)
S. France (green)
Sample sizes and diversity indices for the major mtDNA lineages of Zootoca vivipara
NW Spain (yellow)
0.816 ± 0.035
0.0033 ± 0.0019
NC Spain (red)
0.596 ± 0.067
0.0007 ± 0.0005
NE Spain (blue)
0.333 ± 0.065
0.0003 ± 0.0003
S. France (green)
0.868 ± 0.028
0.0022 ± 0.0013
AFLP structure and diversity
Genetic diversity of AFLP loci in four major Zootoca vivipara genetic clusters
SE (H e )
NW Spain (yellow)
NC Spain (red)
NE Spain (blue)
S. France (green)
Genetic differentiation among Zootoca vivipara lineages based on AFLP data
NW Spain (yellow)
NC Spain (red)
NE Spain (blue)
S. France (green)
None of the 34 loci showed significant departures from neutrality in the BayScan 2.0 analysis across all comparisons, although one of them (locus “15c194”) was identified as an outlier in the K = 2 comparison, and another one (locus “19c84”) in the K = 3 and K = 4 comparisons (Additional file 1: Figure S2). Excluding these two loci from the dataset yielded results for the Structure analysis and the PCoA analysis in Genalex that were indistinguishable from those using all loci. This is due to the fact that most markers in these comparisons had FST values above 0.20 (mean across-loci FST-K2 = 0.27, SD = 0.01; Mean FST-K3 = 0.25, SD = 0.01; Mean FST-K4 = 0.23, SD = 0.02) and thus most markers contributing to the pattern of genetic structure were selectively neutral.
Introgression at contact zones
Using the geographic distribution of the main four mtDNA clades and the comparable AFLP-defined genetic units corresponding to K = 4, we detected localities with mixed haplotypes (for mtDNA) or assignment probability values (for ALFPs) in two main areas within the region, one corresponding to the contact zone between blue and green lineages across the Pyrenees, and one between blue and red lineages in Navarre, just east of the Basque Country (Figure 2).
At a second mtDNA contact zone detected between the red and blue clades at the Ibañeta site in the Spanish side of the western Pyrenees (Figure 2A), the AFLP data revealed that all individuals had a higher probability of assignment to the blue cluster (Figure 2B). A more detailed examination of the 10 individuals genotyped from this site shows that six of them had >90% probability of assignment to the blue cluster and less than 10% to the red cluster, and the remaining four individuals had 86-48% probability of assignment to the blue cluster and 10-33% to the red cluster (Figure 6B).
A final discrepancy between the two datasets pertains to the assignment of the Beret site (BER), near the border in the central Pyrenees, where the six individuals sampled carry southern France mtDNA haplotypes (green clade), yet their AFLP profiles cluster with the NE Spain group (blue cluster) (Figure 2), a situation similar to that of many southern France individuals.
Phylogeography and mtDNA lineage history
Our mtDNA dataset revealed four major clades and a remarkable degree of geographic structuring among Zootoca vivipara lineages. The high variation found in the ND2 gene relative to cyt-b, together with a more thorough geographic sampling allowed us to identify several distinct and geographically segregated clades in northern Spain, which help explain their historical interactions with the divergent clade found in southern France and clarify the phylogeographic history of the region.
In line with previous phylogeographic studies on the species, our study shows that populations of Z. vivipara at the Ibero-Pyrenean zone do not show an actual “suture zone” between Iberian and northern European lineages [8, 24], in sharp contrast to patterns found in other species, like the grasshopper Chorthippus parallelus or the warbler Phylloscopus collybita. We also find no evidence of Iberian clades expanding northward to colonize northern European latitudes, as documented for toads , shrews , hedgehogs , or bears . Instead, the oviparous clades of Z. vivipara on opposite sides of the Pyrenees form a monophyletic clade with respect to all other European lineages, indicating a relatively long time in isolation in this region. However, our mtDNA data do support a major role for the Pyrenees in isolating populations, and the two main divergent lineages (depicted by node “1” in Figure 1A) appear to have diverged on either side of this major barrier within the last one million years, since around the mid Pleistocene. The Spanish clade further differentiated within the last 500 K years into at least three subclades corresponding to separate mountainous regions, namely the western and eastern sections of the Cantabrian Mountains, and the south-facing slope of the Pyrenees. Further structure is also indicated by the small clade in the western-most end of the range, formed by two haplotypes (XJ and X1J) sampled in the Serra do Xistral, Galicia, although further sampling will be necessary to confirm this pattern.
Overall, the general pattern of allopatric differentiation in isolated highlands is consistent with the species ecology, as in the southern part of its distribution the common lizard is found mainly in humid bogs and associated habitats [31, 32], which in Spain are often found in montane areas. Other lizards associated with humid areas, such as members of the genus Iberolacerta, show a deep phylogenetic break between the NW Spain form (I. monticola) and the Pyrenean taxa (I. bonnali, I. aranica and I. aurelioi), suggesting similar patterns of isolation as in Z. vivipara albeit at longer time scales [33, 34]. In contrast, species that do not specialize on humid habitats have shown a contrasting lack of structure across similar regions. For instance, populations of Lacerta schreiberi along northern Spain are genetically uniform and expanded from a refuge in central Portugal since the last glacial maximum , a pattern also shown by Podarcis bocagei and the salamander Chioglossa lusitanica.
The pattern of distinct Z. vivipara lineages associated with different mountain areas along the northern Iberian peninsula is consistent with a pattern of “refugia-within-refugia” documented for a number of other organisms [7, 38, 39], and is consistent with a climatically dynamic history which combined with a complex topography has given rise to allopatrically differentiated populations and lineages. In this context, Z. vivipara lineages are younger than those reported in other taxa, such as Psammodromus hispanicus, Alytes obstetricans, or Lissotriton boscai, and more similar in age to those in Alytes cisternasii.
AFLP structure: contrasting patterns with mtDNA
The pattern of spatial variation in genome-wide AFLP markers revealed both important similarities and differences compared to that of mtDNA sequence. Both datasets revealed the presence of four genetic units that are geographically concordant, yet relative divergence among groups differs in both datasets. AFLP variation forms two main genetic groups that separate NW Spain from the rest, while the two main mtDNA clades separate Spain from France. This difference suggests that at a large geographic scale different neutral and/or selective factors have shaped diversity and structure in both types of markers, yet identifying the specific cause is beyond the scope of our dataset. Potential causes range from neutral factors such as random sorting of loci or chromosomal rearrangements, to the role of natural selection. Even though the AFLP loci we used appear to be largely neutral, these markers are likely to be spread randomly across the nuclear genome, and are therefore more likely to be associated with regions under selection than mtDNA markers. If this were the case, we would expect to find a stronger association between traits under selection and AFLP structure. Available phenotypic data on Z. vivipara in the region provides some support for this hypothesis, and analysis of morphological traits has revealed marked phenotypic differentiation of Cantabrian populations with respect to those in the Basque Country and Pyrenees . This pattern is driven largely by the shape of the head-shield scales, a trait with unknown fitness value, although other traits related to coloration known to be important in related lacertids [44, 45] are yet to be analysed in Z. vivipara. A closer association between traits under selection and nuclear instead of mitochondrial neutral makers has been documented in other studies. For example, in a contact zone between central European newts, Babik et al.  document a tighter correlation between nuclear DNA and morphological hybrid indices than between mtDNA and morphology. Also, a recent study in the Ensatina salamander complex of North America shows that divergence in neutral nuclear markers is better correlated with reproductive isolation than mtDNA . Additional ecological and phenotypic data will be necessary to properly test the hypothesis that neutral nuclear markers show a better association to adaptive variation than mitochondrial markers.
Contact zone dynamics and patterns of introgression
Our phylogeographic results reveal the presence of two independent contact zones between lineages: one in southern France across the high Pyrenees (between the blue and green clades in Figure 2) and one in the foothills of the southwestern Pyrenees (between the red and blue clades). An additional contact zone in the Basque lowlands between the red and green clades was not detected here due to limited sampling in that area, but a population with haplotypes from both clades has been reported by Heulin et al. (2011) at Moura de Montrol sites. Patterns of introgression in this area have been studied in detail among lineages of Corthippus grasshoppers  and Phylloscopus birds [26, 49], both forming suture zones among old lineages. Additional spatial sampling will be necessary to properly compare these suture zones to the secondary contact dynamics among the relatively young lineages of Z. vivipara. Another potential contact zone not detected in the present study corresponds to that found in north-central Spain, between the yellow and red lineages, and further sampling would be necessary to confirm its existence.
Relatively good sampling in the central Pyrenees allows us to infer the demographic and introgression history of lineages there. Different lines of evidence (phylogeny, haplotype frequencies and expansion test) support the hypothesis that the presence of NE Spanish (blue) haplotypes in southern France is due to a recent population expansion from the Spanish side of the Pyrenees, which led to the colonization of southern France across the high Pyrenees. The marked predominance of haplotype “AA” there suggests that this incursion took place recently, probably during or since the last glacial maximum, as indicated by time estimates from the mismatch distribution. The relative geographic congruence among mtDNA and AFLP genetic groups allows us to compare patterns of introgression of maternal and bi-parentally inherited loci. Genome-wide AFLP data revealed limited introgression of nDNA into France, with individuals of mixed origin being restricted to the upper Ossau river drainage and the Spain-France border, suggesting that the contact zone among lineages is located within 10 or 15 km of the border. In contrast, Spanish mtDNA haplotypes are found at least tens of kilometres further than nDNA. Interestingly, it appears that Spanish haplotypes expanded into France (mostly towards the W and NW), they were then trapped at river valleys towards the east (i.e. lower Aspe river and upper Ossau river valleys), giving rise to a sharp contact zone there. This apparent barrier to gene flow in mtDNA but not nDNA at river valleys could be due indirectly to Haldane’s rule , which predicts a lower fitness for hybrids of the heterogametic sex, females in the case of Z. vivipara. A more direct role of selection through differential survivorship of alternative haplogroups has been previously suggested by Heulin et al. (2011), although recent evidence from the same study area in Gabas-piet (Heulin, unpublished data) indicates that further replication across years will be necessary to obtain conclusive evidence. Several studies in diverse organisms have shown that mtDNA can introgress further than nDNA [20, 46, 51–54], and spread relatively quickly. An alternative explanation for the shifted positions of the mtDNA and nDNA clines is that following the initial expansion, the contact zone has moved back south, leaving behind Spanish mtDNA haplotypes in its wake [55, 56]. Distinguishing among these alternative scenarios will require further geographic and molecular sampling, including co-dominant nuclear markers.
Additional fine-scale sampling at contact zones among lineages will also be needed to understand introgression dynamics at other areas, such as the Beret site, where all individuals carry French mtDNA and Spanish nDNA, or the contact zone in the western Pyrenees between the blue and red clades detected at the Ibañeta site, where individuals show blue central-Pyrenean nDNA, yet half of them carry red mtDNA haplotypes. Patterns at these individual sites are similar to those found at some sites in southern France, and a more thorough geographic context provided by finer-scale sampling will be needed to help determine whether they belong to active contact zones with ongoing introgression or instead represent populations left in the genetic wake of a contact zone that shifted away.
Our results show that lineages of Z. vivipara in northern Iberia and the Ibero-Pyrenean suture zone show marked levels of differentiation that is congruent with major topographical features in the region. Our data reveal that mountain ranges do not necessarily represent impassable barriers for this species, as demonstrated by a recent high-elevation colonization event across the Pyrenees that has led to secondary contact and partial introgression among divergent lineages. In contrast, we find indirect evidence of contact zones in lowland areas, in the absence of obvious barriers to dispersal. Our results underscore the importance of multilocus datasets in reconstructing the evolutionary history of independent intraspecific lineages , and understanding the dynamics of introgression upon secondary contact. The degree of mixing among Z. vivipara lineages in secondary contact is similar to that found among different lizard species , suggesting that Z. vivipara lineages in the region may in fact represent separate species. Further research using both molecular and phenotypic data will help determine the extent to which reproductive isolation is underway despite the presence of gene flow . Suture zones are important areas for the study of lineage formation and speciation, as they provide a geographic context for across-taxa comparisons. However, our results on Z. vivipara show that even major biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa depending on their demographic histories and ecological traits. Combining intraspecific studies of genetic variation at the lineage level with ecological data will be essential to better understand phylogeographic differences among taxa and advance our understanding of diversification mechanisms.
Specimens of Z. vivipara were captured at 49 localities across the species range in Northern Spain and Southern France (Table 1). The tail tip of each individual captured was clipped and stored in 95% ethanol at −20°C. Individuals were weighed, measured, photographed, and then released at the site of capture. The capture and handling of lizards complied with national and international ethical guidelines and was conducted under the scientific collecting permits issued by Xunta de Galicia, Junta de Castilla y León, Gobierno del Principado de Asturias, Gobierno de Navarra, Gobierno de Aragón and Generalitat de Catalunya in Spain, and Parc National des Pyrénées in France.
MtDNA sequencing and analysis
Genomic DNA was extracted from tail clippings using a Qiagen™ DNeasy Kit and following the manufacturer’s protocol. The following regions of the mitochondrial DNA were amplified: a 427 bp fragment of the cytochrome b gene (including 21 bp from the adjacent Glu-tRNA) using primers MVZ04 and MVZ05 ; and 964 bp of the NADH dehydrogenase subunit 2 gene (ND2) using primers MetF6 and AsnR2 . Detailed PCR amplification conditions and sequencing details are available in the Additional file 1.
Sequences were automatically aligned using Sequencher 4.1. (GeneCodes, Ann Arbor, Michigan) and variable sites were checked visually for accuracy. Coding gene sequences were unambiguously translated into their amino acid sequence and no double peaks were observed in the chromatographs, suggesting sequences were of mitochondrial origin and not nuclear copies. For most analysis, the two gene fragments (cyt-b and ND2) were concatenated into a single sequence of 1370 bp. Haplotype and nucleotide diversity indices were calculated in DnaSP v5.
We constructed a haplotype network using a maximum parsimony algorithm as implemented in the program TCS 1.21 , setting minimum branch probability at 95%. We used jModelTest to determine the model of sequence evolution for each marker (HKY + G for ND2 and HKY + I for Cyt-b) and obtained an optimal partitioning scheme for the dataset (by gene and by codon position) using the sotware package PartitionFinder . We then constructed a phylogeny of mtDNA haplotypes using Bayesian analysis with MrBayes 3.2.2 (http://mrbayes.csit.fsu.edu/) , and ran two simultaneous parallel runs of four Markov chains each (3 heated and one cold) for 4 million generations and sampled every 100 generations. The first 25% of the trees were excluded as burn-in, and a consensus topology was obtained from the remaining 60,002 samples (30,001 per run). We confirmed MCMC chain convergence using Tracer v1.5 , and output parameters such as an average PSRF value of 1.00, and an average value of 0.0035 for the standard deviation of the split frequencies between simultaneous runs.
To estimate divergence times among lineages we used a coalescence approach that uses Bayesian inference and MCMC simulations to generate posterior probability values for divergence times as implemented in the program Beast. We partitioned our dataset by gene and codon position, and ran the analysis with unlinked substitution and clock models for each partition, yet generated a single tree from both. We conducted preliminary runs using an uncorrelated relaxed lognormal clock to account for rate heterogeneity among lineages, but obtained a value of zero for the ucld.stdev parameter, indicating that our data are clock-like. We therefore used the strict clock in the final analysis, as well as a coalescent model of diversification, and a UPGMA starting tree. After optimizing the priors with preliminary runs, we conducted two final runs of 20 million generations, sampled every 2000 steps. Chain convergence and burn-in were determined with the program Tracer v1.5 , and all ESS values in the final runs were above 600. As a prior for the mutation rate we used a value of 0.0228 subst./site/Myr (0.0114 s/s/Myr per lineage) for Cyt b, a rate estimated from a comprehensive phylogenetic study of several squamate groups  that used as time calibration points the known age of the Canary Islands and the opening of the Strait of Gibraltar 5.3 Mya, well known geological events that have been associated with speciation events in several reptile and amphibian species [70, 71]. We set a lognormal distribution for the prior with a mean of 0.0114 and a log standard deviation of 0.40, so that values of 0.005 and 0.020 s/s/Myr per lineage corresponded to the 5% and 95% quantiles, respectively, as these values encompass the range of rates used for estimating divergence times in other studies on lizards [33, 35, 38, 72].
We tested for past sudden changes in effective population size using Fu’s test of neutrality , which detects departures from neutrality in scenarios characterized by an excess of rare alleles and young mutations in sequences not subjected to recombination. Significant large negative values of F s (generated with Arlequin 3.5) indicate an excess of recent mutations and reject population stasis (Fu, 1997). To estimate the time elapsed since such a sudden population expansion, we used the parameter τ obtained from the distribution of pairwise differences among Cyt b haplotypes, as τ = 2ut, where t is the time elapsed between initial and current population sizes and u = 2μk, where μ is the mutation rate (0.01 substitutions/site/My per lineage) and k is the length of the sequence . We estimated the parameter τ using a generalized non-linear least-square approach  and computed confidence intervals for an alpha level of 0.05 by a parametric bootstrap method based on 100 replicates as implemented in Arlequin 3.5 .
AFLP profiling and analysis
Amplified fragment length polymorphism (AFLP) profiles were generated using the laboratory protocol described in Milá et al. , which was modified slightly from Vos et al.. In brief, whole genomic DNA was digested with restriction enzymes EcoRI and MseI (Tru9) and fragments were ligated to oligonucleotide adapters with T4 DNA ligase. A random sub-sample of fragments was obtained through a pre-selective amplification reaction using primers ET and MC, followed by three selective amplifications using primer pairs ETAG/MCGA, ETAG/MCGT, and ETAG/MCTA, with the E primer fluorescently labelled. Ten pairs of selective amplification primers were tested, but only those producing sufficient loci and repeatable and unambiguously scorable profiles were selected for the analysis. Selectively amplified fragments were run in an ABI 3730XL genetic analyser with a LIZ500 size standard. Peaks were visualized using Genemapper 3.7 and scored manually by a single observer (BM). Only unambiguously scorable loci and individuals were included in the analysis and peaks found in less than 2% of individuals were excluded. Methodological error rate was assessed by running a subset of 10 individuals twice from the pre-selective amplification step. The average per-locus genotyping error rate for the AFLP loci selected, measured as recommended by Bonin et al. , was low at 1.2%.
To determine whether AFLP markers were selectively neutral, we conducted an outlier analysis using the program BayeScan 2.0, which uses differences in allele frequencies between populations and a reversible-jump MCMC algorithm to estimate the locus-specific posterior probability for a model including selection vs. a neutral model . We used 10,000 iterations, a thinning interval of 100, 40 pilot runs, an additional burn-in of 5,000, and a value of 10 for the prior odds for the neutral model .
We estimated allelic frequencies using Zhivotovsky’s  Bayesian method with uniform prior distributions and assuming Hardy-Weinberg genotypic proportions. Genetic diversity (H e ) was based on allele frequencies calculated using the method by Lynch and Milligan  as implemented in the program Aflp-surv v. 1.0 . A matrix of pairwise population F ST values using the F ST analogue ΦPT was calculated with Genalex 6.0 . ΦPT is calculated as VAP/(VAP + VWP) where VAP is the variance among populations and VWP is the variance within populations. Probability values of pairwise ΦPT were based on 10,000 permutations.
To assess genetic structure among samples we conducted a principal coordinate analysis (PCoA; Orloci ) on a genetic distance matrix generated from the binary presence-absence matrix as implemented in Genalex 6.0. .We further examined patterns of population structure using the Bayesian assignment probability test in the program Structure 3.1 . This program uses a Bayesian approach to generate posterior probabilities of assignment of individuals to each of a given number of groups (K). As recommended for dominant markers, we applied a model of no admixture with correlated allele frequencies , and used the Locprior setting, which incorporates sampling locality as a prior without affecting the value of the optimal K . We run five different chains of 1,000,000 steps for each value of K after a burn-in of 100,000 steps. The optimal value of K was estimated by examining the mean Ln likelihood values for different K values from 1 to 10, and by following the method by Evanno et al.  as implemented in Structure Harvester.
Availability of supporting data
Sequences used in this study have been deposited in GenBank under accessions for ND2 ( KF593862-KF593891) and Cyt b (KF593892-KF593905 and KF593907-KF593908) haplotypes.
Comments by three anonymous referees and the subject editor improved previous versions of the manuscript. We are grateful to G. Perron (Parc National des Pyrénées) for providing authorizations to collect tissue samples in the Ossau Valley. BM was funded by an I3P postdoctoral grant from the Spanish Research Council (CSIC), and PSF by a grant from the Spanish Ministry of Education and Science (Programa Ramón y Cajal, RYC-2003-006136). BH was funded by the French National Research Center (CNRS). YSG was financially supported by the European Community’s programme ’Structuring the European Research Area’ under SYNTHESYS at the Museo Nacional de Ciencias Naturales (CSIC), contract number ES-TAF 5305. The project was financed by two grants from the Spanish Ministry of Education and Science (CGL2005-01187, CGL2008-01522) and the Swiss National Science Foundation (PPOOP3_128375) to PSF. Support to cover the publication fee was provided by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
- Avise JC: Phylogeography: the history and formation of species. 2000, Cambridge, Massachusetts: Harvard University PressGoogle Scholar
- Barton N, Hewitt G: Analysis of hybrid zones. Ann Rev Ecol Syst. 1985, 16: 113-148. 10.1146/annurev.es.16.110185.000553.View ArticleGoogle Scholar
- Hewitt GM: Genetic consequeneces of climatic oscillations in the Quaternary. Phil Trans R Soc Lond B. 2004, 359: 183-195. 10.1098/rstb.2003.1388.View ArticleGoogle Scholar
- Hewitt GM: The genetic legacy of the Quaternary ice ages. Nature (London). 2000, 405: 907-913. 10.1038/35016000.View ArticleGoogle Scholar
- Taberlet P, Fumagalli L, Wust-Saucy AG, Cosson JF: Comparative phylogeography and postglacial colonization routes in Europe. Mol Ecol. 1998, 7: 453-464. 10.1046/j.1365-294x.1998.00289.x.PubMedView ArticleGoogle Scholar
- Joger U, Fritz U, Guicking D, Kalyabina-Hauf S, Nagy ZT, Wink M: Phylogeography of western Palaearctic reptiles - Spatial and temporal speciation patterns. Zool Anz. 2007, 246: 293-313. 10.1016/j.jcz.2007.09.002.View ArticleGoogle 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, Netherlands: Springer, 155-188.View ArticleGoogle Scholar
- Surget-Groba Y, Heulin B, Guillaume C-P, Thorpe RS, Kupriyanova L, Vogrin N, Maslak R, Mazzotti S, Venczel M, Ghira I, et al: Intraspecific phylogeography of Lacerta vivipara and the evolution of viviparity. Mol Phylogenet Evol. 2001, 18: 449-459. 10.1006/mpev.2000.0896.PubMedView ArticleGoogle Scholar
- Heulin B, Surget-Groba Y, Guiller A, Guillaume C-P, Deunff J: Comparisons of mitochondrial DNA (mtDNA) sequences (16S rRNA gene) between oviparous and viviparous strains of Lacerta vivipara: a preliminary study. Mol Ecol. 1999, 8 (10): 1627-1631. 10.1046/j.1365-294x.1999.00746.x.PubMedView ArticleGoogle Scholar
- Guillaume C-P, Heulin B, Arrayago MJ, Bea A, Braña F: Refuge areas and suture zones in the Pyrenean and Cantabrian regions: geographic variation of the female MPI sex-linked alleles among oviparous populations of the lizard Lacerta (Zootoca) vivipara. Ecography. 2000, 23: 3-10. 10.1111/j.1600-0587.2000.tb00255.x.View ArticleGoogle Scholar
- Edwards SV, Beerli P: Perspective: gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies. Evolution. 2000, 54: 1839-1854.PubMedGoogle Scholar
- Zink RM, Barrowclough GF: Mitochondrial DNA under siege in avian phylogeography. Mol Ecol. 2008, 17: 2107-2121. 10.1111/j.1365-294X.2008.03737.x.PubMedView ArticleGoogle Scholar
- Moore WS: Inferring phylogenies from mtDNA variation: Mitochondrial-gene trees versus nuclear-gene trees. Evolution. 1995, 49: 718-726. 10.2307/2410325.View ArticleGoogle Scholar
- Excoffier L, Foll M, Petit RJ: Genetic consequences of range expansions. Ann Rev Ecol Syst. 2009, 40: 481-501. 10.1146/annurev.ecolsys.39.110707.173414.View ArticleGoogle Scholar
- Milá B, Toews DPL, Smith TB, Wayne RK: A cryptic contact zone between divergent mitochondrial DNA lineages in southwestern North America supports past introgressive hybridization in the yellow-rumped warbler complex (Aves: Dendroica coronata). Biol J Linn Soc. 2011, 103: 696-706. 10.1111/j.1095-8312.2011.01661.x.View ArticleGoogle Scholar
- McGuire JA, Linkem CW, Koo MS, Hutchison DW, Lappin AK, Orange DI, Lemos-Espinal J, Riddle BR, Jaeger JR, Crandall K: Mitochondrial introgression and incomplete lineage sorting through space and time: phylogenetics of crotaphytid lizards. Evolution. 2007, 61: 2879-2897. 10.1111/j.1558-5646.2007.00239.x.PubMedView ArticleGoogle Scholar
- Funk DJ, Omland KE: Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Ann Rev Ecol Syst. 2003, 34: 397-423. 10.1146/annurev.ecolsys.34.011802.132421.View ArticleGoogle Scholar
- Bazin E, Glémin S, Galtier N: Population size does not influence mitochondrial genetic diversity in animals. Science. 2006, 28: 570-572.View ArticleGoogle Scholar
- Dowling DK, Friberg U, Lindell J: Evolutionary implications of non-neutral mitochondrial genetic variation. Trends Ecol Evol. 2008, 23: 546-554. 10.1016/j.tree.2008.05.011.PubMedView ArticleGoogle Scholar
- Ballard JWO, Whitlock MC: The incomplete natural history of mitochondria. Mol Ecol. 2004, 13: 729-744. 10.1046/j.1365-294X.2003.02063.x.PubMedView ArticleGoogle Scholar
- Knowles LL: Statistical phylogeography. Ann Rev Ecol Syst. 2009, 40: 593-612. 10.1146/annurev.ecolsys.38.091206.095702.View ArticleGoogle Scholar
- Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129: 555-562.PubMed CentralPubMedGoogle Scholar
- Pritchard JK, Wen W: Documentation for Structure Software: Version 2. 2004, Available at: http://pritch.bsd.uchicago.eduGoogle Scholar
- Heulin B, Surget-Groba Y, Sinervo B, Miles D, Guiller A: Dynamics of haplogroup frequencies and survival rates in a contact zone of two mtDNA lineages of the lizard Lacerta vivipara. Ecography. 2011, 34: 436-447. 10.1111/j.1600-0587.2010.06540.x.View ArticleGoogle Scholar
- Cooper SJ, Ibrahim KM, Hewitt G: Postglacial expansion and genome subdivision in the European grasshopper Chorthippus parallelus. Mol Ecol. 1995, 4: 49-60. 10.1111/j.1365-294X.1995.tb00191.x.PubMedView ArticleGoogle Scholar
- Helbig AJ, Salomon M, Bensch S, Seibold I: Male-biased gene flow across an avian hybrid zone: evidence from mitochondrial and microsatellite DNA. J Evol Biol. 2001, 14: 277-287. 10.1046/j.1420-9101.2001.00273.x.View ArticleGoogle Scholar
- Gonçalves H, Martínez-Solano I, Ferrand N, García-París M: Conflicting phylogenetic signal of nuclear vs mitochondrial DNA markers in midwife toads (Anura, Discoglossidae, Alytes): deep coalescence or ancestral hybridization?. Mol Phylogenet Evol. 2007, 44: 494-500. 10.1016/j.ympev.2007.03.001.PubMedView ArticleGoogle Scholar
- Yannic G, Basset P, Hausser J: A new perspective on the evolutionary history of western European Sorex araneus group revealed by paternal and maternal molecular markers. Mol Phylogenet Evol. 2008, 47: 237-250. 10.1016/j.ympev.2008.01.029.PubMedView ArticleGoogle Scholar
- Santucci F, Emerson BC, Hewitt G: Mitochondrial DNA phylogeography of European hedgehogs. Mol Ecol. 1998, 7: 1163-1172. 10.1046/j.1365-294x.1998.00436.x.PubMedView ArticleGoogle Scholar
- Taberlet P, Bouvet J: Mitochondrial DNA polymorphism, phylogeography, and conservation genetics of the brown bear Ursus arctos in Europe. Proc R Soc Lond B Biol Sci. 1994, B 255: 195-200.View ArticleGoogle Scholar
- Pérez-Mellado V: Lacerta vivipara Jaquin, 1797. Fauna Ibérica, vol 10. Edited by: Ramos MA. 1997, Madrid: Museo Nacional de Ciencias Naturales, CSIC, 232-242.Google Scholar
- Heulin B, Guillaume C-P: Le lézard vivipare. Les reptiles de la France, Belgique, Luxembourg et Suisse. Edited by: Vacher JP, Geniez M. 2010, Paris: Biotope, Publications du Museum National d’Histoire Naturelle, 394-401.Google Scholar
- Crochet P-A, Chaline O, Surget-Groba Y, Debain C, Cheylan M: Speciation in mountains: phylogeography and phylogeny of the rock lizards genus Iberolacerta (Reptilia: Lacertidae). Mol Phylogenet Evol. 2004, 30: 860-866. 10.1016/j.ympev.2003.07.016.PubMedView ArticleGoogle Scholar
- Mouret V, Guillaumet A, Cheylan M, Pottier G, Ferchaud AL, Crochet PA: The legacy of ice ages in mountain species: post-glacial colonization of mountain tops rather than current range fragmentation determines mitochondrial genetic diversity in an endemic Pyrenean rock lizard. J Biogeogr. 2011, 38: 1717-1731. 10.1111/j.1365-2699.2011.02514.x.View ArticleGoogle Scholar
- Paulo OS, Dias C, Bruford MW, Jordan WC, Nichols RA: The persistence of Pliocene populations through the Pleistocene climatic cycles: evidence from the phylogeography of an Iberian lizard. Proc R Soc Lond B. 2001, 268: 1625-1630. 10.1098/rspb.2001.1706.View ArticleGoogle Scholar
- Pinho C, Harris DJ, Ferrand N: Contrasting patterns of population subdivision and historical demography in three western Mediterranean lizard species inferred from mitochondrial DNA variation. Mol Ecol. 2007, 16: 1191-1205. 10.1111/j.1365-294X.2007.03230.x.PubMedView ArticleGoogle Scholar
- Alexandrino J, Froufe E, Arntzen JW, Ferrand N: Genetic dubdivision, glacial refucia and postglacial recolonization in the golden-striped salamander, Chioglossa lusitanica (Amphibia: Urodela). Mol Ecol. 2000, 9: 771-781. 10.1046/j.1365-294x.2000.00931.x.PubMedView ArticleGoogle Scholar
- Miraldo A, Hewitt GM, Paulo OS, Emerson BC: Phylogeography and demographic history of Lacerta lepida in the Iberian Peninsula: multiple refugia, range expansions and secondary contact zones. BMC Evol Biol. 2011, 11: 170-10.1186/1471-2148-11-170.PubMed CentralPubMedView ArticleGoogle Scholar
- Fitze PS, González-Jimena V, San-José LM, San Mauro D, Aragón P, Suárez T, Zardoya R: Integrative analyses of speciation and divergence in Psammodromus hispanicus (Squamata: Lacertidae). BMC Evol Biol. 2011, 11: 347-10.1186/1471-2148-11-347.PubMed CentralPubMedView ArticleGoogle Scholar
- Martínez-Solano I, Gonçalves HA, Arntzen JW, García-París M: Phylogenetic relationships and biogeography of midwife toads (Discoglossidae: Alytes). J Biogeogr. 2004, 31: 603-618. 10.1046/j.1365-2699.2003.01033.x.View ArticleGoogle Scholar
- Martínez-Solano I, Teixeira D, Buckley D, García-París M: Mitochondrial DNA phylogeography of Lissotriton boscai (Caudata, Salamandridae): evidence for old, multiple refugia in an Iberian endemic. Mol Ecol. 2006, 15: 3375-3388. 10.1111/j.1365-294X.2006.03013.x.PubMedView ArticleGoogle Scholar
- Gonçalves H, Martínez-Solano I, Pereira RJ, Carvalho B, García-París M, Ferrand N: High levels of population subdivision in a morphologically conserved Mediterranean toad (Alytes cisternasii) result from recent, multiple refugia: evidence from mtDNA, microsatellites and nuclear genealogies. Mol Ecol. 2009, 18: 5143-5160. 10.1111/j.1365-294X.2009.04426.x.PubMedView ArticleGoogle Scholar
- Arribas OJ: Morphological variability of the Cantabro-Pyrenean populations of Zootoca vivipara (Jaquin, 1787) with description of a new subspecies. Herpetozoa. 2009, 21: 123-146.Google Scholar
- Stuart-Fox D, Godinho R, Goüy de Bellocq J, Irwin NR, Brito JC, Moussalli A, Široký P, Hugall AF, Baird SJE: Variation in phenotype, parasite load and male competitive ability across a cryptic hybrid zone. PLOS One. 2009, 4: e5677-10.1371/journal.pone.0005677.PubMed CentralPubMedView ArticleGoogle Scholar
- Nunes VL, Miraldo A, Beaumont MA, Butlin RK, Paulo OS: Association of Mc1r variants with ecologically relevant phenotypes in the European ocellated lizard, Lacerta lepida. J Evol Biol. 2011, 24: 2289-2298. 10.1111/j.1420-9101.2011.02359.x.PubMedView ArticleGoogle Scholar
- Babik W, Szymura JM, Rafiński J: Nuclear markers, mitochondrial DNA and male secondary sexual traits variation in a newt hybrid zone (Triturus vulgaris x T. montandoni). Mol Ecol. 2003, 12: 1913-1930. 10.1046/j.1365-294X.2003.01880.x.PubMedView ArticleGoogle Scholar
- Pereira RJ, Monahan WB, Wake DB: Predictors for reproductive isolation in a ring species complex following genetic and ecological divergence. BMC Evol Biol. 2011, 11: 194-10.1186/1471-2148-11-194.PubMed CentralPubMedView ArticleGoogle Scholar
- Buño I, Torroja E, López-Fernández C, Butlin RK, Hewitt GM, Gosálvez J: A hybrid zone between two subspecies of the grasshopper Chorthippus parallelus along the Pyrenees: the west end. Heredity. 1994, 73: 625-634. 10.1038/hdy.1994.170.View ArticleGoogle Scholar
- Salomon M, Hemin Y: Song variation in the chiffchaffs (Phylloscopus collybita) of the western Pyrenees – the contact zone btween the collybita and brehmii forms. Ethology. 1992, 92: 265-282.View ArticleGoogle Scholar
- Haldane J: Sex ratio and unisexual sterility in hybrid animals. J Genet. 1922, 12: 101-109. 10.1007/BF02983075.View ArticleGoogle Scholar
- Arntzen JW, Wallis GP: Restricted gene flow in a moving hybrid zone of the newts Triturus cristatus and T. marmoratus in western France. Evolution. 1991, 45: 805-826. 10.2307/2409691.View ArticleGoogle Scholar
- Martinsen GD, Whitham TG, Turek RJ, Keim P: Hybrid populations selectively filter gene introgression between species. Evolution. 2001, 55: 1325-1335.PubMedView ArticleGoogle Scholar
- Rohwer S, Bermingham E, Wood C: Plumage and mitochondrial DNA haplotype variation across a moving hybrid zone. Evolution. 2001, 55: 405-422.PubMedView ArticleGoogle Scholar
- Chan KMA, Levin SA: Leaky prezygotic isolation and porous genomes: rapid introgression of maternally inherited DNA. Evolution. 2005, 59: 720-729.PubMedView ArticleGoogle Scholar
- Krosby M, Rohwer S: A 2000 km genetic wake yields evidence for northern glacial refugia and hybrid zone movement in a pair of songbirds. Proc R Soc B. 2009, 276: 615-621. 10.1098/rspb.2008.1310.PubMed CentralPubMedView ArticleGoogle Scholar
- Sequeira F, Alexandrino J, Rocha S, Arntzen W, Ferrand N: Genetic exchange across a hybrid zone within the Iberian endemic golden-striped salamander, Chioglossa lusitanica. Mol Ecol. 2005, 14: 245-254.PubMedView ArticleGoogle Scholar
- Godinho R, Crespo EG, Ferrand N: The limits of mtDNA phylogeography: complex patterns of population fragmentation, expansion and admixture in the highly structured Iberian lizard Lacerta schreiberi are only revealed by the use of nuclear markers. Mol Ecol. 2008, 17: 4670-4683. 10.1111/j.1365-294X.2008.03929.x.PubMedView ArticleGoogle Scholar
- Pinho C, Kaliontzopoulou A, Carretero MA, Harris DJ, Ferrand N: Genetic admixture between the Iberian endemic lizards Podarcis bocagei and Podarcis carbonelli: evidence for limited natural hybridization and a bimodal hybrid zone. J Zool Syst Evol Res. 2009, 47: 368-377. 10.1111/j.1439-0469.2009.00532.x.View ArticleGoogle Scholar
- Mallet J: Hybridization as an invasion of the genome. Trends Ecol Evol. 2005, 20: 229-237. 10.1016/j.tree.2005.02.010.PubMedView ArticleGoogle Scholar
- Smith MF, Patton JL: Variation in mitochondrial cytochrome b sequences in natural populations of South American akodontine rodents (Muridae: Sigmodontinae). Mol Biol Evol. 1991, 8: 85-103.PubMedGoogle Scholar
- Macey JR, Larson A, Ananjeva NB, Fang ZL, Papenfuss TJ: Two novel gene orders and the role of light-strand replication in rearrangement of the vertebrate mitochondrial genome. Mol Biol Evol. 1997, 14: 91-104. 10.1093/oxfordjournals.molbev.a025706.PubMedView ArticleGoogle Scholar
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.PubMedView ArticleGoogle Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.PubMedView ArticleGoogle Scholar
- Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.PubMedView ArticleGoogle Scholar
- Lanfear R, Calcott B, Ho SYW, Guindon S: Combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012, 29: 1695-1701. 10.1093/molbev/mss020.PubMedView ArticleGoogle Scholar
- Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP: Bayesian inference of phylogeny and its impact on evolutionary biology. Science. 2001, 294: 2310-2314. 10.1126/science.1065889.PubMedView ArticleGoogle Scholar
- Rambaut D, Drummond A: MCMC Trace Analysis Package (version 1.4). 2007, Available at: http://tree.bio.ed.ac.uk/software/tracerGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.PubMed CentralPubMedView ArticleGoogle Scholar
- Carranza S, Arnold EN: A review of the geckos of the genus Hemidactylus (Squamata: Gekkonidae) from Oman based on morphology, mitochondrial and nuclear data, with descriptions of eight new species. Zootaxa. 2012, 3378: 1-95.Google Scholar
- Carranza S, Romano A, Arnold EN, Sotgiu G: Biogeography and evolution of European cave salamanders, Hydromantes (Urodela: Plethodontidae), inferred from mtDNA sequences. J Biogeogr. 2008, 35: 724-738. 10.1111/j.1365-2699.2007.01817.x.View ArticleGoogle Scholar
- Cox SC, Carranza S, Brown RP: Divergence times and colonization of the Canary Islands by Gallotia lizards. Mol Phylogenet Evol. 2010, 56: 747-757. 10.1016/j.ympev.2010.03.020.PubMedView ArticleGoogle Scholar
- Carranza S, Arnold EN, Mateo JA, Geniez M: Relationships and evolution of the North African geckos, Geckonia and Tarentola (Reptilia: Gekkonidae), based on mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol. 2002, 23: 244-256. 10.1016/S1055-7903(02)00024-6.PubMedView ArticleGoogle Scholar
- Fu YX: Statistical neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.PubMed CentralPubMedGoogle Scholar
- Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9: 552-569.PubMedGoogle Scholar
- Schneider S, Excoffier L: Estimation of demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics. 1999, 152: 1079-1089.PubMed CentralPubMedGoogle 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: 564-567. 10.1111/j.1755-0998.2010.02847.x.PubMedView ArticleGoogle Scholar
- Milá B, Carranza S, Guillaume O, Clobert J: Marked genetic structuring and extreme dispersal limitation in the Pyrenean brook newt Calotriton asper (Amphibia: Salamandridae) revealed by genome-wide AFLP but not mtDNA. Mol Ecol. 2010, 19: 108-120. 10.1111/j.1365-294X.2009.04441.x.PubMedView ArticleGoogle Scholar
- Vos P, Hogers R, Bleeker M, Reijans M, van de Lee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M, et al: AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res. 1995, 23: 4407-4414. 10.1093/nar/23.21.4407.PubMed CentralPubMedView ArticleGoogle Scholar
- Bonin A, Bellemain E, Bronken Eidesen P, Pompanon F, Brochmann C, Taberlet P: How to track and assess genotyping errors in population genetics studies. Mol Ecol. 2004, 13: 3261-3273. 10.1111/j.1365-294X.2004.02346.x.PubMedView ArticleGoogle Scholar
- Foll M, Gaggiotti OE: A genome scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics. 2008, 180: 977-993. 10.1534/genetics.108.092221.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhivotovsky LA: Estimating population structure in diploids with multilocus dominant DNA markers. Mol Ecol. 1999, 8: 907-913. 10.1046/j.1365-294x.1999.00620.x.PubMedView ArticleGoogle Scholar
- Lynch M, Milligan BG: Analysis of population genetic structure with RAPD markers. Mol Ecol. 1994, 3: 91-99. 10.1111/j.1365-294X.1994.tb00109.x.PubMedView ArticleGoogle Scholar
- Vekemans X, Beauwens T, Lemaire M, Roldán-Ruíz I: Data from amplified fragment length polymorphism (AFLP) markers show indication of size homoplasy and of a relationship between degree of homoplasy and fragment size. Mol Ecol. 2002, 11: 139-151. 10.1046/j.0962-1083.2001.01415.x.PubMedView ArticleGoogle Scholar
- Peakall R, Smouse PE: Genalex 6: genetic analysis in excel. Population genetic software for research and teaching. Mol Ecol Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar
- Orloci L: Multivariate analysis in vegetation research. 1978, The Hague: Dr W. Junk B. VGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Hubisz MJ, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Mol Ecol Resour. 2009, 9: 1322-1332. 10.1111/j.1755-0998.2009.02591.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedView ArticleGoogle Scholar
- Earl D, vonHoldt BM: Structure Harvester: a website and program for visualizing structure output and implementing the Evanno method. Conserv Gen Resour. 2012, 4: 359-361. 10.1007/s12686-011-9548-7.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. 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.