Tracking Natal Dispersal in a Coastal Population of a Migratory Songbird Using Feather Stable Isotope (δ2H, δ34S) Tracers

Adult birds tend to show high fidelity to their breeding territory or disperse over relatively short distances. Gene flow among avian populations is thus expected to occur primarily through natal dispersal. Although natal dispersal is a critical demographic process reflecting the area over which population dynamics take place, low recapture rates of birds breeding for the first time have limited our ability to reliably estimate dispersal rates and distances. Stable isotope approaches can elucidate origins of unmarked birds and so we generated year- and age-specific δ²H and δ³⁴S feather isoscapes (ca. 180 000 km²) of coastal-breeding Ovenbirds (Seiurus aurocapilla) and used bivariate probability density functions to assign the likely natal areas of 35 males recruited as first-year breeders into a population located in northwestern New Brunswick, Canada. Most individuals (80–94% depending on the magnitude of an age correction factor used; i.e. 28–33 out of 35) were classified as residents (i.e. fledged within our study area) and estimated minimum dispersal distances of immigrants were between 40 and 240 km. Even when considering maximum dispersal distances, the likely origin of most first-year breeders was<200 km from our study area. Our method identified recruitment into our population from large geographic areas with relatively few samples whereas previous mark-recapture based methods have required orders of magnitude more individuals to describe dispersal at such geographic scales. Natal dispersal movements revealed here suggest the spatial scale over which many population processes are taking place and we suggest that conservation plans aiming to maintain populations of Ovenbirds and ecologically-similar species should consider management units within 100 or at most 200 km of target breeding populations.     

‘Heritability’ of dispersal propensity in a patchy population

Although dispersal is often considered to be a plastic, condition-dependent trait with low heritability, growing evidence supports medium to high levels of dispersal heritability. Obtaining unbiased estimates of dispersal heritability in natural populations nevertheless remains crucial to understand the evolution of dispersal strategies and their population consequences. Here we show that dispersal propensity (i.e. the probability of dispersal between habitat patches) displays a significant heritability in the collared flycatcher Ficedula albicollis, as estimated by within-family resemblance when accounting for environmental factors. Offspring of dispersing mothers or fathers had a higher propensity to disperse to a new habitat patch themselves. The effect of parental dispersal status was additional to that of local habitat quality, as measured by local breeding population size and success, confirming previous results about condition-dependent dispersal in this population. The estimated levels of heritability varied between 0.30±0.07 and 0.47±0.10, depending on parent–offspring comparisons made and correcting for a significant assortative mating with respect to dispersal status. Siblings also displayed a significant resemblance in dispersal propensity. These results suggest that variation in between-patch natal dispersal in the collared flycatcher is partly genetically determined, and we discuss ways to quantify this genetic basis and its implications.     

Predation and infanticide influence ideal free choice by a parrot occupying heterogeneous tropical habitats

The ideal free distribution (IFD) predicts that organisms will disperse to sites that maximize their fitness based on availability of resources. Habitat heterogeneity underlies resource variation and influences spatial variation in demography and the distribution of populations. We relate nest site productivity at multiple scales measured over a decade to habitat quality in a box-nesting population of Forpus passerinus (green-rumped parrotlets) in Venezuela to examine critical IFD assumptions. Variation in reproductive success at the local population and neighborhood scales had a much larger influence on productivity (fledglings per nest box per year) than nest site or female identity. Habitat features were reliable cues of nest site quality. Nest sites with less vegetative cover produced greater numbers of fledglings than sites with more cover. However, there was also a competitive cost to nesting in high-quality, low-vegetative cover nest boxes, as these sites experienced the most infanticide events. In the lowland local population, water depth and cover surrounding nest sites were related with F. passerinus productivity. Low vegetative cover and deeper water were associated with lower predation rates, suggesting that predation could be a primary factor driving habitat selection patterns. Parrotlets also demonstrated directional dispersal. Pairs that changed nest sites were more likely to disperse from poor-quality nest sites to high-quality nest sites rather than vice versa, and juveniles were more likely to disperse to, or remain in, the more productive of the two local populations. Parrotlets exhibited three characteristics fundamental to the IFD: habitat heterogeneity within and between local populations, reliable habitat cues to productivity, and active dispersal to sites of higher fitness.     

Long-Distance Dispersal via Ocean Currents Connects Omani Clownfish Populations throughout Entire Species Range

Dispersal is a crucial ecological process, driving population dynamics and defining the structure and persistence of populations. Measuring demographic connectivity between discreet populations remains a long-standing challenge for most marine organisms because it involves tracking the movement of pelagic larvae. Recent studies demonstrate local connectivity of reef fish populations via the dispersal of planktonic larvae, while biogeography indicates some larvae must disperse 100–1000 s kilometres. To date, empirical measures of long-distance dispersal are lacking and the full scale of dispersal is unknown. Here we provide the first measure of long-distance dispersal in a coral reef fish, the Omani clownfish Amphiprion omanensis, throughout its entire species range. Using genetic assignment tests we demonstrate bidirectional exchange of first generation migrants, with subsequent social and reproductive integration, between two populations separated by over 400 km. Immigration was 5.4% and 0.7% in each region, suggesting a biased southward exchange, and matched predictions from a physically-coupled dispersal model. This rare opportunity to measure long-distance dispersal demonstrates connectivity of isolated marine populations over distances of 100 s of kilometres and provides a unique insight into the processes of biogeography, speciation and adaptation.     

Sex in space: population dynamic consequences

Sex, so important in the reproduction of bigametic species, is nonetheless often ignored in explorations of the dynamics of populations. Using a growth model of dispersal-coupled populations we can keep track of fluctuations in numbers of females and males. The sexes may differ from each other in their ability to disperse and their sensitivity to population density. As a further complication, the breeding system is either monogamous or polygamous. We use the harmonic mean birth function to account for sex-ratio-dependent population growth in a Moran–Ricker population renewal process. Incorporating the spatial dimension stabilizes the dynamics of populations with monogamy as the breeding system, but does not stabilize the population dynamics of polygamous species. Most notably, in populations coupled with dispersal, where the sexes differ in their dispersal ability there are rarely stable and equal sex ratios. Rather, a two-point cycle, four-point cycle and eventually complex behaviour of sex-ratio dynamics will emerge with increasing birth rates. Monogamy often leads to less noisy sex-ratio dynamics than polygamy. In our model, the sex-ratio dynamics of coupled populations differ from those of an isolated population system, where a stable 50:50 sex ratio is achievable with equal density-dependence costs for females and males. When sexes match in their dispersal ability, population dynamics and sex-ratio dynamics of coupled populations collapse to those of isolated populations.     

No Evidence for Activity Adjustment in Response to Increased Density in Daphnia magna

Increased population density may lead to a decrease in energy available for growth and reproduction via effects on the activity level of individuals. Whilst this may be of particular importance for organisms that compete for defendable resources and/or have a high frequency of social interactions, it is less obvious how individual activity should covary with population density when food resources are not defendable or direct interactions among individuals are negligible. Based on observations that there is a general negative relationship between population density and metabolism it has been suggested that organisms actively reduce activity under increased density to accommodate an expected decrease in food availability. However, in the absence of direct activity measurements the validity of this hypothesis is unclear. Here we test for such anticipatory adjustments of activity levels in the planktonic cladoceran Daphnia magna Straus, a filter feeder whose food resources are not defendable, meaning that density responses can be evaluated in the absence of direct interactions. We tested for changes in activity in response to two separate density cues, one being the direct physical and visual stimuli resulting from being in the vicinity of conspecifics (‘direct density experiment’), and the other being the detection of olfactory cues in their environment (‘olfactory cue experiment’). Ten genetically distinct clones were used to evaluate the potential for genetic variation in these responses. Our measures of activity were highly repeatable, and there was significant variation in activity among clones. Furthermore, this clonal variation was consistent in the ‘direct density’ and ‘olfactory cue’ experiments. The estimated evolvability of the trait (1.3–3.2%) was within the range typically observed in behavioural traits. However, there was no indication that the activity level of individuals respond to population density, either directly to actual density or to olfactory cues representing high density. In this case, the energetic influence of density on population dynamics is sufficiently described by effects on per capita resource availability.     

Limited dispersal in mobile hunter–gatherer Baka Pygmies

Hunter–gatherer Pygmies from Central Africa are described as being extremely mobile. Using neutral genetic markers and population genetics theory, we explored the dispersal behaviour of the Baka Pygmies from Cameroon, one of the largest Pygmy populations in Central Africa. We found a strong correlation between genetic and geographical distances: a pattern of isolation by distance arising from limited parent–offspring dispersal. Our study suggests that mobile hunter–gatherers do not necessarily disperse over wide geographical areas.     

Intrinsic differences in dispersal between populations of gastropods separated by a few metres: Evidence from reciprocal experimental transplantation

Differences between populations in dispersive behaviour can strongly influence population and community structure and have important implications for evolution. Differences in dispersal can be caused by intrinsic differences among populations or by different extrinsic cues. In a semi exposed shore in Co. Wicklow, Ireland, intertidal gastropods, Littorina littorea moved long distances in areas with low natural population density compared to areas with high natural density but were not influenced by density per se. A reciprocal experimental transplantation distinguished between intrinsic and extrinsic factors as causes of differences in dispersive behaviour. Snails from the area of high density were transplanted to the area of low density and vice versa. Comparisons were also carried out with controls (disturbed and undisturbed individuals) in each area. Dispersal of transplanted snails was compared with that of snails translocated within each area. Mean distances displaced and percentage dispersal were monitored after two days period. To test temporal generality the experiment was done twice. The two experiments yielded different outcomes. The majority of evidence supported the intrinsic model: transplanted snails dispersed differently from controls in their new area and similarly to controls in their area of origin. However, in one of the experiments there was some evidence suggesting influence of extrinsic factors or an interaction between extrinsic and intrinsic factors. Further research is required to identify whether the observed intrinsic differences are genetic or due to internal conditions that have been modified by the animal's present and/or past environment. This study reinforces the value of repeated experimental transplantation to characterize factors causing differences in behaviour.     

Probability of successful larval dispersal declines fivefold over 1 km in a coral reef fish

A central question of marine ecology is, how far do larvae disperse? Coupled biophysical models predict that the probability of successful dispersal declines as a function of distance between populations. Estimates of genetic isolation-by-distance and self-recruitment provide indirect support for this prediction. Here, we conduct the first direct test of this prediction, using data from the well-studied system of clown anemonefish (Amphiprion percula) at Kimbe Island, in Papua New Guinea. Amphiprion percula live in small breeding groups that inhabit sea anemones. These groups can be thought of as populations within a metapopulation. We use the x- and y-coordinates of each anemone to determine the expected distribution of dispersal distances (the distribution of distances between each and every population in the metapopulation). We use parentage analyses to trace recruits back to parents and determine the observed distribution of dispersal distances. Then, we employ a logistic model to (i) compare the observed and expected dispersal distance distributions and (ii) determine the relationship between the probability of successful dispersal and the distance between populations. The observed and expected dispersal distance distributions are significantly different (p < 0.0001). Remarkably, the probability of successful dispersal between populations decreases fivefold over 1 km. This study provides a framework for quantitative investigations of larval dispersal that can be applied to other species. Further, the approach facilitates testing biological and physical hypotheses for the factors influencing larval dispersal in unison, which will advance our understanding of marine population connectivity.     

Individual small in‐stream barriers contribute little to strong local population genetic structure five strictly aquatic macroinvertebrate taxa

Water flow in river networks is frequently regulated by man‐made in‐stream barriers. These obstacles can hinder dispersal of aquatic organisms and isolate populations leading to the loss of genetic diversity. Although millions of small in‐stream barriers exist worldwide, their impact on dispersal of macroinvertebrates remains unclear. Therefore, we, therefore, assessed the effects of such barriers on the population structure and effective dispersal of five macroinvertebrate species with strictly aquatic life cycles: the amphipod crustacean Gammarus fossarum (clade 11), three snail species of the Ancylus fluviatilis species complex and the flatworm Dugesia gonocephala. We studied populations at nine weirs and eight culverts (3 pipes, 5 tunnels), built 33–109 years ago, mainly in the heavily fragmented catchment of the river Ruhr (Sauerland, Germany). To assess fragmentation and barrier effects, we generated genome‐wide SNP data using ddRAD sequencing and evaluated clustering, differentiation between populations up‐ and downstream of each barrier and effective migration rates among sites and across barriers. Additionally, we applied population genomic simulations to assess expected differentiation patterns under different gene flow scenarios. Our data show that populations of all species are highly isolated at regional and local scales within few kilometers. While the regional population structure likely results from historical processes, the strong local differentiation suggests that contemporary dispersal barriers exist. However, we identified significant barrier effects only for pipes (for A. fluviatilis II and III) and few larger weirs (>1.3 m; for D. gonocephala). Therefore, our data suggest that most small in‐stream barriers can probably be overcome by all studied taxa frequently enough to prevent fragmentation. However, it remains to be tested if the strong local differentiation is a result of a cumulative effect of small barriers, or if larger in‐stream barriers, land use, chemical pollution, urbanization, or a combination of these factors impede gene flow.     

Historical and contemporary gene dispersal in wild carrot (Daucus carota ssp. carota) populations

Background and Aims

Wild carrot is the ancestor of cultivated carrot and is the most important gene pool for carrot breeding. Transgenic carrot may be released into the environment in the future. The aim of the present study was to determine how far a gene can disperse in wild carrot populations, facilitating risk assessment and management of transgene introgression from cultivated to wild carrots and helping to design sampling strategies for germplasm collections.

Methods

Wild carrots were sampled from Meijendel and Alkmaar in The Netherlands and genotyped with 12 microsatellite markers. Spatial autocorrelation analyses were used to detect spatial genetic structures (SGSs). Historical gene dispersal estimates were based on an isolation by distance model. Mating system and contemporary pollen dispersal were estimated using 437 offspring of 20 mothers with different spatial distances and a correlated paternity analysis in the Meijendel population.

Key Results

Significant SGSs are found in both populations and they are not significantly different from each other. Combined SGS analysis indicated significant positive genetic correlations up to 27 m. Historical gene dispersal σ_g and neighbourhood size N_b were estimated to be 4–12 m [95 % confidence interval (CI): 3–25] and 42–73 plants (95 % CI: 28–322) in Meijendel and 10–31 m (95 % CI: 7–∞) and 57–198 plants (95 % CI: 28–∞) in Alkmaar with longer gene dispersal in lower density populations. Contemporary pollen dispersal follows a fat-tailed exponential-power distribution, implying pollen of wild carrots could be dispersed by insects over long distance. The estimated outcrossing rate was 96 %.

Conclusions

SGSs in wild carrots may be the result of high outcrossing, restricted seed dispersal and long-distance pollen dispersal. High outcrossing and long-distance pollen dispersal suggest high frequency of transgene flow might occur from cultivated to wild carrots and that they could easily spread within and between populations.     

Dispersal of male ortolan buntings away from areas with low female density and a severely male-biased sex ratio

Dispersal is expected to enhance individual fitness, and individuals should thus disperse from areas with poor conditions to areas with more favourable conditions. Few studies have compared conditions before and after dispersal of the same individuals, and in birds little is known about the effects of sex ratio and female density on male dispersal decisions. In this study we examined various fitness-related parameters that adult male ortolan buntings, Emberiza hortulana, might use as cues in their decisions to disperse and settle. The study population has a strongly male-biased sex ratio. Using pairwise comparisons of pre- and post-dispersal conditions, we found that males moved from areas with low female density and a severely male-biased sex ratio to areas with higher female density and a less male-biased sex ratio. Male density and male age structure did not affect male dispersal and settlement. The sex ratio of the pre-dispersal sites was below the population average, but post-dispersal sites were not significantly better than the population average. This suggests that dispersal was triggered by poor conditions, whereas settlement may have been unrelated to the conditions at the new site. In the year of dispersal, males that undertook dispersal were less successful at acquiring a female than males that remained faithful to the site, but in subsequent years there was no difference. We suggest that dispersal might be adaptive under more natural conditions with a more balanced sex ratio. These results show that male ortolan buntings abandon areas with low female density and a severely male-biased sex ratio, and thereby end up in better places, although probably not better than the population average.     

Integrodifference models for persistence in fragmented habitats

Integrodifference models of growth and dispersal are analyzed on finite domains to investigate the effects of emigration, local growth dynamics and habitat heterogeneity on population persistence. We derive the bifurcation structure for a range of population dynamics and present an approximation that allows straighforward calculation of the equilibrium populations in terms of local growth dynamics and dispersal success rates. We show how population persistence in a heterogeneous environment depends on the scale of the heterogeneity relative to the organism's characteristic dispersal distance. When organisms tend to disperse only a short distance, population persistence is dominated by local conditions in high quality patches, but when dispersal distance is relatively large, poor quality habitat exerts a greater influence.     

Spatial Variation in Population Structure and Its Relation to Movement and the Potential for Dispersal in a Model Intertidal Invertebrate

Dispersal, the movement of an individual away from its natal or breeding ground, has been studied extensively in birds and mammals to understand the costs and benefits of movement behavior. Whether or not invertebrates disperse in response to such attributes as habitat quality or density of conspecifics remains uncertain, due in part to the difficulties in marking and recapturing invertebrates. In the upper Bay of Fundy, Canada, the intertidal amphipod Corophium volutator swims at night around the new or full moon. Furthermore, this species is regionally widespread across a large spatial scale with site-to-site variation in population structure. Such variation provides a backdrop against which biological determinants of dispersal can be investigated. We conducted a large-scale study at nine mudflats, and used swimmer density, sampled using stationary plankton nets, as a proxy for dispersing individuals. We also sampled mud residents using sediment cores over 3 sampling rounds (20–28 June, 10–17 July, 2–11 August 2010). Density of swimmers was most variable at the largest spatial scales, indicating important population-level variation. The smallest juveniles and large juveniles or small adults (particularly females) were consistently overrepresented as swimmers. Small juveniles swam at most times and locations, whereas swimming of young females decreased with increasing mud presence of young males, and swimming of large juveniles decreased with increasing mud presence of adults. Swimming in most stages increased with density of mud residents; however, proportionally less swimming occurred as total mud resident density increased. We suggest small juveniles move in search of C. volutator aggregations which possibly act as a proxy for better habitat. We also suggest large juveniles and small adults move if potential mates are limiting. Future studies can use sampling designs over large spatial scales with varying population structure to help understand the behavioral ecology of movement, and dispersal in invertebrate taxa.     

Restricting Prey Dispersal Can Overestimate the Importance of Predation in Trophic Cascades

Predators can affect prey populations and, via trophic cascades, predators can indirectly impact resource populations (2 trophic levels below the predator) through consumption of prey (density-mediated indirect effects; DMIEs) and by inducing predator-avoidance behavior in prey (trait-mediated indirect effects; TMIEs). Prey often employ multiple predator-avoidance behaviors, such as dispersal or reduced foraging activity, but estimates of TMIEs are usually on individual behaviors. We assessed direct and indirect predator effects in a mesocosm experiment using a marine food chain consisting of a predator (toadfish – Opsanus tau), prey (mud crab - Panopeus herbstii) and resource (ribbed mussel – Geukensia demissa). We measured dispersal and foraging activity of prey separately by manipulating both the presence and absence of the predator, and whether prey could or could not disperse into a predator-free area. Consumption of prey was 9 times greater when prey could not disperse, probably because mesocosm boundaries increased predator capture success. Although predator presence did not significantly affect the number of crabs that emigrated, the presence of a predator decreased resource consumption by prey, which resulted in fewer resources consumed for each prey that emigrated in the presence of a predator, and reduced the overall TMIE. When prey were unable to disperse, TMIEs on mussel survival were 3 times higher than the DMIEs. When prey were allowed to disperse, the TMIEs on resource survival increased to 11-times the DMIEs. We found that restricting the ability of prey to disperse, or focusing on only one predator-avoidance behavior, may be underestimating TMIEs. Our results indicate that the relative contribution of behavior and consumption in food chain dynamics will depend on which predator-avoidance behaviors are allowed to occur and measured.     

O father where art thou? Paternity analyses in a natural population of the haploid–diploid seaweed Chondrus crispus

The link between life history traits and mating systems in diploid organisms has been extensively addressed in the literature, whereas the degree of selfing and/or inbreeding in natural populations of haploid–diploid organisms, in which haploid gametophytes alternate with diploid sporophytes, has been rarely measured. Dioecy has often been used as a proxy for the mating system in these organisms. Yet, dioecy does not prevent the fusion of gametes from male and female gametophytes originating from the same sporophyte. This is likely a common occurrence when spores from the same parent are dispersed in clumps and recruit together. This pattern of clumped spore dispersal has been hypothesized to explain significant heterozygote deficiency in the dioecious haploid–diploid seaweed Chondrus crispus. Fronds and cystocarps (structures in which zygotes are mitotically amplified) were sampled in two 25 m² plots located within a high and a low intertidal zone and genotyped at 5 polymorphic microsatellite loci in order to explore the mating system directly using paternity analyses. Multiple males sired cystocarps on each female, but only one of the 423 paternal genotypes corresponded to a field-sampled gametophyte. Nevertheless, larger kinship coefficients were detected between males siring cystocarps on the same female in comparison with males in the entire population, confirming restricted spermatial and clumped spore dispersal. Such dispersal mechanisms may be a mode of reproductive assurance due to nonmotile gametes associated with putatively reduced effects of inbreeding depression because of the free-living haploid stage in C. crispus.     

Evidence that disease-induced population decline changes genetic structure and alters dispersal patterns in the Tasmanian devil

Infectious disease has been shown to be a major cause of population declines in wild animals. However, there remains little empirical evidence on the genetic consequences of disease-mediated population declines, or how such perturbations might affect demographic processes such as dispersal. Devil facial tumour disease (DFTD) has resulted in the rapid decline of the Tasmanian devil, Sarcophilus harrisii, and threatens to cause extinction. Using 10 microsatellite DNA markers, we compared genetic diversity and structure before and after DFTD outbreaks in three Tasmanian devil populations to assess the genetic consequences of disease-induced population decline. We also used both genetic and demographic data to investigate dispersal patterns in Tasmanian devils along the east coast of Tasmania. We observed a significant increase in inbreeding (F_IS pre/post-disease −0.030/0.012, P<0.05; relatedness pre/post-disease 0.011/0.038, P=0.06) in devil populations after just 2–3 generations of disease arrival, but no detectable change in genetic diversity. Furthermore, although there was no subdivision apparent among pre-disease populations (θ=0.005, 95% confidence interval (CI) −0.003 to 0.017), we found significant genetic differentiation among populations post-disease (θ=0.020, 0.010–0.027), apparently driven by a combination of selection and altered dispersal patterns of females in disease-affected populations. We also show that dispersal is male-biased in devils and that dispersal distances follow a typical leptokurtic distribution. Our results show that disease can result in genetic and demographic changes in host populations over few generations and short time scales. Ongoing management of Tasmanian devils must now attempt to maintain genetic variability in this species through actions designed to reverse the detrimental effects of inbreeding and subdivision in disease-affected populations.     

Negative density-dependent dispersal emerges from the joint evolution of density- and body condition-dependent dispersal strategies

Empirical studies have documented both positive and negative density-dependent dispersal, yet most theoretical models predict positive density dependence as a mechanism to avoid competition. Several hypotheses have been proposed to explain the occurrence of negative density-dependent dispersal, but few of these have been formally modeled. Here, we developed an individual-based model of the evolution of density-dependent dispersal. This model is novel in that it considers the effects of density on dispersal directly, and indirectly through effects on individual condition. Body condition is determined mechanistically, by having juveniles compete for resources in their natal patch. We found that the evolved dispersal strategy was a steep, increasing function of both density and condition. Interestingly, although populations evolved a positive density-dependent dispersal strategy, the simulated metapopulations exhibited negative density-dependent dispersal. This occurred because of the negative relationship between density and body condition: high density sites produced low-condition individuals that lacked the resources required for dispersal. Our model, therefore, generates the novel hypothesis that observed negative density-dependent dispersal can occur when high density limits the ability of organisms to disperse. We suggest that future studies consider how phenotype is linked to the environment when investigating the evolution of dispersal.     

Molecular Markers Reveal Limited Population Genetic Structure in a North American Corvid,Clark’s Nutcracker (Nucifraga columbiana)

The genetic impact of barriers and Pleistocene glaciations on high latitude resident species has not been widely investigated. The Clark’s nutcracker is an endemic North American corvid closely associated with Pinus-dominated forests. The nutcracker’s encompasses known barriers to dispersal for other species, and glaciated and unglaciated areas. Clark’s nutcrackers also irruptively disperse long distances in search of pine seed crops, creating the potential for gene flow among populations. Using the highly variable mitochondrial DNA control region, seven microsatellite loci, and species distribution modeling, we examined the effects of glaciations and dispersal barriers on population genetic patterns and population structure of nutcrackers. We sequenced 900 bp of mitochondrial control region for 169 individuals from 15 populations and analysed seven polymorphic microsatellite loci for 13 populations across the Clark’s nutcracker range. We used species distribution modeling and a range of phylogeographic analyses to examine evolutionary history. Clark’s nutcracker populations are not highly differentiated throughout their range, suggesting high levels of gene flow among populations, though we did find some evidence of isolation by distance and peripheral isolation. Our analyses suggested expansion from a single refugium after the last glacial maximum, but patterns of genetic diversity and paleodistribution modeling of suitable habitat were inconclusive as to the location of this refugium. Potential barriers to dispersal (e.g. mountain ranges) do not appear to restrict gene flow in Clark’s nutcracker, and postglacial expansion likely occurred quickly from a single refugium located south of the ice sheets.