Revisiting the Iberian honey bee (Apis mellifera iberiensis) contact zone: maternal and genome‐wide nuclear variations provide support for secondary contact from historical refugia

Dissecting diversity patterns of organisms endemic to Iberia has been truly challenging for a variety of taxa, and the Iberian honey bee is no exception. Surveys of genetic variation in the Iberian honey bee are among the most extensive for any honey bee subspecies. From these, differential and complex patterns of diversity have emerged, which have yet to be fully resolved. Here, we used a genome‐wide data set of 309 neutrally tested single nucleotide polymorphisms (SNPs), scattered across the 16 honey bee chromosomes, which were genotyped in 711 haploid males. These SNPs were analysed along with an intergenic locus of the mtDNA, to reveal historical patterns of population structure across the entire range of the Iberian honey bee. Overall, patterns of population structure inferred from nuclear loci by multiple clustering approaches and geographic cline analysis were consistent with two major clusters forming a well‐defined cline that bisects Iberia along a northeastern–southwestern axis, a pattern that remarkably parallels that of the mtDNA. While a mechanism of primary intergradation or isolation by distance could explain the observed clinal variation, our results are more consistent with an alternative model of secondary contact between divergent populations previously isolated in glacial refugia, as proposed for a growing list of other Iberian taxa. Despite current intense honey bee management, human‐mediated processes have seemingly played a minor role in shaping Iberian honey bee genetic structure. This study highlights the complexity of the Iberian honey bee patterns and reinforces the importance of Iberia as a reservoir of Apis mellifera diversity.     

Molecular biogeography and host relations of a parasitoid fly

Successful geographic range expansion by parasites and parasitoids may also require host range expansion. Thus, the evolutionary advantages of host specialization may trade off against the ability to exploit new host species encountered in new geographic regions. Here, we use molecular techniques and confirmed host records to examine biogeography, population divergence, and host flexibility of the parasitoid fly, Ormia ochracea (Bigot). Gravid females of this fly find their cricket hosts acoustically by eavesdropping on male cricket calling songs; these songs vary greatly among the known host species of crickets. Using both nuclear and mitochondrial genetic markers, we (a) describe the geographical distribution and subdivision of genetic variation in O. ochracea from across the continental United States, the Mexican states of Sonora and Oaxaca, and populations introduced to Hawaii; (b) demonstrate that the distribution of genetic variation among fly populations is consistent with a single widespread species with regional host specialization, rather than locally differentiated cryptic species; (c) identify the more‐probable source populations for the flies introduced to the Hawaiian islands; (d) examine genetic variation and substructure within Hawaii; (e) show that among‐population geographic, genetic, and host song distances are all correlated; and (f) discuss specialization and lability in host‐finding behavior in light of the diversity of cricket songs serving as host cues in different geographically separate populations.     

Increased genetic structuring of isolated Salamandra salamandra populations (Caudata: Salamandridae) at the margins of the Carpathian Mountains

The traditional southern Pleistocene refugia hypothesis in Europe has lately been challenged for several animal and plant species. The Carpathian Basin, especially at the marginal regions, is one of the recently recognized biodiversity hotspots in Europe. Marginal populations are prone to have lower genetic diversity and higher genetic differentiation than central populations. Here, we examined one mitochondrial DNA fragment (D‐loop) and nine nuclear (microsatellite) loci to describe the genetic diversity and phylogeographical pattern of fire salamander (Salamandra salamandra) populations in the Carpathian Basin with focusing on the southern margins of the Western Carpathians, where isolated populations of this species are present. Analyses of microsatellites indicated reduced genetic diversity for most of the isolated populations. Based on the mitochondrial DNA, only two haplotypes were found, whereas the analyses with the nuclear markers revealed a more recent genetic split between Western (Alpine) and Eastern (Carpathian) populations, and separated the Apuseni Mountains population (part of the Western Carpathians). Using approximate Bayesian computation analyses, we identified the most probable colonization scenario for the isolated North Hungarian Carpathian Basin populations. The split between isolated salamander populations from the central populations in the Carpathian Mountains dates back to the beginning of the Late Pleistocene, while the split between most of the Hungarian populations can be associated with the Last Glacial Maximum. We found evidence for long‐time isolation between the marginal Carpathian Basin and central populations. Our results also show that S. salamandra survived glacial periods in the temperate forests of north‐east Pannonia (North Hungarian Mountains), confirming that the Carpathian Basin served as important northerly refugia during the Pleistocene climatic oscillations.     

Environmental gradients shape the genetic structure of two medicinal Salvia species in Jordan

• Environmental gradients, and particularly climatic variables, exert a strong influence on plant distribution and, potentially, population genetic diversity and differentiation. Differences in water availability can cause among‐population variation in ecological processes and can thus interrupt populations’ connectivity and isolate them environmentally. The present study examines the effect of environmental heterogeneity on plant populations due to environmental isolation unrelated to geographic distance.
• Using AFLP markers, we analyzed genetic diversity and differentiation among 12 Salvia spinosa populations and 13 Salvia syriaca populations from three phytogeographical regions (Mediterranean, Irano‐Turanian and Saharo‐Arabian) representing the extent of the species’ geographic range in Jordan. Differences in geographic location and climate were considered in the analyses.
• For both species, flowering phenology varied among populations and regions. Irano‐Turanian and Saharo‐Arabian populations had higher genetic diversity than Mediterranean populations, and genetic diversity increased significantly with increasing temperature. Genetic diversity in Salvia syriaca was affected by population size, while genetic diversity responded to drought in S. spinosa. For both species, high levels of genetic differentiation were found as well as two well‐supported phytogeographical groups of populations, with Mediterranean populations clustering in one group and the Irano‐Turanian and Saharo‐Arabian populations in another. Genetic distance was significantly correlated to environmental distance, but not to geographic distance.
• Our data indicate that populations from moist vs. arid environments are environmentally isolated, where environmental gradients affect their flowering phenology, limit gene flow and shape their genetic structure. We conclude that environmental heterogeneity may act as driver for the observed variation in genetic diversity.

     

Pathways of spread of the introduced ascidian Styela clava (Tunicata) in Northern Europe, as revealed by microsatellite markers

Styela clava, a solitary ascidian native to the NW Pacific, has become a conspicuous member of fouling communities in NW European waters. As its natural dispersal appears to be limited, the wide distribution of S. clava along coasts within its introduced range may be attributed to secondary spread assisted by human activities. Here, we used six microsatellite loci to examine the genetic diversity and extent of gene flow among S. clava populations in its European introduced range. Samples were collected from 21 populations within Europe (N = 808), 4 populations within the USA and two populations within the native range (Japan). Large variation in genetic diversity was observed among the European populations but were not explained either by the geographic distance from the first introduction area (i.e. Plymouth, UK) nor by the time elapsed since the introduction. No founder effect was observed in the introduced populations, except possibly in Puget Sound (USA). At least two different introductions occurred in Europe, identified as distinct genetic clusters: northern Danish populations (resembling one Japanese population), and the rest of Europe; a sample from Shoreham (England) possibly represents a third introduction. In North America, the population from the Atlantic was genetically similar to the majority of European populations, suggesting a European origin for populations on this seaboard, while populations from the Pacific coast were genetically similar to the same Japanese population as the Danish populations.     

Nesting habits influence population genetic structure of a bee living in anthropogenic disturbance

While most organisms are negatively affected by anthropogenic disturbance, a few species thrive in landscapes altered by humans. Typically, native bees are negatively impacted by anthropogenic environmental change, including habitat alteration and climate change. Here, we investigate the population structure of the eastern carpenter bee Xylocopa virginica, a generalist pollinator with a broad geographic range spanning eastern North America. Eastern carpenter bees now nest almost exclusively in artificial wooden structures, linking their geographic distribution and population structure to human activities and disturbance. To investigate the population structure of these bees, we sampled females from 16 different populations from across their range. Nine species‐specific microsatellite loci showed that almost all populations are genetically distinct, but with high levels of genetic diversity and low levels of inbreeding overall. Broadly speaking, populations clustered into three distinct genetic groups: a northern group, a western group and a core group. The northern group had low effective population sizes, decreased genetic variability and the highest levels of inbreeding in the data set, suggesting that carpenter bees may be expanding their range northward. The western group was genetically distinct, but lacked signals of a recent range expansion. Climatic data showed that summer and winter temperatures explained a significant amount of the genetic differentiation seen among populations, while precipitation did not. Our results indicate that X. virginica may be one of the rare ‘anthrophilic’ species that thrive in the face of anthropogenic disturbance.     

Glacial refugia and postglacial expansion of the alpine–prealpine plant species Polygala chamaebuxus

The shrubby milkwort (Polygala chamaebuxus L.) is widely distributed in the Alps, but occurs also in the lower mountain ranges of Central Europe such as the Franconian Jura or the Bohemian uplands. Populations in these regions may either originate from glacial survival or from postglacial recolonization. In this study, we analyzed 30 populations of P. chamaebuxus from the whole distribution range using AFLP (Amplified Fragment Length Polymorphism) analysis to identify glacial refugia and to illuminate the origin of P. chamaebuxus in the lower mountain ranges of Central Europe. Genetic variation and the number of rare fragments within populations were highest in populations from the central part of the distribution range, especially in the Southern Alps (from the Tessin Alps and the Prealps of Lugano to the Triglav Massiv) and in the middle part of the northern Alps. These regions may have served, in accordance with previous studies, as long‐term refugia for the glacial survival of the species. The geographic pattern of genetic variation, as revealed by analysis of molecular variance, Bayesian cluster analysis and a PopGraph genetic network was, however, only weak. Instead of postglacial recolonization from only few long‐term refugia, which would have resulted in deeper genetic splits within the data set, broad waves of postglacial expansion from several short‐term isolated populations in the center to the actual periphery of the distribution range seem to be the scenario explaining the observed pattern of genetic variation most likely. The populations from the lower mountain ranges in Central Europe were more closely related to the populations from the southwestern and northern than from the nearby eastern Alps. Although glacial survival in the Bohemian uplands cannot fully be excluded, P. chamaebuxus seems to have immigrated postglacially from the southwestern or central‐northern parts of the Alps into these regions during the expansion of the pine forests in the early Holocene.     

Population genetic structure and demographic history of the lone star tick,Amblyomma americanum (Ixodida: Ixodidae): New evidence supporting old records

Range expansions are a potential outcome of changes in habitat suitability, which commonly result as a consequence of climate change. Hypotheses on such changes in the geographic distribution of a certain species can be evaluated using population genetic structure and demography. In this study we explore the population genetic structure, genetic variability, demographic history of, and habitat suitability for Amblyomma americanum, a North American tick species that is a known vector of several pathogenic microorganisms. We used a double digestion restriction site‐associated DNA sequencing technique (dd‐RAD seq) and discovered 8,181 independent single nucleotide polymorphisms (SNPs) in 189 ticks from across the geographic range of the species. Genetic diversity was low, particularly when considering the broad geographic range of this species. The edge populations were less diverse than populations belonging to the historic range, possibly indicative of a range expansion, but this hypothesis was not statistically supported by a test based on genetic data. Nonetheless, moderate levels of population structure and substructure were detected between geographic regions. For New England, demographic and species distribution models support a scenario where A. americanum was present in more northern locations in the past, underwent a bottleneck, and subsequently recovered. These results are consistent with a hypothesis that this species is re‐establishing in this area, rather than one focused on range expansion from the south. This hypothesis is consistent with old records describing the presence of A. americanum in the northeastern US in the early colonial period.     

The phylogeographical clade trade: tracing the impact of human-mediated dispersal on the colonization of northern Europe by the oak gallwasp Andricus kollari

Human dispersal of organisms is an important process modifying natural patterns of biodiversity. Such dispersal generates new patterns of genetic diversity that overlie natural phylogeographical signatures, allowing discrimination between alternative dispersal mechanisms. Here we use allele frequency and DNA sequence data to distinguish between alternative scenarios (unassisted range expansion and long range introduction) for the colonization of northern Europe by an oak-feeding gallwasp, Andricus kollari. Native to Mediterranean latitudes from Portugal to Iran, this species became established in northern Europe following human introduction of a host plant, the Turkey oak Quercus cerris. Colonization of northern Europe is possible through three alternative routes: (i) unassisted range expansion from natural populations in the Iberian Peninsula; (ii) unassisted range expansion from natural populations in Italy and Hungary; or (iii) descent from populations imported to the UK as trade goods from the eastern Mediterranean in the 1830s. We show that while populations in France were colonized from sources in Italy and Hungary, populations in the UK and neighbouring parts of coastal northern Europe encompass allozyme and sequence variation absent from the known native range. Further, these populations show demographic signatures expected for large stable populations, rather than signatures of rapid population growth from small numbers of founders. The extent and spatial distribution of genetic diversity in the UK suggests that these A. kollari populations are derived from introductions of large numbers of individuals from each of two genetically divergent centres of diversity in the eastern Mediterranean. The strong spatial patterning in genetic diversity observed between different regions of northern Europe, and between sites in the UK, is compatible with leptokurtic models of population establishment.     

Reproductive isolation and environmental adaptation shape the phylogeography of mountain pine beetle (Dendroctonus ponderosae)

Chromosomal rearrangement can be an important mechanism driving population differentiation and incipient speciation. In the mountain pine beetle (MPB, Dendroctonus ponderosae), deletions on the Y chromosome that are polymorphic among populations are associated with reproductive incompatibility. Here, we used RAD sequencing across the entire MPB range in western North America to reveal the extent of the phylogeographic differences between Y haplotypes compared to autosomal and X‐linked loci. Clustering and geneflow analyses revealed three distinct Y haplogroups geographically positioned within and on either side of the Great Basin Desert. Despite close geographic proximity between populations on the boundaries of each Y haplogroup, there was extremely low Y haplogroup mixing among populations, and gene flow on the autosomes was reduced across Y haplogroup boundaries. These results are consistent with a previous study suggesting that independent degradation of a recently evolved neo‐Y chromosome in previously isolated populations causes male sterility or inviability among Y haplotype lineages. Phylogeographic results supported historic contraction of MPB into three separate Pleistocene glacial refugia followed by postglacial range expansion and secondary contact. Distinct sets of SNPs were statistically associated with environmental data among the most genetically distinct sets of geographic populations. This finding suggests that the process of adaptation to local climatic conditions is influenced by population genetic structure, with evidence for largely independent evolution in the most genetically isolated Y haplogroup.     

Population decline is accompanied by loss of genetic diversity in the Lesser Grey Shrike Lanius minor

The Lesser Grey Shrike has suffered successive declines in population size and a marked contraction of its breeding range since the early 20th century, largely because of long‐term agricultural intensification. This has resulted in a severely fragmented distribution in Western Europe, with isolated breeding nuclei in Spain, France and Italy and a more continuous distribution in Eastern Europe and Asia. Using a combination of nuclear and mitochondrial markers, we assessed the genetic structure and diversity of Lesser Grey Shrike populations from Western Europe, Central Europe and Asia. There was significant genetic differentiation among three major regional groups, one European and two Asian. Genetic diversity measures were lowest in the smallest and most marginal Spanish population. Limited genetic diversity, combined with rapid population decline, suggests the Spanish population may face extinction in the near future.     

Demographic history of a recent invasion of house mice on the isolated Island of Gough

Island populations provide natural laboratories for studying key contributors to evolutionary change, including natural selection, population size and the colonization of new environments. The demographic histories of island populations can be reconstructed from patterns of genetic diversity. House mice (Mus musculus) inhabit islands throughout the globe, making them an attractive system for studying island colonization from a genetic perspective. Gough Island, in the central South Atlantic Ocean, is one of the remotest islands in the world. House mice were introduced to Gough Island by sealers during the 19th century and display unusual phenotypes, including exceptionally large body size and carnivorous feeding behaviour. We describe genetic variation in Gough Island mice using mitochondrial sequences, nuclear sequences and microsatellites. Phylogenetic analysis of mitochondrial sequences suggested that Gough Island mice belong to Mus musculus domesticus, with the maternal lineage possibly originating in England or France. Cluster analyses of microsatellites revealed genetic membership for Gough Island mice in multiple coastal populations in Western Europe, suggesting admixed ancestry. Gough Island mice showed substantial reductions in mitochondrial and nuclear sequence variation and weak reductions in microsatellite diversity compared with Western European populations, consistent with a population bottleneck. Approximate Bayesian computation (ABC) estimated that mice recently colonized Gough Island (~100 years ago) and experienced a 98% reduction in population size followed by a rapid expansion. Our results indicate that the unusual phenotypes of Gough Island mice evolved rapidly, positioning these mice as useful models for understanding rapid phenotypic evolution.     

Phylogeography of the socially polymorphic sweat bee Halictus rubicundus (Hymenoptera: Halictidae)

The evolution of sociality in insects holds a central place in evolutionary theory. By examining the phylogenetic patterns of solitary and social behavior and how they correlate with ecological variables, we may identify factors important in the evolution of sociality. In this study, we investigated historical and biogeographical patterns of sociality in a socially polymorphic bee species (one that demonstrates both social and solitary nesting behavior). This unique system allows for a more powerful examination of evolutionary transitions in sociality than interspecific studies of obligately social and solitary species. We conducted a phylogenetic analysis among populations of the halictine bee Halictus rubicundus and then identified relationships among mitochondrial DNA sequence data, sociality, environmental conditions at the nesting site, and geographic location of populations of this species. Within North America, populations of H. rubicundus expressing social and solitary behavior belong to different genetic lineages. Sociality is also correlated with at least one environmental variable used in this study. Taken together, the results support the predictions for genetic control of sociality, but they are still consistent with social behavior at some level being determined by the environmental conditions at the nesting site.     

Investigating past range dynamics for a weed of cultivation,Silene vulgaris

Since the last glacial maximum (LGM), many plant and animal taxa have expanded their ranges by migration from glacial refugia. Weeds of cultivation may have followed this trend or spread globally following the expansion of agriculture or ruderal habitats associated with human‐mediated disturbance. We tested whether the range expansion of the weed Silene vulgaris across Europe fit the classical model of postglacial expansion from southern refugia, or followed known routes of the expansion of human agricultural practices. We used species distribution modeling to predict spatial patterns of postglacial expansion and contrasted these with the patterns of human agricultural expansion. A population genetic analysis using microsatellite loci was then used to test which scenario was better supported by spatial patterns of genetic diversity and structure. Genetic diversity was highest in southern Europe and declined with increasing latitude. Locations of ancestral demes from genetic cluster analysis were consistent with areas of predicted refugia. Species distribution models showed the most suitable habitat in the LGM on the southern coasts of Europe. These results support the typical postglacial northward colonization from southern refugia while refuting the east‐to‐west agricultural spread as the main mode of expansion for S. vulgaris. We know that S. vulgaris has recently colonized many regions (including North America and other continents) through human‐mediated dispersal, but there is no evidence for a direct link between the Neolithic expansion of agriculture and current patterns of genetic diversity of S. vulgaris in Europe. Therefore, the history of range expansion of S. vulgaris likely began with postglacial expansion after the LGM, followed by more recent global dispersal by humans.     

Patterns of genetic diversity in the polymorphic ground snake (Sonora semiannulata)

We evaluated the genetic diversity of a snake species with color polymorphism to understand the evolutionary processes that drive genetic structure across a large geographic region. Specifically, we analyzed genetic structure of the highly polymorphic ground snake, Sonora semiannulata, (1) among populations, (2) among color morphs (3) at regional and local spatial scales, using an amplified fragment length polymorphism dataset and multiple population genetic analyses, including F_ST-based and clustering analytical techniques. Based upon these methods, we found that there was moderate to low genetic structure among populations. However, this diversity was not associated with geographic locality at either spatial scale. Similarly, we found no evidence for genetic divergence among color morphs at either spatial scale. These results suggest that despite dramatic color polymorphism, this phenotypic diversity is not a major driver of genetic diversity within or among populations of ground snakes. We suggest that there are two mechanisms that could explain existing genetic diversity in ground snakes: recent range expansion from a genetically diverse founder population and current or recent gene flow among populations. Our findings have further implications for the types of color polymorphism that may generate genetic diversity in snakes.     

Phylogeographic insights into the invasion history and secondary spread of the signal crayfish in Japan

Successful invasion by nonindigenous species is often attributed to high propagule pressure, yet some foreign species become widespread despite showing reduced genetic variation due to founder effects. The signal crayfish (Pacifastacus leniusculus) is one such example, where rapid spread across Japan in recent decades is believed to be the result of only three founding populations. To infer the history and explore the success of this remarkable crayfish invasion, we combined detailed phylogeographical and morphological analyses conducted in both the introduced and native ranges. We sequenced 16S mitochondrial DNA of signal crayfish from across the introduced range in Japan (537 samples, 20 sites) and the native range in western North America (700 samples, 50 sites). Because chela size is often related to aggressive behavior in crayfish, and hence, their invasion success, we also measured chela size of a subset of specimens in both introduced and native ranges. Genetic diversity of introduced signal crayfish populations was as high as that of the dominant phylogeographic group in the native range, suggesting high propagule pressure during invasion. More recently established crayfish populations in Japan that originated through secondary spread from one of the founding populations exhibit reduced genetic diversity relative to older populations, probably as a result of founder effects. However, these newer populations also show larger chela size, consistent with expectations of rapid adaptations or phenotypic responses during the invasion process. Introduced signal crayfish populations in Japan originate from multiple source populations from a wide geographic range in the native range of western North America. A combination of high genetic diversity, especially for older populations in the invasive range, and rapid adaptation to colonization, manifested as larger chela in recent invasions, likely contribute to invasion success of signal crayfish in Japan.     

Mitochondrial and microsatellite DNA markers reveal a Balkan origin for the highly invasive horse-chestnut leaf miner Cameraria ohridella (Lepidoptera, Gracillariidae)

Biological invasions usually start with a small number of founder individuals. These founders are likely to represent a small fraction of the total genetic diversity found in the source population. Our study set out to trace genetically the geographical origin of the horse-chestnut leafminer, Cameraria ohridella , an invasive microlepidopteran whose area of origin is still unkown. Since its discovery in Macedonia 25 years ago, this insect has experienced an explosive westward range expansion, progressively colonizing all of Central and Western Europe. We used cytochrome oxidase I sequences (DNA barcode fragment) and a set of six polymorphic microsatellites to assess the genetic variability of C. ohridella populations, and to test the hypothesis that C. ohridella derives from the southern Balkans (Albania, Macedonia and Greece). Analysis of mtDNA of 486 individuals from 88 localities allowed us to identify 25 geographically structured haplotypes. In addition, 480 individuals from 16 populations from Europe and the southern Balkans were genotyped for 6 polymorphic microsatellite loci. High haplotype diversity and low measures of nucleotide diversities including a significantly negative Tajima's D indicate that C. ohridella has experienced rapid population expansion during its dispersal across Europe. Both mtDNA and microsatellites show a reduction in genetic diversity of C. ohridella populations sampled from artificial habitats (e.g. planted trees in public parks, gardens, along roads in urban or sub-urban areas) across Europe compared with C. ohridella sampled in natural stands of horse-chestnuts in the southern Balkans. These findings suggest that European populations of C. ohridella may indeed derive from the southern Balkans.     

Predicting the success of an invader: Niche shift versus niche conservatism

Invasive species can encounter environments different from their source populations, which may trigger rapid adaptive changes after introduction (niche shift hypothesis). To test this hypothesis, we investigated whether postintroduction evolution is correlated with contrasting environmental conditions between the European invasive and source ranges in the Asian tiger mosquito Aedes albopictus. The comparison of environmental niches occupied in European and source population ranges revealed more than 96% overlap between invasive and source niches, supporting niche conservatism. However, we found evidence for postintroduction genetic evolution by reanalyzing a published ddRADseq genomic dataset from 90 European invasive populations using genotype–environment association (GEA) methods and generalized dissimilarity modeling (GDM). Three loci, among which a putative heat‐shock protein, exhibited significant allelic turnover along the gradient of winter precipitation that could be associated with ongoing range expansion. Wing morphometric traits weakly correlated with environmental gradients within Europe, but wing size differed between invasive and source populations located in different climatic areas. Niche similarities between source and invasive ranges might have facilitated the establishment of populations. Nonetheless, we found evidence for environmental‐induced adaptive changes after introduction. The ability to rapidly evolve observed in invasive populations (genetic shift) together with a large proportion of unfilled potential suitable areas (80%) pave the way to further spread of Ae. albopictus in Europe.     

Genetic consequences of interglacial isolation in a steppe bird

In response to climate changes that have occurred during Pleistocene glacial cycles, taxa associated to steppe vegetation might have followed a pattern of historical evolution in which isolation and fragmentation of populations occurred during the short interglacials and expansion events occurred during the long glacial periods, in contrast to the pattern described for temperate species. Here, we use molecular genetic data to evaluate this idea in a steppe bird with Palaearctic distribution, the little bustard (Tetrax tetrax). Overall, extremely low genetic diversity and differentiation was observed among eight little bustard populations distributed in Spain and France. Mismatch distribution analyses showed that most little bustard populations expanded during cooling periods previous to, and just after, the last interglacial period (127,000-111,000 years before present), when steppe habitats were widespread across Europe. Coalescent-based methods suggested that glacial expansions have resulted in substantial admixture in Western Europe due to the existence of different interglacial refugia. Our results are consistent with a model of evolution and genetic consequences of Pleistocene cycles with low between-population genetic differentiation as a result of short-term isolation periods during interglacials and long-term exchange during glacial periods.