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.     

Genomic evidence for convergent evolution of a key trait underlying divergence in island birds

Reproductive isolation can be initiated by changes in one or a few key traits that prevent random mating among individuals in a population. During the early stages of speciation, when isolation is often incomplete, there will be a heterogeneous pattern of differentiation across regions of the genome between diverging populations, with loci controlling these key traits appearing the most distinct as a result of strong diversifying selection. In this study, we used Illumina‐sequenced ddRAD tags to identify genomewide patterns of differentiation in three recently diverged island populations of the Monarcha castaneiventris flycatcher of the Solomon Islands. Populations of this species have diverged in plumage colour, and these differences in plumage colour, in turn, are used in conspecific recognition and likely important in reproductive isolation. Previous candidate gene sequencing identified point mutations in MC1R and ASIP, both known pigmentation genes, to be associated with the difference in plumage colour between islands. Here, we show that background levels of genomic differentiation based on over 70,000 SNPs are extremely low between populations of distinct plumage colour, with no loci reaching the level of differentiation found in either candidate gene. Further, we found that a phylogenetic analysis based on these SNPs produced a taxonomy wherein the two melanic populations appear to have evolved convergently, rather than from a single common ancestor, in contrast to their original classification as a single subspecies. Finally, we found evidence that the pattern of low genomic differentiation is the result of both incomplete lineage sorting and gene flow between populations.     

Genetic population structure and introgression in Anopheles dirus mosquitoes in South-east Asia

Genetic structure and species relationships were studied in three closely related mosquito species, Anopheles dirus A, C and D in Thailand using 11 microsatellite loci and compared with previous mitochondrial DNA (mtDNA) data on the same populations. All three species were well differentiated from each other at the microsatellite loci. Given the almost complete absence of mtDNA differentiation between An. dirus A and D, this endorses the previous suggestion of mtDNA introgression between these species. The high degree of differentiation between the northern and southern population of An. dirus C (RST = 0.401), in agreement with mtDNA data, is suggestive of incipient species. The lack of genetic structure indicated by microsatellites in four populations of An. dirus A across northern Thailand also concurs with mtDNA data. However, in An. dirus D a limited but significant level of structure was detected by microsatellites over ~400 km in northern Thailand, whereas the mtDNA detected no population differentiation over a much larger area (>1200 km). There is prior evidence for population expansion in the mtDNA. If this is due to a selective sweep originating in An. dirus D, the microsatellite data may indicate greater barriers to gene flow within An. dirus D than in species A. Alternatively, there may have been historical introgression of mtDNA and subsequent demographic expansion which occurred first in An. dirus D so enabling it to accumulate some population differentiation. In the latter case the lack of migration-drift equilibrium precludes the inference of absolute or relative values of gene flow in An. dirus A and D.     

Population genetic structure of savannah elephants in Kenya: conservation and management implications

We investigated population genetic structure and regional differentiation among African savannah elephants in Kenya using mitochondrial and microsatellite markers. We observed mitochondrial DNA (mtDNA) nucleotide diversity of 1.68% and microsatellite variation in terms of average number of alleles, expected and observed heterozygosities in the total study population of 10.20, 0.75, and 0.69, respectively. Hierarchical analysis of molecular variance of mtDNA variation revealed significant differentiation among the 3 geographical regions studied (F(CT) = 0.264; P < 0.05) and a relatively lower differentiation among populations within regions (F(SC) = 0.218; P < 0.0001). Microsatellite variation significantly differentiated among populations within regions (F(SC) = 0.019; P < 0.0001) but not at the regional levels (F(CT) = 0.000; P > 0.500). We attribute the high differentiation at the mitochondrial genome to the matrilineal social structure of elephant populations, female natal philopatry, and probably ancient vicariance. Lack of significant regional differentiation at the nuclear loci vis-a-vis strong differences at mtDNA loci between regions is likely the effect of subsequent homogenization through male-mediated gene flow. Our results depicting 3 broad regional mtDNA groups and the observed population genetic differentiation as well as connectivity patterns should be incorporated in the planning of future management activities such as translocations.     

Genetic divergence among native trout Salmo trutta populations from southern Balkans based on mitochondrial DNA and microsatellite variation

The genetic structure and the phylogenetic relationships among five Balkan populations of trout Salmo trutta that have been classified earlier into five different taxa were studied, using microsatellite and mitochondrial DNA (mtDNA) analyses. The pattern of population differentiation observed at microsatellites differed to that depicted by mtDNA variation, yet both methods indicated a very strong partitioning of the genetic variation among sampling locations. Results thus suggest that conservation strategies should be directed towards preserving the genetic integrity and uniqueness of each population.     

Population genetic structure of the African elephant in Uganda based on variation at mitochondrial and nuclear loci: evidence for male-biased gene flow

A drastic decline has occurred in the size of the Uganda elephant population in the last 40 years, exacerbated by two main factors; an increase in the size of the human population and poaching for ivory. One of the attendant consequences of such a decline is a reduction in the amount of genetic diversity in the surviving populations due to increased effects of random genetic drift. Information about the amount of genetic variation within and between the remaining populations is vital for their future conservation and management. The genetic structure of the African elephant in Uganda was examined using nucleotide variation of mitochondrial control region sequences and four nuclear microsatellite loci in 72 individuals from three localities. Eleven mitochondrial DNA (mtDNA) haplotypes were observed, nine of which were geographically localized. We found significant genetic differentiation between the three populations at the mitochondrial locus while three out of the four microsatellite loci differentiated KV and QE, one locus differentiated KV and MF and no loci differentiated MF and QE. Expected heterozygosity at the four loci varied between 0.51 and 0.84 while nucleotide diversity at the mitochondrial locus was 1.4%. Incongruent patterns of genetic variation within and between populations were revealed by the two genetic systems, and we have explained these in terms of the differences in the effective population sizes of the two genomes and male-biased gene flow between populations.     

Extensive phenotypic diversification coexists with little genetic divergence and a lack of population structure in the White Wagtail subspecies complex (Motacilla alba)

Geographically clustered phenotypes often demonstrate consistent patterns in molecular markers, particularly mitochondrial DNA (mtDNA) traditionally used in phylogeographic studies. However, distinct evolutionary trajectories among traits and markers can lead to their discordance. First, geographic structure in phenotypic traits and nuclear molecular markers can be co‐aligned but inconsistent with mtDNA (mito‐nuclear discordance). Alternatively, phenotypic variation can have little to do with patterns in neither mtDNA nor nuclear markers. Disentangling between these distinct patterns can provide insight into the role of selection, demography and gene flow in population divergence. Here, we examined a previously reported case of strong inconsistency between geographic structure in mtDNA and plumage traits in a widespread polytypic bird species, the White Wagtail (Motacilla alba). We tested whether this pattern is due to mito‐nuclear discordance or discrepancy between morphological evolution and both nuclear and mtDNA markers. We analysed population differentiation and structure across six out of nine commonly recognized subspecies using 17 microsatellite loci and a combination of microsatellites and plumage indices in a comprehensively sampled region of a contact between two subspecies. We did not find support for the mito‐nuclear discordance hypothesis: nuclear markers indicated a subtle signal of genetic clustering only partially consistent with plumage groups, similar to previous findings that relied on mtDNA. We discuss evolutionary factors that could have shaped the intricate patterns of phenotypic diversification in the White wagtail and the role that repeated selection on plumage ‘hotspots’ and hybridization may have played.     

Spatial structure of lemming populations (Dicrostonyx groenlandicus) fluctuating in density

The pattern and scale of the genetic structure of populations provides valuable information for the understanding of the spatial ecology of populations, including the spatial aspects of density fluctuations. In the present paper, the genetic structure of periodically fluctuating lemmings (Dicrostonyx groenlandicus) in the Canadian Arctic was analysed using mitochondrial DNA (mtDNA) control region sequences and four nuclear microsatellite loci. Low genetic variability was found in mtDNA, while microsatellite loci were highly variable in all localities, including localities on isolated small islands. For both genetic markers the genetic differentiation was clear among geographical regions but weaker among localities within regions. Such a pattern implies gene flow within regions. Based on theoretical calculations and population census data from a snap-trapping survey, we argue that the observed genetic variability on small islands and the low level of differentiation among these islands cannot be explained without invoking long distance dispersal of lemmings over the sea ice. Such dispersal is unlikely to occur only during population density peaks.     

Sympatric divergence and clinal variation in multiple coloration traits of Ficedula flycatchers

Geographic variation in phenotypes plays a key role in fundamental evolutionary processes such as local adaptation, population differentiation and speciation, but the selective forces behind it are rarely known. We found support for the hypothesis that geographic variation in plumage traits of the pied flycatcher Ficedula hypoleuca is explained by character displacement with the collared flycatcher Ficedula albicollis in the contact zone. The plumage traits of the pied flycatcher differed strongly from the more conspicuous collared flycatcher in a sympatric area but increased in conspicuousness with increasing distance to there. Phenotypic differentiation (P_ST) was higher than that in neutral genetic markers (F_ST), and the effect of geographic distance remained when statistically controlling for neutral genetic differentiation. This suggests that a cline created by character displacement and gene flow explains phenotypic variation across the distribution of this species. The different plumage traits of the pied flycatcher are strongly to moderately correlated, indicating that they evolve non‐independently from each other. The flycatchers provide an example of plumage patterns diverging in two species that differ in several aspects of appearance. The divergence in sympatry and convergence in allopatry in these birds provide a possibility to study the evolutionary mechanisms behind the highly divergent avian plumage patterns.     

Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci

Two hundred and thirty-six mitochondrial DNA nucleotide sequences were used in combination with polymorphism at four nuclear microsatellite loci to assess the amount and distribution of genetic variation within and between African savannah elephants. They were sampled from 11 localities in eastern, western and southern Africa. In the total sample, 43 haplotypes were identified and an overall nucleotide diversity of 2.0% was observed. High levels of polymorphism were also observed at the microsatellite loci both at the level of number of alleles and gene diversity. Nine to 14 alleles per locus across populations and 44 alleles in the total sample were found. The gene diversity ranged from 0.51 to 0.72 in the localities studied. An analysis of molecular variance showed significant genetic differentiation between populations within regions and also between regions. The extent of subdivision between populations at the mtDNA control region was approximately twice as high as shown by the microsatellite loci (mtDNA F(ST) = 0.59; microsatellite R(ST) = 0.31). We discuss our results in the light of Pleistocene refugia and attribute the observed pattern to population divergence in allopatry accompanied by a recent population admixture following a recent population expansion.     

Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation

We performed a combined analysis of mitochondrial DNA (mtDNA) and microsatellite loci among lake whitefish (Coregonus clupeaformis) populations in order to assess the levels of congruence between both types of markers in defining patterns of genetic structuring, introgressive hybridization and inferring population origins in the hybrid zone of the St. John River basin. A second objective was to test the hypothesis that secondary contact between glacial lineages always resulted in the occurrence of sympatric dwarf and normal whitefish ecotypes. Fish were sampled from 35 populations and polymorphism was screened at mtDNA and six microsatellite loci for a total of 688 and 763 whitefish, respectively. Four lakes harbouring a single whitefish population of normal ecotype admixed with mtDNA haplotypes of different lineages were found. This confirmed that secondary contact between whitefish evolutionary lineages did not always result in the persistence of reproductively isolated ecotypes. Microsatellites further supported the definition of distinct glacial lineages by identifying lineage-specific allelic size groups. They also further supported the hypothesis that ecotypes originated from either a single founding lineage (sympatric divergence) or following secondary contacts between lineages (allopatric divergence), depending on the lake. In general, however, the pattern of population differentiation and introgressive hybridization observed at microsatellites was in sharp contrast with that depicted by mtDNA variation. Both factorial correspondence analysis and analysis of admixture proportion revealed a much more pronounced pattern of introgressive hybridization than depicted by mtDNA analyses. Variable levels of introgression indicated that environmental differences may be as important as the historical contingency of secondary contact in explaining the persistence of sympatric ecotypes and the differential pattern of introgressive hybridization among lakes. Whitefish populations from the St. John River basin hybrid zone represent a rare illustration of a continuum of both morphological and genetic differentiation within a given taxon, spanning from complete introgression to possibly complete reproductive isolation, depending on lakes. Thus, each lake may be viewed as a different temporal snapshot taken throughout the gradual process of speciation.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.     

Relative information content of polymorphic microsatellites and mitochondrial DNA for inferring dispersal and population genetic structure in the olive sea snake, Aipysurus laevis

Polymorphic microsatellites are widely considered more powerful for resolving population structure than mitochondrial DNA (mtDNA) markers, particularly for recently diverged lineages or geographically proximate populations. Weaker population subdivision for biparentally inherited nuclear markers than maternally inherited mtDNA may signal male-biased dispersal but can also be attributed to marker-specific evolutionary characteristics and sampling properties. We discriminated between these competing explanations with a population genetic study on olive sea snakes, Aipysurus laevis. A previous mtDNA study revealed strong regional population structure for A. laevis around northern Australia, where Pleistocene sea-level fluctuations have influenced the genetic signatures of shallow-water marine species. Divergences among phylogroups dated to the Late Pleistocene, suggesting recent range expansions by previously isolated matrilines. Fine-scale population structure within regions was, however, poorly resolved for mtDNA. In order to improve estimates of fine-scale genetic divergence and to compare population structure between nuclear and mtDNA, 354 olive sea snakes (previously sequenced for mtDNA) were genotyped for five microsatellite loci. F statistics and Bayesian multilocus genotype clustering analyses found similar regional population structure as mtDNA and, after standardizing microsatellite F statistics for high heterozygosities, regional divergence estimates were quantitatively congruent between marker classes. Over small spatial scales, however, microsatellites recovered almost no genetic structure and standardized F statistics were orders of magnitude smaller than for mtDNA. Three tests for male-biased dispersal were not significant, suggesting that recent demographic expansions to the typically large population sizes of A. laevis have prevented microsatellites from reaching mutation-drift equilibrium and local populations may still be diverging.     

Phylogeography of the Temminck’s Stint (Calidris temminckii): historical vicariance but little present genetic structure in a regionally endangered Palearctic wader

Aim We study the population differentiation and phylogeography of the Temminck’s Stint (Calidris temminckii). Specifically, we seek signs of past and present population size changes and dispersal events and evaluate management and conservation unit status of the populations. We also study the possibility of introgression as the origin of two mitochondrial DNA (mtDNA) lineages found and estimate the divergence time of the lineages. Location Northern Eurasia. Methods We analysed 583 bp of mtDNA control region domains I and II and 11 microsatellite loci from 13 localities throughout the breeding range. In addition, we used mitochondrial cytochrome c oxidase subunit I (COI), a barcoding gene, to search for signs of introgression. Results More population differentiation was found from microsatellites than from mtDNA, although differentiation was weak in both markers. Signs of past population growth were observed, in addition to more recent decline in some areas. Both control region and COI sequences revealed two maternal lineages coexisting in Fennoscandia and in north‐west Siberia. No signs of introgression were detected. Lineage divergence time was estimated to have occurred during the glacial periods of Pleistocene. Main conclusions Slight differences in mtDNA and microsatellite differentiation and diversity may reflect different features – such as the mutation rate and effective population size – of the markers used, or female‐biased dispersal pattern and high male site‐fidelity of the species. The coexistence of the two mitochondrial lineages is most likely a consequence of post‐glacial mixing of two refugial Pleistocene populations. Based on genetic information alone, global conservation concerns are not imminent. However, fast decline of a marginal Bothnian Bay population and the smallness and remoteness of a Central Yakutian population warrant conservation actions.     

High similarity of genetic population structure in the false clown anemonefish (<Emphasis Type="Italic">Amphiprion ocellaris</Emphasis>) found in microsatellite and mitochondrial control region analysis

Many studies, using various marker systems, have been conducted on the genetic population structure of marine organisms to reveal connectivity among locations and dispersal capabilities. Although mitochondrial sequence markers are widely used, their accuracy is controversially discussed in the context of small scale population genetic discrimination. In the present study, the genetic population structure of the False Clown Anemonefish (Amphiprion ocellaris) in the Indo-Malay Archipelago was revealed by screening six microsatellite loci. Results were congruent to previous mitochondrial control region results, with three major genetic breaks within the Indo-Malay Archipelago. Similar to the mitochondrial DNA (mtDNA) analysis, microsatellite data showed a correlation of genetic structure to historical ocean basin separation during Pleistocene sea level low stands, geographic distance, and dominant current patterns. However, microsatellite divergences are not as deep as the mtDNA divergence, suggesting either that admixture of mtDNA lineages is slower than that of nuclear microsatellites, providing a rather historic picture of separation, or the stronger differentiation signal is due to lower effective population sizes presented by mtDNA. As well, the microsatellite analysis did not give a better resolution on the small scale as expected. This study showed that depending on the genetic markers used, different stages of population separation might be illuminated.     

Microsatellite analysis supports mitochondrial phylogeography of the hellbender (Cryptobranchus alleganiensis)

We investigated genetic diversity of the hellbender (Cryptobranchus alleganiensis) throughout its range in the eastern US using nuclear markers and compared our results to a previously published mitochondrial analysis. A variety of nuclear markers, including protein-coding gene introns and microsatellites were tested but only microsatellites were variable enough for population level analysis. Microsatellite loci showed moderate among population sharing of alleles, in contrast to the reciprocal monophyly exhibited by mitochondrial DNA. However, analyses using F-statistics and Bayesian clustering algorithms showed considerable population subdivision and clustered hellbender populations into the same major groups as the mtDNA. The microsatellites combined with the mtDNA data suggest that gene flow is severely restricted or non-existent among eight major groups, and potentially among populations (rivers) within groups. The combined mtDNA and microsatellite data suggest that the currently recognized hellbender subspecies are paraphyletic. We suggest that the eight independent groups identified in our study should be managed as such, rather than basing conservation decisions on the two named subspecies of hellbender.     

Using 2 Genetic Markers to Discriminate Among Canada Goose Populations in Ohio

Abstract: Canada goose (Branta canadensis) harvest management depends on reliable estimates of harvest composition, and established genetic methods provide an alternative to traditional methods. We expanded upon previous genetic studies by comparing the utility of 6 nuclear microsatellite loci and mitochondrial (mtDNA) control region sequences for discriminating among giant (B. c. maxima) and interior (B. c. interior) populations in Ohio (USA) Canada goose harvests at both individual and population levels. Subspecies and populations exhibited greater differentiation in mtDNA (F_ST = 0.202) than microsatellites (F_ST = 0.021), as would be expected based on differences in effective population size. Neither microsatellites nor mtDNA alone were sufficient for estimating harvest composition at the subspecies or population level in simulations and empirical blind tests using individuals of known origin; however, a combined microsatellite + mtDNA dataset yielded accurate and precise harvest derivations at the subspecies level. Both population-level mixed stock analysis and individual-level assignment tests provided accurate results, but a large proportion of birds could not be assigned with confidence at the individual level. We applied mixed stock analysis and the combined microsatellite + mtDNA dataset to Ohio's 2003–2004 harvest and found that interior populations accounted for 4.9% (95% CI = 1.7–8.0%) of the statewide early season and 9.3% (95% CI = 6.9–11.6%) of the regular and late-season harvested sample. These results suggest that maximum likelihood harvest derivations are highly dependent on the choice of genetic markers. Studies should only employ markers that exhibit sufficient variation and have been shown through simulations and empirical testing to accurately discriminate among the subspecies or management populations of interest.     

Ancestral polymorphisms in genetic markers obscure detection of evolutionarily distinct populations in the endangered Florida grasshopper sparrow (Ammodramus savannarum floridanus)

Genetic analyses of bird subspecies designated as conservation units can address whether they represent units with independent evolutionary histories and provide insights into the evolutionary processes that determine the degree to which they are genetically distinct. Here we use mitochondrial DNA control region sequence and six microsatellite DNA loci to examine phylogeographical structure and genetic differentiation among five North American grasshopper sparrow (Ammodramus savannarum) populations representing three subspecies, including a population of the endangered Florida subspecies (A. s. floridanus). This federally listed taxon is of particular interest because it differs phenotypically from other subspecies in plumage and behaviour and has also undergone a drastic decline in population size over the past century. Despite this designation, we observed no phylogeographical structure among populations in either marker: mtDNA haplotypes and microsatellite genotypes from floridanus samples did not form clades that were phylogenetically distinct from variants found in other subspecies. However, there was low but significant differentiation between Florida and all other populations combined in both mtDNA (FST = 0.069) and in one measure of microsatellite differentiation (theta = 0.016), while the non-Florida populations were not different from each other. Based on analyses of mtDNA variation using a coalescent-based model, the effective sizes of these populations are large (approximately 80,000 females) and they have only recently diverged from each other (< 26,000 ybp). These populations are probably far from genetic equilibrium and therefore the lack of phylogenetic distinctiveness of the floridanus subspecies and minimal genetic differentiation is due most probably to retained ancestral polymorphism. Finally, levels of variation in Florida were similar to other populations supporting the idea that the drastic reduction in population size which has occurred within the last 100 years has not yet had an impact on levels of variation in floridanus. We argue that despite the lack of phylogenetic distinctiveness of floridanus genotypes the observed genetic differentiation and previously documented phenotypic differences justify continued designation of this subspecies as a protected population segment.     

Contrasting patterns of mitochondrial and microsatellite population structure in fragmented populations of greater prairie-chickens

Greater prairie-chickens (Tympanuchus cupido pinnatus) were once found throughout the tallgrass prairie of midwestern North America but over the last century these prairies have been lost or fragmented by human land use. As a consequence, many current populations of prairie-chickens have become isolated and small. This fragmentation of populations is expected to lead to reductions in genetic variation as a result of random genetic drift and a decrease in gene flow. As expected, we found that genetic variation at both microsatellite DNA and mitochondrial DNA (mtDNA) markers was reduced in smaller populations, particularly in Wisconsin. There was relatively little range-wide geographical structure (FST) when we examined mtDNA haplotypes but there was a significant positive relationship between genetic (FST) and geographical distance (isolation by distance). In contrast, microsatellite DNA loci revealed significant geographical structure (FST) and a weak effect of isolation by distance throughout the range. These patterns were much stronger when populations with reduced levels of genetic variability (Wisconsin) were removed from the analyses. This suggests that the effects of genetic drift were stronger than gene flow at microsatellite loci, whereas these forces were in range-wide equilibrium at mtDNA markers. These differences between the two molecular markers may be explained by a larger effective population size (Ne) for mtDNA, which is expected in species such as prairie-chickens that have female-biased dispersal and high levels of polygyny. Our results suggest that historic populations of prairie-chickens were once interconnected by gene flow but current populations are now isolated. Thus, maintaining gene flow may be important for the long-term persistence of prairie-chicken populations.