Climatic drivers of leaf traits and genetic divergence in the tree Annona crassiflora: a broad spatial survey in the Brazilian savannas

The Cerrado is the largest South American savanna and encompasses substantial species diversity and environmental variation. Nevertheless, little is known regarding the influence of the environment on population divergence of Cerrado species. Here, we searched for climatic drivers of genetic (nuclear microsatellites) and leaf trait divergence in Annona crassiflora, a widespread tree in the Cerrado. The sampling encompassed all phytogeographic provinces of the continuous area of the Cerrado and included 397 individuals belonging to 21 populations. Populations showed substantial genetic and leaf trait divergence across the species' range. Our data revealed three spatially defined genetic groups (eastern, western and southern) and two morphologically distinct groups (eastern and western only). The east‐west split in both the morphological and genetic data closely mirrors previously described phylogeographic patterns of Cerrado species. Generalized linear mixed effects models and multiple regression analyses revealed several climatic factors associated with both genetic and leaf trait divergence among populations of A. crassiflora. Isolation by environment (IBE) was mainly due to temperature seasonality and precipitation of the warmest quarter. Populations that experienced lower precipitation summers and hotter winters had heavier leaves and lower specific leaf area. The southwestern area of the Cerrado had the highest genetic diversity of A. crassiflora, suggesting that this region may have been climatically stable. Overall, we demonstrate that a combination of current climate and past climatic changes have shaped the population divergence and spatial structure of A. crassiflora. However, the genetic structure of A. crassiflora reflects the biogeographic history of the species more strongly than leaf traits, which are more related to current climate.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

Living in two worlds: Evolutionary mechanisms act differently in the native and introduced ranges of an invasive plant

Identifying the factors that influence spatial genetic structure among populations can provide insights into the evolution of invasive plants. In this study, we used the common reed (Phragmites australis), a grass native in Europe and invading North America, to examine the relative importance of geographic, environmental (represented by climate here), and human effects on population genetic structure and its changes during invasion. We collected samples of P. australis from both the invaded North American and native European ranges and used molecular markers to investigate the population genetic structure within and between ranges. We used path analysis to identify the contributions of each of the three factors—geographic, environmental, and human‐related—to the formation of spatial genetic patterns. Genetic differentiation was observed between the introduced and native populations, and their genetic structure in the native and introduced ranges was different. There were strong effects of geography and environment on the genetic structure of populations in the native range, but the human‐related factors manifested through colonization of anthropogenic habitats in the introduced range counteracted the effects of environment. The between‐range genetic differences among populations were mainly explained by the heterogeneous environment between the ranges, with the coefficient 2.6 times higher for the environment than that explained by the geographic distance. Human activities were the primary contributor to the genetic structure of the introduced populations. The significant environmental divergence between ranges and the strong contribution of human activities to the genetic structure in the introduced range suggest that invasive populations of P. australis have evolved to adapt to a different climate and to human‐made habitats in North America.     

Ecological connectivity assessment in a strongly structured fire salamander (Salamandra salamandra) population

Small populations are more prone to extinction if the dispersal among them is not adequately maintained by ecological connections. The degree of isolation between populations could be evaluated measuring their genetic distance, which depends on the respective geographic (isolation by distance, IBD) and/or ecological (isolation by resistance, IBR) distances. The aim of this study was to assess the ecological connectivity of fire salamander Salamandra salamandra populations by means of a landscape genetic approach. The species lives in broad‐leaved forest ecosystems and is particularly affected by fragmentation due to its habitat selectivity and low dispersal capability. We analyzed 477 biological samples collected in 47 sampling locations (SLs) in the mainly continuous populations of the Prealpine and Eastern foothill lowland (PEF) and 10 SLs in the fragmented populations of the Western foothill (WF) lowland of Lombardy (northern Italy). Pairwise genetic distances (Chord distance, D_C) were estimated from allele frequencies of 16 microsatellites loci. Ecological distances were calculated using one of the most promising methodology in landscape genetics studies, the circuit theory, applied to habitat suitability maps. We realized two habitat suitability models: one without barriers (EcoD) and a second one accounting for the possible barrier effect of main roads (EcoD_b). Mantel tests between distance matrices highlighted how the Log‐D_C in PEF populations was related to log‐transformed geographic distance (confirming a prevalence of IBD), while it was explained by the Log‐EcoD, and particularly by the Log‐EcoD_b, in WF populations, even when accounting for the confounding effect of geographic distance (highlighting a prevalence of IBR). Moreover, we also demonstrated how considering the overall population, the effect of Euclidean or ecological distances on genetic distances acting at the level of a single group (PEF or WF populations) could not be detected, when population are strongly structured.     

Characterizing genic and nongenic molecular markers: comparison of microsatellites and SNPs

The implications of transitioning to single nucleotide polymorphism (SNPs) from microsatellite markers (MSs) have been investigated in a number of population genetics studies, but the effect of genomic location on the amount of information each type of marker reveals has not been explored in detail. We developed novel SNP markers flanking 1 kb regions of 13 genic (within gene or <1 kb away from gene) and 13 nongenic (>10 kb from annotated gene) MSs in the threespine stickleback genome to obtain comparable data for both types of markers. We analysed patterns of genetic diversity and divergence on various geographic scales after converting the SNP loci within each genomic region into haplotypes. Marker type (SNP haplotype or MS) and location (genic or nongenic) significantly affected most estimates of population diversity and divergence. Between‐lineage divergence was significantly higher in SNP haplotypes (genic and nongenic), however, within‐lineage divergence was similar between marker types. Most divergence and diversity measures were uncorrelated between markers, except for population differentiation which was correlated between MSs and SNP haplotypes (both genic and nongenic). Broad‐scale population structure and assignment were similarly resolved by both marker types, however, only the MSs were able to delimit fine‐scale population structuring, particularly when genic and nongenic markers were combined. These results demonstrate that estimates of genetic variability and differentiation among populations can be strongly influenced by marker type, their genomic location in relation to genes and by the interaction of these two factors. This highlights the importance of having an awareness of the inherent strengths and limitations associated with different molecular tools to select the most appropriate methods for accurately addressing various ecological and evolutionary questions.     

Genetic structuring in a Neotropical palm analyzed through an Andean orogenesis‐scenario

Andean orogenesis has driven the development of very high plant diversity in the Neotropics through its impact on landscape evolution and climate. The analysis of the intraspecific patterns of genetic structure in plants would permit inferring the effects of Andean uplift on the evolution and diversification of Neotropical flora. In this study, using microsatellite markers and Bayesian clustering analyses, we report the presence of four genetic clusters for the palm Oenocarpus bataua var. bataua which are located within four biogeographic regions in northwestern South America: (a) Chocó rain forest, (b) Amotape‐Huancabamba Zone, (c) northwestern Amazonian rain forest, and (d) southwestern Amazonian rain forest. We hypothesize that these clusters developed following three genetic diversification events mainly promoted by Andean orogenic events. Additionally, the distinct current climate dynamics among northwestern and southwestern Amazonia may maintain the genetic diversification detected in the western Amazon basin. Genetic exchange was identified between the clusters, including across the Andes region, discarding the possibility of any cluster to diversify as a distinct intraspecific variety. We identified a hot spot of genetic diversity in the northern Peruvian Amazon around the locality of Iquitos. We also detected a decrease in diversity with distance from this area in westward and southward direction within the Amazon basin and the eastern Andean foothills. Additionally, we confirmed the existence and divergence of O. bataua var. bataua from var. oligocarpus in northern South America, possibly expanding the distributional range of the latter variety beyond eastern Venezuela, to the central and eastern Andean cordilleras of Colombia. Based on our results, we suggest that Andean orogenesis is the main driver of genetic structuring and diversification in O. bataua within northwestern South America.     

Exploring the role of Micronesian islands in the maintenance of coral genetic diversity in the Pacific Ocean

Understanding how genetic diversity is maintained across patchy marine environments remains a fundamental problem in marine biology. The Coral Triangle, located in the Indo‐West Pacific, is the centre of marine biodiversity and has been proposed as an important source of genetic diversity for remote Pacific reefs. Several studies highlight Micronesia, a scattering of hundreds of small islands situated within the North Equatorial Counter Current, as a potentially important migration corridor. To test this hypothesis, we characterized the population genetic structure of two ecologically important congeneric species of reef‐building corals across greater Micronesia, from Palau to the Marshall Islands. Genetic divergences between islands followed an isolation‐by‐distance pattern, with Acropora hyacinthus exhibiting greater genetic divergences than A. digitifera, suggesting different migration capabilities or different effective population sizes for these closely related species. We inferred dispersal distance using a biophysical larval transport model, which explained an additional 15–21% of the observed genetic variation compared to between‐island geographical distance alone. For both species, genetic divergence accumulates and genetic diversity diminishes with distance from the Coral Triangle, supporting the hypothesis that Micronesian islands act as important stepping stones connecting the central Pacific with the species‐rich Coral Triangle. However, for A. hyacinthus, the species with lower genetic connectivity, immigration from the subequatorial Pacific begins to play a larger role in shaping diversity than input from the Coral Triangle. This work highlights the enormous dispersal potential of broadcast‐spawning corals and identifies the biological and physical drivers that influence coral genetic diversity on a regional scale.     

Spatial population genetic structure of a bacterial parasite in close coevolution with its host

Knowledge of a species’ population genetic structure can provide insight into fundamental ecological and evolutionary processes including gene flow, genetic drift and adaptive evolution. Such inference is of particular importance for parasites, as an understanding of their population structure can illuminate epidemiological and coevolutionary dynamics. Here, we describe the population genetic structure of the bacterium Pasteuria ramosa, a parasite that infects planktonic crustaceans of the genus Daphnia. This system has become a model for investigations of host–parasite interactions and represents an example of coevolution via negative frequency‐dependent selection (aka “Red Queen” dynamics). To sample P. ramosa, we experimentally infected a panel of Daphnia hosts with natural spore banks from the sediments of 25 ponds throughout much of the species range in Europe and western Asia. Using 12 polymorphic variable number tandem repeat loci (VNTR loci), we identified substantial genetic diversity, both within and among localities, that was structured geographically among ponds. Genetic diversity was also structured among host genotypes within ponds, although this pattern varied by locality, with P. ramosa at some localities partitioned into distinct host‐specific lineages, and other localities where recombination had shuffled genetic variation among different infection phenotypes. Across the sample range, there was a pattern of isolation by distance, and principal components analysis coupled with Procrustes rotation identified congruence between patterns of genetic variation and geography. Our findings support the hypothesis that Pasteuria is an endemic parasite coevolving closely with its host. These results provide important context for previous studies of this model system and inform hypotheses for future research.     

Increased genetic divergence between two closely related fir species in areas of range overlap

Because of introgressive hybridization, closely related species can be more similar to each other in areas of range overlap (parapatry or sympatry) than in areas where they are geographically isolated from each other (allopatry). Here, we report the reverse situation based on nuclear genetic divergence between two fir species, Abies chensiensis and Abies fargesii, in China, at sites where they are parapatric relative to where they are allopatric. We examined genetic divergence across 126 amplified fragment length polymorphism (AFLP) markers in a set of 172 individuals sampled from both allopatric and parapatric populations of the two species. Our analyses demonstrated that AFLP divergence was much greater between the species when comparisons were made between parapatric populations than between allopatric populations. We suggest that selection in parapatry may have largely contributed to this increased divergence.     

The role of geography and ecology in shaping repeated patterns of morphological and genetic differentiation between European minnows (Phoxinus phoxinus) from the Pyrenees and the Alps

Neutral and selective processes can drive repeated patterns of evolution in different groups of populations experiencing similar ecological gradients. In this paper, we used a combination of nuclear and mitochondrial DNA markers, as well as geometric morphometrics, to investigate repeated patterns of morphological and genetic divergence of European minnows in two mountain ranges: the Pyrenees and the Alps. European minnows (Phoxinus phoxinus) are cyprinid fish inhabiting most freshwater bodies in Europe, including those in different mountain ranges that could act as major geographical barriers to gene flow. We explored patterns of P. phoxinus phenotypic and genetic diversification along a gradient of altitude common to the two mountain ranges, and tested for isolation by distance (IBD), isolation by environment (IBE) and isolation by adaptation (IBA). The results indicated that populations from the Pyrenees and the Alps belong to two well differentiated, reciprocally monophyletic mtDNA lineages. Substantial genetic differentiation due to geographical isolation within and between populations from the Pyrenees and the Alps was also found using rapidly evolving AFLPs markers (isolation by distance or IBD), as well as morphological differences between mountain ranges. Also, morphology varied strongly with elevation and so did genetic differentiation to a lower extent. Despite moderate evidence for IBE and IBA, and therefore of repeated evolution, substantial population heterogeneity was found at the genetic level, suggesting that selection and population specific genetic drift act in concert to affect genetic divergence.     

Intercontinental long‐distance seed dispersal across the Mediterranean Basin explains population genetic structure of a bird‐dispersed shrub

Long‐distance dispersal (LDD) is a pivotal process for plants determining their range of distribution and promoting gene flow among distant populations. Most fleshy‐fruited species rely on frugivorous vertebrates to disperse their seeds across the landscape. While LDD events are difficult to record, a few ecological studies have shown that birds move a sizeable number of ingested seeds across geographic barriers, such as sea straits. The foraging movements of migrant frugivores across distant populations, including those separated by geographic barriers, creates a constant flow of propagules that in turn shapes the spatial distributions of the genetic variation in populations. Here, we have analysed the genetic diversity and structure of 74 populations of Pistacia lentiscus, a fleshy‐fruited shrub widely distributed in the Mediterranean Basin, to elucidate whether the Mediterranean Sea acts as a geographic barrier or alternatively whether migratory frugivorous birds promote gene flow among populations located on both sides of the sea. Our results show reduced genetic distances among populations, including intercontinental populations, and they show a significant genetic structure across an eastern‐western axis. These findings are consistent with known bird migratory routes that connect the European and African continents following a north‐southwards direction during the fruiting season of many fleshy‐fruited plants. Further, approximate Bayesian analysis failed to explain the observed patterns as a result of historical population migrations at the end of Last Glacial Maximum. Therefore, anthropic and/or climatic changes that would disrupt the migratory routes of frugivorous birds might have genetic consequences for the plant species they feed upon.     

Ecological adaptation and reproductive isolation in sympatry: genetic and phenotypic evidence for native host races of Rhagoletis pomonella

Ecological speciation with gene flow may be an important mode of diversification for phytophagous insects. The recent shift of Rhagoletis pomonella from its native host downy hawthorn (Crataegus mollis) to introduced apple (Malus domestica) in the northeastern United States is a classic example of sympatric host race formation. Here, we test whether R. pomonella has similarly formed host races on four native Crataegus species in the southern United States: western mayhaw (C. opaca), blueberry hawthorn (C. brachyacantha), southern red hawthorn (C. mollis var. texana) and green hawthorn (C. viridis). These four southern hosts differ from each other in their fruiting phenology and in the volatile compounds emitted from the surface of their fruits. These two traits form the basis of ecological reproductive isolation between downy hawthorn and apple flies in the north. We report evidence from microsatellite population surveys and eclosion studies supporting the existence of genetically differentiated and partially reproductively isolated host races of southern hawthorn flies. The results provide an example of host shifting and ecological divergence involving native plants and imply that speciation with gene flow may be commonly initiated in Rhagoletis when ecological opportunity presents itself.     

Evolutionary insights into Rhinolophus episcopus (Chiroptera,Rhinolophidae) in China: Isolation by distance,environment, or sensory system?

Evolutionary processes can be influenced by several factors, such as geographic isolation, environmental selection, and sensory variation. For most nocturnal bats, echolocation is the primary sensory system used to prey and communicate, and plays important roles in chiropteran diversification and evolution. Understanding the relative contribution of geography, the environment, and this sensory system to population genetic divergence can elucidate the processes involved in bat incipient speciation and evolution. In this study, we collected spatial and environmental information, echolocation calls, as well as the previously published genetic data (six microsatellite loci and the mitochondrial cytochrome b gene) of widely distributed Rhinolophus episcopus populations to test three hypotheses for nuclear and mitochondrial divergence (isolation by distance, isolation by environment, and isolation by sensory variation) and unveil the factors that drive intraspecific genetic differentiation. The moderate level of nuclear differentiation was correlated with geographic/spatial distance and acoustic variation, whereas the relatively high level of mitochondrial differentiation was mainly associated with acoustic divergence. No significant correlation was observed between genetic divergence and environmental variables. Among the three factors, acoustic divergence explained the highest percentage of both nuclear and mitochondrial divergence. Thus, our results indicate that sensory variation may have played important roles in driving population isolation early in bat speciation, which is consistent with the hypothesis of isolation by sensory variation. Our study emphasizes the need to consider more factors, especially sensory traits, and combine multiple statistical methods in landscape genetic studies to test their potential contributions to driving population divergence.     

The challenge of separating signatures of local adaptation from those of isolation by distance and colonization history: The case of two white pines

Accurately detecting signatures of local adaptation using genetic‐environment associations (GEAs) requires controlling for neutral patterns of population structure to reduce the risk of false positives. However, a high degree of collinearity between climatic gradients and neutral population structure can greatly reduce power, and the performance of GEA methods in such case is rarely evaluated in empirical studies. In this study, we attempted to disentangle the effects of local adaptation and isolation by environment (IBE) from those of isolation by distance (IBD) and isolation by colonization from glacial refugia (IBC) using range‐wide samples in two white pine species. For this, SNPs from 168 genes, including 52 candidate genes for growth and phenology, were genotyped in 133 and 61 populations of Pinus strobus and P. monticola, respectively. For P. strobus and using all 153 SNPs, climate (IBE) did not significantly explained among‐population variation when controlling for IBD and IBC in redundancy analyses (RDAs). However, 26 SNPs were significantly associated with climate in single‐locus GEA analyses (Bayenv2 and LFMM), suggesting that local adaptation took place in the presence of high gene flow. For P. monticola, we found no evidence of IBE using RDAs and weaker signatures of local adaptation using GEA and F_ST outlier tests, consistent with adaptation via phenotypic plasticity. In both species, the majority of the explained among‐population variation (69 to 96%) could not be partitioned between the effects of IBE, IBD, and IBC. GEA methods can account differently for this confounded variation, and this could explain the small overlap of SNPs detected between Bayenv2 and LFMM. Our study illustrates the inherent difficulty of taking into account neutral structure in natural populations and the importance of sampling designs that maximize climatic variation, while minimizing collinearity between climatic gradients and neutral structure.     

Patterns of reproductive isolation within and between two Lygaeus species characterized by sexual conflicts over mating

Theory suggests that, under some circumstances, sexual conflict over mating can lead to divergent sexually antagonistic coevolution among populations for traits associated with mating, and that this can promote reproductive isolation and hence speciation. However, sexual conflict over mating may also select for traits (e.g. male willingness to mate) that enhance gene flow between populations, limiting population divergence. In the present study, we compare pre‐ and post‐mating isolation within and between two species characterized by male–female conflict over mating rate. We quantify sexual isolation among five populations of the seed bug Lygaeus equestris collected from Italy and Sweden, and two replicates of a population of the sister‐species Lygaeus simulans, also collected from Italy. We find no evidence of reproductive isolation amongst populations of L. equestris, suggesting that sexual conflict over mating has not led to population divergence in relevant mating traits in L. equestris. However, there was strong asymmetric pre‐mating isolation between L. equestris and L. simulans: male L. simulans were able to mate successfully with female L. equestris, whereas male L. equestris were largely unable to mate with female L. simulans. We found little evidence for strong post‐mating isolation between the two species, however, with hybrid F₂ offspring being produced. Our results suggest that sexual conflict over mating has not led to population divergence, and indeed perhaps supports the contrary theoretical prediction that male willingness to mate may retard speciation by promoting gene flow.     

Postglacial dispersal patterns and mitochondrial genetic structure of the Pyrenean desman (Galemys pyrenaicus) in the northwestern region of the Iberian Peninsula

The genetic structure of small semiaquatic animals may be influenced by dispersal across both rivers and land. The relative importance of these two modes of dispersal may vary across different species and with ecological conditions and evolutionary periods. The Pyrenean desman (Galemys pyrenaicus) is an endemic mammal of the Iberian Peninsula with a strong phylogeographic structure and semiaquatic habits, thus making it an ideal model to study the effects of river and overland dispersal on its genetic structure. Thanks to different types of noninvasive samples, we obtained an extensive sampling of the Pyrenean desman from the northwestern region of the Iberian Peninsula and sequenced two mitochondrial DNA fragments. We then analyzed, using an isolation‐by‐distance approach, the correlation between phylogenetic distances and geographical distances measured along both river networks and land to infer the relative importance of river and overland dispersal. We found that the correlations in the whole area and in a large basin were consistent with an effect of overland dispersal, which may be due to the postglacial colonization of new territories using terrestrial corridors and, possibly, a more extensive fluvial network that may have been present during the Holocene. However, in a small basin, likely to be less influenced by the impact of ancient postglacial dispersal, the correlations suggested significant overall effects of both overland and river dispersal, as expected for a semiaquatic mammal. Therefore, different scales and geographical regions reflect different aspects of the evolutionary history and ecology of this semiaquatic species using this isolation‐by‐distance method. The results we obtained may have crucial implications for the conservation of the Pyrenean desman because they reinforce the importance of interbasin dispersal for this species in the studied area and the need to protect the whole riverine ecosystem, including rivers, upland streams and terrestrial corridors between basins.     

The influence of landscape configuration and environment on population genetic structure in a sedentary passerine: insights from loci located in different genomic regions

The study of the factors structuring genetic variation can help to infer the neutral and adaptive processes shaping the demographic and evolutionary trajectories of natural populations. Here, we analyse the role of isolation by distance (IBD), isolation by resistance (IBR, defined by landscape composition) and isolation by environment (IBE, estimated as habitat and elevation dissimilarity) in structuring genetic variation in 25 blue tit (Cyanistes caeruleus) populations. We typed 1385 individuals at 26 microsatellite loci classified into two groups by considering whether they are located into genomic regions that are actively (TL; 12 loci) or not (NTL; 14 loci) transcribed to RNA. Population genetic differentiation was mostly detected using the panel of NTL. Landscape genetic analyses showed a pattern of IBD for all loci and the panel of NTL, but genetic differentiation estimated at TL was only explained by IBR models considering high resistance for natural vegetation and low resistance for agricultural lands. Finally, the absence for IBE suggests a lack of divergent selection pressures associated with differences in habitat and elevation. Overall, our study shows that markers located in different genomic regions can yield contrasting inferences on landscape‐level patterns of realized gene flow in natural populations.     

Range‐wide analysis of genetic structure in a widespread,highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad‐scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analysed 10 microsatellite loci from over 1900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual‐ and population‐based analyses, we found evidence for three main genetic lineages, one corresponding to the ‘mule deer’ morphological type and two to the ‘black‐tailed deer’ type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles, and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation‐by‐distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation and management.