The biogeography of mitochondrial and nuclear discordance in animals

Combining nuclear (nuDNA) and mitochondrial DNA (mtDNA) markers has improved the power of molecular data to test phylogenetic and phylogeographic hypotheses and has highlighted the limitations of studies using only mtDNA markers. In fact, in the past decade, many conflicting geographic patterns between mitochondrial and nuclear genetic markers have been identified (i.e. mito-nuclear discordance). Our goals in this synthesis are to: (i) review known cases of mito-nuclear discordance in animal systems, (ii) to summarize the biogeographic patterns in each instance and (iii) to identify common drivers of discordance in various groups. In total, we identified 126 cases in animal systems with strong evidence of discordance between the biogeographic patterns obtained from mitochondrial DNA and those observed in the nuclear genome. In most cases, these patterns are attributed to adaptive introgression of mtDNA, demographic disparities and sex-biased asymmetries, with some studies also implicating hybrid zone movement, human introductions and Wolbachia infection in insects. We also discuss situations where divergent mtDNA clades seem to have arisen in the absence of geographic isolation. For those cases where foreign mtDNA haplotypes are found deep within the range of a second taxon, data suggest that those mtDNA haplotypes are more likely to be at a high frequency and are commonly driven by sex-biased asymmetries and/or adaptive introgression. In addition, we discuss the problems with inferring the processes causing discordance from biogeographic patterns that are common in many studies. In many cases, authors presented more than one explanation for discordant patterns in a given system, which indicates that likely more data are required. Ideally, to resolve this issue, we see important future work shifting focus from documenting the prevalence of mito-nuclear discordance towards testing hypotheses regarding the drivers of discordance. Indeed, there is great potential for certain cases of mitochondrial introgression to become important natural systems within which to test the effect of different mitochondrial genotypes on whole-animal phenotypes.     

Estimating the temporal and spatial extent of gene flow among sympatric lizard populations (genus Sceloporus) in the southern Mexican highlands

Interspecific gene flow is pervasive throughout the tree of life. Although detecting gene flow between populations has been facilitated by new analytical approaches, determining the timing and geography of hybridization has remained difficult, particularly for historical gene flow. A geographically explicit phylogenetic approach is needed to determine the overlap of ancestral populations. In this study, we performed population genetic analyses, species delimitation, simulations and a recently developed approach of species tree diffusion to infer the phylogeographic history, timing and geographic extent of gene flow in lizards of the Sceloporus spinosus group. The two species in this group, S. spinosus and S. horridus, are distributed in eastern and western portions of Mexico, respectively, but populations of these species are sympatric in the southern Mexican highlands. We generated data consisting of three mitochondrial genes and eight nuclear loci for 148 and 68 individuals, respectively. We delimited six lineages in this group, but found strong evidence of mito‐nuclear discordance in sympatric populations of S. spinosus and S. horridus owing to mitochondrial introgression. We used coalescent simulations to differentiate ancestral gene flow from secondary contact, but found mixed support for these two models. Bayesian phylogeography indicated more than 60% range overlap between ancestral S. spinosus and S. horridus populations since the time of their divergence. Isolation–migration analyses, however, revealed near‐zero levels of gene flow between these ancestral populations. Interpreting results from both simulations and empirical data indicate that despite a long history of sympatry among these two species, gene flow in this group has only recently occurred.     

Introgression of mtDNA in Urosaurus lizards: historical and ecological processes

Introgression of mtDNA appears common in animals, but the implications of acquiring a novel mitochondrial genome are not well known. This study investigates mito‐genome introgression between the lizard species Urosaurus graciosus, a thermal specialist, and U. ornatus, a species that occupies a wider range of thermal environments. As ectotherms, their metabolic rate is strongly influenced by the thermal environment; with mitochondria being linked to metabolic rates, overall energy budgets could be impacted by introgression. I use mitochondrial gene trees, inferred from Bayesian analyses of Cyt‐B and ND1 gene sequences, along with morphology and microsatellites from nineteen populations of these two species to address if the direction and location of mito‐nuclear discordance match predictions of introgression resulting from past population expansions. MtDNA is expected to move from resident species into expanding or invading species. Second, does having a heterospecific form of mitochondria impact body size, a trait strongly associated with fitness? Multiple independent introgression events of historic origin were detected. All introgression was unidirectional with U. ornatus‐type mtDNA found in U. graciosus parental type individuals. This result was consistent with population expansions detected in U. graciosus but not U. ornatus. Females with heterospecific mtDNA were significantly smaller than homospecific forms, and heterospecific males had a different relationship of body mass to body length than those with homospecific mtDNA. These changes indicate a potential selective disadvantage for individuals with heterospecific mitochondria and are consistent with the theoretical expectation that deleterious alleles are more likely to persist in expanding populations.     

Distinct sources of gene flow produce contrasting population genetic dynamics at different range boundaries of a Choristoneura budworm

Populations are often exposed to multiple sources of gene flow, but accounts are lacking of the population genetic dynamics that result from these interactions or their effects on local evolution. Using a genomic clines framework applied to 1,195 single nucleotide polymorphisms, we documented genomewide, locus‐specific patterns of introgression between Choristoneura occidentalis biennis spruce budworms and two ecologically divergent relatives, C. o. occidentalis and Choristoneura fumiferana, that it interacts with at alternate boundaries of its range. We observe contrasting hybrid indexes between the two hybrid zones, no overlap in “gene‐flow outliers” (clines showing relatively extreme extents or rates of locus‐specific introgression) and variable linkage disequilibrium among those outliers. At the same time, correlated genomewide rates of introgression between zones suggest the presence of processes common to both boundaries. These findings highlight the contrasting population genetic dynamics that can occur at separate frontiers of a single population, while also suggesting that shared patterns may frequently accompany cases of divergence‐with‐gene‐flow that involve a lineage in common. Our results point to potentially complex evolutionary outcomes for populations experiencing multiple sources of gene flow.     

Evolutionary history of a desert perennial Arnebia szechenyi (Boraginaceae): Intraspecific divergence,regional expansion and asymmetric gene flow

The complex interactions of historical, geological and climatic events on plant evolution have been an important research focus for many years. However, the role of desert formation and expansion in shaping the genetic structures and demographic histories of plants occurring in arid areas has not been well explored. In the present study, we investigated the phylogeography of Arnebia szechenyi, a desert herb showing a near-circular distribution surrounding the Tengger Desert in Northwest China. We measured genetic diversity of populations using three maternally inherited chloroplast DNA (cpDNA) fragments and seven bi-paternally inherited nuclear DNA (nDNA) loci that were sequenced from individuals collected from 16 natural populations across its range and modelled current and historical potential habitats of the species. Our data indicated a considerably high level of genetic variation within A. szechenyi and noteworthy asymmetry in historical migration from the east to the west. Moreover, two nuclear genetic groups of populations were revealed, corresponding to the two geographic regions separated by the Tengger Desert. However, analysis of cpDNA data did not show significant geographic structure. The most plausible explanation for the discrepancy between our findings based on cpDNA and nDNA data is that A. szechenyi populations experienced long periods of geographic isolation followed by range expansion, which would have promoted generalized recombination of the nuclear genome. Our findings further highlight the important role that the Tengger Desert, together with the Helan Mountains, has played in the evolution of desert plants and the preservation of biodiversity in arid Northwest China.     

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

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

Stark sexual display divergence among jumping spider populations in the face of gene flow

Gene flow can inhibit evolutionary divergence by eroding genetic differences between populations. A current aim in speciation research is to identify conditions in which selection overcomes this process. We focused on a state of limited differentiation, asking whether selection enables divergence with gene flow in a set of Habronattus americanus jumping spider populations that exhibit three distinct male sexual display morphs. We found that each population is at high frequency or fixed for a single morph. These strong phenotypic differences contrast with low divergence at 210 AFLP markers, suggesting selection has driven or maintains morph divergence. Coinciding patterns of isolation by distance and ‘isolation by phenotype’ (i.e. increased genetic divergence among phenotypically contrasting populations) across the study area support several alternative demographic hypotheses for display divergence, each of which entails gene flow. Display‐associated structure appears broadly distributed across the genome and the markers producing this pattern do not stand out from background levels of differentiation. Overall, the results suggest selection can promote stark sexual display divergence in the face of gene flow among closely related populations.     

Phylogeography of the sergeants Abudefduf sexfasciatus and A. vaigiensis reveals complex introgression patterns between two widespread and sympatric Indo‐West Pacific reef fishes

On evolutionary timescales, sea level oscillations lead to recurrent spatio‐temporal variation in species distribution and population connectivity. In this situation, applying classical concepts of biogeography is challenging yet necessary to understand the mechanisms underlying biodiversity in highly diverse marine ecosystems such as coral reefs. We aimed at studying the outcomes of such complex biogeographic dynamics on reproductive isolation by sampling populations across a wide spatial range of a species‐rich fish genus: the sergeants (Pomacentridae: Abudefduf). We generated a mutlilocus data set that included ten morpho‐species from 32 Indo‐West Pacific localities. We observed a pattern of mito‐nuclear discordance in two common and widely distributed species: Abudefduf sexfasciatus and Abudefduf vaigiensis. The results showed three regional sublineages (Indian Ocean, Coral Triangle region, western Pacific) in A. sexfasciatus (0.6–1.5% divergence at cytb). The other species, A. vaigiensis, is polyphyletic and consists of three distinct genetic lineages (A, B and C) (9% divergence at cytb) whose geographic ranges overlap. Although A. vaigiensis A and A. sexfasciatus were found to be distinct based on nuclear information, A. vaigiensis A was found to be nested within A. sexfasciatus in the mitochondrial gene tree. A. sexfasciatus from the Coral Triangle region and A. vaigiensis A were not differentiated from each other at the mitochondrial locus. We then used coalescent‐based simulation to characterize a spatially widespread but weak gene flow between the two species. We showed that these fishes are good candidates to investigate the evolutionary complexity of the discrepancies between phenotypic and genetic similarity in closely related species.     

Incongruence between mtDNA and nuclear data in the freshwater mussel genus Cyprogenia (Bivalvia: Unionidae) and its impact on species delineation

Accurately identifying species is a crucial step for developing conservation strategies for freshwater mussels, one of the most imperiled faunas in North America. This study uses genetic data to re‐examine species delineation in the genus Cyprogenia. Historically, Cyprogenia found west of the Mississippi River have been ascribed to Cyprogenia aberti (Conrad 1850 ), and those east of the Mississippi River were classified as Cyprogenia stegaria (Rafinesque 1820). Previous studies using mitochondrial DNA sequences indicated that C. aberti and C. stegaria were not reciprocally monophyletic groups, suggesting the need for systematic revision. We generated a novel dataset consisting of 10 microsatellite loci and combined it with sequence data from the mitochondrial ND1 gene for 223 Cyprogenia specimens. Bayesian analysis of the ND1 nucleotide sequences identified two divergent clades that differ by 15.9%. Members of these two clades occur sympatrically across most sampling locations. In contrast, microsatellite genotypes support recognition of three allopatric clusters defined by major hydrologic basins. The divergent mitochondrial lineages are highly correlated with the color of the conglutinate lures used by mussels to attract and infest host fishes, and tests for selection at the ND1 locus were positive. We infer that the incongruence between mtDNA and microsatellite data in Cyprogenia may be the result of a combination of incomplete lineage sorting and balancing selection on lure color. Our results provide further evidence that mitochondrial markers are not always neutral with respect to selection, and highlight the potential problems of relying on a single‐locus‐marker for delineating species.     

Linked vs unlinked markers: multilocus microsatellite haplotype-sharing as a tool to estimate gene flow and introgression

We have explored the use of multilocus microsatellite haplotypes to study introgression from cultivated (Malus domestica) into wild apple (Malus sylvestris), and to study gene flow among remnant populations of M. sylvestris. A haplotype consisted of alleles at microsatellite loci along one chromosome. As destruction of haplotypes through recombination occurs much faster than loss of alleles due to genetic drift, the lifespan of a multilocus haplotype is much shorter than that of the underlying alleles. When different populations share the same haplotype, this may indicate recent gene flow between populations. Similarly, haplotypes shared between two species would be a strong signal for introgression. As the expected lifespan of a haplotype depends on the strength of the linkage, the length [in centiMorgans (cM)] of the haplotype shared contains information on the number of generations passed. This application of shared haplotypes is distinct from using haplotype-sharing to detect association between markers and a certain trait. We inferred haplotypes for four to eight microsatellite loci on Linkage Group 10 of apple from genotype data using the program phase, and then identified those haplotypes shared between populations and species. Compared with a Bayesian analysis of unlinked microsatellite loci using the program structure, haplotype-sharing detected a partially different set of putative hybrids. Cultivated haplotypes present in M. sylvestris were short (< 1.5 cM), indicating that introgression had taken place many generations ago, except for two Belgian plants that contained a haplotype of 47.1 cM, indicating recent introgression. In the estimation of gene flow, F(ST) based on unlinked loci indicated small (0.032-0.058) but statistically significant differentiation between some populations only. However, various M. sylvestris haplotypes were shared in nearly all pairwise comparisons of populations, and their length indicated recent gene flow. Hence, all Dutch populations should be considered as one conservation unit. The added value of using sharing of multilocus microsatellite haplotypes as a source of population genetic information is discussed.     

Recent divergence with gene flow in Tennessee cave salamanders (Plethodontidae: Gyrinophilus) inferred from gene genealogies

Cave organisms occupy a special place in evolutionary biology because convergent morphologies of many species demonstrate repeatability in evolution even as they obscure phylogenetic relationships. The origin of specialized cave-dwelling species also raises the issue of the relative importance of isolation vs. natural selection in speciation. Two alternative hypotheses describe the origin of subterranean species. The 'climate-relict' model proposes allopatric speciation after populations of cold-adapted species become stranded in caves due to climate change. The 'adaptive-shift' model proposes parapatric speciation driven by divergent selection between subterranean and surface habitats. Our study of the Tennessee cave salamander complex shows that the three nominal forms (Gyrinophilus palleucus palleucus, G. p. necturoides, and G. gulolineatus) arose recently and are genealogically nested within the epigean (surface-dwelling) species, G. porphyriticus. Short branch lengths and discordant gene trees were consistent with a complex history involving gene flow between diverging forms. Results of coalescent-based analysis of the distribution of haplotypes among groups reject the allopatric speciation model and support continuous or recurrent genetic exchange during divergence. These results strongly favour the hypothesis that Tennessee cave salamanders originated from spring salamanders via divergence with gene flow.     

Genome‐wide single‐nucleotide polymorphism data reveal cryptic species within cryptic freshwater snail species—The case of the Ancylus fluviatilis species complex

DNA barcoding utilizes short standardized DNA sequences to identify species and is increasingly used in biodiversity assessments. The technique has unveiled an unforeseeably high number of morphologically cryptic species. However, if speciation has occurred relatively recently and rapidly, the use of single gene markers, and especially the exclusive use of mitochondrial markers, will presumably fail in delimitating species. Therefore, the true number of biological species might be even higher. One mechanism that can result in rapid speciation is hybridization of different species in combination with polyploidization, that is, allopolyploid speciation. In this study, we analyzed the population genetic structure of the polyploid freshwater snail Ancylus fluviatilis, for which allopolyploidization was postulated as a speciation mechanism. DNA barcoding has already revealed four cryptic species within A. fluviatilis (i.e., A. fluviatilis s. str., Ancylus sp. A–C), but early allozyme data even hint at the presence of additional cryptic lineages in Central Europe. We combined COI sequencing with high‐resolution genome‐wide SNP data (ddRAD data) to analyze the genetic structure of A. fluviatilis populations in a Central German low mountain range (Sauerland). The ddRAD data results indicate the presence of three cryptic species within A. fluviatilis s. str. occurring in sympatry and even syntopy, whereas mitochondrial sequence data only support the existence of one species, with shared haplotypes between species. Our study hence points to the limitations of DNA barcoding when dealing with organismal groups where speciation is assumed to have occurred rapidly, for example, through the process of allopolyploidization. We therefore emphasize that single marker DNA barcoding can underestimate the true species diversity and argue in strong favor of using genome‐wide data for species delimitation in such groups.     

Phylogeography and lineage‐specific patterns of genetic diversity and molecular evolution in a group of North American skinks

Geography influences the evolutionary trajectory of species by mediating opportunities for hybridization, gene flow, demographic shifts and adaptation. We sought to understand how geography and introgression can generate species‐specific patterns of genetic diversity by examining phylogeographical relationships in the North American skink species Plestiodon multivirgatus and P. tetragrammus (Squamata: Scincidae). Using a multilocus dataset (three mitochondrial genes, four nuclear genes; a total of 3455 bp) we discovered mito‐nuclear discordance, consistent with mtDNA introgression. We further tested for evidence of species‐wide mtDNA introgression by using comparisons of genetic diversity, selection tests and extended Bayesian skyline analyses. Our findings suggest that P. multivirgatus acquired its mitochondrial genome from P. tetragrammus after their initial divergence. This putative species‐wide mitochondrial capture was further evidenced by statistically indistinguishable substitution rates between mtDNA and nDNA in P. multivirgatus. This rate discrepancy was observed in P. multivirgatus but not P. tetragrammus, which has important implications for studies that combine mtDNA and nDNA sequences when inferring time since divergence between taxa. Our findings suggest that by facilitating opportunities for interspecific introgression, geography can alter the course of molecular evolution between recently diverged lineages.     

Adaptation with gene flow across the landscape in a dune sunflower

Isolation by adaptation increases divergence at neutral loci when natural selection against immigrants reduces the rate of gene flow between different habitats. This can occur early in the process of adaptive divergence and is a key feature of ecological speciation. Despite the ability of isolation by distance (IBD) and other forms of landscape resistance to produce similar patterns of neutral divergence within species, few studies have used landscape genetics to control for these other forces. We have studied the divergence of Helianthus petiolaris ecotypes living in active sand dunes and adjacent non-dune habitat, using landscape genetics approaches, such as circuit theory and multiple regression of distance matrices, in addition to coalescent modelling. Divergence between habitats was significant, but not strong, and was shaped by IBD. We expected that increased resistance owing to patchy and unfavourable habitat in the dunes would contribute to divergence. Instead, we found that landscape resistance models with lower resistance in the dunes performed well as predictors of genetic distances among subpopulations. Nevertheless, habitat class remained a strong predictor of genetic distance when controlling for isolation by resistance and IBD. We also measured environmental variables at each site and confirmed that specific variables, especially soil nitrogen and vegetation cover, explained a greater proportion of variance in genetic distance than did landscape or the habitat classification alone. Asymmetry in effective population sizes and numbers of migrants per generation was detected using coalescent modelling with Bayesian inference, which is consistent with incipient ecological speciation being driven by the dune habitat.     

Isolation with asymmetric gene flow during the nonsynchronous divergence of dry forest birds

Dry forest bird communities in South America are often fragmented by intervening mountains and rainforests, generating high local endemism. The historical assembly of dry forest communities likely results from dynamic processes linked to numerous population histories among codistributed species. Nevertheless, species may diversify in the same way through time if landscape and environmental features, or species ecologies, similarly structure populations. Here we tested whether six co‐distributed taxon pairs that occur in the dry forests of the Tumbes and Marañón Valley of northwestern South America show concordant patterns and modes of diversification. We employed a genome reduction technique, double‐digest restriction site‐associated DNA sequencing, and obtained 4407–7186 genomewide SNPs. We estimated demographic history in each taxon pair and inferred that all pairs had the same best‐fit demographic model: isolation with asymmetric gene flow from the Tumbes into the Marañón Valley, suggesting a common diversification mode. Overall, we also observed congruence in effective population size (N_e) patterns where ancestral N_e were 2.9–11.0× larger than present‐day Marañón Valley populations and 0.3–2.0× larger than Tumbesian populations. Present‐day Marañón Valley N_e was smaller than Tumbes. In contrast, we found simultaneous population isolation due to a single event to be unlikely as taxon pairs diverged over an extended period of time (0.1–2.9 Ma) with multiple nonoverlapping divergence periods. Our results show that even when populations of codistributed species asynchronously diverge, the mode of their differentiation can remain conserved over millions of years. Divergence by allopatric isolation due to barrier formation does not explain the mode of differentiation between these two bird assemblages; rather, migration of individuals occurred before and after geographic isolation.     

Frequent mitochondrial gene introgression among high elevation Tibetan megophryid frogs revealed by conflicting gene genealogies

Historical mitochondrial introgression causes differences between a species' mitochondrial gene genealogy and its nuclear gene genealogy, making tree-based species delineation ambiguous. Using sequence data from one mitochondrial gene (cytochrome b ) and three nuclear genes (introns), we examined the evolutionary history of four high elevation Tibetan megophryid frog species, Scutiger boulengeri , Scutiger glandulatus , Scutiger mammatus and Scutiger tuberculatus . The three nuclear genes shared a similar history but the mitochondrial gene tree suggested a drastically different evolutionary scenario. The conflicts between them were explained by multiple episodes of mitochondrial introgression events via historical interspecific hybridization. 'Foreign' mitochondrial genomes might have been fixed in populations and extended through a large portion of the species' distribution. Some hybridization events were probably as old as 10 Myr, while others were recent. An F₁ hybrid was also identified. Historical hybridization events among the four species appeared to be persistent and were not restricted to the period of Pleistocene glaciation, as in several other well-studied cases. Furthermore, hybridization involved several species and occurred in multiple directions, and there was no indication of one mitochondrial genome being superior to others. In addition, incomplete lineage sorting resulting from budding speciation may have also explained some discrepancies between the mitochondrial DNA and nuclear gene trees. Combining all evidences, the former ' Scutiger mammatus ' appeared to be two species, including a new species. With the availability of a wide range of highly variable nuclear gene markers, we recommend using a combination of mitochondrial gene and multiple nuclear genes to reveal a complete species history.     

Divergence with gene flow in the rock-dwelling cichlids of Lake Malawi

Within the past two million years, more than 450 species of haplochromine cichlids have diverged from a single common ancestor in Lake Malawi. Several factors have been implicated in the diversification of this monophyletic clade, including changes in lake level and low levels of gene flow across limited geographic scales. The objectives of this study were to determine the effect of recent lake-level fluctuations on patterns of allelic diversity in the genus Metriaclima, to describe the patterns of population structure within this genus, and to identify barriers to migration. This was accomplished through an analysis of allele frequencies at four microsatellite loci. Twelve populations spanning four species within Metriaclima were surveyed. The effect of lake-level fluctuations can be seen in the reduced genetic diversity of the most recently colonized sites; however, genetic diversity is not depressed at the species level. Low levels of population structure exist among populations, yet some gene flow persists across long stretches of inhospitable habitat. No general barrier to migration was identified. The results of this study are interpreted with respect to several speciation models. Divergence via population bottlenecks is unlikely due to the large allelic diversity observed within each species. Genetic drift and microallopatric divergence are also rejected because some gene flow does occur between adjacent populations. However, the reduced levels of gene flow between populations does suggest that minor changes in the selective environment could cause the divergence of populations.     

Assessment of gene flow across a hybrid zone in red-tailed chipmunks (Tamias ruficaudus)

The role of hybridization in animal speciation is controversial and recent research has challenged the long-standing criterion of complete reproductive isolation to attain species status. The speciation-with-gene-flow model posits that the genome is semi-permeable and hybridization may be a phase in the process of divergence. Here, we apply these concepts to a previously identified zone of mtDNA introgression between the two strongly morphologically differentiated subspecies of red-tailed chipmunk ( Tamias ruficaudus ) in the US Inland Northwest. Using multilocus genotype data from the southern, older contact zone, we demonstrate that neutral gene flow is unusually low between the subspecies across the Lochsa River. This is geographically congruent with the discontinuity in bacular morphology, indicating that the cline of mitochondrial DNA (mtDNA) haplotypes is displaced. Furthermore, we elucidate the evolutionary forces responsible by testing hypotheses of lineage sorting and hybridization. We determined that introgressive hybridization is the cause of mtDNA/morphology incongruence because there are non-zero levels of migration and gene flow. Although our estimate of the age of the hybrid zone has wide credibility intervals, the hybridization events occurred in the Late Pleistocene and the divergence occurred in the Middle Pleistocene. Finally, we assessed substructure within and adjacent to the hybrid zone and found that the hybrid zone constitutes a set of populations that are genetically differentiated from parental sets of populations; therefore, hybridization in this system is not likely an evolutionary sink, but has generated novel combinations of genotypes.