Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Speciation occurs when populations diverge and become reproductively isolated from each other. Natural selection is commonly accepted to play a large role in this process, and it has been widely assumed that reproductive isolation often results as a by‐product of divergence driven by adaptation in allopatry. When such populations come into secondary contact, reinforcement can act to strengthen reproductive isolation, but the frequency and importance of this process are still unknown. Here, we explored genomic signatures of selection in allopatry and sympatry for loci associated with reproductive isolation using a natural primate hybrid zone. By analysing reduced‐representation sequencing data, we quantified admixture and population structure across a howler monkey hybrid zone and examined the relationship between locus‐specific differentiation and introgression. We detected extensive admixture that was mostly limited to the narrow contact zone. Loci with reduced introgression into the heterospecific genomic background (the pattern expected for loci associated with reproductive isolation due to selection against hybrids) were significantly more differentiated between allopatric parental populations than loci with neutral and increased introgression, supporting the hypothesis that reproductive isolation is a by‐product of divergence in allopatry. Further, loci with reduced introgression showed greater differentiation in sympatry than in allopatry, suggesting a role for reinforcement. Thus, our results reflect multiple forms of selection that have shaped reproductive isolation in this system. We conclude that reproductive isolation may have initially been driven by divergence in allopatry, but later reinforced by divergent selection in sympatry.     

Strong postmating reproductive isolation in Mimulus section Eunanus

Postmating reproductive isolation can help maintain species boundaries when premating barriers to reproduction are incomplete. The strength and identity of postmating reproductive barriers are highly variable among diverging species, leading to questions about their genetic basis and evolutionary drivers. These questions have been tackled in model systems but are less often addressed with broader phylogenetic resolution. In this study we analyse patterns of genetic divergence alongside direct measures of postmating reproductive barriers in an overlooked group of sympatric species within the model monkeyflower genus, Mimulus. Within this Mimulus brevipes species group, we find substantial divergence among species, including a cryptic genetic lineage. However, rampant gene discordance and ancient signals of introgression suggest a complex history of divergence. In addition, we find multiple strong postmating barriers, including postmating prezygotic isolation, hybrid seed inviability and hybrid male sterility. M. brevipes and M. fremontii have substantial but incomplete postmating isolation. For all other tested species pairs, we find essentially complete postmating isolation. Hybrid seed inviability appears linked to differences in seed size, providing a window into possible developmental mechanisms underlying this reproductive barrier. While geographic proximity and incomplete mating isolation may have allowed gene flow within this group in the distant past, strong postmating reproductive barriers today have likely played a key role in preventing ongoing introgression. By producing foundational information about reproductive isolation and genomic divergence in this understudied group, we add new diversity and phylogenetic resolution to our understanding of the mechanisms of plant speciation.     

Gene flow-dependent genomic divergence between Anopheles gambiae M and S forms

Anopheles gambiae sensu stricto exists as two often-sympatric races termed the M and S molecular forms, characterized by fixed differences at an X-linked marker. Extreme divergence between M and S forms at pericentromeric "genomic islands" suggested that selection on variants therein could be driving interform divergence in the presence of ongoing gene flow, but recent work has detected much more widespread genomic differentiation. Whether such genomic islands are important in reproductive isolation or represent ancestral differentiation preserved by low recombination is currently unclear. A critical test of these competing hypotheses could be provided by comparing genomic divergence when rates of recent introgression vary. We genotyped 871 single nucleotide polymorphisms (SNPs) in A. gambiae sensu stricto from locations of M and S sympatry and allopatry, encompassing the full range of observed hybridization rates (0-25%). M and S forms were readily partitioned based on genomewide SNP variation in spite of evidence for ongoing introgression that qualitatively reflects hybridization rates. Yet both the level and the heterogeneity of genomic divergence varied markedly in line with levels of introgression. A few genomic regions of differentiation between M and S were common to each sampling location, the most pronounced being two centromere-proximal speciation islands identified previously but with at least one additional region outside of areas expected to exhibit reduced recombination. Our results demonstrate that extreme divergence at genomic islands does not simply represent segregating ancestral polymorphism in regions of low recombination and can be resilient to substantial gene flow. This highlights the potential for islands comprising a relatively small fraction of the genome to play an important role in early-stage speciation when reproductive isolation is limited.     

A new multiplex SNP genotyping assay for detecting hybridization and introgression between the M and S molecular forms of Anopheles gambiae

The M and S forms of Anopheles gambiae have been the subject of intense study, but are morphologically indistinguishable and can only be identified using molecular techniques. PCR‐based assays to distinguish the two forms have been designed and applied widely. However, the application of these assays towards identifying hybrids between the two forms, and backcrossed hybrids in particular, has been problematic as the currently available diagnostic assays are based on single locus and/or are located within a multicopy gene. Here, we present an alternative genotyping method for detecting hybridization and introgression between M and S molecular forms based on a multilocus panel of single‐nucleotide polymorphisms (SNPs) fixed between the M and S forms. The panel of SNPs employed is located in so‐called islands of divergence leading us to describe this method as the ‘Divergence Island SNP’ (DIS) assay. We show this multilocus SNP genotyping approach can robustly and accurately detect F1 hybrids as well as backcrossed individuals.     

No evidence for biased co-transmission of speciation islands in Anopheles gambiae

Genome-scale scans have revealed highly heterogeneous levels of divergence between closely related taxa in many systems. Generally, a small number of regions show high differentiation, with the rest of the genome showing no or only low levels of divergence. These patterns have been interpreted as evidence for ongoing speciation-with-gene-flow, with introgression homogenizing the whole genome except loci involved in reproductive isolation. However, as the number of selected loci increases, the probability of introgression at unselected loci decreases unless there is a transmission ratio distortion causing an over-representation of specific combinations of alleles. Here we examine the transmission of three 'speciation islands' that contain fixed differences between the M and S forms of the mosquito, Anopheles gambiae. We made reciprocal crosses between M and S parents and genotyped over 2000 F(2) individuals, developing a hierarchical likelihood model to identify specific genotypes that are under- or over-represented among the recombinant offspring. Though our overall results did not match the expected number of F(2) genotypes, we found no biased co-transmission among M or S alleles in the three islands. Our likelihood model did identify transmission ratio distortion at two of the three islands, but this distortion was small (approx. 3%) and in opposite directions for the two islands. We discuss how our results impinge on hypotheses of current gene flow between M and S and ongoing speciation-with-gene-flow in this system.     

Sequence diversity, reproductive isolation and species concepts in Saccharomyces

Using the biological species definition, yeasts of the genus Saccharomyces sensu stricto comprise six species and one natural hybrid. Previous work has shown that reproductive isolation between the species is due primarily to sequence divergence acted upon by the mismatch repair system and not due to major gene differences or chromosomal rearrangements. Sequence divergence through mismatch repair has also been shown to cause partial reproductive isolation among populations within a species. We have surveyed sequence variation in populations of Saccharomyces sensu stricto yeasts and measured meiotic sterility in hybrids. This allows us to determine the divergence necessary to produce the reproductive isolation seen among species. Rather than a sharp transition from fertility to sterility, which may have been expected, we find a smooth monotonic relationship between diversity and reproductive isolation, even as far as the well-accepted designations of S. paradoxus and S. cerevisiae as distinct species. Furthermore, we show that one species of Saccharomyces--S. cariocanus--differs from a population of S. paradoxus by four translocations, but not by sequence. There is molecular evidence of recent introgression from S. cerevisiae into the European population of S. paradoxus, supporting the idea that in nature the boundary between these species is fuzzy.     

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.     

Differential patterns of introgression across the X chromosome in a hybrid zone between two species of house mice

A complete understanding of the speciation process requires the identification of genomic regions and genes that confer reproductive barriers between species. Empirical and theoretical research has revealed two important patterns in the evolution of reproductive isolation in animals: isolation typically arises as a result of disrupted epistatic interactions between multiple loci and these disruptions map disproportionately to the X chromosome. These patterns suggest that a targeted examination of natural gene flow between closely related species at X-linked markers with known positions would provide insight into the genetic basis of speciation. We take advantage of the existence of genomic data and a well-documented European zone of hybridization between two species of house mice, Mus domesticus and M. musculus, to conduct such a survey. We evaluate patterns of introgression across the hybrid zone for 13 diagnostic X-linked loci with known chromosomal positions using a maximum likelihood model. Interlocus comparisons clearly identify one locus with reduced introgression across the center of the hybrid zone, pinpointing a candidate region for reproductive isolation. Results also reveal one locus with high frequencies of M. domesticus alleles in populations on the M. musculus side of the zone, suggesting the possibility that positive selection may act to drive the spread of alleles from one species on to the genomic background of the other species. Finally, cline width and cline center are strongly positively correlated across the X chromosome, indicating that gene flow of the X chromosome may be asymmetrical. This study highlights the utility of natural populations of hybrids for mapping speciation genes and suggests that the middle of the X chromosome may be important for reproductive isolation between species of house mice.     

Using differential introgression in hybrid zones to identify genomic regions involved in speciation

Hybrids between species provide information about the evolutionary processes involved in divergence. In addition to creating hybrids in the laboratory, biologists can take advantage of natural hybrid zones to understand the factors that shape gene flow between divergent lineages. In the early stages of speciation, most regions of the genome continue to flow freely between populations. Alternatively, the subset of the genome that confers reproductive barriers between nascent species is expected to reject introgression. Now enabled by advances in genomics, this perspective is motivating detailed comparisons of gene flow across genomic regions in hybrid zones. Here, I review methods for measuring and interpreting introgression at multiple loci in hybrid zones, focusing on the problem of identifying loci that contribute to reproductive isolation. Emerging patterns from multi-locus studies of hybrid zones are highlighted, including remarkable variance in introgression across the genome. Although existing methods have been useful, there is scope for development of new analytical approaches that better connect differential patterns of gene flow in hybrid zones with current knowledge of speciation mechanisms. I outline future prospects for differential introgression studies on a genomic scale.     

Quantified reproductive isolation in Heliconius butterflies: Implications for introgression and hybrid speciation

Heliconius butterflies have become a model for the study of speciation with gene flow. For adaptive introgression to take place, there must be incomplete barriers to gene exchange that allow interspecific hybridization and multiple generations of backcrossing. The recent publication of estimates of individual components of reproductive isolation between several species of butterflies in the Heliconius melpomene–H. cydno clade allowed us to calculate total reproductive isolation estimates for these species. According to these estimates, the butterflies are not as promiscuous as has been implied. Differences between species are maintained by intrinsic mechanisms, while reproductive isolation of geographical races within species is mainly due to allopatry. We discuss the implications of this strong isolation for basic aspects of the hybrid speciation with introgression hypothesis.     

SIGNIFICANT ROLE FOR HISTORICAL EFFECTS IN THE EVOLUTION OF REPRODUCTIVE ISOLATION: EVIDENCE FROM PATTERNS OF INTROGRESSION BETWEEN THE CYPRINID FISHES,LUXILUS CORNUTUS AND LUXILUS CHRYSOCEPHALUS

Samples of Luxilus cornutus, Luxilus chrysocephalus, and their hybrids were collected along hypothesized routes of dispersal from Pleistocene refugia to examine the significance of geographic variation in patterns of introgression between these species. Patterns of allozyme and mitochondrial DNA (mtDNA) variation were generally consistent with those from previous studies. Tests of Hardy-Weinberg equilibrium revealed significant deficiencies of heterozygotes in all samples, indicating some form of reproductive isolation. Mitochondrial DNAs of each species were not equally represented in F₁ hybrids; however, this bias was eliminated when the two largest samples were excluded from the analysis. Backcross hybrids exhibited biased mtDNA introgression, as samples from Lake Erie (eastern) and Lake Michigan (western) drainages showed significant excesses of mtDNAs from L. chrysocephalus and L. cornutus, respectively, relative to frequencies of diagnostic allozyme markers. The extent and direction of allozyme and mtDNA introgression was quantified by calculating isolation index values from morphologically “pure” individuals of each species from each locality. Analysis of variance of these measures identified limited introgression of allozyme variants with no geographic pattern, but significant differences in direction of mtDNA introgression between drainages (i.e., postglacial dispersal route). Association between patterns of mtDNA introgression and dispersal route across the latitudinal width of the contact zone is best explained by genetic divergence during past isolation of ancestral populations from these drainages. These results identify a significant role for historical effects in the evolution of reproductive isolation and the process of speciation.     

Ecological zones rather than molecular forms predict genetic differentiation in the malaria vector Anopheles gambiae s.s. in Ghana

The malaria mosquito Anopheles gambiae s.s. is rapidly becoming a model for studies on the evolution of reproductive isolation. Debate has centered on the taxonomic status of two forms (denoted M and S) within the nominal taxon identified by point mutations in the X-linked rDNA region. Evidence is accumulating that there are significant barriers to gene flow between these forms, but that the barriers are not complete throughout the entire range of their distribution. We sampled populations from across Ghana and southern Burkina Faso, West Africa, from areas where the molecular forms occurred in both sympatry and allopatry. Neither Bayesian clustering methods nor F(ST)-based analysis of microsatellite data found differentiation between the M and S molecular forms, but revealed strong differentiation among different ecological zones, irrespective of M/S status and with no detectable effect of geographical distance. Although no M/S hybrids were found in the samples, admixture analysis detected evidence of contemporary interform gene flow, arguably most pronounced in southern Ghana where forms occur sympatrically. Thus, in the sampled area of West Africa, lack of differentiation between M and S forms likely reflects substantial introgression, and ecological barriers appear to be of greater importance in restricting gene flow.     

GENETIC LEAKAGE AFTER ADAPTIVE AND NONADAPTIVE DIVERGENCE IN THE ENSATINA ESCHSCHOLTZII RING SPECIES

The salamander Ensatina eschscholtzii is an example of a ring species in which extant intermediate stages of terminal forms have a nearly continuous range, offering replicated interactions at several stages of divergence. We employ a greatly expanded allozyme database and individual-based analyses to separate the effects of divergence time and gene flow to evaluate how gradual divergence of populations around the ring contributes to the development of reproductive isolation. Despite the high degree of genetic ( D ≤ 0.39) and ecomorphological divergence observed in secondary contacts around the ring, reproductive isolation or rare hybridization is observed only at the terminus of the ring. Instead, in the secondary contacts sampled around the ring, hybrids are common and reproductively successful, enabling genetic leakage between parental genomes and the potential for genetic merger. Nevertheless, genetic admixture is geographically broad (<100 km) only in contacts between ecomorphologically similar populations (within subspecies). When divergence is accompanied by alternative patterns of adaptive divergence (between subspecies), zones of intergradation are narrower and affect populations only locally (>8 km). Diversification and consequent genetic interactions in Ensatina reveal a continuum between populations, ecological races, and species, where polytypic traits and high genetic differentiation are maintained without reproductive isolation.     

Maintenance of species boundaries in a Neotropical radiation of Begonia

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by‐product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later‐generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Ecological divergence plays an important role in strong but complex reproductive isolation in campions (Silene)*

New species arise through the evolution of reproductive barriers between formerly interbreeding lineages. Yet, comprehensive assessments of potential reproductive barriers, which are needed to make inferences on processes driving speciation, are only available for a limited number of systems. In this study, we estimated individual and cumulative strengths of seven prezygotic and six postzygotic reproductive barriers between the recently diverged taxa Silene dioica (L.) Clairv. and S. latifolia Poiret using both published and new data. A combination of multiple partial reproductive barriers resulted in near‐complete reproductive isolation between S. dioica and S. latifolia, consistent with earlier estimates of gene flow between the taxa. Extrinsic barriers associated with adaptive ecological divergence were most important, while intrinsic postzygotic barriers had moderate individual strength but contributed only little to total reproductive isolation. These findings are in line with ecological divergence as driver of speciation. We further found extensive variation in extrinsic reproductive isolation, ranging from sites with very strong selection against migrants and hybrids to intermediate sites where substantial hybridization is possible. This situation may allow for, or even promote, heterogeneous genetic divergence.     

Multiple gene genealogies reveal asymmetrical hybridization and introgression among strongylocentrotid sea urchins

The evolution of incompatibilities between eggs and sperm is thought to play important roles in establishing and maintaining reproductive isolation among species of broadcast-spawning marine invertebrates. However, the effectiveness of gametic isolation in initiating the speciation process and/or in limiting the introgression of genes among species at later stages of divergence remains largely unknown. In the present study, we collected DNA sequence data from five loci in four species of Strongylocentrotus sea urchins ( S. droebachiensis , S. pallidus , S. purpuratus , and S. franciscanus ) to test whether the susceptibility of S. droebachiensis eggs to fertilization by heterospecific sperm results in gene flow between species. Despite the potential for introgression, a small but statistically significant signal of introgression was observed only between the youngest pair of sister taxa ( S. pallidus and S. droebachiensis ) that was strongly asymmetrical (from the former into the latter). No significant gene flow was observed for either S. purpuratus or S. franciscanus despite the ability of their sperm to readily fertilize the eggs of S. droebachiensis . Our results demonstrate that asymmetrical gamete compatibilities in strongylocentrotids can give rise to asymmetrical patterns of introgression but suggest that gamete traits alone cannot be responsible for maintaining species integrities. The genetic boundaries between strongylocentrotid urchin species in the northeast Pacific appear to be related to postzygotic isolating mechanisms that scale with divergence times and not intrinsic gametic incompatibilities per se .     

The genetic mosaic suggests a new role for hitchhiking in ecological speciation

Early in ecological speciation, the genomically localized effects of divergent selection cause heterogeneity among loci in divergence between incipient species. We call this pattern of genomic variability in divergence the 'genetic mosaic of speciation'. Previous studies have used F(ST) outliers as a way to identify divergently selected genomic regions, but the nature of the relationship between outlier loci and quantitative trait loci (QTL) involved in reproductive isolation has not yet been quantified. Here, we show that F(ST) outliers between a pair of incipient species are significantly clustered around QTL for traits that cause ecologically based reproductive isolation. Around these key QTL, extensive 'divergence hitchhiking' occurs because reduced inter-race mating and negative selection decrease the opportunity for recombination between chromosomes bearing different locally adapted QTL alleles. Divergence hitchhiking is likely to greatly increase the opportunity for speciation in populations that are sympatric, regardless of whether initial divergence was sympatric or allopatric. Early in ecological speciation, analyses of population structure, gene flow or phylogeography based on different random or arbitrarily chosen neutral markers should be expected to conflict--only markers in divergently selected genomic regions will reveal the evolutionary history of adaptive divergence and ecologically based reproductive isolation. Species retain mosaic genomes for a very long time, and gene exchange in hybrid zones can vary dramatically among loci. However, in hybridizing species, the genomic regions that affect ecologically based reproductive isolation are difficult to distinguish from regions that have diverged for other reasons.     

A biogeographic genetic approach for testing the role of reinforcement: the case of Drosophila pseudoobscura and D. persimilis

Abstract.— The role of reinforcement in speciation can be explained by two distinct models. In model I, two diverged populations hybridize and produce fertile hybrids that successfully backcross (hybridization with gene flow). In model II, two populations hybridize but succeeding backcrosses are unproductive (hybridization without gene flow). Using Drosophila persimilis and D. pseudoobscura , we have tested model I by comparing the extent of heterospecific introgression in sympatric versus allopatric populations. We show that certain expectations of this particular model of reinforcement, which is based on hybridization and gene flow between divergent populations after secondary contact, are not realized in these two species. The evidence consists of the similarity of genetic distances as well as proportions of unique/rare alleles between sympatric and allopatric heterospecific populations and a negative correlation between genetic distance and geographical distance between heterospecific populations, which suggests ecological differentiation. This approach in quantifying differential gene flow has important consequences to studies that compare sympatric and allopatric isolation using genetic distance. Following model I, one would expect a pattern of higher prezygotic isolation in sympatric species compared to allopatric species of the same genetic distance simply as a result of an underestimation of genetic distance due to introgression between sympatric populations. We suggest more parsimonious explanations such as reinforcement without genetic exchange (model II) and ecological differentiation, which require high levels of preexisting reproductive isolation between populations.     

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y‐chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least‐divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi‐absence of Y‐chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non‐neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.