Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

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.     

Clinal genetic variation and isolation by distance in the European eel Anguilla anguilla (L.)

The genetic variability and structure of the European eel ( Anguilla anguilla L.) in populations throughout Europe was reassessed using 15 allozymic loci, seven of which were polymorphic. Seven sites were sampled on a latitudinal gradient across the natural continental range, extending from southern France to southern Norway. Heterozygosity ( H _e = 0.05) and level of polymorphism (P = 0.43) were comparable to other marine fish. Populations were poorly differentiated ( G _ST = 0.014, F _ST = 0.002), which is not surprising considering the high dispersal capability of the European eel. However, a significant geographical cline was detected at two alleles ( IDH-1 * 100 and GPI-1 * 110 ), and genetic distances ( D _CE) were concordant with geographical coastal distances. Mantel tests, pairwise F _ST's and multidimensional scaling analyses identify three distinct groups: Northern Europe, Western Europe and the Mediterranean Sea. We propose that the clinal genetic structure in the European eel may be due to (1) isolation by distance (as recently detected with microsatellites), (2) temporal reproductive separation, (3) post-larval selective forces, (4) contact between formerly separated groups or (5) some combination thereof.  © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 77 , 509–521     

Conservation genetics of Hebe speciosa (Plantaginaceae) an endangered New Zealand shrub

All extant populations of Hebe speciosa (Plantaginaceae), a threatened endemic New Zealand shrub, were investigated using the amplified fragment length polymorphism technique (AFLP). Genetic diversity indices varied significantly among geographical regions and were positively correlated with population size. Among-population genetic differentiation was high (mean pairwise Φ_ST = 0.47), implying complex historical relationships between disjunct populations and negligible contemporary gene flow. Southern populations exhibited extremely low genetic diversity relative to those found in Northland, suggesting that these populations may be more recent in origin. Patterns of genetic relationship among some populations indicate pre-European Māori dispersal and cultivation. The three northernmost populations were found to contain the majority of the species' remaining genetic diversity and thus, should be a focus for future conservation management. Some southern sites may also be culturally significant as evidence of Māori trade and cultivation of Hebe speciosa .  © 2005 The Linnean Society of London, Botanical Journal of the Linnean Society , 2005, 149 , 229–239.     

Population genetic structure of the butterfly Melitaea didyma (Nymphalidae) along a northern distribution range border

The population genetic structure of the butterfly Melitaea didyma was studied along the northern distribution range border in Central Germany by means of allozyme electrophoresis. Individuals were sampled from a total of 21 habitat patches from four regions, and two provinces. Sampling was designed to estimate local vs. regional differentiation. High levels of variability were found, H _e= 0.14–0.21. The mean expected sample heterozygosity from one region, Mosel, was significantly lower than from the Hammelburg region, H _e= 0.17 and 0.19, respectively. Two hierarchical levels of genetic differentiation were found. Within regions individuals sampled from different patches behaved as belonging to one population with high levels of gene flow (Hammelburg F _ST= 0.015, Mosel F _ST= 0.044), though local isolation barriers did create a substructuring of these populations. The inbreeding coefficients, F _IS, were constant over all sample levels, suggesting a similar distribution of habitat patches within regions. Between regions gene flow was limited. An isolation by distance analysis indicated that the hierarchical structure, at the provincial level, may be breaking down due to isolation of regional populations. A more general observation was that the sampling design may greatly have influenced the estimation of genetic differentiation. Depending on which samples were included, overall F _ST estimates ranged from 0.059–0.090.     

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.     

Genetic structure and population dynamics of a heteroecious plant pathogen Melampsora larici-epitea in short-rotation coppice willow plantations

Complex life strategies are common among plant pathogens belonging to rust fungi ( Uredinales ). The heteroecious willow rust Melampsora larici-epitea produces five spore stages and alternates on larch ( Larix ). To shed light on the epidemiology of this pathogen, amplified fragment length polymorphisms (AFLPs) were used to determine the genetic diversity and genetic structure of rust samples collected from coppice willow ( Salix ) plantations at three UK sites (LA, CA and MC) over three sampling dates (September 2000, July 2001 and September 2001). Of the total of 819 isolates, 465 were unique AFLP phenotypes and there was a shift in genotype diversity between the two seasons (0.67 in 2000 and 0.87–0.89 in 2001). No phenotypes were common between the two seasons within a site, suggesting that the rust did not overwinter as an asexual stage within plantations. A temporal analysis detected large amounts of genetic drift ( F _S = 0.15–0.26) between the two seasons and very small effective population sizes ( N _e  =  2–3) within sites. These results all point to a new colonization of the plantations by the rust in the second season (2001). The F _ST-analogue values were Φ_CT = 0.121, Weir and Cockerham's θ = 0.086 and the Bayesian estimate θ^B = 0.087–0.096. The results suggest that the sources of inoculum were somewhat localized and the same sources were mainly responsible for disease epidemics in LA and CA over the two seasons. The relatively low F _ST-values among sites (0.055–0.13) suggest the existence of significant gene flow among the three sampled sites.     

The phylogeographical and management implications of genetic population structure in the imperiled snuffbox mussel, Epioblasma triquetra (Bivalvia: Unionidae)

Seven populations of the imperiled snuffbox mussel, Epioblasma triquetra , were sampled from across the central basin of North America. Samples were genotyped using 15 microsatellite DNA loci, and maternal history was inferred using mitochondrial DNA (mtDNA) cytochrome  c oxidase subunit-I (COI) sequences. Populations in the Clinch and St Francis rivers were quite distinct in their mtDNA. The population in the St Francis River had a unique, fixed haplotype. Among a suite of haplotypes, the population in the Clinch River had two unique haplotypes of common ancestry. The other populations were dominated by a common haplotype, which also occurred in the Clinch River population. Analysis of DNA microsatellites revealed much greater divergences and showed significant genetic structure between populations in the formerly glaciated regions. Divergence has occurred between the populations, as evidenced by moderate to high fixation indices ( F _ST and R _ST values) and nearly perfect assignment tests. These results indicate the occurrence of three glacial refugia for E. triquetra : the Tennessee River, rivers south of the Ozark Crest, and the lower Ohio River drainage near the confluence with the Mississippi. Populations in the lower Ohio River were likely to be responsible for the postglacial reinvasion into formerly glaciated regions, and into the upper Tennessee River drainage. The population of the St Francis River may constitute a distinct taxonomic entity. Conservation efforts, if necessary for this imperiled species, should not mix populations.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 93 , 371–384.     

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.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Comparative genetic analyses of historical and contemporary collections highlight contrasting demographic histories for the bumble bees Bombus pensylvanicus and B. impatiens in Illinois

Direct comparison of genetic patterns between museum specimens and contemporary collections can be a powerful approach for detecting recent demographic changes. Using microsatellite markers, we examined historical and contemporary genetic variation from an apparently declining bumble bee species, Bombus pensylvanicus , and from a stable species, Bombus impatiens , in central Illinois. For each species, we genotyped specimens from the Illinois Natural History Survey collected from three populations between 1969–1972 and from a resurvey of the same areas conducted in 2008. Population structure in B . pensylvanicus increased markedly over the last four decades (from θ_ST = 0.001 to 0.027) while no structure was detected in B . impatiens for either time period (θ_ST = –0.006 to –0.003). Changes in genetic diversity were not significant for either species, although small reductions were observed for B . pensylvanicus in all three populations. Coalescent simulations incorporating both contemporary and historical samples suggest that this small change is not surprising for recent population declines, as large reductions in genetic diversity were only apparent under the most severe bottleneck scenarios. These results demonstrate how comparisons of genetic patterns between temporal periods and species can help elucidate potential threats to population health and suggest several strategies that might be useful in the conservation of B . pensylvanicus in the Midwestern USA.     

Marked genetic divergence revealed by allozymes among populations of the guppy Poecilia reticulata (Poeciliidae), in Trinidad

Populations of the guppy Poecilia reticulata from six locations in N. Trinidad were examined by starch gel electrophoresis to estimate their degree of genetic divergence. Variability at seven enzyme-coding loci demonstrated that populations differed markedly in allele frequencies with some allelic substitution between sites (for 23 loci, Nei's mean genetic identity, Ī ranged from 0.869–1.00; the coefficient of gene differentiation, G _ST= 0.086; and the absolute differentiation between populations, _m= 0.044). There was a good correspondence between degree of physical isolation and extent of genetic differentiation, as exemplified by the partitioning of gene diversity (Nei's gene diversity analysis = 66% between river basins; 32% within river basins; 2% within rivers), though there was considerable variability in the contribution of individual loci. Most populations were in Hardy-Weinberg equilibrium, and any significant deviations were due exclusively to heterozygote deficiencies. The patterns of population differentiation are discussed in relation to probable past geological and historical events, and present-day evolutionary forces. Notable among these are previous continental land connections, periodic immigration from N.E. South America, and post-colonization events including differential predation, sexual selection, apparent restricted vagility and small effective population sizes.     

Genetic Variation and Clonal Diversity of the Two Divergent Types of Clonal Populations of Leymus chinensis Tzvel on the Song Liao Steppe in the West of Northeastern China

The genetic variation and clonal diversity of two divergent types (grey-green and yellow-green) of clonal populations of Leymus chinensis Tzvel at 14 loci were compared. Total gene diversity (HT) and the coefficient of genetic differentiation (GST) were all higher for the yellow-green type (HT = 0.270; GST =0.186) than for the grey-green type (HT = 0.250; GST = 0.157) of L. chinensis. Rare alleles usually occurred as heterozygotes rather than homozygotes and significant deviations from Hardy-Weinberg equilibrium were found only at a few loci. This indicated that these two types of populations were mainly out-crossing. Clonal diversity, evenness of clones, and mean clone size were not significantly different between the two types. We found that differences between the clone size and genetic variation of the yellow-green type of populations occurred with different climate and habitat population groups. However, for the grey-green type of populations, these genetic variations decreased under conditions of different climate and habitat population groups.     

POPULATION STRUCTURE IN AN INSHORE CETACEAN REVEALED BY MICROSATELLITE AND mtDNA ANALYSIS: BOTTLENOSE DOLPHINS (TURSIOPS SP.) IN SHARK BAY, WESTERN AUSTRALIA

We examined population substructure of bottlenose dolphins ( Tursiops sp). in Shark Bay, Western Australia, using 10 highly polymorphic microsatellite loci, and mitochondrial DNA (mtDNA). For microsatellite analysis, 302 different animals were sampled from seven localities throughout the bay. Analysis of genetic differentiation between sampling localities showed a significant correlation between the number of migrants ( Nm ) calculated from F _ST, R _ST and private alleles, and distance between localities–a pattern of isolation-by-distance. For mtDNA, 220 individuals from all seven localities were sequenced for a 351 base pair fragment of the control region, resulting in eight haplotypes, with two distinct clusters of haplotypes. Values of F _ST and (φ)_ST for mtDNA yielded statistically significant differences, mostly between localities that were not adjacent to each other, suggesting female gene flow over a scale larger than the sampled localities. We also observed a significant correlation between the number of female migrants calculated from F _ST and φ_ST and the distance of sampling localities. Our results indicate that dispersal in female dolphins in Shark Bay is more restricted than that of males.     

Microsatellite variation in cod Gadus morhua throughout its geographic range

Previous genetic studies using neutral markers such as allozymes, mtDNA and minisatellite loci have demonstrated varying amounts of population structure in cod Gadus morhua throughout the Atlantic. Microsatellite loci, which are potentially the most informative of presently available neutral genetic markers, have been applied extensively within western and eastern Atlantic areas but not on a range-wide basis. In the present study, six microsatellite DNA loci were used to screen cod samples from nine locations throughout the geographic range from the Scotian Shelf in the West Atlantic to the Barents and Baltic Seas in the east. Overall F _ST value was 0·03 ( P = < 0·001) across all samples. Statistically significant population differences over all loci combined were evident between more geographically distant samples, using either heterogeneity tests or F _ST analysis, with at least one locus showing significant differences between all samples (prior to Bonferroni correction). A significant correlation was observed between genetic and geographical distance, suggesting a higher level of historical and contemporary gene flow between adjacent populations than more distant populations. Samples from either end of the geographic range (Scotian Shelf and Baltic Sea) were particularly distinct when analysed using the STRUCTURE programme and also showed a high level of self-assignment when individuals of either the Scotian Shelf or Baltic Sea were tested against the entire data set. The present microsatellite study demonstrates a high level of geographic population structure between the western Atlantic, middle and eastern Atlantic and Baltic Sea, and thus, the findings should be useful in devising overall management and conservation strategies for the species.     

Genetic diversity and genetic structure of populations of Ranunculus japonicus Thunb. (Ranunculaceae)

Abstract In order to clarify the genetic diversity and population structure of Ranunculus japonicus , allozymic analysis was conducted on 60 populations in southwestern Japan. Considerable genetic variati ons were detected among the populations of R. japonicus . The genetic diversities within species ( H _es = 0.215) and within populations ( H _ep = 0.172) were slightly higher than those of other perennial herbs with widespread distribution and outcrossing plants. Significantly higher values of fixation index were detected in some populations, which might have arisen from restricted mating partners. The majority of genetic variation (approx. 80%) resided within a population and a moderate level of genetic differentiation ( G _ST = 0.203) was observed among populations. The F _ST value (0.203) suggests the existence of a substantial population structure in this species. The highly significant correlation between geographic distance and F _ST values indicates that isolation by distance has played an important role in the construction of the genetic structure of this species.     

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.     

Population genetic structure and effective population size of ayu (Plecoglossus altivelis), an amphidromous fish

The temporal and spatial population genetic structure of ayu Plecoglossus altivelis (Salmoniformes: Plecoglossidae), an amphidromous fish, was examined using analysis of variation at six microsatellite DNA loci. Intracohort genetic diversities, as measured by the number of alleles and heterozygosity, were similar among six cohorts (2001–2006) within a population (Nezugaseki River), with the mean number of alleles per cohort ranging from 11·0 to 12·5 and the expected heterozygosity ranging from 0·74 to 0·77. Intrapopulational genetic diversities were also similar across the three studied populations along the 50 km coast, with the mean number of alleles and the expected heterozygosity ranging from 11·33 to 11·67 and from 0·75 to 0·76, respectively. The authors observed only one significant difference in pair-wise population differentiation ( F _ST-value) between the cohorts within a population and among three populations. Estimates of the effective population size ( N _e) based on maximum-likelihood method yielded small values (ranging from 94·8 to 135·5), whereas census population size ranged from c. 4800 to 24 000. As a result, the ratio of annual effective population sizes to census population size ( N _e/ N ) ranged from 0·004 to 0·023. These estimates of N _e/ N agree more closely with estimates for marine fishes than that of the larger estimates for freshwater fishes. The present study suggests that ayu which is highly fecund and shows low survival during the early life stages is also characterized by having low value of N _e/ N , similar to marine species with a pelagic life cycle.     

Microsatellite data reveals weak population substructuring in Copadichromis sp.'virginalis kajose', a demersal cichlid from Lake Malawi, Africa

Small but significant differences were found in allele frequencies among five populations (overall F _ST estimate (θ)=0·004, P=0·006; overall R _ST estimate (RHO)=0·019, P <0·00001) of the demersal cichlid Copadichromis sp.'virginalis kajose', collected from five locations in Lake Malawi. Pairwise F _ST estimates revealed significant differences between the most southerly population (Cape Maclear), and the three most northerly populations (Mbamba Bay, Metangula and Chilola). Pairwise R _ST estimates also revealed significant differences between some populations, but no geographical pattern was discernible. There was no evidence of isolation by distance using either the shortest straight-line distance between samples, or the distance around the shoreline following a 50 m depth contour. F _ST estimates were considerably lower than found in previous studies on the mbuna (rock-dwelling species), but higher than those found in a study of three pelagic cichlid species from Lake Malawi. Substructuring in C. sp.'virginalis kajose' appears to be on a similar scale to the Atlantic cod.     

Gene flow and melanism in garter snakes revisited: a comparison of molecular markers and island vs. coalescent models

Within populations, the stochastic effect of genetic drift and deterministic effect of natural selection are potentially weakened or altered by gene flow among populations. The influence of gene flow on Lake Erie populations of the common garter snake has been of particular interest because of a discontinuous colour pattern polymorphism (striped vs. melanistic) that is a target of natural selection. We reassessed the relative contributions of gene flow and genetic drift using genetic data and population size estimates. We compared all combinations of two marker systems and two analytical approaches to the estimation of gene flow rates: allozymes (data previously published), microsatellite DNA (new data), the island model ( F _ST-based approach), and a coalescence-based approach. For the coalescence approach, mutation rates and sampling effects were also investigated. While the two markers produced similar results, gene flow based on F _ST was considerably higher (Nm > 4) than that from the coalescence-based method (Nm < 1). Estimates of gene flow are likely to be inflated by lack of migration-drift equilibrium and changing population size. Potentially low rates of gene flow (Nm < 1), small population size at some sites, and positive correlations of number of microsatellite DNA alleles and island size and between M , mean ratio of number of alleles to range in allele size, and island size suggest that in addition to selection, random genetic drift may influence colour pattern frequencies. © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 79, 389–399.