Experimental evaluation of the usefulness of microsatellite DNA for detecting demographic bottlenecks

Evolutionary and conservation biologists often use molecular markers to evaluate whether populations have experienced demographic bottlenecks that resulted in a loss of genetic variation. We evaluated the utility of microsatellites for detection of recent, severe bottlenecks and compared the amounts of genetic diversity lost in bottlenecks of different sizes. In experimental mesocosms, we established replicate populations by releasing 1, 2, 4 or 8 pairs of the western mosquitofish, Gambusia affinis (Poeciliidae). Using eight polymorphic microsatellite loci, we quantified seven indices of genetic diversity or change that have been used to assess the effects of demographic bottlenecks on populations. We compared indices for the experimentally bottlenecked populations to those for the source population and examined differences between populations established with different numbers of founders. Direct count heterozygosity and the proportion of polymorphic loci were not very sensitive to genetic changes that resulted from the experimental bottlenecks. Heterozygosity excess and expected heterozygosity were useful to varying degrees in the detection of bottlenecks. Allelic diversity and temporal variance in allele frequencies were most sensitive to genetic changes that resulted from the bottlenecks, and the temporal variance method was slightly more correlated with bottleneck size than was allelic diversity. Based on comparisons to a previous study with allozymes, heterozygosity, temporal variance in allele frequencies and allelic diversity, but not proportion of polymorphic loci, appear to be more sensitive to demographic bottlenecks when quantified using microsatellites. We found that analysis of eight highly polymorphic loci was sufficient to detect a recent demographic bottleneck and to obtain an estimate of the magnitude of bottleneck severity.     

Bottlenecks and rescue effects in a fluctuating population of golden-mantled ground squirrels (Spermophilus lateralis)

Mammal species characterized by highly fluctuating populations often maintain genetic diversity in response to frequent demographic bottlenecks, suggesting the ameliorating influence of life history and behavioral factors. Immigration in particular is expected to promote genetic recovery and is hypothesized to be the most likely process maintaining genetic diversity in fluctuating mammal populations. Most demographic bottlenecks have been inferred retrospectively, and direct analysis of a natural population before, during, and after a bottleneck is rare. Using a continuous 10-year dataset detailing the complete demographic and genetic history of a fluctuating population of golden-mantled ground squirrels (Spermophilus lateralis), we analyzed the genetic consequences of a 4-year demographic bottleneck that reduced the population to seven adult squirrels, and we evaluated the potential “rescue effect” of immigration. Analysis of six microsatellite loci revealed that, while a decline in allelic richness was observed during the bottleneck, there was no observed excess of heterozygosity, a characteristic bottleneck signature, and no evidence for heterozygote deficiency during the recovery phase. In addition, we found no evidence for inbreeding depression during or after the bottleneck. By identifying immigrants and analyzing their demographic and genetic contributions, we found that immigration promoted demographic recovery and countered the genetic effects of the bottleneck, especially the loss of allelic richness. Within 3 years both population size and genetic variation had recovered to pre-bottleneck levels, supporting the role of immigration in maintaining genetic variation during bottleneck events in fluctuating populations. Our analyses revealed considerable variation among analytical techniques in their ability to detect genetic bottlenecks, suggesting that caution is warranted when evaluating bottleneck events based on one technique.     

Reliability of genetic bottleneck tests for detecting recent population declines

The identification of population bottlenecks is critical in conservation because populations that have experienced significant reductions in abundance are subject to a variety of genetic and demographic processes that can hasten extinction. Genetic bottleneck tests constitute an appealing and popular approach for determining if a population decline has occurred because they only require sampling at a single point in time, yet reflect demographic history over multiple generations. However, a review of the published literature indicates that, as typically applied, microsatellite-based bottleneck tests often do not detect bottlenecks in vertebrate populations known to have experienced declines. This observation was supported by simulations that revealed that bottleneck tests can have limited statistical power to detect bottlenecks largely as a result of limited sample sizes typically used in published studies. Moreover, commonly assumed values for mutation model parameters do not appear to encompass variation in microsatellite evolution observed in vertebrates and, on average, the proportion of multi-step mutations is underestimated by a factor of approximately two. As a result, bottleneck tests can have a higher probability of 'detecting' bottlenecks in stable populations than expected based on the nominal significance level. We provide recommendations that could add rigor to inferences drawn from future bottleneck tests and highlight new directions for the characterization of demographic history.     

Island colonization and founder effects: the invasion of the Guadeloupe islands by ship rats (Rattus rattus)

The stepwise colonization of islands within an archipelago is typically punctuated by successive founder effects, with each newly founded population being a subsample of the gene pool of the source island. Thus, the genetic signature of successive bottlenecks should be detected when analysing the genetic structure between islands of an archipelago. To test this prediction, we investigated introduced ship rat populations, Rattus rattus (Linnaeus, 1758), in the Guadeloupe Archipelago. Three different methods, commonly named the heterozygosity excess, the mode-shift indicator and the M ratio method, were used to detect bottlenecks from genetic data obtained with eight microsatellite markers on Guadeloupe and two neighbouring islands, Petite-Terre and Fajou. Moreover, a recent eradication failure on Fajou allowed us to test the accuracy of the methods in an 'experimental-like' situation. The results indicate that rats were introduced on Guadeloupe first, which then became the source population for independent secondary colonization of Fajou and Petite-Terre. Moreover, the heterozygosity excess and the mode-shift indicator only detected bottlenecks for the recent colonization of Petite-Terre and the eradication failure on Fajou. However, bottlenecks were detected for all the populations using the M ratio method. This could be interpreted as the remaining signature of the early introduction of the ship rat in the archipelago.     

Low genetic variation support bottlenecks in Scandinavian red deer

Loss of genetic variation from genetic drift during population bottlenecks has been shown for many species. Red deer (Cervus elaphus) may have been exposed to bottlenecks due to founder events during postglacial colonisation in the early Holocene and during known population reductions in the eighteenth and nineteenth centuries. In this study, we assess loss of genetic variation in Scandinavian red deer due to potential bottlenecks by comparing microsatellite (n = 14) and mitochondrial DNA variation in the Norwegian and Swedish populations with the Scottish, Lithuanian and Hungarian populations. Bottlenecks are also assessed from the M ratio of populations, heterozygosity excess and from hierarchical Bayesian analyses of their demographic history. Strong genetic drift and differentiation was identified in both Scandinavian populations. Microsatellite variation was lower in both Scandinavian populations compared with the other European populations and mitochondrial DNA variation was especially low in the Swedish population where only one unique haplotype was observed. Loss of microsatellite alleles was demonstrated by low M ratios in all populations except the Hungarian. M ratios’ were especially low in the Scandinavian populations, indicating additional or more severe bottlenecks. Heterozygosity excess compared with the expectation from the number of observed microsatellite alleles suggested a recent bottleneck of low severity in the Norwegian population. Hierarchical Bayesian coalescent analyses consistently yielded estimates of a large ancestral and a small current population size in all investigated European populations and suggested the onset of population decline to be between 5,000 and 10,000 years ago, which coincide well with postglacial colonisation.     

Microsatellite analysis of a population crash and bottleneck in the Mauna Kea silversword, Argyroxiphium sandwicense ssp. sandwicense (Asteraceae), and its implications for reintroduction

The Mauna Kea silversword, Argyroxiphium sandwicense ssp. sandwicense, has experienced both a severe population crash associated with an increase in alien ungulate populations on Mauna Kea, and a population bottleneck associated with reintroduction. In this paper, we address the genetic consequences of both demographic events using eight microsatellite loci. The population crash was not accompanied by a significant reduction in number of alleles or heterozygosity. However, the population bottleneck was accompanied by significant reductions in observed number of alleles, effective number of alleles, and expected heterozygosity, though not in observed heterozygosity. The effective size of the population bottleneck was calculated using both observed heterozygosities and allele frequency variances. Both methods corroborated the historical census size of the population bottleneck of at most three individuals. The results suggest that: (i) small populations, even those that result from severe reductions in historical population size and extent, are not necessarily genetically depauperate; and (ii) species reintroduction plans need to be conceived and implemented carefully, with due consideration to the genetic impact of sampling for reintroduction.     

Genetic analysis of microsatellite polymorphism in the Elliot's Pheasant (Syrmaticus ellioti) in China

In this study, we reported the population genetic analyses in the Elliot's Pheasant(Syrnaticus ellioti) using seven polymorphism microsatellite loci based on 105 individuals from 4 geographical populations. Departures from Hardy-Weinberg equilibrium were found in four geographical populations. The average number of alleles was 8.86, with a total of 62 alleles across 7 loci; observed heterozygosity (HO) was generally low and the average number was 0.504. For the seven microsatellite loci, the polymorphism information content ranged from 0.549 to 0.860, with an average number 0.712. Population bottlenecks of the four geographical populations were tested by infinite allele mutation model, step-wise mutation model and two-phase mutation model, which found that each population had experienced bottleneck effect during the recent period. Fst analysis across all geographical populations indicated that the genetic differentiaton between the Guizhou geographical population and the Hunan geographical population was highly significant (P<0.001), a finding supported by the far genetic relationship showed by the neighbor-joining tree of four geographical populations based on Nei's unbiased genetic distances. Using hierarchical analysis of molecular variance (Guizhou geographical population relative to all others pooled), we found a low level of the genetic variation among geographical populations and that between groups. However, differences among populations relative to the total sample explained most of the genetic variance (92.84%), which was significant.     

Founder effect and bottleneck signatures in an introduced,insular population of elk

The population of elk (Cervus elaphus roosevelti) inhabiting Afognak Island, Alaska, USA arose from an introduction of 8 individuals from an established population in Washington, USA in 1929, and recently peaked at approximately 1,400 individuals. We examined indices of diversity for 15 microsatellite loci in the Afognak population and compared them to levels in the parent population to determine effects of translocation and demography on genetic variation. The Afognak population differed significantly (P < 0.0001) from the source population in both allele and genotype frequencies. Allelic richness, number of private alleles and multilocus heterozygosity, but not percent loci polymorphic, were significantly lower in Afognak elk. Mean inbreeding coefficients within Afognak (f = 0.019) and source (f = −0.006) populations did not differ significantly from zero. Despite the demographic bottleneck, no evidence of a genetic bottleneck in the Afognak population was detected using a test for heterozygosity excess or mode shift of allele frequencies. Simulations indicated that rapid population growth after the translocation resulted in heterozygosity excess for only 8 years. Conversely, a statistic testing for a bottleneck signature in the ratio of allele number to allele size range (M-ratio) was significant for both the Afognak and source populations, suggesting that the Afognak population had effectively undergone serial bottlenecks. Nonetheless, Afognak failed to show a smaller M-ratio than the parent population, suggesting a failure of that statistic to detect the bottleneck associated with introduction. We show that a severe bottleneck followed by rapid population growth may be undetectable using available tests.     

The genetic consequences of a demographic bottleneck in an introduced biological control insect

Population bottlenecks may result in the loss of genetic diversity, with potentially negative consequences for species of interest in conservation biology, including rare species, invasive species and biological control agents. We examined mtDNA sequence data and four variable microsatellite loci (SSRs) in the melaleuca psyllid Boreioglycaspis melaleucae, which was introduced from Australia to Florida as a biological control agent of the invasive plant Melaleuca quinquenervia. We sampled psyllids in the native and introduced ranges as well as individuals stored from the original founding population. There was a clear loss of mtDNA haplotype diversity, as well as a loss of rare microsatellite alleles, in the introduced range. However, there was little genetic differentiation between the home and introduced ranges, and no evidence for a genetic bottleneck based on an analysis of heterozygosity with the microsatellite markers. Overall, the data showed that the demographic bottleneck had a limited effect on the genetics of populations in the new range.     

How much can parentage analyses tell us about precapture dispersal?

Estimating rates of movement among populations is never simple, and where young animals cannot all be captured at their birth sites, traditional field methods potentially underestimate dispersal rates. Genetic assignment tests appear to hold promise for detecting 'precapture' dispersal, and recent evidence suggests that even on the scale of dispersal between populations, genetic parentage analyses can also be informative. Herein, we examine the performance of both types of analysis with data from a 17-year study of dispersal in banner-tailed kangaroo rats Dipodomys spectabilis. We compare estimates of precapture dispersal from (i) the commonly used parentage analysis program cervus (ii) a pedigree-reconstruction program, MasterBayes, that combines genetic with spatial and other nongenetic information and (iii) genetic assignment procedures implemented by the program geneclass2, with (iv) rates of dispersal observed through recapture of a subset of animals initially marked shortly after weaning. geneclass2 estimates a larger proportion of precapture dispersers than MasterBayes, but both approaches as well as those based on field data alone, suggest that approximately 10% of adults in local populations are immigrants and that interpopulation dispersal is slightly female-biassed. All genetic procedures detect precapture dispersal between populations, but dispersers identified by MasterBayes are particularly compatible with what is independently known about body mass at dispersal, dispersal distance and distance between parents. Parentage analyses have considerable potential to infer the value of this otherwise elusive demographic parameter when most candidate parents can be genotyped and when nongenetic information, especially the distance separating candidate mothers and fathers, can be incorporated into the procedure.     

Selection maintains MHC diversity through a natural population bottleneck

A perceived consequence of a population bottleneck is the erosion of genetic diversity and concomitant reduction in individual fitness and evolutionary potential. Although reduced genetic variation associated with demographic perturbation has been amply demonstrated for neutral molecular markers, the effective management of genetic resources in natural populations is hindered by a lack of understanding of how adaptive genetic variation will respond to population fluctuations, given these are affected by selection as well as drift. Here, we demonstrate that selection counters drift to maintain polymorphism at a major histocompatibility complex (MHC) locus through a population bottleneck in an inbred island population of water voles. Before and after the bottleneck, MHC allele frequencies were close to balancing selection equilibrium but became skewed by drift when the population size was critically low. MHC heterozygosity generally conformed to Hardy-Weinberg expectations except in one generation during the population recovery where there was a significant excess of heterozygous genotypes, which simulations ascribed to strong differential MHC-dependent survival. Low allelic diversity and highly skewed frequency distributions at microsatellite loci indicated potent genetic drift due to a strong founder affect and/or previous population bottlenecks. This study is a real-time examination of the predictions of fundamental evolutionary theory in low genetic diversity situations. The findings highlight that conservation efforts to maintain the genetic health and evolutionary potential of natural populations should consider the genetic basis for fitness-related traits, and how such adaptive genetic diversity will vary in response to both the demographic fluctuations and the effects of selection.     

Evaluating ring-tailed lemurs (Lemur catta) from southwestern Madagascar for a genetic population bottleneck

In light of historical and recent anthropogenic influences on Malagasy primate populations, in this study ring-tailed lemur (Lemur catta) samples from two sites in southwestern Madagascar, Beza Mahafaly Special Reserve (BMSR) and Tsimanampetsotsa National Park (TNP), were evaluated for the genetic signature of a population bottleneck. A total of 45 individuals (20 from BMSR and 25 from TNP) were genotyped at seven microsatellite loci. Three methods were used to evaluate these populations for evidence of a historical bottleneck: M-ratio, mode-shift, and heterozygosity excess tests. Three mutation models were used for heterozygosity excess tests: the stepwise mutation model (SMM), two-phase model (TPM), and infinite allele model (IAM). M-ratio estimations indicated a potential bottleneck in both populations under some conditions. Although mode-shift tests did not strongly indicate a population bottleneck in the recent historical past when samples from all individuals were included, a female-only analysis indicated a potential bottleneck in TNP. Heterozygosity excess was indicated under two of the three mutation models (IAM and TPM), with TNP showing stronger evidence of heterozygosity excess than BMSR. Taken together, these results suggest that a bottleneck may have occurred among L. catta in southwestern Madagascar in the recent past. Given knowledge of how current major stochastic climatic events and human-induced change can negatively impact extant lemur populations, it is reasonable that comparable events in the historical past could have caused a population bottleneck. This evaluation additionally functions to highlight the continuing environmental and anthropogenic challenges faced by lemurs in southwestern Madagascar.     

Genetic analysis of a documented population bottleneck: introduced Bennett's wallabies (Macropus rufogriseus rufogriseus) in New Zealand

Few bottlenecks of wild populations are sufficiently well-documented to constitute models for testing theories about the impact of bottlenecks on genetic variation, and subsequent population persistence. Relevant details of the Bennett's wallaby (Macropus rufogriseus rufogriseus) introduction into New Zealand were recorded (founder number, source and approximate bottleneck duration) and suggest this may provide a rare opportunity to examine the efficacy of tests designed to detect recent bottlenecks in wild populations. We first assessed the accuracy of historic accounts of the introduction using genetic diversity detected in mitochondrial DNA (mtDNA) and at five microsatellite loci. Phylogenetic analyses of mtDNA D-loop sequence haplotypes were consistent with the reported origin of the founders as Tasmania, rather than one of the Bass Strait islands in which Bennett's wallabies are also found. Microsatellite allele frequencies from the Tasmanian source population were then used to seed bottleneck simulations encompassing varying sizes and numbers of generations, in order to assess the severity of bottleneck consistent with diversity observed in the New Zealand population. The results suggested that the founder number was unlikely to have been as small as the three animals suggested by the account of the introduction. Nonetheless, the bottleneck was probably severe; in the range of three to five pairs of wallabies for one to three generations. It resulted in significantly reduced levels of allelic diversity and heterozygosity relative to the source population. This bottleneck is only detectable under the infinite allele model (IAM) and not under the stepwise mutation model (SMM) or the two-phase model (TPM), and possible explanations for this are discussed.     

Biochemical and DNA markers yield strikingly different results regarding variability and differentiation of roe deer (Capreolus capreolus, Artiodactyla: Cervidae) populations from northern Germany

Three mainland and two island roe deer ( Capreolus capreolus ) populations with a total sample size of 105 individuals from Schleswig–Holstein, northern Germany, were analysed with regard to genetic variability within and differentiation among populations as revealed by eight allozyme loci known to be polymorphic in roe deer, eight microsatellite loci and 404 bp of the mitochondrial control region. Surprisingly, the allozymes were completely monomorphic, but microsatellite and control region variability were high. Hypotheses as to demographic reasons for the variability patterns found, including bottlenecks, founder effects and translocations, are put forward. There were no statistically significant differences between the island and the mainland populations in terms of genetic variability as measured by expected heterozygosity, inbreeding coefficient and allelic richness. The correlations of the various variability indices were not statistically significant after Bonferroni correction. Nevertheless, there was a clear tendency for differentiation indices to yield concordant results for microsatellite and mitochondrial markers.     

Genetic consequences of a demographic bottleneck in the Scandinavian arctic fox

Demographic bottlenecks can result in a loss of genetic variation due to the bottleneck effect and subsequent genetic drift. The arctic fox population in Scandinavia went through a severe demographic bottleneck in the early 20th century, and is today classified as critically endangered. In this study, we investigated the pre-bottleneck genetic variation in Scandinavia and compared it to modern samples from Scandinavia and North Russia. Variation in the mtDNA control region and five microsatellite loci was examined through ancient DNA analysis on museum specimens. The microsatellite data from the museum specimens was further used to simulate the expected effect of the bottleneck. The arctic foxes in Scandinavia have lost approximately 25% of the microsatellite alleles and four out of seven mtDNA haplotypes. The results also suggest that the genetic differentiation between North Russia and Scandinavia has doubled over the last 100 years. However, the level of heterozygosity was significantly higher than expected from the simulations. This highlights both the advantage of using museum specimens and the importance of generating specific predictions in conservation genetics.     

A strong genetic footprint of the re‐introduction history of Alpine ibex (Capra ibex ibex)

A population’s neutral genetic variation is a composite of its size, degree of isolation and demographic history. Bottlenecks and founder events increase genetic drift, leading to the loss of genetic variation and increased genetic differentiation among populations. Gene flow has the opposite effects. Thus, gene flow can override the genetic patterns caused by founder events. Using 37 microsatellite loci, we investigated the effects of serial bottlenecks on genetic variation and differentiation among 42 Alpine ibex populations in Switzerland with known re‐introduction histories. We detected a strong footprint of re‐introduction events on contemporary genetic structure, with re‐introduction history explaining a substantial part of the genetic differentiation among populations. As a result of the translocation of a considerable number of individuals from the sole formerly surviving population in northern Italy, most of the genetic variation of the ancestral population is now present in the combined re‐introduced Swiss populations. However, re‐introductions split up the genetic variation among populations, such that each contemporary Swiss population showed lower genetic variation than the ancestral population. As expected, serial bottlenecks had different effects on the expected heterozygosity (He) and standardized number of alleles (sNa). While loss of sNa was higher in the first bottlenecks than in subsequent ones, He declined to a similar degree with each bottleneck. Thus, genetic drift was detected with each bottleneck, even when no loss of sNa was observed. Overall, more than a hundred years after the beginning of this successful re‐introduction programme, re‐introduction history was the main determinant of today’s genetic structure.     

Gene dispersal and outbreeding in a philopatric mammal

Extensive mark-recapture data from banner-tailed kangaroo rats, Dipodomys spectabilis, have shown that both males and females are highly philopatric and suggest the possibility of close inbreeding. However, indirect analyses based on genetic structure appear to contradict direct observations, suggesting longer dispersal distances. Using microsatellite genotypes from most members of a banner-tailed kangaroo rat population during five successive breeding seasons, we ask how relatedness is influenced by dispersal and how it in turn influences mating patterns. The data confirm that, because of philopatry, neighbours are often close relatives. However, patterns of parentage also show that the average distance between mates is large relative to natal dispersal distances and larger than the average distance between nearest opposite-sexed neighbours. Females' mates were often not their nearest male neighbour and many were less related than the nearest male neighbour. We detected multiple paternity in some females' litters; both sexes produce offspring with multiple mates within and between breeding seasons. At the population level, heterozygosities were high and estimates of F were low, indicating that levels of inbreeding were low. Using individual inbreeding coefficients of all juveniles to estimate their parents' relatedness, we found that parental relatedness was significantly lower than relatedness between nearest opposite-sexed adult neighbours. Thus in philopatric populations, long breeding forays can cause genes to move further than individuals disperse, and polyandry may serve to reduce relatedness between mates.     

Detecting past population bottlenecks using temporal genetic data

Population bottlenecks wield a powerful influence on the evolution of species and populations by reducing the repertoire of responses available for stochastic environmental events. Although modern contractions of wild populations due to human-related impacts have been documented globally, discerning historic bottlenecks for all but the most recent and severe events remains a serious challenge. Genetic samples dating to different points in time may provide a solution in some cases. We conducted serial coalescent simulations to assess the extent to which temporal genetic data are informative regarding population bottlenecks. These simulations demonstrated that the power to reject a constant population size hypothesis using both ancient and modern genetic data is almost always higher than that based solely on modern data. The difference in power between the modern and temporal DNA approaches depends significantly on effective population size and bottleneck intensity and less significantly on sample size. The temporal approach provides more power in cases of genetic recovery (via migration) from a bottleneck than in cases of demographic recovery (via population growth). Choice of genetic region is critical, as mutation rate heavily influences the extent to which temporal sampling yields novel information regarding the demographic history of populations.     

Population genetics of the endangered Crowned Solitary Eagle (<Emphasis Type="Italic">Buteogallus coronatus</Emphasis>) in South America

The Crowned Solitary Eagle (Buteogallus coronatus) is one of the rarest and most severely threatened birds of prey in the Neotropical region. We studied levels of neutral genetic diversity, population structure, and the demographic history of the species using 55 contemporary samples covering a large fraction of the species range, which were genotyped at 17 microsatellite loci. Our results indicated genetic homogeneity across the sampled regions, which may be explained by a high dispersal capability of Crowned Solitary Eagles resulting in high gene flow or relatively recent population expansion. Further demographic tests revealed that the species has experienced a recent demographic reduction, but inbreeding was not detected. The existing connectivity between geographically separated populations may have buffered the negative effects of the demographic bottleneck. Alternatively, the demographic reduction may be too recent to detect a genetic signature due to the long generation time of the species. Potential conservation strategies, including the possibility of translocations of individuals, are discussed.     

Genetic Structure of Earthworm Populations at a Regional Scale: Inferences from Mitochondrial and Microsatellite Molecular Markers in Aporrectodea icterica (Savigny 1826)

Despite the fundamental role that soil invertebrates (e.g. earthworms) play in soil ecosystems, the magnitude of their spatial genetic variation is still largely unknown and only a few studies have investigated the population genetic structure of these organisms. Here, we investigated the genetic structure of seven populations of a common endogeic earthworm (Aporrectodea icterica) sampled in northern France to explore how historical species range changes, microevolutionary processes and human activities interact in shaping genetic variation at a regional scale. Because combining markers with distinct modes of inheritance can provide extra, complementary information on gene flow, we compared the patterns of genetic structure revealed using nuclear (7 microsatellite loci) and mitochondrial markers (COI). Both types of markers indicated low genetic polymorphism compared to other earthworm species, a result that can be attributed to ancient bottlenecks, for instance due to species isolation in southern refugia during the ice ages with subsequent expansion toward northern Europe. Historical events can also be responsible for the existence of two divergent, but randomly interbreeding mitochondrial lineages within all study populations. In addition, the comparison of observed heterozygosity among microsatellite loci and heterozygosity expected under mutation-drift equilibrium suggested a recent decrease in effective size in some populations that could be due to contemporary events such as habitat fragmentation. The absence of relationship between geographic and genetic distances estimated from microsatellite allele frequency data also suggested that dispersal is haphazard and that human activities favour passive dispersal among geographically distant populations.