Differentiation with drift: a spatio-temporal genetic analysis of Galápagos mockingbird populations (Mimus spp.)

Small and isolated island populations provide ideal systems to study the effects of limited population size, genetic drift and gene flow on genetic diversity. We assessed genetic diversity within and differentiation among 19 mockingbird populations on 15 Galápagos islands, covering all four endemic species, using 16 microsatellite loci. We tested for signs of drift and gene flow, and used historic specimens to assess genetic change over the last century and to estimate effective population sizes. Within-population genetic diversity and effective population sizes varied substantially among island populations and correlated strongly with island size, suggesting that island size serves as a good predictor for effective population size. Genetic differentiation among populations was pronounced and increased with geographical distance. A century of genetic drift did not change genetic diversity on an archipelago-wide scale, but genetic drift led to loss of genetic diversity in small populations, especially in one of the two remaining populations of the endangered Floreana mockingbird. Unlike in other Galápagos bird species such as the Darwin''s finches, gene flow among mockingbird populations was low. The clear pattern of genetically distinct populations reflects the effects of genetic drift and suggests that Galápagos mockingbirds are evolving in relative isolation.     

Long-term temporal changes of genetic composition in brown trout (Salmo trutta L.) populations inhabiting an unstable environment

The genetic structure of brown trout (Salmo trutta) populations inhabiting rivers on the island of Bornholm in the Baltic Sea was studied on a spatial and temporal scale. Low water levels in the rivers during the summer period are assumed to have a significant impact on the persistence of local populations, possibly resulting in a metapopulation structure. Extinctions may, however, also be buffered by a remnant strategy, whereby juveniles escape to river outlets during periods of drought. We compared polymorphism at seven microsatellite DNA loci in contemporary and past samples collected from 1944 to 1997. A principal component analysis, a hierarchical gene diversity analysis and assignment tests showed that the genetic composition of populations was not temporally stable, and that temporal genetic differentiation was much stronger than spatial differentiation. Genetic variability was high and stable over time. Effective population sizes (Ne) and migration rate (m) were estimated using a maximum-likelihood-based implementation of the temporal method. Ne estimates were low (ranging from 8.3 to 22.9) and estimates of m were high (between 0.23 and 0.99), in contrast to other Danish trout populations inhabiting larger and more environmentally stable rivers (Ne between 39.2 and 289.9 and m between 0.01 and 0.09). Thus, the observed spatio-temporal patterns of genetic differentiation can be explained by drift in small persisting populations, where levels of genetic variation are maintained by strong gene flow. However, observations of rivers devoid of trout suggested that population turnover also takes place. We suggest that Bornholm trout represent a metapopulation where the genetic structure primarily reflects strong drift and gene flow, combined with occasional extinction-recolonization events.     

Effective population sizes and temporal stability of genetic structure in Rana pipiens, the northern leopard frog

Although studies of population genetic structure are very common, whether genetic structure is stable over time has been assessed for very few taxa. The question of stability over time is particularly interesting for frogs because it is not clear to what extent frogs exist in dynamic metapopulations with frequent extinction and recolonization, or in stable patches at equilibrium between drift and gene flow. In this study we collected tissue samples from the same five populations of leopard frogs, Rana pipiens, over a 22-30 year time interval (11-15 generations). Genetic structure among the populations was very stable, suggesting that these populations were not undergoing frequent extinction and colonization. We also estimated the effective size of each population from the change in allele frequencies over time. There exist few estimates of effective size for frog populations, but the data available suggest that ranid frogs may have much larger ratios of effective size (Ne) to census size (Nc) than toads (bufonidae). Our results indicate that R. pipiens populations have effective sizes on the order of hundreds to at most a few thousand frogs, and Ne/Nc ratios in the range of 0.1-1.0. These estimates of Ne/Nc are consistent with those estimated for other Rana species. Finally, we compared the results of three temporal methods for estimating Ne. Moment and pseudolikelihood methods that assume a closed population gave the most similar point estimates, although the moment estimates were consistently two to four times larger. Wang and Whitlock's new method that jointly estimates Ne and the rate of immigration into a population (m) gave much smaller estimates of Ne and implausibly large estimates of m. This method requires knowing allele frequencies in the source of immigrants, but was thought to be insensitive to inexact estimates. In our case the method may have failed because we did not know the true source of immigrants for each population. The method may be more sensitive to choice of source frequencies than was previously appreciated, and so should be used with caution if the most likely source of immigrants cannot be identified clearly.     

Long-term effective population sizes, temporal stability of genetic composition and potential for local adaptation in anadromous brown trout (Salmo trutta) populations

We examined the long-term temporal (1910s to 1990s) genetic variation at eight microsatellite DNA loci in brown trout (Salmo trutta L) collected from five anadromous populations in Denmark to assess the long-term stability of genetic composition and to estimate effective population sizes (Ne). Contemporary and historical samples consisted of tissue and archived scales, respectively. Pairwise thetaST estimates, a hierarchical analysis of molecular variance (amova) and multidimensional scaling analysis of pairwise genetic distances between samples revealed much closer genetic relationships among temporal samples from the same populations than among samples from different populations. Estimates of Ne, using a likelihood-based implementation of the temporal method, revealed Ne >or= 500 in two of three populations for which we have historical data. A third population in a small (3 km) river showed Ne >or= 300. Assuming a stepping-stone model of gene flow we considered the relative roles of gene flow, random genetic drift and selection to assess the possibilities for local adaptation. The requirements for local adaptation were fulfilled, but only adaptations resulting from strong selection were expected to occur at the level of individual populations. Adaptations resulting from weak selection were more likely to occur on a regional basis, i.e. encompassing several populations. Ne appears to have declined recently in at least one of the studied populations, and the documented recent declines of many other anadromous brown trout populations may affect the persistence of local adaptation.     

Loss of genetic variation in hatchery-reared Indian major carp, Catla catla

The hypothesis that effective population sizes are low in hatchery-reared catla ( Catla catla ) from Bangladesh, possibly leading to inbreeding and loss of variation, was tested. The study was based on analysis of seven microsatellite loci in three samples of hatchery-reared catla and four samples representing wild populations. Pair-wise estimates of genetic differentiation between samples were low between wild samples (θ ranging from 0·012 to 0·034), but high between hatchery samples (θ ranging from 0·153 to 0·185), suggesting strong genetic drift in hatcheries. Genetic variation, both in terms of expected heterozygosity and allelic richness, was significantly lower in hatchery samples than in samples of wild catla. Application of a method for reconstructing families among offspring without parental genetic data showed that the hatchery samples consisted of very few half- and full-sib families, whereas the wild samples consisted of a high number of families, suggesting that most individuals were unrelated. Finally, estimation of the effective number of parents ( N _b) in the largest sample of hatchery fish confirmed that effective population size was low ( N _b= 14·9 for multiallelic loci and N _b= 10·6 if alleles were pooled into two composite alleles). The results show that low effective population sizes leading to loss of variation and possibly inbreeding depression should be a matter of serious concern in aquaculture production of catla.     

Stability of genetic structure and effective population size inferred from temporal changes of microsatellite DNA polymorphisms in the land snail Helix aspersa (Gastropoda: Helicidae)

Temporal evolution of genetic variability may have far-reaching consequences for a diverse array of evolutionary processes. Within the polders of the Bay of Mont-Saint-Michel (France), populations of the land snail Helix aspersa are characterized by a metapopulation structure with occasional extinction processes resulting from farming practices. A temporal survey of genetic structure in H . aspersa was carried out using variability at four microsatellite loci, in ten populations sampled two years apart. Levels of within-population genetic variation, as measured by allelic richness, H _e or F _is , did not change over time and similar levels of population differentiation were demonstrated for both sampling years. The extent of genetic differentiation between temporal samples of the same population established (i) a stable structure for six populations, and (ii) substantial genetic changes for four populations. Using classical F -statistics and a maximum likelihood method, estimates of the effective population size ( N _e) illustrated a mixture of stable populations with high N _e, and unstable populations characterized by very small N _e estimates (of 5–11 individuals). Owing to human disturbances, intermittent gene flow and genetic drift are likely to be the predominant evolutionary processes shaping the observed genetic structure. However, the practice of multiple matings and sperm storage is likely to provide a reservoir of variability, minimizing the eroding genetic effects of population size reduction and increasing the effective population size.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 89–102.     

Geographic patterns of genetic differentiation and plumage colour variation are different in the pied flycatcher (Ficedula hypoleuca)

The pied flycatcher is one of the most phenotypically variable bird species in Europe. The geographic variation in phenotypes has often been attributed to spatial variation in selection regimes that is associated with the presence or absence of the congeneric collared flycatcher. Spatial variation in phenotypes could however also be generated by spatially restricted gene flow and genetic drift. We examined the genetic population structure of pied flycatchers across the breeding range and applied the phenotypic Q _ST ( P _ST)– F _ST approach to detect indirect signals of divergent selection on dorsal plumage colouration in pied flycatcher males. Allelic frequencies at neutral markers were found to significantly differ among populations breeding in central and southern Europe whereas northerly breeding pied flycatchers were found to be one apparently panmictic group of individuals. Pairwise differences between phenotypic ( P _ST) and neutral genetic distances ( F _ST) were positively correlated after removing the most differentiated Spanish and Swiss populations from the analysis, suggesting that genetic drift may have contributed to the observed phenotypic differentiation in some parts of the pied flycatcher breeding range. Differentiation in dorsal plumage colouration however greatly exceeded that observed at neutral genetic markers, which indicates that the observed pattern of phenotypic differentiation is unlikely to be solely maintained by restricted gene flow and genetic drift.     

Historical fragmentation of islands and genetic drift in populations of Galápagos lava lizards (Microlophus albemarlensis complex)

The formation of islands following a rise in sea level at the end of Pleistocene is expected to disrupt the equilibrium between genetic drift and gene flow in species with limited ability to disperse. Here, we test the hypothesis that genetic drift in isolation has caused the differentiation of Galápagos lava lizards ( Microlophus albemarlensis complex) found on 12 islets that are likely to have been connected to a larger island, Isla Santa Cruz, during the late Pleistocene. Using 11 microsatellite loci, screened on 524 individuals from 17 localities distributed among and within 15 islands, we found marked differences in allelic richness and heterozygosity. Genetic differentiation was strong (global F _ST = 0.44), with pairwise differences found among populations on islets being larger than differences among three localities sampled within Isla Santa Cruz. As expected under a scenario of drift in isolation, there was a positive correlation of genetic diversity with island size, no relationship between genetic and geographical distance and a strong negative correlation between heterozygosity and measures of genetic differentiation. We conclude that seawater is a significant barrier to gene flow in lava lizards on this timescale. Our results suggest that the shallow diversification of the M. albemarlensis complex is not due to recent gene flow and that genetic drift may have played a substantial role in observed patterns of phenotypic variation among islands.     

Low population genetic differentiation in the Orchidaceae: implications for the diversification of the family

A leading hypothesis for the immense diversity of the Orchidaceae is that skewed mating success and small, disjunct populations lead to strong genetic drift and switches between adaptive peaks. This mechanism is only possible under conditions of low gene flow that lead to high genetic differentiation among populations. We tested whether orchids typically exhibit high levels of population genetic differentiation by conducting a meta‐analysis to compare mean levels of population genetic differentiation (F_ST) between orchids and other diverse families and between rare and common orchids. Compared with other families, the Orchidaceae is typically characterized by relatively low genetic differentiation among populations (mean F_ST = 0.146) at allozyme loci. Rare terrestrial orchids showed higher population genetic differentiation than common orchids, although this value was still lower than the mean for most plant families. All lines of evidence suggest that orchids are typically characterized by low levels of population genetic differentiation, even in species with naturally disjunct populations. As such, we found no strong evidence that genetic drift in isolated populations has played a major role in the diversification of the Orchidaceae. Further research into the diversification of the family needs to unravel the relative roles of biotic and environmental selective pressures in the speciation of orchids.     

Microsatellites reveal high levels of gene flow among populations of the California squid Loligo opalescens

Information on the extent of genetic differentiation among populations of the squid Loligo opalescens is crucial for the conservation of this commercially utilized species. We analysed six highly variable microsatellite loci in 11 collections of L. opalescens from different locations and spawning seasons to estimate the relative influence of two major evolutionary forces, gene flow and genetic drift. Microsatellite allele frequency patterns suggest that gene flow prevents population differentiation in L. opalescens. Tests for genetic differentiation showed homogeneity of the samples with an overall FST/RST of 0.0028/-0.0013. Genetic uniformity among samples from different year classes indicates that allele frequency patterns in L. opalescens are relatively stable over time. However, a more complete and detailed picture of fine-scale allele frequency shifts in this species will require a systematic microsatellite analysis of local populations over consecutive spawning cycles.     

The Rift Valley complex as a barrier to gene flow for Anopheles gambiae in Kenya

Recent studies of Anopheles gambiae, the principal mosquito vector of malaria in Africa, suggested that the eastern Rift Valley and its surrounding areas act as a barrier to gene flow. To quantify the unique effect of these areas on gene flow, we measured genetic variation within and between populations from each side of the Rift. Low differentiation was measured between populations on each side of the Rift (mean FST < 0.008, mean RST < 0.002). However, high differentiation was measured across the Rift (mean FST = 0.104; mean RST = 0.032). Genetic diversity within populations was lower in eastern populations, suggesting that the effective population sizes (Ne) of these populations were lower than those of western populations. We partitioned the overall differentiation across the Rift into three factors: variation in Ne between populations contributed 7-20%; distance contributed 10-30%, and the remainder, corresponding to the unique effect of the Rift was 50-80%. The Rift's effect was highly significant based on FST. The greater sensitivity of FST in measuring differentiation indicated that drift and not mutation generated the differences between populations. Restricted gene exchange across several hundred kilometers on the face of intense human transportation implies that active mosquito dispersal is the major form of migration, and that migration is a multistep process, where step length is relatively short.     

Comparison of quantitative and molecular genetic variation of native vs. invasive populations of purple loosestrife (Lythrum salicaria L., Lythraceae)

Study of adaptive evolutionary changes in populations of invasive species can be advanced through the joint application of quantitative and population genetic methods. Using purple loosestrife as a model system, we investigated the relative roles of natural selection, genetic drift and gene flow in the invasive process by contrasting phenotypical and neutral genetic differentiation among native European and invasive North American populations ( Q _ST −  F _ST analysis). Our results indicate that invasive and native populations harbour comparable levels of amplified fragment length polymorphism variation, a pattern consistent with multiple independent introductions from a diverse European gene pool. However, it was observed that the genetic variation reduced during subsequent invasion, perhaps by founder effects and genetic drift. Comparison of genetically based quantitative trait differentiation ( Q _ST) with its expectation under neutrality ( F _ST) revealed no evidence of disruptive selection ( Q _ST >  F _ST) or stabilizing selection ( Q _ST <  F _ST). One exception was found for only one trait (the number of stems) showing significant sign of stabilizing selection across all populations. This suggests that there are difficulties in distinguishing the effects of nonadaptive population processes and natural selection. Multiple introductions of purple loosestrife may have created a genetic mixture from diverse source populations and increased population genetic diversity, but its link to the adaptive differentiation of invasive North American populations needs further research.     

Widespread gene flow and high genetic variability in populations of water voles Arvicola terrestris in patchy habitats

Theory predicts that the impact of gene flow on the genetic structure of populations in patchy habitats depends on its scale and the demographic attributes of demes (e.g. local colony sizes and timing of reproduction), but empirical evidence is scarce. We inferred the impact of gene flow on genetic structure among populations of water voles Arvicola terrestris that differed in average colony sizes, population turnover and degree of patchiness. Colonies typically consisted of few reproducing adults and several juveniles. Twelve polymorphic microsatellite DNA loci were examined. Levels of individual genetic variability in all areas were high ( H _O= 0.69–0.78). Assignments of juveniles to parents revealed frequent dispersal over long distances. The populations showed negative F _IS values among juveniles, F _IS values around zero among adults, high F _ST values among colonies for juveniles, and moderate, often insignificant, F _ST values for parents. We inferred that excess heterozygosity within colonies reflected the few individuals dispersing from a large area to form discrete breeding colonies. Thus pre-breeding dispersal followed by rapid reproduction results in a seasonal increase in differentiation due to local family groups. Genetic variation was as high in low-density populations in patchy habitats as in populations in continuous habitats used for comparison. In contrast to most theoretical predictions, we found that populations living in patchy habitats can maintain high levels of genetic variability when only a few adults contribute to breeding in each colony, when the variance of reproductive success among colonies is likely to be low, and when dispersal between colonies exceeds nearest-neighbour distances.     

Plant species diversity and composition of experimental grasslands affect genetic differentiation of Lolium perenne populations

Contrasting hypotheses exist about the relationship between plant species diversity and genetic diversity. However, experimental data of species diversity effects on genetic differentiation among populations are lacking. To address this, Lolium perenne was sown with an equal number of seeds in 78 experimental grasslands (Jena Experiment) varying in species richness (1, 2, 4, 8 to 16) and functional group richness and composition (1-4; grasses, legumes, small herbs, tall herbs). Population sizes were determined 4years after sowing, and single-nucleotide polymorphism (SNP) DNA markers based on bulk samples of up to 100 individuals per population were applied. Genetic distances between the field populations and the initially sown seed population increased with sown species richness. The degree of genetic differentiation from the original seed population was largely explained by actual population sizes, which suggests that genetic drift was the main driver of differentiation. Weak relationships among relative allele frequencies and species diversity or actual population sizes, and a positive correlation between actual population sizes and expected heterozygosity also supported the role of genetic drift. Functional composition had additional effects on genetic differentiation of L. perenne populations, indicating a selection because of genotype-specific interactions with other species. Our study supports that genetic diversity is likely to be lower in plant communities with a higher number of interspecific competitors. Negative effects of species richness on population sizes may increase the probability of genetic drift, and selection because of genotype-specific interactions depending on species and genotypic community composition may modulate this relationship.     

Long‐term ‘islands’ in the landscape: low gene flow,effective population size and genetic divergence in the shrub Hakea oldfieldii (Proteaceae)

The genetic structure of disjunct populations is determined by founding genetic properties, demographic processes, gene flow, drift and local selection. We aim to identify the genetic consequences of natural population disjunction at regional and local scales in Hakea oldfieldii using nuclear and plastid markers to investigate long‐term effective population sizes and gene flow, and patterns of diversity and divergence, among populations. Regional divergence was significant as shown by a consistent pattern in principal coordinates, neighbor‐joining and Bayesian analyses, but divergence at the local level was also significant with localized distribution of plastid haplotypes and populations clustering separately in Bayesian analyses. Historical, recent and first‐generation gene flow was low, suggesting that recent habitat fragmentation has not reduced gene migration significantly. Genetic bottlenecks were detected in three populations. Long‐term effective population size was significantly correlated with the number of alleles/locus and observed heterozygosity, but not with census population size, suggesting that the loss of diversity is associated with long‐term changes rather than recent fragmentation. Inbreeding coefficients were significant in only three populations, suggesting that the loss of diversity is linked to drift and bottlenecks associated with demographic processes (local extinction by fires) rather than inbreeding. Historical disjunction as a result of specific ecological requirements, contraction of habitats following drying during the Pleistocene, low gene flow and changes in population size are likely to have been important forces driving divergence through isolation by distance and drift. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 319–334.     

Genetic estimates of contemporary effective population size: what can they tell us about the importance of genetic stochasticity for wild population persistence?

Genetic stochasticity due to small population size contributes to population extinction, especially when population fragmentation disrupts gene flow. Estimates of effective population size ( N _e) can therefore be informative about population persistence, but there is a need for an assessment of their consistency and informative relevance. Here we review the body of empirical estimates of N _e for wild populations obtained with the temporal genetic method and published since Frankham's (1995 ) review. Theoretical considerations have identified important sources of bias for this analytical approach, and we use empirical data to investigate the extent of these biases. We find that particularly model selection and sampling require more attention in future studies.
We report a median unbiased N _e estimate of 260 (among 83 studies) and find that this median estimate tends to be smaller for populations of conservation concern, which may therefore be more sensitive to genetic stochasticity. Furthermore, we report a median N _e/ N ratio of 0.14, and find that this ratio may actually be higher for small populations, suggesting changes in biological interactions at low population abundances. We confirm the role of gene flow in countering genetic stochasticity by finding that N _e correlates strongest with neutral genetic metrics when populations can be considered isolated. This underlines the importance of gene flow for the estimation of N _e, and of population connectivity for conservation in general. Reductions in contemporary gene flow due to ongoing habitat fragmentation will likely increase the prevalence of genetic stochasticity, which should therefore remain a focal point in the conservation of biodiversity.     

Microsatellite DNA analysis of northern pike (Esox lucius L.) populations: insights into the genetic structure and demographic history of a genetically depauperate species

The northern pike Esox lucius L. is a freshwater fish exhibiting pronounced population subdivision and low genetic variability. However, there is limited knowledge on phylogeographical patterns within the species, and it is not known whether the low genetic variability reflects primarily current low effective population sizes or historical bottlenecks. We analysed six microsatellite loci in ten populations from Europe and North America. Genetic variation was low, with the average number of alleles within populations ranging from 2.3 to 4.0 per locus. Genetic differentiation among populations was high (overall θ_ST = 0.51; overall ρ_ST = 0.50). Multidimensional scaling analysis of genetic distances between populations and spatial analysis of molecular variance suggested a single phylogeographical race within the sampled populations from northern Europe, whereas North American and southern European populations were highly distinct. A population from Ireland was monomorphic at all loci, presumably reflecting founder events associated with introduction of the species to the island in the sixteenth century. Bayesian analysis of demographic parameters showed differences in θ (a product of effective population size and mutation rate) among populations from large and small water bodies, but the relative differences in θ were smaller than expected, which could reflect population subdivision within the larger water bodies. Finally, the analyses showed drastic population declines on a time scale of several thousand years within European populations, which we ascribe to either glacial bottlenecks or postglacial founder events.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 91–101.     

Genetic differentiation in pikeperch (Sander lucioperca): the relative importance of gene flow, drift and common history

To investigate the relative importance of homogenizing factors, such as gene flow, and diversifying factors, such as drift, genetic variation in pikeperch ( Sander lucioperca ) in two Fennoscandian regions (North and South) was analysed with microsatellites. Allelic richness and the degree of differentiation were significantly higher in the North ( F _ST= 0·20) than in the South ( F _ST= 0·064). In northern areas, assignments of genotypes were almost exclusively to the population of origin, but in southern areas, the proportion of correct assignments was significantly lower. Most samples exhibited significant heterozygote deficits, and the level of relatedness was higher than expected from randomness. These combined results suggest that there has been more gene flow between populations in southern areas than in northern areas, where the importance of genetic drift has been greater. Effective population sizes were small ( c. 100) and did not differ between areas. The effect of a common history appears minor, and thus processes such as genetic drift and gene flow have been more influential in shaping the patterns of genetic diversity in this species.     

Local adaptation and ecological differentiation under selection,migration, and drift in Arabidopsis lyrata*

How the balance between selection, migration, and drift influences the evolution of local adaptation has been under intense theoretical scrutiny. Yet, empirical studies that relate estimates of local adaptation to quantification of gene flow and effective population sizes have been rare. Here, we conducted a reciprocal transplant trial, a common garden trial, and a whole‐genome‐based demography analysis to examine these effects among Arabidopsis lyrata populations from two altitudinal gradients in Norway. Demography simulations indicated that populations within the two gradients are connected by gene flow (0.1 < 4N_em < 11) and have small effective population sizes (N_e < 6000), suggesting that both migration and drift can counteract local selection. However, the three‐year field experiments showed evidence of local adaptation at the level of hierarchical multiyear fitness, attesting to the strength of differential selection. In the lowland habitat, local superiority was associated with greater fecundity, while viability accounted for fitness differences in the alpine habitat. We also demonstrate that flowering time differentiation has contributed to adaptive divergence between these locally adapted populations. Our results show that despite the estimated potential of gene flow and drift to hinder differentiation, selection among these A. lyrata populations has resulted in local adaptation.     

Genetic structure of continental and island populations of the Mediterranean endemic Cyclamen balearicum (Primulaceae)

Cyclamen balearicum is a self-compatible perennial herb endemic to the western Mediterranean Basin. This species occurs in five geographically isolated terrestrial islands in southern France and on four Balearic islands. In this study, we compare genetic variability and differentiation within and among 11 terrestrial island populations and 17 true island populations. Of nine readable enzyme loci, five were polymorphic in both terrestrial and true islands. F statistics showed a significant heterozygote deficiency in all populations, probably due to high levels of autonomous selfing, restricted gene flow, and subsequent genetic drift. Genetic diversity was higher in terrestrial islands than on the Balearic islands, suggesting that the Balearic islands were colonized when they were in contact with the continent. Population differentiation was greater among terrestrial islands (Fst = 0.417 and Gst = 0.344) than among true islands (Fst = 0.112 and Gst = 0.093). Furthermore, differentiation among populations on the Basses Cévennes terrestrial island was greater (Fst = 0.254) than among populations on the true island of Mallorca (Fst = 0.163). The greater genetic differentiation among terrestrial islands could have been caused by genetic bottlenecks associated with changes in climate and human land use that may have reduced population sizes more severely in terrestrial islands in southern France than on the Balearic islands.