Steep clines within a highly permeable genome across a hybrid zone between two subspecies of the European rabbit

Maintenance of genetic distinction in the face of gene flow is an important aspect of the speciation process. Here, we provide a detailed spatial and genetic characterization of a hybrid zone between two subspecies of the European rabbit. We examined patterns of allele frequency change for 22 markers located on the autosomes, X‐chromosome, Y‐chromosome and mtDNA in 1078 individuals sampled across the hybrid zone. While some loci revealed extremely wide clines (w ≥ 300 km) relative to an estimated dispersal of 1.95–4.22 km/generation, others showed abrupt transitions (w ≈ 10 km), indicating localized genomic regions of strong selection against introgression. The subset of loci showing steep clines had largely coincident centers and stepped changes in allele frequency that did not co‐localize with any physical barrier or ecotone, suggesting that the rabbit hybrid zone is a tension zone. The steepest clines were for X‐ and Y‐chromosome markers. Our results are consistent with previous inference based on DNA sequence variation of individuals sampled in allopatry in suggesting that a large proportion of each genome has escaped the overall barrier to gene flow in the middle of the hybrid zone. These results imply an old history of hybridization and high effective gene flow and anticipate that isolation factors should often localize to small genomic regions.     

Y not introgress? Insights into the genetics of speciation in European rabbits

Hybridization between genetically divergent populations is common in nature. By comparing the rate of gene flow throughout the genome, loci that impede genetic mixing, and therefore underlie reproductive isolation, can be identified, providing key insights into the process of speciation. In a previous issue of Molecular Ecology , Geraldes et al . (2008) report the geographical distribution of Y chromosome lineages in wild populations of rabbits from across the Iberian Peninsula and southern France. Y chromosomes showed a high level of differentiation between rabbit subspecies on either side of a hybrid zone, despite clear evidence for gene flow at other loci. This provocative result suggests a role for the Y chromosome in reproductive isolation, and adds to a growing list of nascent species with mosaic patterns of differentiation across their genomes.     

Contrasting patterns of introgression at X-linked loci across the hybrid zone between subspecies of the European rabbit (Oryctolagus cuniculus)

Hybrid zones provide an excellent opportunity for studying the consequences of genetic changes between closely related taxa. Here we investigate patterns of genetic variability and gene flow at four X-linked loci within and between the two subspecies of European rabbit (Oryctolagus cuniculus cuniculus and O. c. algirus). Two of these genes are located near the centromere and two are located near the telomeres. We observed a deep split in the genealogy of each gene with the root located along the deepest branch in each case, consistent with the evolution of these subspecies in allopatry. The two centromeric loci showed low levels of variability, high levels of linkage disequilibrium, and little introgression between subspecies. In contrast, the two telomeric loci showed high levels of variability, low levels of linkage disequilibrium, and considerable introgression between subspecies. These data are consistent with suppression of recombination near the centromere of the rabbit X chromosome. These observations support a view of speciation where genomic incompatibilities at different loci in the genome create localized differences in levels of gene flow between nascent species.     

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.     

Discordant Phylogeographic Patterns between the Y Chromosome and Mitochondrial DNA in the House Mouse: Selection on the Y Chromosome?

We have compared patterns of geographic variation and molecular divergence of mitochondrial DNA (mtDNA) and Y chromosome over the range of the different subspecies of Mus musculus. MtDNA was typed for 305 nucleotides in the control region, the Y chromosome for 834 base pairs (bp) in Zfy introns and 242 bp in Sry, a Zfy2 18-bp deletion, and two microsatellites. Apparent discrepancies exist between the distributions of the lineages of mtDNA and of the two major Y-chromosome lineages thus defined: some subspecies share the same mtDNA lineage but have different Y-chromosome lineages or vice versa. One microsatellite reveals a geographically clustered variation inside the distribution of each Y-chromosome lineage, showing that new Y-chromosome variants can rapidly spread locally. The two major Y-chromosome lineages have a divergence time only about one fourth of that between mtDNA lineages. Although this recent coalescence of the Y chromosomes between subspecies could partly be due to a lower ancestral polymorphism of the Y chromosome, it suggests that secondary introgression after the radiation of the subspecies might have occurred. There is evidence that the differentiation of the Y-chromosome lineages contributes to partial reproductive isolation between subspecies, and patterns of molecular evolution suggest that selection has played a role in the rapid spread across subspecies.     

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.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

Counterselection on sex chromosomes in the Mus musculus European hybrid zone

The extent to which alleles can disperse across a hybrid zone depends on the selection they are subjected to in the hybrid genetic background or, for those that are selectively neutral, on their ability to escape from the unfavourable environment by recombination. Three markers spanning a 45 cM segment in the center of the X chromosome were used to investigate the degree to which selection against X chromosome linked genes helps to maintain the barrier to gene flow in the hybrid zone between Mus musculus domesticus and M. m. musculus in Denmark. The introgression of all the sex chromosome specific markers was more limited than that of the autosomal enzymes (Idh1, Amy, Gpd1, Pgm1, Es1, Es2, Mpi, Np1, Es10, Sod1) and the mitochondrial DNA. The cline for DXPas2, which is in the center of the X chromosome, is extremely steep and shows that certain genes located in this region are strongly selected against in the hybrid background. The clines of the other two X-linked markers, Hprt and DXPas1, and of the Y chromosome are not as abrupt and all three have similar asymmetric introgression patterns. Although the musculus variants appear to behave in much the same way as those of the autosomal genes, the domesticus variants do not introgress. The results show that X-linked and to a lesser extent Y-linked genes are more strongly selected against in the hybrid genome than the mitochondrial genome or the different autosomal loci. This suggests that co-adapted gene systems involving the sex chromosomes may play an important role in the hybrid breakdown between the two subspecies.     

SPECIATION IN THE EUROPEAN RABBIT (ORYCTOLAGUS CUNICULUS): ISLANDS OF DIFFERENTIATION ON THE X CHROMOSOME AND AUTOSOMES

Studies of gene flow between closely related taxa can provide insight into the genetic basis of speciation. To evaluate the importance of the X chromosome in reproductive isolation between subspecies of the European rabbit and to study the genomic scale over which islands of differentiation extend, we resequenced a total of 34 loci distributed along the X chromosome and chromosome 14. Previous studies based on few markers suggested that loci in centromeric regions were highly differentiated between rabbit subspecies, whereas loci in telomeric regions were less differentiated. Here, we confirmed this finding but also discovered remarkable variation in levels of differentiation among loci, with F_ST values from nearly 0 to 1. Analyses using isolation‐with‐migration models suggest that this range appears to be largely explained by differential levels of gene flow among loci. The X chromosome was significantly more differentiated than the autosomes. On chromosome 14, differentiation decayed very rapidly at increasing distances from the centromere, but on the X chromosome distinct islands of differentiation encompassing several megabases were observed both at the centromeric region and along the chromosome arms. These findings support the idea that the X chromosome plays an important role in reproductive isolation between rabbit subspecies. These results also demonstrate the mosaic nature of the genome at species boundaries.     

The uncharacterized gene 1700093K21Rik and flanking regions are correlated with reproductive isolation in the house mouse,Mus musculus

Reproductive barriers exist between the house mouse subspecies, Mus musculus musculus and M. m. domesticus, members of the Mus musculus species complex, primarily as a result of hybrid male infertility, and a hybrid zone exists where their ranges intersect in Europe. Using single nucleotide polymorphisms (SNPs) diagnostic for the two taxa, the extent of introgression across the genome was previously compared in these hybrid populations. Sixty-nine of 1316 autosomal SNPs exhibited reduced introgression in two hybrid zone transects suggesting maladaptive interactions among certain loci. One of these markers is within a region on chromosome 11 that, in other studies, has been associated with hybrid male sterility of these subspecies. We assessed sequence variation in a 20 Mb region on chromosome 11 flanking this marker, and observed its inclusion within a roughly 150 kb stretch of DNA showing elevated sequence differentiation between the two subspecies. Four genes are associated with this genomic subregion, with two entirely encompassed. One of the two genes, the uncharacterized 1700093K21Rik gene, displays distinguishing features consistent with a potential role in reproductive isolation between these subspecies. Along with its expression specifically within spermatogenic cells, we present various sequence analyses that demonstrate a high rate of molecular evolution of this gene, as well as identify a subspecies amino acid variant resulting in a structural difference. Taken together, the data suggest a role for this gene in reproductive isolation.     

Comparison of Y chromosome and mtDNA phylogenies leads to unique inferences of macaque evolutionary history

We report here the results of one of the first analyses to use male-specific nuclear markers in elucidating primate phylogenetic relationships at the intrageneric level. Two closely linked Y chromosome markers, TSPY and SRY, were sequenced for a total of 3100 bases. Forty-four macaques, representing 18 of the 19 recognized species, were sequenced for the full 3.1 kb, as was 1 individual from each of the following outgroup genera: Papio, Theropithecus, Mandrillus, Allenopithecus,Cercopithecus, Trachypithecus, Presbytis, and Homo. In contrast to recent mtDNA phylogenies, Y chromosome loci support four monophyletic species groups, including a sinica group containing M. arctoides-a classification largely congruent with those of Fooden and Delson. Comparison of mtDNA and Y chromosome phylogenies highlight (1) a potential hybrid origin of Macaca arctoides from M. fascicularis and proto-M. assamensis/thibetana and (2) cases of mitochondrial paraphyly in macaque species whose Y chromosome lineages are monophyletic-a probable evolutionary consequence of philopatric females vs dispersing males. These results raise the question of whether a phylogenetic tree should be a topology of species origins or a depiction of more current species relationships, including subsequent episodes of introgression.     

Genetic analysis of autosomal and X-linked markers across a mouse hybrid zone

In this paper, we present results of the first comprehensive study of the introgression of both autosomal and sex-chromosome markers across the central European portion of the hybrid zone between two house mouse subspecies, Mus musculus musculus and M. m. domesticus. More than 1800 individuals sampled from 105 sites were analyzed with a set of allozyme loci (hopefully representing neutral or nearly neutral markers) and X-linked loci (which are assumed to be under selection). The zone center is best modeled as a single straight line independent of fine-scale local geographic or climatic conditions, being maintained by a balance between dispersal and selection against hybrids. The width (w) of the multilocus autosomal cline was estimated as 9.6 km whereas the estimate for the compound X-chromosome cline was about 4.6 km only. As the former estimate is comparable to that of the Danish portion of the zone (assumed to be much younger than the central European one), zone width does not appear to be related to its age. The strength (B) of the central barrier was estimated as about 20 km; with dispersal (sigma) of about 1 km/gen(1/2), this means effective selection (s*) is approximately 0.06-0.09 for autosomal loci and about 0.25 for X-linked loci. The number of loci under selection was estimated as N= 56-99 for autosomes and about 380 for X-linked loci. Finally, we highlight some potential pitfalls in hybrid zone analyses and in comparisons of different transects. We suggest that conclusions about parts of the mouse genome involved in reproductive isolation and speciation should be drawn with caution and that analytical approaches always providing some estimates should not be used without due care regarding the support or confidence of such estimates, especially if conclusions are based on the difference between these estimates. Finally, we recommend that analysis in two-dimensional space, dense sampling, and rigorous treatment of data, including inspection of likelihood profiles, are essential for hybrid zone studies.     

Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation

We performed a combined analysis of mitochondrial DNA (mtDNA) and microsatellite loci among lake whitefish (Coregonus clupeaformis) populations in order to assess the levels of congruence between both types of markers in defining patterns of genetic structuring, introgressive hybridization and inferring population origins in the hybrid zone of the St. John River basin. A second objective was to test the hypothesis that secondary contact between glacial lineages always resulted in the occurrence of sympatric dwarf and normal whitefish ecotypes. Fish were sampled from 35 populations and polymorphism was screened at mtDNA and six microsatellite loci for a total of 688 and 763 whitefish, respectively. Four lakes harbouring a single whitefish population of normal ecotype admixed with mtDNA haplotypes of different lineages were found. This confirmed that secondary contact between whitefish evolutionary lineages did not always result in the persistence of reproductively isolated ecotypes. Microsatellites further supported the definition of distinct glacial lineages by identifying lineage-specific allelic size groups. They also further supported the hypothesis that ecotypes originated from either a single founding lineage (sympatric divergence) or following secondary contacts between lineages (allopatric divergence), depending on the lake. In general, however, the pattern of population differentiation and introgressive hybridization observed at microsatellites was in sharp contrast with that depicted by mtDNA variation. Both factorial correspondence analysis and analysis of admixture proportion revealed a much more pronounced pattern of introgressive hybridization than depicted by mtDNA analyses. Variable levels of introgression indicated that environmental differences may be as important as the historical contingency of secondary contact in explaining the persistence of sympatric ecotypes and the differential pattern of introgressive hybridization among lakes. Whitefish populations from the St. John River basin hybrid zone represent a rare illustration of a continuum of both morphological and genetic differentiation within a given taxon, spanning from complete introgression to possibly complete reproductive isolation, depending on lakes. Thus, each lake may be viewed as a different temporal snapshot taken throughout the gradual process of speciation.     

Introgression at differentially aged hybrid zones in red-tailed chipmunks

Hybrid zones allow us to investigate the maintenance and the break down of reproductive isolation; they are a window into the speciation process. Tamias ruficaudus (red-tailed chipmunk) has a roughly ring-like distribution in the Inland Northwest and includes two morphologically well-differentiated subspecies, T. r. ruficaudus (in the eastern portion of its range) and T. r. simulans (in the western portion). These taxa meet at two contact zones: the Lochsa River in Idaho and 200 km to the north, near Whitefish, Montana. The Lochsa Zone is encompassed within the Clearwater River Drainage, which has been proposed as a glacial refugium for many taxa throughout the Pleistocene, whereas the Whitefish Zone was under the Cordilleran ice sheet during the most recent glacial maxima approximately 10,000 years ago. Mitochondrial DNA introgression has been documented at both contact zones, yet the subspecies remain significantly distinct with respect to bacular morphology and no intermediate morphologies have ever been reported. Here, we elucidate differentiation and introgression of the nuclear genome using ten microsatellite loci and compare findings to previously described mitochondrial DNA haplotype distribution and introgression. We found significant substructure in the nuclear data; each subspecies is divided into at least two genetically distinct demes. At the Lochsa contact zone, individuals restricted to the mtDNA zone of introgression form a distinct deme at microsatellite loci whereas in the younger, Whitefish contact zone, there is no hybrid-zone specific group. The genetic distances of the demes within these two subspecies indicate recent northward expansion.     

Genetic analysis of the hybridization zone between two subspecies Mus musculus domesticus and Mus musculus musculus in Bulgaria

The hybrid zone between the two subspecies of mice Mus musculus domesticus and Mus musculus musculus, which has been studied extensively in Denmark, crosses Europe to the Black Sea through the Alps and the Balkans. Two hundred and seventy-nine animals were captured in 22 localities along a transect across the Balkans. The animals were characterized for seven diagnostic nuclear loci by protein electrophoresis and by restriction pattern analysis of their mitochondrial DNA. The nuclear data show a sharp transition between the two subspecies, most of the variations in allele frequencies (from 0.9 to 0.1) occurring within a 36-km section of the transect. The introgression varies from one locus to the other and is more pronounced, in terms of distance, in M. m. musculus territory. Mitochondrial DNA introgression is important but occurs in one direction only, i.e. from M. m. musculus to M. m. domesticus, while a cytoplasmic transfer from M. m. domesticus to M. m. musculus has been reported. A previous study showed that no Y chromosome introgression occurs. The different behaviour of these three types of markers could be due to the interaction between selection against hybrid genomes and meiotic recombination. Objectively, it would appear that the genes that can introgress are neutral or nearly so and have been separated from deleterious genes they were linked to by recombination. This could explain the differential introgression between autosomal loci. The mitochondrial and Y chromosomes undergo no or very little recombination and each is transmitted as a whole. Their degree of introgression is thus indicative of the intensity of selection resulting from the amount of functional differentiation between the two taxa, which seems to be strong for the Y chromosome and weak for mitochondrial DNA. We propose that the asymmetry of nuclear introgression is due to different population structures. As M. m. musculus is relatively less structured, the rapid spreading of introgressed genes would be favoured. Such a scheme, however, can hardly account for the unidirectionality of the mitochondrial flow, which could be due to sex-dependent behaviour.     

A hybrid zone with coincident clines for autosomal and sex-specific markers in the Sorex araneus group

In hybrid zones, endogenous counter-selection of hybrids is usually first expressed as reduced fertility or viability in hybrids of the heterogametic sex, a mechanism known as Haldane's rule. This phenomenon often leads to a differential of gene flow between sex-linked markers. Here, we address the possibility of a differential gene flow for Y chromosome, mtDNA and autosomal markers across the hybrid zone between the genetically and chromosomally well-differentiated species Sorex antinorii and Sorex araneus race Vaud. Intermarker comparison clearly revealed coincidental centre and very abrupt clines for all three types of markers. The overall level of genetic differentiation between the two species must be strong enough to hinder asymmetric introgression. Cyto-nuclear mismatches were also observed in the centre of hybrid zone. The significantly lower number of mismatches observed in males than in females possibly results from Y chromosome-mtDNA interactions. Results are compared with those previously reported in another hybrid zone between S. antinori and S. araneus race Cordon.     

The evolving male: spinner dolphin (Stenella longirostris) ecotypes are divergent at Y chromosome but not mtDNA or autosomal markers

The susceptibility of the Y chromosome to sexual selection may make this chromosome an important player in the formation of reproductive isolating barriers, and ultimately speciation. Here, we investigate the role of the Y chromosome in phenotypic divergence and reproductive isolation of spinner dolphin (Stenella longirostris) ecotypes. This species contains six known ecotypes (grouped into four subspecies) that exhibit striking differences in morphology, habitat and mating system, despite having adjacent or overlapping ranges and little genetic divergence at previously studied mtDNA and autosomal markers. We examined the phylogeographic structure for all six ecotypes across the species range (n = 261, 17 geographic locations) using DNA sequences from three Y chromosome markers, two maternally inherited mitochondrial (mtDNA) markers, and a biparentally inherited autosomal intron. mtDNA and autosomal analyses revealed low divergence (most Φ_ST values <0.1) between ecotypes and geographic regions, concordant with previous studies. In contrast, Y intron analyses revealed fixed differences amongst the three most phenotypically divergent groups: S. l. longirostris vs. S. l. roseiventris vs. combined S. l. orientalis/S. l. centroamericana/Tres Marias ecotypes). Another ecotype (whitebelly), previously postulated to be a hybrid between the two phenotypically most divergent ecotypes, had Y haplotypes from both putative parent ecotypes, supporting a hybrid designation. Reduced introgression of the Y chromosome has previously been observed in other organisms ranging from insects to terrestrial mammals, and here we demonstrate this phenomenon in a marine mammal with high dispersal capabilities. These results indicate that reduced introgression of the Y chromosome occurs in a wide taxonomic range of organisms and support the growing body of evidence that rapid evolution of the Y chromosome is important in evolutionary diversification.     

Asymmetric and differential gene introgression at a contact zone between two highly divergent lineages of field voles (Microtus agrestis)

Secondary contact zones have the potential to shed light on the mode and rate at which reproductive isolation accumulates during allopatric speciation. We investigated the population genetics of a contact zone between two highly divergent lineages of field voles (Microtus agrestis) in the Swiss Jura mountains. To shed light on the processes underlying introgression, we used maternally, paternally, and bi-parentally inherited markers. Though the two lineages maintained a strong genetic structure, we found some hybrids and evidence of gene flow. The extent of introgression varied with the mode of inheritance, being highest for mtDNA and absent for the Y chromosome. In addition, introgression was asymmetric, occurring only from the Northern to the Southern lineage. Both patterns seem parsimoniously explained by neutral processes linked to differences in effective sizes and sex-biased dispersal rates. The lineage with lower effective population size was also the more introgressed, and the mode-of-inheritance effect correlated with the male-biased dispersal rate of microtine rodents. We cannot exclude, however, that Haldane's effect contributed to the latter, as we found a marginally significant deficit in males (the heterogametic sex) among hybrids. We propose a possible demographic scenario to account for the patterns documented, and empirical extensions to further investigate this contact zone.     

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.